High-Fat Diet-Induced Obesity Increases Brain Mitochondrial Complex I and Lipoxidation-Derived Protein Damage.

Obesity is a risk factor for highly prevalent age-related neurodegenerative diseases, the pathogenesis of whichinvolves mitochondrial dysfunction and protein oxidative damage. Lipoxidation, driven by high levels of peroxidizable unsaturated fatty acids and low antioxidant protection of the brain, stands out as a significant risk factor. To gain information on the relationship between obesity and brain molecular damage, in a porcine model of obesity we evaluated (1) the level of mitochondrial respiratory chain complexes, as the main source of free radical generation, by Western blot; (2) the fatty acid profile by gas chromatography; and (3) the oxidative modification of proteins by mass spectrometry. The results demonstrate a selectively higher amount of the lipoxidation-derived biomarker malondialdehyde-lysine (MDAL) (34% increase) in the frontal cortex, and positive correlations between MDAL and LDL levels and body weight. No changes were observed in brain fatty acid profile by the high-fat diet, and the increased lipid peroxidative modification was associated with increased levels of mitochondrial complex I (NDUFS3 and NDUFA9 subunits) and complex II (flavoprotein). Interestingly, introducing n3 fatty acids and a probiotic in the high-fat diet prevented the observed changes, suggesting that dietary components can modulate protein oxidative modification at the cerebral level and opening new possibilities in neurodegenerative diseases' prevention.

A compendium of genetic regulatory effects across pig tissues.

The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.

Preweaning piglet survival on commercial farms.

Preweaning piglet mortality (PWM), a trait highly related to litter size, is one of the main concerns associated with productive efficiency and animal welfare in commercial pig farms. The objectives of this work were to study piglet survival at the farm level, to establish a survival rate (SR) as a target indicator to be improved, and to model it based on other reproductive parameters. Analyzed data corresponded to 580 Spanish commercial farms with a total inventory of 809,768 sows. These farms showed a mean SR of 85.70% piglets born alive (BA), which decreased to 81.81% when total piglets born (TB) were considered. The SR was strongly associated with prolificacy (P < 0.01), the parities with the highest prolificacy being those that had the lowest SR. Thus, the highest correlations were for the SR of piglets BA in the third and fourth parities (r = -0.460 and r = -0.452, respectively, P < 0.01), and for the SR of piglets TB in the fourth parity (r = -0.546, P < 0.01), which was the one with the highest prolificacy. The values corresponding to the quartile of farms with the highest SR within the most productive farms were established as targets to be improved, which were ≥88.5% of piglets BA and 83.2% of piglets TB. Nevertheless, the direct associations shown between the piglet's survival and prolificacy and other productive factors, such as the age of piglets at weaning, the farrowings per sow and year and the farrowing interval, suggest the convenience of modeling the risk of PWM on farms to have its own target of survival index to be improved.

Sow prolificacy and preweaning piglet mortality have increased parallelly on commercial farms. This loss of piglets is a concern of efficiency and animal welfare, and it requires the improvement of piglet survival by reducing the number of stillborn piglets and piglet mortality during lactation, paying particular attention to hyperprolific sows (≥15 total piglets born per litter). Data from 580 commercial farms with an average inventory of 809,768 sows have been analyzed to propose two predictive models based on several reproductive parameters and two survival rate targets with the aim of reducing this problem, which are ≥88.5% of piglets born alive and ≥83.2% of total piglets born.

A meta-analysis of genetic and phenotypic diversity of European local pig breeds reveals genomic regions associated with breed differentiation for production traits.

BACKGROUND: Intense selection of modern pig breeds has resulted in genetic improvement of production traits while the performance of local pig breeds has remained lower. As local pig breeds have been bred in extensive systems, they have adapted to specific environmental conditions, resulting in a rich genotypic and phenotypic diversity. This study is based on European local pig breeds that have been genetically characterized using DNA-pool sequencing data and phenotypically characterized using breed level phenotypes related to stature, fatness, growth, and reproductive performance traits. These data were analyzed using a dedicated approach to detect signatures of selection linked to phenotypic traits in order to uncover potential candidate genes that may underlie adaptation to specific environments. RESULTS: Analysis of the genetic data of European pig breeds revealed four main axes of genetic variation represented by the Iberian and three modern breeds (i.e. Large White, Landrace, and Duroc). In addition, breeds clustered according to their geographical origin, for example French Gascon and Basque breeds, Italian Apulo Calabrese and Casertana breeds, Spanish Iberian, and Portuguese Alentejano breeds. Principal component analysis of the phenotypic data distinguished the larger and leaner breeds with better growth potential and reproductive performance from the smaller and fatter breeds with low growth and reproductive efficiency. Linking the signatures of selection with phenotype identified 16 significant genomic regions associated with stature, 24 with fatness, 2 with growth, and 192 with reproduction. Among them, several regions contained candidate genes with possible biological effects on stature, fatness, growth, and reproductive performance traits. For example, strong associations were found for stature in two regions containing, respectively, the ANXA4 and ANTXR1 genes, for fatness in a region containing the DNMT3A and POMC genes and for reproductive performance in a region containing the HSD17B7 gene. CONCLUSIONS: In this study on European local pig breeds, we used a dedicated approach for detecting signatures of selection that were supported by phenotypic data at the breed level to identify potential candidate genes that may have adapted to different living environments and production systems.

Multi-Criteria Evaluation Model of Management for Weaned Piglets and Its Relations with Farm Performance and Veterinary Medicine Consumption.

Weaned piglets, being immature, demand careful handling to mitigate post-weaning stress in order to avoid immunosuppression and the use of antimicrobials to palliate the effects of disease outbreaks due to poor management. The objective of this work is to design a quick scan calculator or multi-criteria evaluation model of management for weaned piglets, founded on 10 critical indices covering post-weaning management aspects based on hygienic measures and management of facilities and animals. These include pre-weaning handling, batch management, biosecurity, water management, feed management, health program, stockmen training, temperature, ventilation, and floor conditions and density to relate handling and hygiene practices with farm performance and the consumption of veterinary medication. Each index carries a maximum score of ten, with evaluations derived from different management factors that make up each index (from three to eight factors were evaluated per index). Their cumulative score reflects the degree of adequacy of on-farm management. Therefore, a perfectly managed farm would achieve 100 points. The calculator underwent testing on 23 intensive farms with a total population of close to 16,000 sows and more than 400,000 weaned piglets, revealing the highest mean scores in floor conditions and density (8.03 out of 10) and pre-weaning handling and health programs (6.87 and 6.28, respectively). Conversely, the lowest scores corresponded to temperature, ventilation, water management, and stockmen training (4.08, 4.32, 4.81, and 4.93, respectively). The assessed farms averaged a global score of 56.12 out of 100 (from 37.65 to 76.76). The calculator's global score correlated with key post-weaning productivity and piglet health indicators, such as the feed conversion ratio, mortality rate, and piglet production cost, with r values of -0.442, -0.437, and -0.435, respectively (p < 0.05). Additionally, it negatively correlated with medication costs per piglet (r = -0.414; p < 0.05) and positively with annual farm productivity (r = 0.592; p < 0.01). To enhance management, hygiene, and prevention, farms should prioritize addressing indices with the lowest scores, thereby reducing medication consumption and enhancing productivity and health outcomes. Additionally, this quick scan calculator can be used for benchmarking purposes.

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.

Assessing the relationship between the in silico predicted consequences of 97 missense mutations mapping to 68 genes related to lipid metabolism and their association with porcine fatness traits.

In general, the relationship between the predicted functional consequences of missense mutations mapping to genes known to be involved in human diseases and the severity of disease manifestations is weak. In this study, we tested in pigs whether missense single nucleotide polymorphisms (SNPs), predicted to have consequences on the function of genes related to lipid metabolism are associated with lipid phenotypes. Association analysis demonstrated that nine out of 72 nominally associated SNPs were classified as "highly" or "very highly consistent" in silico-predicted functional mutations and did not show association with lipid traits expected to be affected by inactivation of the corresponding gene. Although the lack of endophenotypes and the limited sample size of certain genotypic classes might have limited to some extent the reach of the current study, our data indicate that present-day bioinformatic tools have a modest ability to predict the impact of missense mutations on complex phenotypes.

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.

Simultaneous placement of fecal and urinary continent channel stomas in the umbilicus: Single-center experience.

INTRODUCTION: In patients with urinary continent channel (UCC) and Malone Antegrade Continent Enema (MACE) procedures, two separate abdominal stomas are needed. The umbilicus is a preferred site for single channel stomas given the ability to conceal the stoma. However, there are no studies describing outcomes of both stomas being created in the umbilicus. We aimed to describe our experience in patients who underwent UCC and MACE stomas both placed in the umbilicus. METHODS: A retrospective review from 2009 to 2020 was performed in our institution for patients who underwent the creation of UCC and MACE stomas simultaneously in the umbilicus. The variation in the technique involves two V-skin shaped flaps in the umbilicus; the MACE and UCC stomas are delivered from both flaps and placed at the right and left side respectively. Patients with greater than 3 months of follow-up were included in the study. RESULTS: There were 17 patients identified with the median age of 13.5 years and a median follow-up of 32.8 months. The mean BMI percentile was 89.5%. Monti technique and split appendix with cecal extension were utilized in 8 (47.1%) and 7 (41.2%) patients respectively and 13 (76.5%) patients required concurrent urological procedures. All channel-related complications occurred within a mean time of 15.7 months. Skin-level stenosis in the MACE occurred in 5 (29.4%) events, and all were successfully managed by placing an indwelling catheter for up to two weeks. There were 2 (11.8%) complications related to UCC, which required subfascial minor surgical revision. The rate of patients with symptomatic UTI decreased 35.3% postoperatively, and no new onset of UTI occurred in patients without a prior history of UTI. During follow-up, all patients remained dry between CIC, however one had occasional leakage related to delay in catheterization. Total fecal continence was achieved in 14 (82.3%) patients. Additionally, 3 (16.6%) patients experienced improvement in fecal continence with sporadic soiling episodes. COMMENTS: Placement of UCC and MACE stomas in the umbilicus demonstrate a percentage of complication of 7/34 (20.6%) with only 2 patients requiring surgical intervention, comparable to the standard. UTI rate decreased in patients with a prior history of UTI. We believe the patients' perspective and degree of satisfaction will fully determine the benefits of this technique. CONCLUSIONS: Simultaneous UCC and MACE stomas placed at the umbilicus showed good functional outcomes and similar complication rates to traditional approach where stomas were placed separately in the abdominal wall.

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.

Admixture and breed traceability in European indigenous pig breeds and wild boar using genome-wide SNP data.

Preserving diversity of indigenous pig (Sus scrofa) breeds is a key factor to (i) sustain the pork chain (both at local and global scales) including the production of high-quality branded products, (ii) enrich the animal biobanking and (iii) progress conservation policies. Single nucleotide polymorphism (SNP) chips offer the opportunity for whole-genome comparisons among individuals and breeds. Animals from twenty European local pigs breeds, reared in nine countries (Croatia: Black Slavonian, Turopolje; France: Basque, Gascon; Germany: Schwabisch-Hällisches Schwein; Italy: Apulo Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano, Sarda; Lithuania: Indigenous Wattle, White Old Type; Portugal: Alentejana, Bísara; Serbia: Moravka, Swallow-Bellied Mangalitsa; Slovenia: Krskopolje pig; Spain: Iberian, Majorcan Black), and three commercial breeds (Duroc, Landrace and Large White) were sampled and genotyped with the GeneSeek Genomic Profiler (GGP) 70 K HD porcine genotyping chip. A dataset of 51 Wild Boars from nine countries was also added, summing up to 1186 pigs (~ 49 pigs/breed). The aim was to: (i) investigate individual admixture ancestries and (ii) assess breed traceability via discriminant analysis on principal components (DAPC). Albeit the mosaic of shared ancestries found for Nero Siciliano, Sarda and Moravka, admixture analysis indicated independent evolvement for the rest of the breeds. High prediction accuracy of DAPC mark SNP data as a reliable solution for the traceability of breed-specific pig products.

Glaucoma Cascade Screening in a High Risk Afro-Caribbean Haitian Population: A Pilot Study.

PRCIS: Glaucoma cascade screening in first-degree relatives (FDRs) of young Haitian glaucoma patients had high yield for diagnosing manifest and suspected glaucoma in 30.8% of those screened despite modest participation. PURPOSE: To evaluate the outcomes of glaucoma cascade screening in FDRs (parents, siblings, and offspring) of Haitian juvenile open-angle glaucoma (JOAG) patients. PATIENTS AND METHODS: Consecutive index patients (Haitians with JOAG) were identified, and the number/type of FDRs residing in South Florida were recorded. These FDRs were invited for free glaucoma screening, which included a comprehensive ophthalmic exam, gonioscopy, automated visual field testing and optical coherence tomographic analysis of the retinal nerve fiber layers. FDR characteristics and clinical findings from screening are reported. RESULTS: A total of 77 FDRs were invited, 26 (33.8%) agreed to undergo screening (18 females, 9 males), which revealed 2 (7.7%) with manifest glaucoma (mean age 77.5 y; one of whom was previously unaware of his glaucoma diagnosis), 6 (23.1%) with suspected glaucoma (mean age 29.8±18.3 y), and 18 (69.2%) without manifest or suspected glaucoma (mean age 37.2±21.8 y). Siblings of index patients were least likely to participate in cascade glaucoma screening when compared with index patients' parents or offspring. FDR eyes with manifest glaucoma had significantly worse best-corrected visual acuities, higher intraocular pressures, thinner central corneal thicknesses, and thinner circumferential papillary retinal nerve fiber layer thicknesses than those without glaucoma. CONCLUSION: Glaucoma cascade screening of Haitian JOAG patients' FDRs revealed that 30.8% had suspected or manifest glaucoma. Future efforts centered on provider-initiated recruitment and improving public glaucoma awareness and education may increase screening participation.

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.

Leveraging host-genetics and gut microbiota to determine immunocompetence in pigs.

BACKGROUND: The gut microbiota influences host performance playing a relevant role in homeostasis and function of the immune system. The aim of the present work was to identify microbial signatures linked to immunity traits and to characterize the contribution of host-genome and gut microbiota to the immunocompetence in healthy pigs. RESULTS: To achieve this goal, we undertook a combination of network, mixed model and microbial-wide association studies (MWAS) for 21 immunity traits and the relative abundance of gut bacterial communities in 389 pigs genotyped for 70K SNPs. The heritability (h2; proportion of phenotypic variance explained by the host genetics) and microbiability (m2; proportion of variance explained by the microbial composition) showed similar values for most of the analyzed immunity traits, except for both IgM and IgG in plasma that was dominated by the host genetics, and the haptoglobin in serum which was the trait with larger m2 (0.275) compared to h2 (0.138). Results from the MWAS suggested a polymicrobial nature of the immunocompetence in pigs and revealed associations between pigs gut microbiota composition and 15 of the analyzed traits. The lymphocytes phagocytic capacity (quantified as mean fluorescence) and the total number of monocytes in blood were the traits associated with the largest number of taxa (6 taxa). Among the associations identified by MWAS, 30% were confirmed by an information theory network approach. The strongest confirmed associations were between Fibrobacter and phagocytic capacity of lymphocytes (r = 0.37), followed by correlations between Streptococcus and the percentage of phagocytic lymphocytes (r = -0.34) and between Megasphaera and serum concentration of haptoglobin (r = 0.26). In the interaction network, Streptococcus and percentage of phagocytic lymphocytes were the keystone bacterial and immune-trait, respectively. CONCLUSIONS: Overall, our findings reveal an important connection between gut microbiota composition and immunity traits in pigs, and highlight the need to consider both sources of information, host genome and microbial levels, to accurately characterize immunocompetence in pigs.

Variability in porcine microRNA genes and its association with mRNA expression and lipid phenotypes.

BACKGROUND: Mature microRNAs (miRNAs) play an important role in repressing the expression of a wide range of mRNAs. The presence of polymorphic sites in miRNA genes and their corresponding 3'UTR binding sites can disrupt canonical conserved miRNA-mRNA pairings, and thus modify gene expression patterns. However, to date such polymorphic sites in miRNA genes and their association with gene expression phenotypes and complex traits are poorly characterized in pigs. RESULTS: By analyzing whole-genome sequences from 120 pigs and wild boars from Europe and Asia, we identified 285 single nucleotide polymorphisms (SNPs) that map to miRNA loci, and 109,724 SNPs that are located in predicted 7mer-m8 miRNA binding sites within porcine 3'UTR. In porcine miRNA genes, SNP density is reduced compared with their flanking non-miRNA regions. By sequencing the genomes of five Duroc boars, we identified 12 miRNA SNPs that were subsequently genotyped in their offspring (N = 345, Lipgen population). Association analyses of miRNA SNPs with 38 lipid-related traits and hepatic and muscle microarray expression phenotypes recorded in the Lipgen population were performed. The most relevant detected association was between the genotype of the rs319154814 (G/A) SNP located in the apical loop of the ssc-miR-326 hairpin precursor and PPP1CC mRNA levels in the liver (q-value = 0.058). This result was subsequently confirmed by qPCR (P-value = 0.027). The rs319154814 (G/A) genotype was also associated with several fatty acid composition traits. CONCLUSIONS: Our findings show a reduced variability of porcine miRNA genes, which is consistent with strong purifying selection, particularly in the seed region that plays a critical role in miRNA binding. Although it is generally assumed that SNPs mapping to the seed region are those with the most pronounced consequences on mRNA expression, we show that a SNP mapping to the apical region of ssc-miR-326 is significantly associated with hepatic mRNA levels of the PPP1CC gene, one of its predicted targets. Although experimental confirmation of such an interaction is reported in humans but not in pigs, this result highlights the need to further investigate the functional effects of miRNA polymorphisms that are located outside the seed region on gene expression in pigs.

kernInt: A Kernel Framework for Integrating Supervised and Unsupervised Analyses in Spatio-Temporal Metagenomic Datasets.

The advent of next-generation sequencing technologies allowed relative quantification of microbiome communities and their spatial and temporal variation. In recent years, supervised learning (i.e., prediction of a phenotype of interest) from taxonomic abundances has become increasingly common in the microbiome field. However, a gap exists between supervised and classical unsupervised analyses, based on computing ecological dissimilarities for visualization or clustering. Despite this, both approaches face common challenges, like the compositional nature of next-generation sequencing data or the integration of the spatial and temporal dimensions. Here we propose a kernel framework to place on a common ground the unsupervised and supervised microbiome analyses, including the retrieval of microbial signatures (taxa importances). We define two compositional kernels (Aitchison-RBF and compositional linear) and discuss how to transform non-compositional beta-dissimilarity measures into kernels. Spatial data is integrated with multiple kernel learning, while longitudinal data is evaluated by specific kernels. We illustrate our framework through a single point soil dataset, a human dataset with a spatial component, and a previously unpublished longitudinal dataset concerning pig production. The proposed framework and the case studies are freely available in the kernInt package at https://github.com/elies-ramon/kernInt.

A gene co-association network regulating gut microbial communities in a Duroc pig population.

BACKGROUND: Analyses of gut microbiome composition in livestock species have shown its potential to contribute to the regulation of complex phenotypes. However, little is known about the host genetic control over the gut microbial communities. In pigs, previous studies are based on classical "single-gene-single-trait" approaches and have evaluated the role of host genome controlling gut prokaryote and eukaryote communities separately. RESULTS: In order to determine the ability of the host genome to control the diversity and composition of microbial communities in healthy pigs, we undertook genome-wide association studies (GWAS) for 39 microbial phenotypes that included 2 diversity indexes, and the relative abundance of 31 bacterial and six commensal protist genera in 390 pigs genotyped for 70 K SNPs. The GWAS results were processed through a 3-step analytical pipeline comprised of (1) association weight matrix; (2) regulatory impact factor; and (3) partial correlation and information theory. The inferred gene regulatory network comprised 3561 genes (within a 5 kb distance from a relevant SNP-P < 0.05) and 738,913 connections (SNP-to-SNP co-associations). Our findings highlight the complexity and polygenic nature of the pig gut microbial ecosystem. Prominent within the network were 5 regulators, PRDM15, STAT1, ssc-mir-371, SOX9 and RUNX2 which gathered 942, 607, 588, 284 and 273 connections, respectively. PRDM15 modulates the transcription of upstream regulators of WNT and MAPK-ERK signaling to safeguard naive pluripotency and regulates the production of Th1- and Th2-type immune response. The signal transducer STAT1 has long been associated with immune processes and was recently identified as a potential regulator of vaccine response to porcine reproductive and respiratory syndrome. The list of regulators was enriched for immune-related pathways, and the list of predicted targets includes candidate genes previously reported as associated with microbiota profile in pigs, mice and human, such as SLIT3, SLC39A8, NOS1, IL1R2, DAB1, TOX3, SPP1, THSD7B, ELF2, PIANP, A2ML1, and IFNAR1. Moreover, we show the existence of host-genetic variants jointly associated with the relative abundance of butyrate producer bacteria and host performance. CONCLUSIONS: Taken together, our results identified regulators, candidate genes, and mechanisms linked with microbiome modulation by the host. They further highlight the value of the proposed analytical pipeline to exploit pleiotropy and the crosstalk between bacteria and protists as significant contributors to host-microbiome interactions and identify genetic markers and candidate genes that can be incorporated in breeding program to improve host-performance and microbial traits. Video Abstract.