Search

Abstract Background The composition and distribution of fatty acids (FA) are important factors determining the quality, flavor, and nutrient value of meat. In addition, FAs synthesized in the body participate in energy metabolism and are involved in different regulatory pathways in the form of signaling molecules or by acting as agonist or antagonist ligands of different nuclear receptors. Finally, synthesis and catabolism of FAs affect adaptive immunity by regulating lymphocyte metabolism. The present study performed genome-wide association studies using FA profiles of blood, liver, backfat and muscle from 432 commercial Duroc pigs. Results Twenty-five genomic regions located on 15 Sus scrofa chromosomes (SSC) were detected. Annotation of the quantitative trait locus (QTL) regions identified 49 lipid metabolism-related candidate genes. Among these QTLs, four were identified in more than one tissue. The ratio of C20:4n-6/C20:3n-6 was associated with the region on SSC2 at 7.56–14.26 Mb for backfat, liver, and muscle. Members of the fatty acid desaturase gene cluster (FADS1, FADS2, and FADS3) are the most promising candidate genes in this region. Two QTL regions on SSC14 (103.81–115.64 Mb and 100.91–128.14 Mb) were identified for FA desaturation in backfat and muscle. In addition, two separate regions on SSC9 at 0 – 14.55 Mb and on SSC12 at 0–1.91 Mb were both associated with the same multiple FA traits for backfat, with candidate genes involved in de novo FA synthesis and triacylglycerol (TAG) metabolism, such as DGAT2 and FASN. The ratio C20:0/C18:0 was associated with the region on SSC5 at 64.84–78.32 Mb for backfat. Furthermore, the association of the C16:0 content with the region at 118.92–123.95 Mb on SSC4 was blood specific. Finally, candidate genes involved in de novo lipogenesis regulate T cell differentiation and promote the generation of palmitoleate, an adipokine that alleviates inflammation. Conclusions Several SNPs and candidate genes were associated with lipid metabolism in blood, liver, backfat, and muscle. These results contribute to elucidating the molecular mechanisms implicated in the determination of the FA profile in different pig tissues and can be useful in selection programs that aim to improve health and energy metabolism in pigs.

Stress significantly affects the health, welfare, and productivity of farm animals. We performed a longitudinal study to evaluate stress’s effects on pig performance, feeding behavior, and fecal microbiota composition. This study involved 64 Duroc pigs during the fattening period, divided into two experimental groups: a stress group (n = 32) and a control group (n = 32). Stressed groups had less space and were mixed twice during the experiment. We monitored body weight, feed efficiency, feeding behavior, and fecal microbiota composition. Compared to the control group, the stressed pigs exhibited reduced body weight, feed efficiency, fewer feeder visits, and longer meal durations. In the fecal microbiota, resilience was observed, with greater differences between groups when sampling was closer to the stressful stimulus. Stressed pigs showed an increase in opportunistic bacteria, such as Streptococcus, Treponema and members of the Erysipelotrichaceae family, while control pigs had more butyrate- and propionate-producing genera like Anaerobutyricum, Coprococcus and HUN007. Our findings confirm that prolonged stress negatively impacts porcine welfare, behavior, and performance, and alters their gut microbiota. Specific microorganisms identified could serve as non-invasive biomarkers for stress, potentially informing both animal welfare and similar gut-brain axis mechanisms relevant to human research.

Organic livestock farming is committed to high environmental and animal welfare standards, although pathologies such as post-weaning diarrhoea (PWD) may appear. The main objective of this study was to assess nutritional strategies to prevent PWD in organic piglets. A total of 134 weaned piglets were fed one of three diets: high crude protein (17.8%, HCP), low crude protein (16.8%, LCP), and low crude protein supplemented with liquid whey (LCP+W). Piglets were assessed weekly for four weeks on the following parameters: diarrhoea incidence, additional health parameters, average daily gain, and behaviour. Faecal samples were taken to analyse the intestinal microbiota composition. Data were analysed using LMM and GLMM models and Shannon and Whittaker indexes. No significant effect of diet on diarrhoea incidence was found, but the LCP+W diet increased average daily gain. Pigs fed the LCP+W diet presented a lower percentage of drinking and negative social behaviour compared with the HCP diet, and LCP pigs presented higher exploration compared with HCP. In addition, LCP+W piglets showed a higher abundance of the beneficial genus Frisingicoccus. Although liquid whey did not reduce diarrhoea incidence, the benefits found in growth, microbiota composition, and reduced negative social behaviour indicate that it could be an optimal supplement to organic diets.

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.

ABSTRACT 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.

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.

Abstract Background The effect of the production environment and different management practices in rabbit cecal microbiota remains poorly understood. While previous studies have proved the impact of the age or the feed composition, research in the breeding farm and other animal management aspects, such as the presence of antibiotics in the feed or the level of feeding, is still needed. Characterization of microbial diversity and composition of growing rabbits raised under different conditions could help better understand the role these practices play in cecal microbial communities and how it may result in different animal performance. Results Four hundred twenty-five meat rabbits raised in two different facilities, fed under two feeding regimes (ad libitum or restricted) with feed supplemented or free of antibiotics, were selected for this study. A 16S rRNA gene-based assessment through the MiSeq Illumina sequencing platform was performed on cecal samples collected from these individuals at slaughter. Different univariate and multivariate approaches were conducted to unravel the influence of the different factors on microbial alpha diversity and composition at phylum, genus and OTU taxonomic levels. The animals raised in the facility harboring the most stable environmental conditions had greater, and less variable, microbial richness and diversity. Bootstrap univariate analyses of variance and sparse partial least squares-discriminant analyses endorsed that farm conditions exerted an important influence on rabbit microbiota since the relative abundances of many taxa were found differentially represented between both facilities at all taxonomic levels characterized. Furthermore, only five OTUs were needed to achieve a perfect classification of samples according to the facility where animals were raised. The level of feeding and the presence of antibiotics did not modify the global alpha diversity but had an impact on some bacteria relative abundances, albeit in a small number of taxa compared with farm, which is consistent with the lower sample classification power according to these factors achieved using microbial information. Conclusions This study reveals that factors associated with the farm effect and other management factors, such as the presence of antibiotics in the diet or the feeding level, modify cecal microbial communities. It highlights the importance of offering a controlled breeding environment that reduces differences in microbial cecal composition that could be responsible for different animal performance.

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.

Abstract Background To date, the molecular mechanisms that underlie residual feed intake (RFI) in pigs are unknown. Results from different genome-wide association studies and gene expression analyses are not always consistent. The aim of this research was to use machine learning to identify genes associated with feed efficiency (FE) using transcriptomic (RNA-Seq) data from pigs that are phenotypically extreme for RFI. Methods RFI was computed by considering within-sex regression on mean metabolic body weight, average daily gain, and average backfat gain. RNA-Seq analyses were performed on liver and duodenum tissue from 32 high and 33 low RFI pigs collected at 153 d of age. Machine-learning algorithms were used to predict RFI class based on gene expression levels in liver and duodenum after adjusting for batch effects. Genes were ranked according to their contribution to the classification using the permutation accuracy importance score in an unbiased random forest (RF) algorithm based on conditional inference. Support vector machine, RF, elastic net (ENET) and nearest shrunken centroid algorithms were tested using different subsets of the top rank genes. Nested resampling for hyperparameter tuning was implemented with tenfold cross-validation in the outer and inner loops. Results The best classification was obtained with ENET using the expression of 200 genes in liver [area under the receiver operating characteristic curve (AUROC): 0.85; accuracy: 0.78] and 100 genes in duodenum (AUROC: 0.76; accuracy: 0.69). Canonical pathways and candidate genes that were previously reported as associated with FE in several species were identified. The most remarkable pathways and genes identified were NRF2-mediated oxidative stress response and aldosterone signalling in epithelial cells, the DNAJC6, DNAJC1, MAPK8, PRKD3 genes in duodenum, and melatonin degradation II, PPARα/RXRα activation, and GPCR-mediated nutrient sensing in enteroendocrine cells and SMOX, IL4I1, PRKAR2B, CLOCK and CCK genes in liver. Conclusions ML algorithms and RNA-Seq expression data were found to provide good performance for classifying pigs into high or low RFI groups. Classification was better with gene expression data from liver than from duodenum. Genes associated with FE in liver and duodenum tissue that can be used as predictive biomarkers for this trait were identified.

Abstract Background The molecular basis of compensatory growth in monogastric animals has not yet been fully explored. Herewith, in this study we aim to determine changes in the pig skeletal muscle transcriptome profile during compensatory growth following a feed restriction period. A RNA-Seq experiment was performed with a total of 24 females belonging to a Duroc commercial line. Half of the animals received either a restricted (RE) or ad libitum (AL) diet during the first fattening period (60–125 d of age). After that, all gilts were fed ad libitum for a further ~30 d until the age of ~155 d, when animals were slaughtered and samples of gluteus medius muscle were harvested to perform RNA-Seq analyses and intramuscular fat content determination. Results During the period following food restriction, RE animals re-fed ad libitum displayed compensatory growth, showed better feed conversion rate and tended to deposit more subcutaneous fat than AL fed animals. Animals were slaughtered in the phase of accelerated growth, when RE animals had not completely compensated the performance of AL group, showing lower live and carcass weights. At intramuscular level, RE gilts showed a higher content of polyunsaturated fatty acids during the compensatory growth phase. The comparison of RE and AL expression profiles allowed the identification of 86 (ǀlog2Fold-Changeǀ > 1, padj

Abstract Background The genetic basis of muscle fat deposition in pigs is not well known. So far, we have only identified a limited number of genes involved in the absorption, transport, storage and catabolism of lipids. Such information is crucial to interpret, from a biological perspective, the results of genome-wide association analyses for intramuscular fat content and composition traits. Herewith, we have investigated how the ingestion of food changes gene expression in the gluteus medius muscle of Duroc pigs. Results By comparing the muscle mRNA expression of fasted pigs (T0) with that of pigs sampled 5 h (T1) and 7 h (T2) after food intake, we have detected differential expression (DE) for 148 (T0-T1), 520 (T0-T2) and 135 (T1-T2) genes (q-value of 1.5). Many of these DE genes were transcription factors, suggesting that we have detected the coordinated response of the skeletal muscle to nutrient supply. We also found DE genes with a dual role in oxidative stress and angiogenesis (THBS1, THBS2 and TXNIP), two biological processes that are probably activated in the post-prandial state. Finally, we have identified several loci playing a key role in the modulation of circadian rhythms (ARNTL, PER1, PER2, BHLHE40, NR1D1, SIK1, CIART and CRY2), a result that indicates that the porcine muscle circadian clock is modulated by nutrition. Conclusion We have shown that hundreds of genes change their expression in the porcine skeletal muscle in response to nutrient intake. Many of these loci do not have a known metabolic role, a result that suggests that our knowledge about the genetic basis of muscle energy homeostasis is still incomplete.

To gain insight into the importance of carefully selecting the sampling area for intestinal microbiota studies, cecal and fecal microbial communities of Caldes meat rabbit were characterized. The animals involved in the study were divided in two groups according to the feed intake level they received during the fattening period; ad libitum (n = 10) or restricted to 75% of ad libitum intake (n = 11). Cecum and internal hard feces were sampled from sacrificed animals. Assessment of bacterial and archaeal populations was performed by means of Illumina sequencing of 16S rRNA gene amplicons in a MiSeq platform. A total of 596 operational taxonomic units (OTUs) were detected using QIIME software. Taxonomic assignment revealed that microbial diversity was dominated by phyla Firmicutes (76.42%), Tenericutes (7.83%), and Bacteroidetes (7.42%); kingdom Archaea was presented at low percentage (0.61%). No significant differences were detected between sampling origins in microbial diversity or richness assessed using two alpha-diversity indexes: Shannon and the observed number of OTUs. However, the analysis of variance at genus level revealed a higher presence of genera Clostridium, Anaerofustis, Blautia, Akkermansia, rc4-4, and Bacteroides in cecal samples. By contrast, genera Oscillospira and Coprococcus were found to be overrepresented in feces, suggesting that bacterial species of these genera would act as fermenters at the end of feed digestion process. At the lowest taxonomic level, 83 and 97 OTUs in feces and cecum, respectively, were differentially represented. Multivariate statistical assessment revealed that sparse partial least squares discriminant analysis (sPLS-DA) was the best approach for this purpose. Interestingly, the majority of the most discriminative OTUs selected by sPLS-DA were found to be differentially represented between sampling origins in univariate analysis. Our study provides evidence that the choice of intestinal sampling area is relevant due to important differences in some taxa’s relative abundance that have been revealed between rabbits’ cecal and fecal microbiota. An appropriate sampling intestinal area should be chosen in each microbiota assessment.

Abstract Background The pig gut microbiome harbors thousands of species of archaea, bacteria, viruses and eukaryotes such as protists and fungi. However, since the majority of published studies have been focused on prokaryotes, little is known about the diversity, host-genetic control, and contributions to host performance of the gut eukaryotic counterparts. Here we report the first study that aims at characterizing the diversity and composition of gut commensal eukaryotes in pigs, exploring their putative control by host genetics, and analyzing their association with piglets body weight. Results Fungi and protists from the faeces of 514 healthy Duroc pigs of two sexes and two different ages were characterized by 18S and ITS ribosomal RNA gene sequencing. The pig gut mycobiota was dominated by yeasts, with a high prevalence and abundance of Kazachstania spp. Regarding protists, representatives of four genera (Blastocystis, Neobalantidium, Tetratrichomonas and Trichomitus) were predominant in more than the 80% of the pigs. Heritabilities for the diversity and abundance of gut eukaryotic communities were estimated with the subset of 60d aged piglets (N = 390). The heritabilities of α-diversity and of the abundance of fungal and protists genera were low, ranging from 0.15 to 0.28. A genome wide association study reported genetic variants related to the fungal α-diversity and to the abundance of Blastocystis spp. Annotated candidate genes were mainly associated with immunity, gut homeostasis and metabolic processes. Additionally, we explored the association of gut commensal eukaryotes with piglet body weight. Our results pointed to a positive contribution of fungi from the Kazachstania genus, while protists displayed both positive (Blastocystis and Entamoeba) and negative (Trichomitus) associations with piglet body weight. Conclusions Our results point towards a minor and taxa specific genetic control over the diversity and composition of the pig gut eukaryotic communities. Moreover, we provide evidences of the associations between piglets’ body weight after weaning and members from the gut fungal and protist eukaryote community. Overall, this study highlights the relevance of considering, along with that of bacteria, the contribution of the gut eukaryote communities to better understand host-microbiome association and their role on pig performance, welfare and health.