Novel SNPs of the Bovine PRLR Gene Associated with Milk Production Traits

The single nucleotide polymorphisms (SNPs) within exon 10 of the prolactin receptor gene (PRLR) were detected in Chinese Holstein cows using polymerase chain reaction-single-strand conformation polymorphism (PCR-SSCP) and DNA sequencing methods, and their genetic effects on milk production traits were evaluated in this study. Two newly detected SNPs (g.9206G→A and g.9681C→T) caused amino acid variations E378K and A536V, respectively, which were then preliminarily predicted at the topological level. Statistical results indicated that the two SNPs were significantly associated with milk yields, and cows with the combined genotype GGCC showed superior milk performance. A putative phosphorylation site was identified at residue 378K ([ST]-×-[RK]), which offers a partial explanation for the associations. These results suggest that the two novel SNPs within exon 10 of the PRLR gene associated with milk production traits are useful genetic markers in a selection program for Holstein dairy cattle.     

Genome wide CNV analysis reveals additional variants associated with milk production traits in Holsteins

Background

Milk production is an economically important sector of global agriculture. Much attention has been paid to the identification of quantitative trait loci (QTL) associated with milk, fat, and protein yield and the genetic and molecular mechanisms underlying them. Copy number variation (CNV) is an emerging class of variants which may be associated with complex traits.

Results

In this study, we performed a genome-wide association between CNVs and milk production traits in 26,362 Holstein bulls and cows. A total of 99 candidate CNVs were identified using Illumina BovineSNP50 array data, and association tests for each production trait were performed using a linear regression analysis with PCA correlation. A total of 34 CNVs on 22 chromosomes were significantly associated with at least one milk production trait after false discovery rate (FDR) correction. Some of those CNVs were located within or near known QTL for milk production traits. We further investigated the relationship between associated CNVs with neighboring SNPs. For all 82 combinations of traits and CNVs (less than 400 kb in length), we found 17 cases where CNVs directly overlapped with tag SNPs and 40 cases where CNVs were adjacent to tag SNPs. In 5 cases, CNVs located were in strong linkage disequilibrium with tag SNPs, either within or adjacent to the same haplotype block. There were an additional 20 cases where CNVs did not have a significant association with SNPs, suggesting that the effects of those CNVs were probably not captured by tag SNPs.

Conclusion

We conclude that combining CNV with SNP analyses reveals more genetic variations underlying milk production traits than those revealed by SNPs alone.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-683) contains supplementary material, which is available to authorized users.     

Genome wide association studies for milk production traits in Chinese Holstein population

Genome-wide association studies (GWAS) based on high throughput SNP genotyping technologies open a broad avenue for exploring genes associated with milk production traits in dairy cattle. Motivated by pinpointing novel quantitative trait nucleotide (QTN) across Bos Taurus genome, the present study is to perform GWAS to identify genes affecting milk production traits using current state-of-the-art SNP genotyping technology, i.e., the Illumina BovineSNP50 BeadChip. In the analyses, the five most commonly evaluated milk production traits are involved, including milk yield (MY), milk fat yield (FY), milk protein yield (PY), milk fat percentage (FP) and milk protein percentage (PP). Estimated breeding values (EBVs) of 2,093 daughters from 14 paternal half-sib families are considered as phenotypes within the framework of a daughter design. Association tests between each trait and the 54K SNPs are achieved via two different analysis approaches, a paternal transmission disequilibrium test (TDT)-based approach (L1-TDT) and a mixed model based regression analysis (MMRA). In total, 105 SNPs were detected to be significantly associated genome-wise with one or multiple milk production traits. Of the 105 SNPs, 38 were commonly detected by both methods, while four and 63 were solely detected by L1-TDT and MMRA, respectively. The majority (86 out of 105) of the significant SNPs is located within the reported QTL regions and some are within or close to the reported candidate genes. In particular, two SNPs, ARS-BFGL-NGS-4939 and BFGL-NGS-118998, are located close to the DGAT1 gene (160bp apart) and within the GHR gene, respectively. Our findings herein not only provide confirmatory evidences for previously findings, but also explore a suite of novel SNPs associated with milk production traits, and thus form a solid basis for eventually unraveling the causal mutations for milk production traits in dairy cattle.     

Association between polymorphisms in the SLC27A1 gene and milk production traits in Chinese Holstein cattle

We assessed SLC27A1, a candidate gene for milk production traits in Chinese Holstein cattle. DNA was extracted from the blood of 48 top Chinese Holstein Cattle selected according to phenotypic character and mixed into DNA pool for SNP detection. We tested blood samples of these cattle for SNPs in exon 3 and the 3'-flanking region of the SLC27A1 gene by using polymerase chain reaction-single-stranded conformation polymorphism (PCR-SSCP) and DNA sequencing. We found 2 polymorphic sites: 112T>C, a synonymous mutation, in exon 3 (SNP(1)), and 64G>A in the 3'-UTR (SNP(2)). We also determined the genotypes of 330 Chinese Holstein cattle by using PCR-restriction fragment length polymorphism (RFLP). We found 3 genotypes each at SNP(1) (TT, TC, and CC) and SNP(2) (GG, GA, and AA). The association among the different genotypes at these 2 sites and milk production traits was analyzed using a least-squares procedure. The results showed that cows with genotype CC had higher milk yields than those with genotype TC (0.01 < p < 0.05). No significant difference was detected among the 3 SNP(2) genotypes in terms of milk production traits. Our results provide evidence that the C allele have potential effects on milk yield trait.     

Three novel SNPs of the bovine Tf gene in Chinese native cattle and their associations with milk production traits

Transferrin (Tf) is a β-globulin protein that transports iron ions in mammalian cells. It contributes to innate immunity to microbial pathogens, primarily by limiting microbial access to iron. Thus, polymorphisms present in bovine Tf could potentially underlie inherited differences in mastitis resistance and milk production traits. We detected three novel single-nucleotide polymorphisms of the Tf gene in Chinese native cattle by screening for genetic variation of Tf in 751 individuals of three Chinese cattle breeds, namely China Holstein, Luxi Yellow and Bohai Black, using PCR-RFLP and DNA sequencing techniques. The three new SNPs, g.-1748G>A ss250608649, g.13942T>C ss250608650, and g.14037A>G ss250608651, had allele frequencies of 85.9, 86.3 and 92.5%, 64.5, 73.3 and 65.0%, and 67.6, 73.7 and 60.0%, respectively. SNP g.-1748G>A was located in the 5' flanking region of Tf. SNP g.14037A>G was located in intron 8 of Tf. SNP g.13942T>C, located in exon 8 of Tf, was a synonymous mutation (TTA > CTA), encoding a leucine (326 aa) in the Tf protein. Associations of the Tf SNPs with milk traits were also analyzed. Significant (P < 0.05) relationships among the Tf polymorphisms, somatic cell scores (SCS), and milk productive traits were observed. Cows with genotypes TT (g.13942T>C), GG (g.-1748G>A) and AG (g.14037A>G) had a lower SCS and higher protein levels and 305-day milk yield. Nineteen combinations of different haplotypes from the three SNPs were identified in Chinese Holstein cattle. The haplotype combination ATA/GCA, GCA/GCA and GCG/ GTA was dominant in cows with a lower SCS, a higher protein level and a higher 305-day milk yield, respectively. Moreover, the gene expression level of Tf was higher in mastitis-affected mammary tissues than in normal mammary tissues. These results suggest that the Tf gene affects milk production, as well as mastitis-resistance traits, in Chinese Holsteins.     

A Post-GWAS Replication Study Confirming the PTK2 Gene Associated with Milk Production Traits in Chinese Holstein

Our initial genome-wide association study (GWAS) demonstrated that two SNPs (ARS-BFGL-NGS-33248, UA-IFASA-9288) within the protein tyrosine kinase 2 (PTK2) gene were significantly associated with milk production traits in Chinese Holstein dairy cattle. To further validate if the statistical evidence provided in GWAS were true-positive findings, a replication study was performed herein through genotype-phenotype associations. The two tested SNPs were found to show significant associations with milk production traits, which confirmed the associations observed in the original study. Specifically, SNPs lying in the PTK2 gene were also detected by sequencing 14 unrelated sires in Chinese Holsteins and a total of thirty-three novel SNPs were identified. Thirteen out of these identified SNPs were genotyped and tested for association with milk production traits in an independent resource population. After Bonferroni correction for multiple testing, twelve SNPs were statistically significant for more than two milk production traits. Analyses of pairwise D’ measures of linkage disequilibrium (LD) between all SNPs were also explored. Two haplotype blocks were inferred and the association study at haplotype level revealed similar effects on milk production traits. In addition, the RNA expression analyses revealed that a non-synonymous coding SNP (g.4061098T>G) was involved in the regulation of gene expression. Thus the findings presented here provide strong evidence for associations of PTK2 variants with dairy production traits and may be applied in Chinese Holstein breeding program.     

Associations between newly discovered polymorphisms in the Bos taurusgrowth hormone receptor gene and performance traits in Holstein–Friesian dairy cattle

Variations in the growth hormone receptor (GHR) gene sequence are associated with performance traits in cattle. For example, the single nucleotide polymorphism (SNP) F279Y in transmembrane exon 8 has a strong association with milk yield. In this study, 32 previously unreported, putative novel SNPs (31 in the 5′ non‐coding region) were identified by resequencing ～19 kb of the GHR gene in genomic DNA from 22 cattle of multiple breeds. A population of 848 Holstein–Friesian AI sires was subsequently genotyped for the 32 putative novel SNPs and seven published SNPs (including F279Y, one in exon 1A promoter and five in exon 10). Associations between each segregating SNP and genetic merit for performance were quantified in the 848 Holstein–Friesians using weighted animal linear mixed models. Six of the published SNPs and seven of the novel SNPs were associated with at least one of the traits – milk yield, fat yield, protein yield, fat percentage, protein percentage, somatic cell score, calving interval, survival and growth and size traits. Even when the allelic substitution effect (P < 0.001) of F279Y was accounted for, the allelic substitution effect of one of the novel SNPs (GHR4.2) in the 5′ non‐coding region of GHR was associated with a lactation milk yield of 37.46 kg (P < 0.001). GHR4.2 and F279Y were not in linkage disequilibrium (r² = 0.00, D’ = 0.04) in the 848 Holstein–Friesians, indicating that their association with milk yield was independent.     

Genome Wide Analysis of Fertility and Production Traits in Italian Holstein Cattle

A genome wide scan was performed on a total of 2093 Italian Holstein proven bulls genotyped with 50K single nucleotide polymorphisms (SNPs), with the objective of identifying loci associated with fertility related traits and to test their effects on milk production traits. The analysis was carried out using estimated breeding values for the aggregate fertility index and for each trait contributing to the index: angularity, calving interval, non-return rate at 56 days, days to first service, and 305 day first parity lactation. In addition, two production traits not included in the aggregate fertility index were analysed: fat yield and protein yield. Analyses were carried out using all SNPs treated separately, further the most significant marker on BTA14 associated to milk quality located in the DGAT1 region was treated as fixed effect. Genome wide association analysis identified 61 significant SNPs and 75 significant marker-trait associations. Eight additional SNP associations were detected when SNP located near DGAT1 was included as a fixed effect. As there were no obvious common SNPs between the traits analyzed independently in this study, a network analysis was carried out to identify unforeseen relationships that may link production and fertility traits.     

Identification of SNPs in Interferon Gamma,Interleukin-22, and Their Receptors and Associations with Health and Production-Related Traits in Canadian Holstein Bulls

Genetic variants in a number of immunoregulatory genes have been previously associated with health and production traits in dairy cattle. Therefore, in the following study, the genes coding interferon gamma (IFNG), IFNG receptor 1 and 2 domains, interleukin-22 (IL22), and IL22 receptor alpha 1, were investigated for single nucleotide polymorphisms (SNPs) in Holstein bulls. These SNPs, along with SNPs previously identified in IL10, IL10 receptor, and transforming growth factor beta 1 (TGFB1) genes, were evaluated for statistical associations to estimated breeding values for milk somatic cell score (SCS), a trait highly correlated to mastitis incidence, and various production-related traits, including milk yield, protein yield, fat yield, and lactation persistency. While no significant associations were found between these SNPs and SCS, SNPs in IL10 receptor beta subunit showed a significant effect on protein yield and lactation persistency. While there is evidence that IL10 plays an important role during lactation, it is also likely that the effects of SNPs in IL10 receptor beta subunit on protein yield and lactation persistency are due to linkage disequilibrium with a neighboring QTL.     

Identification of SNPs in interferon gamma, interleukin-22, and their receptors and associations with health and production-related traits in Canadian Holstein bulls

Genetic variants in a number of immunoregulatory genes have been previously associated with health and production traits in dairy cattle. Therefore, in the following study, the genes coding interferon gamma (IFNG), IFNG receptor 1 and 2 domains, interleukin-22 (IL22), and IL22 receptor alpha 1, were investigated for single nucleotide polymorphisms (SNPs) in Holstein bulls. These SNPs, along with SNPs previously identified in IL10, IL10 receptor, and transforming growth factor beta 1 (TGFB1) genes, were evaluated for statistical associations to estimated breeding values for milk somatic cell score (SCS), a trait highly correlated to mastitis incidence, and various production-related traits, including milk yield, protein yield, fat yield, and lactation persistency. While no significant associations were found between these SNPs and SCS, SNPs in IL10 receptor beta subunit showed a significant effect on protein yield and lactation persistency. While there is evidence that IL10 plays an important role during lactation, it is also likely that the effects of SNPs in IL10 receptor beta subunit on protein yield and lactation persistency are due to linkage disequilibrium with a neighboring QTL.     

中国荷斯坦牛POU1F1基因与PRL基因的多态性及其聚合效应对产奶性状的影响

文章采用DNA测序、PCR-RFLP和CRS-PCR技术对979头中国荷斯坦牛POU1F1基因与PRL基因进行研究,发现了3个新SNPs,分别是POU1F1基因第二外显子G1178C、PRL基因5侧翼区A906G和A1134G。采用SAS统计软件GLM程序,利用最小二乘法拟合线性模型,分析基因多态性与产奶性状的关系。结果表明:POU1F1基因1178位点GC基因型在产奶量、乳蛋白量、乳脂量方面均为优良基因型。PRL基因5侧翼区906位点AG基因型在产奶量方面为优良基因型,1134位点不同基因型产奶性状差异不显著。对PRL基因5侧翼区的906位点和POU1F1基因的1178位点进行基因互作分析,结果在乳脂率、乳蛋白率、产奶量、乳蛋白量和乳脂量方面各基因型组合之间均未观察到显著差异,说明基因聚合效应并不是单基因效应的简单相加,基因聚合效应在分子育种中具有更重要的意义。     

Detection of loci affecting milk production and health traits in an elite US Holstein population using microsatellite markers

Quantitative trait loci (QTL) affecting health and milk production traits were studied in seven large half-sib US Holstein families by using the granddaughter design. Genotyping for 16 markers was completed and marker allele differences within and pooled-across families were analysed. Potential QTL were identified for somatic cell score (SCS), fat yield, fat percentage, protein yield and protein percentage. Three markers (BM203, BM4505 and BM2078) were associated with significant effects for different traits and, after further analysis, may be useful in marker-assisted selection in specific families. Comparisons between these data and previously identified QTL support the location of a QTL for milk yield and protein yield on chromosome 21.     

Genetic effects of stearoyl-coenzyme A desaturase (SCD) polymorphism on milk production traits in the Chinese dairy population

Stearoyl-CoA desaturase (SCD) is a multifunctional complex enzyme important in the cellular biosynthesis of fatty acids. The present study was to investigate the association of the SCD gene with milk production traits in dairy cattle. Two single nucleotide polymorphisms (SNPs) (g.6926A>G and g.8646A>G) in introns 3 and 4, and three SNPs (g.10153A>G, g.10213T>C and g.10329C>T) in exon 5 were identified with pooled DNA sequencing and genotyped using matrix-assisted laser desorption/ionization time of flight mass spectrometry assay in 752 Chinese Holstein cows. Polymorphism g.10329C>T was predicted to result in an amino acid replacement from alanine to valine in the SCD protein. With a mixed animal model, the significant associations of the five SNPs with 305-day milk, fat and protein yields and protein percentage were determined. We further demonstrated cows with heterozygous genotypes (A/G or C/T) had highest 305 day milk yield, fat yield, protein yield and lowest protein percentage. Heterozygous cows with genotype AG at the g.6926A>G locus showed the greatest milk yield (P < 0.0001), fat yield (P < 0.0001) and protein yield (P < 0.0001) among other heterozygous genotypes at any of the loci. Dominance effects of all identified SNPs on milk, fat and protein yields and protein percentage were significant. Moreover, significant allele substitution effects at g.6926A>G locus on milk yield and at g.10213T>C on protein yield were observed. Five-locus haplotypes and strong linkage disequilibrium (D' > 0.9) between the five SNPs were also observed. The results suggest that identified polymorphisms could be potential genetic markers to improve the production performance of Chinese Holstein.     

Polymorphisms of mitochondrial ATPase 8/6 genes and association with milk production traits in holstein cows

The maternal effect has been widely proposed to affect the production traits in domestic animals. However, the sequence polymorphisms of mitochondrial DNA (mtDNA) and association with milk production traits in Holstein cows have remained unclear. In this study, we investigated the single nucleotide polymorphisms (SNPs) of mtDNA ATPase 8/6 genes and association with four milk production traits of interest in 303 Holstein cows. A total of 18 SNPs were detected among the 842?bp fragment of ATPase 8/6 genes, which determined six haplotypes of B. taurus (H1-H4) and B. indicus (H5-H6). The mixed model analysis revealed that there was significant association between haplotype and 305-day milk yield (MY). The highest MY was observed in haplotype H4. However, we did not detect statistically significant differences among haplotypes for the traits of milk fat (MF), milk protein (MP), and somatic cell count (SC). The overall haplotype diversity and nucleotide diversity of ATPase 8/6 genes were 0.563?±?0.030 and 0.00609?±?0.00043, respectively. The results suggested that mitochondrial ATPase 8/6 genes could be potentially used as molecular marker to genetically improve milk production in Holstein cows.     

Revisiting the quantitative trait loci for milk production traits on BTA6

A parallel association study was performed in two independent cattle populations based on 41 validated, targeted single nucleotide polymorphisms (SNPs) and four microsatellite markers to re-evaluate the multiple quantitative trait loci (QTL) architecture for milk performance on bovine chromosome 6 (BTA6). Two distinct QTL located in the vicinity of the middle region of BTA6, but differing unambiguously regarding their effects on milk composition and yield traits were validated in the German Holstein population. A highly significant association of the protein variant ABCG2 p.Tyr581Ser with milk composition traits reconfirmed the causative molecular relevance of the ABCG2 gene in QTL region 1, whereas in QTL region 2, significant and tentative associations between gene variants RW070 and RW023 (located in the promoter region and exon 9 of the PPARGC1A gene for milk yield traits) were detected. For the German Fleckvieh population, only RW023 showed a tentative association with milk yield traits, whereas those loci with significant effects in German Holsteins (ABCG2 p.Tyr581Ser, RW070) showed fixed alleles. Even though our new data highlight two variants in the PPARGC1A gene (RW023, RW070) in QTL region 2, based on the results of our study, currently no unequivocal conclusion about the causal background of this QTL affecting milk yield traits can be drawn. Notably, the German Holstein and Fleckvieh populations, known for their divergent degree of dairy type, differ substantially in their allele frequencies for the growth-associated NCAPG p.Ile442Met locus.     

Nucleotide sequence polymorphism within exon 4 of the bovine prolactin gene and its associations with milk performance traits

Prolactin plays an important regulatory function in mammary gland development, milk secretion, and expression of milk protein genes. Hence the PRL gene is a potential quantitative trait locus and genetic marker of production traits in dairy cattle. We analysed the sequence of the PRL gene to investigate whether mutations in this sequence might be responsible for quantitative variations in milk yield and composition. Using SSCP and direct sequencing, we detected six single-nucleotide polymorphisms within a 294-bp prolactin gene fragment involving exon 4. All detected mutations were silent with respect to the amino acid sequence of the protein. PCR-RFLP genotyping of SNP 8398 R (RsaI) was used to assess allele frequencies in 186 Black-and-White cows (0.113 and 0.887 for A and G, respectively) and in 138 Jersey cows (0.706 and 0.294 for A and G, respectively). Black-and-White cows with genotype AG showed the highest milk yield, while cows with genotype GG showed the highest fat content.     

Polymorphism in the exon 4 of β-lactoglobulin variant B precursor gene and its association with milk traits and protein structure in Chinese Holstein

β-lactoglobulin (β-LG) is the major whey protein in the milk. In order to investigate the polymorphism of β-LG variant B precursor (β-LG B*: GenBank accession no. DQ489319) gene and its effects on the milk traits, the single-strand conformation polymorphism method (PCR-SSCP) were adopted to analyze polymorphism between 5229th and 5476th bp in the β-LG B* gene in Chinese Holstein. Four genotypes were found (AA, AB, AC and ABC) and 3 single nucleotide polymorphisms (SNPs) were detected (g.5239C>A, g.5240A>C, g.5305C>T and mix type g.5305C/T) in the exon 4 of β-LG B* gene. It was also found that the protein contents of AB, AC and ABC dairy cows were higher than AA (P < 0.05), and AC cows were the highest among them. Three SNPs (g.5239C>A, g.5240A>C and g.5305C>T) might affect the milk trait and all of them were high polymorphism (0.5 < PIC < 1.0). In further researches, the three SNPs also caused amino acid change (Asp>Glu, Thr>Pro and Ala>Val) respectively, and the spatial secondary and tertiary structure forecasting result also showed that single amino acid change influence protein spatial structure change in Chinese Holstein. Taken together, it is suggested that these SNPs change β-LG B* gene structure and expression. The polymorphism possibly holds the secret of milk protein and fat contents in the milk of Chinese Holstein.     

鸭生长激素基因内含子2、3多态性分析

根据鸭生长激素基因内含子2、3的序列设计5对引物,利用PCR-SSCP方法对北京鸭、西湖野鸭、金定鸭、山麻鸭、荆江鸭、绍兴鸭等6个鸭品种进行了单核苷酸多态性分析, 并检测其多态性。结果共发现8个突变位点, 其中内含子2有7个: 2593处C-T, 2770处G-A, 2813处T-A, 2829处C-A, 2894处C-T, 2896处T-C,3100处C-G; 内含子3有1个: 3270处A-G。统计结果显示, 这8个变异位点的基因型频率分布与品种有关, 在这些基因座的变异水平上, 北京鸭和绍兴鸭表现出了相当的品种保守性, 本研究所检测到的这些基因座可能与鸭的生产性能有关。