Gut microbiome and resistome characterization of pigs treated with commonly used post-weaning diarrhea treatments.

BACKGROUND: The global burden of antimicrobial resistance demands additional measures to ensure the sustainable and conscious use of antimicrobials. For the swine industry, the post-weaning period is critical and for many years, antimicrobials have been the most effective strategy to control and treat post-weaning related infections. Among them, post-weaning diarrhea causes vast economic losses, as it severely compromises piglets' health and growth performance. In this study, 210 piglets were transferred from a farm with recurrent cases of post-weaning diarrhea to an experimental farm and divided into six different treatment groups to determine the effect of the different treatments on the growth performance and survival, the microbiome, and the resistome in a cross-sectional and longitudinal study. The different treatments included antimicrobials trimethoprim/sulfamethoxazole, colistin, and gentamicin, an oral commercial vaccine, a control with water acidification, and an untreated control. An extra group remained at the farm of origin following the implemented amoxicillin routine treatment. A total of 280 fecal samples from pigs at four different sampling times were selected for metagenomics: before weaning-treatment at the farm of origin, and three days, two weeks, and four weeks post-treatment. RESULTS: The control group with water acidification showed a reduced death risk in the survival analyses and non-significant differences in average daily weight gain in comparison to the antibiotic-treated groups. However, the growth-promoting effect among antibiotic-treated groups was demonstrated when comparing against the untreated control group at the experimental farm. After four weeks of treatment, diversity indexes revealed significantly decreased diversity for the untreated control and the group that remained at the farm of origin treated with amoxicillin. For this last group, impaired microbial diversity could be related to the continuous amoxicillin treatment carried out at the farm. Analysis of the resistome showed that both gentamicin and amoxicillin treatments significantly contributed to the emergence of resistance, while trimethoprim/sulphonamide and colistin did not, suggesting that different treatments contribute differently to the emergence of resistance. CONCLUSIONS: Overall, this shotgun longitudinal metagenomics analysis demonstrates that non-antibiotic alternatives, such as water acidification, can contribute to reducing the emergence of antimicrobial resistance without compromising pig growth performance and gut microbiome.

Identification of candidate regulatory genes for intramuscular fatty acid composition in pigs by transcriptome analysis.

BACKGROUND: Intramuscular fat (IMF) content and its fatty acid (FA) composition are typically controlled by several genes, each with a small effect. In the current study, to pinpoint candidate genes and putative regulators involved in FA composition, we performed a multivariate integrative analysis between intramuscular FA and transcriptome profiles of porcine longissimus dorsi (LD) muscle. We also carried out a combination of network, regulatory impact factor (RIF), in silico prediction of putative target genes, and functional analyses to better support the biological relevance of our findings. RESULTS: For this purpose, we used LD RNA-Seq and intramuscular FA composition profiles of 129 Iberian × Duroc backcrossed pigs. We identified 378 correlated variables (13 FA and 365 genes), including six FA (C20:4n-6, C18:2n-6, C20:3n-6, C18:1n-9, C18:0, and C16:1n-7) that were among the most interconnected variables in the predicted network. The detected FA-correlated genes include genes involved in lipid and/or carbohydrate metabolism or in regulation of IMF deposition (e.g., ADIPOQ, CHUK, CYCS, CYP4B1, DLD, ELOVL6, FBP1, G0S2, GCLC, HMGCR, IDH3A, LEP, LGALS12, LPIN1, PLIN1, PNPLA8, PPP1R1B, SDR16C5, SFRP5, SOD3, SNW1, and TFRC), meat quality (GALNT15, GOT1, MDH1, NEU3, PDHA1, SDHD, and UNC93A), and transport (e.g., EXOC7 and SLC44A2). Functional analysis highlighted 54 over-represented gene ontology terms, including well-known biological processes and pathways that regulate lipid and carbohydrate metabolism. RIF analysis suggested a pivotal role for six transcription factors (CARHSP1, LBX1, MAFA, PAX7, SIX5, and TADA2A) as putative regulators of gene expression and intramuscular FA composition. Based on in silico prediction, we identified putative target genes for these six regulators. Among these, TADA2A and CARHSP1 had extreme RIF scores and present novel regulators in pigs. In addition, the expression of TADA2A correlated (either positively or negatively) with C20:4n-6, C18:2n-6, C20:3n-6, C18:1n-9, and that of CARHSP1 correlated (positively) with the C16:1n-7 lipokine. We also found that these two transcription factors share target genes that are involved in lipid metabolism (e.g., GOT1, PLIN1, and TFRC). CONCLUSIONS: This integrative analysis of muscle transcriptome and intramuscular FA profile revealed valuable information about key candidate genes and potential regulators for FA and lipid metabolism in pigs, among which some transcription factors are proposed to control gene expression and modulate FA composition differences.

Recent Advances in Enteric Methane Mitigation and the Long Road to Sustainable Ruminant Production.

Mitigation of methane emission, a potent greenhouse gas, is a worldwide priority to limit global warming. A substantial part of anthropogenic methane is emitted by the livestock sector, as methane is a normal product of ruminant digestion. We present the latest developments and challenges ahead of the main efficient mitigation strategies of enteric methane production in ruminants. Numerous mitigation strategies have been developed in the last decades, from dietary manipulation and breeding to targeting of methanogens, the microbes that produce methane. The most recent advances focus on specific inhibition of key enzymes involved in methanogenesis. But these inhibitors, although efficient, are not affordable and not adapted to the extensive farming systems prevalent in low- and middle-income countries. Effective global mitigation of methane emissions from livestock should be based not only on scientific progress but also on the feasibility and accessibility of mitigation strategies.

Correction: Mixed management in growing and finishing pigs: Differences between gender and their impacts on behavior, growth performance, and physiological parameters.

[This corrects the article DOI: 10.1371/journal.pone.0284481.].

Mutations on a conserved distal enhancer in the porcine C-reactive protein gene impair its expression in liver.

C-reactive protein (CRP) is an evolutionary highly conserved protein. Like humans, CRP acts as a major acute phase protein in pigs. While CRP regulatory mechanisms have been extensively studied in humans, little is known about the molecular mechanisms that control pig CRP gene expression. The main goal of the present work was to study the regulatory mechanisms and identify functional genetic variants regulating CRP gene expression and CRP blood levels in pigs. The characterization of the porcine CRP proximal promoter region revealed a high level of conservation with both cow and human promoters, sharing binding sites for transcription factors required for CRP expression. Through genome-wide association studies and fine mapping, the most associated variants with both mRNA and protein CRP levels were localized in a genomic region 39.3 kb upstream of CRP. Further study of the region revealed a highly conserved putative enhancer that contains binding sites for several transcriptional regulators such as STAT3, NF-kB or C/EBP-ß. Luciferase reporter assays showed the necessity of this enhancer-promoter interaction for the acute phase induction of CRP expression in liver, where differences in the enhancer sequences significantly modified CRP activity. The associated polymorphisms disrupted the putative binding sites for HNF4α and FOXA2 transcription factors. The high correlation between HNF4α and CRP expression levels suggest the participation of HNF4α in the regulatory mechanism of porcine CRP expression through the modification of its binding site in liver. Our findings determine, for the first time, the relevance of a distal regulatory element essential for the acute phase induction of porcine CRP in liver and identify functional polymorphisms that can be included in pig breeding programs to improve immunocompetence.

Copy Number Variation on ABCC2-DNMBP Loci Affects the Diversity and Composition of the Fecal Microbiota in Pigs.

Genetic variation in the pig genome partially modulates the composition of porcine gut microbial communities. Previous studies have been focused on the association between single nucleotide polymorphisms (SNPs) and the gut microbiota, but little is known about the relationship between structural variants and fecal microbial traits. The main goal of this study was to explore the association between porcine genome copy number variants (CNVs) and the diversity and composition of pig fecal microbiota. For this purpose, we used whole-genome sequencing data to undertake a comprehensive identification of CNVs followed by a genome-wide association analysis between the estimated CNV status and the fecal bacterial diversity in a commercial Duroc pig population. A CNV predicted as gain (DUP) partially harboring ABCC2-DNMBP loci was associated with richness (P = 5.41 × 10-5, false discovery rate [FDR] = 0.022) and Shannon α-diversity (P = 1.42 × 10-4, FDR = 0.057). The in silico predicted gain of copies was validated by real-time quantitative PCR (qPCR), and its segregation, and positive association with the richness and Shannon α-diversity of the porcine fecal bacterial ecosystem was confirmed in an unrelated F1 (Duroc × Iberian) cross. Our results advise the relevance of considering the role of host-genome structural variants as potential modulators of microbial ecosystems and suggest the ABCC2-DNMBP CNV as a host-genetic factor for the modulation of the diversity and composition of the fecal microbiota in pigs. IMPORTANCE A better understanding of the environmental and host factors modulating gut microbiomes is a topic of greatest interest. Recent evidence suggests that genetic variation in the pig genome partially controls the composition of porcine gut microbiota. However, since previous studies have been focused on the association between single nucleotide polymorphisms and the fecal microbiota, little is known about the relationship between other sources of genetic variation, like the structural variants and microbial traits. Here, we identified, experimentally validated, and replicated in an independent population a positive link between the gain of copies of ABCC2-DNMBP loci and the diversity and composition of pig fecal microbiota. Our results advise the relevance of considering the role of host-genome structural variants as putative modulators of microbial ecosystems and open the possibility of implementing novel holobiont-based management strategies in breeding programs for the simultaneous improvement of microbial traits and host performance.

On the holobiont 'predictome' of immunocompetence in pigs.

BACKGROUND: Gut microbial composition plays an important role in numerous traits, including immune response. Integration of host genomic information with microbiome data is a natural step in the prediction of complex traits, although methods to optimize this are still largely unexplored. In this paper, we assess the impact of different modelling strategies on the predictive capacity for six porcine immunocompetence traits when both genotype and microbiota data are available. METHODS: We used phenotypic data on six immunity traits and the relative abundance of gut bacterial communities on 400 Duroc pigs that were genotyped for 70 k SNPs. We compared the predictive accuracy, defined as the correlation between predicted and observed phenotypes, of a wide catalogue of models: reproducing kernel Hilbert space (RKHS), Bayes C, and an ensemble method, using a range of priors and microbial clustering strategies. Combined (holobiont) models that include both genotype and microbiome data were compared with partial models that use one source of variation only. RESULTS: Overall, holobiont models performed better than partial models. Host genotype was especially relevant for predicting adaptive immunity traits (i.e., concentration of immunoglobulins M and G), whereas microbial composition was important for predicting innate immunity traits (i.e., concentration of haptoglobin and C-reactive protein and lymphocyte phagocytic capacity). None of the models was uniformly best across all traits. We observed a greater variability in predictive accuracies across models when microbiability (the variance explained by the microbiome) was high. Clustering microbial abundances did not necessarily increase predictive accuracy. CONCLUSIONS: Gut microbiota information is useful for predicting immunocompetence traits, especially those related to innate immunity. Modelling microbiome abundances deserves special attention when microbiability is high. Clustering microbial data for prediction is not recommended by default.

Mixed management in growing and finishing pigs: Differences between gender and their impacts on behavior, growth performance, and physiological parameters.

Mixing, a common management strategy used to regroup pigs, has been reported to impair individual performance and affect pig welfare because of the establishment of a new social hierarchy after regrouping. In this study we aimed to determine whether mixing management (non-mixed vs. mixed) and gender (gilts vs. barrows) affect the social and non-social behavior, performance, and physiological parameters of pigs. A total of 96 growing pigs (48 barrows and 48 females) were separated into two treatments: control (CT)-pigs that were mixed once during the growing-finishing period; and social stress (SS)-pigs that were mixed thrice during the growing-finishing period. We recorded social and non-social behaviors, injury score, performance, and physiological parameters during the experimental period. Data were grouped by the period, based on each mix performed, and overall values. The statistical analysis performed considered gender and treatment. For treatment, during period-II and III, the SS group presented the highest frequency of agonistic interactions (AI), stayed longer lying laterally (LL) and sternly (LS), and explored more enrichment material (ER) than the CT group. Furthermore, SS pigs presented the highest injury score in the ear, head, and middle and posterior regions. Compared to the females, the barrows spent more time at the electronic feed station and initiated most of the agonistic interactions during period-II, and they presented a higher injury score for the ear and head regions during period-III. In conclusion, repeated regrouping significantly affected social and feeding behavior without severely altering performance and physiological parameters. Furthermore, different patterns of social and feeding behavior, agonistic interactions, and injury scores between barrows and females were observed. This study provides an understanding of the impact of mixing management and gender differences on pigs, and this knowledge can be used to improve swine productivity and welfare.

Genetic architecture of innate and adaptive immune cells in pigs.

Pig industry is facing new challenges that make necessary to reorient breeding programs to produce more robust and resilient pig populations. The aim of the present work was to study the genetic determinism of lymphocyte subpopulations in the peripheral blood of pigs and identify genomic regions and biomarkers associated to them. For this purpose, we stained peripheral blood mononuclear cells to measure ten immune-cell-related traits including the relative abundance of different populations of lymphocytes, the proportions of CD4+ T cells and CD8+ T cells, and the ratio of CD4+/CD8+ T cells from 391 healthy Duroc piglets aged 8 weeks. Medium to high heritabilities were observed for the ten immune-cell-related traits and significant genetic correlations were obtained between the proportion of some lymphocytes populations. A genome-wide association study pointed out 32 SNPs located at four chromosomal regions on pig chromosomes SSC3, SSC5, SSC8, and SSCX as significantly associated to T-helper cells, memory T-helper cells and γδ T cells. Several genes previously identified in human association studies for the same or related traits were located in the associated regions, and were proposed as candidate genes to explain the variation of T cell populations such as CD4, CD8A, CD8B, KLRC2, RMND5A and VPS24. The transcriptome analysis of whole blood samples from animals with extreme proportions of γδ T, T-helper and memory T-helper cells identified differentially expressed genes (CAPG, TCF7L1, KLRD1 and CD4) located into the associated regions. In addition, differentially expressed genes specific of different T cells subpopulations were identified such as SOX13 and WC1 genes for γδ T cells. Our results enhance the knowledge about the genetic control of lymphocyte traits that could be considered to optimize the induction of immune responses to vaccines against pathogens. Furthermore, they open the possibility of applying effective selection programs for improving immunocompetence in pigs and support the use of the pig as a very reliable human biomedical model.

Hierarchy Establishment in Growing Finishing Pigs: Impacts on Behavior, Growth Performance, and Physiological Parameters.

In recent years, an increased number of studies have dealt with the analysis of social dominance related to animal behavior, physiology, and performance. This study aimed to investigate whether hierarchical ranking affects the coping style, non-social behavior during open field and novel object tests, performance, and physiological parameters of pigs. A total of 48 growing pigs (24 barrows and 24 females) were mixed three times during the growing-finishing period. The social and non-social behaviors of pigs were directly noted, and three behavioral tests were performed during the experimental period. Performance and physiological parameters were also recorded. Statistical analysis considered hierarchical classification (dominant vs. intermediary vs. subordinate) and p-values ≤ 0.05 were considered significant. After three regroupings, the pigs in different hierarchical classifications showed no change in hair cortisol values and open-field and novel object tests. Mean corpuscular hemoglobin concentration and leukocyte values increased in intermediary pigs, and the lowest counts were found in pigs classified as dominants. Furthermore, dominant pigs visited the feeder more but spent shorter time there compared to subordinate and intermediary pigs. Our results suggest that hierarchical classification influenced feeding behavior and physiological parameters without affecting cortisol values and growth performance, demonstrating a possible compensation skill.

Global analysis of the association between pig muscle fatty acid composition and gene expression using RNA-Seq.

Fatty acids (FAs) play an essential role as mediators of cell signaling and signal transduction, affecting metabolic homeostasis and determining meat quality in pigs. However, FAs are transformed by the action of several genes, such as those encoding desaturases and elongases of FAs in lipogenic tissues. The aim of the current work was to identify candidate genes, biological processes, and pathways involved in the modulation of intramuscular FA profile from longissimus dorsi muscle. FA profile by gas chromatography of methyl esters and gene expression by RNA-Seq were determined in 129 Iberian × Duroc backcrossed pigs. An association analysis between the muscle transcriptome and its FA profile was performed, followed by a concordance and functional analysis. Overall, a list of well-known (e.g., PLIN1, LEP, ELOVL6, SC5D, NCOA2, ACSL1, MDH1, LPL, LGALS12, TFRC, GOT1, and FBP1) and novel (e.g., TRARG1, TANK, ENSSSCG00000011196, and ENSSSCG00000038429) candidate genes was identified, either in association with specific or several FA traits. Likewise, several of these genes belong to biological processes and pathways linked to energy, lipid, and carbohydrate metabolism, which seem determinants in the modulation of FA compositions. This study can contribute to elucidate the complex relationship between gene expression and FA profile in pig muscle.

Disentangling the causal relationship between rabbit growth and cecal microbiota through structural equation models.

BACKGROUND: The effect of the cecal microbiome on growth of rabbits that were fed under different regimes has been studied previously. However, the term "effect" carries a causal meaning that can be confounded because of potential genetic associations between the microbiome and production traits. Structural equation models (SEM) can help disentangle such a complex interplay by decomposing the effect on a production trait into direct host genetics effects and indirect host genetic effects that are exerted through microbiota effects. These indirect effects can be estimated via structural coefficients that measure the effect of the microbiota on growth while the effects of the host genetics are kept constant. In this study, we applied the SEM approach to infer causal relationships between the cecal microbiota and growth of rabbits fed under ad libitum (ADGAL) or restricted feeding (ADGR). RESULTS: We identified structural coefficients that are statistically different from 0 for 138 of the 946 operational taxonomic units (OTU) analyzed. However, only 15 and 38 of these 138 OTU had an effect greater than 0.2 phenotypic standard deviations (SD) on ADGAL and ADGR, respectively. Many of these OTU had a negative effect on both traits. The largest effects on ADGR were exerted by an OTU that is taxonomically assigned to the Desulfovibrio genus (- 1.929 g/d, CSS-normalized OTU units) and by an OTU that belongs to the Ruminococcaceae family (1.859 g/d, CSS-normalized OTU units). For ADGAL, the largest effect was from OTU that belong to the S24-7 family (- 1.907 g/d, CSS-normalized OTU units). In general, OTU that had a substantial effect had low to moderate estimates of heritability. CONCLUSIONS: Disentangling how direct and indirect effects act on production traits is relevant to fully describe the processes of mediation but also to understand how these traits change before considering the application of an external intervention aimed at changing a given microbial composition by blocking/promoting the presence of a particular microorganism.

Modeling microRNA-driven post-transcriptional regulation using exon-intron split analysis in pigs.

The contribution of microRNAs (miRNAs) to mRNA post-transcriptional regulation has often been explored by the post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes in mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Göttingen minipigs. Exon-intron split analysis of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of miRNAs upregulated in fed pigs. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, exon-intron split analysis provides a valuable view of the regulation of gene expression, complementary to canonical differential expression analyses.

Supplementation of xylo-oligosaccharides to suckling piglets promotes the growth of fiber-degrading gut bacterial populations during the lactation and nursery periods.

Modulating early-life microbial colonization through xylo-oligosacharides (XOS) supplementation represents an opportunity to accelerate the establishment of fiber-degrading microbial populations and improve intestinal health. Ninety piglets from 15 litters were orally administered once a day from d7 to d27 of lactation with either 5 mL of water (CON) or 5 mL of a solution containing 30 to 60 mg of XOS (XOS). Supplementation ceased at weaning (d28) when all piglets were fed the same commercial pre-starter diet. Growth performance did not differ between treatments during the experimental period (d7 to d40). Piglet's fecal microbiota (n = 30) shifted significantly from the end of lactation (d27) to nursery period (d40) exhibiting an increase in microbial alpha diversity. Animals supplemented with XOS showed higher richness and abundance of fiber-degrading bacteria and short-chain fatty acid (SCFA) production at d27 and d40. Additionally, the predicted abundance of the pyruvate to butanoate fermentation pathway was increased in the XOS group at d40. These results show that supplementation of XOS to lactating piglets promotes fiber-degrading bacterial populations in their hindgut. Moreover, differences observed in the nursery period suggest that XOS can influence the microbiota in the long-term.

Use of Bayes factors to evaluate the effects of host genetics, litter and cage on the rabbit cecal microbiota.

BACKGROUND: The rabbit cecum hosts and interacts with a complex microbial ecosystem that contributes to the variation of traits of economic interest. Although the influence of host genetics on microbial diversity and specific microbial taxa has been studied in several species (e.g., humans, pigs, or cattle), it has not been investigated in rabbits. Using a Bayes factor approach, the aim of this study was to dissect the effects of host genetics, litter and cage on 984 microbial traits that are representative of the rabbit microbiota. RESULTS: Analysis of 16S rDNA sequences of cecal microbiota from 425 rabbits resulted in the relative abundances of 29 genera, 951 operational taxonomic units (OTU), and four microbial alpha-diversity indices. Each of these microbial traits was adjusted with mixed linear and zero-inflated Poisson (ZIP) models, which all included additive genetic, litter and cage effects, and body weight at weaning and batch as systematic factors. The marginal posterior distributions of the model parameters were estimated using MCMC Bayesian procedures. The deviance information criterion (DIC) was used for model comparison regarding the statistical distribution of the data (normal or ZIP), and the Bayes factor was computed as a measure of the strength of evidence in favor of the host genetics, litter, and cage effects on microbial traits. According to DIC, all microbial traits were better adjusted with the linear model except for the OTU present in less than 10% of the animals, and for 25 of the 43 OTU with a frequency between 10 and 25%. On a global scale, the Bayes factor revealed substantial evidence in favor of the genetic control of the number of observed OTU and Shannon indices. At the taxon-specific level, significant proportions of the OTU and relative abundances of genera were influenced by additive genetic, litter, and cage effects. Several members of the genera Bacteroides and Parabacteroides were strongly influenced by the host genetics and nursing environment, whereas the family S24-7 and the genus Ruminococcus were strongly influenced by cage effects. CONCLUSIONS: This study demonstrates that host genetics shapes the overall rabbit cecal microbial diversity and that a significant proportion of the taxa is influenced either by host genetics or environmental factors, such as litter and/or cage.

Identification of transcriptional regulatory variants in pig duodenum, liver, and muscle tissues.

BACKGROUND: In humans and livestock species, genome-wide association studies (GWAS) have been applied to study the association between variants distributed across the genome and a phenotype of interest. To discover genetic polymorphisms affecting the duodenum, liver, and muscle transcriptomes of 300 pigs from 3 different breeds (Duroc, Landrace, and Large White), we performed expression GWAS between 25,315,878 polymorphisms and the expression of 13,891 genes in duodenum, 12,748 genes in liver, and 11,617 genes in muscle. RESULTS: More than 9.68 × 1011 association tests were performed, yielding 14,096,080 significantly associated variants, which were grouped in 26,414 expression quantitative trait locus (eQTL) regions. Over 56% of the variants were within 1 Mb of their associated gene. In addition to the 100-kb region upstream of the transcription start site, we identified the importance of the 100-kb region downstream of the 3'UTR for gene regulation, as most of the cis-regulatory variants were located within these 2 regions. We also observed 39,874 hotspot regulatory polymorphisms associated with the expression of 10 or more genes that could modify the protein structure or the expression of a regulator gene. In addition, 2 motifs (5'-GATCCNGYGTTGCYG-3' and a poly(A) sequence) were enriched across the 3 tissues within the neighboring sequences of the most significant single-nucleotide polymorphisms in each cis-eQTL region. CONCLUSIONS: The 14 million significant associations obtained in this study are publicly available and have enabled the identification of expression-associated cis-, trans-, and hotspot regulatory variants within and across tissues, thus shedding light on the molecular mechanisms of regulatory variations that shape end-trait phenotypes.

A Co-Association Network Analysis Reveals Putative Regulators for Health-Related Traits in Pigs.

In recent years, the increase in awareness of antimicrobial resistance together with the societal demand of healthier meat products have driven attention to health-related traits in livestock production. Previous studies have reported medium to high heritabilities for these traits and described genomic regions associated with them. Despite its genetic component, health- and immunity-related traits are complex and its study by association analysis with genomic markers may be missing some information. To analyse multiple phenotypes and gene-by-gene interactions, systems biology approaches, such as the association weight matrix (AWM), allows combining genome wide association study results with network inference algorithms. The present study aimed to identify gene networks, key regulators and candidate genes associated to immunocompetence in pigs by integrating multiple health-related traits, enriched for innate immune phenotypes, using the AWM approach. The co-association network analysis unveiled a network comprised of 3,636 nodes (genes) and 451,407 edges (interactions), including a total of 246 regulators. From these, five genes (ARNT2, BRMS1L, MED12L, SUPT3H and TRIM25) were selected as key regulators as they were associated with the maximum number of genes with the minimum overlapping (1,827 genes in total). The five regulators were involved in pathways related to immunity such as lymphocyte differentiation and activation, platelet activation and degranulation, megakaryocyte differentiation, FcγR-mediated phagocytosis and response to nitric oxide, among others, but also in immunometabolism. Furthermore, we identified genes co-associated with the key regulators previously reported as candidate genes (e.g., ANGPT1, CD4, CD36, DOCK1, PDE4B, PRKCE, PTPRC and SH2B3) for immunity traits in humans and pigs, but also new candidate ones (e.g., ACSL3, CXADR, HBB, MMP12, PTPN6, WLS) that were not previously described. The co-association analysis revealed new regulators associated with health-related traits in pigs. This approach also identified gene-by-gene interactions and candidate genes involved in pathways related to cell fate and metabolic and immune functions. Our results shed new light in the regulatory mechanisms involved in pig immunity and reinforce the use of the pig as biomedical model.

The value of gut microbiota to predict feed efficiency and growth of rabbits under different feeding regimes.

Gut microbiota plays an important role in nutrient absorption and could impact rabbit feed efficiency. This study aims at investigating such impact by evaluating the value added by microbial information for predicting individual growth and cage phenotypes related to feed efficiency. The dataset comprised individual average daily gain and cage-average daily feed intake from 425 meat rabbits, in which cecal microbiota was assessed, and their cage mates. Despite microbiota was not measured in all animals, consideration of pedigree relationships with mixed models allowed the study of cage-average traits. The inclusion of microbial information into certain mixed models increased their predictive ability up to 20% and 46% for cage-average feed efficiency and individual growth traits, respectively. These gains were associated with large microbiability estimates and with reductions in the heritability estimates. However, large microbiabililty estimates were also obtained with certain models but without any improvement in their predictive ability. A large proportion of OTUs seems to be responsible for the prediction improvement in growth and feed efficiency traits, although specific OTUs taxonomically assigned to 5 different phyla have a higher weight. Rabbit growth and feed efficiency are influenced by host cecal microbiota, thus considering microbial information in models improves the prediction of these complex phenotypes.