DOI: https://doi.org/10.21203/rs.3.rs-754325/v2
Background: Subcutaneous fat and intramuscular fat (IMF) deposition are closely related to meat production and pork quality. The Dingyuan pig is a local pig breed in Anhui Province, China, that has great potential for fat deposition. Individuals with extreme subcutaneous fat and intramuscular fat content can be found in this breed, which provides a good study system for investigating the molecular mechanisms regulating these two types of fat deposit.
Results: In this study, we used RNA-Seq and tandem mass tags-based proteomics to analyze the key pathways and genes that specifically regulate subcutaneous fat and intramuscular fat deposition in Dingyuan pigs. We identified 191 differentially expressed genes (DEGs) and 61 differentially abundant proteins (DAPs) in the high backfat thickness (HBF) and low backfat thickness (LBF) groups. In the high intramuscular fat and low intramuscular fat groups, we found 85 DEGs and 12 DAPs. The gene ontology and KEGG pathway enrichment analysis showed that the DEGs and DAPs in the backfat groups were mainly involved in various metabolic pathways, such as those related to carbohydrates, amino acids, esters, and fatty acids, whereas the DEGs and DAPs of the IMF groups were involved in a wide range of signaling pathways, including metabolic pathways, the insulin pathway, ketone body synthesis and degradation, longevity, and some disease-related pathways. Among the genes related to the metabolic pathways of carbohydrates, amino acids, esters, and fatty acids, we found 26 candidate genes that specifically regulate subcutaneous fat deposition and 7 genes that specifically regulate IMF deposition in Dingyuan pigs.
Conclusion: Our data show that subcutaneous fat deposition and IMF deposition are regulated by the same genes, but there are also genes that specifically regulate these two fat depositions. Our data provide insights into the mechanisms of pig fat deposition.
China's pig breeds are rich in genetic resources and have good meat quality. However, they also suffer from problems such as a slow growth rate, high fat deposits, and low lean meat percentages [1]. Commercial pig breeds have improved lean meat percentages and reduced subcutaneous fat deposits; however, they simultaneously have low intramuscular fat content and poor meat quality [2–5]. Porcine subcutaneous fat is quantified by backfat thickness, which is negatively correlated with lean meat percentage and positively correlated with intramuscular fat [6]. The intramuscular fat (IMF) content plays a decisive role in the quality and flavor of pork [7]. Due to the positive genetic correlation between the two fat traits [8], increasing the intramuscular fat content while reducing the deposition of subcutaneous fat has always been a challenge in pig breeding.
The Dingyuan pig is an indigenous pig breed distributed in Dingyuan County, Anhui Province, China. It has the advantages of strong disease resistance, high intramuscular fat, and good meat quality, but also has a low lean meat percentage and relatively slow growth rate. The slaughter test of a large sample size of Dingyuan pigs revealed that some individuals had thinner subcutaneous fat even though their growth rate and IMF content were similar to those of other individuals. Therefore, we believe that this population may harbor key genes that specifically regulate subcutaneous and intramuscular fat deposition, providing a good study system to investigate the molecular mechanisms regulating tissue-level differences in fat deposition.
RNA sequencing (RNA-Seq) is a transcript quantification technology widely used to measure gene expression in tissues. Although RNA-Seq has been used to analyze the transcriptome profiles of subcutaneous and intramuscular fat in pigs, most of these studies have been conducted using different pig breeds [9, 10]. Some researchers have used pigs of the same breed with different phenotypes to study the mechanism of subcutaneous or intramuscular fat deposition [6, 11]; however, there are few reports on the simultaneous study of subcutaneous and intramuscular fat. In this study, we used RNA-Seq to study gene expression in subcutaneous and intramuscular fat simultaneously. We also used proteome analysis technology based on tandem mass tags (TMT) because it has improved proteome coverage and more reliable peptide identification and quantification, which is suitable for analyzing thousands of proteins in complex biological samples [12]. Such integrated transcriptome and proteome analyses can provide insights into key functional genes and their regulatory mechanisms between different types of fat deposits in pigs.
In this study, we collected tissues of subcutaneous fat and the longissimus dorsi (LD) muscle from pigs of similar age and body weight with divergent backfat thickness and intramuscular fat content. Transcriptomic and proteomic data were obtained using RNA-Seq and TMT to identify the key genes and pathways that specifically regulate subcutaneous fat and intramuscular fat deposition in Dingyuan pigs.
The phenotypic traits of the pigs used in this study are listed in Table 1. The backfat thickness and dressing percentage of the HBF group were significantly higher than those of the LBF group (P < 0.05). However, the body weight, intramuscular fat content, and lean meat percentage were not significantly different between the two groups (P > 0.05). The IMF content was significantly higher in the HIMF group than in the LIMF group (P < 0.05).
Groups |
HBF |
LBF |
HIMF |
LIMF |
---|---|---|---|---|
Weight (kg) |
88.17 ± 5.06a |
90.57 ± 3.55a |
86.17 ± 3.29a |
89.67 ± 2.04a |
IMF (%) |
6.88 ± 2.68a |
4.20 ± 2.05a |
6.30 ± 0.56a |
2.97 ± 0.95c |
Backfat (thickness/cm) |
5.53 ± 0.49a |
2.53 ± 0.42c |
5.10 ± 1.11a |
3.67 ± 1.12a |
Dressing percentage |
78.33 ± 0.01a |
72.33 ± 0.01c |
78.10 ± 2.13a |
74.53 ± 3.17a |
Lean meat percentage |
46.01 ± 5.03a |
50.94 ± 2.35a |
47.21 ± 3.33a |
50.42 ± 3.81a |
The total number of raw reads obtained by sequencing the 12 samples was 125.21 G and the number of clean reads was 120.87 G with an average proportion of 96.11% (94.28–97.59%; Table S2). The transcriptome data were compared to the reference genome (Sscrofa11.1) and the proportion of mapped reads was > 95% (95.37–96.66%), where an average of 85.8% reads were mapped within exons, 9.2% mapped to introns, and 5.0% mapped to intergenic regions (Fig. 1A).
After quantifying the expression of protein-coding genes using FPKM analysis, we compared the expression patterns of protein-coding genes in different samples. A total of 17,302 genes (FPKM > 0.05) were obtained in the LD muscle and back subcutaneous fat tissue of Dingyuan pigs, and 13081 genes were co-expressed in the four groups (Fig. 1B). The gene expression distributions of protein-coding genes were similar in all samples from different tissues (Fig. 1C).
A total of 191 genes were differentially expressed between the HBF and LBF groups (Table S3), of which 63 were upregulated and 128 were downregulated in HBF (Fig. 2A). The 191 DEGs were enriched in 83 significantly enriched GO terms, which mainly involved ion binding, metal ion binding, cation binding, cell differentiation, small molecule metabolic processes, and other related biological processes (Fig. 2B, Table S4A). KEGG enrichment analysis showed that the 21 significantly enriched pathways were primarily involved in metabolism, fatty acid degradation, fatty acid metabolism, notch signaling pathway, and hypertrophic cardiomyopathy, among others (Fig. 2C). Moreover, fatty acid biosynthesis and the AMP-activated protein kinase (AMPK), PPAR, MAPK, PI3K-Akt, adipocytokine, and FoxO signaling pathways—all related to lipogenesis—were also enriched (Table S4b).
A total of 85 DEGs were identified between the HIMF and LIMF groups, including 31 upregulated and 54 downregulated genes in HIMF (Fig. 3A, Table S5). The GO analysis identified 120 significantly enriched terms (Table S6a), of which the top 20 most enriched genes were mainly associated with catalytic activity, cell population proliferation, growth factor receptor binding, cytoskeleton organization, and cellular lipid metabolic process, among other functions (Fig. 3B). The results of KEGG pathway analysis showed that the DEGs were involved 81 pathways and 18 were significantly enriched (Fig. 3C, Table S6b). Among these, many were related to fat formation and metabolism, such as the FoxO, adipocytokine, PPAR, and AMPK signaling pathways.
To identify the genes that specifically regulate subcutaneous and intramuscular fat deposition, we compared the DEGs from the two tissues and found 3 overlapping DEGs—KRT80, PRKAG3, and SLC7A5. PRKAG3 is considered a candidate gene involved in subcutaneous and intramuscular fat deposition.
We selected 22 DEGs with specific differences in subcutaneous fat as candidate genes, including genes related to fatty acid biosynthesis and degradation (FASN, ACADSB, and GCDH), amino acid metabolism (ADHFE1, TKT, and ACAT1), propanoate metabolism (ECHDC1 and ACSS3), and carbon metabolism (ALDH1L1 and ALDH6A1), as well as 12 genes (AACS, SERPINE1, PPARD, UBD, UCP2, SERPINE1, FBP1, CA3, KLF2, PFKFB1, ACP5, PRG4) associated with fat cell differentiation, butanoate metabolism, the PPAR signaling pathway, regulation of gluconeogenesis, and glucose metabolism (Table 2).
Gene name |
log2FC in subcutaneous fat |
FDR in subcutaneous fat |
log2FC in the LD muscle |
FDR in the LD muscle |
Gene function |
---|---|---|---|---|---|
FASN |
-2.96285 |
0.000001 |
0.47826 |
0.67607 |
Metabolic pathways, fatty acid metabolism, AMPK signaling pathway, insulin signaling pathway, fatty acid biosynthesis |
CA3 |
-2.26031 |
0.004719 |
-0.84586 |
0.57883 |
Metabolic pathways, nitrogen metabolism, cellular anatomical entity, cellular process, metabolic process, binding, intracellular |
TKT |
-1.05613 |
0.000520 |
0.21110 |
0.97801 |
Metabolic pathways, carbon metabolism, amino acid biosynthesis, pentose phosphate pathway, alpha-amino acid metabolic process, sulfur compound metabolic process, positive regulation of potassium ion transport |
KLF2 |
1.41477 |
0.000055 |
0.16371 |
0.99102 |
White fat cell differentiation |
ACAT1 |
-1.01960 |
0.000071 |
-0.43082 |
0.95335 |
Reproductive structure development, sulfur amino acid biosynthetic process, metabolic pathways, valine, leucine, and isoleucine degradation, propanoate metabolism, carbon metabolism, pyruvate metabolism |
ALDH1L1 |
-2.35834 |
8.61E-08 |
-0.90571 |
0.91253 |
One carbon pool by folate, One Carbon Metabolism, Folate Metabolism |
ECHDC1 |
-1.43454 |
5.01E-09 |
-0.02139 |
0.99586 |
Metabolic pathways, propanoate metabolism, metabolic process, catalytic activity |
ADHFE1 |
-1.85280 |
2.19E-08 |
-0.50819 |
0.92831 |
Reproduction, alpha-amino acid metabolic process, sulfur compound metabolic process, reproductive structure development, regulation of muscle system process, citric acid cycle, respiratory electron transport, pyruvate metabolism |
ACADSB |
-1.18737 |
0.000048 |
-0.29033 |
0.95987 |
Metabolic pathways, valine, leucine, and isoleucine degradation, fatty acid degradation, fatty acid metabolism |
ALDH6A1 |
-1.17949 |
0.004352 |
-0.79455 |
0.54191 |
Metabolic pathways, valine, leucine, and isoleucine degradation, propanoate metabolism, carbon metabolism, inositol phosphate metabolism, beta-alanine metabolism |
AACS |
-2.27971 |
2.28E-10 |
-0.68441 |
0.73266 |
Metabolic pathways, valine, leucine, and isoleucine degradation, butanoate metabolism |
SERPINE1 |
1.63599 |
0.003084 |
-0.99366 |
0.77249 |
Adipogenesis, blood clotting cascade, complement and coagulation cascades |
PPARD |
1.30507 |
0.000548 |
-1.19111 |
0.23847 |
Pathways in cancer, PPAR signaling pathway, acute myeloid leukemia, Wnt signaling pathway, ion binding, metal ion binding, small molecule metabolic process |
UBD |
1.02161 |
0.001181 |
0.07256 |
0.99490 |
Proteasome binding, protein ubiquitination, positive regulation of apoptotic process |
UCP2 |
1.44530 |
0.002745 |
0.27950 |
0.97436 |
Anatomical structure morphogenesis, reproductive structure development, cardiac muscle tissue development, cell junction organization, muscle cell development |
FBP1 |
-1.78347 |
0.004800 |
-3.33968 |
/ |
Negative regulation of glycolytic process, regulation of gluconeogenesis, fructose 6-phosphate metabolic process, negative regulation of cell growth |
ACSS3 |
-1.19635 |
0.009090 |
-0.73653 |
0.62802 |
Metabolic pathways, propanoate metabolism |
PFKFB1 |
-3.26636 |
6.02E-06 |
0.90415 |
0.26914 |
Regulation of glycolysis by fructose 2,6-bisphosphate metabolism, metabolism, glycolysis, glucose metabolism, focal adhesion-PI3K-Akt-mTOR-signaling pathway, carbohydrate metabolism |
GCDH |
-1.16760 |
2.19E-08 |
-0.19610 |
0.98930 |
Metabolic pathways, fatty acid degradation, lysine degradation, tryptophan metabolism |
ACP5 |
2.07205 |
0.007519 |
0.92431 |
0.95335 |
Metabolism, metabolism of water-soluble vitamins and cofactors, NAD phosphorylation and dephosphorylation |
PRG4 |
3.61705 |
0.525778 |
-1.38560 |
1.94832 |
Phospholipase-C Pathway, ERK signaling, integrin pathway, MAPK signaling |
LPL |
1.00349 |
0.002139 |
1.15750 |
0.10242 |
Fatty acid β-oxidation, adipogenesis, PPAR signaling pathway, lipoprotein metabolism, triacylglyceride synthesis |
PRKAG3 |
-2.73175 |
0.007801 |
1.87501 |
0.00149 |
Longevity regulating pathway, AMPK signaling pathway, apelin signaling pathway, insulin signaling pathway, oxytocin signaling pathway, non-alcoholic fatty liver disease, tight junction |
RETREG1 |
-0.73548 |
0.403354 |
-2.72111 |
0.00000 |
nucleolus, endoplasmic reticulum, Golgi apparatus |
PRKAG2 |
-0.50079 |
0.068174 |
-2.57251 |
0.00403 |
Vitamin digestion and absorption, thermogenesis, lipid metabolism |
SMPDL3A |
0.20620 |
0.868544 |
-3.52920 |
0.00000 |
Vitamin digestion and absorption, sphingomyelin phosphodiesterase activity, sphingomyelin metabolic process, cellular lipid metabolic process, membrane lipid catabolic process, sphingolipid catabolic process, phospholipid catabolic process |
IRS2 |
-0.03815 |
0.971924 |
-1.78695 |
0.00030 |
Adipogenesis genes, focal adhesion-PI3K-Akt-mTOR-signaling pathway, erythropoietin activates phosphoinositide-3-kinase, IL-13 signaling pathway, signaling by type 1 insulin-like growth factor 1 receptor |
BDH1 |
-0.05789 |
0.975117 |
-1.92965 |
0.00878 |
Vitamin digestion and absorption, butanoate metabolism, metabolic pathways, ketone body catabolism, lipid metabolism |
PPARA |
-0.58082 |
0.502387 |
-1.59100 |
0.00067 |
Vitamin digestion and absorption, cAMP signaling pathway, Hepatitis C |
GK |
0.06790 |
0.966532 |
-1.98564 |
0.00415 |
Vitamin digestion and absorption, glycerolipid metabolism, metabolic pathways |
LEP |
0.71676 |
0.328614 |
3.16540 |
0.00023 |
AMP-activated protein kinase Signaling, adipocytokine signaling pathway, peptide hormone metabolism, cytokine-cytokine receptor interaction |
In the LD muscle, following GO and KEGG analysis, we selected 8 genes mainly involved in fat metabolism-related pathways as candidate genes for the specific regulation of intramuscular fat deposition. These included the fat synthesis gene IRS2 (insulin receptor substrate 2), the transcription factor PPARA (peroxisome proliferator activated receptor alpha), the adipocyte secretory product LEP (leptin), genes related to fatty acid β-oxidation, the lipid metabolism-related gene RETREG1 (reticulophagy regulator 1), GK (glycerol kinase), BDH1 (3-hydroxybutyrate dehydrogenase 1), SMPDL3A (sphingomyelin phosphodiesterase acid like 3A), and PRKAG2 (protein kinase AMP-activated non-catalytic subunit gamma 2). The expression levels of these candidate genes in the two tissues are shown in Fig. 4.
A total of 29,519 peptides were obtained from the MS analysis of the subcutaneous fat tissue, and 4,099 proteins were identified. In contrast, 14,619 peptides and 2,178 proteins were obtained from the LD muscle. Among the identified proteins, 78.3% of those in the subcutaneous fat tissue (Fig. S1A) and 82.3% of those in the LD muscle (Fig. S1B) were represented by 1–10 peptides, and the molecular weight of the proteins mainly ranged from 10 to 80 kDa (Fig. S1C, Fig. S1D). The coefficient of variation of > 90% of the four groups of replicates was less than 30%, indicating that our experimental samples had good biological reproducibility (Fig. S2).
In total, 62 DAPs were identified between the HBF and LBF groups, of which 22 were upregulated and 40 were downregulated (Fig. S3, Table S7a). The heatmap of the DAPs showed that the expression levels had good repeatability within each group (Fig. S4). The 62 DAPs were enriched in 69 GO terms (Table S7b). Approximately half of the 20 most significantly enriched GO terms were related to lipid and fatty acid metabolism, including the fatty acid catabolic process, fatty acid oxidation, fatty acid metabolic process, lipid oxidation, lipid metabolic process, cellular lipid metabolic process, and lipid catabolic process, among others (Fig. 5A, Table S7b). The enriched pathways indicated by KEGG analysis mainly included fatty acid degradation, fatty acid metabolism, and the PPAR signaling pathway and metabolic pathways (Fig. 5B, Table S7c). The 8 DAPs (ACAA2, ACAT1, ACOX1, CPT1A, ACSL4, SDHD, and IDH3A) were mainly associated with these pathways. In addition, 4 genes—MMUT, PCCB, HMGCL, and ALDH6A1—were involved in the synthesis and degradation of ketone bodies and propanoate metabolism, and are also considered important candidate proteins (Table 3).
Protein name |
Fold change (HBF/LBF) in the proteome |
P-value |
Fold change (HIMF/LIMF) in the proteome |
P-value |
Functional analysis |
---|---|---|---|---|---|
ACAA2 |
0.81683 |
0.02263 |
0.84136 |
0.36089 |
Propanoate metabolism, fatty acid degradation, fatty acid metabolism |
ACAT1 |
0.70871 |
0.04268 |
0.95727 |
0.87357 |
Propanoate metabolism, fatty acid degradation, fatty acid metabolism, synthesis and degradation of ketone bodies |
ALDH6A1 |
0.81354 |
0.00886 |
0.81433 |
0.08507 |
Valine, leucine, and isoleucine degradation, propanoate metabolism, carbon metabolism, beta-alanine metabolism, metabolic pathways |
ACOX1 |
0.80644 |
0.02469 |
/ |
/ |
Propanoate metabolism, fatty acid degradation, fatty acid metabolism, PPAR signaling pathway |
PCCB |
0.78891 |
0.04731 |
0.92308 |
0.41421 |
Propanoate metabolism |
MMUT |
0.73611 |
0.01595 |
0.93050 |
0.34846 |
Propanoate metabolism |
IDH3A |
0.83206 |
0.00272 |
0.97174 |
0.79603 |
Carbon metabolism, metabolic pathways, citric acid cycle |
HMGCL |
0.78701 |
0.04835 |
/ |
/ |
Valine, leucine, and isoleucine degradation, synthesis and degradation of ketone bodies, peroxisome, metabolic pathways, butanoate metabolism |
CPT1A |
1.37812 |
0.03484 |
/ |
/ |
Fatty acid degradation, fatty acid metabolism, PPAR signaling pathway, adipocytokine signaling pathway, thermogenesis |
SDHD |
0.77096 |
0.00692 |
0 |
0 |
Carbon metabolism, metabolic pathways, citric acid cycle, thermogenesis, Alzheimer’s disease |
ACSL4 |
1.20345 |
0.00644 |
/ |
/ |
Fatty acid degradation, fatty acid metabolism, peroxisome, metabolic pathways, PPAR signaling pathway, adipocytokine signaling pathway |
Twelve DAPs were found in the HIMF vs. LIMF groups, of which 4 were upregulated and 8 were downregulated (Table S8a). The biological processes associated with the DAPs included developmental processes, multicellular organismal processes, protein modification processes, cellular processes, biological regulation, metabolic processes, and other biological processes (Table S8b). In total, 14 signal pathways enriched by KEGG analysis (Table S8c), including the PI3K-Akt signaling pathway, regulation of the actin cytoskeleton, and metabolic pathways, among others. TRIM55 (tripartite motif containing 55; involved in cellular protein metabolism), PTGR2 (prostaglandin reductase 2; involved in prostaglandin metabolism), and UBD (ubiquitin D; involved in cellular protein metabolism) were considered to be a candidate gene for IMF deposition.
On integrating the 191 DEGs and 62 DAPs from comparisons of the HBF and LBF groups, we found 8 overlapping genes (ACAT1, ALDH6A1, ISLR, HSDL2, MCCC2, IVD, EPHX1, and PRG4). Moreover, all 8 genes showed the same expression trends in terms of mRNA and protein. The PRG4 gene was highly expressed in the HBF group, whereas the remaining 7 genes were low in HBF. Based on the DEGs and DAPs, a protein-protein interaction (PPI) network was established for the Sus scrofa database using the STRING v.10.0 online software (Fig. S5). We found that IVD, ALDH6A1, ACSS3, ECDH1, ACADSB, ACAA2, ACSS1, HSDL2, HMGCL, PCCB, MCCC2, PCCA, AACS, ACOX1, and ACAT1 played pivotal roles in the network.
For the LD muscle tissue, among the 85 DEGs and 12 DAPs, only one overlapping gene—LMOD2 (leiomodin 2)—was found. This gene is mainly involved in the processes of myofibril assembly and muscle structure development. The PPI network diagram shows that the DEGs and DAPs in the LD muscle group were mainly centered on the INS gene (Fig. S6) and interacted with actin-related proteins and regulators such as ACTR10, ACTR1A, PFN2, and ENAH, whereas those in the other group interacted with insulin- and insulin receptor-related genes and substrate genes such as INSR, IRS1, and fat regulatory genes LEP and PPARA.
We selected ten DEGs (ALDH1L1, FASN, LPL, IGFBP5, ACAT1, INSR, RETREG1, FZD10, PPARA, PRKAG3) to test the validity of the RNA-Seq and TMT-based proteomic results. For this, we performed RT-qPCR using RNA samples from fat and muscle tissues of Dingyuan pigs. The expression levels of these genes in each group are shown in Fig. S1. The 10 genes selected were differentially expressed among the groups, and the mRNA and protein expression trends of these 10 genes were concordant with those obtained by RT-qPCR. Further, fold changes of all 10 genes in the qRT-PCR and RNA-Seq showed the same trends (Fig. S2). The results indicated that the DEGs identified by RNA-Seq and the DAPs identified by TMT-based proteomics were reliable and efficient.
Dingyuan pigs have high subcutaneous and intramuscular fat contents. The population is a native breed that has not undergone strong artificial selection, and contains individuals with divergent backfat thickness and intramuscular fat content. Comparing gene expression between individuals with divergent traits in the same population can reduce noise due to different genetic backgrounds. The clean reads obtained from the transcriptome mapped to > 94.2% of the genome, which is higher than the value reported in previous studies on the pig backfat and LD muscle transcriptome [9, 13–15]. The current results show that our transcriptome sequencing results are of high quality. Through integrated transcriptomic and proteomic analysis, we found 8 overlapping genes in the subcutaneous fat group and only 1 overlapping gene in the LD muscle group. Due to the existence of multiple post-transcriptional regulatory mechanisms, similar results have been reported in many studies [16–19]. This prompted us to analyze integrated transcriptomic and proteomic data to better identify functional genes and their regulatory pathways.
Adipose tissue is not only an important energy store but is also essential for normal carbohydrate and lipid homeostasis. In addition to storing energy, white adipose tissue can also secrete enzymes and transcription factors, and are important secretory organs in the body [20, 21]. On comparing the KEGG pathways enriched by DEGs in the subcutaneous fat and the LD muscle, we found that the hypertrophic cardiomyopathy signaling pathway was the only one common to both groups. Most of the DEG-enriched pathways identified in subcutaneous fat were metabolic pathways, including the metabolic processes of various amino acids, esters, fatty acids, and other intermediate products in glucose metabolism. The DEGs in the LD muscle group participated in a wide range of signaling pathways, including insulin, ketone body synthesis and degradation, longevity, and some disease-related pathways. Our results are consistent with those of another study which reported that pathways related to adipogenesis, lipolysis, glycolysis, and fatty acid oxidation were significantly downregulated in intramuscular adipocytes compared those in adipocytes from other locations [22, 23]. Previous studies have also reported the differences in the regulatory mechanisms of subcutaneous and intramuscular fat production [24].
We separately analyzed the signaling pathways involved in the genes upregulated and downregulated in the HBF group and found that the highly expressed DEGs in the LBF group were mainly involved in fatty acid degradation and metabolism, pyruvate metabolism, pentose phosphate pathway, and other pathways related to lipid metabolism. The DEGs more highly expressed in the HBF group had functional categories related to inflammation and immunity, such as viral myocarditis, lipid and atherosclerosis, the AGE-RAGE signaling pathway in diabetic complications, and natural killer cell-mediated cytotoxicity. Human obesity is known to be associated with many diseases, including diabetes, dyslipidemia, cardiovascular disease, and cancer [25]. As an animal model for studying human obesity, pigs exhibit obesity complications similar to those of humans [26]. Therefore, it is particularly important to identify the different mechanisms of regulation of intramuscular and subcutaneous fat such that IMF content can be increased without increasing backfat deposition.
Candidate genes affecting the production of subcutaneous adipose tissue and their regulatory changes
Among the 8 genes identified in the subcutaneous fat group, MCCC2 and IVD are involved in leucine metabolism. Based on their functions, the genes ACAT1, ALDH6A1, and PRG4 were selected as candidate proteins involved in regulating subcutaneous fat deposition in pigs. In addition to these common genes, we also identified proteins involved in fatty acid metabolism (ACAA2, ACOX1, CPT1A, and ACSL4), propionate metabolism (ALDH6A1, PCCB, MMUT), ketone body metabolism (HMGCL), and the citric acid cycle (IDH3A); the genes coding these proteins may be candidate genes for regulating subcutaneous fat in pigs. Our transcriptome proteome had very few overlapping genes, and the results indicate that there is post-transcriptional regulation of fat deposition in Dingyuan pigs. In addition, due to the limitations of sequencing technology, the number of proteins identified by the proteome was far lower than that in the transcriptome. Moreover, the number of proteins detected in the two tissues was very different, and many proteins were not detected. Therefore, we believe that DEGs that exist only in the transcriptome can also be used as important candidate genes.
In our study, PRKAG3 was the only DEG found to participate in fat metabolism in the transcriptome of both tissues. The PRKAG3 gene encodes the AMPK regulatory subunit γ3 [27]. AMPK activity is associated with increased glucose uptake and fatty acid oxidation, as well as inhibition of glycogen synthase activity and fatty acid synthesis [28]. Two missense mutations in the pig PRKAG3 gene—A595G (Ile199Val) and G154A (Gly52Ser)—are related to the IMF content, water-holding, and meat quality of pigs [29], and the AMPKγ3R200Q mutation results in a lack of AMP dependence and elevated basal activity of AMPK [28]. The current research on pig PRKAG3 is focused on IMF. Based on the FPKM value of the transcriptome, we found that PRKAG3 expression is high in muscle tissue and very low in fat tissue; therefore, the role of this gene in pig backfat deposition remains to be verified.
Our results showed that in the HBF and HIMF groups, most of the candidate genes that regulate the biological process of adipogenesis were downregulated. Lipoprotein lipase (LPL) is one of the key enzymes involved in lipid metabolism. It is the rate-limiting enzyme for the catabolism of triglycerides and can remove TG-rich, very low-density lipoproteins and chylomicrons from the blood [30, 31]. Fatty acid synthase (FASN) is a key enzyme in the de novo fatty acid synthesis pathway and cell substrate energy metabolism [32]. Acetoacetyl-CoA synthetase (AACS) is an enzyme that provides acetyl units for the biosynthesis of cholesterol and fatty acids [33]. AACS gene expression in white adipose tissue is lower in Zucker fatty rats than in lean rats; however, in high-fat diet induced obese rats, the expression of this gene is increased [34]. ALDH1L1 and ALDH6A1 belong to the aldehyde dehydrogenases family. The ALDH6A1 gene is involved in the metabolism of a variety of aliphatic and aromatic aldehydes produced by various endogenous and exogenous precursors, and is also an important candidate gene that affects drip loss [35]. Moreover, ALDH1L1 is an enzyme that metabolizes folate [36]. At present, there are almost no reports on the relationship between these two genes and fat deposition. However, based on the various metabolic pathways they participate in and their differential expression in the HBF and LBF groups, we speculate that they may play a role in the subcutaneous fat deposition process in pigs.
ADHFE1 expression is highly enriched in adipose tissue and other highly metabolically active tissues, and its expression is closely related to the phenotype of mature adipocytes in vivo and in vitro [37]. The candidate genes related to subcutaneous adipogenesis that were upregulated in the HBF group included KLF2, ACP5, SERPINE1, PPARD, UBD, and UCP2. KLF2 inhibits fat formation by negatively regulating adipocyte differentiation, but does not affect the commitment of multipotent stem cells to the preadipocytic lineage [38]. ACP5, which is associated with excessive fat deposition and the development of atherosclerosis in pigs [39]. SERPINE1 (also known as PAI-1) is mainly produced in white adipose tissue and increases the cardiovascular risk in obese and diabetic patients, whereas its deficiency can prevent obesity and metabolic dysfunction [40, 41]. The peroxisome proliferator-activated receptor delta gene (PPARD) is a key regulator of lipid metabolism, and its haplotype is associated with pig backfat thickness [42]. Studies in C2C12 myoblasts have shown that overexpression of PPARD inhibits myotube formation and enhances adipocyte differentiation [43]. In mice, knockout of the UBD gene prevents the development of age-related obesity while prolonging lifespan and vitality [44]. UCP2 participates in the oxidation of mitochondrial substrates and plays an important role in metabolic pathways [45, 46]. Studies have reported that some SNPs in UCP2 are related to obesity and type 2 diabetes [47–49]. Among these differential genes, genes that are positively related to adipogenesis were down-regulated (e.g., FASN), and some were upregulated (e.g., ACP5), indicating that fat metabolism in the body is a dynamic and balanced process.
Candidate genes affecting the production of IMF and their regulatory changes
In the transcriptome of the IMF group, we found many enriched pathways related to lipid metabolism. Based on these pathways, we identified genes that may regulate intramuscular fat deposition. RETREG1, also known as pFAM134B, promotes the accumulation of subcutaneous fat in pigs by increasing the mRNA levels of PPARγ, FAS, and ACC, and reducing the levels of ATGL and HSL [50]. PRKAG2 is an important regulator of cardiac metabolism [51]; an SNP of PRKAG2 is associated with insulin sensitivity and is essential for obesity, triglycerides, and HDL cholesterol [52]. SMPDL3A—a recently identified phosphodiesterase—plays a role in cholesterol and fat metabolism by regulating the liver X-receptor [53]. The insulin receptor substrate 2 (IRS2) is the first downstream effector of the insulin receptor. Studies have found that the absence of IRS2 in myeloid cells improves glucose homeostasis and enhances resistance to metabolic dysfunction induced by a high-fat diet [54]. In addition, the proteins LEP [55], GK [56], and PPARA [57]—reported to be associated with obesity or adipose-related biological processes—were also included in our results.
The purpose of our experiment was to identify genes that specifically regulate subcutaneous and intramuscular fat in pig with extreme backfat and IMF. Interestingly, some of the DEGs and DAPs identified in the backfat groups in this study—namely, IGFBP5, PTGR2, PPARD, GK, LEP, PPARA, and UBD—have also been reported as candidate genes regulating intramuscular fat in other studies. For instance, UBD expression was upregulated in our HBF transcriptome and in the proteome of the HIMF group. Studies have shown that UBD has a positive regulatory effect on intramuscular and subcutaneous fat, and that inhibition of UBD in intramuscular and subcutaneous adipocytes results in the inhibition of cell proliferation and lipid droplet production [58].
From our results, it is worth mentioning that the RETREG1 gene was downregulated in the HIMF group. Although it has been suggested that this gene has an inhibitory effect on intramuscular fat deposition in Dingyuan pigs, studies have shown that it also promotes subcutaneous fat deposition. The opposing roles of genes in the two types of adipose tissue are beneficial for the selection of high-quality pigs, and their functions can be explored through further studies.
Excluding the genes that co-regulate intramuscular and subcutaneous fat, we suggest that in Dingyuan pigs, the FASN, LPL, AACS, TKT, ALDH6A1, ALDH1L1, KLF2, ACP5, SERPINE1, UCP2, FBP1, ACSS3, PFKFB1, ECHDC1, ACADSB, GCDH, ACAA2, ACOX1, CPT1A, ACSL4, ALDH6A1, PCCB, MMUT, HMGCL, IDH3A, PRG4, and CTSL genes may regulate subcutaneous fat deposition, whereas PRKAG2, SMPDL3A, IRS2, BDH1, TRIM55, ACSL4, and CTP1A may regulate intramuscular fat deposition.
The Dingyuan pigs were raised in the Ankang Agriculture and Animal Husbandry Company, Dingyuan County, Anhui Province, with consistent feeding conditions. Over 200 unmated sows that were 10 months old and showed normal growth were selected and their body weight and live backfat thickness were measured using a platform scale and a B-ultrasound instrument (SSD-500V, ALOKA, Japan). Based on the measurements, 12 individuals with similar body weights were selected such that six had extremely high live backfat thickness and six had extremely low live backfat thickness. The pigs were slaughtered to measure the carcass weight, backfat thickness, lean meat percentage, and meat quality traits. The LD muscle and back subcutaneous fat tissues were collected, immediately frozen in liquid nitrogen, and stored at -80°C. Animals were reared and handled in accordance with the Guide for the Care and Use of Laboratory Animals in China. All experimental guidelines were approved by the Committee on the Ethics of Animal Experiments of China Agricultural University (permit number: SKLAB-2012-04-07).
A fresh sample of the LD muscle tissue was cut into pieces, the connective tissue and fascia were removed, and the remaining tissue was homogenized. The LD muscle samples were then dried to a constant weight in an oven and the dried samples were ground into powder. The Soxhlet extraction method was used to determine IMF content in muscle tissue.
Based on backfat thickness and IMF content, subcutaneous fat samples were divided into high backfat thickness (HBF, n = 3) and low backfat thickness (LBF, n = 3), and LD muscle samples were divided into the high intramuscular fat (HIMF, n = 3) and low intramuscular fat (LIMF, n = 3) groups.
Total RNA was extracted from the LD muscle and back subcutaneous fat tissues using the TRIzol method. RNase-free DNase was added to total RNA to remove DNA contamination. The integrity of RNA was checked by 1.2% agarose gel electrophoresis and the quality and amount of total RNA were determined using a Nanodrop 2000 instrument (Thermo Fisher Scientific, Waltham, MA, USA). The RNA integrity number (RIN) values of all samples were above 7, and 2.5 µg of total RNA was used to construct the cDNA library. The Illumina NovaSeq6000 platform (Illumina, San Diego, CA, USA) was used for transcriptome sequencing, and the sequencing read length was paired-end 150 bp. From the raw data obtained by sequencing, we filtered out low-quality sequences (Q-value < 19), sequences contaminated by linkers, and reads with an N content > 5% to obtain clean data.
Clean reads were mapped to the reference pig (Sus scrofa) genome version scrofa.Sscrofa11.1 (ftp://ftp.ensembl.org/pub/release-99/fasta/sus_scrofa/dna/Sus_scrofa.Sscrofa11.1.dna.toplevel.fa) using the HISAT2 (V2.0.5) software for alignment. We used the Cufflinks (V2.2.1) software to calculate fragments per kilobase million mapped reads (FPKM) values to represent gene expression levels. Genes with an average FPKM value > 0.05 (including all individuals in one group) were defined as being expressed for that group, and those with an average FPKM value < 0.05 were defined as being non-expressed in the group. Here, we used DESeq2 (V1.25.9) to analyze the differences between the four groups of Dingyuan pig samples. Differentially expressed genes (DEGs) were defined as those with false discovery rate (FDR) < 0.01 and absolute log2 of the fold change (|log2FC|) ≥ 1.
The samples used for the proteomic analyses were the same as those used for RNA-SEq. An appropriate amount of sample tissue was added to the protein lysis solution (Mammal Tissue Total Protein Extraction Kit AP0601-50, Invent Biotechnologies, Inc., Plymouth, MN, USA). After sonication, the samples were placed on ice for 20 min, centrifuged at 10,000 × g at 4°C for 30 min, transferred to a new tube, and stored at -80°C. Protein concentrations of the samples were determined using a BCA kit (Thermo Fisher Scientific), and protein quality was detected by SDS-PAGE. The peptides were labeled with TMT isobaric tags (Thermo Fisher Scientific) at room temperature for 1 h.
The TMT-labeled samples were fractionated using a high-pH peptide fractionation protocol. Each sample was fractionated into 60 fractions using a step gradient of 0–95% acetonitrile. Next, the 60 fractions were pooled into 10 fractions in a non-contiguous manner. Each sample was separated using a high-performance liquid chromatography (HPLC) system. The chromatographic column was balanced with 95% solution A (98% H2O, pH10), loaded onto the mass spectrometer pre-column using the autosampler, and then separated by an analytical column. The separated samples were analyzed using a Q Exactive Plus mass spectrometer (Thermo Fisher Scientific).
After mass spectrometry analysis, the Proteome Discoverer 2.4 software (Thermo Fisher Scientific) was used to process all the original files against the Sus_scrofa.Sscrofa11.1.pep.all.fa database. The peptide confidence was set to a high level (q-value < 0.01) for peptide filtering. Quantification of experimental bias was set to normalize the total peptide amount. FC ≥ 1.2, FC ≤ 0.833, and P-value ≤ 0.05 were set as thresholds for identifying differentially abundant proteins (DAPs).
Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the DEGs was performed using the KOBAS 3.0 enrichment analysis tool (https://kobas.cbi.pku.edu.cn/anno_iden.php). Functional analysis of DAPs was conducted using the UniProt ftp database and the KEGG PATHWAY database. The Pathway Maps tool was used to enrich the pathways, and P-values were calculated using R software based on a hypergeometric distribution, with the default database being used as the background.
To confirm the DEGs identified in RNA-Seq and the DAPs identified in TMT-based proteomics, qRT-PCR was performed on 10 genes to measure their expression levels in different groups. The primer sequences used are listed in Table S1. Quantitative real-time PCR (qPCR) reactions were performed using the SYBR Green qPCR SuperMix (Transgen, Beijing, China) using a Bio-Rad CFX96 System (Bio-Rad, Hercules, CA, USA) with a reaction volume of 20 µL. A cDNA pool of all samples was used for calibration and the experiment was performed in triplicate for each sample. Gene expression levels were calculated using the 2-ΔΔCt method.
The differences in slaughter traits and gene expression levels were analyzed using the t-test in SPSS (version 21.0; IBM Corp. Released 2012. IBM SPSS Statistics for Windows, Armonk, NY, USA). Graphs were prepared using GraphPad Prism (version 7), and data are presented as the mean ± standard error.
DAP, differentially abundant protein; DEG, differentially expressed gene; FDR, false discovery rate; FPKM, fragments per kilobase million mapped reads; GO, gene ontology; HBF, high backfat thickness group; HIMF, high intramuscular fat group; IMF, intramuscular fat; KEGG, Kyoto Encyclopedia of Genes and Genomes; LBF, low backfat thickness group; LD, longissimus dorsi; LIMF, low intramuscular fat group; log2FC, log2 fold change; qPCR, quantitative real-time PCR; TMT, tandem mass tags.
Ethics approval and consent to participate
All animal work was conducted following the guidelines for the care and use of experimental animals. The Animal Welfare Committee of the State Key Laboratory for Agro-Biotechnology of the China Agricultural University approved all procedures for animal care (approval number, SKLAB‐2012‐04‐07). The study was carried out in compliance with the ARRIVE guidelines.
Consent for publication
Not applicable.
Availability of data and material
All RNA-Seq data were deposited in the Gene Expression Omnibus under accession number GSE179457 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE179457). The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD027061.The datasets supporting the conclusions of this article are included within the article and its additional files.
Competing interests
The authors declare no conflict of interest.
Funding
This work was supported by the National Natural Science Foundation of China (32060736), the Yunnan Science and Technology Project (No. 22018FG001-003), and Central Government Guides Local Science and Technology (No. YDZX20195400004426).
Authors' contributions
P.Z., and H.Z. designed the research; P.Z., Y.Z., H.Z., Y.W., and Q.L. collected samples, extracted RNAs, and performed qRT-PCR verification experiments; P.Z. performed most of the bioinformatics studies; P.Z., H.Z., and B.Z. wrote and revised the manuscript; and all authors reviewed the manuscript. The authors read and approved the final manuscript.
Acknowledgements
We thanks Annoroad (Beijing, China) and Beijing Bangfei Bioscience &Technology Co., Ltd (Beijing, China) for the sequencing services and bioinformatic analysis.