Genome-wide associations and detection of potential candidate genes for direct genetic and maternal genetic effects influencing dairy cattle body weight at different ages

Maternal effects and its importance in the genetic evaluations of preweaning live weight traits of beef cattle. A review

Maternal effects in cattle genetics are defined as the causal influence of the phenotype or maternal genotype on the offspring's phenotype by effects occurring when the genetic and environmental characteristics of the mother influence the phenotype of the offspring beyond the direct inheritance of genes. Its relevance has been strongly described in genetic models focused on the genetic improvement of preweaning traits in cow-calf beef cattle production systems. Here, basic concepts and the importance of maternal effects when using linear and animal model procedures for genetic evaluations of growth and live-weight traits in beef cattle are reviewed and discussed. A brief history of estimation methods from classical studies to recent studies used for the development of animal models for studying maternal effects is also provided. Some important biometric concepts for maternal effect estimation are described, and the antagonism between direct genetic effects and maternal effects, its biological basis, and sources of error in the estimation of direct genetic and maternal covariance are discussed. Finally, some genomic perspectives are presented.

Genome-Wide Association Study of Conformation Traits in Brazilian Holstein Cattle

The linear conformation of animals exerts an influence on health, reproduction, production, and welfare, in addition to longevity, which directly affects the profitability of milk-producing farms. The objectives of this study were (1) to perform genome-wide association studies (GWASs) of conformation traits, namely the Rump, Feet and Legs, Mammary System, Dairy Strength, and Final Classification traits, and (2) to identify genes and related pathways involved in physiological processes associated with conformation traits in Brazilian Holstein cattle. Phenotypic and genotypic data from 2339 Holstein animals distributed across the states of Rio Grande do Sul, Paraná, São Paulo, and Minas Gerais were used. The genotypic data were obtained with a 100 K SNP marker panel. The single-step genome-wide association study (ssGWAS) method was employed in the analyses. Genes close to a significant SNP were identified in an interval of 100 kb up- and downstream using the Ensembl database available in the BioMart tool. The DAVID database was used to identify the main metabolic pathways and the STRING program was employed to create the gene regulatory network. In total, 36 significant SNPs were found on 15 chromosomes; 27 of these SNPs were linked to genes that may influence the traits studied. Fourteen genes most closely related to the studied traits were identified, as well as four genes that showed interactions in important metabolic pathways such as myogenesis, adipogenesis, and angiogenesis. Among the total genes, four were associated with myogenesis (TMOD2, TMOD3, CCND2, and CTBP2), three with angiogenesis (FGF23, FGF1, and SCG3), and four with adipogenesis and body size and development (C5H12orf4, CCND2, EMILIN1, and FGF6). These results contribute to a better understanding of the biological mechanisms underlying phenotypic variability in conformation traits in Brazilian Holstein cattle.

Mendelian randomization analysis of 34,497 German Holstein cows to infer causal associations between milk production and health traits

BACKGROUND

Claw diseases and mastitis represent the most important health issues in dairy cattle with a frequently mentioned connection to milk production. Although many studies have aimed at investigating this connection in more detail by estimating genetic correlations, they do not provide information about causality. An alternative is to carry out Mendelian randomization (MR) studies using genetic variants to investigate the effect of an exposure on an outcome trait mediated by genetic variants. No study has yet investigated the causal association of milk yield (MY) with health traits in dairy cattle. Hence, we performed a MR analysis of MY and seven health traits using imputed whole-genome sequence data from 34,497 German Holstein cows. We applied a method that uses summary statistics and removes horizontal pleiotropic variants (having an effect on both traits), which improves the power and unbiasedness of MR studies. In addition, genetic correlations between MY and each health trait were estimated to compare them with the estimates of causal effects that we expected.

RESULTS

All genetic correlations between MY and each health trait were negative, ranging from - 0.303 (mastitis) to - 0.019 (digital dermatitis), which indicates a reduced health status as MY increases. The only non-significant correlation was between MY and digital dermatitis. In addition, each causal association was negative, ranging from - 0.131 (mastitis) to - 0.034 (laminitis), but the number of significant associations was reduced to five nominal and two experiment-wide significant results. The latter were between MY and mastitis and between MY and digital phlegmon. Horizontal pleiotropic variants were identified for mastitis, digital dermatitis and digital phlegmon. They were located within or nearby variants that were previously reported to have a horizontal pleiotropic effect, e.g., on milk production and somatic cell count.

CONCLUSIONS

Our results confirm the known negative genetic connection between health traits and MY in dairy cattle. In addition, they provide new information about causality, which for example points to the negative energy balance mediating the connection between these traits. This knowledge helps to better understand whether the negative genetic correlation is based on pleiotropy, linkage between causal variants for both trait complexes, or indeed on a causal association.

Including microbiome information in a multi-trait genomic evaluation: a case study on longitudinal growth performance in beef cattle

BACKGROUND

Growth rate is an important component of feed conversion efficiency in cattle and varies across the different stages of the finishing period. The metabolic effect of the rumen microbiome is essential for cattle growth, and investigating the genomic and microbial factors that underlie this temporal variation can help maximize feed conversion efficiency at each growth stage.

RESULTS

By analysing longitudinal body weights during the finishing period and genomic and metagenomic data from 359 beef cattle, our study demonstrates that the influence of the host genome on the functional rumen microbiome contributes to the temporal variation in average daily gain (ADG) in different months (ADG1, ADG2, ADG3, ADG4). Five hundred and thirty-three additive log-ratio transformed microbial genes (alr-MG) had non-zero genomic correlations (rg) with at least one ADG-trait (ranging from |0.21| to |0.42|). Only a few alr-MG correlated with more than one ADG-trait, which suggests that a differential host-microbiome determinism underlies ADG at different stages. These alr-MG were involved in ribosomal biosynthesis, energy processes, sulphur and aminoacid metabolism and transport, or lipopolysaccharide signalling, among others. We selected two alternative subsets of 32 alr-MG that had a non-uniform or a uniform rg sign with all the ADG-traits, regardless of the rg magnitude, and used them to develop a microbiome-driven breeding strategy based on alr-MG only, or combined with ADG-traits, which was aimed at shaping the rumen microbiome towards increased ADG at all finishing stages. Combining alr-MG information with ADG records increased prediction accuracy of genomic estimated breeding values (GEBV) by 11 to 22% relative to the direct breeding strategy (using ADG-traits only), whereas using microbiome information, only, achieved lower accuracies (from 7 to 41%). Predicted selection responses varied consistently with accuracies. Restricting alr-MG based on their rg sign (uniform subset) did not yield a gain in the predicted response compared to the non-uniform subset, which is explained by the absence of alr-MG showing non-zero rg at least with more than one of the ADG-traits.

CONCLUSIONS

Our work sheds light on the role of the microbial metabolism in the growth trajectory of beef cattle at the genomic level and provides insights into the potential benefits of using microbiome information in future genomic breeding programs to accurately estimate GEBV and increase ADG at each finishing stage in beef cattle.

Multitrait meta-analyses identify potential candidate genes for growth-related traits in Holstein heifers

Understanding the underlying pleiotropic relationships among growth and body size traits is important for refining breeding strategies in dairy cattle for optimal body size and growth rate. Therefore, we performed single-trait GWAS for monthly-recorded body weight (BW), hip height, body length, and chest girth from birth to 12 mo of age in Holstein animals, followed by stepwise multiple regression of independent or lowly-linked markers from GWAS loci using conditional and joint association analyses (COJO). Subsequently, we conducted a multitrait meta-analysis to detect pleiotropic markers. Based on the single-trait GWAS, we identified 170 significant SNPs, in which 59 of them remained significant after the COJO analyses. The most significant SNP, located at BTA7:3,676,741, explained 2.93% of the total phenotypic variance for BW6 (BW at 6 mo of age). We identified 17 SNPs with potential pleiotropic effects based on the multitrait meta-analyses, which resulted in 3 additional SNPs in comparison to those detected based on the single-trait GWAS. The identified quantitative trait loci regions overlap with genes known to influence human growth-related traits. According to positional and functional analyses, we proposed HMGA2, HNF4G, MED13L, BHLHE40, FRZB, DMP1, TRIB3, and GATAD2A as important candidate genes influencing the studied traits. The combination of single-trait GWAS and meta-analyses of GWAS results improved the efficiency of detecting associated SNPs, and provided new insights into the genetic mechanisms of growth and development in Holstein cattle.

Genomic signatures of selection, local adaptation and production type characterisation of East Adriatic sheep breeds

BACKGROUND

The importance of sheep breeding in the Mediterranean part of the eastern Adriatic has a long tradition since its arrival during the Neolithic migrations. Sheep production system is extensive and generally carried out in traditional systems without intensive systematic breeding programmes for high uniform trait production (carcass, wool and milk yield). Therefore, eight indigenous Croatian sheep breeds from eastern Adriatic treated here as metapopulation (EAS), are generally considered as multipurpose breeds (milk, meat and wool), not specialised for a particular type of production, but known for their robustness and resistance to certain environmental conditions. Our objective was to identify genomic regions and genes that exhibit patterns of positive selection signatures, decipher their biological and productive functionality, and provide a "genomic" characterization of EAS adaptation and determine its production type.

RESULTS

We identified positive selection signatures in EAS using several methods based on reduced local variation, linkage disequilibrium and site frequency spectrum (eROHi, iHS, nSL and CLR). Our analyses identified numerous genomic regions and genes (e.g., desmosomal cadherin and desmoglein gene families) associated with environmental adaptation and economically important traits. Most candidate genes were related to meat/production and health/immune response traits, while some of the candidate genes discovered were important for domestication and evolutionary processes (e.g., HOXa gene family and FSIP2). These results were also confirmed by GO and QTL enrichment analysis.

CONCLUSIONS

Our results contribute to a better understanding of the unique adaptive genetic architecture of EAS and define its productive type, ultimately providing a new opportunity for future breeding programmes. At the same time, the numerous genes identified will improve our understanding of ruminant (sheep) robustness and resistance in the harsh and specific Mediterranean environment.

Effects of prenatal heat stress on birth weight and birth weight genetic parameters in German Holstein calves

The aim of this study was to infer the effects of heat stress (HS) during late gestation of dams on phenotypes and on direct and maternal genetic parameters for birth weight (BiW). We considered 171,221 Holstein calves kept in 56 large-scale co-operator herds. For a clear separation of maternal effects, only calves from dams with at least 3 offspring were included in the analyses. The genotype data set comprised 41,143 SNPs from 1,883 Holstein bulls. Temperature-humidity indices (THI) during the last 8 wk of gestation were calculated in each herd to reflect prenatal HS. A further prenatal HS descriptor was the first principal component (PC1) from principal component analysis considering the daily THI during the last 56 d of gestation. Regression coefficients of BiW on prenatal THI during the last 12 wk of gestation and PC1 were estimated in 13 consecutive phenotypic analyses. The strongest BiW decline was -0.63 kg per standardized THI, identified during 50 to 56 d before birth. A reaction norm model with weekly prenatal THI or PC1 nested within maternal genetic and maternal permanent environmental effects was defined to infer maternal sensitivity in response to prenatal THI alterations. Direct BiW heritabilities were close to 0.33 in the course of prenatal THI. Maternal BiW heritabilities marginally increased from 0.07 to 0.08 with increasing THI. Genetic correlations between maternal genetic effects at maximum HS levels and remaining THI were larger than 0.95, indicating the absence of genotype by time-lagged HS interactions. In contrast, maternal permanent environmental correlations between BiW at prenatal THI indicating HS with BiW at remaining THI substantially declined with increasing THI distances. Hence, from a herd management perspective, avoiding HS during the dry period of the dams will contribute to a slight increase in fetus growth.

Genome-Wide Association Study of Maternal Genetic Effects on Intramuscular Fat and Fatty Acid Composition in Rabbits

Maternal genetic effects (MGE) could affect meat quality traits such as intramuscular fat (IMF) and its fatty acid composition. However, it has been scarcely studied, especially in rabbits. The objectives of the present study were, first, to assess the importance of MGE on intramuscular fat and fatty acid composition by applying a Bayesian maternal animal model in two rabbit lines divergently selected for IMF. The second objective was to identify genomic regions and candidate genes of MGE that are associated with the traits of these offspring, using Bayesian methods in a Genome Wide Association Study (GWAS). Quantitative analyses were performed using data from 1982 rabbits, and 349 animals from the 9th generation and 76 dams of the 8th generation with 88,512 SNPs were used for the GWAS. The studied traits were IMF, saturated fatty acids (total SFA, C14:0; myristic acid, C16:0; palmitic acid and C18:0; stearic acid), monounsaturated fatty acids (total MUFA, C16:1n-7; palmitoleic acid and C18:1n-9; oleic acid), polyunsaturated fatty acids (total PUFA, C18:2n-6; linoleic acid, C18:3n-3; α-linolenic acid and C20:4n-6; arachidonic acid), MUFA/SFA and PUFA/SFA. The proportion of phenotypic variance explained by the maternal genetic effect ranged from 8 to 22% for IMF, depending on the model. For fatty acid composition, the proportion of phenotypic variance explained by maternal genetic effects varied from 10% (C18:0) to 46% (MUFA) in a model including both direct and additive maternal genetic effects, together with the common litter effect as a random variable. In particular, there were significant direct maternal genetic correlations for C16:0, C18:1n9, C18:2n6, SFA, MUFA, and PUFA with values ranging from -0.53 to -0.89. Relevant associated genomic regions were located on the rabbit chromosomes (OCU) OCU1, OCU5 and OCU19 containing some relevant candidates (TANC2, ACE, MAP3K3, TEX2, PRKCA, SH3GL2, CNTLN, RPGRIP1L and FTO) related to lipid metabolism, binding, and obesity. These regions explained about 1.2 to 13.9% of the total genomic variance of the traits studied. Our results showed an important maternal genetic effect on IMF and its fatty acid composition in rabbits and identified promising candidate genes associated with these traits.

Genome-wide identification and evolutionary analysis of the FGF gene family in buffalo

Fibroblast growth factors (FGFs) are important polypeptide growth factors that play a critical role in many developmental processes, including differentiation, cell proliferation, and migration in mammals. This study employs in silico analyses to characterize the FGF gene family in buffalo, investigating their genome-wide identification, physicochemical properties, and evolutionary patterns. For this purpose, genomic and proteomic sequences of buffalo, cattle, goat, and sheep were retrieved from NCBI database. We identified a total of 22 FGF genes in buffalo. Physicochemical properties observed through ProtParam tool showed notable features of these proteins including in-vitro instability, thermostability, hydrophilicity, and basic nature. Phylogenetic analysis grouped 22 identified genes into nine sub-families based on evolutionary relationships. Additionally, analysis of gene structure, motif patterns, and conserved domains using TBtools revealed the remarkable conservation of this gene family across selected species throughout the course of evolution. Comparative amino acid analysis performed through ClustalW demonstrated significant conservation between buffalo and cattle FGF proteins. Mutational analysis showed three non-synonymous mutations at positions R103 > G, P7 > L, and E98 > Q in FGF4, FGF6, and FGF19, respectively in buffalo. Duplication events revealed only one segmental duplication (FGF10/FGF22) in buffalo and two in cattle (FGF10/FGF22 and FGF13/FGF13-like) with Ka/Ks values <1 indicating purifying selection pressure for these duplications. Comparison of protein structures of buffalo, goat, and sheep exhibited more similarities in respective structures. In conclusion, our study highlights the conservation of the FGF gene family in buffalo during evolution. Furthermore, the identified non-synonymous mutations may have implications for the selection of animals with better performance.Communicated by Ramaswamy H. Sarma.

Longitudinal genome-wide association analysis using a single-step random regression model for height in Japanese Holstein cattle

Growth traits, such as body weight and height, are essential in the design of genetic improvement programs of dairy cattle due to their relationship with feeding efficiency, longevity, and health. We investigated genomic regions influencing height across growth stages in Japanese Holstein cattle using a single-step random regression model. We used 72,921 records from birth to 60 mo of age with 4,111 animals born between 2000 and 2016. The analysis included 1,410 genotyped animals with 35,319 single nucleotide polymorphisms, consisting of 883 females with records and 527 bulls, and 30,745 animals with pedigree information. A single genomic region at the 58.4 megabase pair on chromosome 18 was consistently identified across 6 age points of 10, 20, 30, 40, 50, and 60 mo after multiple testing corrections for the significance threshold. Twelve candidate genes, previously reported for longevity and gestation length, were found near the identified genomic region. Another location near the identified region was also previously associated with body conformation, fertility, and calving difficulty. Functional Gene Ontology enrichment analysis suggested that the candidate genes regulate dephosphorylation and phosphatase activity. Our findings show that further study of the identified candidate genes will contribute to a better understanding of the genetic basis of height in Japanese Holstein cattle.

Invited review: Good practices in genome-wide association studies to identify candidate sequence variants in dairy cattle

Genotype data from dairy cattle selection programs have greatly facilitated GWAS to identify variants related to economic traits. Results can enhance the accuracy of genomic prediction, analyze more complex models that go beyond additive effects, elucidate the genetic architecture of a trait, and finally, decipher the underlying biology of traits. The entire process, comprising data generation, quality control, statistical analyses, interpretation of association results, and linking results to biology should be designed and executed to minimize the generation of false-positive and false-negative associations and misleading links to biological processes. This review aims to provide general guidelines for data analysis that address data quality control, association tests, adjustment for population stratification, and significance evaluation to improve the reliability of conclusions. We also provide guidance on post-GWAS strategy and the interpretation of results. These guidelines are tailored to dairy cattle, which are characterized by long-range linkage disequilibrium, large half-sib families, and routinely collected phenotypes, requiring different approaches than those applied in human GWAS. We discuss common limitations and challenges that have been overlooked in the analysis and interpretation of GWAS to identify candidate sequence variants in dairy cattle.

Identification of Candidate Genes and Functional Pathways Associated with Body Size Traits in Chinese Holstein Cattle Based on GWAS Analysis

Body size is one of the most economically important traits of dairy cattle, as it is significantly associated with cow longevity, production, health, fertility, and environmental adaptation. The identification and application of genetic variants using a novel genetic approach, such as genome-wide association studies (GWASs), may give more insights into the genetic architecture of complex traits. The identification of genes, single nucleotide polymorphisms (SNPs), and pathways associated with the body size traits may offer a contribution to genomic selection and long-term planning for selection in dairy cows. In this study, we performed GWAS analysis to identify the genetic markers and genes associated with four body size traits (body height, body depth, chest width, and angularity) in 1000 Chinese Holstein cows. We performed SNPs genotyping in 1000 individuals, based on the GeneSeek Genomic Profiler Bovine 100 K. In total, we identified 11 significant SNPs in association with body size traits at the threshold of Bonferroni correction (5.90 × 10−7) using the fixed and random model circulating probability unification (FarmCPU) model. Several genes within 200 kb distances (upstream or downstream) of the significant SNPs were identified as candidate genes, including MYH15, KHDRBS3, AIP, DCC, SQOR, and UBAP1L. Moreover, genes within 200 kb of the identified SNPs were significantly enriched (p ≤ 0.05) in 25 Gene Ontology terms and five Kyoto Encyclopedia of Genes and Genomes pathways. We anticipate that these results provide a foundation for understanding the genetic architecture of body size traits. They will also contribute to breeding programs and genomic selection work on Chinese Holstein cattle.

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.

Direct genetic effects, maternal genetic effects, and maternal genetic sensitivity on prenatal heat stress for calf diseases and corresponding genomic loci in German Holsteins

The aim of this study was to infer the effects of heat stress (HS) of dams during late gestation on direct and maternal genetic parameters for pneumonia (PNEU, 112,563 observations), diarrhea (DIAR, 176,904 observations), and omphalitis (OMPH, 176,872 observations) in Holstein calves kept in large-scale co-operator herds. The genotype dataset included 41,135 SNPs from 19,247 male and female cattle. Temperature-humidity indices (THI) during the last 8 wk of pregnancy were calculated, using the climate data from the nearest public weather station for each herd. Heat load effects were considered for average weekly THI larger than 60. Phenotypically, regression coefficients of calf diseases on prenatal THI during the last 8 wk of gestation were estimated in 8 consecutive runs. The strongest detrimental effects of prenatal HS on PNEU and DIAR were identified for the last week of pregnancy (wk 1). Thus, only wk 1 was considered in ongoing genetic and genomic analyses. In an advanced model considering prenatal HS, random regression coefficients on THI in wk 1 nested within maternal genetic effects (maternal slope effects for heat load) were considered as parameters to infer maternal sensitivity in response to prenatal THI alterations. Direct heritabilities from the advanced model ranged from 0.10 (THI 60) to 0.08 (THI 74) for PNEU and were close to 0.16 for DIAR. Maternal heritabilities for PNEU increased from 0.03 to 0.10 along the THI gradient. For DIAR, the maternal heritability was largest (0.07) at the minimum THI (THI = 60) and decreased to 0.05 at THI 74. Genetic correlations smaller than 0.80 for PNEU and DIAR recorded at THI 60 with corresponding diseases at THI 74 indicated genotype by climate interactions for maternal genetic effects. Genome-wide associations studies were performed using de-regressed proofs of genotyped sires for direct genetic, maternal genetic, and maternal slope effects. Thirty suggestive and 2 significant SNPs were identified from the GWAS. Forty-three genes located close to the suggestive SNPs (±100 kb) were annotated as potential candidate genes. Three biological processes were inferred on the basis of the these genes, addressing the negative regulation of the viral life cycle, innate immune response, and protein ubiquitination. Hence, the genetics of prenatal heat stress mechanisms are associated with immune physiology and disease resistance mechanisms.

Estimation of direct and maternal genetic effects and annotation of potential candidate genes for weight and meat quality traits in a genotyped outdoor dual-purpose cattle breed

With regard to potential applications of genomic selection in small numbered breeds, we evaluated genomic models and focused on potential candidate gene annotations for weight and meat quality traits in the local Rotes Höhenvieh (RHV) breed. Traits included 6,003 birth weights (BWT), 5,719 200 d-weights (200dw), 4,594 365 d-weights (365dw), and 547 records for intramuscular fat content (IMF). A total of 581,304 SNP from 370 genotyped cattle with phenotypic records were included in genomic analyses. Model evaluations focused on single- and multiple-trait models with direct and with direct and maternal genetic effects. Genetic relationship matrices were based on pedigree (A-matrix), SNP markers (G-matrix), or both (H-matrix). Genome-wide association studies (GWASs) were carried out using linear mixed models to identify potential candidate genes for the traits of interest. De-regressed proofs (DRP) for direct and maternal genetic components were used as pseudo-phenotypes in the GWAS. Accuracies of direct breeding values were higher from models based on G or on H compared to A. Highest accuracies (> 0.89) were obtained for IMF with multiple-trait models using the G-matrix. Direct heritabilities with maternal genetic effects ranged from 0.62 to 0.66 for BWT, from 0.45 to 0.55 for 200dW, from 0.40 to 0.44 for 365dW, and from 0.48 to 0.75 for IMF. Maternal heritabilities for BWT, 200dW, and 365dW were in a narrow range from 0.21 to 0.24, 0.24 to 0.27, and 0.21 to 0.25, respectively, and from 0.25 to 0.65 for IMF. Direct genetic correlations among body weight traits were positive and favorable, and very similar from different models but showed a stronger variation with 0.31 (A), −0.13 (G), and 0.45 (H) between BWT and IMF. In gene annotations, we identified 6, 3, 1, and 6 potential candidate genes for direct genetic effect on BWT, 200dW, 365dW, and IMF traits, respectively. Regarding maternal genetic effects, four (SHROOM3, ZNF609, PECAM1, and TEX2) and two (TMEM182 and SEC11A) genes were detected as potential candidate genes for BWT and 365dW, respectively. Potential candidate genes for maternal effect on IMF were GRHL2, FGA, FGB, and CTNNA3. As the most important finding from a practical breeding perspective, a small number of genotyped RHV cattle enabled accurate breeding values for high heritability IMF.

Genome-wide associations for heat stress response suggest potential candidate genes underlying milk fatty acid composition in dairy cattle

Contents of milk fatty acids (FA) display remarkable alterations along climatic gradients. Detecting candidate genes underlying such alterations might be beneficial for the exploration of climate sensitivity in dairy cattle. Consequently, we aimed on the definition of FA heat stress indicators, considering FA breeding values in response to temperature-humidity index (THI) alterations. Indicators were used in GWAS, in ongoing gene annotations and for the estimation of chromosome-wide variance components. The phenotypic data set consisted of 39,600 test-day milk FA records from 5,757 first-lactation Holstein dairy cows kept in 16 large-scale German cooperator herds. The FA traits were C18:0, polyunsaturated fatty acids (PUFA), saturated fatty acids (SFA), and unsaturated fatty acids (UFA). After genotype quality control, 40,523 SNP markers from 3,266 cows and 930 sires were considered. Meteorological data from the weather station in closest herd distance were used for the calculation of maximum hourly daily THI, which were allocated to 10 different THI classes. The same FA from 3 stages of lactation were considered as different, but genetically correlated traits. Consequently, a 3-trait reaction norm model was used to estimate genetic parameters and breeding values for FA along THI classes, considering either pedigree (A) or genomic (G) relationship matrices. De-regressed proofs and genomic estimated breeding values at the intermediate THI class 5 and at the extreme THI class 10 were used as pseudophenotypes in ongoing genomic analyses for thermoneutral (TNC) and heat stress conditions (HSC), respectively. The differences in de-regressed proofs and in genomic estimated breeding values from both THI classes were pseudophenotypes for heat stress response (HSR). Genetic correlations between the same FA under TNC and HSC were smallest in the first lactation stage and ranged from 0.20 for PUFA to 0.87 for SFA when modeling with the A matrix, and from 0.35 for UFA to 0.86 for SFA when modeling with the G matrix. In the first lactation stage, larger additive genetic variances under HSC compared with TNC indicate climate sensitivity for C18:0, PUFA, and UFA. Climate sensitivity was also reflected by pronounced chromosome-wide genetic variances for HSR of PUFA and UFA in the first stage of lactation. For all FA under TNC, HSC, and HSR, quite large genetic variance proportions were explained by BTA14. In GWAS, 30 SNP (within or close to 38 potential candidate genes) overlapped for HSR of the different FA. One unique potential candidate gene (AMFR) was detected for HSR of PUFA, 15 for HSR of SFA (ADGRB1, DENND3, DUSP16, EFR3A, EMP1, ENSBTAG00000003838, EPS8, MGP, PIK3C2G, STYK1, TMEM71, GSG1, SMARCE1, CCDC57, and FASN) and 3 for HSR of UFA (ENSBTAG00000048091, PAEP, and EPPK1). The identified unique genes play key roles in milk FA synthesis and are associated with disease resistance in dairy cattle. The results suggest consideration of FA in combination with climatic responses when inferring genetic mechanisms of heat stress in dairy cows.

Genome Wide Association Study of Beef Traits in Local Alpine Breed Reveals the Diversity of the Pathways Involved and the Role of Time Stratification

Knowledge of the genetic architecture of key growth and beef traits in livestock species has greatly improved worldwide thanks to genome-wide association studies (GWAS), which allow to link target phenotypes to Single Nucleotide Polymorphisms (SNPs) across the genome. Local dual-purpose breeds have rarely been the focus of such studies; recently, however, their value as a possible alternative to intensively farmed breeds has become clear, especially for their greater adaptability to environmental change and potential for survival in less productive areas. We performed single-step GWAS and post-GWAS analysis for body weight (BW), average daily gain (ADG), carcass fleshiness (CF) and dressing percentage (DP) in 1,690 individuals of local alpine cattle breed, Rendena. This breed is typical of alpine pastures, with a marked dual-purpose attitude and good genetic diversity. Moreover, we considered two of the target phenotypes (BW and ADG) at different times in the individuals' life, a potentially important aspect in the study of the traits' genetic architecture. We identified 8 significant and 47 suggestively associated SNPs, located in 14 autosomal chromosomes (BTA). Among the strongest signals, 3 significant and 16 suggestive SNPs were associated with ADG and were located on BTA10 (50-60 Mb), while the hotspot associated with CF and DP was on BTA18 (55-62 MB). Among the significant SNPs some were mapped within genes, such as SLC12A1, CGNL1, PRTG (ADG), LOC513941 (CF), NLRP2 (CF and DP), CDC155 (DP). Pathway analysis showed great diversity in the biological pathways linked to the different traits; several were associated with neurogenesis and synaptic transmission, but actin-related and transmembrane transport pathways were also represented. Time-stratification highlighted how the genetic architectures of the same traits were markedly different between different ages. The results from our GWAS of beef traits in Rendena led to the detection of a variety of genes both well-known and novel. We argue that our results show that expanding genomic research to local breeds can reveal hitherto undetected genetic architectures in livestock worldwide. This could greatly help efforts to map genomic complexity of the traits of interest and to make appropriate breeding decisions.

Genetic and Genome-Wide Association Analysis of Yearling Weight Gain in Israel Holstein Dairy Calves

Yearling weight gain in male and female Israeli Holstein calves, defined as 365 × ((weight - 35)/age at weight) + 35, was analyzed from 814,729 records on 368,255 animals from 740 herds recorded between 1994 and 2021. The variance components were calculated based on valid records from 2008 through 2017 for each sex separately and both sexes jointly by a single-trait individual animal model analysis, which accounted for repeat records on animals. The analysis model also included the square root, linear, and quadratic effects of age at weight. Heritability and repeatability were 0.35 and 0.71 in the analysis of both sexes and similar in the single sex analyses. The regression of yearling weight gain on birth date in the complete data set was -0.96 kg/year. The complete data set was also analyzed by the same model as the variance component analysis, including both sexes and accounting for differing variance components for each sex. The genetic trend for yearling weight gain, including both sexes, was 1.02 kg/year. Genetic evaluations for yearling weight gain was positively correlated with genetic evaluations for milk, fat, protein production, and cow survival but negatively correlated with female fertility. Yearling weight gain was also correlated with the direct effect on dystocia, and increased yearling weight gain resulted in greater frequency of dystocia. Of the 1749 Israeli Holstein bulls genotyped with reliabilities >50%, 1445 had genetic evaluations. As genotyping of these bulls was performed using several single nucleotide polymorhphism (SNP) chip platforms, we included only those markers that were genotyped in >90% of the tested cohort. A total of 40,498 SNPs were retained. More than 400 markers had significant effects after permutation and correction for multiple testing (p_nominal < 1 × 10^-8). Considering all SNPs simultaneously, 0.69 of variance among the sires' transmitting ability was explained. There were 24 markers with coefficients of determination for yearling weight gain >0.04. One marker, BTA-75458-no-rs on chromosome 5, explained ≈6% of the variance among the estimated breeding values for yearling weight gain. ARS-BFGL-NGS-39379 had the fifth largest coefficient of determination in the current study and was also found to have a significant effect on weight at an age of 13-14 months in a previous study on Holsteins. Significant genomic effects on yearling weight gain were mainly associated with milk production quantitative trait loci, specifically with kappa casein metabolism.

Heat stress during late pregnancy and postpartum influences genetic parameter estimates for birth weight and weight gain in dual-purpose cattle offspring generations

Impact of direct heat stress (HS) on genetic parameter estimates, i.e., HS close to the trait recording date, was verified in several previous studies conducted in dairy and beef cattle populations. The aim of the present study was to analyze the impact of time-lagged HS at different recording periods during late pregnancy (a.p.) and postpartum (p.p.) on genetic parameter estimates for birth weight (BWT) and weight gain traits (200 d- and 365 d-weight gain (200dg, 365dg)) in offspring of the dual-purpose cattle breed "Rotes Höhenvieh" (RHV). Furthermore, we estimated genetic correlations within traits across time-lagged climatic indicators, in order to proof possible genotype by environment interactions (G×E). Trait recording included 5,434 observations for BWT, 3,679 observations for 200dg and 2,998 observations for 365dg. Time-lagged climatic descriptors were classes for the mean temperature humidity index (mTHI) and number of HS days (nHS) from the following periods: 7 d-period a.p. (BWT), 56 d-period a.p., and 56 d-period p.p. (200dg and 365dg). Genetic parameters were estimated via 2-trait animal models, i.e., defining the same trait in different climatic environments as different traits. Genetic variances and heritabilities for all traits increased with increasing mTHI- and nHS-classes for all recording periods, indicating pronounced genetic differentiation with regard to time-lagged in utero HS and HS directly after birth. Similarly, in low mTHI- and nHS-classes indicating cold stress, genetic variances, and heritabilities were larger than for temperate climates. Genetic correlations substantially smaller than 0.80 indicating G × E were observed when considering same traits from mTHI- and nHS-classes in greater distance. Estimated breeding values (EBV) of the 10 most influential sires with the largest number of offspring records fluctuated across mTHI- and nHS-classes. Correlations between sire EBV for same traits from distant climatic classes confirmed the genetic correlation estimates. Sires displaying stable EBV with climatic alterations were also identified. Selection of those sires might contribute to improved robustness in the RHV outdoor population genetically.

Single-Nucleotide Polymorphism Variations Associated With Specific Genes Putatively Identified Enhanced Genetic Predisposition for 305-Day Milk Yield in the Girolando Crossbreed

Milk production phenotypes are the main focus of genetic selection in dairy herds, and although there are many genes identified as related to the biology of these traits in pure breeds, little is known about crossbreed animals. This study aimed to identify potential genes associated with the 305-day milk yield in 337 crossbreed Gir × Holstein (Girolando) animals. Milk production records were genotyped for 45,613 single-nucleotide polymorphisms (SNPs). This dataset was used for a genome-wide association study (GWAS) using the 305-day milk yield adjusted for the fixed effects of herd and year and linear and quadratic effects of age at calving (in days) and calving factor averaged per animal. Genes within the significant SNPs were retrieved from the Bos taurus ARS-UCD1.2 assembly (bosTau9) for gene ontology analysis. In summary, the GWAS identified 52 SNPs associated [p ≤ 10^–4, false discovery rate (FDR) = 8.77%] with milk production, including NUB1 and SLC24A2, which were previously described as related to milk production traits in cattle. The results suggest that SNPs associated mainly with NUB1 and SLC24A2 could be useful to understand milk production in Girolando and used as predictive markers for selecting genetic predisposition for milk yield in Girolando.

Maternal and direct genetic parameters for tail length, tail lesions, and growth traits in pigs

Tail length and tail lesions are the major triggers for tail biting in pigs. Against this background, 2 datasets were analyzed to estimate genetic parameters for tail characteristics and growth traits. Dataset 1 considered measurements for trait tail length (T-LEN) and for the growth traits birth weight (BW), weaning weight (WW), postweaning weight (PWW), and average daily gain (ADG) from 9,348 piglets. Piglets were born in the period from 2015 to 2018 and kept on the university Gießen research station. Dataset 2 included 4,943 binary observations from 1,648 pigs from the birth years 2016 to 2019 for tail lesions (T-LES) as indicators for nail necrosis, tail abnormalities, or tail biting. T-LES were recorded at 30 ± 7 d after entry for rearing (T-Les-1), at 50 ± 7 d after entry for rearing (end of the rearing period, T-LES-2), and 130 ± 20 d after entry for rearing (end of fattening period, T-LES-3). Genetic statistical model evaluation for dataset 1 based on Akaike's information criterion and likelihood ration tests suggested multiple-trait animal models considering covariances between direct and maternal genetic effects. The direct heritability for T-LEN was 0.42 (±0.03), indicating the potential for genetic selection on short tails. The maternal genetic heritability for T-LEN was 0.05 (±0.04), indicating the influence of uterine characteristics on morphological traits. The negative correlation between direct and maternal effects for T-LEN of -0.35 (±0.13), as well as the antagonistic relationships (i.e., positive direct genetic correlations in the range from 0.03 to 0.40) between T-LEN with the growth traits BW, WW, PWW, and ADG, complicate selection strategies and breeding goal definitions. The correlations between direct effects for T-LEN and maternal effects for breeding goal traits, and vice versa, were positive but associated with a quite large SE. The heritability for T-LES when considering the 3 repeated measurements was 0.23 (±0.04) from the linear (repeatability of 0.30) and 0.21 (±0.06; repeatability of 0.29) from the threshold model. The breeding value correlations between T-LES-3 with breeding values from the repeatability models were quite large (0.74 to 0.90), suggesting trait lesion recording at the end of the rearing period. To understand all genetic mechanisms in detail, ongoing studies are focusing on association analyses between T-LEN and T-LES, and the identification of tail biting from an actor's perspective.

Population demographic history and population structure for Pakistani Nili-Ravi breeding bulls based on SNP genotyping to identify genomic regions associated with male effects for milk yield and body weight

The domestic Nili-Ravi water buffalo (Bubalus bubalis) is the best dairy animal contributing 68% to total milk production in Pakistan. In this study, we identified genome-wide single nucleotide polymorphisms (SNPs) to estimate various population genetic parameters such as diversity, pairwise population differentiation, linkage disequilibrium (LD) distribution and for genome-wide association study for milk yield and body weight traits in the Nili-Ravi dairy bulls that they may pass on to their daughters who are retained for milking purposes. The genotyping by sequencing approach revealed 13,039 reference genome-anchored SNPs with minor allele frequency of 0.05 among 167 buffalos. Population structure analysis revealed that the bulls were grouped into two clusters (K = 2), which indicates the presence of two different lineages in the Pakistani Nili-Ravi water buffalo population, and we showed the extent of admixture of these two lineages in our bull collection. LD analysis revealed 4169 significant SNP associations, with an average LD decay of 90 kb for these buffalo genome. Genome-wide association study involved a multi-locus mixed linear model for milk yield and body weight to identify genome-wide male effects. Our study further illustrates the utility of the genotyping by sequencing approach for identifying genomic regions to uncover additional demographic complexity and to improve the complex dairy traits of the Pakistani Nili-Ravi water buffalo population that would provide the lot of economic benefits to dairy industry.

GWAS-Based Identification of New Loci for Milk Yield, Fat, and Protein in Holstein Cattle

High-yield and high-quality of milk are the primary goals of dairy production. Understanding the genetic architecture underlying these milk-related traits is beneficial so that genetic variants can be targeted toward the genetic improvement. In this study, we measured five milk production and quality traits in Holstein cattle population from China. These traits included milk yield, fat, and protein. We used the estimated breeding values as dependent variables to conduct the genome-wide association studies (GWAS). Breeding values were estimated through pedigree relationships by using a linear mixed model. Genotyping was carried out on the individuals with phenotypes by using the Illumina BovineSNP150 BeadChip. The association analyses were conducted by using the fixed and random model Circulating Probability Unification (FarmCPU) method. A total of ten single-nucleotide polymorphisms (SNPs) were detected above the genome-wide significant threshold (p < 4.0 × 10-7), including six located in previously reported quantitative traits locus (QTL) regions. We found eight candidate genes within distances of 120 kb upstream or downstream to the associated SNPs. The study not only identified the effect of DGAT1 gene on milk fat and protein, but also discovered novel genetic loci and candidate genes related to milk traits. These novel genetic loci would be an important basis for molecular breeding in dairy cattle.

Estimating Conformational Traits in Dairy Cattle With DeepAPS: A Two-Step Deep Learning Automated Phenotyping and Segmentation Approach

Assessing conformation features in an accurate and rapid manner remains a challenge in the dairy industry. While recent developments in computer vision has greatly improved automated background removal, these methods have not been fully translated to biological studies. Here, we present a composite method (DeepAPS) that combines two readily available algorithms in order to create a precise mask for an animal image. This method performs accurately when compared with manual classification of proportion of coat color with an adjusted R² = 0.926. Using the output mask, we are able to automatically extract useful phenotypic information for 14 additional morphological features. Using pedigree and image information from a web catalog (www.semex.com), we estimated high heritabilities (ranging from h² = 0.18–0.82), indicating that meaningful biological information has been extracted automatically from imaging data. This method can be applied to other datasets and requires only a minimal number of image annotations (∼50) to train this partially supervised machine-learning approach. DeepAPS allows for the rapid and accurate quantification of multiple phenotypic measurements while minimizing study cost. The pipeline is available at https://github.com/lauzingaretti/deepaps.

Weighted Single-Step Genome-Wide Association Study for Growth Traits in Chinese Simmental Beef Cattle

Improving the genetic process of growth traits is one of the major goals in the beef cattle industry, as it can increase meat production and reduce the cost of raising animals. Although several quantitative trait loci affecting growth traits in beef cattle have been identified, the genetic architecture of these economically important traits remains elusive. This study aims to map single nucleotide polymorphisms (SNPs) and genes associated with birth weight (BW), yearling weight (YW), average daily gain from birth to yearling (BYADG), and body weight at the age of 18 months (18MW) in a Chinese Simmental beef cattle population using a weighted, single-step, genome-wide association study (wssGWAS). Phenotypic and pedigree data from 6022 animals and genotypes from 744 animals (596,297 SNPs) were used for an association analysis. The results showed that 66 genomic windows explained 1.01-20.15% of the genetic variance for the four examined traits, together with the genes near the top SNP within each window. Furthermore, the identified genomic windows (>1%) explained 50.56%, 57.71%, 61.78%, and 37.82% of the genetic variances for BW, YW, BYADG, and 18MW, respectively. Genes with potential functions in muscle development and regulation of cell growth were highlighted as candidates for growth traits in Simmental cattle (SQOR and TBCB for BW, MYH10 for YW, RLF for BYADG, and ARHGAP31 for 18MW). Moreover, we found 40 SNPs that had not previously been identified as being associated with growth traits in cattle. These findings will further advance our understanding of the genetic basis for growth traits and will be useful for the molecular breeding of BW, YW, BYADG, and 18MW in the context of genomic selection in beef cattle.

Multi-breed genome-wide association studies across countries for electronically recorded behavior traits in local dual-purpose cows

Basic bovine behavior is a crucial parameter influencing cattle domestication. In addition, behavior has an impact on cattle productivity, welfare and adaptation. The aim of the present study was to infer quantitative genetic and genomic mechanisms contributing to natural dual-purpose cow behavior in grazing systems. In this regard, we genotyped five dual-purpose breeds for a dense SNP marker panel from four different European countries. All cows from the across-country study were equipped with the same electronic recording devices. In this regard, we analyzed 97,049 longitudinal sensor behavior observations from 319 local dual-purpose cows for rumination, feeding, basic activity, high active, not active and ear temperature. According to the specific sensor behaviors and following a welfare protocol, we computed two different welfare indices. For genomic breed characterizations and multi-breed genome-wide association studies, sensor traits and test-day production records were merged with 35,826 SNP markers per cow. For the estimation of variance components, we used the pedigree relationship matrix and a combined similarity matrix that simultaneously included both pedigree and genotypes. Heritabilities for feeding, high active and not active were in a moderate range from 0.16 to 0.20. Estimates were very similar from both relationship matrix-modeling approaches and had quite small standard errors. Heritabilities for the remaining sensor traits (feeding, basic activity, ear temperature) and welfare indices were lower than 0.09. Five significant SNPs on chromosomes 11, 17, 27 and 29 were associated with rumination, and two different SNPs significantly influenced the sensor traits “not active” (chromosome 13) and “feeding” (chromosome 23). Gene annotation analyses inferred 22 potential candidate genes with a false discovery rate lower than 20%, mostly associated with rumination (13 genes) and feeding (8 genes). Mendelian randomization based on genomic variants (i.e., the instrumental variables) was used to infer causal inference between an exposure and an outcome. Significant regression coefficients among behavior traits indicate that all specific behavioral mechanisms contribute to similar physiological processes. The regression coefficients of rumination and feeding on milk yield were 0.10 kg/% and 0.12 kg/%, respectively, indicating their positive influence on dual-purpose cow productivity. Genomically, an improved welfare behavior of grazing cattle, i.e., a higher score for welfare indices, was significantly associated with increased fat and protein percentages.