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

Cold exposure affects glucose metabolism, lipid droplet deposition and mitophagy in skeletal muscle of newborn goats

Cold exposure is a common stressor for newborn goats. Skeletal muscle plays an important role in maintaining whole-body homeostasis of glucose and lipid metabolism. However, the molecular mechanisms underlying regulation of skeletal muscle of newborn goats by cold exposure remains unclear. In this study, we found a significant increase (P < 0.01) in serum glucagon levels after 24 h of cold exposure (COLD, 6°C), while glucose and insulin concentrations were significantly decreased (P < 0.01) compared to room temperature (RT, 25°C). Additionally, we found that cold exposure reduced glycogen content (P < 0.01) in skeletal muscle. Pathway enrichment analysis revealed that cold exposure activated skeletal muscle glucose metabolism pathways (including insulin resistance and the insulin signaling pathway) and mitophagy-related pathways. Cold exposure up-regulated the expression of genes involved in fatty acid and triglyceride synthesis, promoting skeletal muscle lipid deposition. Notably, cold exposure induced mitophagy in skeletal muscle.

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.

Identification of genomic regions, genetic variants and gene networks regulating candidate genes for lipid metabolism in pig muscle

The intramuscular fat content and fatty acid composition of porcine meat have a significant impact on its quality and nutritional value. This research aimed to investigate the expression of 45 genes involved in lipid metabolism in the longissimus dorsi muscle of three experimental pig backcrosses, with a 25% of Iberian background. To achieve this objective, we conducted an expression Genome-Wide Association Study (eGWAS) using gene expression levels in muscle measured by high-throughput real-time qPCR for 45 target genes and genotypes from the PorcineSNP60 BeadChip or Axiom Porcine Genotyping Array and 65 single nucleotide polymorphisms (SNPs) located in 20 genes genotyped by a custom-designed Taqman OpenArray in a cohort of 354 animals. The eGWAS analysis identified 301 eSNPs associated with 18 candidate genes (ANK2, APOE, ARNT, CIITA, CPT1A, EGF, ELOVL6, ELOVL7, FADS3, FASN, GPAT3, NR1D2, NR1H2, PLIN1, PPAP2A, RORA, RXRA and UCP3). Three cis-eQTL (expression quantitative trait loci) were identified for GPAT3, RXRA, and UCP3 genes, which indicates that a genetic polymorphism proximal to the same gene is affecting its expression. Furthermore, 24 trans-eQTLs were detected, and eight candidate regulatory genes were located in these genomic regions. Additionally, two trans-regulatory hotspots in Sus scrofa chromosomes 13 and 15 were identified. Moreover, a co-expression analysis performed on 89 candidate genes and the fatty acid composition revealed the regulatory role of four genes (FABP5, PPARG, SCD, and SREBF1). These genes modulate the levels of α-linolenic, arachidonic, and oleic acids, as well as regulating the expression of other candidate genes associated with lipid metabolism. The findings of this study offer novel insights into the functional regulatory mechanism of genes involved in lipid metabolism, thereby enhancing our understanding of this complex biological process.

Characterization of novel loci controlling seed oil content in Brassica napus by marker metabolite-based multi-omics analysis

BACKGROUND

Seed oil content is an important agronomic trait of Brassica napus (B. napus), and metabolites are considered as the bridge between genotype and phenotype for physical traits.

RESULTS

Using a widely targeted metabolomics analysis in a natural population of 388 B. napus inbred lines, we quantify 2172 metabolites in mature seeds by liquid chromatography mass spectrometry, in which 131 marker metabolites are identified to be correlated with seed oil content. These metabolites are then selected for further metabolite genome-wide association study and metabolite transcriptome-wide association study. Combined with weighted correlation network analysis, we construct a triple relationship network, which includes 21,000 edges and 4384 nodes among metabolites, metabolite quantitative trait loci, genes, and co-expression modules. We validate the function of BnaA03.TT4, BnaC02.TT4, and BnaC05.UK, three candidate genes predicted by multi-omics analysis, which show significant impacts on seed oil content through regulating flavonoid metabolism in B. napus.

CONCLUSIONS

This study demonstrates the advantage of utilizing marker metabolites integrated with multi-omics analysis to dissect the genetic basis of agronomic traits in crops.

CREG1 stimulates AMPK phosphorylation and glucose uptake in skeletal muscle cells

The cellular repressor of adenovirus early region 1A-stimulated gene 1 (CREG1) is a secreted glycoprotein involved in cell differentiation and energy metabolism. It also binds to insulin-like growth factor 2 receptor (IGF2R), a protein implicated in muscle regeneration. However, whether CREG1 regulates the regeneration and metabolism of skeletal muscles via IGF2R remains unclear. This study investigates the role of CREG1 in skeletal muscle regeneration and glucose uptake in C2C12 myotubes and a cardiotoxin (CTX)-induced mouse skeletal muscle regeneration model. CTX-treated skeletal muscle showed significantly higher levels of IGF2R, CREG1, phospho-AMPKα Thr¹⁷², and GLUT4 proteins. Similarly, treatment of myotubes with CREG1 also stimulated AMPKα phosphorylation and GLUT4 expression. CREG1-induced AMPKα phosphorylation and 2DG uptake in myotubes were suppressed by IGF2R knockdown and Compound C, an AMPK inhibitor. These results suggest that CREG1 stimulates glucose uptake in skeletal muscles partially through AMPK activation. Hence, CREG1 plays an essential role in muscle regeneration by affecting glucose metabolism in skeletal muscles.

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.

An insight into the runs of homozygosity distribution and breed differentiation in Mangalitsa pigs

Mangalitsa pigs exhibit three distinct coat color patterns based on which they are described as Red, Blond, and Swallow-bellied. The current study investigated genome-wide diversity and selection signatures in the three breeds using fixation index, runs of homozygosity and population structure analyses. The analyses were originally based on quality-controlled data on 77 Mangalitsa animals from Germany, including 23 Blond, 30 Swallow-bellied and 24 Red Mangalitsa genotyped with a customized version of the ProcineSNP60 v2 Genotyping Bead Chip. Also, 20 Hungarian Mangalitsa genotypes were included as outgroup data for comparison. Estimates of observed heterozygosity were 0.27, 0.28, and 0.29, and inbreeding coefficients estimated based on runs of homozygosity were 24.11%, 20.82%, and 16.34% for Blond, Swallow-bellied and Red Mangalitsa, respectively. ROH islands were detected in all breeds, however, none of these were shared amongst them. The KIF16B gene previously reported to play a role in synaptic signaling was found in a ROH island (SSC17: 16–26) in Swallow-bellied Mangalitsa. The same gene was found to harbor a significantly differentiated SNP (MARC0032380) while contrasting either Blond or Red to Swallow-belied Mangalitsa. In the Red Mangalitsa, some ROH islands were associated with genes that play a role in meat quality traits, i.e., ABCA12, VIL1, PLSCR5, and USP37. Our population structure analysis highlighted a separation of the three breeds, but also showed the closest relatedness between Red and Blond Mangalitsa pigs. Findings of this study improve our understanding of the diversity in the three breeds of Mangalitsa pigs.

Comparative Transcriptomic Analysis of mRNAs, miRNAs and lncRNAs in the Longissimus dorsi Muscles between Fat-Type and Lean-Type Pigs

In pigs, meat quality and production are two important traits affecting the pig industry and human health. Compared to lean-type pigs, fat-type pigs contain higher intramuscular fat (IMF) contents, better taste and nutritional value. To uncover genetic factors controlling differences related to IMF in pig muscle, we performed RNA-seq analysis on the transcriptomes of the Longissimus dorsi (LD) muscle of Laiwu pigs (LW, fat-type pigs) and commercial Duroc × Landrace × Yorkshire pigs (DLY, lean-type pigs) at 150 d to compare the expression profiles of mRNA, miRNA and lncRNA. A total of 225 mRNAs, 12 miRNAs and 57 lncRNAs were found to be differentially expressed at the criteria of |log2(foldchange)| > 1 and q < 0.05. The mRNA expression of LDHB was significantly higher in the LD muscle of LW compared to DLY pigs with log2(foldchange) being 9.66. Using protein interaction prediction method, we identified more interactions of estrogen-related receptor alpha (ESRRA) associated with upregulated mRNAs, whereas versican (VCAN) and proenkephalin (PENK) were associated with downregulated mRNAs in LW pigs. Integrated analysis on differentially expressed (DE) mRNAs and miRNAs in the LD muscle between LW and DLY pigs revealed two network modules: between five upregulated mRNA genes (GALNT15, FKBP5, PPARGC1A, LOC110258214 and LOC110258215) and six downregulated miRNA genes (ssc-let-7a, ssc-miR190-3p, ssc-miR356-5p, ssc-miR573-5p, ssc-miR204-5p and ssc-miR-10383), and between three downregulated DE mRNA genes (IFRD1, LOC110258600 and LOC102158401) and six upregulated DE miRNA genes (ssc-miR1379-3p, ssc-miR1379-5p, ssc-miR397-5p, ssc-miR1358-5p, ssc-miR299-5p and ssc-miR1156-5p) in LW pigs. Based on the mRNA and ncRNA binding site targeting database, we constructed a regulatory network with miRNA as the center and mRNA and lncRNA as the target genes, including GALNT15/ssc-let-7a/LOC100523888, IFRD1/ssc-miR1379-5p/CD99, etc., forming a ceRNA network in the LD muscles that are differentially expressed between LW and DLY pigs. Collectively, these data may provide resources for further investigation of molecular mechanisms underlying differences in meat traits between lean- and fat-type pigs.

Transcriptome-metabolome analysis reveals how sires affect meat quality in hybrid sheep populations

Crossbreeding improves and enhances meat quality and is widely used in sheep production; however, the molecular mechanisms underlying the meat quality of various crossbred sheep remain unknown. In this study, male Southdown, Suffolk and Hu sheep were crossbred with female Hu sheep, and the transcriptomes and metabolomes of the longissimus dorsi muscle of the F1 generation were sequenced to explore how different sire breeds affect meat quality. The results showed that 631 differentially expressed genes and 119 significantly altered metabolites contributed to muscle development characteristics and meat quality-related diversity (P < 0.05). These genes and metabolites were significantly enriched in lipid metabolism pathways, including arachidonic acid metabolism and PPAR signaling. Several candidate genes were associated with muscle growth, such as MYLK3, MYL10, FIGN, MYH8, MYOM3, LMCD1, and FLRT1. Among these, MYH8 and MYL10 participated in regulating muscle growth and development and were correlated with meat quality-related fatty acid levels (|r| > 0.5 and p < 0.05). We selected mRNA from four of these genes to verify the accuracy of the sequencing data via qRT-PCR. Our findings provide further insight into the key genes and metabolites involved in muscle growth and meat quality in hybrid sheep populations.

Evidence for and localization of proposed causative variants in cattle and pig genomes

BACKGROUND

This paper reviews the localization of published potential causative variants in contemporary pig and cattle reference genomes, and the evidence for their causality. In spite of the difficulties inherent to the identification of causative variants from genetic mapping and genome-wide association studies, researchers in animal genetics have proposed putative causative variants for several traits relevant to livestock breeding.

RESULTS

For this review, we read the literature that supports potential causative variants in 13 genes (ABCG2, DGAT1, GHR, IGF2, MC4R, MSTN, NR6A1, PHGK1, PRKAG3, PLRL, RYR1, SYNGR2 and VRTN) in cattle and pigs, and localized them in contemporary reference genomes. We review the evidence for their causality, by aiming to separate the evidence for the locus, the proposed causative gene and the proposed causative variant, and report the bioinformatic searches and tactics needed to localize the sequence variants in the cattle or pig genome.

CONCLUSIONS

Taken together, there is usually good evidence for the association at the locus level, some evidence for a specific causative gene at eight of the loci, and some experimental evidence for a specific causative variant at six of the loci. We recommend that researchers who report new potential causative variants use referenced coordinate systems, show local sequence context, and submit variants to repositories.

Identification of Distant Regulatory Elements Using Expression Quantitative Trait Loci Mapping for Heat-Responsive Genes in Oysters

Many marine ectotherms, especially those inhabiting highly variable intertidal zones, develop high phenotypic plasticity in response to rapid climate change by modulating gene expression levels. Herein, we examined the regulatory architecture of heat-responsive gene expression plasticity in oysters using expression quantitative trait loci (eQTL) analysis. Using a backcross family of Crassostrea gigas and its sister species Crassostrea angulata under acute stress, 56 distant regulatory regions accounting for 6–26.6% of the gene expression variation were identified for 19 heat-responsive genes. In total, 831 genes and 164 single nucleotide polymorphisms (SNPs) that could potentially regulate expression of the target genes were screened in the eQTL region. The association between three SNPs and the corresponding target genes was verified in an independent family. Specifically, Marker13973 was identified for heat shock protein (HSP) family A member 9 (HspA9). Ribosomal protein L10a (RPL10A) was detected approximately 2 kb downstream of the distant regulatory SNP. Further, Marker14346-48 and Marker14346-85 were in complete linkage disequilibrium and identified for autophagy-related gene 7 (ATG7). Nuclear respiratory factor 1 (NRF1) was detected approximately 3 kb upstream of the two SNPs. These results suggested regulatory relationships between RPL10A and HSPA9 and between NRF1 and ATG7. Our findings indicate that distant regulatory mutations play an important role in the regulation of gene expression plasticity by altering upstream regulatory factors in response to heat stress. The identified eQTLs provide candidate biomarkers for predicting the persistence of oysters under future climate change scenarios.

Candidate gene screening for lipid deposition using combined transcriptomic and proteomic data from Nanyang black pigs

Background

Lower selection intensities in indigenous breeds of Chinese pig have resulted in obvious genetic and phenotypic divergence. One such breed, the Nanyang black pig, is renowned for its high lipid deposition and high genetic divergence, making it an ideal model in which to investigate lipid position trait mechanisms in pigs. An understanding of lipid deposition in pigs might improve pig meat traits in future breeding and promote the selection progress of pigs through modern molecular breeding techniques. Here, transcriptome and tandem mass tag-based quantitative proteome (TMT)-based proteome analyses were carried out using longissimus dorsi (LD) tissues from individual Nanyang black pigs that showed high levels of genetic variation.

Results

A large population of Nanyang black pigs was phenotyped using multi-production trait indexes, and six pigs were selected and divided into relatively high and low lipid deposition groups. The combined transcriptomic and proteomic data identified 15 candidate genes that determine lipid deposition genetic divergence. Among them, FASN, CAT, and SLC25A20 were the main causal candidate genes. The other genes could be divided into lipid deposition-related genes (BDH2, FASN, CAT, DHCR24, ACACA, GK, SQLE, ACSL4, and SCD), PPARA-centered fat metabolism regulatory factors (PPARA, UCP3), transcription or translation regulators (SLC25A20, PDK4, CEBPA), as well as integrin, structural proteins, and signal transduction-related genes (EGFR).

Conclusions

This multi-omics data set has provided a valuable resource for future analysis of lipid deposition traits, which might improve pig meat traits in future breeding and promote the selection progress in pigs, especially in Nanyang black pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07764-2.

Proteomic and lipidomic analyses reveal saturated fatty acids, phosphatidylinositol, phosphatidylserine, and associated proteins contributing to intramuscular fat deposition

Intramuscular fat (IMF) content is an important factor in porcine meat quality. Previous studies have screened multiple candidate genes related to IMF deposition, but the lipids that affect IMF deposition and their lipid-protein network remain unknown. In this study, we performed proteomic and lipidomic analyses of the longissimus dorsi (LD) muscle from high-IMF (IMFH) and low-IMF (IMF-L) groups of Xidu black pigs. Eighty-eight proteins and 143 lipids were differentially abundant between the groups. The differentially abundant proteins were found to be involved in cholesterol metabolism, the PPAR signaling pathway, and ferroptosis. The triacylglycerols (TAGs) upregulated in the IMF-H group were mainly shown to be synthesized by saturated fatty acids (SFAs), while the downregulated TAGs were mainly synthesized by polyunsaturated fatty acids (PUFAs). All differentially abundant phosphatidylinositols (PIs) and phosphatidylserines (PSs) were found to be upregulated in the IMF-H group. A correlation analysis of the proteomic and lipidomic revealed candidate proteins (APOA4, VDAC3, PRNP, CTSB, GSPT1) related to TAG, PI, and PS lipids. These results revealed differences in proteins and lipids between the IMF-H and IMF-L groups, which represent new candidate proteins and lipids that should be investigated to determine the molecular mechanisms controlling IMF deposition in pigs. SIGNIFICANCE: Intramuscular fat (IMF) is a key factor affecting meat quality, and meat with a higher IMF content can have a better flavor. In this study, proteomic results show that the ferroptosis pathway, including the PRNP, VDAC3 and CP proteins, affects IMF deposition. Lipid composition is the key factor affecting IMF deposition, but there are few reports on this. In this study, through lipidomic analysis, we suggest that saturated fatty acid (SFA), phosphatidylinositol (PI), and phosphatidylserine (PS) may contribute to IMF deposition. A correlation analysis reveals the potential regulatory network between lipids and proteins. This study clarifies the difference in protein and lipid compositions in longissimus dorsi (LD) muscle with high and low IMF contents. This information suggests that it would be beneficial to increase the intramuscular fat content of pork not only from a genetic perspective but also from a nutritional perspective.