The Resilient Dairy Genome Project-A general overview of methods and objectives related to feed efficiency and methane emissions.

The Resilient Dairy Genome Project (RDGP) is an international large-scale applied research project that aims to generate genomic tools to breed more resilient dairy cows. In this context, improving feed efficiency and reducing greenhouse gases from dairy is a high priority. The inclusion of traits related to feed efficiency (e.g., dry matter intake [DMI]) or greenhouse gases (e.g., methane emissions [CH4]) relies on available genotypes as well as high quality phenotypes. Currently, 7 countries (i.e., Australia, Canada, Denmark, Germany, Spain, Switzerland, and United States) contribute with genotypes and phenotypes including DMI and CH4. However, combining data are challenging due to differences in recording protocols, measurement technology, genotyping, and animal management across sources. In this study, we provide an overview of how the RDGP partners address these issues to advance international collaboration to generate genomic tools for resilient dairy. Specifically, we describe the current state of the RDGP database, data collection protocols in each country, and the strategies used for managing the shared data. As of February 2022, the database contains 1,289,593 DMI records from 12,687 cows and 17,403 CH4 records from 3,093 cows and continues to grow as countries upload new data over the coming years. No strong genomic differentiation between the populations was identified in this study, which may be beneficial for eventual across-country genomic predictions. Moreover, our results reinforce the need to account for the heterogeneity in the DMI and CH4 phenotypes in genomic analysis.

Genetic analysis of production traits and body size measurements and their relationships with metabolic diseases in German Holstein cattle.

This study sheds light on the genetic complexity and interplay of production, body size, and metabolic health in dairy cattle. Phenotypes for body size-related traits from conformation classification (130,166 animals) and production (101,562 animals) of primiparous German Holstein cows were available. Additionally, 21,992, 16,641, and 7,096 animals were from herds with recordings of the metabolic diseases ketosis, displaced abomasum, and milk fever in first, second, and third lactation. Moreover, all animals were genotyped. Heritabilities of traits and genetic correlations between all traits were estimated and GWAS were performed. Heritability was between 0.240 and 0.333 for production and between 0.149 and 0.368 for body size traits. Metabolic diseases were lowly heritable, with estimates ranging from 0.011 to 0.029 in primiparous cows, from 0.008 to 0.031 in second lactation, and from 0.037 to 0.052 in third lactation. Production was found to have negative genetic correlations with body condition score (BCS; -0.279 to -0.343) and udder depth (-0.348 to -0.419). Positive correlations were observed for production and body depth (0.138-0.228), dairy character (DCH) (0.334-0.422), and stature (STAT) (0.084-0.158). In first parity cows, metabolic disease traits were unfavorably correlated with production, with genetic correlations varying from 0.111 to 0.224, implying that higher yielding cows have more metabolic problems. Genetic correlations of disease traits in second and third lactation with production in primiparous cows were low to moderate and in most cases unfavorable. While BCS was negatively correlated with metabolic diseases (-0.255 to -0.470), positive correlations were found between disease traits and DCH (0.269-0.469) as well as STAT (0.172-0.242). Thus, the results indicate that larger and sharper animals with low BCS are more susceptible to metabolic disorders. Genome-wide association studies revealed several significantly associated SNPs for production and conformation traits, confirming previous findings from literature. Moreover, for production and conformation traits, shared significant signals on Bos taurus autosome (BTA) 5 (88.36 Mb) and BTA 6 (86.40 to 87.27 Mb) were found, implying pleiotropy. Additionally, significant SNPs were observed for metabolic diseases on BTA 3, 10, 14, 17, and 26 in first lactation and on BTA 2, 6, 8, 17, and 23 in third lactation. Overall, this study provides important insights into the genetic basis and interrelations of relevant traits in today's Holstein cattle breeding programs, and findings may help to improve selection decisions.

A genomic assessment of the correlation between milk production traits and claw and udder health traits in Holstein dairy cattle.

Claw diseases and mastitis represent the most important disease traits in dairy cattle with increasing incidences and a frequently mentioned connection to milk yield. Yet, many studies aimed to detect the genetic background of both trait complexes via fine-mapping of quantitative trait loci. However, little is known about genomic regions that simultaneously affect milk production and disease traits. For this purpose, several tools to detect local genetic correlations have been developed. In this study, we attempted a detailed analysis of milk production and disease traits as well as their interrelationship using a sample of 34,497 50K genotyped German Holstein cows with milk production and claw and udder disease traits records. We performed a pedigree-based quantitative genetic analysis to estimate heritabilities and genetic correlations. Additionally, we generated GWAS summary statistics, paying special attention to genomic inflation, and used these data to identify shared genomic regions, which affect various trait combinations. The heritability on the liability scale of the disease traits was low, between 0.02 for laminitis and 0.19 for interdigital hyperplasia. The heritabilities for milk production traits were higher (between 0.27 for milk energy yield and 0.48 for fat-protein ratio). Global genetic correlations indicate the shared genetic effect between milk production and disease traits on a whole genome level. Most of these estimates were not significantly different from zero, only mastitis showed a positive one to milk (0.18) and milk energy yield (0.13), as well as a negative one to fat-protein ratio (-0.07). The genomic analysis revealed significant SNPs for milk production traits that were enriched on Bos taurus autosome 5, 6, and 14. For digital dermatitis, we found significant hits, predominantly on Bos taurus autosome 5, 10, 22, and 23, whereas we did not find significantly trait-associated SNPs for the other disease traits. Our results confirm the known genetic background of disease and milk production traits. We further detected 13 regions that harbor strong concordant effects on a trait combination of milk production and disease traits. This detailed investigation of genetic correlations reveals additional knowledge about the localization of regions with shared genetic effects on these trait complexes, which in turn enables a better understanding of the underlying biological pathways and putatively the utilization for a more precise design of breeding schemes.

Microbial signature inferred from genomic breeding selection on milk urea concentration and its relation to proxies of nitrogen-utilization efficiency in Holsteins.

Increasing the nitrogen-utilization efficiency (NUE) of dairy cows by breeding selection would offer advantages from nutritional, environmental, and economic perspectives. Because data collection of NUE phenotypes is not feasible in large cow cohorts, the cow individual milk urea concentration (MU) has been suggested as an indicator trait. Considering the symbiotic interplay between dairy cows and their rumen microbiome, individual MU was thought to be influenced by host genetics and by the rumen microbiome, the latter in turn being partly attributed to host genetics. To enhance our knowledge of MU as an indicator trait for NUE, we aimed to identify differential abundant rumen microbial genera between Holstein cows with divergent genomic breeding values for MU (GBVMU; GBVHMU vs. GBVLMU, where H and L indicate high and low MU phenotypes, respectively). The microbial genera identified were further investigated for their correlations with MU and 7 additional NUE-associated traits in urine, milk, and feces in 358 lactating Holsteins. Statistical analysis of microbial 16S rRNA amplicon sequencing data revealed significantly higher abundances of the ureolytic genus Succinivibrionaceae UCG-002 in GBVLMU cows, whereas GBVHMU animals hosted higher abundances of Clostridia unclassified and Desulfovibrio. The entire discriminating ruminal signature of 24 microbial taxa included a further 3 genera of the Lachnospiraceae family that revealed significant correlations to MU values and were therefore proposed as considerable players in the GBVMU-microbiome-MU axis. The significant correlations of Prevotellaceae UCG-003, Anaerovibrio, Blautia, and Butyrivibrio abundances with MU measurements, milk nitrogen, and N content in feces suggested their contribution to genetically determined N-utilization in Holstein cows. The microbial genera identified might be considered for future breeding programs to enhance NUE in dairy herds.

Single-step SNPBLUP evaluation in six German beef cattle breeds.

The implementation of genomic selection for six German beef cattle populations was evaluated. Although the multiple-step implementation of genomic selection is the status quo in most national dairy cattle evaluations, the breeding structure of German beef cattle, coupled with the shortcoming and complexity of the multiple-step method, makes single step a more attractive option to implement genomic selection in German beef cattle populations. Our objective was to develop a national beef cattle single-step genomic evaluation in five economically important traits in six German beef cattle populations and investigate its impact on the accuracy and bias of genomic evaluations relative to the current pedigree-based evaluation. Across the six breeds in our study, 461,929 phenotyped and 14,321 genotyped animals were evaluated with a multi-trait single-step model. To validate the single-step model, phenotype data in the last 2 years were removed in a forward validation study. For the conventional and single-step approaches, the genomic estimated breeding values of validation animals and other animals were compared between the truncated and the full evaluations. The correlation of the GEBVs between the full and truncated evaluations in the validation animals was slightly higher in the single-step evaluation. The regression of the full GEBVs on truncated GEBVs was close to the optimal value of 1 for both the pedigree-based and the single-step evaluations. The SNP effect estimates from the truncated evaluation were highly correlated with those from the full evaluation, with values ranging from 0.79 to 0.94. The correlation of the SNP effect was influenced by the number of genotyped animals shared between the full and truncated evaluations. The regression coefficients of the SNP effect of the full evaluation on the truncated evaluation were all close to the expected value of 1, indicating unbiased estimates of the SNP markers for the production traits. The Manhattan plot of the SNP effect estimates identified chromosomal regions harbouring major genes for muscling and body weight in breeds of French origin. Based on the regression intercept and slope of the GEBVs of validation animals, the single-step evaluation was neither inflated nor deflated across the six breeds. Overall, the single-step model resulted in a more accurate and stable evaluation. However, due to the small number of genotyped individuals, the single-step method only provided slightly better results when compared to the pedigree-based method.

Mendelian sampling covariability of marker effects and genetic values.

BACKGROUND: Measures of the expected genetic variability among full-sibs are of practical relevance, such as in the context of mating decisions. An important application field in animal and plant breeding is the selection and allocation of mates when large or small amounts of genetic variability among offspring are desired, depending on user-specific goals. Estimates of the Mendelian sampling variance can be obtained by simulating gametes from parents with known diplotypes. Knowledge of recombination rates and additive marker effects is also required. In this study, we aimed at developing an exact method that can account for both additive and dominance effects. RESULTS: We derived parent-specific covariance matrices that exactly quantify the within-family (co-)variability of additive and dominance marker effects. These matrices incorporate prior knowledge of the parental diplotypes and recombination rates. When combined with additive marker effects, they allow the exact derivation of the Mendelian sampling (co-)variances of (estimated) breeding values for several traits, as well for the aggregate genotype. A comparative analysis demonstrated good average agreement between the exact values and the simulation results for a practical dataset (74,353 German Holstein cattle). CONCLUSIONS: The newly derived method is suitable for calculating the exact amount of intra-family variation of the estimated breeding values and genetic values (comprising additive and dominance effects).

Prediction of expected genetic variation within groups of offspring for innovative mating schemes.

BACKGROUND: Experience from progeny-testing indicates that the mating of popular bull sires that have high estimated breeding values with excellent dams does not guarantee the production of offspring with superior breeding values. This is explained partly by differences in the standard deviation of gamete breeding values (SDGBV) between animals at the haplotype level. The SDGBV depends on the variance of the true effects of single nucleotide polymorphisms (SNPs) and the degree of heterozygosity. Haplotypes of 58 035 Holstein animals were used to predict and investigate expected SDGBV for fat yield, protein yield, somatic cell score and the direct genetic effect for stillbirth. RESULTS: Differences in SDGBV between animals were detected, which means that the groups of offspring of parents with low SDGBV will be more homogeneous than those of parents with high SDGBV, although the expected mean breeding values of the progeny will be the same. SDGBV was negatively correlated with genomic and pedigree inbreeding coefficients and a small loss of SDGBV over time was observed. Sires that had relatively low mean gamete breeding values but high SDGBV had a higher probability of producing extremely positive offspring than sires that had a high mean gamete breeding value and low SDGBV. CONCLUSIONS: An animal's SDGBV can be estimated based on genomic information and used to design specific genomic mating plans. Estimated SDGBV are an additional tool for mating programs, which allows breeders to identify and match mating partners using specific haplotype information.

Ruminal background of predisposed milk urea (MU) concentration in Holsteins.

Efforts to reduce nitrogen (N) emissions are currently based on the optimization of dietary- N supply at average herd N requirements. The implementation of the considerable individual differences and predispositions in N- use efficiency and N- excretion in breeding programs is hampered by the difficulty of data collection. Cow individual milk urea (MU) concentration has been proposed as an easy-to-measure surrogate trait, but recent studies questioned its predictive power. Therefore, a deeper understanding of the biological mechanisms underlying predisposed higher (HMUg) or lower (LMUg) MU concentration in dairy cows is needed. Considering the complex N- metabolism in ruminants, the distinction between HMUg and LMUg could be based on differences in (i) the rumen microbial community, (ii) the host-specific transcription processes in the rumen villi, and (iii) the host-microbe interaction in the rumen. Therefore, rumen fluid and rumen epithelial samples from 10 HMUg and 10 LMUg cows were analyzed by 16S sequencing and HiSeq sequencing. In addition, the effect of dietary-N reduction on ruminal shifts was investigated in a second step. In total, 10 differentially abundant genera (DAG) were identified between HMUg and LMUg cows, elucidating greater abundances of ureolytic Succinivibrionaceae_UCG-002 and Ruminococcaceae_unclassified in LMUg animals and enhanced occurrences of Butyvibrio in HMUg cows. Differential expression analysis revealed genes of the bovine Major Histocompatibility Complex (BOLA genes) as well as MX1, ISG15, and PRSS2 displaying candidates of MU predisposition that further attributed to enhanced immune system activities in LMUg cows. A number of significant correlations between microbial genera and host transcript abundances were uncovered, including strikingly positive correlations of BOLA-DRA transcripts with Roseburia and Lachnospiraceae family abundances that might constitute particularly prominent microbial-host interplays of MU predisposition. The reduction of feed-N was followed by 18 DAG in HMUg and 19 DAG in LMUg, depicting pronounced interest on Shuttleworthia, which displayed controversial adaption in HMUg and LMUg cows. Lowering feed-N further elicited massive downregulation of immune response and energy metabolism pathways in LMUg. Considering breeding selection strategies, this study attributed information content to MU about predisposed ruminal N-utilization in Holstein-Friesians.

Long-Term Impact of Optimum Contribution Selection Strategies on Local Livestock Breeds with Historical Introgression Using the Example of German Angler Cattle.

The long-term performance of different selection strategies was evaluated via simulation using the example of a local cattle breed, German Angler cattle. Different optimum contribution selection (OCS) approaches to maximize genetic gain were compared to a reference scenario without selection and truncation selection. The kinships and migrant contribution (MC) were estimated from genomic data. Truncation selection achieved the highest genetic gain but decreased diversity considerably at native alleles. It also caused the highest increase in MCs. Traditional OCS, which only constrains kinship, achieved almost the same genetic gain but also caused a small increase of MC and remarkably reduced the diversity of native alleles. When MC was required not to increase and the increase of kinship at native alleles was restricted, the MC levels and the diversity at native alleles were well managed, and the genetic gain was only slightly reduced. However, genetic progress was substantially lower in the scenario that aimed to recover the original genetic background. Truncation selection and traditional OCS selection both reduce the genetic originality of breeds with historical introgression. The inclusion of MC and kinship at native alleles as additional constraints in OCS showed great potential for conservation. Recovery of the original genetic background is possible but requires many generations of selection and reduces the genetic progress in performance traits. Hence, constraining MCs at their current values can be recommended to avoid further reduction of genetic originality.

The mammalian cervical vertebrae blueprint depends on the T (brachyury) gene.

A key common feature all but three known mammalian genera is the strict seven cervical vertebrae blueprint, suggesting the involvement of strong conserving selection forces during mammalian radiation. This is further supported by reports indicating that children with cervical ribs die before they reach reproductive age. Hypotheses were put up, associating cervical ribs (homeotic transformations) to embryonal cancer (e.g., neuroblastoma) or ascribing the constraint in cervical vertebral count to the development of the mammalian diaphragm. Here, we describe a spontaneous mutation c.196A > G in the Bos taurus T gene (also known as brachyury) associated with a cervical vertebral homeotic transformation that violates the fundamental mammalian cervical blueprint, but does not preclude reproduction of the affected individual. Genome-wide mapping, haplotype tracking within a large pedigree, resequencing of target genome regions, and bioinformatic analyses unambiguously confirmed the mutant c.196G allele as causal for this previously unknown defect termed vertebral and spinal dysplasia (VSD) by providing evidence for the mutation event. The nonsynonymous VSD mutation is located within the highly conserved T box of the T gene, which plays a fundamental role in eumetazoan body organization and vertebral development. To our knowledge, VSD is the first unequivocally approved spontaneous mutation decreasing cervical vertebrae number in a large mammal. The spontaneous VSD mutation in the bovine T gene is the first in vivo evidence for the hypothesis that the T protein is directly involved in the maintenance of the mammalian seven-cervical vertebra blueprint. It therefore furthers our knowledge of the T-protein function and early mammalian notochord development.