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

Exploring the value of genomic predictions to simultaneously improve production potential and resilience of farmed animals

Sustainable livestock production requires that animals have a high production potential but are also highly resilient to environmental challenges. The first step to simultaneously improve these traits through genetic selection is to accurately predict their genetic merit. In this paper, we used simulations of sheep populations to assess the effect of genomic data, different genetic evaluation models and phenotyping strategies on prediction accuracies and bias for production potential and resilience. In addition, we also assessed the effect of different selection strategies on the improvement of these traits. Results show that estimation of both traits greatly benefits from taking repeated measurements and from using genomic information. However, the prediction accuracy for production potential is compromised, and resilience estimates tends to be upwards biased, when families are clustered in groups even when genomic information is used. The prediction accuracy was also found to be lower for both traits, resilience and production potential, when the environment challenge levels are unknown. Nevertheless, we observe that genetic gain in both traits can be achieved even in the case of unknown environmental challenge, when families are distributed across a large range of environments. Simultaneous genetic improvement in both traits however greatly benefits from the use of genomic evaluation, reaction norm models and phenotyping in a wide range of environments. Using models without the reaction norm in scenarios where there is a trade-off between resilience and production potential, and phenotypes are collected from a narrow range of environments may result in a loss for one trait. The study demonstrates that genomic selection coupled with reaction-norm models offers great opportunities to simultaneously improve productivity and resilience of farmed animals even in the case of a trade-off.

Improving Genomic Prediction Accuracy in the Chinese Holstein Population by Combining with the Nordic Holstein Reference Population

The size of the reference population is critical in order to improve the accuracy of genomic prediction. Indeed, improving genomic prediction accuracy by combining multinational reference populations has proven to be effective. In this study, we investigated the improvement of genomic prediction accuracy in seven complex traits (i.e., milk yield; fat yield; protein yield; somatic cell count; body conformation; feet and legs; and mammary system conformation) by combining the Chinese and Nordic Holstein reference populations. The estimated genetic correlations between the Chinese and Nordic Holstein populations are high with respect to protein yield, fat yield, and milk yield-whereby these correlations range from 0.621 to 0.720-and are moderate with respect to somatic cell count (0.449), but low for the three conformation traits (which range from 0.144 to 0.236). When utilizing the joint reference data and a two-trait GBLUP model, the genomic prediction accuracy in the Chinese Holsteins improves considerably with respect to the traits with moderate-to-high genetic correlations, whereas the improvement in Nordic Holsteins is small. When compared with the single population analysis, using the joint reference population for genomic prediction in younger animals, results in a 2.3 to 8.1 percent improvement in accuracy. Meanwhile, 10 replications of five-fold cross-validation were also implemented in order to evaluate the performance of joint genomic prediction, thereby resulting in a 1.6 to 5.2 percent increase in accuracy. With respect to joint genomic prediction, the bias was found to be quite low. However, for traits with low genetic correlations, the joint reference data do not improve the prediction accuracy substantially for either population.

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.

Genetic analysis of disease resilience of wean-to-finish pigs under a natural disease challenge model using reaction norms

Background

Disease resilience is the ability to maintain performance across environments with different disease challenge loads (CL). A reaction norm describes the phenotypes that a genotype can produce across a range of environments and can be implemented using random regression models. The objectives of this study were to: (1) develop measures of CL using growth rate and clinical disease data recorded under a natural polymicrobial disease challenge model; and (2) quantify genetic variation in disease resilience using reaction norm models.

Methods

Different CL were derived from contemporary group effect estimates for average daily gain (ADG) and clinical disease phenotypes, including medical treatment rate (TRT), mortality rate, and subjective health scores. Resulting CL were then used as environmental covariates in reaction norm analyses of ADG and TRT in the challenge nursery and finisher, and compared using model loglikelihoods and estimates of genetic variance associated with CL. Linear and cubic spline reaction norm models were compared based on goodness-of-fit and with multi-variate analyses, for which phenotypes were separated into three traits based on low, medium, or high CL.

Results

Based on model likelihoods and estimates of genetic variance explained by the reaction norm, the best CL for ADG in the nursery was based on early ADG in the finisher, while the CL derived from clinical disease traits across the nursery and finisher was best for ADG in the finisher and for TRT in the nursery and across the nursery and finisher. With increasing CL, estimates of heritability for nursery and finisher ADG initially decreased, then increased, while estimates for TRT generally increased with CL. Genetic correlations for ADG and TRT were low between high versus low CL, but high for close CL. Linear reaction norm models fitted the data significantly better than the standard genetic model without genetic slopes, while the cubic spline model fitted the data significantly better than the linear reaction norm model for most traits. Reaction norm models also fitted the data better than multi-variate models.

Conclusions

Reaction norm models identified genotype-by-environment interactions related to disease CL. Results can be used to select more resilient animals across different levels of CL, high-performance animals at a given CL, or a combination of these.

The theory on and software simulating large-scale genomic data for genotype-by-environment interactions

Background

With the emphasis on analysing genotype-by-environment interactions within the framework of genomic selection and genome-wide association analysis, there is an increasing demand for reliable tools that can be used to simulate large-scale genomic data in order to assess related approaches.

Results

We proposed a theory to simulate large-scale genomic data on genotype-by-environment interactions and added this new function to our developed tool GPOPSIM. Additionally, a simulated threshold trait with large-scale genomic data was also added. The validation of the simulated data indicated that GPOSPIM2.0 is an efficient tool for mimicking the phenotypic data of quantitative traits, threshold traits, and genetically correlated traits with large-scale genomic data while taking genotype-by-environment interactions into account.

Conclusions

This tool is useful for assessing genotype-by-environment interactions and threshold traits methods.

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.

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.