Evaluation of sire by environment interactions for growth rate and backfat depth using reaction norm models in pigs

Validation of reaction norm breeding values for robustness in Australian sheep

BACKGROUND

There can be variation between animals in how stable their genetic merit is across different environments due to genotype-by-environment (G×E) interactions. This variation could be used in breeding programs to select robust genotypes that combine high overall performance with stable genetic ranking across environments. There have been few attempts to validate breeding values for robustness in livestock, although this is a necessary step towards their implementation in selection decisions. The objective of this study was to validate breeding values for the robustness of body weight across different growth environments that were estimated using reaction norm models in sheep data.

RESULTS

Using threefold cross-validation for the progeny of 337 sires, the average correlation between single-step breeding values for the reaction norm slope and the realised robustness of progeny across different growth environments was 0.21. The correlation between breeding values for the reaction slope estimated independently in two different datasets linked by common sires was close to the expected correlation based on theory.

CONCLUSIONS

Slope estimated breeding values (EBV) obtained using reaction norm models were predictive of the phenotypic robustness of progeny across different environments and were consistent for sires with progeny in two different datasets. Selection based on reaction norm EBV could be used to increase the robustness of a population to environmental variation.

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.

Estimation of macro- and micro-genetic environmental sensitivity in unbalanced datasets

Genotype-by-environment interaction is caused by variation in genetic environmental sensitivity (GES), which can be subdivided into macro- and micro-GES. Macro-GES is genetic sensitivity to macro-environments (definable environments often shared by groups of animals), while micro-GES is genetic sensitivity to micro-environments (individual environments). A combined reaction norm and double hierarchical generalised linear model (RN-DHGLM) allows for simultaneous estimation of base genetic, macro- and micro-GES effects. The accuracy of variance components estimated using a RN-DHGLM has been explicitly studied for balanced data and recommendation of a data size with a minimum of 100 sires with at least 100 offspring each have been made. In the current study, the data size (numbers of sires and progeny) and structure requirements of the RN-DHGLM were investigated for two types of unbalanced datasets. Both datasets had a variable number of offspring per sire, but one dataset also had a variable number of offspring within macro-environments. The accuracy and bias of the estimated macro- and micro-GES effects and the estimated breeding values (EBVs) obtained using the RN-DHGLM depended on the data size. Reasonably accurate and unbiased estimates were obtained with data containing 500 sires with 20 offspring or 100 sires with 50 offspring, regardless of the data structure. Variable progeny group sizes, alone or in combination with an unequal number of offspring within macro-environments, had little impact on the dispersion of the EBVs or the bias and accuracy of variance component estimation, but resulted in lower accuracies of the EBVs. Compared to genetic correlations of zero, a genetic correlation of 0.5 between base genetic, macro- and micro-GES components resulted in a slight decrease in the percentage of replicates that converged out of 100 replicates, but had no effect on the dispersion and accuracy of variance component estimation or the dispersion of the EBVs. The results show that it is possible to apply the RN-DHGLM to unbalanced datasets to obtain estimates of variance due to macro- and micro-GES. Furthermore, the levels of accuracy and bias of variance estimates when analysing macro- and micro-GES simultaneously are determined by average family size, with limited impact from variability in family size and/or cohort size. This creates opportunities for the use of field data from populations with unbalanced data structures when estimating macro- and micro-GES.

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.

Reaction norm analysis of pig growth using environmental descriptors based on alternative traits

Contemporary group (CG) estimates of different phenotypes have not been widely explored for pigs. The objective of this study was to extend the traits used to derive environmental descriptors of the growing pig, to include CG estimates of early growth between birth and start of feed intake test (EADG), growth during feed intake test (TADG), lifetime growth (ADG), daily feed intake (DFI), backfat (BF) and muscle depth (MD). Pedigree and performance records (n = 7,746) from a commercial Australian piggery were used to derive environmental descriptors based on CG estimates of these six traits. The CG estimates of growth traits described different aspects of the environment from the CG estimates of carcass traits (r < 0.10). These definitions of the environment then were used in reaction norm analysis of growth, to evaluate sire-by-environment interaction (Sire × E) for growth. The most appropriate reaction norm model to evaluate Sire × E for growth was dependent on the environmental descriptor used. If the trait used to derive an environmental descriptor was distinctly different from growth (e.g., BF and MD), CG as an additional random effect was required in the model. If not included, inflated common litter effect and sire intercept variance suggest there was unaccounted environmental variability. There was no significant Sire × E using any of the definitions of the environment, with estimated variance in sire slopes largest when environments were defined by BF ( σ ^ bi 2  = 97 ± 83 (g/day)² ), followed by environments defined by DFI ( σ ^ bi 2  = 39 ± 101 (g/day)² ). While there appears to be differences in ability to detect Sire × E, improved data structure is required to better assess these environmental descriptors based on alternative traits. The ideal trait, or combination of traits, used to derive environmental descriptors may be unique for individual herds. Therefore, multiple phenotypes should be further explored for the evaluation of Sire × E for growth in the pig.

Macro-environmental sensitivity for growth rate in Danish Duroc pigs is under genetic control

The aim of this study was to examine (i) the genetic variation in macro-environmental sensitivity (macro-ES) for ADG in Danish Duroc pigs, (ii) the genetic heterogeneity among sexes, and (iii) residual variance heterogeneity among herds. Record of ADG for 32,297 boars (19 herds) and 42,724 gilts (16 herds) was used for analysis. The data were provided by the National Danish Pig Research Centre. The analysis was performed by fitting univariate reaction norm models with the herd-year-month on test (HYM) effect as environmental covariates and herd-specific residual variance for boars and gilts separately under a Bayesian setting. The environmental covariate was inferred simultaneously with other parameters of the model. Gibbs sampling was used to sample model dispersion and location parameters. The posterior means and highest posterior density intervals of the additive genetic variance, genetic correlations for ADG, and heritability were calculated over the continuous environmental range of -3σh to +3σh (SD of the HYM effect). The coheritability of ADG at the average environmental level and ADG in the environments along the -3σh to +3σh environmental gradient were also calculated. The analysis showed significant variation in macro-ES, revealing genotype by environment interactions (G × E) for ADG. The presence of G × E resulted in changes in additive genetic variance and heritability across the -3σh to +3σh range. The genetic correlations were high and positive between ADG in environments differing by 1σh units or less and decreased to moderately positive between ADG in the extreme environments in both sexes. The coheritability of ADG in the environment at the average level and the -3σh environment for boars were greater than the heritability in the environment at the average level, while it was less for gilts. The coheritability of ADG in the environment at the average level and the +3σh environment for boars was less than heritability in the environment at the average level, while it was either the same or greater for gilts, depending on the residual variance. Boars had larger additive genetic and residual variances than gilts. Heterogeneity of residual variances across herds was shown for both sexes. In conclusion, this study shows the presence of macro-ES, genetic variance heterogeneity among sexes for ADG in pigs, and residual variance heterogeneity across herds.

Contemporary group estimates adjusted for climatic effects provide a finer definition of the unknown environmental challenges experienced by growing pigs

Environmental descriptors derived from mean performances of contemporary groups (CGs) are assumed to capture any known and unknown environmental challenges. The objective of this paper was to obtain a finer definition of the unknown challenges, by adjusting CG estimates for the known climatic effects of monthly maximum air temperature (MaxT), minimum air temperature (MinT) and monthly rainfall (Rain). As the unknown component could include infection challenges, these refined descriptors may help to better model varying responses of sire progeny to environmental infection challenges for the definition of disease resilience. Data were recorded from 1999 to 2013 at a piggery in south-east Queensland, Australia (n = 31,230). Firstly, CG estimates of average daily gain (ADG) and backfat (BF) were adjusted for MaxT, MinT and Rain, which were fitted as splines. In the models used to derive CG estimates for ADG, MaxT and MinT were significant variables. The models that contained these significant climatic variables had CG estimates with a lower variance compared to models without significant climatic variables. Variance component estimates were similar across all models, suggesting that these significant climatic variables accounted for some known environmental variation captured in CG estimates. No climatic variables were significant in the models used to derive the CG estimates for BF. These CG estimates were used to categorize environments. There was no observable sire by environment interaction (Sire×E) for ADG when using the environmental descriptors based on CG estimates on BF. For the environmental descriptors based on CG estimates of ADG, there was significant Sire×E only when MinT was included in the model (p = .01). Therefore, this new definition of the environment, preadjusted by MinT, increased the ability to detect Sire×E. While the unknown challenges captured in refined CG estimates need verification for infection challenges, this may provide a practical approach for the genetic improvement of disease resilience.

A sandwich ELISA for porcine alpha-1 acid glycoprotein (pAGP, ORM-1) and further demonstration of its use to evaluate growth potential in newborn pigs

A simple, reproducible sandwich, ELISA was developed to measure porcine alpha-1 acid glycoprotein (pAGP, ORM-1) in pig plasma. Porcine AGP isolated from serum was purchased and a polyclonal antisera was prepared in rabbits using the whole pAGP molecule as immunogen. The antiserum was affinity purified, and a portion of the purified antibody fraction was labeled with horseradish peroxidase. Porcine AGP protein was used as a standard, whereas commercially available buffers and reagents were utilized throughout the assay. The assay was specific for pAGP, had a lower limit of detection of 3.2 ng/mL, and could be used to quantify pAGP in plasma or serum. Using this ELISA, we corroborated our previous findings obtained by RID assay, which demonstrated that the AGP concentration in newborn piglets is negatively associated with preweaning growth rate. The current data were obtained using piglets from a different geographical location and genetic background and showed that elevated AGP at birth was associated with reduced preweaning growth rate (P < 0.001, r = 0.433, n = 19 litters). In addition, litters with a greater average AGP at birth were at a growth disadvantage compared with litters with reduced average AGP plasma concentrations (P < 0.001, r = 0.708, n = 19 litters). Litter average plasma AGP was a better predictor of litter preweaning growth rate than average litter birth weight. The data represent further support for using perinatal AGP concentrations as a tool to identify potential slower growing pigs and as a plasma biomarker for predicting litter growth rate.

What does the ‘closed herd’ really mean for Australian breeding companies and their customers?

The perception that the genetic background of the Australian pig population is limiting for genetic improvement of commercial pigs in Australia is considered in the context of well established theory combined with practical evidence. The diversity of pig breeds used in modern commercial pig-breeding programs is diminished worldwide relative to all the pig breeds available. Australia is no different in this respect. The use of predominantly three main breeds (Large White, Landrace, Duroc) and synthetic lines, with contributions from other minor breeds to form the basis of a cross-breeding system for commercial pig production is well established internationally. The Australian concern of relatively small founder populations is potentially of relevance, from a theoretical perspective, for (1) the prevalence of defects or the presence of desirable alleles, and (2) the loss of genetic variation or increase in inbreeding depression resulting from increased inbreeding in closed nucleus lines, potentially reducing response to selection. However, rates of response achieved in Australian herds are generally commensurate with the performance recording and selection emphasis applied, and do not appear to be unduly restricted. Moreover, favourable alleles present in unrepresented breeds are frequently present in the three major breeds elsewhere, and therefore would be expected to be present within the Australian populations. Wider testing would provide confirmation of this. Comparison of estimates of effective population size of Australian populations with experimental selection lines overseas (e.g. INRA) or other intensely selected species (e.g. Holstein cattle) suggest adequate genetic diversity to achieve ongoing genetic improvement in the Australian pig industry. However, fitness traits should be included in breeding goals. What remains to be seen is whether novel phenotypes or genotypes are required to meet future challenges, which might be imposed by changes in the environment (e.g. climate change, disease) or market needs. Given probable overlap in genetic merit across Australian and foreign populations for unselected attributes, we suggest that sufficient genetic resources are already present in Australian herds to continue commercial progress within existing Australian populations that have adapted to Australian conditions.