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.

Identifying miRNA-mRNA regulatory networks on extreme n-6/n-3 polyunsaturated fatty acid ratio expression profiles in porcine skeletal muscle.

Omega-3 (n-3) and omega-6 (n-6) polyunsaturated fatty acids (PUFAs) are essential fatty acids with antagonistic inflammatory functions that play vital roles in metabolic health and immune response. Current commercial swine diets tend to over-supplement with n-6 PUFAs, which may increase the likelihood of developing inflammatory diseases and affect the overall well-being of the animals. However, it is still poorly understood how n-6/n-3 PUFA ratios affect the porcine transcriptome expression and how messenger RNAs (mRNAs) and microRNAs (miRNAs) might regulate biological processes related to PUFA metabolism. On account of this, we selected a total of 20 Iberian × Duroc crossbred pigs with extreme values for n-6/n-3 FA ratio (10 high vs 10 low), and longissimus dorsi muscle samples were used to identify differentially expressed mRNAs and miRNAs. The observed differentially expressed mRNAs were associated to biological pathways related to muscle growth and immunomodulation, while the differentially expressed microRNAs (ssc-miR-30a-3p, ssc-miR-30e-3p, ssc-miR-15b and ssc-miR-7142-3p) were correlated to adipogenesis and immunity. Relevant miRNA-to-mRNA regulatory networks were also predicted (i.e., mir15b to ARRDC3; mir-7142-3p to METTL21C), and linked to lipolysis, obesity, myogenesis, and protein degradation. The n-6/n-3 PUFA ratio differences in pig skeletal muscle revealed genes, miRNAs and enriched pathways involved in lipid metabolism, cell proliferation and inflammation.

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.

Meta-Analysis and Validation of a Colorectal Cancer Risk Prediction Model Using Deep Sequenced Fecal Metagenomes.

The gut microbiome is a potential modifiable risk factor for colorectal cancer (CRC). We re-analyzed all eight previously published stool sequencing data and conducted an MWAS meta-analysis. We used cross-validated LASSO predictive models to identify a microbiome signature for predicting the risk of CRC and precancerous lesions. These models were validated in a new study, Colorectal Cancer Screening (COLSCREEN), including 156 participants that were recruited in a CRC screening context. The MWAS meta-analysis identified 95 bacterial species that were statistically significantly associated with CRC (FDR < 0.05). The LASSO CRC predictive model obtained an area under the receiver operating characteristic curve (aROC) of 0.81 (95%CI: 0.78−0.83) and the validation in the COLSCREEN dataset was 0.75 (95%CI: 0.66−0.84). This model selected a total of 32 species. The aROC of this CRC-trained model to predict precancerous lesions was 0.52 (95%CI: 0.41−0.63). We have identified a signature of 32 bacterial species that have a good predictive accuracy to identify CRC but not precancerous lesions, suggesting that the identified microbes that were enriched or depleted in CRC are merely a consequence of the tumor. Further studies should focus on CRC as well as precancerous lesions with the intent to implement a microbiome signature in CRC screening programs.

Phytogenic Compounds Supplemented to Gestating Hyperprolific Sows Affects the Gut Health-Related Gene Expression and Histological Responses in Neonate Piglets.

This research aims to determine whether a specific blend of phytogenic compounds (BPC) supplemented in gestating hyperprolific sow diets can promote prenatal maternal effects in terms of piglet gut function and morphology. Twenty-eight (Landrace × Yorkshire) gilts and sows (parity 0 to 7) were randomly distributed by parity number and body weight into two dietary treatments: unsupplemented Control (CON) (n = 14) or CON diet supplemented with 1 g/kg feed of BPC during gestation (n = 14). The BPC supplementation during gestation of sows downregulated the neonate piglets' jejunal genes involved in oxidation (SOD2) and nutrient transport (SLC16A1/MCT1, SLC11A2/DMT1, and SLC39A/ZIP4), while IFNG and CLDN4 related to immune response and barrier function, respectively, were upregulated (q < 0.10). In addition, the jejunal villus height and the ratio of the villus height to crypt depth tended to increase (p < 0.10), while goblet cell volume density was higher (p < 0.05) in BPC compared to CON. In conclusion, dietary supplementation of BPC in gestating diets for hyperprolific sows influences neonatal histomorphology and expression of genes related to the intestinal function and health.

Expression analysis of porcine miR-33a/b in liver, adipose tissue and muscle and its potential role in fatty acid metabolism.

mir-33a and mir-33b are co-transcribed with the SREBF2 and SREBF1 transcription factors, respectively. The main role of SREBF1 is the regulation of genes involved in fatty acid metabolism, while SREBF2 regulates genes participating in cholesterol biosynthesis and uptake. Our objective was to study the expression of both miR-33a and miR-33b, together with their host SREBF genes, in liver, adipose tissue and muscle to better understand the role of miR-33a/b in the lipid metabolism of pigs. In our study, the expression of miR-33a, miR-33b and SREBF2 in liver, adipose tissue, and muscle was studied in 42 BC1_LD (25% Iberian x 75% Landrace backcross) pigs by RT-qPCR. In addition, the expression of in-silico predicted target genes and fatty acid composition traits were correlated with the miR-33a/b expression. We observed different tissue expression patterns for both miRNAs. In adipose tissue and muscle a high correlation between miR-33a and miR-33b expression was found, whereas a lower correlation was observed in liver. The expression analysis of in-silico predicted target-lipid related genes showed negative correlations between miR-33b and CPT1A expression in liver. Conversely, positive correlations between miR-33a and PPARGC1A and USF1 gene expression in liver were observed. Lastly, positive and negative correlations between miR-33a/b expression and saturated fatty acid (SFA) and polyunsaturated fatty acid (PUFA) content, respectively, were identified. Overall, our results suggested that both miRNAs are differentially regulated and have distinct functions in liver, in contrast to muscle and adipose tissue. Furthermore, the correlations between miR-33a/b expression both with the expression of in-silico predicted target-lipid related genes and with fatty acid composition, opens new avenues to explore the role of miR33a/b in the regulation of lipid metabolism.

Porcine Digestible Peptides (PDP) in Weanling Diets Regulates the Expression of Genes Involved in Gut Barrier Function, Immune Response and Nutrient Transport in Nursery Pigs.

This study was conducted to investigate the effects of dietary supplementation of porcine digestible peptides (PDP), spray-dried plasma (SDP), or a combination of both, on growth performance and the expression of genes related to intestinal function of weaned pigs. A total of 180 piglets (trial 1) and 198 piglets (trial 2) were used to evaluate the partial substitution of soybean ingredients with 2% SDP or 2% PDP (trial 1), and with 3% SDP or the combination of 1% SDP and 2% PDP (SDP-PDP; trial 2) during the pre-starter period (0-14 days). The gene expression of 56 genes was quantified in a qPCR platform in jejunum and ileum samples obtained from piglets 14 d after weaning (trial 2). Piglets fed SDP, PDP and SDP-PDP had a higher body weight (BW), average daily gain (ADG) and feed efficiency (G:F) than the soybean control on day 14 (p < 0.05). In addition, the combination of SDP and PDP upregulated ten genes in jejunum samples (p < 0.05) related to intestinal function. More research is needed to confirm that gene expression upregulation by PDP in combination with SDP has an impact on intestinal function and to elucidate its underlying mechanisms.

Identification of strong candidate genes for backfat and intramuscular fatty acid composition in three crosses based on the Iberian pig.

Meat quality has an important genetic component and can be modified by the fatty acid (FA) composition and the amount of fat contained in adipose tissue and muscle. The present study aimed to find genomic regions associated with the FA composition in backfat and muscle (longissimus dorsi) in 439 pigs with three different genetic backgrounds but having the Iberian breed in common. Genome-wide association studies (GWAS) were performed between 38,424 single-nucleotide polymorphisms (SNPs) covering the pig genome and 60 phenotypic traits related to backfat and muscle FA composition. Nine significant associated regions were found in backfat on the Sus scrofa chromosomes (SSC): SSC1, SSC2, SSC4, SSC6, SSC8, SSC10, SSC12, and SSC16. For the intramuscular fat, six significant associated regions were identified on SSC4, SSC13, SSC14, and SSC17. A total of 52 candidate genes were proposed to explain the variation in backfat and muscle FA composition traits. GWAS were also reanalysed including SNPs on five candidate genes (ELOVL6, ELOVL7, FADS2, FASN, and SCD). Regions and molecular markers described in our study may be useful for meat quality selection of commercial pig breeds, although several polymorphisms were breed-specific, and further analysis would be needed to evaluate possible causal mutations.

Identification of eQTLs associated with lipid metabolism in Longissimus dorsi muscle of pigs with different genetic backgrounds.

Intramuscular fat content and its fatty acid composition affect porcine meat quality and its nutritional value. The present work aimed to identify genomic variants regulating the expression in the porcine muscle (Longissimus dorsi) of 45 candidate genes for lipid metabolism and fatty acid composition in three experimental backcrosses based on the Iberian breed. Expression genome-wide association studies (eGWAS) were performed between the muscle gene expression values, measured by real-time quantitative PCR, and the genotypes of 38,426 SNPs distributed along all chromosomes. The eGWAS identified 186 eSNPs located in ten Sus scrofa regions and associated with the expression of ACSM5, ACSS2, ATF3, DGAT2, FOS and IGF2 (FDR < 0.05) genes. Two expression quantitative trait loci (eQTLs) for IGF2 and ACSM5 were classified as cis-acting eQTLs, suggesting a mutation in the same gene affecting its expression. Conversely, ten eQTLs showed trans-regulatory effects on gene expression. When the eGWAS was performed for each backcross independently, only three common trans-eQTL regions were observed, indicating different regulatory mechanisms or allelic frequencies among the breeds. In addition, hotspot regions regulating the expression of several genes were detected. Our results provide new data to better understand the functional regulatory mechanisms of lipid metabolism genes in muscle.

Analysis of porcine IGF2 gene expression in adipose tissue and its effect on fatty acid composition.

IGF2:g.3072G>A polymorphism has been described as the causal mutation of a maternally imprinted QTL for muscle growth and fat deposition in pigs. The objective of the current work was to study the association between the IGF2:g.3072G>A polymorphism and the IGF2 gene expression and its effect on fatty acid composition in adipose tissue in different pig genetic backgrounds. A cis-eQTL region associated with the IGF2 mRNA expression in adipose tissue was identified in an eGWAS with 355 animals. The IGF2 gene was located in this genomic interval and IGF2g.3072G>A was the most significant SNP, explaining a 25% of the gene expression variance. Significant associations between IGF2:g.3072G>A polymorphism and oleic (C18:1(n-9); p-value = 4.18x10-07), hexadecanoic (C16:1(n-9); p-value = 4.04x10-07), linoleic (C18:2(n-6); p-value = 6.44x10-09), α-linoleic (C18:3(n-3); p-value = 3.30x10-06), arachidonic (C20:4(n-6); p-value = 9.82x10-08) FAs and the MUFA/PUFA ratio (p-value = 2.51x10-9) measured in backfat were identified. Animals carrying the A allele showed an increase in IGF2 gene expression and higher PUFA and lower MUFA content. However, in additional studies was observed that there could be other proximal genetic variants affecting FA composition in adipose tissue. Finally, no differences in the IGF2 gene expression in adipose tissue were found between heterozygous animals classified according to the IGF2:g.3072G>A allele inherited from the father (APGM or AMGP). However, pyrosequencing analysis revealed that there is imprinting of the IGF2 gene in muscle and adipose tissues, with stronger differences among the paternally and maternally inherited alleles in muscle. Our results suggested that IGF2:g.3072G>A polymorphism plays an important role in the regulation of IGF2 gene expression and can be involved in the fatty acid composition in adipose tissue. In both cases, further studies are still needed to deepen the mechanism of regulation of IGF2 gene expression in adipose tissue and the IGF2 role in FA composition.

Indel detection from Whole Genome Sequencing data and association with lipid metabolism in pigs.

The selection in commercial swine breeds for meat-production efficiency has been increasing among the past decades, reducing the intramuscular fat content, which has changed the sensorial and technological properties of pork. Through processes of natural adaptation and selective breeding, the accumulation of mutations has driven the genetic divergence between pig breeds. The most common and well-studied mutations are single-nucleotide polymorphisms (SNPs). However, insertions and deletions (indels) usually represents a fifth part of the detected mutations and should also be considered for animal breeding. In the present study, three different programs (Dindel, SAMtools mpileup, and GATK) were used to detect indels from Whole Genome Sequencing data of Iberian boars and Landrace sows. A total of 1,928,746 indels were found in common with the three programs. The VEP tool predicted that 1,289 indels may have a high impact on protein sequence and function. Ten indels inside genes related with lipid metabolism were genotyped in pigs from three different backcrosses with Iberian origin, obtaining different allelic frequencies on each backcross. Genome-Wide Association Studies performed in the Longissimus dorsi muscle found an association between an indel located in the C1q and TNF related 12 (C1QTNF12) gene and the amount of eicosadienoic acid (C20:2(n-6)).

Association between the pig genome and its gut microbiota composition.

The gut microbiota has been evolving with its host along the time creating a symbiotic relationship. In this study, we assess the role of the host genome in the modulation of the microbiota composition in pigs. Gut microbiota compositions were estimated through sequencing the V3-V4 region of the 16S rRNA gene from rectal contents of 285 pigs. A total of 1,261 operational taxonomic units were obtained and grouped in 18 phyla and 101 genera. Firmicutes (45.36%) and Bacteroidetes (37.47%) were the two major phyla obtained, whereas at genus level Prevotella (7.03%) and Treponema (6.29%) were the most abundant. Pigs were also genotyped with a high-throughput method for 45,508 single nucleotide polymorphisms that covered the entire pig genome. Subsequently, genome-wide association studies were made among the genotypes of these pigs and their gut microbiota composition. A total of 52 single-nucleotide polymorphisms distributed in 17 regions along the pig genome were associated with the relative abundance of six genera; Akkermansia, CF231, Phascolarctobacterium, Prevotella, SMB53, and Streptococcus. Our results suggest 39 candidate genes that may be modulating the microbiota composition and manifest the association between host genome and gut microbiota in pigs.

Characterization of bacterial microbiota compositions along the intestinal tract in pigs and their interactions and functions.

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.

Expression analysis of candidate genes for fatty acid composition in adipose tissue and identification of regulatory regions.

The aim of this work was to study the genetic basis of the backfat expression of lipid-related genes associated with meat quality traits in pigs. We performed a genome-wide association study with the backfat gene expression measured in 44 genes by qPCR and the PorcineSNP60 BeadChip genotypes in 115 Iberian x Landrace backcross animals. A total of 193 expression-associated SNPs located in 19 chromosomal regions were associated with expression levels of ACSM5, ELOVL6, FABP4, FADS2, and SLC27A4 genes. Three expression quantitative trail loci (eQTLs) corresponding to ACSM5, FABP4, and FADS2 were classified as cis-acting eQTLs, whereas the remaining 16 eQTLs have trans-regulatory effects. Remarkably, a SNP in the ACSM5 promoter region and a SNP in the 3'UTR region of FABP4 were the most associated polymorphisms with the ACSM5 and FABP4 expression levels, respectively. Moreover, relevant lipid-related genes mapped in the trans-eQTLs regions associated with the ACSM5, FABP4, FADS2, and SLC27A4 genes. Interestingly, a trans-eQTL hotspot on SSC13 regulating the gene expression of ELOVL6, ELOLV5, and SCD, three important genes implicated in the elongation and desaturation of fatty acids, was identified. These findings provide new data to further understand the functional regulatory mechanisms implicated in the variation of fatty acid composition in pigs.