Genotype-by-environment interaction of fertility traits in Danish Holstein cattle using a single-step genomic reaction norm model

Review: Opportunities and challenges for the genetic selection of dairy calf disease traits

Interest in dairy cow health continues to grow as we better understand health's relationship with production potential and animal welfare. Over the past decade, efforts have been made to incorporate health traits into national genetic evaluations. However, they have focused on the mature cow, with calf health largely being neglected. Diarrhoea and respiratory disease comprise the main illnesses with regard to calf health. Conventional methods to control calf disease involve early separation of calves from the dam and housing calves individually. However, public concern regarding these methods, and growing evidence that these methods may negatively impact calf development, mean the dairy industry may move away from these practices. Genetic selection may be a promising tool to address these major disease issues. In this review, we examined current literature for enhancing calf health through genetics and discussed alternative approaches to improve calf health via the use of epidemiological modelling approaches, and the potential of indirectly selecting for improved calf health through improving colostrum quality. Heritability estimates on the observed scale for diarrhoea ranged from 0.03 to 0.20, while for respiratory disease, estimates ranged from 0.02 to 0.24. The breadth in these ranges is due, at least in part, to differences in disease prevalence, population structure, data editing and models, as well as data collection practices, which should be all considered when comparing literature values. Incorporation of epidemiological theory into quantitative genetics provides an opportunity to better determine the level of genetic variation in disease traits, as it accounts for disease transmission among contemporaries. Colostrum intake is a major determinant of whether a calf develops either respiratory disease or diarrhoea. Colostrum traits have the advantage of being measured and reported on a continuous scale, which removes the issues classically associated with binary disease traits. Overall, genetic selection for improved calf health is feasible. However, to ensure the maximum response, first steps by any industry members should focus efforts on standardising recording practices and encouragement of uploading information to genetic evaluation centres through herd management software, as high-quality phenotypes are the backbone of any successful breeding programme.

Review: Ruminant heat-stress terminology

With increasing climate variability, there is a rise in the exposure to, and incidence of, ruminant heat stress (HS), increasing the requirement for focused research. As such, precise terminology is crucial to maintain effective communication and knowledge advancement. Despite this, several key terms are currently defined inconsistently, leading to confusion and misinterpretation. This paper examines the historical and contemporary use of the terms 'resistance', 'tolerance', 'resilience', and 'susceptibility' across various disciplines, revealing significant ambiguities that hinder both research and practice. Through this comprehensive review, we propose new definitions for each term as they are used relating to HS, with a focus on ruminant production. Proposed definitions align with current scientific understanding, providing a robust framework for future research and application. As further research is conducted, we hope these definitions can be improved through the inclusion of quantitative measures which align with these classifications. This present review provides definition clarity for common heat abatement terminology, enabling consistency and from this, progress in the field to ameliorate HS for ruminants.

CFAP58 is involved in the sperm head shaping and flagellogenesis of cattle and mice

CFAP58 is a testis-enriched gene that plays an important role in the sperm flagellogenesis of humans and mice. However, the effect of CFAP58 on bull semen quality and the underlying molecular mechanisms involved in spermatogenesis remain unknown. Here, we identified two single-nucleotide polymorphisms (rs110610797, A>G and rs133760846, G>T) and one indel (g.-1811_ g.-1810 ins147bp) in the promoter of CFAP58 that were significantly associated with semen quality of bulls, including sperm deformity rate and ejaculate volume. Moreover, by generating gene knockout mice, we found for the first time that the loss of Cfap58 not only causes severe defects in the sperm tail, but also affects the manchette structure, resulting in abnormal sperm head shaping. Cfap58 deficiency causes an increase in spermatozoa apoptosis. Further experiments confirmed that CFAP58 interacts with IFT88 and CCDC42. Moreover, it may be a transported cargo protein that plays a role in stabilizing other cargo proteins, such as CCDC42, in the intra-manchette transport/intra-flagellar transport pathway. Collectively, our findings reveal that CFAP58 is required for spermatogenesis and provide genetic markers for evaluating semen quality in cattle.

Uncovering genetic parameters and environmental influences on fertility, milk production, and quality in autochthonous Reggiana cattle

Reggiana is a local cattle breed from northern Italy known for its rusticity and profitability, due to the production of branded Parmigiano Reggiano cheese. To ensure the persistence of such profitability in the long term, an adequate breeding program is required. To this aim, in the present study we estimate the genetic parameters of the main productive and reproductive traits, and we evaluate the effect of genotype by environment interaction (GxE) on these traits using 2 environmental covariates: (1) productivity and (2) temperature-humidity index (THI). Milk, fat, protein, and casein yield were considered as daily production traits, whereas protein, fat, casein percentage, casein index, and somatic cell score were considered as milk quality traits. Finally, reproductive traits such as the number of inseminations, days open, calving interval, and calving-to-first-insemination interval were evaluated. Reggiana cattle produce an average of 19 kg of milk per day with 3.7% fat and 3.4% protein content and have excellent fertility parameters. Compared with other breeds, they have slightly lower heritability for production and quality for production traits (e.g., 0.12 [0.09; 0.15] for milk yield), but similar heritability for fertility traits. Milk, protein, and fat daily yields are highly correlated but negatively correlated with the percentage of protein, fat, and casein, whereas fertility traits have an unfavorable genetic correlation with daily production traits. When considering productivity, a consistent amount of variability due to GxE was observed for all daily production traits, somatic cell count, and casein index. A modest amount of GxE was observed for fertility parameters, while the percentage of solid content showed almost no GxE effect. A similar situation occurred when considering the THI, but no GxE interaction was observed for reproduction traits. In conclusion, this study provides useful information for the implementation of accurate selection plans in this local breed, accounting for environmental plasticity measured through the consistent GxE interaction observed.

A simulation-based assessment of the efficiency of QTL mapping under environment and genotype x environment interaction effects

The objective of this simulation-based study was to assess how genes, environments, and genotype x environment (GxE) interaction affect the quantitative trait loci (QTL) mapping efficiency. The simulation software performed 50 samplings of 300 recombinant inbred lines (RILs) from a F2, which were assessed in six environments. The RILs were genotyped for 977 single nucleotide polymorphisms (SNP) and phenotyped for grain yield. The average SNP density was 2 cM. We defined six QTLs and 190 minor genes. The trait heritability ranged from 30 to 80%. We fitted the single QTL model and the multiple QTL model on multiple phenotypes. The environment and complex GxE interaction effects led to a low correlation between the QTL heritability and power. The single- and across-environment analyses allowed all QTLs be declared, with an average power of 28 to 100%. In the across-environment analysis, five QTLs showed average power in the range 46 to 82%. Both models provided a good control of the false positive rate (6%, on average) and a precise localization of the QTLs (bias of 2 cM, on average). The QTL power in each environment has a high positive correlation with the range between QTL genotypes for the sum of the additive, environment, and GxE interaction effects (0.76 to 0.96). The uncertainty about the magnitude and sign of the environment and GxE interaction effects makes QTL mapping in multi-environment trials unpredictable. Unfortunately, this uncertainty has no solution because the geneticist has no control over the magnitude and sign of the environment and GxE interaction effects. However, the single- and across-environment analyses are efficient even under a low correlation between QTL heritability and power.

Organic Milk Production and Dairy Farming Constraints and Prospects under the Laws of the European Union

In recent years, there has been rapid development in organic farming. When choosing organic livestock products, consumers are guided by the conviction that animals are provided with the highest welfare standards and access to pasture. The purpose of this article was to trace the principles of organic farming prevailing in the EU with regard to milk production and cattle breeding. The principles of organic production are universal and their application is not limited to certified farms. Organic certification is intended to assure the consumer of the quality and method of production. Due to additional requirements imposed by law, organic cows are usually kept in better welfare conditions compared to conventional cattle, but this is not the rule. The altered taste and texture of organic milk and its products compared to conventional products mainly depends on the presence of pasture greens in the cows' diet. Therefore, milk from conventionally kept, pasture-grazed cows may have similar characteristics and composition. Organic farms tend to have lower milk yields compared to conventional farms due to the lower consumption of concentrate feed. In the future, it is expected that the proportion of land that is unsuitable for the production of crops for human consumption will increasingly be used for cow grazing.

Parent-offspring genotyped trios unravelling genomic regions with gametic and genotypic epistatic transmission bias on the cattle genome

Several biological mechanisms affecting the sperm and ova fertility and viability at developmental stages of the reproductive cycle resulted in observable transmission ratio distortion (i.e., deviation from Mendelian expectations). Gene-by-gene interactions (or epistasis) could also potentially cause specific transmission ratio distortion patterns at different loci as unfavorable allelic combinations are under-represented, exhibiting deviation from Mendelian proportions. Here, we aimed to detect pairs of loci with epistatic transmission ratio distortion using 283,817 parent-offspring genotyped trios (sire-dam-offspring) of Holstein cattle. Allelic and genotypic parameterization for epistatic transmission ratio distortion were developed and implemented to scan the whole genome. Different epistatic transmission ratio distortion patterns were observed. Using genotypic models, 7, 19 and 6 pairs of genomic regions were found with decisive evidence with additive-by-additive, additive-by-dominance/dominance-by-additive and dominance-by-dominance effects, respectively. Using the allelic transmission ratio distortion model, more insight was gained in understanding the penetrance of single-locus distortions, revealing 17 pairs of SNPs. Scanning for the depletion of individuals carrying pairs of homozygous genotypes for unlinked loci, revealed 56 pairs of SNPs with recessive epistatic transmission ratio distortion patterns. The maximum number of expected homozygous offspring, with none of them observed, was 23. Finally, in this study, we identified candidate genomic regions harboring epistatic interactions with potential biological implications in economically important traits, such as reproduction.

The effect of high-density genotypic data and different methods on joint genomic prediction: A case study in large white pigs

Joint genomic prediction (GP) is an attractive method to improve the accuracy of GP by combining information from multiple populations. However, many factors can negatively influence the accuracy of joint GP, such as differences in linkage disequilibrium phasing between single nucleotide polymorphisms (SNPs) and causal variants, minor allele frequencies and causal variants' effect sizes across different populations. The objective of this study was to investigate whether the imputed high-density genotype data can improve the accuracy of joint GP using genomic best linear unbiased prediction (GBLUP), single-step GBLUP (ssGBLUP), multi-trait GBLUP (MT-GBLUP) and GBLUP based on genomic relationship matrix considering heterogenous minor allele frequencies across different populations (wGBLUP). Three traits, including days taken to reach slaughter weight, backfat thickness and loin muscle area, were measured on 67 276 Large White pigs from two different populations, for which 3334 were genotyped by SNP array. The results showed that a combined population could substantially improve the accuracy of GP compared with a single-population GP, especially for the population with a smaller size. The imputed SNP data had no effect for single population GP but helped to yield higher accuracy than the medium-density array data for joint GP. Of the four methods, ssGLBUP performed the best, but the advantage of ssGBLUP decreased as more individuals were genotyped. In some cases, MT-GBLUP and wGBLUP performed better than GBLUP. In conclusion, our results confirmed that joint GP could be beneficial from imputed high-density genotype data, and the wGBLUP and MT-GBLUP methods are promising for joint GP in pig breeding.

Heteroscedastic Reaction Norm Models Improve the Assessment of Genotype by Environment Interaction for Growth, Reproductive, and Visual Score Traits in Nellore Cattle

The assessment of the presence of genotype by environment interaction (GxE) in beef cattle is very important in tropical countries with diverse climatic conditions and production systems. The present study aimed to assess the presence of GxE by using different reaction norm models for eleven traits related to growth, reproduction, and visual score in Nellore cattle. We studied five reaction norm models (RNM), fitting a linear model considering homoscedastic residual variance (RNM_homo), and four models considering heteroskedasticity, being linear (RNM_hete), quadratic (RNM_quad), linear spline (RNM_l-l), and quadratic spline (RNM_q-q). There was the presence of GxE for age at first calving (AFC), scrotal circumference (SC), weaning to yearling weight gain (WYG), and yearling weight (YW). The best models were RNM_l-l for YW and RNM_q-q for AFC, SC, and WYG. The heritability estimates for RNM_l-l ranged from 0.07 to 0.20, 0.42 to 0.61, 0.24 to 0.42, and 0.47 to 0.63 for AFC, SC, WYG, and YW, respectively. The heteroskedasticity in reaction norm models improves the assessment of the presence of GxE for YW, WYG, AFC, and SC. Additionally, the trajectories of reaction norms for these traits seem to be affected by a non-linear component, and selecting robust animals for these traits is an alternative to increase production and reduce environmental sensitivity.

G × EBLUP: A novel method for exploring genotype by environment interactions and genomic prediction

Genotype by environment (G × E) interaction is fundamental in the biology of complex traits and diseases. However, most of the existing methods for genomic prediction tend to ignore G × E interaction (GEI). In this study, we proposed the genomic prediction method G × EBLUP by considering GEI. Meanwhile, G × EBLUP can also detect the genome-wide single nucleotide polymorphisms (SNPs) subject to GEI. Using comprehensive simulations and analysis of real data from pigs and maize, we showed that G × EBLUP achieved higher efficiency in mapping GEI SNPs and higher prediction accuracy than the existing methods, and its superiority was more obvious when the GEI variance was large. For pig and maize real data, compared with GBLUP, G × EBLUP showed improvement by 3% in the prediction accuracy for backfat thickness, while our findings indicated that the trait of days to 100 kg of pig was not affected by GEI and G × EBLUP did not improve the accuracy of genomic prediction for the trait. A significant advantage was observed for G × EBLUP in maize; the prediction accuracy was improved by ∼5.0 and 7.7% for grain weight and water content, respectively. Furthermore, G × EBLUP was not influenced by the number of environment levels. It could determine a favourable environment using SNP Bayes factors for each environment, implying that it is a robust and useful method for market-specific animal and plant breeding. We proposed G × EBLUP, a novel method for the estimation of genomic breeding value by considering GEI. This method identified the genome-wide SNPs that were susceptible to GEI and yielded higher genomic prediction accuracies and lower mean squared error compared with the GBLUP method.

Genomic analysis of the slope of the reaction norm for body weight in Australian sheep

Background

Selection of livestock based on their robustness or sensitivity to environmental variation could help improve the efficiency of production systems, particularly in the light of climate change. Genetic variation in robustness arises from genotype-by-environment (G × E) interactions, with genotypes performing differently when animals are raised in contrasted environments. Understanding the nature of this genetic variation is essential to implement strategies to improve robustness. In this study, our aim was to explore the genetics of robustness in Australian sheep to different growth environments using linear reaction norm models (RNM), with post-weaning weight records of 22,513 lambs and 60 k single nucleotide polymorphisms (SNPs). The use of scale-corrected genomic estimated breeding values (GEBV) for the slope to account for scale-type G × E interactions was also investigated.

Results

Additive genetic variance was observed for the slope of the RNM, with genetic correlations between low- and high-growth environments indicating substantial re-ranking of genotypes (0.44–0.49). The genetic variance increased from low- to high-growth environments. The heritability of post-weaning body weight ranged from 0.28 to 0.39. The genetic correlation between intercept and slope of the reaction norm for post-weaning body weight was low to moderate when based on the estimated (co)variance components but was much higher when based on back-solved SNP effects. An initial analysis suggested that a region on chromosome 11 affected both the intercept and the slope, but when the GEBV for the slope were conditioned on the GEBV for the intercept to remove the effect of scale-type G × E interactions on SNP effects for robustness, a single genomic region on chromosome 7 was found to be associated with robustness. This region included genes previously associated with growth traits and disease susceptibility in livestock.

Conclusions

This study shows a significant genetic variation in the slope of RNM that could be used for selecting for increased robustness of sheep. Both scale-type and rank-type G × E interactions contributed to variation in the slope. The correction for scale effects of GEBV for the slope should be considered when analysing robustness using RNM. Overall, robustness appears to be a highly polygenic trait.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00734-6.

Candidate Genes in Bull Semen Production Traits: An Information Approach Review

Semen quality plays a crucial role in the successful implementation of breeding programs, especially where artificial insemination (AI) is practiced. Bulls with good semen traits have good fertility and can produce a volume of high semen per ejaculation. The aim of this review is to use an information approach to highlight candidate genes and their relation to bull semen production traits. The use of genome-wide association studies (GWAS) has been demonstrated to be successful in identifying genomic regions and individual variations associated with production traits. Studies have reported over 40 genes associated with semen traits using Illumina BeadChip single-nucleotide polymorphism (SNPs).

Genome-wide association study of beef bull semen attributes

Background

Cattle production is dependent upon fertility because it results in producing offspring to offset production costs. A number of semen attributes are believed to affect fertility and are frequently measured as part of routine breeding soundness exams or semen collection procedures. The objective of this study was to perform a single-step genome-wide association study (ssGWAS) for beef bull semen attributes. Beef bull fertility phenotypes including volume (VOL), concentration (CONC), number of spermatozoa (NSP), initial motility (IMot), post-thaw motility (PTMot), three-hour post-thaw motility (3HRPTMot), percentage of normal spermatozoa (%NORM), primary abnormalities (PRIM), and secondary abnormalities (SEC) were obtained from two artificial insemination (AI) centers. A total of 1819 Angus bulls with 50,624 collection records were used for ssGWAS. A five-generation pedigree was obtained from the American Angus Association and consisted of 6521 sires and 17,136 dams. Genotypes on 1163 bulls were also obtained from the American Angus Association and utilized in ssGWAS.

Results

A multi-trait animal model was used for the estimation of single nucleotide polymorphism (SNP) effects. Significant SNP were those with a -log10P-value threshold greater than 4.0. Volume, CONC, NSP, IMot, PTMot, 3HRPTMot, %NORM, PRIM, and SEC have five, three, six, seven, two, six, six, and two genome-wide significant SNP, respectively.

Conclusions

Several significant SNP were determined to be near or within quantitative trait loci (QTL) associated with beef bull semen attributes. In addition, genes associated with fertility were found to contain or be near the significant SNP found in the study. The results indicate there are regions of the genome that impact fertility, proving inclusion of genomic information into genetic evaluation should be advantageous for genetic improvement of male fertility traits.

Definition of Environmental Variables and Critical Periods to Evaluate Heat Tolerance in Large White Pigs Based on Single-Step Genomic Reaction Norms

Properly quantifying environmental heat stress (HS) is still a major challenge in livestock breeding programs, especially as adverse climatic events become more common. The definition of critical periods and climatic variables to be used as the environmental gradient is a key step for genetically evaluating heat tolerance (HTol). Therefore, the main objectives of this study were to define the best critical periods and environmental variables (ENV) to evaluate HT and estimate variance components for HT in Large White pigs. The traits included in this study were ultrasound backfat thickness (BFT), ultrasound muscle depth (MDP), piglet weaning weight (WW), off-test weight (OTW), interval between farrowing (IBF), total number of piglets born (TNB), number of piglets born alive (NBA), number of piglets born dead (NBD), number of piglets weaned (WN), and weaning to estrus interval (IWE). Seven climatic variables based on public weather station data were compared based on three criteria, including the following: (1) strongest G×E estimate as measured by the slope term, (2) ENV yielding the highest theoretical accuracy of the genomic estimated breeding values (GEBV), and (3) variable yielding the highest distribution of GEBV per ENV. Relative humidity (for BFT, MDP, NBD, WN, and WW) and maximum temperature (for OTW, TNB, NBA, IBF, and IWE) are the recommended ENV based on the analyzed criteria. The acute HS (average of 30 days before the measurement date) is the critical period recommended for OTW, BFT, and MDP in the studied population. For WN, WW, IBF, and IWE, a period ranging from 34 days prior to farrowing up to weaning is recommended. For TNB, NBA, and NBD, the critical period from 20 days prior to breeding up to 30 days into gestation is recommended. The genetic correlation values indicate that the traits were largely (WN, WW, IBF, and IWE), moderately (OTW, TNB, and NBA), or weakly (MDP, BFT, and NBD) affected by G×E interactions. This study provides relevant recommendations of critical periods and climatic gradients for several traits in order to evaluate HS in Large White pigs. These observations demonstrate that HT in Large White pigs is heritable, and genetic progress can be achieved through genetic and genomic selection.

Genetic and Genomic Analyses of Service Sire Effect on Female Reproductive Traits in Holstein Cattle

Fertility and reproductive performance are key drivers of dairy farm profitability. Hence, reproduction traits have been included in a large majority of worldwide dairy cattle selection indexes. The reproductive traits are lowly heritable but can be improved through direct genetic selection. However, most scientific studies and dairy cattle breeding programs have focused solely on the genetic effects of the dam (GED) on reproductive performance and, therefore, ignored the contribution of the service sire in the phenotypic outcomes. This study aimed to investigate the service sire effects on female reproductive traits in Holstein cattle from a genomic perspective. Genetic parameter estimation and genome-wide association studies (GWAS) were performed for the genetic effect of service sire (GESS) on conception rate (CR), 56-day non-return rate (NRR56), calving ease (CE), stillbirth (SB), and gestation length (GL). Our findings indicate that the additive genetic effects of both sire and dam contribute to the phenotypic variance of reproductive traits measured in females (0.0196 vs. 0.0109, 0.0237 vs. 0.0133, 0.0040 vs. 0.0289, 0.0782 vs. 0.0083, and 0.1024 vs. 0.1020 for GESS and GED heritability estimates for CR, NRR56, CE, SB, and GL, respectively), and these two genetic effects are positively correlated for SB (0.1394) and GL (0.7871). Interestingly, the breeding values for GESS on insemination success traits (CR and NRR56) are unfavorably and significantly correlated with some production, health, and type breeding values (ranging from -0.449 to 0.274), while the GESS values on calving traits (CE, SB, and GL) are usually favorably associated with those traits (ranging from -0.493 to 0.313). One hundred sixty-two significant single-nucleotide polymorphisms (SNPs) and their surrounding protein-coding genes were identified as significantly associated with GESS and GED, respectively. Six genes overlapped between GESS and GED for calving traits and 10 genes overlapped between GESS for success traits and calving traits. Our findings indicate the importance of considering the GESS when genetically evaluating the female reproductive traits in Holstein cattle.

Genome-Wide Associations for Microscopic Differential Somatic Cell Count and Specific Mastitis Pathogens in Holstein Cows in Compost-Bedded Pack and Cubicle Farming Systems

Simple Summary

New free walk housing systems such as compost-bedded pack barns might positively influence animal welfare. However, udder health can be negatively affected due to the microbial environment in the pack. Udder health depends on many factors, such as the environment, the feed, the pathogen species, and the genetic mechanisms of the cow’s immune system. For a more precise evaluation of udder health, we examined novel traits including specific mastitis pathogens and differential somatic cell fractions in milk. In order to identify possible candidate genes for udder health, a genome-wide association study, including single-nucleotide polymorphisms (SNP) by housing system interactions (compost-bedded pack barn and conventional cubicle barn), was performed. We identified two potential candidate genes for the interaction effect in relation to udder health. The identified potential candidate gene HEMK1 (HemK methyltransferase family member 1) is involved in immune system development, and CHL1 (cell adhesion molecule L1 like) has an immunosuppressive effect during stress conditions. The results suggest housing system-specific breeding strategies in order to improve udder health in compost-bedded pack and conventional cubicle barns.

Abstract

The aim of the present study was to detect significant SNP (single-nucleotide polymorphism) effects and to annotate potential candidate genes for novel udder health traits in two different farming systems. We focused on specific mastitis pathogens and differential somatic cell fractions from 2198 udder quarters of 537 genotyped Holstein Friesian cows. The farming systems comprised compost-bedded pack and conventional cubicle barns. We developed a computer algorithm for genome-wide association studies allowing the estimation of main SNP effects plus consideration of SNPs by farming system interactions. With regard to the main effect, 35 significant SNPs were detected on 14 different chromosomes for the cell fractions and the pathogens. Six SNPs were significant for the interaction effect with the farming system for most of the udder health traits. We inferred two possible candidate genes based on significant SNP interactions. HEMK1 plays a role in the development of the immune system, depending on environmental stressors. CHL1 is regulated in relation to stress level and influences immune system mechanisms. The significant interactions indicate that gene activity can fluctuate depending on environmental stressors. Phenotypically, the prevalence of mastitis indicators differed between systems, with a notably lower prevalence of minor bacterial indicators in compost systems.

Genotype-by-environment interactions for reproduction, body composition, and growth traits in maternal-line pigs based on single-step genomic reaction norms

BACKGROUND

There is an increasing need to account for genotype-by-environment (G × E) interactions in livestock breeding programs to improve productivity and animal welfare across environmental and management conditions. This is even more relevant for pigs because selection occurs in high-health nucleus farms, while commercial pigs are raised in more challenging environments. In this study, we used single-step homoscedastic and heteroscedastic genomic reaction norm models (RNM) to evaluate G × E interactions in Large White pigs, including 8686 genotyped animals, for reproduction (total number of piglets born, TNB; total number of piglets born alive, NBA; total number of piglets weaned, NW), growth (weaning weight, WW; off-test weight, OW), and body composition (ultrasound muscle depth, MD; ultrasound backfat thickness, BF) traits. Genetic parameter estimation and single-step genome-wide association studies (ssGWAS) were performed for each trait.

RESULTS

The average performance of contemporary groups (CG) was estimated and used as environmental gradient in the reaction norm analyses. We found that the need to consider heterogeneous residual variance in RNM models was trait dependent. Based on estimates of variance components of the RNM slope and of genetic correlations across environmental gradients, G × E interactions clearly existed for TNB and NBA, existed for WW but were of smaller magnitude, and were not detected for NW, OW, MD, and BF. Based on estimates of the genetic variance explained by the markers in sliding genomic windows in ssGWAS, several genomic regions were associated with the RNM slope for TNB, NBA, and WW, indicating specific biological mechanisms underlying environmental sensitivity, and dozens of novel candidate genes were identified. Our results also provided strong evidence that the X chromosome contributed to the intercept and slope of RNM for litter size traits in pigs.

CONCLUSIONS

We provide a comprehensive description of G × E interactions in Large White pigs for economically-relevant traits and identified important genomic regions and candidate genes associated with GxE interactions on several autosomes and the X chromosome. Implementation of these findings will contribute to more accurate genomic estimates of breeding values by considering G × E interactions, in order to genetically improve the environmental robustness of maternal-line pigs.

Genetic Features of Reproductive Traits in Bovine and Buffalo: Lessons From Bovine to Buffalo

Bovine and buffalo are important livestock species that have contributed to human lives for more than 1000 years. Improving fertility is very important to reduce the cost of production. In the current review, we classified reproductive traits into three categories: ovulation, breeding, and calving related traits. We systematically summarized the heritability estimates, molecular markers, and genomic selection (GS) for reproductive traits of bovine and buffalo. This review aimed to compile the heritability and genome-wide association studies (GWASs) related to reproductive traits in both bovine and buffalos and tried to highlight the possible disciplines which should benefit buffalo breeding. The estimates of heritability of reproductive traits ranged were from 0 to 0.57 and there were wide differences between the populations. For some specific traits, such as age of puberty (AOP) and calving difficulty (CD), the majority beef population presents relatively higher heritability than dairy cattle. Compared to bovine, genetic studies for buffalo reproductive traits are limited for age at first calving and calving interval traits. Several quantitative trait loci (QTLs), candidate genes, and SNPs associated with bovine reproductive traits were screened and identified by candidate gene methods and/or GWASs. The IGF1 and LEP pathways in addition to non-coding RNAs are highlighted due to their crucial relevance with reproductive traits. The distribution of QTLs related to various traits showed a great differences. Few GWAS have been performed so far on buffalo age at first calving, calving interval, and days open traits. In addition, we summarized the GS studies on bovine and buffalo reproductive traits and compared the accuracy between different reports. Taken together, GWAS and candidate gene approaches can help to understand the molecular genetic mechanisms of complex traits. Recently, GS has been used extensively and can be performed on multiple traits to improve the accuracy of prediction even for traits with low heritability, and can be combined with multi-omics for further analysis.

Genotype-by-environment interaction in Holstein heifer fertility traits using single-step genomic reaction norm models

Background

The effect of heat stress on livestock production is a worldwide issue. Animal performance is influenced by exposure to harsh environmental conditions potentially causing genotype-by-environment interactions (G × E), especially in highproducing animals. In this context, the main objectives of this study were to (1) detect the time periods in which heifer fertility traits are more sensitive to the exposure to high environmental temperature and/or humidity, (2) investigate G × E due to heat stress in heifer fertility traits, and, (3) identify genomic regions associated with heifer fertility and heat tolerance in Holstein cattle.

Results

Phenotypic records for three heifer fertility traits (i.e., age at first calving, interval from first to last service, and conception rate at the first service) were collected, from 2005 to 2018, for 56,998 Holstein heifers raised in 15 herds in the Beijing area (China). By integrating environmental data, including hourly air temperature and relative humidity, the critical periods in which the heifers are more sensitive to heat stress were located in more than 30 days before the first service for age at first calving and interval from first to last service, or 10 days before and less than 60 days after the first service for conception rate. Using reaction norm models, significant G × E was detected for all three traits regarding both environmental gradients, proportion of days exceeding heat threshold, and minimum temperature-humidity index. Through single-step genome-wide association studies, PLAG1, AMHR2, SP1, KRT8, KRT18, MLH1, and EOMES were suggested as candidate genes for heifer fertility. The genes HCRTR1, AGRP, PC, and GUCY1B1 are strong candidates for association with heat tolerance.

Conclusions

The critical periods in which the reproductive performance of heifers is more sensitive to heat stress are trait-dependent. Thus, detailed analysis should be conducted to determine this particular period for other fertility traits. The considerable magnitude of G × E and sire re-ranking indicates the necessity to consider G × E in dairy cattle breeding schemes. This will enable selection of more heat-tolerant animals with high reproductive efficiency under harsh climatic conditions. Lastly, the candidate genes identified to be linked with response to heat stress provide a better understanding of the underlying biological mechanisms of heat tolerance in dairy cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07496-3.

Climate sensitivity of milk production traits and milk fatty acids in genotyped Holstein dairy cows

The aim of this study was the evaluation of climate sensitivity via genomic reaction norm models [i.e., to infer cow milk production and milk fatty acid (FA) responses on temperature-humidity index (THI) alterations]. Test-day milk traits were recorded between 2010 and 2016 from 5,257 first-lactation genotyped Holstein dairy cows. The cows were kept in 16 large-scale cooperator herds, being daughters of 344 genotyped sires. The longitudinal data consisted of 47,789 test-day records for the production traits milk yield (MY), fat yield (FY), and protein yield (PY), and of 20,742 test-day records for 6 FA including C16:0, C18:0, saturated fatty acids (SFA), unsaturated fatty acids (UFA), monounsaturated fatty acids (MUFA), and polyunsaturated fatty acids (PUFA). After quality control of the genotypic data, 41,057 SNP markers remained for genomic analyses. Meteorological data from the weather station in closest herd distance were used for the calculation of maximum hourly daily THI. Genomic reaction norm models were applied to estimate genetic parameters in a single-step approach for production traits and FA in dependency of THI at different lactation stages, and to evaluate the model stability. In a first evaluation strategy (New_sire), all phenotypic records from daughters of genotyped sires born after 2010 were masked, to mimic a validation population. In the second strategy (New_env), only daughter records of the new sires recorded in the most extreme THI classes were masked, aiming at predicting sire genomic estimated breeding values (GEBV) under heat stress conditions. Model stability was the correlation between GEBV of the new sires in the reduced data set with respective GEBV estimated from all phenotypic data. Among all test-day production traits, PY responded as the most sensitive to heat stress. As observed for the remaining production traits, genetic variances were quite stable across THI, but genetic correlations between PY from temperate climates with PY from extreme THI classes dropped to 0.68. Genetic variances in dependency of THI were very similar for C16:0 and SFA, indicating marginal climatic sensitivity. In the early lactation stage, genetic variances for C18:0, MUFA, PUFA, and UFA were significantly larger in the extreme THI classes compared with the estimates under thermoneutral conditions. For C18:0 and MUFA, PUFA, and UFA in the middle THI classes, genetic correlations in same traits from the early and the later lactation stages were lower than 0.50, indicating strong days in milk influence. Interestingly, within lactation stages, genetic correlations for C18:0 and UFA recorded at low and high THI were quite large, indicating similar genetic mechanisms under stress conditions. The model stability was improved when applying the New_env instead of New_sire strategy, especially for FA in the first stage of lactation. Results indicate moderately accurate genomic predictions for milk traits in extreme THI classes when considering phenotypic data from a broad range of remaining THI. Phenotypically, thermal stress conditions contributed to an increase of UFA, suggesting value as a heat stress biomarker. Furthermore, the quite large genetic variances for UFA at high THI suggest the consideration of UFA in selection strategies for improved heat stress resistance.

Breed Ancestry, Divergence, Admixture, and Selection Patterns of the Simbra Crossbreed

In this study, we evaluated an admixed South African Simbra crossbred population, as well as the Brahman (Indicine) and Simmental (Taurine) ancestor populations to understand their genetic architecture and detect genomic regions showing signatures of selection. Animals were genotyped using the Illumina BovineLD v2 BeadChip (7K). Genomic structure analysis confirmed that the South African Simbra cattle have an admixed genome, composed of 5/8 Taurine and 3/8 Indicine, ensuring that the Simbra genome maintains favorable traits from both breeds. Genomic regions that have been targeted by selection were detected using the linkage disequilibrium-based methods iHS and Rsb. These analyses identified 10 candidate regions that are potentially under strong positive selection, containing genes implicated in cattle health and production (e.g., TRIM63, KCNA10, NCAM1, SMIM5, MIER3, and SLC24A4). These adaptive alleles likely contribute to the biological and cellular functions determining phenotype in the Simbra hybrid cattle breed. Our data suggested that these alleles were introgressed from the breed's original indicine and taurine ancestors. The Simbra breed thus possesses derived parental alleles that combine the superior traits of the founder Brahman and Simmental breeds. These regions and genes might represent good targets for ad-hoc physiological studies, selection of breeding material and eventually even gene editing, for improved traits in modern cattle breeds. This study represents an important step toward developing and improving strategies for selection and population breeding to ultimately contribute meaningfully to the beef production industry.

Genotype by environment interaction for Holstein cattle populations using autoregressive and within- and across-country multi-trait reaction norms test-day models

The progenies of international bulls in diverse climatic conditions and management levels may lead to different expressions of their genetic potential resulting in a re-ranking of these bulls. Therefore, evaluate the presence of genotype by environment interaction (G×E) within and across countries is important to guide the decision-making on alternative selection strategies. Thus, a two-step reaction norm (RN) approach was used to investigate the presence of G×E in Portuguese and Brazilian Holstein cattle. In step 1, we performed a within-country genetic evaluation using an autoregressive model to obtain precorrected phenotypes and environmental gradients (herd test-day solutions, HTD levels). In step 2, the precorrected phenotypes were considered as two distinct traits in a bi-trait RN model to estimate variance components across HTD levels, genetic correlation between HTD levels in Portugal and Brazil, and RN of the estimated breeding values. Additionally, the genetic correlation between countries using a bi-trait random regression (RR) sire model was obtained. In step 1, genetic additive variance for milk yield (MY) in Portugal was 14.1% higher than in Brazil. For somatic cell score (SCS), the genetic additive variance in Portugal was 12.7% lower than in Brazil. Although similar heritability estimates for SCS were observed in both countries, MY heritabilities were 0.31 for Portugal and 0.23 for Brazil. Genetic correlations (SD) between both countries obtained using RR sire model were 0.78 (0.051) for MY and 0.75 (0.062) for SCS. In step 2, MY genetic correlations among HTD levels within countries were higher than 0.94 for Portugal and 0.98 for Brazil. Somatic cell score genetic correlations among HTD levels ranged from 0.70 to 0.99 for Portugal and from 0.84 to 0.99 for Brazil. The average (SD) of genetic correlation estimates between Portuguese and Brazilian HTD levels were 0.74 (0.009) for MY and 0.57 (0.060) for SCS. These results suggest the presence of G×E for MY and SCS of Holstein cattle between both countries. Although there was no indication of G×E between Brazilian herd environments, the low genetic correlation for SCS indicates potential re-ranking of bulls between extreme environmental gradient in Portugal. Overall, the results of this study evidence the importance of national and international genetic evaluation systems to assist dairy farmers in the selection of the best genotypes to obtain the expected returns from investments in imported semen and to realize genetic progress in dairy populations under local environmental conditions.

Genotype-by-environment interactions at the trait level and total merit index level for milk production and functional traits in Brown Swiss cattle

The production environments of the German-Austrian Brown Swiss population show a wide range due to differences in topography, landscapes, local climates, and different farm management systems. Extensive production systems such as organic farming have become increasingly popular in recent decades because of interest in sustainability and consumer preferences. Compared with conventional farmers, organic farmers put more weight on fitness traits. Besides the official total merit index (TMI), a selection index applying relative economic weights (REWs) suitable for organic production systems is provided for Brown Swiss cattle in Germany. The aim of the study was to investigate genotype-by-environment interactions (GxE) for milk production traits and functional traits (including longevity, fertility traits, and calving traits) in a sample of the German-Austrian Brown Swiss population housed in Baden-Wuerttemberg (southern Germany) by applying bivariate and random regression sire models. For bivariate analyses, the production environment was binary classified by farm management system (organic and conventional) and altitude of farm location (above or below 800 m above sea level (ASL)). Milk energy yields (MEY) obtained from herd effects were used as continuously scaled environmental descriptor in the reaction norm approach. The TMIs for sires were calculated based on breeding values estimated with different models and environment-specific REWs to determine possible GxE at TMI levels and rerankings of sires. In bivariate analyses, genetic correlations at the trait level were high and ranged from r_g = 0.99 (calving to first insemination, cystic ovaries, and maternal stillbirth rate) to r_g = 0.79 (first insemination to conception for altitude). Except for the latter, no severe GxE were found at the trait level using the bivariate models. Fat yield was the only trait showing minor GxE in the reaction norm model approach. Investigating the environmental sensitivity at the TMI level revealed rank correlations between the different environment-specific TMIs that were close to unity, implying no severe reranking effects. The results show no need to account for different environments in Brown Swiss cattle breeding programs.

The superiority of multi-trait models with genotype-by-environment interactions in a limited number of environments for genomic prediction in pigs

Background

Different production systems and climates could lead to genotype-by-environment (G × E) interactions between populations, and the inclusion of G × E interactions is becoming essential in breeding decisions. The objective of this study was to investigate the performance of multi-trait models in genomic prediction in a limited number of environments with G × E interactions.

Results

In total, 2,688 and 1,384 individuals with growth and reproduction phenotypes, respectively, from two Yorkshire pig populations with similar genetic backgrounds were genotyped with the PorcineSNP80 panel. Single- and multi-trait models with genomic best linear unbiased prediction (GBLUP) and BayesC π were implemented to investigate their genomic prediction abilities with 20 replicates of five-fold cross-validation. Our results regarding between-environment genetic correlations of growth and reproductive traits (ranging from 0.618 to 0.723) indicated the existence of G × E interactions between these two Yorkshire pig populations. For single-trait models, genomic prediction with GBLUP was only 1.1% more accurate on average in the combined population than in single populations, and no significant improvements were obtained by BayesC π for most traits. In addition, single-trait models with either GBLUP or BayesC π produced greater bias for the combined population than for single populations. However, multi-trait models with GBLUP and BayesC π better accommodated G × E interactions, yielding 2.2% – 3.8% and 1.0% – 2.5% higher prediction accuracies for growth and reproductive traits, respectively, compared to those for single-trait models of single populations and the combined population. The multi-trait models also yielded lower bias and larger gains in the case of a small reference population. The smaller improvement in prediction accuracy and larger bias obtained by the single-trait models in the combined population was mainly due to the low consistency of linkage disequilibrium between the two populations, which also caused the BayesC π method to always produce the largest standard error in marker effect estimation for the combined population.

Conclusions

In conclusion, our findings confirmed that directly combining populations to enlarge the reference population is not efficient in improving the accuracy of genomic prediction in the presence of G × E interactions, while multi-trait models perform better in a limited number of environments with G × E interactions.

Genomic prediction of growth traits for pigs in the presence of genotype by environment interactions using single-step genomic reaction norm model

Economically important traits are usually complex traits influenced by genes, environment and genotype-by-environment (G × E) interactions. Ignoring G × E interaction could lead to bias in the estimation of breeding values and selection decisions. A total of 1,778 pigs were genotyped using the PorcineSNP80 BeadChip. The existence of G × E interactions was investigated using a single-step reaction norm model for growth traits of days to 100 kg (AGE) and backfat thickness adjusted to 100 kg (BFT), based on a pedigree-based relationship matrix (A) or a genomic-pedigree joint relationship matrix (H). In the reaction norm model, the herd-year-season effect was measured as the environmental variable (EV). Our results showed no G × E interactions for AGE, but for BFT. For both AGE and BFT, the genomic reaction norm model (H) produced more accurate predictions than the conventional reaction norm model (A). For BFT, the accuracies were greater based on the reaction norm model than those based on the reduced model without exploiting G × E interaction, with EV ranging from 0.5 to 1, and accuracy increasing by 3.9% and 4.6% in the reaction norm model based on A and H matrices, respectively, while reaction norm model yielded approximately 8.4% and 7.9% lower accuracy for EVs ranging from 0 to 0.4, based on A and H matrices, respectively. In addition, for BFT, the highest accuracy was obtained in the BJLM6 farm for realizing directional selection. This study will help to apply G × E interactions to practical genomic selection.

Genetic Parameters and Genome-Wide Association Studies for Anti-Müllerian Hormone Levels and Antral Follicle Populations Measured After Estrus Synchronization in Nellore Cattle

Reproductive efficiency plays a major role in the long-term sustainability of livestock industries and can be improved through genetic and genomic selection. This study aimed to estimate genetic parameters (heritability and genetic correlation) and identify genomic regions and candidate genes associated with anti-Müllerian hormone levels (AMH) and antral follicle populations measured after estrous synchronization (AFP) in Nellore cattle. The datasets included phenotypic records for 1099 and 289 Nellore females for AFP and AMH, respectively, high-density single nucleotide polymorphism (SNP) genotypes for 944 animals, and 4129 individuals in the pedigree. The heritability estimates for AMH and AFP were 0.28 ± 0.07 and 0.30 ± 0.09, and the traits were highly and positively genetically correlated (rG = 0.81 ± 0.02). These findings indicated that these traits can be improved through selective breeding, and substantial indirect genetic gains are expected by selecting for only one of the two traits. A total of 31 genomic regions were shown to be associated with AMH or AFP, and two genomic regions located on BTA1 (64.9-65.0 Mb and 109.1-109.2 Mb) overlapped between the traits. Various candidate genes were identified to be potentially linked to important biological processes such as ovulation, tissue remodeling, and the immune system. Our findings support the use of AMH and AFP as indicator traits to genetically improve fertility rates in Nellore cattle and identify better oocyte donors.

Genomics of Heat Tolerance in Reproductive Performance Investigated in Four Independent Maternal Lines of Pigs

Improving swine climatic resilience through genomic selection has the potential to minimize welfare issues and increase the industry profitability. The main objective of this study was to investigate the genetic and genomic determinism of tolerance to heat stress in four independent purebred populations of swine. Three female reproductive traits were investigated: total number of piglets born (TNB), number of piglets born alive (NBA) and average birth weight (ABW). More than 80,000 phenotypic and 12,000 genotyped individuals were included in this study. Genomic random-regression models were fitted regressing the phenotypes of interest on a set of 95 environmental covariates extracted from public weather station records. The models yielded estimates of (genomic) reactions norms for individual pigs, as indicator of heat tolerance. Heat tolerance is a heritable trait, although the heritabilities are larger under comfortable than heat-stress conditions (larger than 0.05 vs. 0.02 for TNB; 0.10 vs. 0.05 for NBA; larger than 0.20 vs. 0.10 for ABW). TNB showed the lowest genetic correlation (-38%) between divergent climatic conditions, being the trait with the strongest impact of genotype by environment interaction, while NBA and ABW showed values slightly negative or equal to zero reporting a milder impact of the genotype by environment interaction. After estimating genetic parameters, a genome-wide association study was performed based on the single-step GBLUP method. Heat tolerance was observed to be a highly polygenic trait. Multiple and non-overlapping genomic regions were identified for each trait based on the genomic breeding values for reproductive performance under comfortable or heat stress conditions. Relevant regions were found on chromosomes (SSC) 1, 3, 5, 6, 9, 11, and 12, although there were important regions across all autosomal chromosomes. The genomic region located on SSC9 appears to be of particular interest since it was identified for two traits (TNB and NBA) and in two independent populations. Heat tolerance based on reproductive performance indicators is a heritable trait and genetic progress for heat tolerance can be achieved through genetic or genomic selection. Various genomic regions and candidate genes with important biological functions were identified, which will be of great value for future functional genomic studies.