Fatty acid composition of beef is associated with exonic nucleotide variants of the gene encoding FASN

Genetic associations of fatty acid composition with exonic single nucleotide polymorphisms (SNPs) in the gene encoding fatty acid synthase (FASN) were examined using 513 Korean cattle. All five individual SNPs of g.12870 T>C, g.13126 T>C, g.15532 C>A, g.16907 T>C and g.17924 G>A were associated with a variety of fatty acid compositions and further with marbling score (P < 0.05). Their genotypes of CC, TT, AA, TT, and GG were associated with increased monounsaturated fatty acids and with decreased saturated fatty acids (P < 0.05). The genotypes at all the SNPs also increased marbling score (P < 0.05). Further genetic associations with fatty acid composition suggested that homozygous genotype with the haplotype of ATG at g.15532, g.16907, and g.17924 in a linkage disequilibrium block increased monounsaturated fatty acids and marbling score (P < 0.05). We concluded that the five exonic SNPs of g.12870, g.13126, g.15532, g.16907, and g.17924 in the FASN gene could change fatty acid contents. Their genotypes of CC, TT, AA, TT, and GG and haplotype of ATG at g.15532, g.16907, and g.17924 were recommended for genetic improvement of beef quality.     

A post‐GWAS confirming the SCD gene associated with milk medium‐ and long‐chain unsaturated fatty acids in Chinese Holstein population

The stearoyl‐CoA desaturase (delta‐9‐desaturase) gene encodes a key enzyme in the cellular biosynthesis of monounsaturated fatty acids. In our initial genome‐wide association study (GWAS) of Chinese Holstein cows, 19 SNPs fell in a 1.8‐Mb region (20.3–22.1 Mb) on chromosome 26 underlying the SCD gene and were highly significantly associated with C14:1 or C14 index. The aims of this study were to verify whether the SCD gene has significant genetic effects on milk fatty acid composition in dairy cattle. By resequencing the entire coding region of the bovine SCD gene, a total of six variations were identified, including three coding variations (g.10153G>A, g.10213T>C and g.10329C>T) and three intronic variations (g.6926A>G, g.8646G>A and g.16158G>C). The SNP in exon 3, g.10329C>T, was predicted to result in an amino acid replacement from alanine (GCG) to valine (GTG) in the SCD protein. An association study for 16 milk fatty acids using 346 Chinese Holstein cows with accurate phenotypes and genotypes was performed using the mixed animal model with the proc mixed procedure in sas 9.2. All six detected SNPs were revealed to be associated with six medium‐ and long‐chain unsaturated fatty acids (P = 0.0457 to P < 0.0001), specifically for C14:1 and C14 index (P = 0.0005 to P < 0.0001). Subsequently, strong linkage disequilibrium (D′ = 0.88–1.00) was observed among all six SNPs in SCD and the five SNPs (rs41623887, rs109923480, rs42090224, rs42092174 and rs42091426) within the 1.8‐Mb region identified in our previous GWAS, indicating that the significant association of the SCD gene with milk fatty acid content traits reduced the observed significant 1.8‐Mb chromosome region in GWAS. Haplotype‐based analysis revealed significant associations of the haplotypes encompassing the six SCD SNPs and one SNP (rs109923480) in a GWAS with C14:1, C14 index, C16:1 and C16 index (P = 0.0011 to P < 0.0001). In summary, our findings provide replicate evidence for our previous GWAS and demonstrate that variants in the SCD gene are significantly associated with milk fatty acid composition in dairy cattle, which provides clear evidence for an increased understanding of milk fatty acid synthesis and enhances opportunities to improve milk‐fat composition in dairy cattle.     

Superior single nucleotide polymorphisms that contribute to two main routes of the fatty acid synthesis pathway in Korean cattle

The fatty acid composition of adipose tissue in livestock has been recognized as an important carcass trait that affects meat quality. To determine the overall positive effect on fatty acid synthesis, we investigated PPARγ, SREBPs, FABP4, SCD, and FASN genes because SNP and/or SNP combinations were strongly affecting each fatty acid synthesis pathway. We screened their SNPs, SNP combinations, and genotype. Furthermore, we studied the economic traits C18:1, monounsaturated fatty acid (MUFA), and marbling score (MS). To enhance the accuracy of the predictive genetic effect, we applied a statistical adjustment model excluding environmental factors. We also selected superior genes, gene combinations, and genotypes for each economic trait by using multifactor dimensionality reduction. Based on the results, g.3977-325 T>C (CC) was selected as the best genotype of the SNP and the combination (g.10153 A>G, g.3977-325 T>C) (AACC, AGCC, GGCC), (g.3690 G>A, g.3977-325 T>C) (GGCC, GACC), (g.3977-325 T>C, g.25670 C>T) (CCCC, CCCT, CCTT) (genotypes of the best SNP combination) was selected as the best gene combination for C18:1 and MUFA. In addition, g.3977-325 T>C (CC) was selected as the best genotype of the SNP and the combination (g.3977-325 T>C, g.6974 G>A) (TTAA, CTAA, CCAA, CCGA, CCGG) as the best SNP combination for MS. We integrated genes that were selected as excellent SNP and SNP combinations affecting two main routes of the fatty acid synthesis pathway and identified the best and final SNP, SNP combination, and superior genotype.     

Single nucleotide polymorphisms,haplotypes and combined genotypes of lactoferrin gene and their associations with mastitis in Chinese Holstein cattle

Lactoferrin (Lf) is naturally produced by the mammary gland, having biological functions of antibacterial and anti-inflammatory activities. To investigate whether the Lf gene is associated with mastitis in dairy cattle, a DNA sequencing approach was used to identify single nucleotide polymorphisms (SNPs) in the gene. Three previously reported SNPs in the 5′ flanking region and one novel SNP in exon1 of Lf gene were identified. A total of 353 individuals from Holstein cattle populations were genotyped for their SNPs using Created Restriction Site PCR (CRS-PCR) and PCR-RFLP methods. Twenty-two and nineteen combinations of three SNPs (g.3440T>G, g.3879_3880insG, and g.4432T>C) and another three SNPs (g.3429G>A, g.3440T>G, g.3879_3880insG) were observed, respectively. The result of haplotype analysis of four SNPs showed that fourteen different haplotypes were identified. Two major haplotypes (GECB and GECA) occurred with a frequency of 22.5 and 18.5% in the study population, respectively. Statistical analyses revealed no significant association between one single SNP of Lf gene and SCS, whereas significant associations between their combined genotypes of three SNPs, haplotype and SCS. Combined genotype EFCDBB and GGEFDD with the lowest SCS were favorable for the mastitis resistance. They may be used as a possible candidate for marker-assisted selection in dairy cattle breeding program.     

Genetic variants of the <Emphasis Type="Italic">FABP4</Emphasis> gene are associated with marbling scores and meat quality grades in Hanwoo (Korean cattle)

This study was aimed to search new genetic variants in the bovine FABP4 gene as molecular markers for meat quality and carcass traits. PCR–RFLP analysis revealed that three SNPs located at nucleotide positions g.2834C>G, g.3533T>A, and g.3691G>A were identified based on a GenBank accession number (NC_007312.4). Sequence analysis revealed that SNPs were located in intron 1 (g.2834C>G) and 2 (g.3533T>A), and an exon 3 (g.3691G>A), showing allele frequencies as 0.592, 0.579, and 0.789, respectively. Genetic variabilities of heterozygosity (He) and polymorphic information contents (PIC) were estimated for g.2834C>G (0.608 and 0.531), g.3533T>A (0.615 and 0.539), and g.3691G>A (0.498 and 0.401) loci, respectively. A SNP located in the exon 3 of FABP4 was characterized and associated with desirable increases of MS (marbling scores) and MG (meat quality grades) in Hanwoo. The statistical analysis revealed that additive effects by GG genotypes in g.3691G>A SNP were significantly greater than AA genotypes in MS and MG traits. These findings suggest that the FABP4g.3691G>A SNP will be a useful candidate locus to maximize economic benefits for cattle populations.     

Identification of novel single nucleotide polymorphisms (SNPs) of the lipoprotein lipase (LPL) gene associated with fatty acid composition in Korean cattle

The lipoprotein lipase (LPL) gene can be considered a functional candidate gene that regulates fatty acid composition. In this study, genetic associations between fatty acid composition and exonic single nucleotide polymorphisms (SNPs) in the LPL gene were examined using 612 Korean cattle. We investigated the relationship between unsaturated fatty acids and five novel SNPs (c.322G>A, c.329A>T, c.527T>G, c.988C>T and c.1591G>A), and confirmed that three polymorphic SNPs (c.322G>A, c.329A>T and c.1591G>A) were associated with fatty acid composition. Korean cattle with an AA genotype of c.322G>A, c.329A>T, and GA genotype of c.1591G>A had higher levels of monounsaturated fatty acids and carcass traits (P < 0.05). Our findings confirmed that three novel SNPs we identified in the LPL gene can affect fatty acid composition and carcass traits. Therefore, selection for AA and GA genotypes should be recommended to genetically improve beef quality and flavor.     

Polymorphisms in calpastatin and mu‐calpain genes are associated with beef iron content

The objective of this study was to assess the association of markers in the calpastatin and mu‐calpain loci with iron in beef cattle muscle. The population consisted of 259 cross‐bred steers from Beefmaster, Brangus, Bonsmara, Romosinuano, Hereford and Angus sires. Total iron and heme iron concentrations were measured. Markers in the calpastatin (referred to as CAST) and mu‐calpain (referred to as CAPN4751) genes were used to assess their association with iron levels. The mean and standard error for iron and heme iron content in the population was 35.6 ± 1.3 μg and 27.1 ± 1.4 μg respectively. Significant associations (P < 0.01) of markers were observed for both iron and heme iron content. For CAST, animals with the CC genotype had higher levels of iron and heme iron in longissimus dorsi muscle. For CAPN4751, individuals with the TT genotype had higher concentrations of iron and heme iron than did animals with the CC and CT genotypes. Genotypes known to be associated with tougher meat were associated with higher levels of iron concentration.     

Gene–gene interactions of fatty acid synthase (FASN) using multifactor-dimensionality reduction method in Korean cattle

We examined the gene–gene interactions of five exonic single nucleotide polymorphisms (SNPs) in the gene encoding fatty acid synthase using 513 Korean cattle and using the model free and the non-parametrical multifactor dimensionality reduction method for the analysis. The five SNPs of g.12870 T>C, g.13126 T>C, g.15532 C>A, g.16907 T>C and g.17924 G>A associated with a variety of fatty acid compositions and marbling score were used in this study. The two-factor interaction between g.13126 T>C and g.15532 C>A had the highest training-balanced among the five-factor models and a testing-balanced accuracy at 70.18 % on C18:1 with a cross-validation consistency of 10 out of 10. Also, the two-factor interaction between g.13126 T>C and g.15532 C>A had the highest testing-balanced accuracy at 68.59 % with a 10 out of 10 cross-validation consistency, than any other models on MUFA. In MS, a single SNP g.15532 C>A had the best accuracy at 58.85 % and the two-factor interaction model g.12870 T>C and g.15532 C>A had the highest testing-balanced accuracy at 64.00 %. The three-factor interaction model g.12870 T>C, g.13126 T>C and g.15532 C>A was recorded as having a high testing-balanced accuracy of 63.24 %, but it was lower than the two-factor interaction model. We used likelihood ratio tests for interaction, and Chi square tests to validate our results, with all tests showing statistical significance. We also compared this with mean scores between the high-risk trait group and low-risk trait group. The genotypes of TTCA, TTAA and TCAA at g.15532 and g.13126 on C18:1, genotypes TTCC, TTCA, TTAA, TCAA CCAA at g.15532 and g.13126 on MUFA and genotypes CCCC, TCCA, CCCA, TTAA, TCAA and CCAA at g.15532 and g.12870 on MS were recommended for the genetic improvement of beef quality.     

The polymorphisms of bovine <Emphasis Type="Italic">VEGF</Emphasis> gene and their associations with growth traits in Chinese cattle

PCR-SSCP and DNA sequencing methods were employed to screen the genetic variation of VEGF gene in 671 individuals belonging to three Chinese indigenous cattle breeds including Nanyang, Jiaxian Red and Qinchuan. Three haplotypes (A, B and C), four observed genotypes (AA, AB, BB and AC) and three new SNPs (6765T>C ss130456744, 6860A>G ss130456745, 6893T>C ss130456746) were detected. The analysis suggested that one SNP (ss130456744) in the bovine VEGF gene had significant effects on birth weight, body weight and heart girth at 6 months old in the Nanyang breed (P < 0.05). The results showed that the SNP (ss130456744) in intron 2 of the VEGF gene is associated with early development and growth of Chinese cattle. These findings raise hope that this polymorphism can be a molecular breeding marker in breeding strategies through marker assisted selection (MAS) in Chinese domestic cattle.     

Genetic association of <Emphasis Type="Italic">Phosphodiesterase 1B</Emphasis> (PDE1B) with carcass traits in Korean cattle

Quantitative trait loci for fat deposition and carcass traits have been identified in the vicinity of the gene encoding phosphodiesterase 1B (PDE1B) on bovine chromosome 5. Therefore, the PDE1B gene can be considered as a positional and functional candidate gene for carcass traits in beef cattle. This study aimed to identify single nucleotide polymorphisms (SNPs) in the PDE1B gene and to evaluate their associations with carcass traits in Korean cattle. Eight SNPs, g.440T>G, g.17122A>G, g.17507A>C, g.17575A>G, g.17607T>C, g.17609C>A, g.17692C>T, and g.17707C>G, were identified in the region ranged from exon 1 to intron 6. Five of them were used for association analysis because of their availability of restriction fragment length polymorphisms. As a result, g.17122A>G in intron 3 was significantly associated with backfat thickness (BFT), and g.17507A>C in exon 5 was associated with longissimus dorsi muscle area (LMA, P < 0.05). Animals with the AG genotype of g.17122 had thicker BFT than those with the AA genotype. Animals with the AA or AC genotype of g.17507A>C had larger LMA than those with the CC genotype. We suggested the PDE1B gene as a candidate gene for carcass traits of beef cattle. Fine mapping would be required for application to marker-assisted selection.     

Identification of single nucleotide polymorphisms in the ovine acetyl-CoA carboxylase-alpha gene

Acetyl-CoA carboxylase (ACACA) is the rate-limiting enzyme in the biosynthesis of palmitic acid and long-chain fatty acids. The dietary intake of palmitic acid, which represents approximately 22% of sheep milk fatty acids, increases low-density lipoprotein (LDL) levels and the risk of developing human cardiovascular diseases. Following the candidate gene approach for improving sheep milk composition, and as a first step in assessing the possible influence of the ovine ACACA gene on milk fatty acid composition and its potential use as an animal genetic model of human atherosclerosis disease, we present here an investigation into the genetic variability of the ovine ACACA gene. We sequenced approximately 6.6 kb of ovine ACACA cDNA, including most of the coding sequence of the protein (except 348 bp), in Spanish Churra sheep. A total of 22 synonymous single nucleotide polymorphisms (SNPs) were identified in the analysed sequence, which were genotyped in a set of eight sheep breeds with different productive aptitudes (dairy, meat and double aptitudes). Two of the SNPs identified, SNP03 (c.1450T>C) and SNP15 (c.5134T>C), which appeared to be breed-specific variations, were situated in the gene sequence coding for the biotin-carboxylase (BC) and acetyl-CoA carboxyltransferase (ACCT) domains of the protein, respectively. Particularly interesting is SNP12 (c.4579G>A), which displayed higher frequencies in the dairy-specialised breeds relative to the meat-producing breeds. Moreover, in the dairy breeds studied, the frequency of this SNP showed a positive correlation with the degree of dairy specialisation. A previously described alternative splicing site (Ser-1200) affecting an important regulatory region of the enzyme was observed in one of the Churra animals. Despite the high genetic variability observed in this gene, none of the identified SNPs caused an amino acid change. However, these polymorphisms could be in linkage disequilibrium with other mutations showing a functional effect on the ACACA enzyme. Hence, the characterisations of the allelic variants reported herein lay the groundwork for evaluation of the potential use of these SNPs as genetic markers of fat content and fatty acid composition in sheep dairy products.     

Association of single nucleotide polymorphisms in the microRNA miR‐1596 locus with residual feed intake in chickens

MicroRNAs are an abundant class of small non‐coding RNAs that regulate gene expression. Genetic variations in microRNA sequences may be associated with phenotype differences by influencing the expression of microRNAs and/or their targets. This study identified two single nucleotide polymorphisms (SNPs) in the genomic region of the microRNA miR‐1596 locus of chicken. Of the two SNPs, one was 95 bp upstream of miR‐1596 (g.5678784A>T) and the other was in the middle of the sequence producing the mature microRNA gga‐miR‐1596‐3p (g.5678944A>G). Genotypic distribution of the two SNPs had large differences among 12 chicken breeds (lines), especially between the fast‐growing commercial lines and the slow‐growing Chinese indigenous breeds for the g.5678784A>T SNP. Only the g.5678784A>T SNP was significantly associated with residual feed intake (RFI) in the F₂ population derived from a fast‐growing and a slow‐growing broiler as well as in the pure Huiyang bearded chicken. The birds with the AA genotype of the g.5678784A>T SNP had lower RFI and higher expression of the mature gga‐miR‐1596‐3p microRNA of miR‐1596 than did those with the other genotypes of the same SNP. We also found that the expression of the mature gga‐miR‐1596‐3p microRNA of miR‐1596 was significantly associated with RFI. These findings suggest that miR‐1596 can become a candidate gene related to RFI, and its genetic variation may contribute to changes in RFI by altering expression levels of the mature gga‐miR‐1596‐3p microRNA in chicken.     

Effects of hydroxy‐delta‐5‐steroid dehydrogenase, 3 beta‐ and steroid delta‐isomerase 1 polymorphisms on fat androstenone level and gene expression in Duroc pigs

A high level of androstenone in porcine adipose tissue is a major factor contributing to boar taint. Porcine hydroxy‐delta‐5‐steroid dehydrogenase, 3 beta‐ and steroid delta‐isomerase 1 (3β‐HSD, also known as HSD3B1) plays a key role in the hepatic metabolism that catalyzes androstenone to β‐androstenol. Therefore, 3β‐HSD is a candidate gene for boar taint. This study aimed to investigate functional 3β‐HSD polymorphisms in Duroc pigs. We found eight single nucleotide polymorphisms (SNPs) in the full‐length porcine 3β‐HSD. Four of the SNPs had restriction enzyme sites, and we genotyped them in 147 uncastrated male Duroc pigs using a polymerase chain reaction–restriction fragment length polymorphism method. Pigs with the GG genotype at the g.165262G>A locus (SNP5) had significantly lower androstenone levels than did those with other genotypes (P = 0.030). SNP5 also was associated with differences in 3β‐HSD mRNA levels: pigs with the GG genotype had higher levels than those with other genotypes (P = 0.019). The SNP5 polymorphism could affect the hepatic catabolism of androstenone and consequently impact androstenone accumulation in the adipose tissue. Therefore, SNP5 in the 3β‐HSD of Duroc pigs could be a useful selective marker for decreasing boar taint.     

Whole‐genome association study for fatty acid composition of oleic acid in Japanese Black cattle

Fatty acid composition, especially oleic acid (C18:1), plays an important role in the eating quality of meat in Japanese Black cattle. Therefore, the objective of this study was to identify loci associated with C18:1 in the intramuscular fat of the trapezius muscles in Japanese Black cattle using the Illumina BovineSNP50 BeadChip whole genome single nucleotide polymorphism (SNP) assay. We also evaluated the relationship between C18:1 and three fatty acid synthesis genes, fatty acid synthase (FASN), stearoyl‐CoA desaturase and sterol regulatory element‐binding protein‐1. In this experiment, we applied a mixed model and Genomic Control approach using selective genotyping to perform a genome‐wide association study. A total of 160 animals (80 animals with higher values and 80 animals with lower values), selected from 3356 animals based on corrected phenotype, were genotyped using the Illumina BovineSNP50 BeadChip and three fatty acid synthesis genes, and the quality of these SNPs was assessed. In this study, a total of 38 955 SNPs, which included SNPs in the three fatty acid synthesis genes, were used, and the estimated inflation factor was 1.06. In the studied population, a total of 32 SNPs, including the FASN gene, had significant effects, and in particular 30 SNPs of all significant SNPs were located between 49 and 55 Mbp on chromosome 19. This study is one of the first genome‐wide association studies for fatty acid composition in a cattle population using the recently released Illumina BovineSNP50 BeadChip.     

Mutations of <Emphasis Type="Italic">MC4R</Emphasis> gene and its association with economic traits in Qinchuan cattle

MC4R belongs to a seven-transmembrane G-protein-coupled receptor which may regulate body composition and insulin action. Many mutations in the MC4R gene are associated with obesity, energy expenditure and serum triglyceride levels in human and animals. Six mutations in the MC4R gene were identified in our study (-293C>G, -193A>T, -192T>G, -129A>G, -84T>C and 1,069C>G). The -129A>G was significantly associated with live weight (LW) (P < 0.05), Cattle with the genotypes AG and GG had higher LW than genotype AA. The 1,069C>G was significantly associated with LW, carcass weight (CW), backfat thickness and marbling score (MS). Cattle with the genotype GG had higher LW, CW and MS than genotype CC; Cattle with the genotypes GG and CG had higher MS than CC. The results suggested that -129A>G and 1,069C>G SNP of the MC4R gene may be useful as a genetic marker for carcass and meat quality traits in Qinchuan cattle.     

Identification of single nucleotide polymorphisms in the bovine follicle‐stimulating hormone receptor and effects of genotypes on superovulatory response traits

In dairy cows, there is evidence that failure to respond to superovulation protocols is a heritable trait. In women, genotyping for the p.N680S single nucleotide polymorphism (SNP) in the follicle‐stimulating hormone receptor (FSHR) gene may help identify poor responders before ovarian stimulation is initiated. Our objectives were to identify SNPs in the coding region of the bovine FSHR gene and to investigate the effect of FSHR genotypes on superovulatory response in Holstein cattle. Sequencing of FSHR exons 1–10 revealed seven SNPs. Three were non‐synonymous mutations (c.337C>G, c.871A>G and c.1973C>G). SNP c.337C>G encodes for a proline‐to‐alanine (p.Pro113Ala) amino acid replacement in the extracellular ligand‐binding domain of the receptor. PCR‐RFLP analyses showed that homozygous GG Holstein cows present a higher percentage of viable embryos, whereas GG and CG animals have less unfertilised oocytes. SNP c.871A>G results in an isoleucine‐to‐valine (p.Ile291Val) modification, and homozygous AA animals present lower embryo yield after superovulatory treatments. SNP c.1973C>G corresponds to a threonine‐to‐serine (p.The658Ser) modification in the intracellular carboxyl‐terminal domain of the FSHR protein, and homozygous GG Holstein cows were associated with a lower embryo yield and a higher percentage of unfertilised oocytes. Our results suggest that specific alleles of the bovine FSHR gene are associated with variations in embryo yield and in the number of unfertilised oocytes.     

Identification of genetic associations of ECHS1 gene with milk fatty acid traits in dairy cattle

Our previous genome‐wide association study identified 83 genome‐wide significant SNPs and 20 novel promising candidate genes for milk fatty acids in Chinese Holstein. Among them, the enoyl‐CoA hydratase, short chain 1 (ECHS1) and enoyl‐CoA hydratase and 3‐hydroxyacyl CoA dehydrogenase (EHHADH) genes were located near two SNPs and one SNP respectively, and they play important roles in fatty acid metabolism pathways. We herein validated whether the two genes have genetic effects on milk fatty acid traits in dairy cattle. By re‐sequencing the full‐length coding region, partially adjacent introns and 3000 bp up/downstream flanking sequences, we identified 12 SNPs in ECHS1: two in exons, four in the 3′ flanking region and six in introns. The g.25858322C>T SNP results in an amino acid replacement from leucine to phenylalanine and changes the secondary structure of the ECHS1 protein, and single‐locus association analysis showed that it was significantly associated with three milk fatty acids (P = 0.0002–0.0013). The remaining 11 SNPs were found to be significantly associated with at least one milk fatty acid (P = <0.0001–0.0040). Also, we found that two haplotype blocks, consisting of nine and two SNPs respectively, were significantly associated with eight milk fatty acids (P = <0.0001–0.0125). However, none of polymorphisms was observed in the EHHADH gene. In conclusion, our findings are the first to indicate that the ECHS1 gene has a significant genetic impact on long‐chain unsaturated and medium‐chain saturated fatty acid traits in dairy cattle, although the biological mechanism is still undetermined and requires further in‐depth validation.     

Association study of single nucleotide polymorphisms in JAK2 and STAT5B genes and their differential mRNA expression with mastitis susceptibility in Chinese Holstein cattle

The JAK–STAT pathway plays a key role in mediating immune responses. The genetic effects of single nucleotide polymorphisms (SNPs) in JAK2 and STAT5B were investigated for serum cytokines, mastitis indicators and productions traits in a population of 468 Chinese Holstein cattle. Pooled DNA sequencing revealed one SNP (BTA8:g.39645396A>G) in JAK2 and two SNPs (BTA19:g.43673888A>G and BTA19:g.43660093T>C) in STAT5B. A fixed effect model considering the effects of SNPs, parity, herd, season and year of calving was used by way of the general linear model procedure of sas . Genotype frequencies of these SNPs in the population were in Hardy–Weinberg equilibrium (P > 0.05). A novel SNP (g.39645396A>G) in JAK2 was predicted to change the amino acid from lysine to asparagine and was significantly associated with the somatic cell count (SCC) and somatic cell score (SCS), whereas g.43673888A>G in STAT5B was significantly associated with SCC, SCS and interleukin‐4 (IL‐4) (P < 0.05). The dominant effect of g.39645396A>G in JAK2 was significant for SCS, and its additive effect was significant for SCC, whereas the dominant effect of g.43673888A>G in STAT5B was significant for SCS and IL‐4 (P < 0.05). The combination of g.39645396A>G in JAK2 and g.43673888A>G in STAT5B showed a significant effect on SCC, SCS, IL‐4 and TNF‐α (P < 0.05). As for mRNA expression analysis, the AA genotype g.39645396A>G and GG genotype g.43673888A>G indicated higher mRNA expression level and were significantly different from other genotypes (P < 0.05). The results imply that JAK2 and STAT5B genes could be useful candidate genes, and the identified polymorphisms might potentially be strong genetic markers for selection of dairy cattle against mastitis development.