Search

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

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.

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.

[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 (ADG) or restricted feeding (ADG)., [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 ADG and ADG, respectively. Many of these OTU had a negative effect on both traits. The largest effects on ADG 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 ADG, 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.

Weaning is a critical period in the life of pigs with repercussions on their health and welfare and on the economy of the swine industry. This study aimed to assess the efect of the commercial early weaning on gut microbiota, intestinal gene expression and serum metabolomic response via an integratedomic approach combining 16S rRNA gene sequencing, the OpenArray gene expression technology and 1 H-NMR spectroscopy. Fourteen piglets from diferent litters were sampled for blood, jejunum tissue and caecal content two days before (− 2d), and three days after (+3d) weaning. A clearly diferential ordination of caecal microbiota was observed. Higher abundances of Roseburia, Ruminococcus, Coprococcus, Dorea and Lachnospira genera in weaned piglets compared to prior to weaning showed the quick microbial changes of the piglets’ gut microbiota. Downregulation of OCLN, CLDN4, MUC2, MUC13, SLC15A1 and SLC13A1 genes, also evidenced the negative impact of weaning on gut barrier and digestive functions. Metabolomic approach pinpointed signifcant decreases in choline, LDL, triglycerides, fatty acids, alanine and isoleucine and increases in 3-hydroxybutyrate after weaning. Moreover, the correlation between microbiota and metabolome datasets revealed the existence of metabolic clusters interrelated to diferent bacterial clusters. Our results demonstrate the impact of weaning stress on the piglet and give insights regarding the associations between gut microbiota and the animal gene activity and metabolic response.

The aim of this study was to determine the possible impact of early socialization and an enriched neonatal environment to improve adaptation of piglets to weaning. We hypothesized that changes in the microbiota colonization process and in their metabolic response and intestinal functionality could help the animals face weaning stress. A total of 48 sows and their litters were allotted into a control (CTR) or an enriched treatment (ENR), in which piglets from two adjacent pens were combined and enriched with toys. The pattern of caecal microbial colonization, the jejunal gene expression, the serum metabolome and the intestinal physiology of the piglets were assessed before (-2 d) and after weaning (+ 3d). A differential ordination of caecal microbiota was observed after weaning. Serum metabolome suggested a reduced energetic metabolism in ENR animals, as evidenced by shifts in triglycerides and fatty acids, VLDL/LDL and creatine regions. The TLR2 gene showed to be downregulated in the jejunum of ENR pigs after weaning. The integration of gene expression, metabolome and microbiota datasets confirmed that differences between barren and enriched neonatal environments were evident only after weaning. Our results suggest that improvements in adaptation to weaning could be mediated by a better response to the post-weaning stress.

[Background]: The role of non-coding RNAs in the porcine muscle metabolism is poorly understood, with few studies investigating their expression patterns in response to nutrient supply. Therefore, we aimed to investigate the changes in microRNAs (miRNAs), long intergenic non-coding RNAs (lincRNAs) and mRNAs muscle expression before and after food intake., [Results]: We measured the miRNA, lincRNA and mRNA expression levels in the gluteus medius muscle of 12 gilts in a fasting condition (AL-T0) and 24 gilts fed ad libitum during either 5 h. (AL-T1, N = 12) or 7 h. (AL-T2, N = 12) prior to slaughter. The small RNA fraction was extracted from muscle samples retrieved from the 36 gilts and sequenced, whereas lincRNA and mRNA expression data were already available. In terms of mean and variance, the expression profiles of miRNAs and lincRNAs in the porcine muscle were quite different than those of mRNAs. Food intake induced the differential expression of 149 (AL-T0/AL-T1) and 435 (AL-T0/AL-T2) mRNAs, 6 (AL-T0/AL-T1) and 28 (AL-T0/AL-T2) miRNAs and none lincRNAs, while the number of differentially dispersed genes was much lower. Among the set of differentially expressed miRNAs, we identified ssc-miR-148a-3p, ssc-miR-22-3p and ssc-miR-1, which play key roles in the regulation of glucose and lipid metabolism. Besides, co-expression network analyses revealed several miRNAs that putatively interact with mRNAs playing key metabolic roles and that also showed differential expression before and after feeding. One case example was represented by seven miRNAs (ssc-miR-148a-3p, ssc-miR-151-3p, ssc-miR-30a-3p, ssc-miR-30e-3p, ssc-miR-421-5p, ssc-miR-493-5p and ssc-miR-503) which putatively interact with the PDK4 mRNA, one of the master regulators of glucose utilization and fatty acid oxidation., [Conclusions]: As a whole, our results evidence that microRNAs are likely to play an important role in the porcine skeletal muscle metabolic adaptation to nutrient availability.

[Background]: Genetic pressure in animal breeding is sparking the interest of breeders for selecting elite boars with higher sperm quality to optimize ejaculate doses and fertility rates. However, the molecular basis of sperm quality is not yet fully understood. Our aim was to identify candidate genes, pathways and DNA variants associated to sperm quality in swine by analysing 25 sperm-related phenotypes and integrating genome-wide association studies (GWAS) and RNA-seq under a systems biology framework., [Results]: By GWAS, we identified 12 quantitative trait loci (QTL) associated to the percentage of head and neck abnormalities, abnormal acrosomes and motile spermatozoa. Candidate genes included CHD2, KATNAL2, SLC14A2 and ABCA1. By RNA-seq, we identified a wide repertoire of mRNAs (e.g. PRM1, OAZ3, DNAJB8, TPPP2 and TNP1) and miRNAs (e.g. ssc-miR-30d, ssc-miR-34c, ssc-miR-30c-5p, ssc-miR-191, members of the let-7 family and ssc-miR-425-5p) with functions related to sperm biology. We detected 6128 significant correlations (P-value ≤ 0.05) between sperm traits and mRNA abundances. By expression (e)GWAS, we identified three trans-expression QTL involving the genes IQCJ, ACTR2 and HARS. Using the GWAS and RNA-seq data, we built a gene interaction network. We considered that the genes and interactions that were present in both the GWAS and RNA-seq networks had a higher probability of being actually involved in sperm quality and used them to build a robust gene interaction network. In addition, in the final network we included genes with RNA abundances correlated with more than four semen traits and miRNAs interacting with the genes on the network. The final network was enriched for genes involved in gamete generation and development, meiotic cell cycle, DNA repair or embryo implantation. Finally, we designed a panel of 73 SNPs based on the GWAS, eGWAS and final network data, that explains between 5% (for sperm cell concentration) and 36% (for percentage of neck abnormalities) of the phenotypic variance of the sperm traits., [Conclusions]: By applying a systems biology approach, we identified genes that potentially affect sperm quality and constructed a SNP panel that explains a substantial part of the phenotypic variance for semen quality in our study and that should be tested in other swine populations to evaluate its relevance for the pig breeding sector.

[Background]: The comparison of expression QTL (eQTL) maps obtained in different tissues is an essential step to understand how gene expression is genetically regulated in a context-dependent manner. In the current work, we have compared the transcriptomic and eQTL profiles of two porcine tissues (skeletal muscle and liver) which typically show highly divergent expression profiles, in 103 Duroc pigs genotyped with the Porcine SNP60 BeadChip (Illumina) and with available microarray-based measurements of hepatic and muscle mRNA levels. Since structural variation could have effects on gene expression, we have also investigated the co-localization of cis-eQTLs with copy number variant regions (CNVR) segregating in this Duroc population., [Results]: The analysis of differential expresssion revealed the existence of 1204 and 1490 probes that were overexpressed and underexpressed in the gluteus medius muscle when compared to liver, respectively (|fold-change| > 1.5, q-value < 0.05). By performing genome scans in 103 Duroc pigs with available expression and genotypic data, we identified 76 and 28 genome-wide significant cis-eQTLs regulating gene expression in the gluteus medius muscle and liver, respectively. Twelve of these cis-eQTLs were shared by both tissues (i.e. 42.8% of the cis-eQTLs identified in the liver were replicated in the gluteus medius muscle). These results are consistent with previous studies performed in humans, where 50% of eQTLs were shared across tissues. Moreover, we have identified 41 CNVRs in a set of 350 pigs from the same Duroc population, which had been genotyped with the Porcine SNP60 BeadChip by using the PennCNV and GADA softwares, but only a small proportion of these CNVRs co-localized with the cis-eQTL signals., [Conclusion]: Despite the fact that there are considerable differences in the gene expression patterns of the porcine liver and skeletal muscle, we have identified a substantial proportion of common cis-eQTLs regulating gene expression in both tissues. Several of these cis-eQTLs influence the mRNA levels of genes with important roles in meat (CTSF) and carcass quality (TAPT1), lipid metabolism (TMEM97) and obesity (MARC2), thus evidencing the practical importance of dissecting the genetic mechanisms involved in their expression.

This study aims identifying candidate genes and pathways associated with feed efficiency (FE) in pigs. Liver and duodenum transcriptomes of 37 gilts showing high and low residual feed intake (RFI) were analysed by RNA-Seq. Gene expression data was explored through differential expression (DE) and weighted gene co-expression network analyses. DE analysis revealed 55 and 112 differentially regulated genes in liver and duodenum tissues, respectively. Clustering genes according to their connectivity resulted in 23 (liver) and 25 (duodenum) modules of genes with a co-expression pattern. Four modules, one in liver (with 444 co-expressed genes) and three in duodenum (gathering 37, 126 and 41 co-expressed genes), were significantly associated with FE indicators. Intra-module analyses revealed tissue-specific candidate genes; 12 of these genes were also identified as DE between individuals with high and low RFI. Pathways enriched by the list of genes showing DE and/or belonging to FE co-expressed modules included response to oxidative stress, inflammation, immune response, lipid metabolism and thermoregulation. Low overlapping between genes identified in duodenum and liver tissues was observed but heat shock proteins were associated to FE in both tissues. Our results suggest tissue-specific rather than common transcriptome regulatory processes associated with FE in pigs.

In this study, we used BovineSNP50 Genotyping BeadChip data to estimate the structure, putative ancestral origin as well as to identify regions with selective sweeps that may have had an important role in the adaptation to tropical conditions of the ‘Charolais de Cuba’ (CHCU) breed. According to a principal component analysis, CHCU samples cluster together with taurine breeds with an estimated 93% of taurus ancestral alleles. Despite the short period since importation, we detected differentiation (Fst = 0.049) between the French Charolaise (CHA) and CHCU. However, CHA breed was the closest breed to CHCU followed by other hybrids breed with a clear CHA origin. Linkage disequilibrium (r2) decay tends to be lower in CHCU compared to CHA probably due to a less intense artificial selection programs of CHCU. Signals of recent adaptation to tropical conditions between CHCU and CHA were identified. Genes mapping within those regions reflect different functions related to immunity, metabolic changes and heat tolerance (CHCU) and muscle development and meat quality (CHA) that may have had an important role in the phenotypic differentiation of these breeds. Further studies will expand our knowledge on the molecular basis of adaptation of cattle to tropical conditions and molecular process associated with meat quality traits.

In addition to its value in meat production, the pig is an interesting animal model for human digestive tract studies due to its physiological similarities. The aim of this study was to describe the microbiome composition, distribution and interaction along the Iberian pig intestinal tract and its role in whole-body energy homeostasis. The V3-V4 region of the 16S rRNA gene was amplified and sequenced from the microbiomes of five gut sections (duodenum, jejunum, ileum, and proximal and distal colon) in thirteen castrated male pigs. A total of 1,669 operational taxonomic units distributed in 179 genera were found among all samples. The two most abundant genera in the small intestine were Lactobacillus and Clostridium, while Prevotella was predominant in the colon. The colon samples were more similar among the pigs and richer in species than the small intestine samples were. In the small intestine, the metagenome prediction pointed to rapid internalization and conversion of the available simple carbohydrates for microbial proliferation and maintenance. In the colon, a competition among anaerobic bacteria for plant polysaccharide degradation to produce short chain fatty acids was found. This study confirms that the energy pathways of the gut microbiome differ along its sections and provides a description of the correlations between genera.

In the present study, liver co-expression networks and expression Genome Wide Association Study (eGWAS) were performed to identify DNA variants and molecular pathways implicated in the functional regulatory mechanisms of meat quality traits in pigs. With this purpose, the liver mRNA expression of 44 candidates genes related with lipid metabolism was analysed in 111 Iberian x Landrace backcross animals. The eGWAS identified 92 eSNPs located in seven chromosomal regions and associated with eight genes: CROT, CYP2U1, DGAT1, EGF, FABP1, FABP5, PLA2G12A, and PPARA. Remarkably, cis-eSNPs associated with FABP1 gene expression which may be determining the C18:2(n-6)/C18:3(n-3) ratio in backfat through the multiple interaction of DNA variants and genes were identified. Furthermore, a hotspot on SSC8 associated with the gene expression of eight genes was identified and the TBCK gene was pointed out as candidate gene regulating it. Our results also suggested that the PI3K-Akt-mTOR pathway plays an important role in the control of the analysed genes highlighting nuclear receptors as the NR3C1 or PPARA. Finally, sex-dimorphism associated with hepatic lipid metabolism was identified with over-representation of female-biased genes. These results increase our knowledge of the genetic architecture underlying fat composition traits.

The pig, Sus scrofa, is a foreign species to the American continent. Although pigs originally introduced in the Americas should be related to those from the Iberian Peninsula and Canary islands, the phylogeny of current creole pigs that now populate the continent is likely to be very complex. Because of the extreme climates that America harbors, these populations also provide a unique example of a fast evolutionary phenomenon of adaptation. Here, we provide a genome wide study of these issues by genotyping, with a 60k SNP chip, 206 village pigs sampled across 14 countries and 183 pigs from outgroup breeds that are potential founders of the American populations, including wild boar, Iberian, international and Chinese breeds. Results show that American village pigs are primarily of European ancestry, although the observed genetic landscape is that of a complex conglomerate. There was no correlation between genetic and geographical distances, neither continent wide nor when analyzing specific areas. Most populations showed a clear admixed structure where the Iberian pig was not necessarily the main component, illustrating how international breeds, but also Chinese pigs, have contributed to extant genetic composition of American village pigs. We also observe that many genes related to the cardiovascular system show an increased differentiation between altiplano and genetically related pigs living near sea level.