Genome-wide association studies for additive and dominance effects for body composition traits in commercial crossbred Piétrain pigs.

Fat depth (FD) and muscle depth (MD) are economically important traits and used to estimate carcass lean content (LMP), which is one of the main breeding objectives in pig breeding programmes. We assessed the genetic architectures of body composition traits for additive and dominance effects in commercial crossbred Piétrain pigs using both 50 K array and sequence genotypes. We first performed a genome-wide association study (GWAS) using single-marker association analysis with a false discovery rate of 0.1. Then, we estimated the additive and dominance effects of the most significant variant in the quantitative trait loci (QTL) regions. It was investigated whether the use of whole-genome sequence (WGS) will improve the QTL detection (both additive and dominance) with a higher power compared with lower density SNP arrays. Our results showed that more QTL regions were detected by WGS compared with 50 K array (n = 54 vs. n = 17). Of the novel associated regions associated with FD and LMP and detected by WGS, the most pronounced peak was on SSC13, situated at ~116-118, 121-127 and 129-134 Mbp. Additionally, we found that only additive effects contributed to the genetic architecture of the analysed traits and no significant dominance effects were found for the tested SNPs at QTL regions, regardless of panel density. The associated SNPs are located in or near several relevant candidate genes. Of these genes, GABRR2, GALR1, RNGTT, CDH20 and MC4R have been previously reported as being associated with fat deposition traits. However, the genes on SSC1 (ZNF292, ORC3, CNR1, SRSF12, MDN1, TSHZ1, RELCH and RNF152) and SSC18 (TTC26 and KIAA1549) have not been reported previously to our best knowledge. Our current findings provide insights into the genomic regions influencing composition traits in Piétrain pigs.

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.

Applying multi-omics data to study the genetic background of bovine respiratory disease infection in feedlot crossbred cattle.

Bovine respiratory disease (BRD) is the most common and costly infectious disease affecting the wellbeing and productivity of beef cattle in North America. BRD is a complex disease whose development is dependent on environmental factors and host genetics. Due to the polymicrobial nature of BRD, our understanding of the genetic and molecular mechanisms underlying the disease is still limited. This knowledge would augment the development of better genetic/genomic selection strategies and more accurate diagnostic tools to reduce BRD prevalence. Therefore, this study aimed to utilize multi-omics data (genomics, transcriptomics, and metabolomics) analyses to study the genetic and molecular mechanisms of BRD infection. Blood samples of 143 cattle (80 BRD; 63 non-BRD animals) were collected for genotyping, RNA sequencing, and metabolite profiling. Firstly, a genome-wide association study (GWAS) was performed for BRD susceptibility using 207,038 SNPs. Two SNPs (Chr5:25858264 and BovineHD1800016801) were identified as associated (p-value <1 × 10-5) with BRD susceptibility. Secondly, differential gene expression between BRD and non-BRD animals was studied. At the significance threshold used (log2FC>2, logCPM>2, and FDR<0.01), 101 differentially expressed (DE) genes were identified. These DE genes significantly (p-value <0.05) enriched several immune responses related functions such as inflammatory response. Additionally, we performed expression quantitative trait loci (eQTL) analysis and identified 420 cis-eQTLs and 144 trans-eQTLs significantly (FDR <0.05) associated with the expression of DE genes. Interestingly, eQTL results indicated the most significant SNP (Chr5:25858264) identified via GWAS was a cis-eQTL for DE gene GPR84. This analysis also demonstrated that an important SNP (rs209419196) located in the promoter region of the DE gene BPI significantly influenced the expression of this gene. Finally, the abundance of 31 metabolites was significantly (FDR <0.05) different between BRD and non-BRD animals, and 17 of them showed correlations with multiple DE genes, which shed light on the interactions between immune response and metabolism. This study identified associations between genome, transcriptome, metabolome, and BRD phenotype of feedlot crossbred cattle. The findings may be useful for the development of genomic selection strategies for BRD susceptibility, and for the development of new diagnostic and therapeutic tools.

Imputation to whole-genome sequence and its use in genome-wide association studies for pork colour traits in crossbred and purebred pigs.

Imputed whole-genome sequence (WGS) has been proposed to improve genome-wide association studies (GWAS), since all causative mutations responsible for phenotypic variation are expected to be present in the data. This approach was applied on a large number of purebred (PB) and crossbred (CB) pigs for 18 pork color traits to evaluate the impact of using imputed WGS relative to medium-density marker panels. The traits included Minolta A*, B*, and L* for fat (FCOL), quadriceps femoris muscle (QFCOL), thawed loin muscle (TMCOL), fresh ham gluteus medius (GMCOL), ham iliopsoas muscle (ICOL), and longissimus dorsi muscle on the fresh loin (FMCOL). Sequence variants were imputed from a medium-density marker panel (61K for CBs and 50K for PBs) in all genotyped pigs using BeagleV5.0. We obtained high imputation accuracy (average of 0.97 for PBs and 0.91 for CBs). GWAS were conducted for three datasets: 954 CBs and 891 PBs, and the combined CBs and PBs. For most traits, no significant associations were detected, regardless of panel density or population type. However, quantitative trait loci (QTL) regions were only found for a few traits including TMCOL Minolta A* and GMCOL Minolta B* (CBs), FMCOL Minolta B*, FMCOL Minolta L*, and ICOL Minolta B* (PBs) and FMCOL Minolta A*, FMCOL Minolta B*, GMCOL Minolta B*, and ICOL Minolta B* (Combined dataset). More QTL regions were identified with WGS (n = 58) relative to medium-density marker panels (n = 22). Most of the QTL were linked to previously reported QTLs or candidate genes that have been previously reported to be associated with meat quality, pH and pork color; e.g., VIL1, PRKAG3, TTLL4, and SLC11A1, USP37. CTDSP1 gene on SSC15 has not been previously associated with meat color traits in pigs. The findings suggest any added value of WGS was only for detecting novel QTL regions when the sample size is sufficiently large as with the Combined dataset in this study. The percentage of phenotypic variance explained by the most significant SNPs also increased with WGS compared with medium-density panels. The results provide additional insights into identification of a number of candidate regions and genes for pork color traits in different pig populations.

Integrative analyses of genomic and metabolomic data reveal genetic mechanisms associated with carcass merit traits in beef cattle.

Improvement of carcass merit traits is a priority for the beef industry. Discovering DNA variants and genes associated with variation in these traits and understanding biological functions/processes underlying their associations are of paramount importance for more effective genetic improvement of carcass merit traits in beef cattle. This study integrates 10,488,742 imputed whole genome DNA variants, 31 plasma metabolites, and animal phenotypes to identify genes and biological functions/processes that are associated with carcass merit traits including hot carcass weight (HCW), rib eye area (REA), average backfat thickness (AFAT), lean meat yield (LMY), and carcass marbling score (CMAR) in a population of 493 crossbred beef cattle. Regression analyses were performed to identify plasma metabolites associated with the carcass merit traits, and the results showed that 4 (3-hydroxybutyric acid, acetic acid, citric acid, and choline), 6 (creatinine, L-glutamine, succinic acid, pyruvic acid, L-lactic acid, and 3-hydroxybutyric acid), 4 (fumaric acid, methanol, D-glucose, and glycerol), 2 (L-lactic acid and creatinine), and 5 (succinic acid, fumaric acid, lysine, glycine, and choline) plasma metabolites were significantly associated with HCW, REA, AFAT, LMY, and CMAR (P-value < 0.1), respectively. Combining the results of metabolome-genome wide association studies using the 10,488,742 imputed SNPs, 103, 160, 83, 43, and 109 candidate genes were identified as significantly associated with HCW, REA, AFAT, LMY, and CMAR (P-value < 1 × 10-5), respectively. By applying functional enrichment analyses for candidate genes of each trait, 26, 24, 26, 24, and 28 significant cellular and molecular functions were predicted for HCW, REA, AFAT, LMY, and CMAR, respectively. Among the five topmost significantly enriched biological functions for carcass merit traits, molecular transport and small molecule biochemistry were two top biological functions associated with all carcass merit traits. Lipid metabolism was the most significant biological function for LMY and CMAR and it was also the second and fourth highest biological function for REA and HCW, respectively. Candidate genes and enriched biological functions identified by the integrative analyses of metabolites with phenotypic traits and DNA variants could help interpret the results of previous genome-wide association studies for carcass merit traits. Our integrative study also revealed additional potential novel genes associated with these economically important traits. Therefore, our study improves understanding of the molecular and biological functions/processes that influence carcass merit traits, which could help develop strategies to enhance genomic prediction of carcass merit traits with incorporation of metabolomic data. Similarly, this information could guide management practices, such as nutritional interventions, with the purpose of boosting specific carcass merit traits.

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.

Breeding for disease resilience: opportunities to manage polymicrobial challenge and improve commercial performance in the pig industry.

Disease resilience, defined as an animal's ability to maintain productive performance in the face of infection, provides opportunities to manage the polymicrobial challenge common in pig production. Disease resilience can deliver a number of benefits, including more sustainable production as well as improved animal health and the potential for reduced antimicrobial use. However, little progress has been made to date in the application of disease resilience in breeding programs due to a number of factors, including (1) confusion around definitions of disease resilience and its component traits disease resistance and tolerance, and (2) the difficulty in characterizing such a complex trait consisting of multiple biological functions and dynamic elements of rates of response and recovery from infection. Accordingly, this review refines the definitions of disease resistance, tolerance, and resilience based on previous studies to help improve the understanding and application of these breeding goals and traits under different scenarios. We also describe and summarize results from a "natural disease challenge model" designed to provide inputs for selection of disease resilience. The next steps for managing polymicrobial challenges faced by the pig industry will include the development of large-scale multi-omics data, new phenotyping technologies, and mathematical and statistical methods adapted to these data. Genome editing to produce pigs resistant to major diseases may complement selection for disease resilience along with continued efforts in the more traditional areas of biosecurity, vaccination and treatment. Altogether genomic approaches provide exciting opportunities for the pig industry to overcome the challenges provided by hard-to-manage diseases as well as new environmental challenges associated with climate change.

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.

Improving Chicken Responses to Glycoconjugate Vaccination Against Campylobacter jejuni.

Campylobacter jejuni is a common cause of diarrheal disease worldwide. Human infection typically occurs through the ingestion of contaminated poultry products. We previously demonstrated that an attenuated Escherichia coli live vaccine strain expressing the C. jejuni N-glycan on its surface reduced the Campylobacter load in more than 50% of vaccinated leghorn and broiler birds to undetectable levels (responder birds), whereas the remainder of the animals was still colonized (non-responders). To understand the underlying mechanism, we conducted three vaccination and challenge studies using 135 broiler birds and found a similar responder/non-responder effect. Subsequent genome-wide association studies (GWAS), analyses of bird sex and levels of vaccine-induced IgY responses did not correlate with the responder versus non-responder phenotype. In contrast, antibodies isolated from responder birds displayed a higher Campylobacter-opsonophagocytic activity when compared to antisera from non-responder birds. No differences in the N-glycome of the sera could be detected, although minor changes in IgY glycosylation warrant further investigation. As reported before, the composition of the microbiota, particularly levels of OTU classified as Clostridium spp., Ruminococcaceae and Lachnospiraceae are associated with the response. Transplantation of the cecal microbiota of responder birds into new birds in combination with vaccination resulted in further increases in vaccine-induced antigen-specific IgY responses when compared to birds that did not receive microbiota transplants. Our work suggests that the IgY effector function and microbiota contribute to the efficacy of the E. coli live vaccine, information that could form the basis for the development of improved vaccines targeted at the elimination of C. jejuni from poultry.

Relationship between indirect genetic effects for growth, environmental enrichment, coping style and sex with the serum metabolome profile of pigs.

Including Indirect Genetic Effects (IGE) in breeding programs to reduce aggression in group housed animals has been proposed. However, the effect of selection for IGE for growth on animal metabolism and physiology is unknown. The purpose of this study was twofold: (1) To investigate the effects of this new breeding method along with two housing (barren and straw), coping style (high and low resisters) and sex (female and castrated males) options on the metabolome profile of pigs. (2) To identify and map biological processes associated with a regrouping test at 9 weeks of age. We used Nuclear Magnetic Resonance to quantify 49 serum metabolites at week 8, 9 and 22. Also, we quantified 3 catecholamines (tyramine, epinephrine, phenylethylamine) and serotonin and three water soluble vitamins (B2, B5 and B7). Overall, no significant differences were observed between negative and positive IGE animals. The magnitude of change (delta) of many metabolites as a response to the regrouping test was significantly affected by IGE, especially that of the amino acids (P < 0.05), being greater in positive IGE pigs. The regrouping test was associated with alteration in glycine, serine and threonine metabolism. In conclusion positive and negative IGE animals respond differently to the regrouping test.

Identification of candidate genes and enriched biological functions for feed efficiency traits by integrating plasma metabolites and imputed whole genome sequence variants in beef cattle.

BACKGROUND: Feed efficiency is one of the key determinants of beef industry profitability and sustainability. However, the cellular and molecular background behind feed efficiency is largely unknown. This study combines imputed whole genome DNA variants and 31 plasma metabolites to dissect genes and biological functions/processes that are associated with residual feed intake (RFI) and its component traits including daily dry matter intake (DMI), average daily gain (ADG), and metabolic body weight (MWT) in beef cattle. RESULTS: Regression analyses between feed efficiency traits and plasma metabolites in a population of 493 crossbred beef cattle identified 5 (L-valine, lysine, L-tyrosine, L-isoleucine, and L-leucine), 4 (lysine, L-lactic acid, L-tyrosine, and choline), 1 (citric acid), and 4 (L-glutamine, glycine, citric acid, and dimethyl sulfone) plasma metabolites associated with RFI, DMI, ADG, and MWT (P-value < 0.1), respectively. Combining the results of metabolome-genome wide association studies using 10,488,742 imputed SNPs, 40, 66, 15, and 40 unique candidate genes were identified as associated with RFI, DMI, ADG, and MWT (P-value < 1 × 10-5), respectively. These candidate genes were found to be involved in some key metabolic processes including metabolism of lipids, molecular transportation, cellular function and maintenance, cell morphology and biochemistry of small molecules. CONCLUSIONS: This study identified metabolites, candidate genes and enriched biological functions/processes associated with RFI and its component traits through the integrative analyses of metabolites with phenotypic traits and DNA variants. Our findings could enhance the understanding of biochemical mechanisms of feed efficiency traits and could lead to improvement of genomic prediction accuracy via incorporating metabolite data.

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.

Genetic parameters of drinking and feeding traits of wean-to-finish pigs under a polymicrobial natural disease challenge.

BACKGROUND: The pork industry faces unprecedented challenges from disease, which increases cost of production and use of antibiotics, and reduces production efficiency, carcass quality, and animal wellbeing. One solution is to improve the overall resilience of pigs to a broad array of common diseases through genetic selection. Behavioral changes in feeding and drinking are usually the very first clinical signs when animals are exposed to stressors such as disease. Changes in feeding and drinking behaviors in diseased pigs may reflect the way they cope with the challenge and, thus, could be used as indicator traits to select for disease resilience. The objectives of this study were to estimate genetic parameters of feeding and drinking traits for wean-to-finish pigs in a natural polymicrobial disease challenge model, to estimate genetic correlations of feeding and drinking traits with growth rate and clinical disease traits, and to develop indicator traits to select for disease resilience. RESULTS: In general, drinking traits had moderate to high estimates of heritability, especially average daily water dispensed, duration, and number of visits (0.44 to 0.58). Similar estimates were observed for corresponding feeding traits (0.35 to 0.51). Most genetic correlation estimates among drinking traits were moderate to high (0.30 to 0.92) and higher than among feeding traits (0 to 0.11). Compared to other drinking traits, water intake duration and number of visits had relatively stronger negative genetic correlation estimates with treatment rate and mortality, especially across the challenge nursery and finisher (- 0.39 and - 0.45 for treatment rate; - 0.20 and - 0.19 for mortality). CONCLUSION: Most of the recorded drinking and feeding traits under a severe disease challenge had moderate to high estimates of heritability, especially for feed or water intake duration and number of visits. Phenotypic and genetic correlations among the recorded feeding traits under disease were generally low but drinking traits showed high correlations with each other. Water intake duration and number of visits are potential indicator traits to select for disease resilience because of their high heritability and had moderate genetic correlations with treatment and mortality rates under severe disease.

Quantitative analysis of the blood transcriptome of young healthy pigs and its relationship with subsequent disease resilience.

BACKGROUND: Disease resilience, which is the ability of an animal to maintain performance under disease, is important for pigs in commercial herds, where they are exposed to various pathogens. Our objective was to investigate population-level gene expression profiles in the blood of 912 healthy F1 barrows at ~ 27 days of age for associations with performance and health before and after their exposure to a natural polymicrobial disease challenge at ~ 43 days of age. RESULTS: Most significant (q < 0.20) associations of the level of expression of individual genes in blood of young healthy pigs were identified for concurrent growth rate and subjective health scores prior to the challenge, and for mortality, a combined mortality-treatment trait, and feed conversion rate after the challenge. Gene set enrichment analyses revealed three groups of gene ontology biological process terms that were related to disease resilience: 1) immune and stress response-related terms were enriched among genes whose increased expression was unfavorably associated with both pre- and post-challenge traits, 2) heme-related terms were enriched among genes that had favorable associations with both pre- and post-challenge traits, and 3) terms related to protein localization and viral gene expression were enriched among genes that were associated with reduced performance and health traits after but not before the challenge. CONCLUSIONS: Gene expression profiles in blood from young healthy piglets provide insight into their performance when exposed to disease and other stressors. The expression of genes involved in stress response, heme metabolism, and baseline expression of host genes related to virus propagation were found to be associated with host response to disease.

Investigating the genetic architecture of disease resilience in pigs by genome-wide association studies of complete blood count traits collected from a natural disease challenge model.

BACKGROUND: Genetic improvement for disease resilience is anticipated to be a practical method to improve efficiency and profitability of the pig industry, as resilient pigs maintain a relatively undepressed level of performance in the face of infection. However, multiple biological functions are known to be involved in disease resilience and this complexity means that the genetic architecture of disease resilience remains largely unknown. Here, we conducted genome-wide association studies (GWAS) of 465,910 autosomal SNPs for complete blood count (CBC) traits that are important in an animal's disease response. The aim was to identify the genetic control of disease resilience. RESULTS: Univariate and multivariate single-step GWAS were performed on 15 CBC traits measured from the blood samples of 2743 crossbred (Landrace × Yorkshire) barrows drawn at 2-weeks before, and at 2 and 6-weeks after exposure to a polymicrobial infectious challenge. Overall, at a genome-wise false discovery rate of 0.05, five genomic regions located on Sus scrofa chromosome (SSC) 2, SSC4, SSC9, SSC10, and SSC12, were significantly associated with white blood cell traits in response to the polymicrobial challenge, and nine genomic regions on multiple chromosomes (SSC1, SSC4, SSC5, SSC6, SSC8, SSC9, SSC11, SSC12, SSC17) were significantly associated with red blood cell and platelet traits collected before and after exposure to the challenge. By functional enrichment analyses using Ingenuity Pathway Analysis (IPA) and literature review of previous CBC studies, candidate genes located nearby significant single-nucleotide polymorphisms were found to be involved in immune response, hematopoiesis, red blood cell morphology, and platelet aggregation. CONCLUSIONS: This study helps to improve our understanding of the genetic basis of CBC traits collected before and after exposure to a polymicrobial infectious challenge and provides a step forward to improve disease resilience.

Proliferation of peripheral blood mononuclear cells from healthy piglets after mitogen stimulation as indicators of disease resilience.

Disease resilience refers to the productivity of an animal under disease. Given the high biosecurity of pig nucleus herds, traits that can be measured on healthy pigs and that are genetically correlated with disease resilience, that is, genetic indicator traits, offer a strategy to select for disease resilience. Our objective was to evaluate mitogen stimulation assays (MSAs) on peripheral blood mononuclear cells (PBMCs) from young healthy pigs as genetic indicators for disease resilience. Data were from a natural disease challenge in which batches of 60 or 75 naïve Yorkshire × Landrace piglets were introduced every 3 wk into a continuous flow barn that was seeded with multiple diseases. In this environment, disease resilience traits, including growth, treatment, and mortality rates, were recorded on 3,136 pigs that were genotyped with a high-density marker panel. PBMCs from 882 of these pigs from 19 batches were isolated from whole blood collected prior to the disease challenge and stimulated with five mitogens: concanavalin A (ConA), phytohemagglutinin (PHA), pokeweed mitogen (PWM), lipopolysaccharide (LPS), and phorbol myristate acetate (PMA). The proliferation of cells was evaluated at 48, 72, and 96 h and compared with unstimulated samples (rest count). Heritabilities of cell proliferation were estimated using a model with batch as a fixed effect and covariates of entry age; rest count; complete blood count proportions of lymphocytes, monocytes, eosinophils, and basophils; and pen, litter, and animal genetics as random effects. Heritability estimates were highest for response to ConA (0.30 ± 0.09, 0.28 ± 0.10, 0.17 ± 0.10, and 0.25 ±0.10 at 48, 72, and 96 h after stimulation and for area under the curve across the three time points, respectively). Estimates were in a similar range for response to PHA and PMA but low for PWM and LPS. Responses to ConA, PHA, and PMA were moderately genetically correlated with several disease resilience traits and in the expected direction, but individual estimates were not significantly different from zero due to large SEs. In conclusion, although validation is needed, MSAss, in particular based on ConA, show promise as genetic indicator traits for disease resilience.

Genomic Heritability and Genome-Wide Association Studies of Plasma Metabolites in Crossbred Beef Cattle.

Metabolites, substrates or products of metabolic processes, are involved in many biological functions, such as energy metabolism, signaling, stimulatory and inhibitory effects on enzymes and immunological defense. Metabolomic phenotypes are influenced by combination of genetic and environmental effects allowing for metabolome-genome-wide association studies (mGWAS) as a powerful tool to investigate the relationship between these phenotypes and genetic variants. The objectives of this study were to estimate genomic heritability and perform mGWAS and in silico functional enrichment analyses for a set of plasma metabolites in Canadian crossbred beef cattle. Thirty-three plasma metabolites and 45,266 single nucleotide polymorphisms (SNPs) were available for 475 animals. Genomic heritability for all metabolites was estimated using genomic best linear unbiased prediction (GBLUP) including genomic breed composition as covariates in the model. A single-step GBLUP implemented in BLUPF90 programs was used to determine SNP P values and the proportion of genetic variance explained by SNP windows containing 10 consecutive SNPs. The top 10 SNP windows that explained the largest genetic variation for each metabolite were identified and mapped to detect corresponding candidate genes. Functional enrichment analyses were performed on metabolites and their candidate genes using the Ingenuity Pathway Analysis software. Eleven metabolites showed low to moderate heritability that ranged from 0.09 ± 0.15 to 0.36 ± 0.15, while heritability estimates for 22 metabolites were zero or negligible. This result indicates that while variations in 11 metabolites were due to genetic variants, the majority are largely influenced by environment. Three significant SNP associations were detected for betaine (rs109862186), L-alanine (rs81117935), and L-lactic acid (rs42009425) based on Bonferroni correction for multiple testing (family wise error rate <0.05). The SNP rs81117935 was found to be located within the Catenin Alpha 2 gene (CTNNA2) showing a possible association with the regulation of L-alanine concentration. Other candidate genes were identified based on additive genetic variance explained by SNP windows of 10 consecutive SNPs. The observed heritability estimates and the candidate genes and networks identified in this study will serve as baseline information for research into the utilization of plasma metabolites for genetic improvement of crossbred beef cattle.

Genomic Analysis of IgG Antibody Response to Common Pathogens in Commercial Sows in Health-Challenged Herds.

Losses due to infectious diseases are one of the main factors affecting productivity in the swine industry, motivating the investigation of disease resilience-related traits for genetic selection. However, these traits are not expected to be expressed in the nucleus herds, where selection is performed. One alternative is to use information from the commercial level to identify and select nucleus animals genetically superior for coping with pathogen challenges. In this study, we analyzed the genetic basis of antibody (Ab) response to common infectious pathogens in health-challenged commercial swine herds as potential indicator traits for disease resilience, including Ab response to influenza A virus of swine (IAV), Mycoplasma hyopneumoniae (MH), porcine circovirus (PCV2), and Actinobacillus pleuropneumoniae (APP; different serotypes). Ab response was measured in blood at entry into gilt rearing, post-acclimation (∼40 days after entering the commercial herd), and parities 1 and 2. Heritability estimates for Ab response to IAV, MH, and PCV2 ranged from 0 to 0.76. Ab response to APP ranged from 0 to 0.40. The genetic correlation (r G ) of Ab response to IAV with MH, PCV2, PRRSV, and APPmean (average Ab responses for all serotypes of APP) were positive (>0.29) at entry. APPmean was negatively correlated with PCV2 and MH at entry and parity 2 but positively correlated with MH at post-acclimation and parity 1. Genomic regions associated with Ab response to different APP serotypes were identified on 13 chromosomes. The region on chromosome 14 (2 Mb) was associated with several serotypes of APP, explaining up to 4.3% of the genetic variance of Ab to APP7 at entry. In general, genomic prediction accuracies for Ab response were low to moderate, except average Ab response to all infectious pathogens evaluated. These results suggest that genetic selection of Ab response in commercial sows is possible, but with variable success depending on the trait and the time-point of collection. Future work is needed to determine genetic correlations of Ab response with disease resilience, reproductive performance, and other production traits.

The genetic basis of natural antibody titers of young healthy pigs and relationships with disease resilience.

BACKGROUND: Disease resilience is the ability to maintain performance under pathogen exposure but is difficult to select for because breeding populations are raised under high health. Selection for resilience requires a trait that is heritable, easy to measure on healthy animals, and genetically correlated with resilience. Natural antibodies (NAb) are important parts of the innate immune system and are found to be heritable and associated with disease susceptibility in dairy cattle and poultry. Our objective was to investigate NAb and total IgG in blood of healthy, young pigs as potential indicator traits for disease resilience. RESULTS: Data were from Yorkshire x Landrace pigs, with IgG and IgM NAb (four antigens) and total IgG measured by ELISA in blood plasma collected ~ 1 week after weaning, prior to their exposure to a natural polymicrobial challenge. Heritability estimates were lower for IgG NAb (0.12 to 0.24, + 0.05) and for total IgG (0.19 + 0.05) than for IgM NAb (0.33 to 0.53, + 0.07) but maternal effects were larger for IgG NAb (0.41 to 0.52, + 0.03) and for total IgG (0.19 + 0.05) than for IgM NAb (0.00 to 0.10, + 0.04). Phenotypically, IgM NAb titers were moderately correlated with each other (average 0.60), as were IgG NAb titers (average 0.42), but correlations between IgM and IgG NAb titers were weak (average 0.09). Phenotypic correlations of total IgG were moderate with NAb IgG (average 0.46) but weak with NAb IgM (average 0.01). Estimates of genetic correlations among NAb showed similar patterns but with small SE, with estimates averaging 0.76 among IgG NAb, 0.63 among IgM NAb, 0.17 between IgG and IgM NAb, 0.64 between total IgG and IgG NAb, and 0.13 between total IgG and IgM NAb. Phenotypically, pigs that survived had slightly higher levels of NAb and total IgG than pigs that died. Genetically, higher levels of NAb tended to be associated with greater disease resilience based on lower mortality and fewer parenteral antibiotic treatments. Genome-wide association analyses for NAb titers identified several genomic regions, with several candidate genes for immune response. CONCLUSIONS: Levels of NAb in blood of healthy young piglets are heritable and potential genetic indicators of resilience to polymicrobial disease.

Genetic analysis of disease resilience in wean-to-finish pigs from a natural disease challenge model.

The objective was to estimate the genetic parameters of performance and resilience of growing pigs under disease. Data were from 3,139 Yorkshire × Landrace wean-to-finish pigs that were exposed to a natural polymicrobial disease challenge that was established by entering naturally infected animals into a nursery barn, targeting various viral and bacterial diseases. The challenge was maintained by entering batches of 60 or 75 healthy nursery pigs every 3 wk in a continuous flow system. Traits analyzed included average daily gain (ADG), feed intake (ADFI) and duration (ADFD); feed conversion ratio (FCR); residual feed intake (RFI); mortality (MOR); number of health treatments (TRT); health scores (HScore); carcass weight (CWT), back fat (CBF) and loin depth (CLD); dressing percentage (DRS); lean yield (LYLD); day-to-day variation in feed intake and duration (VARFI and VARDUR); and the proportion of off-feed days (OFFFI and OFFDUR). Analyses were performed by mixed linear models with genomic relationships. The resilience traits, such as TRT, MOR, and HScore, were lowly heritable (up to 0.15) but had high genetic correlations with each other. Performance traits, such as ADG, ADFI, ADFD, FCR, RFI, and carcass traits, were moderate to highly heritable (0.17 to 0.49). Heritabilities of resilience indicator traits such as OFF and VAR had low to moderate heritabilities (0.08 to 0.23) but were higher when based on duration vs. amount. ADFI had a low genetic correlation with ADFD (0.13). ADG in the challenge nursery had stronger negative genetic correlations with both TRT and MOR than ADG in the finisher (-0.37 to -0.74 vs. -0.15 to -0.56). ADFI and FCR had moderate negative (-0.21 to -0.39) and positive (0.34 to 0.49) genetic correlations, respectively, with TRT and MOR. ADFD and RFI had very low genetic correlations with TRT and MOR. CWT and DRS were moderately negatively correlated with TRT and MOR (-0.33 to -0.59). Resilience indicator traits based on feed intake or duration had moderate to high positive genetic correlations with TRT (0.18 to 0.81) and MOR (0.33 to 0.87). In conclusion, performance and resilience traits under a polymicrobial disease challenge are heritable and can be changed by selection. Phenotypes extracted from feed intake patterns can be used as genetic indicator traits for disease resilience. Most promising is day-to-day variation in intake duration, which had a sizeable heritability (0.23) and favorable genetic correlations with MOR (0.79) and treatment rate (0.20).