Novel SNPs in IL-17F and IL-17A genes associated with somatic cell count in Chinese Holstein and Inner-Mongolia Sanhe cattle

Bovine mastitis is the most common and costly disease of lactating cattle worldwide. Apart from milk somatic cell count (SCC) and somatic cell score (SCS), serum cytokines such as interleukin-17 (IL-17) and interleukin-4 (IL-4) may also be potential indicators for bovine mastitis. The present study was designed to investigate the effects of single nucleotide polymorphisms (SNPs) in bovine IL-17F and IL-17A genes on SCC, SCS and serum cytokines in Chinese Holstein and Inner-Mongolia Sanhe cattle, and to compare the mRNA expression variations of the cows with different genotypes. A total of 464 lactating cows (337 Holstein and 127 Inner-Mongolia Sanhe cattle) were screened for SNPs identification and the data were analyzed using fixed effects of herd, parity, season and year of calving by general linear model procedure. The results revealed that SNP g.24392436C > T in IL-17F and SNP g.24345410A > G in IL-17A showed significant effects on SCC and IL-4 in Holstein (n = 337) and on IL-17 and IL-4 in Sanhe cattle (n = 127). The homozygous GG genotype of SNP g.24345410A > G had significantly higher mRNA expression compared with the heterozygous AG genotype. The results indicate that IL-17F and IL-17A could be powerful candidate genes of mastitis resistance and the significant SNPs might be useful genetic markers against mastitis in both dairy and dual purpose cattle.


Background
Bovine mastitis is characterized by a range of chemical, physical and bacteriological changes in the milk accompanied by pathological changes in the udder tissues [1]. Mastitis is the most costly disease of dairy cattle and reported to cause $2 billion dollars annual losses to the U.S. dairy industry and about $35 billion to the world dairy sector [2,3]. Genetics and environment are the most important factors that contribute to mastitis development.
Milk somatic cell count (SCC) and somatic cell score (SCS) are the most suitable indirect indexes to evaluate the degree of mastitis because these traits are convenient, inexpensive and easy to record [4]. The genetic evaluation and indirect selection of cattle for lower SCC or SCS may reduce the incidence of susceptibility to mastitis [5,6]. The cytokines in serum such as interleukin-4 (IL-4), IL-6, IL-17, tumor necrosis factor-α (TNF-α) and interferon-γ (IFN-γ) also act as indirect indexes in inflammatory conditions [7], which suggests that beside SCC and SCS, serum cytokines could be considered as crucial indicators for bovine mastitis. Previous studies indicated that mutations in IL-17F and IL-17A genes were related with inflammatory conditions e.g. inflammatory bowel disease [7], asthma [8], rheumatoid arthritis [9], ovarian cancer [10], colon cancer [11] and breast cancer [12]. The role of the two genes in disease process in humans and other model animals indicates that IL-17F and IL-17A could be potential candidate genes for mastitis resistance in bovine as well.
IL-17 is a bridge of the innate and adaptive immune systems. The IL-17 family contains six members with IL-17A as the founding member (which was cloned about two decades ago) [13] and IL-17F is the most recently discovered member of the family [14]. From the alignment of the predicted amino acid sequence, it was found that IL-17F and IL-17A share the strongest homology compared to other members of the family and have similar functions to induce inflammatory response [15]. Both IL-17F and IL-17A were reported to be associated with increased risk of certain subtypes of gastric cancer [16]. However, the influence on the pathophysiological features of ulcerative colitis is partially different among IL-17A and IL-17F polymorphisms. A study revealed that IL-17A/-197A allele was significantly associated with chronic relapsing phenotype of ulcerative colitis and the -197A/A homozygote was more frequent in steroid dependent cases, while the IL-17F/7488 T allele was associated with the chronic continuous phenotype [17]. The regulation of IL-17 directly influences the nature of the cellular recruitment at the site of inflammation and it represents a marker of molecule regulating neutrophil and eosinophil infiltration [18].
It was reported that pure dairy breeds e.g. Holstein are more susceptible to mastitis than dual purpose breeds e.g. Simmental and Sanhe cattle [4]. Sanhe cattle is a precious dual purpose breed of northern (Sanhe area) China. The Sanhe area is located in the Inner Mongolia between longitude 117°15'~124°02' east and latitude 47°05'~51°30' north. Its average temperature in winter ranges between −20°C to −31°C and the mean annual temperature is −0.95°C [5]. However, the genetic difference of Sanhe cattle with Holsteins and the genetic effects on mastitis resistance in this breed is rarely reported till date. The present study was designed to 1) evaluate IL-17F and IL-17A as candidate genes for association analysis with SCS as well several serum cytokines in Chinese Holstein and Sanhe cattle, and to 2) identify novel potential SNP markers for bovine mastitis resistance which could be used for udder health improvement.

Cattle population and sample collection
A total of 337 Chinese Holstein cows were randomly selected from three dairy cattle farms located in Northern China (Qiqihar, Tianjin and Shanxi) and 127 Sanhe cattle were randomly collected from Sanhe cattle breeding farm in Hailar, Inner Mongolia (Fig. 1). The cows ranged between parity one and five and were milked thrice a day. The cattle were fed a lactation diet as recommended by the Dairy Association of China for lactating cows.
Blood samples were collected from the jugular vein of all selected cattle in three 9 mL tubes, one for DNA extraction (in EDTA coated tube), the second for RNA extraction (containing TRIzol) and the third for serum isolation by an in-house technician of the farm. The RNA extraction tube was immediately placed at −80°C to avoid any damage to the RNA. For serum isolation, the blood samples were placed at room temperature 30 min to enable blood coagulation and then centrifuged at 3,000 r/min for 10 min to separate serum. The serum samples were then stored at 4°C and sent to Beijing Huaying Biological Technology Research Institute within 24 h to detect the concentration of IL-4, IL-6, IL-17, IFN-γ and TNF-α (Sino-UK, China) using radioimmunoassay. Briefly, each serum sample was initially centrifuged for 5 min at 3,000 r/min at 4°C. Then, 100 μL of the supernatant ("cold" antigen), 100 μL of antibody and 100 μL of radiolabeled antigen (125-I, "hot" antigen) were mixed and then stayed at 4°C for 24 h. Next, 500 μL of separating solution was added in and mixed. The mixture was stayed at room temperature for 20 min and then centrifuged for 25 min at 3,500 r/ min at 4°C. Finally, the supernatant was removed and the radioactivity of the bound antigen remaining in the precipitate was measured by a gamma counter. Fresh milk samples were collected in 50 mL tubes from all the four quarters of each cow's udder and mixed with 0.03 g potassium dichromate, stored at 4°C and then sent to Beijing Dairy Cattle Centre for somatic cell count (SCC) test within 48 h. SCS were converted from SCC according to Rupp and Biochard [19] using a formula SCS = log 2 (SCC/100,000) + 3.

DNA isolation and SNPs identification
Genomic DNA was extracted from whole blood using Tiangen Blood DNA Kit (Tiangen Biotech Co., China) following the manufacturer's instructions. The quantity and quality of DNA were measured using NanoDrop™ ND-2000c Spectrophotometer (Thermo Scientific, Inc.).
All SNPs in bovine IL-17F and IL-17A genes were identified by sequencing polymerase chain reaction (PCR) amplicon using a DNA pool constructed with genomic DNA of thirty randomly selected cattle (50 ng/ uL per sample). Based on the information of the NCBI and UCSC database, primers for 12 SNPs were designed to amplify the two genes fragments (Table 1). Two SNPs, one in IL-17F (g.24392436C > T) and another one in IL-17A (g.24345410 A > G), were arbitrarily chosen and screened in the two cattle populations using SNaPshot assay (ABI Multiplex SNaPshot, USA).

RNA isolation, purification and reverse transcription
The Bioteke RNA Isolation Kit (Bioteke, Beijing) was used to extract total RNA from peripheral blood samples of the Sanhe and Holstein cattle. RNase-Free DNaseSet (Qiagen, Germany) was used to purify RNA ensuring that genomic DNA was discarded. Reverse transcription was carried out with High Capacity cDNA Archive Kit (ABI, USA) according to the manufacturer's protocol. PCR primer sets for bovine IL-17A and IL-4 genes were designed by software oligo6.0, considering the golden rules for real-time PCR (RT-PCR) ( Table 1). The amplification efficiency of these primer pairs was tested by RT-PCR firstly. The mRNA expression of the two genes was normalized against the housekeeping gene glyceraldehyde-3-phosphatedehydrogenase (GAPDH) cDNA in the corresponding samples. The primers sequences of GAPDH are also listed in Table 1.

Quantitative RT-PCR
Quantitative real time polymerase chain reaction (qRT-PCR) was carried out to determine the mRNA expression levels of IL-17A and IL-4 genes. The reactions were performed in a total volume of 20 μL containing 2 μL cDNA, 1 μL each primers, 10 μL SYBR Green Master Mix (Roche, Switzerland), 6 μL nuclease-free water using the following amplification condition: 94°C for 10 min, followed by 44 cycles of 94°C for 15 s, 60°C for 10 s, 72°C for 10 s, and 72°C for 30 s. Fluorescence signals were collected during 60°C step. Mean was derived from the two repeats for each sample. Light Cycler 480 RT-PCR system was used to perform amplification, detection and data analyses.

Statistical analyses
The SCC data were classified into three grades: (I) SCC ≤ 200,000 cells/mL; (II) 200,000 cells/mL< SCC < 500,000 cells/mL; and (III) SCC ≥ 500,000 cells/mL. The influence of SCC levels on five serum cytokines and SCS in the two cattle populations was analyzed using the Table 1 Primer pairs of PCR and real time qRT-PCR used in the present study general linear model procedure of SAS 9.1 using the following model: Where y ij represents SCS or serum concentration of cytokine IL-4, IL-6, IL-17, TNF-α and IFN-γ; μ is overall mean; α i is effect of SCC levels; e is the random error.
Associations of the 12 SNPs with SCS and five serum cytokines were analyzed using the GLM model 2 in the two breeds separately (SAS 9.1): Where y ijkl represents each phenotype; μ is overall mean; α i is the fixed effect of genotype; β j is the fixed effect of the herd, year and season of birth; γ l is the fixed effect of parity; e is the random error.
In model 2, the estimated genotype effects were further divided into additive effect (A) and dominant effect (D). The additive effect was the mean deviation of two homozygous genotypes (Formula 1), and the dominant effect was calculated by the deviation of heterozygous genotype from the mean of two homozygous genotypes (Formula 2) [20].
Where, AA, AB and BB were least square means of genotype AA, AB and BB, respectively. Student t test was used for qRT-PCR analyses for comparing the difference of mRNA expression level of IL-17A and IL-4 between different genotypes of IL-17A.

SNP discovery and genotypes of bovine IL-17 F and IL-17A
In the present study, a total of 12 SNPs (comprising 8 SNPs in IL-17 F and 4 SNPs in IL-17A gene) were revealed by screening the pooled DNA of 30 randomly selected Chinese Holstein and Sanhe cattle ( Table 2). Of the 8 SNPs in IL-17 F, one each was located in exon 4 and intron 3 and six were located in 2 kb promoter region, whereas, in IL-17A, three SNPs were located in exon 3 and one in 2 kb promoter region. Out of the 12 SNPs, two were then genotyped in a total population of 337 Holstein and 127 Sanhe cattle (Table 3). Allele and genotype frequencies and Chi square test χ 2 results are summarized in Table 3. Chi square test (χ 2 ) showed that genotype frequencies of all SNPs in the population were in Hardy-Weinberg equilibrium.

Effect of three SCC grades on SCS and cytokines of Chinese Holstein and Sanhe cattle
The descriptive statistics for SCC, SCS and each serum cytokines were listed in the Additional file 1: Table S1. The effect of SCC grades on SCS and serum cytokines were analyzed in the two cattle population ( Table 4). The results showed that SCC grade had highly significant effect on SCS, IL-17 and IFN-γ in Holstein and on SCS, IL-17 and IL-6 in Sanhe cattle (P < 0.001). In Holstein cattle, the grade 3 SCC was significantly associated with higher values of SCS and IL-17 compared with grad I and II. Whereas, in Sanhe cattle, SCC grade II was significantly associated with higher values of cytokine IL-17 and IL-6 compared to grade 1 and 3 (P < 0.001).

Effects of the SNPs on mastitis indicator traits
The results of association study are shown in Table 5. SNP (g.24392436C > T) in IL-17F showed significant association with SCS in Holstein (P < 0.05) and highly significant association with cytokine IL-17 in Sanhe cattle (P < 0.01). Whereas, the association of SNP (g.24345410A > G) in IL-17A was found significant with cytokine IL-4 in both Holstein and Sanhe cattle (P < 0.05). The CC genotype of SNP (g.24392436C > T) was significantly associated with higher SCS compared to the other genotypes in Holstein, whereas, the TT genotype was significantly associated with higher cytokine IL-17 (P < 0.01) than genotype CC and CT. The AG genotype of IL-17A was significantly associated with higher values of IL-4 than AA genotype in both Holstein and Sanhe cattle (P < 0.01).

The additive and dominant effects of the SNPs
To dissect the genotype effects of the significant SNPs, their additive and dominant effects were calculated using  Table 6). In addition, the dominant effect of the SNP g.24345410A > G of IL-17A was found significant on IL-4 in both Holstein and Sanhe cattle (P < 0.01).

Effects of the combination genotypes on mastitis indicator traits
The effects of the combination genotypes are mentioned in Table 7. The results showed that the combination genotype of SNPs in IL-17F and IL-17A were significantly associated with IL-17 cytokine (P < 0.01) and showed a tendency towards significance for association with SCS (P = 0.054) in Sanhe cattle.

mRNA expression level in different genotypes of IL-17A and its association with IL-4 gene, cytokine IL-4 and SCS
The results of qRT-PCR for mRNA expression of IL-17A and IL-4 gene showed that mRNA expression of IL-17A was significantly higher in genotype GG compared with genotype AG (Fig. 2). Moreover, the AG genotype of IL-17A gene showed significantly higher mRNA expression compared with AA and GG genotype with respect to IL-

Discussion
In the present study, a total of 12 SNPs were identified in IL-17F and IL-17A genes in Holstein and Inner-Mongolia Sanhe cattle, of which 2 SNPs were arbitrarily chosen for genotyping and further screening to evaluate their potential association with mastitis. The two SNPs were found significantly associated with mastitis indicator traits. To the best of our knowledge, this is the first study to examine the associations of polymorphisms in IL-17A and IL-17 F genes with bovine mastitis and to evaluate these genes as prognostic markers for mastitis in dairy cattle.
Genetic polymorphisms in IL-17F and IL-17A were reported to be significantly associated with susceptibility  of breast cancer in human [12]. In a bovine model of tuberculosis, a higher expression of IL-17 gene was reported to be positively associated with bovine tuberculosis suggesting IL-17 as a potential biomarker for prognosis in bovine tuberculosis [21]. In the present study, we found that the SNP in IL-17F was significantly associated with SCS in Holstein cattle and the effect of the combination genotype on SCS showed a tendency towards significance in Sanhe cattle. The SNPs in IL-17 F and IL-17A were significantly associated with serum    in mastitis susceptibility studies in dairy cattle. In addition, although SCS is continuous trait which generally serves as an important indicator for subclinical mastitis, is highly influenced by various environmental factors. Thus, the records of clinical incidence of mastitis should be collected and analyze their association with these two SNPs in IL17F and IL17A.
The cytokine IL-6 is a major player in hematopoiesis as well as in the immune system. It possesses both pro-and anti-inflammatory properties and is a pleiotropic inflammatory cytokine involved in numerous biological functions including hematopoiesis, inflammation, immune regulation and oncogenesis [22]. It was demonstrated that detection of cytokine IL-6 in milk indicated subclinical mastitis earlier than the detection of elevated SCC [23]. The author concluded that the detection of IL-6 in milk could be a reliable prediction marker for subclinical mastitis. In the present study, SNPs in IL-17F and IL-17A genes were non-significantly associated with cytokine IL-6, but were significantly associated with SCS, cytokine IL-17 and IL-4. The epigenetic regulation of gene expression by IL-6 can lead to tumor progression by altering the promoter methylation and the genes regulatory pathways [24]. Noticeably, in the present study the grades of SCC were significantly associated with SCS, cytokines IL-6 and IL-17. The significant association of the two SNPs in these genes with serum cytokines and mastitis indicator trait is a positive clue to consider these genes as candidate genes in mastitis resistance studies.

Conclusion
In conclusion, both IL-17A and IL-17F gene polymorphisms (g.24392436C > T and g.24345410A > G) may provide valuable information for predicting the prognosis of bovine mastitis. However, further studies are recommended to validate both IL-17 gene expression and early lymphocyte activation as biomarker of immune status. The present study provides preliminary findings of the relationship between IL-17A and IL-17F gene with SCC/SCS and cytokines levels. The results infer that IL-17A and IL-17F genes could be crucial modifiers of inflammatory diseases and the SNPs might be useful markers of genetic resistance against bovine mastitis development in both dairy and dual purpose cattle. We suggest further in-depth research using large population size to evaluate the association of these genes with bovine mastitis.

Additional file
Additional file 1: Table S1.