Two missense mutations in exon 9 of caprine PRLR gene were associated with litter size

Guanzhong (n = 321) and Boer (n = 191) goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the coding regions of the prolactin receptor (PRLR) gene by DNA sequencing and PCR‐RFLP. Two SNPs (c.1457G>A and c.1645G>A) were identified that caused amino acid variations p.Ser485Asn and p.Val548Met respectively. Statistical results indicated that the c.1457G>A and c.1645G>A SNPs were significantly associated with litter size in Boer and Guanzhong goat breeds. Further analysis revealed that combined genotype C4 (GGGG) and haplotype G‐G were better than the others for litter size in both goat breeds. These results might contribute to goat genetic resources and breeding.     

Two mutations in the 5′‐flanking region of the KITLG gene are associated with litter size of dairy goats

In this study, Xinong Saanen (SN) and Guanzhong (GZ) dairy goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the 5′‐flanking region of the KITLG gene by DNA sequencing and primer‐introduced restriction analysis–polymerase chain reaction. Two novel SNPs (g.13090G>T and g.13664C>A) were identified (GenBank Accession no. KM658964). Furthermore, g.13090G>T and g.13664C>A loci were closely linked in SN and GZ breeds (r² > 0.33). Association analysis results showed that g.13090G>T and g.13664C>A SNPs significantly affected litter size (P < 0.05). The litter size of individuals with the combined genotype GG/CC from both dairy goat breeds was greater than that of individuals with TT/AA in average parity (P < 0.05). Known biochemical and physiological functions, along with our results, indicated that GG/CC could be used in marker‐assisted selection to choose individuals with greater litter size from both breeds. These results extend the spectrum of genetic variation in the caprine KITLG gene and may contribute to genetic resources and breeding of goats.     

Genetic polymorphisms of the bovine Fatty acid binding protein 4 gene are significantly associated with marbling and carcass weight in Hanwoo (Korean Cattle)

The objective of this study was to investigate an association between polymorphisms in the FABP4 gene and phenotypic variation for marbling and carcass weight (CWT) in a population of Hanwoo steers. We re‐sequenced 4.3 kb of the FABP4 gene region in 24 Hanwoo bulls and identified 16 SNPs and 1 microsatellite polymorphism. Of these 16 SNPs, three SNPs [g.2774G>C (intron I), g.3473A>T (intron II) and g.3631G>A (exon III, creating a p.Met >Val amino acid substitution)] were genotyped in 583 steers to assess their association with carcass traits. The g.3473A allele showed a significant increasing effect on CWT (P = 0.01) and the g.3631G allele was associated with higher marbling score (P = 0.006). One haplotype of these three SNPs (CAG) was significantly associated with CWT (P = 0.02) and marbling score (P = 0.05) and could potentially be of value for marker assisted selection in Hanwoo cattle. The CAG haplotype effect for CWT was larger (11.14 ± 5.03 kg) than the largest single locus effect of g.3473A>T (5.01 ± 2.2 kg).     

Detection of exonic variants within the melanocortin 1 receptor (MC1R) gene in Black Silky,White Leghorn and Golden duckwing Araucana chicken

The melanocortin 1 receptor (MC1R) gene can be considered a candidate functional gene for the pigmentation of plumage color. The aim of this study was to investigate the association between the genotype frequencies of g.69 T>C, g.376 G>A and g.427 A>G SNPs within the MC1R gene in Black silky (O), Golden duckwing Araucana (GA) and White Leghorn (W). The CC and AA genotype frequencies of g.69 T>C and g.427 A>G SNPs in White Leghorn (W) were both 1.000, and the TT genotype frequency of the g.69 T>C SNP in Golden duckwing Araucana (GA) was also 1.000. The GG and AA genotype frequencies of g.376 G>A and g.427 A>G SNPs in Black silky (O) were both 0.100. When a haplotype is observed using a combination of markers, a Golden duckwing Araucana (GA) can especially be distinguished when it is a TAG, TGG and TAA type in the SNP combination of the MC1R gene. In case of the CAA types, only White Leghorn (W) could specifically be distinguished. Therefore, three SNPs in MC1R may provide identification in chicken breeds.     

The CDKN2A/p16INK4a 5′UTR sequence and translational regulation: impact of novel variants predisposing to melanoma

Many variants of uncertain functional significance in cancer susceptibility genes lie in regulatory regions, and clarifying their association with disease risk poses significant challenges. We studied 17 germline variants (nine of which were novel) in the CDKN2A 5′UTR with independent approaches, which included mono and bicistronic reporter assays, Western blot of endogenous protein, and allelic representation after polysomal profiling to investigate their impact on CDKN2A mRNA translation regulation. Two of the novel variants (c.‐27del23, c.‐93‐91delAGG) were classified as causal mutations (score ≥3), along with the c.‐21C>T, c.‐34G>T, and c.‐56G>T, which had already been studied by a subset of assays. The novel c.‐42T>A as well as the previously described c.‐67G>C were classified as potential mutations (score 1 or 2). The remaining variants (c.‐14C>T, c.‐20A>G, c.‐25C>T+c.‐180G>A, c.‐30G>A, c.‐40C>T, c.‐45G>A, c.‐59C>G, c.‐87T>A, c.‐252A>T) were classified as neutral (score 0). In conclusion, we found evidence that nearly half of the variants found in this region had a negative impact on CDKN2A mRNA translation, supporting the hypothesis that 5′UTR can act as a cellular Internal Ribosome Entry Site (IRES) to modulate p16^INK4a translation.     

Identification of single nucleotide polymorphisms in the TNFRSF17 gene and their association with gastrointestinal disorders

TNFRSF17 is preferentially expressed in mature B lymphocytes, and may be important for the development of B cells. TNFRSF17 is selected as a candidate susceptibility gene to IBD pathogenesis by our cDNA microarray analysis, and we showed the specific expression of TNFRSF17 in resting and activated CD19⁺ cells obtained from human blood. We identified four SNPs (g-1729G>A, g.2295T>C, g.2445G>A and g.2493G>A) and one variation site (g.894delT) in the TNFRSF17 gene using direct sequencing analysis. In addition, the association of the genotype and allelic frequencies of these SNPs was studied in healthy controls and in patients with ulcerative colitis (UC) or irritable bowel syndrome (IBS). Although, the genotype and allelic frequencies of these SNPs, in the UC and IBS patients, were not significantly different from those in the healthy controls, the distribution of the AAG, GGA, AGG and AAA haplotypes, of the SNPs (g.-1729G>A, g.2445G> A and g.2493G>A) associated with the TNFRSF17 gene, in the UC patients, were notably different from those of the healthy controls (P = 0.002, 0.002, 4.7E-4 and 3.3E-6, respectively). Moreover, the frequencies of the AAG, AGG, GAG and GAA haplotypes were significantly different in the IBS patients compared to the healthy controls (P = 4.2E-5, 4.4E-17, 1.8E-22 and 1.6E-10, respectively). These results suggest that the haplotypes of the TNFRSF17 polymorphisms might be associated with UC and IBS susceptibility.     

Gender‐Specific Association of a Perilipin Gene Haplotype with Obesity Risk in a White Population

Objective: Perilipin is a class of protein‐coating lipid droplets in adipocytes and steroidogenic cells. Our purpose was to examine the association between common single‐nucleotide polymorphisms (SNPs) at the perilipin (PLIN) locus and obesity, as well as related phenotypes, in unrelated American adults. Research Methods and Procedures: Four PLIN SNPs (PLIN 6209T>C, 11482G>A, 13041A>G, and 14995A>T) were typed in 734 white subjects (373 men and 361 women) attending a residential lifestyle intervention program. The baseline anthropometric and biochemical measures were used. Obesity was defined as BMI ≥ 30 kg/m². Results: Multivariate analysis demonstrated that, in women, two of the SNPs (13041A>G, and 14995A>T) were significantly associated with percentage body fat (p = 0.016 for 13041A>G and p = 0.010 for 14995A>T) and waist circumference (p = 0.020 for 13041A>G and p = 0.045 for 14995A>T). Moreover, haplotype analysis using these two SNPs indicated that haplotypes A/T and G/T were both associated with significantly increased obesity risk (odds ratio = 1.76, 95% confidence interval 1.07 to 2.90 for haplotype A/T, and odds ratio = 1.73, 95% confidence interval 1.06 to 2.82 for haplotype G/T) when compared with haplotype A/A. No significant associations between PLIN variations and obesity were found in men. Discussion: Our data support the hypothesis that the PLIN locus may be a significant genetic determinant for obesity risk in whites and that women are more sensitive to the genetic effects of perilipin than men.     

Sample-to-SNP kit: A reliable,easy and fast tool for the detection of HFE p.H63D and p.C282Y variations associated to hereditary hemochromatosis

Classical hereditary hemochromatosis involves the HFE-gene and diagnostic analysis of the DNA variants HFE p.C282Y (c.845G > A; rs1800562) and HFE p.H63D (c.187C > G; rs1799945). The affected protein alters the iron homeostasis resulting in iron overload in various tissues. The aim of this study was to validate the TaqMan-based Sample-to-SNP protocol for the analysis of the HFE-p.C282Y and p.H63D variants with regard to accuracy, usefulness and reproducibility compared to an existing SNP protocol. The Sample-to-SNP protocol uses an approach where the DNA template is made accessible from a cell lysate followed by TaqMan analysis. Besides the HFE-SNPs other eight SNPs were used as well. These SNPs were: Coagulation factor II-gene F2 c.20210G > A, Coagulation factor V-gene F5 p.R506Q (c.1517G > A; rs121917732), Mitochondria SNP: mt7028 G > A, Mitochondria SNP: mt12308 A > G, Proprotein convertase subtilisin/kexin type 9-gene PCSK9 p.R46L (c.137G > T), Plutathione S-transferase pi 1-gene GSTP1 p.I105V (c313A > G; rs1695), LXR g.-171 A > G, ZNF202 g.-118 G > T. In conclusion the Sample-to-SNP kit proved to be an accurate, reliable, robust, easy to use and rapid TaqMan-based SNP detection protocol, which could be quickly implemented in a routine diagnostic or research facility.     

Genome-wide and local pattern of linkage disequilibrium and persistence of phase for 3 Danish pig breeds

Background

A genome wide association study for litter size in Norwegian White Sheep (NWS) was conducted using the recently developed ovine 50K SNP chip from Illumina. After genotyping 378 progeny tested artificial insemination (AI) rams, a GWAS analysis was performed on estimated breeding values (EBVs) for litter size.

Results

A QTL-region was identified on sheep chromosome 5, close to the growth differentiation factor 9 (GDF9), which is known to be a strong candidate gene for increased ovulation rate/litter size. Sequencing of the GDF9 coding region in the most extreme sires (high and low BLUP values) revealed a single nucleotide polymorphism (c.1111G>A), responsible for a Val→Met substitution at position 371 (V371M). This polymorphism has previously been identified in Belclare and Cambridge sheep, but was not found to be associated with fertility. In our NWS-population the c.1111G>A SNP showed stronger association with litter size than any other single SNP on the Illumina 50K ovine SNP chip. Based on the estimated breeding values, daughters of AI rams homozygous for c.1111A will produce minimum 0.46 - 0.57 additional lambs compared to daughters of wild-type rams.

Conclusion

We have identified a missense mutation in the bioactive part of the GDF9 protein that shows strong association with litter size in NWS. Based on the NWS breeding history and the marked increase in the c.1111A allele frequency in the AI ram population since 1983, we hypothesize that c.1111A allele originate from Finnish landrace imported to Norway around 1970. Because of the widespread use of Finnish landrace and the fact that the ewes homozygous for the c.1111A allele are reported to be fertile, we expect the commercial impact of this mutation to be high.     

Identification of SNPs within the sheep <Emphasis Type="Italic">PROP1</Emphasis> gene and their effects on wool traits

Regarding mutations of PROP1 (Prophet of POU1F1) gene significantly associating with combined pituitary hormone deficiency (CPHD) in human patients and animals, PROP1 gene is a novel important candidate gene for detecting genetic variation and growth, reproduction, metabolism traits selection and breeding. The aim of this study was to detect PROP1 gene mutation of the exon 1–3 and its association with wool traits in 345 Chinese Merino sheep. In this study, on the basis of PCR-SSCP and DNA sequencing methods, ten novel SNPs within the sheep PROP1 gene, namely, AY533708: g.45A > G resulting in Glu15Glu, g.1198A > G, g.1341G > C resulting in Arg63Ser, g.1389G > A resulting in Ala79Ala, g.1402C > T resulting in Leu84Leu, g.1424A > G resulting in Asn91Ser, g.1522C > T, g.1556A > T, g.1574T > C, g.2430C > G were reported. In addition, association analysis showed that three genotypes of P4 fragment were significantly associated with fiber diameter in the analyzed population (P = 0.044). These results strongly suggested that polymorphisms of the PROP1 gene could be a useful molecular marker for sheep breeding and genetics through marker-assisted selection (MAS).     

Single-nucleotide polymorphisms g.151435C>T and g.173057T>C in PRLR gene regulated by bta-miR-302a are associated with litter size in goats

Single-nucleotide polymorphisms (SNPs) located at microRNA-binding sites (miR-SNPs) can affect the expression of genes. This study aimed to identify the miR-SNPs associated with litter size. Guanzhong (n = 321) and Boer (n = 191) goat breeds were used to detect SNPs in the caprine prolactin receptor (PRLR) gene by DNA sequencing, primer-introduced restriction analysis-polymerase chain reaction, and polymerase chain reaction-restriction fragment length polymorphism. Three novel SNPs (g.151435C>T, g.151454A>G, and g.173057T>C) were identified in the caprine PRLR gene. Statistical results indicated that the g.151435C>T and g.173057T>C SNPs were significantly associated with litter size in Guanzhong and Boer goat breeds. Further analysis revealed that combinative genotype C6 (TTAACC) was better than the others for litter size in both goat breeds. Furthermore, the PRLR g.173057T>C polymorphism was predicted to regulate the binding activity of bta-miR-302a. Luciferase reporter gene assay confirmed that 173057C to T substitution disrupted the binding site for bta-miR-302a, resulting in the reduced levels of luciferase. Taken together, these findings suggested that bta-miR-302a can influence the expression of PRLR protein by binding with 3′untranslated region, resulting in that the g.173057T>C SNP had significant effects on litter size.     

Association of the variants and haplotypes in the DOCK7, PCSK9 and GALNT2 genes and the risk of hyperlipidaemia

Little is known about the association between the single nucleotide polymorphisms (SNPs) and haplotypes of the dedicator of cytokinesis 7 (DOCK7), pro‐protein convertase subtilisin/kexin type 9 (PCSK9) and polypeptide N‐acetylgalactosaminyltransferase 2 (GALNT2) and serum lipid traits in the Chinese populations. This study was to determine the association between nine SNPs in the three genes and their haplotypes and hypercholesterolaemia (HCH)/hypertriglyceridaemia (HTG), and to identify the possible gene–gene interactions among these SNPs. Genotyping was performed in 733 HCH and 540 HTG participants. The haplotype of C‐C‐G‐C‐T‐G‐C‐C‐G [in the order of DOCK7 rs1168013 (G>C), rs10889332 (C>T); PCSK9 rs615563 (G>A), rs7552841 (C>T), rs11206517 (T>G); and GALNT2 rs1997947 (G>A), rs2760537 (C>T), rs4846913 (C>A) and rs11122316 (G>A) SNPs] was associated with increased risk of HCH and HTG. The haplotypes of C‐C‐G‐C‐T‐G‐C‐C‐A and G‐C‐G‐T‐T‐G‐T‐C‐G were associated with a reduced risk of HCH and HTG. The haplotypes of G‐C‐G‐C‐T‐G‐C‐C‐A and G‐C‐G‐C‐T‐G‐T‐C‐G were associated with increased risk of HCH. The haplotypes of C‐T‐G‐C‐T‐G‐C‐C‐G, G‐C‐A‐C‐T‐G‐C‐C‐G and G‐C‐G‐C‐T‐G‐C‐C‐A were associated with an increased risk of HTG. The haplotypes of G‐C‐G‐C‐T‐G‐T‐C‐A and G‐C‐G‐T‐T‐G‐T‐C‐G were associated with a reduced risk of HTG. In addition, possible inter‐locus interactions among the DOCK7, PCSK9 and GALNT2 SNPs were also noted. However, further functional studies of these genes are still required to clarify which SNPs are functional and how these genes actually affect the serum lipid levels.     

Significant association between SNPs in the superoxide dismutase 3, extracellular (SOD3) gene and resistance to Aeromonas hydrophila in the freshwater mussel Hyriopsis cumingii

Extracellular superoxide dismutase (SOD3) is a major antioxidant enzyme that protects organs from damage by reactive oxygen species (ROS). In this study, the SOD3 gene was identified and characterized from the freshwater mussel Hyriopsis cumingii (Hc‐SOD3). The cDNA sequence consists of 763 bp, encoding a protein of 208 amino acids. The amino acid sequence possesses two CuZnSOD signature sequences, and amino acids required for binding of Cu (His‐93, ‐95, ‐110 and ‐169) and Zn (His‐110, ‐118, ‐129 and Asp‐132) were conserved in Hc‐SOD3. The Hc‐SOD3 genomic sequence was 9165 bp in length, containing four exons and three introns. Eighteen single nucleotide polymorphisms were detected in the Hc‐SOD3 gene from resistant stock (RS) and susceptible stock (SS) of H. cumingii to Aeromonas hydrophila. The genotype and allele distribution were examined in resistant and susceptible stocks. Among them, a C/G substitution at the g.7994C>G locus and G/C substitution at the g.8087G>C locus were significantly associated with resistance/susceptibility of H. cumingii to A. hydrophila, both in genotype (P = 0.017, P = 0.004 respectively) and allele frequency (P = 0.021, P = 0.006 respectively). Linkage disequilibrium analysis revealed that g.7994C>G, g.8001A>G, g.8035G>A, g.8087G>C and g.8191T>A were in linkage disequilibrium. The results suggest that the two polymorphic loci, g.7994C>G and g.8087G>C, could be potential genetic markers for future molecular selection of strains that are resistant to diseases.     

Association of missense variants in GDF9 with litter size in Entlebucher Mountain dogs

In the past two decades, average litter size (ALS) in Entlebucher Mountain dogs decreased by approximately 0.8 puppies. We conducted a GWAS for ALS using the single-step methodology to take advantage of 1632 pedigree records, 892 phenotypes and 372 genotypes (173 662 markers) for which only 12% of the dogs had both phenotypes and genotypes available. Our analysis revealed associations towards the growth differentiation factor 9 gene (GDF9), which is known to regulate oocyte maturation. The trait heritability was estimated at 43.1%, from which approximately 15% was accountable by the GDF9 locus alone. Therefore, markers flanking GDF9 explained approximately 6.5% of the variance in ALS. Analysis of WGSs revealed two missense substitutions in GDF9, one of which (g.11:21147009G>A) affected a highly conserved nucleotide in vertebrates. The derived allele A was validated in 111 dogs and shown to be associated with decreased ALS (−0.75 ± 0.22 puppies per litter). The variant was further predicted to cause a proline to serine substitution. The affected residue was immediately followed by a six-residue deletion that is fixed in the canine species but absent in non-canids. We further confirmed that the deletion is prevalent in the Canidae family by sequencing three species of wild canids. Since canids uniquely ovulate oocytes at the prophase stage of the first meiotic division, requiring maturation in the oviduct, we conjecture that the amino acid substitution and the six-residue deletion of GDF9 may serve as a model for insights into the dynamics of oocyte maturation in canids.     

Associations of polymorphisms in the promoter I of bovine acetyl‐CoA carboxylase‐α gene with beef fatty acid composition

The objectives of this study were to identify single nucleotide polymorphisms (SNPs) in the promoter I (PI) region of the bovine acetyl‐CoA carboxylase‐α (ACACA) gene and to evaluate the extent to which they were associated with lipid‐related traits. Eight novel SNPs were identified, which were AJ276223:g.2064T>A (SNP1), g.2155C>T (SNP2), g.2203G>T (SNP3), g.2268T>C (SNP4), g.2274G>A (SNP5), g.2340A>G (SNP6), g.2350T>C (SNP7) and g.2370A>G (SNP8). Complete linkage disequilibrium was observed among SNP1, 2, 4, 5, 6 and 8. Phenotypic data were collected from 573 cross‐bred steers with six sire breeds, including Hereford, Angus, Brangus, Beefmaster, Bonsmara and Romosinuano. The genotypes of SNP1/2/4/5/6/8 were significantly associated with adjusted backfat thickness. The genotypes of SNP3 were significantly associated with triacylglycerol (TAG) content and fatty acid composition of longissimus dorsi muscle (LM) in Brangus‐, Romosinuano‐ and Bonsmara‐sired cattle. Cattle with g.2203GG genotype had greater concentrations of TAG, total lipid, total saturated fatty acid and total monounsaturated fatty acid than did cattle with g.2203GT genotype. The genotypes of SNP7 were significantly associated with fatty acid composition of LM. Cattle with genotype g.2350TC had greater amounts of several fatty acids in LM than did cattle with genotype g.2350CC. Our results suggested that the SNPs in the PI region of ACACA gene are associated with variations in the fatty acid contents in LM.     

Distribution and linkage disequilibrium analysis of polymorphisms of <Emphasis Type="Italic">GH1</Emphasis> gene in different populations of pigs associated with body size

Growth hormone (GH) has been considered as a candidate gene for growth and body size in pigs. In this study, polymorphisms of the GH1 gene were evaluated for associations with body size traits in 190 pig individuals. Seventeen single-nucleotide polymorphisms (SNPs) were identified in GH1 gene of the large pig breeds and miniature pig breeds using direct sequencing and genotyped by allele-specific PCR approach. Notably, six (g.237A>G, g.283T>C, g.309A>G, g.318A>G, g.540A>G and g.544A>G) of them were significantly associated with body size, of which three loci (g.283T>C, g.309A>G, g.318A>G) located in the signal-peptide coding region of GH1 gene compose a CGG haplotype for large pigs and TAA haplotype for miniature pigs (P< 0.001), two loci (g.540A>G and g.544A>G) located in the second intron of GH1 gene compose a GG haplotype for large pigs and AA haplotype for miniature pigs (P< 0.001). Our results demonstrate that these SNPs in GH1 gene are associated with the body size of pigs providing genetic basis for pig breeding with the improved economic benefits.     

Association between METTL3 gene polymorphisms and neuroblastoma susceptibility: A nine‐centre case‐control study

Neuroblastoma ranks as the most commonly seen and deadly solid tumour in infancy. The aberrant activity of m⁶A‐RNA methyltransferase METTL3 is involved in human cancers. Therefore, functional genetic variants in the METTL3 gene may contribute to neuroblastoma risk. In the current nine‐centre case‐control study, we aimed to analyse the association between the METTL3 gene single nucleotide polymorphisms (SNPs) and neuroblastoma susceptibility. We genotyped four METTL3 gene SNPs (rs1061026 T>G, rs1061027 C>A, rs1139130 A>G, and rs1263801 G>C) in 968 neuroblastoma patients and 1814 controls in China. We found significant associations between these SNPs and neuroblastoma risk in neither single‐locus nor combined analyses. Interestingly, in the stratified analysis, we observed a significant risk association with rs1061027 AA in subgroups of children ≤ 18 months of age (adjusted OR = 1.87, 95% CI = 1.03‐3.41, P = .040) and females (adjusted OR = 1.86, 95% CI = 1.07‐3.24, P = .028). Overall, we identified a significant association between METTL3 gene rs1061027 C>A polymorphism and neuroblastoma risk in children ≤18 months of age and females. Our findings provide novel insights into the genetic determinants of neuroblastoma.     

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.