Reaction-norm animal model analysis of average daily gain heat tolerance in purebred Duroc pigs

The present study aimed to genetically improve growth performance under high-heat environments by specifically designing a reaction-norm animal model (RNAM) for purebred Duroc pigs in Japan. A total of 54,750 records of average daily gain (ADG) measured for pigs reared at four farms in different prefectures were analyzed. Estimated maximum daily temperatures at the respective farm locations were used to calculate the average cumulative thermal load (TL). The TL values served as an indicator of high-heat environments for pigs. The plausible cumulative period length and threshold temperature for calculating TL were determined to be 8 weeks until just before shipping and 25°C, respectively. Variance components were estimated via RNAM analysis using TL as a linear covariate. The estimated additive genetic variances under both responsive and non-responsive to TL were found to be significant. Moreover, the estimated heritability of ADG ranged from 0.38 to 0.73 for TL values of 0-8. These results suggest that the RNAM developed holds the potential for improving the genetic ability of growth under high-heat environments in pigs.

Single-step genomic predictions for heat tolerance of production yields in US Holsteins and Jerseys

The physiological stress caused by excessive heat affects dairy cattle health and production. This study sought to investigate the effect of heat stress on test-day yields in US Holstein and Jersey cows and develop single-step genomic predictions to identify heat tolerant animals. Data included 12.8 million and 2.1 million test-day records, respectively, for 923,026 Holstein and 153,710 Jersey cows in 27 US states. From 2015 through 2021, test-day records from the first 5 lactations included milk, fat, and protein yields (kg). Cow records were included if they had at least 5 test-day records per lactation. Heat stress was quantified by analyzing the effect of a 5-d hourly average temperature-humidity index (THI5d¯) on observed test-day yields. Using a multiple trait repeatability model, a heat threshold (THI threshold) was determined fowr each breed based on the point that the average adjusted yields started to decrease, which was 69 for Holsteins and 72 for Jerseys. An additive genetic component of general production and heat tolerance production were estimated using a multiple trait reaction norm model and single-step genomic BLUP methodology. Random effects were regressed on a function of 5-d hourly average (THI5d¯) and THI threshold. The proportion of test-day records that occurred on or above the respective heat thresholds was 15% for Holstein and 10% for Jersey. Heritability of milk, fat, and protein yields under heat stress for Holsteins increased, with a small standard error, indicating that the additive genetic component for heat tolerance of these traits was observed. This was not as evident in Jersey traits. For Jersey, the permanent environment explained the same or more of the variation in fat and protein yield under heat stress indicating that nongenetic factors may determine heat tolerance for these Jersey traits. Correlations between the general genetic merit of production (in the absence of heat stress) and heat tolerance genetic merit of production traits were moderate in strength and negative. This indicated that selecting for general genetic merit without consideration of heat tolerance genetic merit of production may result in less favorable performance in hot and humid climates. A general genomic estimated breeding value for genetic merit and a heat tolerance genomic estimated breeding value were calculated for each animal. This study contributes to the investigation of the impact of heat stress on US dairy cattle production yields and offers a basis for the implementation of genomic selection. The results indicate that genomic selection for heat tolerance of production yields is possible for US Holsteins and Jerseys, but a study to validate the genomic predictions should be explored.

Genome-wide scans for selection signatures in indigenous pigs revealed candidate genes relating to heat tolerance

Heat stress is a major problem that constrains pig productivity. Understanding and identifying adaptation to heat stress has been the focus of recent studies, and the identification of genome-wide selection signatures can provide insights into the mechanisms of environmental adaptation. Here, we generated whole-genome re-sequencing data from six Chinese indigenous pig populations to identify genomic regions with selection signatures related to heat tolerance using multiple methods: three methods for intra-population analyses (Integrated Haplotype Score, Runs of Homozygosity and Nucleotide diversity Analysis) and three methods for inter-population analyses (Fixation index (F_ST), Cross-population Composite Likelihood Ratio and Cross-population Extended Haplotype Homozygosity). In total, 1 966 796 single nucleotide polymorphisms were identified in this study. Genetic structure analyses and F_ST indicated differentiation among these breeds. Based on information on the location environment, the six breeds were divided into heat and cold groups. By combining two or more approaches for selection signatures, outlier signals in overlapping regions were identified as candidate selection regions. A total of 163 candidate genes were identified, of which, 29 were associated with heat stress injury and anti-inflammatory effects. These candidate genes were further associated with 78 Gene Ontology functional terms and 30 Kyoto Encyclopedia of Genes and Genomes pathways in enrichment analysis (P < 0.05). Some of these have clear relevance to heat resistance, such as the AMPK signalling pathway and the mTOR signalling pathway. The results improve our understanding of the selection mechanisms responsible for heat resistance in pigs and provide new insights of introgression in heat adaptation.

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.

Systematic review of animal-based indicators to measure thermal, social, and immune-related stress in pigs

The intense nature of pig production has increased the animals’ exposure to stressful conditions, which may be detrimental to their welfare and productivity. Some of the most common sources of stress in pigs are extreme thermal conditions (thermal stress), density and mixing during housing (social stress), or exposure to pathogens and other microorganisms that may challenge their immune system (immune-related stress). The stress response can be monitored based on the animals’ coping mechanisms, as a result of specific environmental, social, and health conditions. These animal-based indicators may support decision making to maintain animal welfare and productivity. The present study aimed to systematically review animal-based indicators of social, thermal, and immune-related stresses in farmed pigs, and the methods used to monitor them. Peer-reviewed scientific literature related to pig production was collected using three online search engines: ScienceDirect, Scopus, and PubMed. The manuscripts selected were grouped based on the indicators measured during the study. According to our results, body temperature measured with a rectal thermometer was the most commonly utilized method for the evaluation of thermal stress in pigs (87.62%), as described in 144 studies. Of the 197 studies that evaluated social stress, aggressive behavior was the most frequently-used indicator (81.81%). Of the 535 publications examined regarding immune-related stress, cytokine concentration in blood samples was the most widely used indicator (80.1%). Information about the methods used to measure animal-based indicators is discussed in terms of validity, reliability, and feasibility. Additionally, the introduction and wide spreading of alternative, less invasive methods with which to measure animal-based indicators, such as cortisol in saliva, skin temperature and respiratory rate via infrared thermography, and various animal welfare threats via vocalization analysis are highlighted. The information reviewed was used to discuss the feasible and most reliable methods with which to monitor the impact of relevant stressors commonly presented by intense production systems on the welfare of farmed pigs.

Effect of ambient temperature on average daily gain of pigs evaluated using public weather data and a plateau-linear regression model

We performed a plateau-linear regression model analysis of the average daily gain (ADG) of pigs on daily average temperature at the end of performance testing (T). Records for performance testing between 30 kg and 105 kg of 2268 purebred Duroc pigs raised at the National Livestock Breeding Center Miyazaki Station were used. Off-farm ambient temperatures were measured at the nearest Automated Meteorological Data Acquisition System station at Kobayashi, Miyazaki (Kobayashi station). A plateau-linear regression equation was obtained in which ADG decreased by 12.6 g for every 1°C when T > 21.1°C. We calculated the expected age in day at the end of testing (D105) using the regression equation obtained and T observed at the Kobayashi station in 2020. The number of days that D105 was prolonged due to higher T was 125 days, corresponding to approximately one third of the year. These results could contribute to planning and management of stable pork production in response to heat in Japan.

Validation of single-step genomic predictions using the linear regression method for milk yield and heat tolerance in a Thai-Holstein population

Background and Aim:

Genomic selection improves accuracy and decreases the generation interval, increasing the selection response. This study was conducted to assess the benefits of using single-step genomic best linear unbiased prediction (ssGBLUP) for genomic evaluations of milk yield and heat tolerance in Thai-Holstein cows and to test the value of old phenotypic data to maintain the accuracy of predictions.

Materials and Methods:

The dataset included 104,150 milk yield records collected from 1999 to 2018 from 15,380 cows. The pedigree contained 33,799 animals born between 1944 and 2016, of which 882 were genotyped. Analyses were performed with and without genomic information using ssGBLUP and BLUP, respectively. Statistics for bias, dispersion, the ratio of accuracies, and the accuracy of estimated breeding values were calculated using the linear regression (LR) method. A partial dataset excluded the phenotypes of the last generation, and 66 bulls were identified as validation individuals.

Results:

Bias was considerable for BLUP (0.44) but negligible (−0.04) for ssGBLUP; dispersion was similar for both techniques (0.84 vs. 1.06 for BLUP and ssGBLUP, respectively). The ratio of accuracies was 0.33 for BLUP and 0.97 for ssGBLUP, indicating more stable predictions for ssGBLUP. The accuracy of predictions was 0.18 for BLUP and 0.36 for ssGBLUP. Excluding the first 10 years of phenotypic data (i.e., 1999-2008) decreased the accuracy to 0.09 for BLUP and 0.32 for ssGBLUP. Genomic information doubled the accuracy and increased the persistence of genomic estimated breeding values when old phenotypes were removed.

Conclusion:

The LR method is useful for estimating accuracies and bias in complex models. When the population size is small, old data are useful, and even a small amount of genomic information can substantially improve the accuracy. The effect of heat stress on first parity milk yield is small.

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.

Genetic Parameters for Tolerance to Heat Stress in Crossbred Swine Carcass Traits

Data for loin and backfat depth, as well as carcass growth of 126,051 three-way crossbred pigs raised between 2015 and 2019, were combined with climate records of air temperature, relative humidity, and temperature-humidity index. Environmental covariates with the largest impact on the studied traits were incorporated in a random regression model that also included genomic information. Genetic control of tolerance to heat stress and the presence of genotype by environment interaction were detected. Its magnitude was more substantial for loin depth and carcass growth, but all the traits studied showed a different impact of heat stress and different magnitude of genotype by environment interaction. For backfat depth, heritability was larger under comfortable conditions (no heat stress), as compared to heat stress conditions. Genetic correlations between extreme values of environmental conditions were lower (∼0.5 to negative) for growth and loin depth. Based on the solutions obtained from the model, sires were ranked on their breeding value for general performance and tolerance to heat stress. Antagonism between overall performance and tolerance to heat stress was moderate. Still, the models tested can provide valuable information to identify genetic material that is resilient and can perform equally when environmental conditions change. Overall, the results obtained from this study suggest the existence of genotype by environment interaction for carcass traits, as a possible genetic contributor to heat tolerance in swine.

Effects of sire line, birth weight and sex on growth performance and carcass traits of crossbred pigs under standardized environmental conditions

A variety of available terminal sire lines makes the choice of terminal sire line complex for the pig producer. Higher birth weights are important for subsequent growth performance and selection for this trait is also necessary in sire lines. The aim was to investigate the effect of sire line, birth weight and gender on growth performance, carcass traits and meat quality. In total 3844 crossbred pigs from Camborough Pig Improvement Company (PIC) dams matched with either a Synthetic (A) or Piétrain (B) sire line were used. Pigs from line A grew faster (p<0.01), showed higher feed intake (p<0.01) and reached a higher final body weight (p≤0.01), but they had a similar efficiency (p=0.179). Leaner carcasses and heavier primal cuts (p<0.001) were observed in pigs from line B. Carcasses from pigs sired by line A had higher meat quality (p<0.001). Males had a higher growth rate (p≤0.05) but had a poorer feed efficiency (p<0.01). Heavier birth weight pigs and females had leaner, higher value carcasses with heavier primal cuts (p<0.001) compared to middle and low birth weight females or males. Sire line by sex interactions was significant for growth (p≤0.05) and carcass traits (p<0.001). Interaction between sire line and birth weight classes were only detected for loin depth (p<0.01). Line A is preferable if the numbers of fatting pigs per fattening place and year should be improved, and line B is an option to increase leanness and carcass primal cuts.

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.

Large-Scale Phenotyping of Livestock Welfare in Commercial Production Systems: A New Frontier in Animal Breeding

Genomic breeding programs have been paramount in improving the rates of genetic progress of productive efficiency traits in livestock. Such improvement has been accompanied by the intensification of production systems, use of a wider range of precision technologies in routine management practices, and high-throughput phenotyping. Simultaneously, a greater public awareness of animal welfare has influenced livestock producers to place more emphasis on welfare relative to production traits. Therefore, management practices and breeding technologies in livestock have been developed in recent years to enhance animal welfare. In particular, genomic selection can be used to improve livestock social behavior, resilience to disease and other stress factors, and ease habituation to production system changes. The main requirements for including novel behavioral and welfare traits in genomic breeding schemes are: (1) to identify traits that represent the biological mechanisms of the industry breeding goals; (2) the availability of individual phenotypic records measured on a large number of animals (ideally with genomic information); (3) the derived traits are heritable, biologically meaningful, repeatable, and (ideally) not highly correlated with other traits already included in the selection indexes; and (4) genomic information is available for a large number of individuals (or genetically close individuals) with phenotypic records. In this review, we (1) describe a potential route for development of novel welfare indicator traits (using ideal phenotypes) for both genetic and genomic selection schemes; (2) summarize key indicator variables of livestock behavior and welfare, including a detailed assessment of thermal stress in livestock; (3) describe the primary statistical and bioinformatic methods available for large-scale data analyses of animal welfare; and (4) identify major advancements, challenges, and opportunities to generate high-throughput and large-scale datasets to enable genetic and genomic selection for improved welfare in livestock. A wide variety of novel welfare indicator traits can be derived from information captured by modern technology such as sensors, automatic feeding systems, milking robots, activity monitors, video cameras, and indirect biomarkers at the cellular and physiological levels. The development of novel traits coupled with genomic selection schemes for improved welfare in livestock can be feasible and optimized based on recently developed (or developing) technologies. Efficient implementation of genetic and genomic selection for improved animal welfare also requires the integration of a multitude of scientific fields such as cell and molecular biology, neuroscience, immunology, stress physiology, computer science, engineering, quantitative genomics, and bioinformatics.

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.

Comparing Alternative Single-Step GBLUP Approaches and Training Population Designs for Genomic Evaluation of Crossbred Animals

As crossbreeding is extensively used in some livestock species, we aimed to evaluate the performance of single-step GBLUP (ssGBLUP) and weighted ssGBLUP (WssGBLUP) methods to predict Genomic Estimated Breeding Values (GEBVs) of crossbred animals. Different training population scenarios were evaluated: (SC1) ssGBLUP based on a single-trait model considering purebred and crossbred animals in a joint training population; (SC2) ssGBLUP based on a multiple-trait model to enable considering phenotypes recorded in purebred and crossbred training animals as different traits; (SC3) WssGBLUP based on a single-trait model considering purebred and crossbred animals jointly in the training population (both populations were used for SNP weights' estimation); (SC4) WssGBLUP based on a single-trait model considering only purebred animals in the training population (crossbred population only used for SNP weights' estimation); (SC5) WssGBLUP based on a single-trait model and the training population characterized by purebred animals (purebred population used for SNP weights' estimation). A complex trait was simulated assuming alternative genetic architectures. Different scaling factors to blend the inverse of the genomic (G -1) and pedigree (A 22 - 1 ) relationship matrices were also tested. The predictive performance of each scenario was evaluated based on the validation accuracy and regression coefficient. The genetic correlations across simulated populations in the different scenarios ranged from moderate to high (0.71-0.99). The scenario mimicking a completely polygenic trait (h Q T L 2 = 0) yielded the lowest validation accuracy (0.12; for SC3 and SC4). The simulated scenarios assuming 4,500 QTLs affecting the trait andh Q T L 2 = h 2 resulted in the greatest GEBV accuracies (0.47; for SC1 and SC2). The regression coefficients ranged from 0.28 (for SC3 assuming polygenic effect) to 1.27 (for SC2 considering 4,500 QTLs). In general, SC3 and SC5 resulted in inflated GEBVs, whereas other scenarios yielded deflated GEBVs. The scaling factors used to combine G -1 andA 22 - 1 had a small influence on the validation accuracies, but a greater effect on the regression coefficients. Due to the complexity of multiple-trait models and WssGBLUP analyses, and a similar predictive performance across the methods evaluated, SC1 is recommended for genomic evaluation in crossbred populations with similar genetic structures [moderate-to-high (0.71-0.99) genetic correlations between purebred and crossbred populations].

Improving accuracy of direct and maternal genetic effects in genomic evaluations using pooled boar semen: a simulation study1

Pooling semen of multiple boars is commonly used in swine production systems. Compared with single boar systems, this technique changes family structure creating maternal half-sib families. The aim of this simulation study was to investigate how pooling semen affects the accuracy of estimating direct and maternal effects for individual piglet birth weight, in purebred pigs. Different scenarios of pooling semen were simulated by allowing the same female to mate from 1 to 6 boars, per insemination, whereas litter size was kept constant (N = 12). In each pooled boar scenario, genomic information was used to construct either the genomic relationship matrix (G) or to reconstruct pedigree in addition to G. Genotypes were generated for 60,000 SNPs evenly distributed across 18 autosomes. From the 5 simulated generations, only animals from generations 3 to 5 were genotyped (N = 36,000). Direct and maternal true breeding values (TBV) were computed as the sum of the effects of the 1,080 QTLs. Phenotypes were constructed as the sum of direct TBV, maternal TBV, an overall mean of 1.25 kg, and a residual effect. The simulated heritabilities for direct and maternal effects were 0.056 and 0.19, respectively, and the genetic correlation between both effects was -0.25. All simulations were replicated 5 times. Variance components and direct and maternal heritability were estimated using average information REML. Predictions were computed via pedigree-based BLUP and single-step genomic BLUP (ssGBLUP). Genotyped littermates in the last generation were used for validation. Prediction accuracies were calculated as correlations between EBV and TBV for direct (accdirect) and maternal (accmat) effects. When boars were known, accdirect were 0.21 (1 boar) and 0.26 (6 boars) for BLUP, whereas for ssGBLUP, they were 0.38 (1 boar) and 0.43 (6 boars). When boars were unknown, accdirect was lower in BLUP but similar in ssGBLUP. For the scenario with known boars, accmat was 0.58 and 0.63 for 1 and 6 boars, respectively, under ssGBLUP. For unknown boars, accmat was 0.63 for 2 boars and 0.62 for 6 boars in ssGBLUP. In general, accdirect and accmat were lower in the single-boar scenario compared with pooled semen scenarios, indicating that a half-sib structure is more adequate to estimate direct and maternal effects. Using pooled semen from multiple boars can help us to improve accuracy of predicting maternal and direct effects when maternal half-sib families are larger than 2.

Whole-genome mapping of quantitative trait loci and accuracy of genomic predictions for resistance to columnaris disease in two rainbow trout breeding populations

Background

Columnaris disease (CD) is an emerging problem for the rainbow trout aquaculture industry in the US. The objectives of this study were to: (1) identify common genomic regions that explain a large proportion of the additive genetic variance for resistance to CD in two rainbow trout (Oncorhynchus mykiss) populations; and (2) estimate the gains in prediction accuracy when genomic information is used to evaluate the genetic potential of survival to columnaris infection in each population.

Methods

Two aquaculture populations were investigated: the National Center for Cool and Cold Water Aquaculture (NCCCWA) odd-year line and the Troutlodge, Inc., May odd-year (TLUM) nucleus breeding population. Fish that survived to 21 days post-immersion challenge were recorded as resistant. Single nucleotide polymorphism (SNP) genotypes were available for 1185 and 1137 fish from NCCCWA and TLUM, respectively. SNP effects and variances were estimated using the weighted single-step genomic best linear unbiased prediction (BLUP) for genome-wide association. Genomic regions that explained more than 1% of the additive genetic variance were considered to be associated with resistance to CD. Predictive ability was calculated in a fivefold cross-validation scheme and using a linear regression method.

Results

Validation on adjusted phenotypes provided a prediction accuracy close to zero, due to the binary nature of the trait. Using breeding values computed from the complete data as benchmark improved prediction accuracy of genomic models by about 40% compared to the pedigree-based BLUP. Fourteen windows located on six chromosomes were associated with resistance to CD in the NCCCWA population, of which two windows on chromosome Omy 17 jointly explained more than 10% of the additive genetic variance. Twenty-six windows located on 13 chromosomes were associated with resistance to CD in the TLUM population. Only four associated genomic regions overlapped with quantitative trait loci (QTL) between both populations.

Conclusions

Our results suggest that genome-wide selection for resistance to CD in rainbow trout has greater potential than selection for a few target genomic regions that were found to be associated to resistance to CD due to the polygenic architecture of this trait, and because the QTL associated with resistance to CD are not sufficiently informative for selection decisions across populations.

Electronic supplementary material

The online version of this article (10.1186/s12711-019-0484-4) contains supplementary material, which is available to authorized users.