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

Effect of SNPs on Litter Size in Swine

Although sows do not directly enter the market, they play an important role in piglet breeding on farms. They consume large amounts of feed, resulting in a significant environmental burden. Pig farms can increase their income and reduce environmental pollution by increasing the litter size (LS) of swine. PCR-RFLP/SSCP and GWAS are common methods to evaluate single-nucleotide polymorphisms (SNPs) in candidate genes. We conducted a systematic meta-analysis of the effect of SNPs on pig LS. We collected and analysed data published over the past 30 years using traditional and network meta-analyses. Trial sequential analysis (TSA) was used to analyse population data. Gene set enrichment analysis and protein-protein interaction network analysis were used to analyse the GWAS dataset. The results showed that the candidate genes were positively correlated with LS, and defects in PCR-RFLP/SSCP affected the reliability of candidate gene results. However, the genotypes with high and low LSs did not have a significant advantage. Current breeding and management practices for sows should consider increasing the LS while reducing lactation length and minimizing the sows' non-pregnancy period as much as possible.

Genome-wide association study considering genotype-by-environment interaction for productive and reproductive traits using whole-genome sequencing in Nellore cattle

BACKGROUND

The genotype-by-environment interaction (GxE) in beef cattle can be investigated using reaction norm models to assess environmental sensitivity and, combined with genome-wide association studies (GWAS), to map genomic regions related to animal adaptation. Including genetic markers from whole-genome sequencing in reaction norm (RN) models allows us to identify high-resolution candidate genes across environmental gradients through GWAS. Hence, we performed a GWAS via the RN approach using whole-genome sequencing data, focusing on mapping candidate genes associated with the expression of reproductive and growth traits in Nellore cattle. For this purpose, we used phenotypic data for age at first calving (AFC), scrotal circumference (SC), post-weaning weight gain (PWG), and yearling weight (YW). A total of 20,000 males and 7,159 females genotyped with 770k were imputed to the whole sequence (29 M). After quality control and linkage disequilibrium (LD) pruning, there remained ∼ 2.41 M SNPs for SC, PWG, and YW and ∼ 5.06 M SNPs for AFC.

RESULTS

Significant SNPs were identified on Bos taurus autosomes (BTA) 10, 11, 14, 18, 19, 20, 21, 24, 25 and 27 for AFC and on BTA 4, 5 and 8 for SC. For growth traits, significant SNP markers were identified on BTA 3, 5 and 20 for YW and PWG. A total of 56 positional candidate genes were identified for AFC, 9 for SC, 3 for PWG, and 24 for YW. The significant SNPs detected for the reaction norm coefficients in Nellore cattle were found to be associated with growth, adaptative, and reproductive traits. These candidate genes are involved in biological mechanisms related to lipid metabolism, immune response, mitogen-activated protein kinase (MAPK) signaling pathway, and energy and phosphate metabolism.

CONCLUSIONS

GWAS results highlighted differences in the physiological processes linked to lipid metabolism, immune response, MAPK signaling pathway, and energy and phosphate metabolism, providing insights into how different environmental conditions interact with specific genes affecting animal adaptation, productivity, and reproductive performance. The shared genomic regions between the intercept and slope are directly implicated in the regulation of growth and reproductive traits in Nellore cattle raised under different environmental conditions.

Using lifestyle information in polygenic modeling of blood pressure traits: a simple method to reduce bias

Complex traits are determined by the effects of multiple genetic variants, multiple environmental factors, and potentially their interaction. Predicting complex trait phenotypes from genotypes is a fundamental task in quantitative genetics that was pioneered in agricultural breeding for selection purposes. However, it has recently become important in human genetics. While prediction accuracy for some human complex traits is appreciable, this remains low for most traits. A promising way to improve prediction accuracy is by including not only genetic information but also environmental information in prediction models. However, environmental factors can, in turn, be genetically determined. This phenomenon gives rise to a correlation between the genetic and environmental components of the phenotype, which violates the assumption of independence between the genetic and environmental components of most statistical methods for polygenic modeling. In this work, we investigated the impact of including 27 lifestyle variables as well as genotype information (and their interaction) for predicting diastolic blood pressure, systolic blood pressure, and pulse pressure in older individuals in UK Biobank. The 27 lifestyle variables were included as either raw variables or adjusted by genetic and other non-genetic factors. The results show that including both lifestyle and genetic data improved prediction accuracy compared to using either piece of information alone. Both prediction accuracy and bias can improve substantially for some traits when the models account for the lifestyle variables after their proper adjustment. Our work confirms the utility of including environmental information in polygenic models of complex traits and highlights the importance of proper handling of the environmental variables.

CYP24A1 is associated with fetal mummification in pigs

Mummified piglets are among the leading causes of fertility loss and severely hamper reproductive performance in pigs. However, the contributions of genomic variation to the emergence of mummified piglets (MUM) have rarely been studied. This study aims to (1) elucidate the genetic architecture of MUM in sows of parity 1 - 3 using a single-step genome-wide association study (ssGWAS). The ssGWAS involved genotyping-by-sequencing of Large White and Landrace pig breeds. (2) Explore the biological role of the candidate genes at the cellular level. A total of 185 and 48 genome-wide significant SNPs are associated with MUM in Large White and Landrace pigs, explaining 0.01-36.52% genetic variance for different significant loci, respectively. All the significant SNPs are parity-specific, and the numerous, consecutive significant loci likely generated the nine significant peaks in different parities. Multiple candidate genes (including CYP24A1, FBXO30, and ARHGEF28) are associated with fetal congenital and maternal diseases. Collectively, CYP24A1 regulation contributes to steady-state levels of embryo development genes. CYP24A1 is involved in reproduction and, immune and gestational disorders. Thus, it is associated with known newborn death traits and MUM in Large White sows. Altogether, these results improve the current understanding of the genetic architecture of MUM and expand the knowledge on genetic variations for selecting against mummified piglets in pig breeding.

Genetic parameters for automatically-measured vaginal temperature, respiration efficiency, and other thermotolerance indicators measured on lactating sows under heat stress conditions

BACKGROUND

Genetic selection based on direct indicators of heat stress could capture additional mechanisms that are involved in heat stress response and enable more accurate selection for more heat-tolerant individuals. Therefore, the main objectives of this study were to estimate genetic parameters for various heat stress indicators in a commercial population of Landrace × Large White lactating sows measured under heat stress conditions. The main indicators evaluated were: skin surface temperatures (SST), automatically-recorded vaginal temperature (TV), respiration rate (RR), panting score (PS), body condition score (BCS), hair density (HD), body size (BS), ear size, and respiration efficiency (Reff).

RESULTS

Traits based on TV presented moderate heritability estimates, ranging from 0.15 ± 0.02 to 0.29 ± 0.05. Low heritability estimates were found for SST traits (from 0.04 ± 0.01 to 0.06 ± 0.01), RR (0.06 ± 0.01), PS (0.05 0.01), and Reff (0.03 ± 0.01). Moderate to high heritability values were estimated for BCS (0.29 ± 0.04 for caliper measurements and 0.25 ± 0.04 for visual assessments), HD (0.25 ± 0.05), BS (0.33 ± 0.05), ear area (EA; 0.40 ± 0.09), and ear length (EL; 0.32 ± 0.07). High genetic correlations were estimated among SST traits (> 0.78) and among TV traits (> 0.75). Similarly, high genetic correlations were also estimated for RR with PS (0.87 ± 0.02), with BCS measures (0.92 ± 0.04), and with ear measures (0.95 ± 0.03). Low to moderate positive genetic correlations were estimated between SST and TV (from 0.25 ± 0.04 to 0.76 ± 0.07). Low genetic correlations were estimated between TV and BCS (from - 0.01 ± 0.08 to 0.06 ± 0.07). Respiration efficiency was estimated to be positively and moderately correlated with RR (0.36 ± 0.04), PS (0.56 ± 0.03), and BCS (0.56 ± 0.05 for caliper measurements and 0.50 ± 0.05 for the visual assessments). All other trait combinations were lowly genetically correlated.

CONCLUSIONS

A comprehensive landscape of heritabilities and genetic correlations for various thermotolerance indicators in lactating sows were estimated. All traits evaluated are under genetic control and heritable, with different magnitudes, indicating that genetic progress is possible for all of them. The genetic correlation estimates provide evidence for the complex relationships between these traits and confirm the importance of a sub-index of thermotolerance traits to improve heat tolerance in pigs.

Genomic prediction and association mapping of maize grain yield in multi-environment trials based on reaction norm models

Genotype-by-environment interaction (GEI) is among the greatest challenges for maize breeding programs. Strong GEI limits both the prediction of genotype performance across variable environmental conditions and the identification of genomic regions associated with grain yield. Incorporating GEI into yield prediction models has been shown to improve prediction accuracy of yield; nevertheless, more work is needed to further understand this complex interaction across populations and environments. The main objectives of this study were to: 1) assess GEI in maize grain yield based on reaction norm models and predict hybrid performance across a gradient of environmental (EG) conditions and 2) perform a genome-wide association study (GWAS) and post-GWAS analyses for maize grain yield using data from 2014 to 2017 of the Genomes to Fields initiative hybrid trial. After quality control, 2,126 hybrids with genotypic and phenotypic data were assessed across 86 environments representing combinations of locations and years, although not all hybrids were evaluated in all environments. Heritability was greater in higher-yielding environments due to an increase in genetic variability in these environments in comparison to the low-yielding environments. GWAS was carried out for yield and five single nucleotide polymorphisms (SNPs) with the highest magnitude of effect were selected in each environment for follow-up analyses. Many candidate genes in proximity of selected SNPs have been previously reported with roles in stress response. Genomic prediction was performed to assess prediction accuracy of previously tested or untested hybrids in environments from a new growing season. Prediction accuracy was 0.34 for cross validation across years (CV0-Predicted EG) and 0.21 for cross validation across years with only untested hybrids (CV00-Predicted EG) when compared to Best Linear Unbiased Prediction (BLUPs) that did not utilize genotypic or environmental relationships. Prediction accuracy improved to 0.80 (CV0-Predicted EG) and 0.60 (CV00-Predicted EG) when compared to the whole-dataset model that used the genomic relationships and the environmental gradient of all environments in the study. These results identify regions of the genome for future selection to improve yield and a methodology to increase the number of hybrids evaluated across locations of a multi-environment trial through genomic prediction.

X-linked genes influence various complex traits in dairy cattle

BACKGROUND

The search for quantitative trait loci (QTL) affecting traits of interest in mammals is frequently limited to autosomes, with the X chromosome excluded because of its hemizygosity in males. This study aimed to assess the importance of the X chromosome in the genetic determinism of 11 complex traits related to milk production, milk composition, mastitis resistance, fertility, and stature in 236,496 cows from three major French dairy breeds (Holstein, Montbéliarde, and Normande) and three breeds of regional importance (Abondance, Tarentaise, and Vosgienne).

RESULTS

Estimates of the proportions of heritability due to autosomes and X chromosome (h²_X) were consistent among breeds. On average over the 11 traits, h²_X=0.008 and the X chromosome explained ~ 3.5% of total genetic variance. GWAS was performed within-breed at the sequence level (~ 200,000 genetic variants) and then combined in a meta-analysis. QTL were identified for most breeds and traits analyzed, with the exception of Tarentaise and Vosgienne and two fertility traits. Overall, 3, 74, 59, and 71 QTL were identified in Abondance, Montbéliarde, Normande, and Holstein, respectively, and most were associated with the most-heritable traits (milk traits and stature). The meta-analyses, which assessed a total of 157 QTL for the different traits, highlighted new QTL and refined the positions of some QTL found in the within-breed analyses. Altogether, our analyses identified a number of functional candidate genes, with the most notable being GPC3, MBNL3, HS6ST2, and DMD for dairy traits; TMEM164, ACSL4, ENOX2, HTR2C, AMOT, and IRAK1 for udder health; MAMLD1 and COL4A6 for fertility; and NRK, ESX1, GPR50, GPC3, and GPC4 for stature.

CONCLUSIONS

This study demonstrates the importance of the X chromosome in the genetic determinism of complex traits in dairy cattle and highlights new functional candidate genes and variants for these traits. These results could potentially be extended to other species as many X-linked genes are shared among mammals.

Genomic-based genetic parameters for resilience across lactations in North American Holstein cattle based on variability in daily milk yield records

Considering the increasing challenges imposed by climate change and the need to improve animal welfare, breeding more resilient animals capable of better coping with environmental disturbances is of paramount importance. In dairy cattle, resilience can be evaluated by measuring the longitudinal occurrences of abnormal daily milk yield throughout lactation. Aiming to estimate genetic parameters for dairy cattle resilience, we collected 5,643,193 daily milk yield records on automatic milking systems (milking robots) and milking parlors across 21,350 lactations 1 to 3 of 11,787 North American Holstein cows. All cows were genotyped with 62,029 SNPs. After determining the best fitting models for each of the 3 lactations, daily milk yield residuals were used to derive 4 resilience indicators: weighted occurrence frequency of yield perturbations (wfPert), accumulated milk losses of yield perturbations (dPert), and log-transformed variance (LnVar) and lag-1 autocorrelation (r_auto) of daily yield residuals. The indicator LnVar presented the highest heritability estimates (±standard error), ranging from 0.13 ± 0.01 in lactation 1 to 0.15 ± 0.02 in lactation 2; the other 3 indicators had relatively lower heritabilities across the 3 lactations (0.01-0.06). Based on bivariate analyses of each resilience indicator across lactations, stronger genetic correlations were observed between lactations 2 and 3 (0.88-0.96) than between lactations 1 and 2 or 3 (0.34-0.88) for dPert, LnVar, and r_auto. For the pairwise comparisons of different resilience indicators within each lactation, dPert had the strongest genetic correlations with wfPert (0.64) and r_auto (0.53) in lactation 1, whereas the correlations in lactations 2 and 3 were more variable and showed relatively high standard errors. The genetic correlation results indicated that different resilience indicators across lactations might capture additional biological mechanisms and should be considered as different traits in genetic evaluations. We also observed favorable genetic correlations of these resilience indicators with longevity and Net Merit index, but further biological validation of these resilience indicators is needed. In conclusion, this study provided genetic parameter estimates for different resilience indicators derived from daily milk yields across the first 3 lactations in Holstein cattle, which will be useful when potentially incorporating these traits in dairy cattle breeding schemes.

A dynamic concept of animal welfare: The role of appetitive and adverse internal and external factors and the animal’s ability to adapt to them

Animal welfare is a multifaceted issue that can be approached from different viewpoints, depending on human interests, ethical assumptions, and culture. To properly assess, safeguard and promote animal welfare, concepts are needed to serve as guidelines in any context the animal is kept in. Several different welfare concepts have been developed during the last half decade. The Five Freedoms concept has provided the basis for developing animal welfare assessment to date, and the Five Domains concept has guided those responsible for safeguarding animal welfare, while the Quality of Life concept focuses on how the individual perceives its own welfare state. This study proposes a modified and extended version of an earlier animal welfare concept - the Dynamic Animal Welfare Concept (DAWCon). Based on the adaptability of the animal, and taking the importance of positive emotional states and the dynamic nature of animal welfare into account, an individual animal is likely in a positive welfare state when it is mentally and physically capable and possesses the ability and opportunity to react adequately to sporadic or lasting appetitive and adverse internal and external stimuli, events, and conditions. Adequate reactions are elements of an animal’s normal behavior. They allow the animal to cope with and adapt to the demands of the (prevailing) environmental circumstances, enabling it to reach a state that it perceives as positive, i.e., that evokes positive emotions. This paper describes the role of internal as well as external factors in influencing welfare, each of which exerts their effects in a sporadic or lasting manner. Behavior is highlighted as a crucial read-out parameter. As most animals under human care are selected for certain traits that may affect their behavioral repertoire it is crucial to have thorough ethograms, i.e., a catalogue of specific behaviors of the species/strain/breed under study. DAWCon highlights aspects that need to be addressed when assessing welfare and may stimulate future research questions.

A Mutation in Endogenous saRNA miR-23a Influences Granulosa Cells Response to Oxidative Stress

Phenotypes are the result of the interaction between the gene and the environment, so the response of individuals with different genotypes to an environment is variable. Here, we reported that a mutation in miR-23a influences granulosa cells (GCs) response to oxidative stress, a common mechanism of environmental factors affecting female reproduction. We showed that nuclear miR-23a is a pro-apoptotic miRNA in porcine GCs through the activation of the transcription and function of NORHA, a long non-coding RNA (lncRNA) induces GC apoptosis and responses to oxidative stress. Mechanistically, miR-23a acts as an endogenous small activating RNA (saRNA) to alter histone modifications of the NORHA promoter through the direct binding to its core promoter. A C > T mutation was identified at −398 nt of the miR-23a core promoter, which created a novel binding site for the transcription factor SMAD4 and recruited the transcription repressor SMAD4 to inhibit miR-23a transcription and function in GCs. Notably, g.−398C > T mutation in the miR-23a promoter reduced GCs response to oxidative stress. In addition, g.−398C > T mutation was significantly associated with sow fertility traits. In short, our findings preliminarily revealed the genetic basis of individual differences in the response to oxidative stress from the perspective of a single mutation and identified miR-23a as a candidate gene for the environmental adaptation to oxidative stress.

Use of Milk Infrared Spectral Data as Environmental Covariates in Genomic Prediction Models for Production Traits in Canadian Holstein

The purpose of this study was to provide a procedure for the inclusion of milk spectral information into genomic prediction models. Spectral data were considered a set of covariates, in addition to genomic covariates. Milk yield and somatic cell score were used as traits to investigate. A cross-validation was employed, making a distinction for predicting new individuals' performance under known environments, known individuals' performance under new environments, and new individuals' performance under new environments. We found an advantage of including spectral data as environmental covariates when the genomic predictions had to be extrapolated to new environments. This was valid for both observed and, even more, unobserved families (genotypes). Overall, prediction accuracy was larger for milk yield than somatic cell score. Fourier-transformed infrared spectral data can be used as a source of information for the calculation of the 'environmental coordinates' of a given farm in a given time, extrapolating predictions to new environments. This procedure could serve as an example of integration of genomic and phenomic data. This could help using spectral data for traits that present poor predictability at the phenotypic level, such as disease incidence and behavior traits. The strength of the model is the ability to couple genomic with high-throughput phenomic information.

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.