Identification of eQTLs using different sets of single nucleotide polymorphisms associated with carcass and body composition traits in pigs

BACKGROUND

Mapping expression quantitative trait loci (eQTLs) in skeletal muscle tissue in pigs is crucial for understanding the relationship between genetic variation and phenotypic expression of carcass traits in meat animals. Therefore, the primary objective of this study was to evaluate the impact of different sets of single nucleotide polymorphisms (SNP), including scenarios removing SNPs pruned for linkage disequilibrium (LD) and SNPs derived from SNP chip arrays and RNA-seq data from liver, brain, and skeletal muscle tissues, on the identification of eQTLs in the Longissimus lumborum tissue, associated with carcass and body composition traits in Large White pigs. The SNPs identified from muscle mRNA were combined with SNPs identified in the brain and liver tissue transcriptomes, as well as SNPs from the GGP Porcine 50 K SNP chip array. Cis- and trans-eQTLs were identified based on the skeletal muscle gene expression level, followed by functional genomic analyses and statistical associations with carcass and body composition traits in Large White pigs.

RESULTS

The number of cis- and trans-eQTLs identified across different sets of SNPs (scenarios) ranged from 261 to 2,539 and from 29 to 13,721, respectively. Furthermore, 6,180 genes were modulated by eQTLs in at least one of the scenarios evaluated. The eQTLs identified were not significantly associated with carcass and body composition traits but were significantly enriched for many traits in the "Meat and Carcass" type QTL. The scenarios with the highest number of cis- (n = 304) and trans- (n = 5,993) modulated genes were the unpruned and LD-pruned SNP set scenarios identified from the muscle transcriptome. These genes include 84 transcription factor coding genes.

CONCLUSIONS

After LD pruning, the set of SNPs identified based on the transcriptome of the skeletal muscle tissue of pigs resulted in the highest number of genes modulated by eQTLs. Most eQTLs are of the trans type and are associated with genes influencing complex traits in pigs, such as transcription factors and enhancers. Furthermore, the incorporation of SNPs from other genomic regions to the set of SNPs identified in the porcine skeletal muscle transcriptome contributed to the identification of eQTLs that had not been identified based on the porcine skeletal muscle transcriptome alone.

Production and meat quality traits of genetically lean immunocastrated pigs naturally divergent for loin tenderness

Warner-Bratzler Shear Force (WBSF) is a quantitative measurement of meat toughness that has great impact on the consumer acceptability of meat. This study was conducted to evaluate growth performance, carcass and meat quality characteristics, and fatty acids profile of longissimus lumborum (LL) intramuscular fat (IMF) of pigs that are genetically divergent for WBSF. Based on WBSF values of the LL from a previous study, 12 immunocastrated male pigs selected from 96 pigs were divided into two groups with high WBSF (53.28 to 42.50 N) and low WBSF (37.27 to 27.79 N). Although high-WBSF pigs tended to have improved (P = 0.08) gain-to-feed ratio, overall performance was similar between WBSF groups. High-WBSF pigs also tended to have higher (P = 0.09) cooling loss and lean percentage as well as decreased (P = 0.08) 10th-rib backfat depth than low-WBSF pigs. Loins from high-WBSF pigs tended to have lower (P = 0.07) IMF content and higher (P = 0.09) cooking loss than low-WBSF pigs. Compared to low-WBSF pigs, IMF of the LL from high-WBSF pigs had lower (P = 0.05) percentage of oleic acid and tended to have a decreased (P = 0.07) percentage of total monounsaturated fatty acids. Loins from pigs with high WBSF tended to have increased (P = 0.09) total polyunsaturated fatty acids (PUFA) content and had higher (P = 0.03) PUFA: saturated fatty acid ratio than low-WBSF pigs. Selecting pigs for pork tenderness could potentially conflict with lean growth efficiency and a healthier fatty acids profile for human consumption.

Brain fatty acid and transcriptome profiles of pig fed diets with different levels of soybean oil

BACKGROUND

The high similarity in anatomical and neurophysiological processes between pigs and humans make pigs an excellent model for metabolic diseases and neurological disorders. Lipids are essential for brain structure and function, and the polyunsaturated fatty acids (PUFA) have anti-inflammatory and positive effects against cognitive dysfunction in neurodegenerative diseases. Nutrigenomics studies involving pigs and fatty acids (FA) may help us in better understanding important biological processes. In this study, the main goal was to evaluate the effect of different levels of dietary soybean oil on the lipid profile and transcriptome in pigs' brain tissue.

RESULTS

Thirty-six male Large White pigs were used in a 98-day study using two experimental diets corn-soybean meal diet containing 1.5% soybean oil (SOY1.5) and corn-soybean meal diet containing 3.0% soybean oil (SOY3.0). No differences were found for the brain total lipid content and FA profile between the different levels of soybean oil. For differential expression analysis, using the DESeq2 statistical package, a total of 34 differentially expressed genes (DEG, FDR-corrected p-value < 0.05) were identified. Of these 34 DEG, 25 are known-genes, of which 11 were up-regulated (log2 fold change ranging from + 0.25 to + 2.93) and 14 were down-regulated (log2 fold change ranging from - 3.43 to -0.36) for the SOY1.5 group compared to SOY3.0. For the functional enrichment analysis performed using MetaCore with the 34 DEG, four pathway maps were identified (p-value < 0.05), related to the ALOX15B (log2 fold change - 1.489), CALB1 (log2 fold change - 3.431) and CAST (log2 fold change + 0.421) genes. A "calcium transport" network (p-value = 2.303e-2), related to the CAST and CALB1 genes, was also identified.

CONCLUSION

The results found in this study contribute to understanding the pathways and networks associated with processes involved in intracellular calcium, lipid metabolism, and oxidative processes in the brain tissue. Moreover, these results may help a better comprehension of the modulating effects of soybean oil and its FA composition on processes and diseases affecting the brain tissue.

RNA-seq transcriptome profiling of pigs' liver in response to diet with different sources of fatty acids

Pigs (Sus scrofa) are an animal model for metabolic diseases in humans. Pork is an important source of fatty acids (FAs) in the human diet, as it is one of the most consumed meats worldwide. The effects of dietary inclusion of oils such as canola, fish, and soybean oils on pig gene expression are mostly unknown. Our objective was to evaluate FA composition, identify changes in gene expression in the liver of male pigs fed diets enriched with different FA profiles, and identify impacted metabolic pathways and gene networks to enlighten the biological mechanisms' variation. Large White male pigs were randomly allocated to one of three diets with 18 pigs in each; all diets comprised a base of corn and soybean meal to which either 3% of soybean oil (SOY), 3% canola oil (CO), or 3% fish oil (FO) was added for a 98-day trial during the growing and finishing phases. RNA sequencing was performed on the liver samples of each animal by Illumina technology for differential gene expression analyses, using the R package DESeq2. The diets modified the FA profile, mainly in relation to polyunsaturated and saturated FAs. Comparing SOY vs. FO, 143 differentially expressed genes (DEGs) were identified as being associated with metabolism, metabolic and neurodegenerative disease pathways, inflammatory processes, and immune response networks. Comparing CO vs. SOY, 148 DEGs were identified, with pathways related to FA oxidation, regulation of lipid metabolism, and metabolic and neurodegenerative diseases. Our results help explain the behavior of genes with differential expression in metabolic pathways resulting from feeding different types of oils in pig diets.

Differential Gene Expression Associated with Soybean Oil Level in the Diet of Pigs

Simple Summary

Findings from the analysis of the pig transcriptome may help to better understand the biological mechanisms that can be modulated by the diet. Thus, the aim of this study was to identify the differentially expressed genes from the skeletal muscle and liver samples of pigs fed diets with two different levels of soybean oil (1.5 or 3%). The FA profile in the tissues was modified by the diet mainly related to monounsaturated (MUFA) and polyunsaturated (PUFA). This nutrigenomics study verified the effect of different levels of soybean oil in the pig diet on the transcriptome profile of skeletal muscle and liver, where the higher level of soybean oil added to the diet led to a higher expression of genes targeting biological processes related to lipid oxidation and consequently to metabolic diseases and inflammation.

Abstract

The aim of this study was to identify the differentially expressed genes (DEG) from the skeletal muscle and liver samples of animal models for metabolic diseases in humans. To perform the study, the fatty acid (FA) profile and RNA sequencing (RNA-Seq) data of 35 samples of liver tissue (SOY1.5, n = 17 and SOY3.0, n = 18) and 36 samples of skeletal muscle (SOY1.5, n = 18 and SOY3.0, n = 18) of Large White pigs were analyzed. The FA profile of the tissues was modified by the diet, mainly those related to monounsaturated (MUFA) and polyunsaturated (PUFA) FA. The skeletal muscle transcriptome analysis revealed 45 DEG (FDR 10%), and the functional enrichment analysis identified network maps related to inflammation, immune processes, and pathways associated with oxidative stress, type 2 diabetes, and metabolic dysfunction. For the liver tissue, the transcriptome profile analysis revealed 281 DEG, which participate in network maps related to neurodegenerative diseases. With this nutrigenomics study, we verified that different levels of soybean oil in the pig diet, an animal model for metabolic diseases in humans, affected the transcriptome profile of skeletal muscle and liver tissue. These findings may help to better understand the biological mechanisms that can be modulated by the diet.

Effect of dietary soybean oil inclusion on liver-related transcription factors in a pig model for metabolic diseases

Dietary fatty acids (FA) are components of the lipids, which contribute to membrane structure, energy input, and biological functions related to cellular signaling and transcriptome regulation. However, the consumers still associate dietary FA with fat deposition and increased occurrence of metabolic diseases such as obesity and atherosclerosis. Previous studies already demonstrated that some fatty acids are linked with inflammatory response, preventing metabolic diseases. To better understand the role of dietary FA on metabolic diseases, for the first time, a study to identify key transcription factors (TF) involved in lipid metabolism and inflammatory response by transcriptome analysis from liver samples of animal models was performed. The key TF were identified by functional enrichment analysis from the list of differentially expressed genes identified in liver samples between 35 pigs fed with 1.5% or 3.0% soybean oil. The functional enrichment analysis detected TF linked to lipid homeostasis and inflammatory response, such as RXRA, EGFR, and SREBP2 precursor. These findings demonstrated that key TF related to lipid metabolism could be modulated by dietary inclusion of soybean oil. It could contribute to nutrigenomics research field that aims to elucidate dietary interventions in animal and human health, as well as to drive food technology and science.

Dietary fat content and fiber type modulate hind gut microbial community and metabolic markers in the pig

Obesity leads to changes in the gut microbial community which contribute to the metabolic dysregulation in obesity. Dietary fat and fiber affect the caloric density of foods. The impact of dietary fat content and fiber type on the microbial community in the hind gut is unknown. Effect of dietary fat level and fiber type on hindgut microbiota and volatile fatty acid (VFA) profiles was investigated. Expression of metabolic marker genes in the gut, adipose tissue and liver was determined. A 2×2 experiment was conducted in pigs fed at two dietary fat levels (5% or 17.5% swine grease) and two fiber types (4% inulin, fermentable fructo-oligosaccharide or 4% solka floc, non-fermentable cellulose). High fat diets (HFD) resulted in a higher (P<0.05) total body weight gain, feed efficiency and back fat accumulation than the low fat diet. Feeding of inulin, but not solka floc, attenuated (P<0.05) the HFD-induced higher body weight gain and fat mass accumulation. Inulin feeding tended to lead to higher total VFA production in the cecum and resulted in a higher (P<0.05) expression of acyl coA oxidase (ACO), a marker of peroxisomal β-oxidation. Inulin feeding also resulted in lower expression of sterol regulatory element binding protein 1c (SREBP-1c), a marker of lipid anabolism. Bacteria community structure characterized by DGGE analysis of PCR amplified 16S rRNA gene fragments showed that inulin feeding resulted in greater bacterial population richness than solka floc feeding. Cluster analysis of pairwise Dice similarity comparisons of the DGGE profiles showed grouping by fiber type but not the level of dietary fat. Canonical correspondence analysis (CCA) of PCR- DGGE profiles showed that inulin feeding negatively correlated with back fat thickness. This study suggests a strong interplay between dietary fat level and fiber type in determining susceptibility to obesity.