Association of UDP-galactose-4-epimerase with milk protein concentration in the Chinese Holstein population

Article information

Asian-Australas J Anim Sci. 2020;33(11):1725-1731
Publication date (electronic) : 2020 February 25
doi : https://doi.org/10.5713/ajas.19.0549
1Shaanxi Key Laboratory of Molecular Biology for Agriculture, College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi 712100, China
2College of Animal Science and Technology, Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing 100193, China
3Mammalian NutriPhysioGenomics, Department of Animal Sciences and Division of Nutritional Sciences, University of Illinois, Urbana, IL 61801, USA
*Corresponding Authors: Shengli Zhang, Tel: +86-010-62733967, Fax: +86-010-62733967, E-mail: zhangslcau@cau.edu.cn. Juan J. Loor, Tel: +1-217-244-5957, Fax: +1-217-333-5044, E-mail: loor@illinois.edu
Received 2019 July 8; Revised 2019 September 29; Accepted 2020 February 14.

Abstract

Objective

An initial RNA-Sequencing study revealed that UDP-galactose-4-epimerase (GALE) was one of the most promising candidates for milk protein concentration in Chinese Holstein cattle. This enzyme catalyzes the interconversion of UDP-galactose and UDP-glucose, an important step in galactose catabolism. To further validate the genetic effect of GALE on milk protein traits, genetic variations were identified, and genotypes-phenotypes associations were performed.

Methods

The entire coding region and the 5′-regulatory region (5′-UTR) of GALE were re-sequenced using pooled DNA of 17 unrelated sires. Association studies for five milk production traits were performed using a mixed linear animal model with a population encompassing 1,027 Chinese Holstein cows.

Results

A total of three variants in GALE were identified, including two novel variants (g.2114 A>G and g.2037 G>A) in the 5′-UTR and one previously reported variant (g.3836 G>C) in an intron. All three single nucleotide polymorphisms (SNPs) were associated with milk yield (p<0.0001), fat yield (p = 0.0006 to <0.0001), protein yield (p = 0.0232 to <0.0001) and protein percentage (p<0.0001), while no significant associations were detected between the SNPs and fat percentage. A strong linkage disequilibrium (D’ = 0.96 to 1.00) was observed among all three SNPs, and a 5 Kb haplotype block involving three main haplotypes with GAG, AGC, and AGG was formed. The results of haplotype association analyses were consistent with the results of single locus association analysis (p<0.0001). The phenotypic variance ratio above 3.00% was observed for milk protein yield that was explained by SNP-g.3836G >C.

Conclusion

Overall, our findings provided new insights into the polymorphic variations in bovine GALE gene and their associations with milk protein concentration. The data indicate their potential uses for marker-assisted breeding or genetic selection schemes.

INTRODUCTION

Milk proteins are important nutrients and milk protein concentration serves as valuable index for milk quality. Dairy industry concerns have driven increasing efforts to improve milk protein concentration [1]. With the development of genomics, bioinformatics and statistical genetics, a single gene or chromosome segments affecting important economic traits can be analyzed [2]. It is possible to improve milk protein concentration through marker assisted selection (MAS) or genomic selection schemes, the challenge, however, is to identify key genes or causal variations affecting milk protein traits [36]. Our previously published research has identified that UDP-galactose-4-epimerase (GALE) was a strong candidate gene for milk protein traits due to its large differential expression (Log2 fold-change = −0.74, q-value = 4.41E-03) in mammary tissues of cows with high and low milk protein percentage [7]. In addition, strong interactions were also observed between GALE and several other genes such as lactalbumin, alpha (LALBA), beta 1,4-galactosyltransferase, polypeptide 1 (B4GALT1), and UDP-glucose 6-dehydrogenase (UGDH) that play important roles in milk composition synthesis (Supplementary Figure S1) [810]. Therefore, based on the biological function and transcriptional effects on milk protein traits, the current study mined to screen the full-length coding regions of the GALE gene for single nucleotide polymorphisms (SNPs) and to evaluate the genetic effects of polymorphisms on milk production traits in a large Chinese Holstein population.

MATERIALS AND METHODS

Ethics statement

Animal handling and sample collection procedures were performed in accordance with protocols approved by the Institutional Animal Care and Use Committee (IACUC) at China Agricultural University (Permit Number: DK996).

Genetic sampling

A total of 1,027 Chinese Holstein cows and their 17 corresponding sires were considered as the study population. Cows were selected from 17 farms in the Beijing Sanyuan Lvhe Dairy Farm Center, where routine standard performance test, i.e. Dairy Herd Improvement system (DHI) has been implemented since 1999. The phenotype observations for all individuals for five milk production traits (305 d milk yield, 305 d protein yield, 305 d fat yield, average 305 d protein percentage and average 305 d fat percentage) were collected for subsequent analyses via the complete DHI data from the Chinese dairy cattle population.

Single nucleotide polymorphism identification and genotyping

Blood samples were collected from 1,027 cows via coccygeal vein to isolate genomic DNA using DP (318) Blood DNA kits (TianGen, Beijing, China). Genomic DNA were also isolated from frozen semen of 17 sires using standard phenol-chloroform procedures and were pooled with 50 ng/μL DNA of each individual to identify variants of the GALE gene. GALE gene is 4,591 bp in length located at BTA2, contains 11 exons and 10 introns encoding 348 amino acids. All exons and their adjacent intronic sequences were targeted for selective amplification by polymerase chain reaction (PCR). A total of 12 pairs of nucleotide primers (Supplementary Table S1) targeting the regions of interest were designed using Primer3 based on the genomic sequence of the bovine GALE gene referring to Bos_taurus_UMD_3.1 assembly (NCBI reference sequence: AC_000159.1).

The PCR amplification was performed in a total volume of 25 μL containing 50 to 100 ng of genomic DNA, 0.5 μL of each primer, 2.5 mM of dNTP mix, 2.5 μL of 10× PCR buffer, and 1 U of Taq DNA Polymerase (Takara Biotechnology Co., Ltd, Dalian, China). The PCR reaction conditions included a pre-denaturation at 95°C for 5 min, followed by 34 cycles of 94°C for 30 s, annealing from 46°C to 56°C for 30 s, 72°C for 30 s, and a final extension at 72°C for 10 min. The PCR products were directly sequenced using the ABI3730xl DNA analyzer (Applied Biosystems, Foster City, CA, USA), and the sequences aligned to the bovine reference sequence (UMD3.1.1) using BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi) to identify potential SNPs.

The details of novel SNPs that were identified in the present study were submitted to dbSNP (http://www.ncbi.nlm.nih.gov/SNP/) and are publicly available (accession numbers from ss1996900612 to ss1996900613). All identified SNPs for subsequent genotyping in the 1,027 Chinese Holstein cows were performed with matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI-TOF MS, Squenom MassARRAY, Bioyong Technologies Inc. Hong Kong, China) assay.

Linkage disequilibrium analysis

Haploview [11] was used to measure pairwise linkage disequilibrium (LD) for all identified SNPs within GALE. Briefly, missing genotypes were first imputed for each individual using the Beagle3.2 software program [12]. Subsequently, the LD blocks were generated with the subject genotype data using the LD coefficient (D’) [13]. A haplotype with a frequency >5% was considered as a distinguishable haplotype, while the haplotypes with relative frequency <5% were pooled into a single group. Haplotype blocks within these SNPs were used to test their associations with milk production traits.

Association analyses

A goodness-of-fit test (Chi-square) was applied to compare the numbers of expected and observed genotypes to test Hardy-Weinberg equilibrium for each identified SNP, 0.05 as the significant threshold value. Association analyses were conducted to estimate the effects of GALE variants on milk production traits based on both single SNP genotypes and the haplotype blocks. The effects of single SNP or haplotypes in GALE on the five milk production traits were analyzed with the mixed procedure of SAS9.3 software (SAS Institute Inc., Cary, NC, USA) using the following mixed linear animal model:

yijklmn=μ+Fi+YSj+Pk+b×M+Gl+αm+eijklmn

where, yijklmn was the phenotypic value of each trait for each cow (n = 1,027 for each trait); μ was the overall mean; Fi was the fixed effect of farm; YSj was the fixed effect of year-season; Pk was the fixed effect of parity; M was the covariate effect of calving month; b was the regression coefficient of M; Gl was the fixed effect corresponding to the genotype of polymorphisms or haplotype (genotypes of SNPs were modelled as 0–1–2, haplotypes were modelled as 0–5); αm was the random polygenic effect, distributed as N (0, Aσa2), with the additive genetic relationship matrix A and the additive genetic variance σa2; A-matrix was constructed by tracing the pedigree back to three generations of 2,312 involved individuals; and eijklmn was the random residual, distributed as N (0, Iσe2), with identity matrix I and residual error variance σe2.

For single SNP and haplotype analyses, the Bonferroni method was adopted to correct for multiple-testing according to the number of SNP loci or haplotype blocks. Associations were considered as significant if a raw p value <0.05/N, where N was the number of SNP loci or haplotype blocks tested in analyses. The additive (a), dominance (d), and allele substitution (α) effects were estimated using the equation from Falconer and Mackay [14], i.e. a = (AA–BB)⁄2, d = AB–(AA+BB)⁄2, and α = a+d(q–p), where AA and BB represented the two homozygous genotypes, AB was heterozygous genotype, and p and q were the allele frequencies of corresponding loci.

The effect of a SNP on a specific trait was measured as the proportion of phenotypic variance of the trait explained by the SNP. The proportion of variance explained by a SNP was calculated as 2p(1-p)α2/σp2, where p was the allele frequency of SNP, α was the average effect of gene substitution calculated based on the linear mixed model, and σp2 was the estimate of phenotypic variance using the complete DHI data of the Chinese dairy cattle population.

RESULTS

A total of three SNPs were discovered in the GALE gene of which two (g.2114A>G and g.2037G>A) in the 5′-UTR are novel (ss1996900612 and ss1996900613). The other SNP (g.3836G>C) previously reported was located in the intronic region (rs211659075) (Table 1). All three SNPs were in Hardy-Weinberg equilibrium (p>0.05, Table 2).

Information for the three identified single nucleotide polymorphisms in UDP-galactose-4-epimerase gene

Genotypic and allelic frequencies and Hardy-Weinberg equilibrium test of three single nucleotide polymorphisms of UDP-galactose-4-epimerase gene in Chinese Holstein cattle

The association results between the identified SNPs in GALE and five milk production traits are presented in Table 3. All three SNPs (g.2114A>G, g.2037G>A and g.3836G>C) were highly associated with milk yield (p<0.0001), fat yield (p = 0.0006 to <0.0001), protein yield (p = 0.0232 to <0.0001), and protein percentage (p<0.0001). No significant associations were observed between the SNPs and milk fat percentage. Greater than 1% of phenotypic variation accounted for by the three SNPs was detected in six significant SNP-trait pairs. Within these pairs, the pairs of g.3836G>C and milk yield, g.3836G>C and milk protein yield and g.3836G>C and milk protein percentage accounted for up to 2.61%, 3.00%, and 1.08% of phenotypic variation, respectively. In addition, significant additive effects, dominant effects and allele substitution effects were observed for the significant related traits (Table 4).

Associations of UDP-galactose-4-epimerase single nucleotide polymorphisms with milk production traits in Chinese Holstein cattle (LSM±SE)

Additive, dominant and allele substitution effects of the three single nucleotide polymorphisms associated with milk production traits of UDP-galactose-4-epimerase in Chinese Holstein

The specific LD results are showed in Supplementary Table S2–S3 and Supplementary Figure S2. Strong LD (D’ = 0.96 to 1.00) was observed between the three identified SNPs. A 5 Kb haplotype block was inferred (Figure 1), and three major haplotypes were constructed: GAG, AGC and AGG, with frequencies of 59.8%, 19.8%, and 15.0%, respectively (Table 5). A total of six genotypes of haplotypes, H1H1, H2H1, H2H2, H2H3, H3H1, and H3H3 (H1 = GAG, H2 = AGC, H3 = AGG) were formed. Haplotype association analysis revealed consistent results with single-locus analysis (Table 6). The haplotypes H1 and H3 had higher milk yield and protein yield than H2 haplotype (Table 6).

Figure 1

The haplotype blocks and pairwise linkage disequilibrium (LD) values (D’) for the UDP-galactose-4-epimerase (GALE) single nucleotide polymorphisms (SNPs). The values within boxes are pair-wise SNP correlation (D’), black boxes without numbers indicate complete LD (D’ = 1).

Main haplotypes and their frequencies observed in GALE gene

Haplotype associations of GALE single nucleotide polymorphisms with milk production traits in Chinese Holstein cattle (LSM±SE)

DISCUSSION

Polymorphisms located in the promoter of a gene may affect transcription by altering transcription factor binding sites or RNA stability [15], which indicated the importance of the two novel SNPs (g.2114A>G and g.2037G>A) identified in the 5′-UTR of GALE. The intronic SNP (g.3836G>C) may have a potential regulatory effect on gene expression, regulation, transcription and mRNA splicing, although it does not hold a sequence encoding a protein [1619]. The greater expression of GALE in mammary tissues of cows with high versus low milk protein percentage [7] agreed with such effect. To our knowledge, this was the first evidence showing significant associations of the GALE gene with milk protein traits in dairy cattle.

From a statistical standpoint, the single SNP association analysis was less powerful than multiple SNPs analysis due to the lack of simultaneous use of multiple SNPs information [20,21]. Thus, the haplotype-based association analysis was further performed to investigate the association of GALE variants with milk production traits in the present study. We observed that the three identified SNPs were associated with milk yield and milk protein traits, which was further confirmed by haplotype-based association analysis. Both single and haplotype association analyses demonstrated that the GALE gene was a promising candidate gene affecting milk yield and protein traits.

Protein GALE is UDP-galactose-4-epimerase, which catalyzes the interconversion of UDP-galactose and UDP-glucose in the final step of the Leloir pathway [22,23], and catalyzes the epimerization of UDP-N-acetylglucosamine to UDP-N-acetylgalactosamine [24,25]. GALE plays critical roles in dietary galactose metabolism, endogenous galactose production, and glycoprotein and glycolipid biosynthesis [26,27]. String interaction network (https://string-db.org/network/9606.ENSP00000363621) revealed that GALE protein interacts with LALBA, UDP-Gal: betaGlcNAc B4GALT1, and UGDH. Among this list, LALBA is the major component of milk protein and a subunit of lactose synthase. As one of the well-studied glycosyltransferases, B4GALT1 is responsible for the synthesis of complex-type N-liked oligosaccharides in many glycoproteins [8]. The association of polymorphisms of the B4GALT1 with milk production traits in Holstein cows has been reported in previously published research [9]. The UGDH gene was suggested to be associated with milk yield and milk composition [10]. Taken altogether, the influences of GALE on milk production and composition are likely due to the interaction of GALE with those known genes.

CONCLUSION

Results in the present study demonstrated the significant genetic effects of GALE on milk protein traits, which is in close agreement with our previous RNA-Seq study. Results also confirmed the phenotypes of milk protein were directly affected by GALE gene at genome level and transcriptional level. Due to the high phenotypic variance ratio, the SNP g.3836G>C in bovine GALE may be the most promising marker implicated in milk protein concentration in dairy cattle and has the capability to be used in MAS. Hence, results lay a preliminary foundation for further identifying the causal mutations affecting milk proteins in follow-up studies.

Supplementary Information

Notes

CONFLICT OF INTEREST

We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.

ACKNOWLEDGMENTS

This study was supported by the Young Talent Fund of University Association for Science and Technology in Shaanxi, China (20170202), the Key Research and Development Plan of Shaanxi Province, China (2018ZDCXL-NY-01-05), the Fundamental Research Funds for the Central Universities (2452018126), and Ph.D research start-up foundation of Northwest A&F University (2452016203).

References

1. Dallas DC, Murray NM, Gan JN. Proteolytic systems in milk: perspectives on the evolutionary function within the mammary gland and the infant. J Mammary Gland Biol Neoplasia 2015;20:133–47. https://doi.org/10.1007/s10911-015-9334-3 .
2. Carta A, Casu S, Salaris S. Invited review: Current state of genetic improvement in dairy sheep. J Dairy Sci 2009;92:5814–33. https://doi.org/10.3168/jds.2009-2479 .
3. Chen HY, Zhang Q, Yin CC, Wang CK, Gong WJ, Mei G. Detection of quantitative trait loci affecting milk production traits on bovine chromosome 6 in a Chinese Holstein population by the daughter design. J Dairy Sci 2006;89:782–90. https://doi.org/10.3168/jds.S0022-0302(06)72140-3 .
4. Gebreyesus G, Lund MS, Janss L, et al. Short communication: Multi-trait estimation of genetic parameters for milk protein composition in the Danish Holstein. J Dairy Sci 2016;99:2863–6. https://doi.org/10.3168/jds.2015-10501 .
5. Kolbehdari D, Wang Z, Grant JR, et al. A whole genome scan to map QTL for milk production traits and somatic cell score in Canadian Holstein bulls. J Anim Breed Genet 2009;126:216–27. https://doi.org/10.1111/j.1439-0388.2008.00793.x .
6. Sanchez MP, Ferrand M, Gele M, et al. Short communication: Genetic parameters for milk protein composition predicted using mid-infrared spectroscopy in the French Montbeliarde, Normande, and Holstein dairy cattle breeds. J Dairy Sci 2017;100:6371–5. https://doi.org/10.3168/jds.2017-12663 .
7. Li C, Cai W, Zhou C, et al. RNA-Seq reveals 10 novel promising candidate genes affecting milk protein concentration in the Chinese Holstein population. Sci Rep 2016;6:26813. https://doi.org/10.1038/srep26813 .
8. Amado M, Almeida R, Schwientek T, Clausen H. Identification and characterization of large galactosyltransferase gene families: galactosyltransferases for all functions. Biochim Biophys Acta Gen Subj 1999;1473:35–53. https://doi.org/10.1016/S0304-4165(99)00168-3 .
9. Shahbazkia HR, Aminlari M, Cravador A. Association of polymorphism of the beta(1, 4)-galactosyltransferase-I gene with milk production traits in Holsteins. Mol Biol Rep 2012;39:6715–21. https://doi.org/10.1007/s11033-012-1495-1 .
10. Xu Q, Mei G, Sun DX, et al. Detection of genetic association and functional polymorphisms of UGDH affecting milk production trait in Chinese Holstein cattle. BMC Genomics 2012;13:590. https://doi.org/10.1186/1471-2164-13-590 .
11. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 2005;21:263–5. https://doi.org/10.1093/bioinformatics/bth457 .
12. Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet 2007;81:1084–97. https://doi.org/10.1086/521987 .
13. Gabriel SB, Schaffner SF, Nguyen H, et al. The structure of haplotype blocks in the human genome. Science 2002;296:2225–9. https://doi.org/10.1126/science.1069424 .
14. Falconer DS, Mackay TFC. Introduction to quantitative genetics 4th ednth ed. New York, USA: Longman Scientific and Technical; 1996.
15. Huang YZ, Li MX, Wang J, et al. A 5′-regulatory region and two coding region polymorphisms modulate promoter activity and gene expression of the growth suppressor gene ZBED6 in cattle. Plos One 2013;8:e79744. https://doi.org/10.1371/journal.pone.0079744 .
16. Nott A, Muslin SH, Moore MJ. A quantitative analysis of intron effects on mammalian gene expression. RNA 2003;9:607–17. https://doi.org/10.1261/rna.5250403 .
17. Park HJ, Lee S, Ju E, Jones JA, Choi I. Alternative transcription of sodium/bicarbonate transporter SLC4A7 gene enhanced by single nucleotide polymorphisms. Physiol Genomics 2017;49:167–76. https://doi.org/10.1152/physiolgenomics.00112.2016 .
18. Sagne C, Marcel V, Amadou A, Hainaut P, Olivier M, Hall J. A meta-analysis of cancer risk associated with the TP53 intron 3 duplication polymorphism (rs17878362): geographic and tumor-specific effects. Cell Death Dis 2013;4:e492. https://doi.org/10.1038/cddis.2013.24 .
19. Visser M, Palstra RJ, Kayser M. Human skin color is influenced by an intergenic DNA polymorphism regulating transcription of the nearby BNC2 pigmentation gene. Hum Mol Genet 2014;23:5750–62. https://doi.org/10.1093/hmg/ddu289 .
20. Akey J, Jin L, Xiong MM. Haplotypes vs single marker linkage disequilibrium tests: what do we gain? Eur J Hum Genet 2001;9:291–300. https://doi.org/10.1038/sj.ejhg.5200619 .
21. Martin ER, Lai EH, Gilbert JR, et al. SNPing away at complex diseases: Analysis of single-nucleotide polymorphisms around APOE in Alzheimer disease. Am J Hum Genet 2000;67:383–94. https://doi.org/10.1086/303003 .
22. Daenzer JM, Sanders RD, Hang D, Fridovich-Keil JL. UDP-galactose 4′-epimerase activities toward UDP-Gal and UDP-GalNAc play different roles in the development of Drosophila melanogaster . PLoS Genet 2012;8:e1002721. https://doi.org/10.1371/journal.pgen.1002721 .
23. Schulz JM, Ross KL, Malmstrom K, Krieger M, Fridovich-Keil JL. Mediators of galactose sensitivity in UDP-galactose 4′-epimerase-impaired mammalian cells. J Biol Chem 2005;280:13493–502. https://doi.org/10.1074/jbc.M414045200 .
24. Roper JR, Guther MLS, Milne KG, Ferguson MAJ. Galactose metabolism is essential for the African sleeping sickness parasite Trypanosoma brucei . Proc Natl Acad Sci USA 2002;99:5884–9. https://doi.org/10.1073/pnas.092669999 .
25. Seo A, Gulsuner S, Pierce S, et al. Inherited thrombocytopenia associated with mutation of UDP-Galactose-4-Epimerase (GALE). Hum Mol Genet 2019;28:133–42. https://doi.org/10.1093/hmg/ddy334 .
26. Sanders RD, Sefton JMI, Moberg KH, Fridovich-Keil JL. UDP-galactose 4′ epimerase (GALE) is essential for development of Drosophila melanogaster . Dis Model Mech 2010;3:628–38. https://doi.org/10.1242/dmm.005058 .
27. Song HB, He M, Cai ZP, et al. UDP-glucose 4-epimerase and-1,4-galactosyltransferase from the oyster Magallana gigas as valuable biocatalysts for the production of galactosylated products. Int J Mol Sci 2018;19:1600. https://doi.org/10.3390/ijms19061600 .

Article information Continued

Figure 1

The haplotype blocks and pairwise linkage disequilibrium (LD) values (D’) for the UDP-galactose-4-epimerase (GALE) single nucleotide polymorphisms (SNPs). The values within boxes are pair-wise SNP correlation (D’), black boxes without numbers indicate complete LD (D’ = 1).

Table 1

Information for the three identified single nucleotide polymorphisms in UDP-galactose-4-epimerase gene

CHR RefSNP SNP locus Alleles Location Position Gene
2 ss1996900612 g.2114A>G A/G 5′-UTR 129705167 GALE
2 ss1996900613 g.2037G>A G/A 5′-UTR 129705244 GALE
2 rs211659075 g.3836G>C G/C Intron-9 129711117 GALE

SNPs, single nucleotide polymorphisms; GALE, UDP-galactose-4-epimerase.

Table 2

Genotypic and allelic frequencies and Hardy-Weinberg equilibrium test of three single nucleotide polymorphisms of UDP-galactose-4-epimerase gene in Chinese Holstein cattle

Position Locus Genotypes N Frequency Allele Frequency Hardy-Weinberg equilibrium χ2 test
5′flanking region ss1996900612 AG 457 0.454 A 0.350 p>0.05
g.2114A>G AA 124 0.123 G 0.650
GG 425 0.422
5′flanking region ss1996900613 AG 491 0.483 A 0.594 p>0.05
g.2037G>A AA 358 0.352 G 0.406
GG 167 0.164
Intron-9 rs211659075 CG 326 0.323 C 0.198 p>0.05
g.3836G>C CC 37 0.037 G 0.802
GG 645 0.640

Table 3

Associations of UDP-galactose-4-epimerase single nucleotide polymorphisms with milk production traits in Chinese Holstein cattle (LSM±SE)

Locus Genotypes Milk yield (kg) Fat yield (kg) Fat percentage (%) Protein yield (kg) Protein percentage (%)
ss 1996900612 AA(124) 10,507±88.95AB 371.88±3.75AB 3.595±0.036 331.02±2.73ab 3.171±0.012A
g.2114A>G AG(457) 10,402±61.90A 364.66±2.59A 3.585±0.025 327.83±1.89a 3.192±0.009A
GG(425) 10,667±63.55B 373.54±2.65B 3.546±0.026 332.53±1.93b 3.136±0.009B
p-value1) <0.0001 0.0006 0.1567 0.0232 <0.0001
Variance 0.05% 0.13% 0.22% 0.02% 0.16%
ss 1996900613 AA(358) 10,768±65.48A 375.79±2.74A 3.543±0.026 335.94±1.99A 3.141±0.009A
g.2037G>A AG(491) 10,326±61.30B 361.05±2.56B 3.572±0.025 324.63±1.87B 3.185±0.009B
GG(167) 10,404±81.06B 368.40±3.40AB 3.602±0.033 327.01±2.48B 3.169±0.011B
p-value1) <0.0001 <0.0001 0.1643 <0.0001 <0.0001
Variance 1.64% 0.14% 0.50% 1.03% 0.30%
rs211659075 CC(37) 9,597.66±141.44A 338.99±5.98A 3.605±0.057 301.50±4.36A 3.166±0.020AB
g.3836G>C CG(326) 10,359±66.80B 363.19±2.80B 3.596±0.027 326.87±2.04B 3.202±0.009A
GG(645) 10,581±59.48C 370.71±2.48C 3.559±0.024 330.91±1.81B 3.154±0.008B
p-value1) <0.0001 <0.0001 0.2294 <0.0001 <0.0001
Variance 2.61% 1.56% 0.08% 3.00% 1.08%

LSM, least square mean; SE, standard error; SNP, single nucleotide polymorphisms.

1)

p-value refers to the results of association analysis between each SNP and milk production traits.

Different letter (small letters, p<0.05; capital letters, p<0.01) superscripts (adjusted value after correction for multiple testing) indicate significant differences among the genotypes. Phenotypicvarianceratio=2p(1-p)α2/σp2.

Table 4

Additive, dominant and allele substitution effects of the three single nucleotide polymorphisms associated with milk production traits of UDP-galactose-4-epimerase in Chinese Holstein

Locus Genetic effect Milk yield Fat yield Fat percentage Protein yield Protein percentage
ss 1996900612 Additive 79.68±42.75 −0.83±1.81 0.0245±0.02 −0.75±1.32 0.0176±0.01**
g.2114A>G Dominant −184.83±56.79** −8.05±2.41** 0.0146±0.02 −3.95±1.76* 0.0384±0.01**
Substitution −24.38±51.78 1.58±2.19 0.0201±0.02 0.43±1.60 0.0061±0.01
ss 1996900613 Additive 182.15±39.80** 3.70±1.69* −0.0295±0.02 4.46±1.23** −0.0140±0.01**
g.2037G>A Dominant −259.98±52.84** −11.05±2.24** −0.0001±0.02 −6.84±1.63** 0.0304±0.01**
Substitution 133.27±43.35** 1.62±1.84 −0.0295±0.02 3.18±1.34* −0.0083±0.01
rs211659075 Additive −491.46±69.94** 15.86±2.96** 0.0229±0.03 −14.70±2.16** 0.0062±0.01
g.3836G>C Dominant 269.59±81.71** 8.34±3.47** 0.0140±0.03 10.66±2.53** 0.0420±0.01**
Substitution −654.07±111.31** −20.89±4.72** 0.0144±0.04 −21.13±3.44** −0.0191±0.02
*

Means differ at p<0.05.

**

Means differ at p<0.01.

Table 5

Main haplotypes and their frequencies observed in GALE gene

GALE haplotypes ss1996900612 A>G ss1996900613 G>A rs211659075 G>C Frequency (%)
GAG G A G 59.8
AGC A G C 19.8
AGG A G G 15.0

GALE, UDP-galactose-4-epimerase; SNPs, single nucleotide polymorphisms.

The Ref number of each SNP can be found in the haplotype Figure 1.

Table 6

Haplotype associations of GALE single nucleotide polymorphisms with milk production traits in Chinese Holstein cattle (LSM±SE)

GALE haplotypes Milk yield (kg) Fat yield (kg) Fat percentage (%) Protein yield (kg) Protein percentage (%)
H1H1(369) 10,779±66.50A 378.16±2.78AC 3.548±0.027 336.27±2.03A 3.140±0.009A
H2H1(243) 10,403±72.46B 362.51±3.04BD 3.558±0.029 328.14±2.21B 3.199±0.010B
H2H2(39) 9,891.26±138.99C 350.80±5.86B 3.611±0.056 311.33±4.27C 3.172±0.020AB
H2H3(61) 10,798±118.84AD 387.90±5.01A 3.628±0.047 341.15±3.65A 3.175±0.016AB
H3H1(181) 10,535±79.06BD 371.77±3.31CD 3.578±0.032 331.72±2.41AB 3.177±0.011B
H3H3(26) 10,986±173.95AD 369.53±7.35ABC 3.405±0.070 342.22±5.36AB 3.137±0.024AB
p-value1) <0.0001 <0.0001 0.0807 <0.0001 <0.0001

GALE, UDP-galactose-4-epimerase; LSM, least square mean; SE, standard error.

1)

p-value refers to the results of association analysis between each haplotype and milk production traits.

Different letter superscripts indicate significant differences among the haplotypes (p<0.01). H1, H2, and H3 represented the types of haplotypes, of these, H1 = GAG, H2 = AGC, H3 = AGG.