Current challenges and future of agricultural genomes to phenomes in the USA.

Dramatic improvements in measuring genetic variation across agriculturally relevant populations (genomics) must be matched by improvements in identifying and measuring relevant trait variation in such populations across many environments (phenomics). Identifying the most critical opportunities and challenges in genome to phenome (G2P) research is the focus of this paper. Previously (Genome Biol, 23(1):1-11, 2022), we laid out how Agricultural Genome to Phenome Initiative (AG2PI) will coordinate activities with USA federal government agencies expand public-private partnerships, and engage with external stakeholders to achieve a shared vision of future the AG2PI. Acting on this latter step, AG2PI organized the "Thinking Big: Visualizing the Future of AG2PI" two-day workshop held September 9-10, 2022, in Ames, Iowa, co-hosted with the United State Department of Agriculture's National Institute of Food and Agriculture (USDA NIFA). During the meeting, attendees were asked to use their experience and curiosity to review the current status of agricultural genome to phenome (AG2P) work and envision the future of the AG2P field. The topic summaries composing this paper are distilled from two 1.5-h small group discussions. Challenges and solutions identified across multiple topics at the workshop were explored. We end our discussion with a vision for the future of agricultural progress, identifying two areas of innovation needed: (1) innovate in genetic improvement methods development and evaluation and (2) innovate in agricultural research processes to solve societal problems. To address these needs, we then provide six specific goals that we recommend be implemented immediately in support of advancing AG2P research.

Identification of the genetic basis of sow pelvic organ prolapse.

Introduction: Pelvic organ prolapse (POP) is one contributor to recent increases in sow mortality that have been observed in some populations and environments, leading to financial losses and welfare concerns. Methods: With inconsistent previous reports, the objective here was to investigate the role of genetics on susceptibility to POP, using data on 30,429 purebred sows, of which 14,186 were genotyped (25K), collected from 2012 to 2022 in two US multiplier farms with a high POP incidence of 7.1% among culled and dead sows and ranging from 2% to 4% of all sows present by parity. Given the low incidence of POP for parities 1 and >6, only data from parities 2 to 6 were retained for analyses. Genetic analyses were conducted both across parities, using cull data (culled for POP versus another reason), and by parity, using farrowing data. (culled for POP versus culled for another reason or not culled). Results and Discussion: Estimates of heritability from univariate logit models on the underlying scale were 0.35 ± 0.02 for the across-parity analysis and ranged from 0.41 ± 0.03 in parity 2 to 0.15 ± 0.07 in parity 6 for the by-parity analyses. Estimates of genetic correlations of POP between parities based on bivariate linear models indicated a similar genetic basis of POP across parities but less similar with increasing distance between parities. Genome wide association analyses revealed six 1 Mb windows that explained more than 1% of the genetic variance in the across-parity data. Most regions were confirmed in several by-parity analyses. Functional analyses of the identified genomic regions showed a potential role of several genes on chromosomes 1, 3, 7, 10, 12, and 14 in susceptibility to POP, including the Estrogen Receptor gene. Gene set enrichment analyses showed that genomic regions that explained more variation for POP were enriched for several terms from custom transcriptome and gene ontology libraries. Conclusion: The influence of genetics on susceptibility to POP in this population and environment was confirmed and several candidate genes and biological processes were identified that can be targeted to better understand and mitigate the incidence of POP.

Genome wide association study of thyroid hormone levels following challenge with porcine reproductive and respiratory syndrome virus.

Introduction: Porcine reproductive and respiratory syndrome virus (PRRSV) causes respiratory disease in piglets and reproductive disease in sows. Piglet and fetal serum thyroid hormone (i.e., T3 and T4) levels decrease rapidly in response to Porcine reproductive and respiratory syndrome virus infection. However, the genetic control of T3 and T4 levels during infection is not completely understood. Our objective was to estimate genetic parameters and identify quantitative trait loci (QTL) for absolute T3 and/or T4 levels of piglets and fetuses challenged with Porcine reproductive and respiratory syndrome virus. Methods: Sera from 5-week-old pigs (N = 1792) at 11 days post inoculation (DPI) with Porcine reproductive and respiratory syndrome virus were assayed for T3 levels (piglet_T3). Sera from fetuses (N = 1,267) at 12 or 21 days post maternal inoculation (DPMI) with Porcine reproductive and respiratory syndrome virus of sows (N = 145) in late gestation were assayed for T3 (fetal_T3) and T4 (fetal_T4) levels. Animals were genotyped using 60 K Illumina or 650 K Affymetrix single nucleotide polymorphism (SNP) panels. Heritabilities, phenotypic correlations, and genetic correlations were estimated using ASREML; genome wide association studies were performed for each trait separately using Julia for Whole-genome Analysis Software (JWAS). Results: All three traits were low to moderately heritable (10%-16%). Phenotypic and genetic correlations of piglet_T3 levels with weight gain (0-42 DPI) were 0.26 ± 0.03 and 0.67 ± 0.14, respectively. Nine significant quantitative trait loci were identified for piglet_T3, on Sus scrofa chromosomes (SSC) 3, 4, 5, 6, 7, 14, 15, and 17, and collectively explaining 30% of the genetic variation (GV), with the largest quantitative trait loci identified on SSC5, explaining 15% of the genetic variation. Three significant quantitative trait loci were identified for fetal_T3 on SSC1 and SSC4, which collectively explained 10% of the genetic variation. Five significant quantitative trait loci were identified for fetal_T4 on SSC1, 6, 10, 13, and 15, which collectively explained 14% of the genetic variation. Several putative immune-related candidate genes were identified, including CD247, IRF8, and MAPK8. Discussion: Thyroid hormone levels following Porcine reproductive and respiratory syndrome virus infection were heritable and had positive genetic correlations with growth rate. Multiple quantitative trait loci with moderate effects were identified for T3 and T4 levels during challenge with Porcine reproductive and respiratory syndrome virus and candidate genes were identified, including several immune-related genes. These results advance our understanding of growth effects of both piglet and fetal response to Porcine reproductive and respiratory syndrome virus infection, revealing factors associated with genomic control of host resilience.

Plasma protein levels of young healthy pigs as indicators of disease resilience.

Selection for disease resilience, which refers to the ability of an animal to maintain performance when exposed to disease, can reduce the impact of infectious diseases. However, direct selection for disease resilience is challenging because nucleus herds must maintain a high health status. A possible solution is indirect selection of indicators of disease resilience. To search for such indicators, we conducted phenotypic and genetic quantitative analyses of the abundances of 377 proteins in plasma samples from 912 young and visually healthy pigs and their relationships with performance and subsequent disease resilience after natural exposure to a polymicrobial disease challenge. Abundances of 100 proteins were significantly heritable (false discovery rate (FDR) <0.10). The abundance of some proteins was or tended to be genetically correlated (rg) with disease resilience, including complement system proteins (rg = -0.24, FDR = 0.001) and IgG heavy chain proteins (rg = -0.68, FDR = 0.22). Gene set enrichment analyses (FDR < 0.2) based on phenotypic and genetic associations of protein abundances with subsequent disease resilience revealed many pathways related to the immune system that were unfavorably associated with subsequent disease resilience, especially the innate immune system. It was not possible to determine whether the observed levels of these proteins reflected baseline levels in these young and visually healthy pigs or were the result of a response to environmental disturbances that the pigs were exposed to before sample collection. Nevertheless, results show that, under these conditions, the abundance of proteins in some immune-related pathways can be used as phenotypic and genetic predictors of disease resilience and have the potential for use in pig breeding and management.

A challenge of selection for disease resilience is that it is difficult to directly select pigs that have greater resilience to multiple diseases in the healthy nucleus herd environment which is essential for breeding programs. A possible alternative is to select an indicator trait or marker that can be measured in a healthy setting, is heritable, and is associated with the genetics of disease resilience. In this study, we investigated plasma protein levels measured on young healthy pigs as indicator traits to select for disease resilience. For this purpose, we used plasma proteome data collected prior to the natural exposure of nursery pigs to multiple diseases, performed phenotypic and genetic quantitative analyses, and investigated their relationships with disease resilience. Our results suggest that plasma protein levels of young healthy pigs have the potential as biomarkers to select for disease resistance.

Exploration of plasma metabolite levels in healthy nursery pigs in response to environmental enrichment and disease resilience.

The purpose of this study was to explore plasma metabolite levels in young healthy pigs and their potential association with disease resilience and estimate genetic and phenotypic correlation with the change in lymphocyte concentration following disease challenge. Plasma samples were collected from 968 healthy nursery pigs over 15 batches at an average of 28 ± 3.23 d of age. Forty-four metabolites were identified and quantified by nuclear magnetic resonance. Pigs were then introduced into a natural disease challenge barn, and were classified into four groups based on the growth rate of each animal in the grow-to-finish phase (GFGR) and treatment rate (TR): resilient (RES), average (MID), susceptible (SUS), and dead (pigs that died before harvest). Blood samples were collected from all pigs before and 2 wk after disease challenge and complete blood count was determined. Environmental enrichment (inedible point source objects) was provided for half of the pigs in seven batches (N = 205) to evaluate its impact on resilience and metabolite concentrations. Concentration of all metabolites was affected by batch, while entry age affected the concentration of 16 metabolites. The concentration of creatinine was significantly lower for pigs classified as "dead" and "susceptible" when compared to "average" (P < 0.05). Pigs that received enrichment had significantly lower concentrations of six metabolites compared with pigs that did not receive enrichment (P ≤ 0.05). Both, group classification and enrichment affected metabolites that are involved in the same pathways of valine, leucine, and isoleucine biosynthesis and degradation. Resilient pigs had higher increase in lymphocyte concentration after disease challenge. The concentration of plasma l-α-aminobutyric acid was significantly negatively genetically correlated with the change in lymphocyte concentration following challenge. In conclusion, creatinine concentration in healthy nursery pigs was lower in pigs classified as susceptible or dead after disease challenge, whilst l-α-aminobutyric may be a genetic biomarker of lymphocyte response after pathogen exposure, and both deserve further investigation. Batch, entry age, and environmental enrichment were important factors affecting the concentration of metabolites and should be taken into consideration in future studies.

The focus of this study was to explore plasma metabolite levels in young healthy pigs and their potential association with health outcome classification following the exposure to a polymicrobial disease challenge. In addition, we explored the effect of the environmental enrichment on metabolite concentrations. Finally, we estimated genetic and phenotypic correlations between metabolites and the magnitude of change in lymphocytes levels following exposure to a polymicrobial disease challenge. We found that concentration of creatinine was lower in pigs that died before marketing, classified as "dead" and susceptible when compared to average group. This indicates that creatinine can be used as an early indicator of death and/or susceptibility of disease in pigs. Providing environmental enrichment affected the concentration of six metabolites and branched chain amino acids index. Such results would be very useful to design environmental enrichment strategies when pigs are challenged by disease in commercial farms. The magnitude of change in lymphocytes level was negatively genetically correlated with l-α-aminobutyric acid. This result indicates that l-αs-aminobutyric acid can be an early indicator of the magnitude of change in lymphocytes level. Such indicator can be collected from nucleus breeding herds in healthy animals and could provide an early biomarker of resilience.

Genetic Analyses of Response of Local Ghanaian Tanzanian Chicken Ecotypes to a Natural Challenge with Velogenic Newcastle Disease Virus.

Newcastle disease is a devastating poultry disease that often causes significant economic losses in poultry in the developing countries of Africa, Asia, as well as South and Central America. Velogenic Newcastle disease virus (NDV) outbreaks are associated with high mortalities, which can threaten household livelihoods, especially in the rural areas, and lead to loss of high-quality proteins in the form of meat and eggs, as well as household purchasing power. In this study, we exposed unvaccinated Ghanaian and Tanzanian chickens of six local ecotypes to velogenic NDV strains, measured NDV response traits, sequenced their DNA on a genotyping-by-sequencing platform, and performed variance component analyses. The collected phenotypes included: growth rates (pre- and post-exposure); lesion scores (gross lesion severity) in the trachea, proventriculus, intestine, and cecal tonsils; natural antibody levels; anti-NDV antibody levels at 7 days post exposure (dpe); tear and cloacal viral load at 2, 4, and 6 dpe; and survival time. Heritability estimates were low to moderate, ranging from 0.11 for average lesion scores to 0.36 for pre-exposure growth rate. Heritability estimates for survival time were 0.23 and 0.27 for the Tanzanian and Ghanaian ecotypes, respectively. Similar heritability estimates were observed when data were analyzed either separately or combined for the two countries. Survival time was genetically negatively correlated with lesion scores and with viral load. Results suggested that response to mesogenic or velogenic NDV of these local chicken ecotypes could be improved by selective breeding. Chickens that are more resilient to velogenic NDV can improve household livelihoods in developing countries.

Identification of circulating metabolites associated with wooden breast and white striping.

Current diagnostic methods for wooden breast and white striping, common breast muscle myopathies of modern commercial broiler chickens, rely on subjective examinations of the pectoralis major muscle, time-consuming microscopy, or expensive imaging technologies. Further research on these disorders would benefit from more quantitative and objective measures of disease severity that can be used in live birds. To this end, we utilized untargeted metabolomics alongside two statistical approaches to evaluate plasma metabolites associated with wooden breast and white striping in 250 male commercial broiler chickens. First, mixed linear modeling was employed to identify metabolites with a significant association with these muscle disorders and found 98 metabolites associated with wooden breast and 44 metabolites associated with white striping (q-value < 0.05). Second, a support vector machine was constructed using stepwise feature selection to determine the smallest subset of metabolites with the highest categorization accuracy for wooden breast. The final support vector machine achieved 94% accuracy using only 6 metabolites. The metabolite 3-methylhistidine, which is often used as an index of myofibrillar breakdown in skeletal muscle, was the top metabolite for both wooden breast and white striping in our mixed linear model and was also the metabolite with highest marginal prediction accuracy (82%) for wooden breast in our support vector machine. Overall, this study identified a candidate set of metabolites for an objective measure of wooden breast or white striping severity in live birds and expanded our understanding of these muscle disorders.

Quality of breeding value predictions from longitudinal analyses, with application to residual feed intake in pigs.

BACKGROUND: An important goal in animal breeding is to improve longitudinal traits. The objective of this study was to explore for longitudinal residual feed intake (RFI) data, which estimated breeding value (EBV), or combination of EBV, to use in a breeding program. Linear combinations of EBV (summarized breeding values, SBV) or phenotypes (summarized phenotypes) derived from the eigenvectors of the genetic covariance matrix over time were considered, and the linear regression method (LR method) was used to facilitate the evaluation of their prediction accuracy. RESULTS: Weekly feed intake, average daily gain, metabolic body weight, and backfat thickness measured on 2435 growing French Large White pigs over a 10-week period were analysed using a random regression model. In this population, the 544 dams of the phenotyped animals were genotyped. These dams did not have own phenotypes. The quality of the predictions of SBV and breeding values from summarized phenotypes of these females was evaluated. On average, predictions of SBV at the time of selection were unbiased, slightly over-dispersed and less accurate than those obtained with additional phenotypic information. The use of genomic information did not improve the quality of predictions. The use of summarized instead of longitudinal phenotypes resulted in predictions of breeding values of similar quality. CONCLUSIONS: For practical selection on longitudinal data, the results obtained with this specific design suggest that the use of summarized phenotypes could facilitate routine genetic evaluation of longitudinal traits.

Application of Bayesian genomic prediction methods to genome-wide association analyses.

BACKGROUND: Bayesian genomic prediction methods were developed to simultaneously fit all genotyped markers to a set of available phenotypes for prediction of breeding values for quantitative traits, allowing for differences in the genetic architecture (distribution of marker effects) of traits. These methods also provide a flexible and reliable framework for genome-wide association (GWA) studies. The objective here was to review developments in Bayesian hierarchical and variable selection models for GWA analyses. RESULTS: By fitting all genotyped markers simultaneously, Bayesian GWA methods implicitly account for population structure and the multiple-testing problem of classical single-marker GWA. Implemented using Markov chain Monte Carlo methods, Bayesian GWA methods allow for control of error rates using probabilities obtained from posterior distributions. Power of GWA studies using Bayesian methods can be enhanced by using informative priors based on previous association studies, gene expression analyses, or functional annotation information. Applied to multiple traits, Bayesian GWA analyses can give insight into pleiotropic effects by multi-trait, structural equation, or graphical models. Bayesian methods can also be used to combine genomic, transcriptomic, proteomic, and other -omics data to infer causal genotype to phenotype relationships and to suggest external interventions that can improve performance. CONCLUSIONS: Bayesian hierarchical and variable selection methods provide a unified and powerful framework for genomic prediction, GWA, integration of prior information, and integration of information from other -omics platforms to identify causal mutations for complex quantitative traits.

Response of three local chicken ecotypes of Ghana to lentogenic and velogenic Newcastle disease virus challenge.

This study was carried out to assess the response of three Ghanaian local chicken ecotypes to LaSota (lentogenic) and virulent field strains of Newcastle disease virus (NDV). Local chickens sampled from the Interior Savannah (IS), Forest (FO) and Coastal Savannah (CS) agro-ecological zones were bred and their offspring were challenged with LaSota NDV at 4 weeks of age. The LaSota challenge was replicated four times with different chicken groups. A total of 1438 chicks comprising 509 Coastal Savannah, 518 Forest and 411 Interior Savannah ecotypes were used. Pre- and post-challenge anti-NDV antibody titre levels were determined via ELISA assays. A second trial was conducted by introducing sick birds infected with virulent NDV to a flock of immunologically naïve chickens at 4 weeks old. Body weights were measured pre- and post-infection. Sex of the chickens was determined using a molecular method. In both trials, there was no significant difference among ecotypes in body weight and growth rate. In the LaSota trial, anti-NDV antibody titre did not differ by ecotype or sex. However, there was a positive linear relationship between body weight and antibody titre. In the velogenic NDV trial, survivability and lesion scores were similar among the three ecotypes. This study confirms that a relatively high dose of LaSota (NDV) challenge has no undesirable effect on Ghanaian local chicken ecotypes. All three Ghanaian local chicken ecotypes were susceptible to velogenic NDV challenge. Resistance to NDV by Ghanaian local chickens appears to be determined more by the individual's genetic makeup than by their ecotype.

Comparison of linkage disequilibrium estimated from genotypes versus haplotypes for crossbred populations.

BACKGROUND: Linkage disequilibrium (LD) is commonly measured based on the squared coefficient of correlation [Formula: see text] between the alleles at two loci that are carried by haplotypes. LD can also be estimated as the [Formula: see text] between unphased genotype dosage at two loci when the allele frequencies and inbreeding coefficients at both loci are identical for the parental lines. Here, we investigated whether [Formula: see text] for a crossbred population (F1) can be estimated using genotype data. The parental lines of the crossbred (F1) can be purebred or crossbred. METHODS: We approached this by first showing that inbreeding coefficients for an F1 crossbred population are negative, and typically differ in size between loci. Then, we proved that the expected [Formula: see text] computed from unphased genotype data is expected to be identical to the [Formula: see text] computed from haplotype data for an F1 crossbred population, regardless of the inbreeding coefficients at the two loci. Finally, we investigated the bias and precision of the [Formula: see text] estimated using unphased genotype versus haplotype data in stochastic simulation. RESULTS: Our findings show that estimates of [Formula: see text] based on haplotype and unphased genotype data are both unbiased for different combinations of allele frequencies, sample sizes (900, 1800, and 2700), and levels of LD. In general, for any allele frequency combination and [Formula: see text] value scenarios considered, and for both methods to estimate [Formula: see text], the precision of the estimates increased, and the bias of the estimates decreased as sample size increased, indicating that both estimators are consistent. For a given scenario, the [Formula: see text] estimates using haplotype data were more precise and less biased using haplotype data than using unphased genotype data. As sample size increased, the difference in precision and biasedness between the [Formula: see text] estimates using haplotype data and unphased genotype data decreased. CONCLUSIONS: Our theoretical derivations showed that estimates of LD between loci based on unphased genotypes and haplotypes in F1 crossbreds have identical expectations. Based on our simulation results, we conclude that the LD for an F1 crossbred population can be accurately estimated from unphased genotype data. The results also apply for other crosses (F2, F3, Fn, BC1, BC2, and BCn), as long as (selected) individuals from the two parental lines mate randomly.

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.

Predictions of the accuracy of genomic prediction: connecting R2, selection index theory, and Fisher information.

BACKGROUND: Deterministic predictions of the accuracy of genomic estimated breeding values (GEBV) when combining information sources have been developed based on selection index theory (SIT) and on Fisher information (FI). These two approaches have resulted in slightly different results when considering the combination of pedigree and genomic information. Here, we clarify this apparent contradiction, both for the combination of pedigree and genomic information and for the combination of subpopulations into a joint reference population. RESULTS: First, we show that existing expressions for the squared accuracy of GEBV can be understood as a proportion of the variance explained. Next, we show that the apparent discrepancy that has been observed between accuracies based on SIT vs. FI originated from two sources. First, the FI referred to the genetic component that is captured by the marker genotypes, rather than the full genetic component. Second, the common SIT-based derivations did not account for the increase in the accuracy of GEBV due to a reduction of the residual variance when combining information sources. The SIT and FI approaches are equivalent when these sources are accounted for. CONCLUSIONS: The squared accuracy of GEBV can be understood as a proportion of the variance explained. The SIT and FI approaches for combining information for GEBV are equivalent and provide identical accuracies when the underlying assumptions are equivalent.

Genome-wide association study of disease resilience traits from a natural polymicrobial disease challenge model in pigs identifies the importance of the major histocompatibility complex region.

Infectious diseases cause tremendous financial losses in the pork industry, emphasizing the importance of disease resilience, which is the ability of an animal to maintain performance under disease. Previously, a natural polymicrobial disease challenge model was established, in which pigs were challenged in the late nursery phase by multiple pathogens to maximize expression of genetic differences in disease resilience. Genetic analysis found that performance traits in this model, including growth rate, feed and water intake, and carcass traits, as well as clinical disease phenotypes, were heritable and could be selected for to increase disease resilience of pigs. The objectives of the current study were to identify genomic regions that are associated with disease resilience in this model, using genome-wide association studies and fine-mapping methods, and to use gene set enrichment analyses to determine whether genomic regions associated with disease resilience are enriched for previously published quantitative trait loci, functional pathways, and differentially expressed genes subject to physiological states. Multiple quantitative trait loci were detected for all recorded performance and clinical disease traits. The major histocompatibility complex region was found to explain substantial genetic variance for multiple traits, including for growth rate in the late nursery (12.8%) and finisher (2.7%), for several clinical disease traits (up to 2.7%), and for several feeding and drinking traits (up to 4%). Further fine mapping identified 4 quantitative trait loci in the major histocompatibility complex region for growth rate in the late nursery that spanned the subregions for class I, II, and III, with 1 single-nucleotide polymorphism in the major histocompatibility complex class I subregion capturing the largest effects, explaining 0.8-27.1% of genetic variance for growth rate and for multiple clinical disease traits. This single-nucleotide polymorphism was located in the enhancer of TRIM39 gene, which is involved in innate immune response. The major histocompatibility complex region was pleiotropic for growth rate in the late nursery and finisher, and for treatment and mortality rates. Growth rate in the late nursery showed strong negative genetic correlations in the major histocompatibility complex region with treatment or mortality rates (-0.62 to -0.85) and a strong positive genetic correlation with growth rate in the finisher (0.79). Gene set enrichment analyses found genomic regions associated with resilience phenotypes to be enriched for previously identified disease susceptibility and immune capacity quantitative trait loci, for genes that were differentially expressed following bacterial or virus infection and immune response, and for gene ontology terms related to immune and inflammatory response. In conclusion, the major histocompatibility complex and other quantitative trait loci that harbor immune-related genes were identified to be associated with disease resilience traits in a large-scale natural polymicrobial disease challenge. The major histocompatibility complex region was pleiotropic for growth rate under challenge and for clinical disease traits. Four quantitative trait loci were identified across the class I, II, and III subregions of the major histocompatibility complex for nursery growth rate under challenge, with 1 single-nucleotide polymorphism in the major histocompatibility complex class I subregion capturing the largest effects. The major histocompatibility complex and other quantitative trait loci identified play an important role in host response to infectious diseases and can be incorporated in selection to improve disease resilience, in particular the identified single-nucleotide polymorphism in the major histocompatibility complex class I subregion.

Multiple trait breeding programs with genotype-by-environment interactions based on reaction norms, with application to genetic improvement of disease resilience.

BACKGROUND: Genotype-by-environment interactions for a trait can be modeled using multiple-trait, i.e. character-state, models, that consider the phenotype as a different trait in each environment, or using reaction norm models based on a functional relationship, usually linear, between phenotype and a quantitative measure of the quality of the environment. The equivalence between character-state and reaction norm models has been demonstrated for a single trait. The objectives of this study were to extend the equivalence of the reaction norm and character-state models to a multiple-trait setting and to both genetic and environmental effects, and to illustrate the application of this equivalence to the design and optimization of breeding programs for disease resilience. METHODS: Equivalencies between reaction norm and character-state models for multiple-trait phenotypes were derived at the genetic and environmental levels, which demonstrates how multiple-trait reaction norm parameters can be derived from multiple-trait character state parameters. Methods were applied to optimize selection for a multiple-trait breeding goal in a target environment based on phenotypes collected in a healthy and disease-challenged environment, and to optimize the environment in which disease-challenge phenotypes should be collected. RESULTS AND CONCLUSIONS: The equivalence between multiple-trait reaction norm and multiple-trait character-state parameters allow genetic improvement for a multiple-trait breeding goal in a target environment to be optimized without recording phenotypes and estimating parameters for the target environment.

Further host-genomic characterization of total antibody response to PRRSV vaccination and its relationship with reproductive performance in commercial sows: genome-wide haplotype and zygosity analyses.

BACKGROUND: The possibility of using antibody response (S/P ratio) to PRRSV vaccination measured in crossbred commercial gilts as a genetic indicator for reproductive performance in vaccinated crossbred sows has motivated further studies of the genomic basis of this trait. In this study, we investigated the association of haplotypes and runs of homozygosity (ROH) and heterozygosity (ROHet) with S/P ratio and their impact on reproductive performance. RESULTS: There was no association (P-value ≥ 0.18) of S/P ratio with the percentage of ROH or ROHet, or with the percentage of heterozygosity across the whole genome or in the major histocompatibility complex (MHC) region. However, specific ROH and ROHet regions were significantly associated (P-value ≤ 0.01) with S/P ratio on chromosomes 1, 4, 5, 7, 10, 11, 13, and 17 but not (P-value ≥ 0.10) with reproductive performance. With the haplotype-based genome-wide association study (GWAS), additional genomic regions associated with S/P ratio were identified on chromosomes 4, 7, and 9. These regions harbor immune-related genes, such as SLA-DOB, TAP2, TAPBP, TMIGD3, and ADORA. Four haplotypes at the identified region on chromosome 7 were also associated with multiple reproductive traits. A haplotype significantly associated with S/P ratio that is located in the MHC region may be in stronger linkage disequilibrium (LD) with the quantitative trait loci (QTL) than the previously identified single nucleotide polymorphism (SNP) (H3GA0020505) given the larger estimate of genetic variance explained by the haplotype than by the SNP. CONCLUSIONS: Specific ROH and ROHet regions were significantly associated with S/P ratio. The haplotype-based GWAS identified novel QTL for S/P ratio on chromosomes 4, 7, and 9 and confirmed the presence of at least one QTL in the MHC region. The chromosome 7 region was also associated with reproductive performance. These results narrow the search for causal genes in this region and suggest SLA-DOB and TAP2 as potential candidate genes associated with S/P ratio on chromosome 7. These results provide additional opportunities for marker-assisted selection and genomic selection for S/P ratio as genetic indicator for litter size in commercial pig populations.

Thyroid hormone suppression in feeder pigs following polymicrobial or porcine reproductive and respiratory syndrome virus-2 challenge.

Thyroid hormones are powerful regulators of growth, development, and basal metabolic rate and can be dysregulated under conditions of severe stress or illness. To understand the role of these hormones in porcine disease response, serum samples were obtained from three batches of nursery-aged pigs (n = 208) exposed to a natural polymicrobial disease challenge with an array of bacterial and viral pathogens. Levels of total thyroxin (T4) and triiodothyronine (T3) assessed in sera by radioimmunoassay, decreased significantly by 14 days post-exposure (DPE). Levels of T3 partially rebounded by 48 DPE, while T4 levels remain depressed. Post-exposure T3 and T4 levels were positively correlated with acute and long-term average daily gain (ADG). Cross-sectional sampling of animals maintained at the high health source farms, showed no equivalent change in either hormone when managed under standard industrial conditions. To further elucidate the effect of porcine reproductive and respiratory syndrome virus (PRRSV)-infection on thyroid hormone levels, archived sera over 42 days post inoculation (DPI) from nursery pigs (N = 190) challenged with one of two PRRSV2 strains by the PRRS Host Genetics Consortium were similarly assessed, with animals selected in a two-by-two design, to investigate biological extremes in ADG and viral load (VL). All animals showed a similar decrease in both thyroid hormones reaching a minimum at 7 DPI and returning to near pre-challenge levels by 42 DPI. Post-challenge T3 and T4 levels were significantly greater in high ADG groups, with no significant association with VL or strain. The results of this study demonstrate porcine susceptibility to thyroid disruption in response to disease challenge and demonstrate a relationship between this response and growth performance.