1
|
Olasege BS, van den Berg I, Haile-Mariam M, Ho PN, Yin Oh Z, Porto-Neto LR, Hayes BJ, Pryce JE, Fortes MRS. Dissecting loci that underpin the genetic correlations between production, fertility, and urea traits in Australian Holstein cattle. Anim Genet 2024; 55:540-558. [PMID: 38885945 DOI: 10.1111/age.13455] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2023] [Revised: 05/09/2024] [Accepted: 05/25/2024] [Indexed: 06/20/2024]
Abstract
Unfavorable genetic correlations between milk production, fertility, and urea traits have been reported. However, knowledge of the genomic regions associated with these unfavorable correlations is limited. Here, we used the correlation scan method to identify and investigate the regions driving or antagonizing the genetic correlations between production vs. fertility, urea vs. fertility, and urea vs. production traits. Driving regions produce an estimate of correlation that is in the same direction as the global correlation. Antagonizing regions produce an estimate in the opposite direction of the global estimates. Our dataset comprised 6567, 4700, and 12,658 Holstein cattle with records of production traits (milk yield, fat yield, and protein yield), fertility (calving interval) and urea traits (milk urea nitrogen and blood urea nitrogen predicted using milk-mid-infrared spectroscopy), respectively. Several regions across the genome drive the correlations between production, fertility, and urea traits. Antagonizing regions were confined to certain parts of the genome and the genes within these regions were mostly involved in preventing metabolic dysregulation, liver reprogramming, metabolism remodeling, and lipid homeostasis. The driving regions were enriched for QTL related to puberty, milk, and health-related traits. Antagonizing regions were mostly related to muscle development, metabolic body weight, and milk traits. In conclusion, we have identified genomic regions of potential importance for dairy cattle breeding. Future studies could investigate the antagonizing regions as potential genomic regions to break the unfavorable correlations and improve milk production as well as fertility and urea traits.
Collapse
Affiliation(s)
- Babatunde S Olasege
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia
- CSIRO Agriculture and Food, Saint Lucia, Queensland, Australia
| | - Irene van den Berg
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, Victoria, Australia
| | - Mekonnen Haile-Mariam
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, Victoria, Australia
- School of Applied Systems Biology, La Trobe University, Bundoora, Victoria, Australia
| | - Phuong N Ho
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, Victoria, Australia
| | - Zhen Yin Oh
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia
| | | | - Ben J Hayes
- The University of Queensland, Queensland Alliance for Agriculture and Food Innovation (QAAFI), Brisbane, Queensland, Australia
| | - Jennie E Pryce
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, Victoria, Australia
- School of Applied Systems Biology, La Trobe University, Bundoora, Victoria, Australia
| | - Marina R S Fortes
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia
- The University of Queensland, Queensland Alliance for Agriculture and Food Innovation (QAAFI), Brisbane, Queensland, Australia
| |
Collapse
|
2
|
Atashi H, Chen Y, Vanderick S, Hubin X, Gengler N. Single-step genome-wide association analyses for milk urea concentration in Walloon Holstein cows. J Dairy Sci 2024; 107:3020-3031. [PMID: 38056570 DOI: 10.3168/jds.2023-23902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Accepted: 11/07/2023] [Indexed: 12/08/2023]
Abstract
The aims of this study were to conduct a single-step genome-wide association to identify genomic regions associated with milk urea (MU) and to estimate genetic correlations between MU and milk yield (MY), milk composition (calcium content [CC], fat percentage [FP], protein percentage [PP], and casein percentage [CNP]), and cheese-making properties (CMP; coagulation time [CT], curd firmness after 30 min from rennet addition [a30], and titratable acidity [TA]). The used data have been collected from 2015 to 2020 on 78,073 first-parity (485,218 test-day records) and 48,766 second-parity (284,942 test-day records) Holstein cows distributed in 671 herds in the Walloon Region of Belgium. Data of 565,533 SNP located on 29 BTA of 6,617 animals (1,712 males) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic BLUP approach. The proportion of the total additive genetic variance explained by windows of 50 consecutive SNPs (with an average size of ∼216 kb) was calculated, and the top-3 genomic regions explaining the largest rate of the genetic variance were considered promising regions and used to identify potential candidate genes. Mean (SD) MU was 25.38 (8.02) mg/dL and 25.03 (8.06) mg/dL in the first and second lactation, respectively. Mean heritability estimates for daily MU were 0.21 and 0.23 for the first and second lactation, respectively. The genetic correlations estimated between MU and MY, milk composition, and CMP were quite low (ranged from -0.10 [CC] to 0.10 [TA] and -0.05 [CT] to 0.13 [TA] for the first and second lactations, respectively). The top-3 regions associated with MU were located from 80.61 to 80.74 Mb on BTA6, 103.26 to 103.41 Mb on BTA11, and 1.59 to 2.15 Mb on BTA14. Genes including KCNT1, MROH1, SHARPIN, TSSK5, CPSF1, HSF1, TONSL, DGAT1, SCX, and MAF1 were identified as positional candidate genes for MU. The findings of this study can be used for a better understanding of the genomic architecture underlying MU in Holstein cattle.
Collapse
Affiliation(s)
- H Atashi
- TERRA Teaching and Research Center, Gembloux Agro-Bio Tech, University of Liège, 5030 Gembloux, Belgium; Department of Animal Science, Shiraz University, 71441-13131 Shiraz, Iran.
| | - Y Chen
- TERRA Teaching and Research Center, Gembloux Agro-Bio Tech, University of Liège, 5030 Gembloux, Belgium
| | - S Vanderick
- TERRA Teaching and Research Center, Gembloux Agro-Bio Tech, University of Liège, 5030 Gembloux, Belgium
| | - X Hubin
- Elevéo asbl Awé Group, 5590 Ciney, Belgium
| | - N Gengler
- TERRA Teaching and Research Center, Gembloux Agro-Bio Tech, University of Liège, 5030 Gembloux, Belgium
| |
Collapse
|
3
|
Buitenhuis AJ, Poulsen NA. Estimation of heritability for milk urea and genetic correlations with milk production traits in 3 Danish dairy breeds. J Dairy Sci 2023; 106:5562-5569. [PMID: 37331871 DOI: 10.3168/jds.2022-22798] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Accepted: 02/17/2023] [Indexed: 06/20/2023]
Abstract
The aim of this study was to estimate genetic parameters for milk urea (MU) content in 3 main Danish dairy breeds. As a part of the Danish milk recording system, milk samples from cows on commercial farms were analyzed for MU concentration (mmol/L) and the percentages of fat and protein. There were 323,800 Danish Holstein, 70,634 Danish Jersey, and 27,870 Danish Red cows sampled with a total of 1,436,580, 368,251, and 133,922 test-day records per breed, respectively, included in the data set. Heritabilities for MU were low to moderate (0.22, 0.18, and 0.24 for the Holstein, Jersey, and Red breeds, respectively). The genetic correlation was close to zero between MU and milk yield in Jersey and Red, and -0.14 for Holstein. The genetic correlations between MU and fat and protein percentages, respectively, were positive for all 3 dairy breeds. Herd-test-day explained 51%, 54%, and 49% of the variation in MU in Holstein, Jersey, and Red, respectively. This indicates that MU levels in milk can be reduced by farm management. The current study shows that there are possibilities to influence MU by genetic selection as well as by farm management.
Collapse
Affiliation(s)
- A J Buitenhuis
- Center for Quantitative Genetics and Genomics, Aarhus University, DK-8000 Aarhus C, Denmark.
| | - N A Poulsen
- Department of Food Science, Aarhus University, DK-8200 Aarhus N, Denmark
| |
Collapse
|
4
|
Jahnel RE, Blunk I, Wittenburg D, Reinsch N. Relationship between milk urea content and important milk traits in Holstein cattle. Animal 2023; 17:100767. [PMID: 37141636 DOI: 10.1016/j.animal.2023.100767] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2022] [Revised: 02/27/2023] [Accepted: 02/28/2023] [Indexed: 05/06/2023] Open
Abstract
Breeding cattle with low nitrogen emissions has been proposed as a countermeasure against eutrophication due to dairy production. Milk urea content (MU) could potentially serve as a new readily measured indicator trait for nitrogen emissions by cows. Therefore, we estimated genetic parameters related to MU and its relationship with other milk traits. We analysed 4 178 735 milk samples collected between January 2008 and June 2019 from 261 866 German Holstein dairy cows during their first, second, and third lactations. Restricted maximum likelihood estimation was conducted using univariate and bivariate random regression sire models in WOMBAT. We obtained moderate average daily heritability estimates for the daily MU of 0.24 in first lactation cows, 0.23 in second lactation cows, and 0.21 in third lactation cows with average daily genetic SDs of 25.16 mg/kg, 24.93 mg/kg, and 23.75 mg/kg, respectively. Averaged over days in milk, the repeatability estimates were low at 0.41 in first, second, and third lactation cows. A strong positive genetic correlation was found between MU and milk urea yield (MUY; 0.72 on average). In addition, 305-day heritabilities were estimated as 0.50, 0.52, and 0.50 in first, second, and third lactation cows, respectively, with genetic correlations of 0.94 or higher for MU in different lactations. By contrast, the averaged estimates of the genetic correlations between MU and other milk traits were low (-0.07 to 0.15). Moderate heritability estimates clearly allow the possible selection for MU, and the near-zero estimates of genetic correlations indicate no risk of undesired correlated selection responses in other milk traits. However, a relationship still needs to be established between MU as an indicator trait and the target trait, defined as total individual nitrogen emissions.
Collapse
Affiliation(s)
- R E Jahnel
- Research Institute for Farm Animal Biology (FBN), Wilhelm-Stahl-Allee 2, 18196 Dummerstorf, Germany
| | - I Blunk
- Research Institute for Farm Animal Biology (FBN), Wilhelm-Stahl-Allee 2, 18196 Dummerstorf, Germany
| | - D Wittenburg
- Research Institute for Farm Animal Biology (FBN), Wilhelm-Stahl-Allee 2, 18196 Dummerstorf, Germany
| | - N Reinsch
- Research Institute for Farm Animal Biology (FBN), Wilhelm-Stahl-Allee 2, 18196 Dummerstorf, Germany.
| |
Collapse
|
5
|
Atashi H, Chen Y, Wilmot H, Vanderick S, Hubin X, Gengler N. Genome-wide association for milk urea concentration in Dual-Purpose Belgian Blue cows. J Anim Breed Genet 2022; 139:710-722. [PMID: 35834354 DOI: 10.1111/jbg.12732] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2021] [Accepted: 06/25/2022] [Indexed: 11/27/2022]
Abstract
The objectives of this study were to estimate genetic parameters and identify genomic regions associated with milk urea concentration (MU) in Dual-Purpose Belgian Blue (DPBB) cows. The data were 29,693 test-day records of milk yield (MY), fat yield (FY), protein yield (PY), fat percentage (FP), protein percentage (PP) and MU collected between 2014 and 2020 on 2498 first parity cows (16,935 test-day records) and 1939 second-parity cows (12,758 test-day records) from 49 herds in the Walloon Region of Belgium. Data of 28,266 single nucleotide polymorphisms (SNP), located on 29 Bos taurus autosomes (BTA), on 1699 animals (639 males and 1060 females) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method using a single chain of 100,000 iterations after a burn-in period of 20,000. SNP solutions were estimated using a single-step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by windows of 25 consecutive SNPs (with an average size of ~2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. The mean (SD) of MU was 22.89 (10.07) and 22.35 (10.07) mg/dl for first and second parity, respectively. The mean (SD) heritability estimates for daily MU were 0.18 (0.01) and 0.22 (0.02), for first and second parity, respectively. The mean (SD) genetic correlations between daily MU and MY, FY, PY, FP and PP were -0.05 (0.09), -0.07 (0.11), -0.03 (0.13), -0.05 (0.08) and -0.03 (0.11) for first parity, respectively. The corresponding values estimated for second parity were 0.02 (0.10), -0.02 (0.09), 0.02 (0.08), -0.08 (0.06) and -0.05 (0.05). The genome-wide association analyses identified three genomic regions (BTA2, BTA3 and BTA13) associated with MU. The identified regions showed contrasting results between parities and among different stages within each parity. This suggests that different groups of candidate genes underlie the phenotypic expression of MU between parities and among different lactation stages within a parity. The results of this study can be used for future implementation and use of genomic evaluation to reduce MU in DPBB cows.
Collapse
Affiliation(s)
- Hadi Atashi
- TERRA Teaching and Research Center, Gembloux Agro-Bio Tech, University of Liège, Gembloux, Belgium.,Department of Animal Science, Shiraz University, Shiraz, Iran
| | - Yansen Chen
- TERRA Teaching and Research Center, Gembloux Agro-Bio Tech, University of Liège, Gembloux, Belgium
| | - Hélène Wilmot
- TERRA Teaching and Research Center, Gembloux Agro-Bio Tech, University of Liège, Gembloux, Belgium.,National Fund for Scientific Research (F.R.S.-FNRS), Brussels, Belgium
| | - Sylvie Vanderick
- TERRA Teaching and Research Center, Gembloux Agro-Bio Tech, University of Liège, Gembloux, Belgium
| | | | - Nicolas Gengler
- TERRA Teaching and Research Center, Gembloux Agro-Bio Tech, University of Liège, Gembloux, Belgium
| |
Collapse
|
6
|
van den Berg I, Ho PN, Nguyen TV, Haile-Mariam M, Luke TDW, Pryce JE. Using mid-infrared spectroscopy to increase GWAS power to detect QTL associated with blood urea nitrogen. Genet Sel Evol 2022; 54:27. [PMID: 35436852 PMCID: PMC9014603 DOI: 10.1186/s12711-022-00719-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2021] [Accepted: 04/05/2022] [Indexed: 11/20/2022] Open
Abstract
Blood urea nitrogen (BUN) is an indicator trait for urinary nitrogen excretion. Measuring BUN level requires a blood sample, which limits the number of records that can be obtained. Alternatively, BUN can be predicted using mid-infrared (MIR) spectroscopy of a milk sample and thus records become available on many more cows through routine milk recording processes. The genetic correlation between MIR predicted BUN (MBUN) and BUN is 0.90. Hence, genetically, BUN and MBUN can be considered as the same trait. The objective of our study was to perform genome-wide association studies (GWAS) for BUN and MBUN, compare these two GWAS and detect quantitative trait loci (QTL) for both traits, and compare the detected QTL with previously reported QTL for milk urea nitrogen (MUN). The dataset used for our analyses included 2098 and 18,120 phenotypes for BUN and MBUN, respectively, and imputed whole-genome sequence data. The GWAS for MBUN was carried out using either the full dataset, the 2098 cows with records for BUN, or 2000 randomly selected cows, so that the dataset size is comparable to that for BUN. The GWAS results for BUN and MBUN were very different, in spite of the strong genetic correlation between the two traits. We detected 12 QTL for MBUN, on bovine chromosomes 2, 3, 9, 11, 12, 14 and X, and one QTL for BUN on chromosome 13. The QTL detected on chromosomes 11, 14 and X overlapped with QTL detected for MUN. The GWAS results were highly sensitive to the subset of records used. Hence, caution is warranted when interpreting GWAS based on small datasets, such as for BUN. MBUN may provide an attractive alternative to perform a more powerful GWAS to detect QTL for BUN.
Collapse
|
7
|
van den Berg I, Ho PN, Nguyen TV, Haile-Mariam M, MacLeod IM, Beatson PR, O'Connor E, Pryce JE. GWAS and genomic prediction of milk urea nitrogen in Australian and New Zealand dairy cattle. Genet Sel Evol 2022; 54:15. [PMID: 35183113 PMCID: PMC8858489 DOI: 10.1186/s12711-022-00707-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2021] [Accepted: 01/31/2022] [Indexed: 11/24/2022] Open
Abstract
Background Urinary nitrogen leakage is an environmental concern in dairy cattle. Selection for reduced urinary nitrogen leakage may be done using indicator traits such as milk urea nitrogen (MUN). The result of a previous study indicated that the genetic correlation between MUN in Australia (AUS) and MUN in New Zealand (NZL) was only low to moderate (between 0.14 and 0.58). In this context, an alternative is to select sequence variants based on genome-wide association studies (GWAS) with a view to improve genomic prediction accuracies. A GWAS can also be used to detect quantitative trait loci (QTL) associated with MUN. Therefore, our objectives were to perform within-country GWAS and a meta-GWAS for MUN using records from up to 33,873 dairy cows and imputed whole-genome sequence data, to compare QTL detected in the GWAS for MUN in AUS and NZL, and to use sequence variants selected from the meta-GWAS to improve the prediction accuracy for MUN based on a joint AUS-NZL reference set. Results Using the meta-GWAS, we detected 14 QTL for MUN, located on chromosomes 1, 6, 11, 14, 19, 22, 26 and the X chromosome. The three most significant QTL encompassed the casein genes on chromosome 6, PAEP on chromosome 11 and DGAT1 on chromosome 14. We selected 50,000 sequence variants that had the same direction of effect for MUN in AUS and MUN in NZL and that were most significant in the meta-analysis for the GWAS. The selected sequence variants yielded a genetic correlation between MUN in AUS and MUN in NZL of 0.95 and substantially increased prediction accuracy in both countries. Conclusions Our results demonstrate how the sharing of data between two countries can increase the power of a GWAS and increase the accuracy of genomic prediction using a multi-country reference population and sequence variants selected based on a meta-GWAS. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-022-00707-9.
Collapse
Affiliation(s)
- Irene van den Berg
- Centre for AgriBioscience, Agriculture Victoria, 5 Ring Road, Bundoora, AgriBioVIC, 3083, Australia.
| | - Phuong N Ho
- Centre for AgriBioscience, Agriculture Victoria, 5 Ring Road, Bundoora, AgriBioVIC, 3083, Australia
| | - Tuan V Nguyen
- Centre for AgriBioscience, Agriculture Victoria, 5 Ring Road, Bundoora, AgriBioVIC, 3083, Australia
| | - Mekonnen Haile-Mariam
- Centre for AgriBioscience, Agriculture Victoria, 5 Ring Road, Bundoora, AgriBioVIC, 3083, Australia
| | - Iona M MacLeod
- Centre for AgriBioscience, Agriculture Victoria, 5 Ring Road, Bundoora, AgriBioVIC, 3083, Australia
| | | | | | - Jennie E Pryce
- Centre for AgriBioscience, Agriculture Victoria, 5 Ring Road, Bundoora, AgriBioVIC, 3083, Australia.,School of Applied Systems Biology, La Trobe University, Bundoora, VIC, 3083, Australia
| |
Collapse
|