Identification of Candidate Genes That Specifically Regulate Subcutaneous and Intramuscular Fat Deposition Using Transcriptomic and Proteomic Profiles in Dingyuan Pigs

Abstract

Subcutaneous fat and intramuscular fat (IMF) deposition are closely related to meat production and pork quality. Dingyuan pig, as a native pig breed in China, low selection leads to obvious genetic and phenotypic differences in the population. Individuals with extreme fat content in the population are ideal models for studying the mechanism of fat deposition. In this study, we used RNA-Seq and tandem mass tags-based (TMT) proteomics to analyze the key pathways and genes that specifically regulate subcutaneous fat and IMF deposition in Dingyuan pigs. We identified 191 differentially expressed genes (DEGs) and 61 differentially abundant proteins (DAPs) in the high and low back fat thickness (HBF, LBF) groups, 85 DEGs and 12 DAPs were obtained in the high and low intramuscular fat (HIMF, LIMF) groups. The functional analysis showed that the DEGs and DAPs in the backfat groups were mainly involved in carbohydrates, amino acids, and fatty acids, whereas the IMF groups were involved in the insulin pathway, longevity, and some disease-related pathways. we found 33 candidate genes that might tissue-specifically lipids deposition for subcutaneous and intramuscular fat. Our research provides theoretical reference materials for the improvement of fat deposition traits of local pig breeds in my country.

Introduction

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 ¹. Commercial pig breeds increase the proportion of lean meat and reduce the deposition of subcutaneous fat, while the reduced intramuscular fat content at the same time results in poor flavor ^2–5. Porcine subcutaneous fat is quantified by backfat thickness, which is negatively correlated with lean meat percentage and positively correlated with intramuscular fat ⁶. The intramuscular fat (IMF) content plays a decisive role in the quality and flavor of pork ⁷. Due to the positive genetic correlation between the two fat traits ⁸, 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 high intramuscular fat, and good meat quality, but also has a low lean meat percentage and relatively slow growth rate. A large population of Dingyuan pigs was phenotyped using multi-production trait indexes, and found that the growth rate and IMF content of some individuals were similar to those of others, but their subcutaneous fat was thinner. 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-Seq has been widely used to analyze the transcriptome profile of pig subcutaneous and intramuscular fat. Jinhua pig and Landrace pig intramuscular fat transcriptome sequencing revealed the differences in fat and protein metabolism at the transcriptional level of the two pigs ⁹. Sequencing of subcutaneous fat and IMF of 100 pigs identified 215 subcutaneous fat genes and 90 IMF genes ¹⁰. The mechanism of subcutaneous fat and IMF deposition has not only been studied in different breeds of pigs, but also in Italian Great White pigs or Duroc pigs ¹¹ of the same breed with different phenotypes ⁶. Dingyuan pig is an excellent indigenous pig breed, but there is no relevant report yet.

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. 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 ¹². 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.

Materials And Methods

Animal ethics compliance

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.

Animals and samples

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 and lean meat percentage. The LD muscle and back subcutaneous fat tissues were collected, immediately frozen in liquid nitrogen, and stored at -80°C.

Measurement of intramuscular fat

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.

Library preparation and sequencing

Total RNA was extracted from the LD muscle and back subcutaneous fat tissues using the TRIzol (Invitrogen, Carlsbad, CA, USA) method. Isolated total RNA was quantified (Nanodrop, ND2000) and quality controlled with typical curves (Agilent, Bioanalyzer 2100). 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 (TruSeq RNA Sample Preparation Guide, Illumina Inc., San Diego, CA). The Illumina NovaSeq6000 platform (Illumina, San Diego, CA, USA) was used for transcriptome sequencing, and the sequencing read length was paired-end 150 bp. The obtained raw data were filtered to clean data with fastp (version 0.12.3) by removing reads containing adapters, low-quality reads (Q-value <19), and reads with an N content > 5%. Clean reads were mapped to the reference pig (Sus scrofa) genome version scrofa.Sscrofa11.1.99 ⁵⁹(http://ftp.ensembl.org/pub/release-99/fasta/sus_scrofa/dna/Sus_scrofa.Sscrofa11.1.dna.toplevel.fa.gz) 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. The DESeq2 (V1.25.9) was used to analyze the differences between groups, and the P-values were adjusted using the Benjamini-Hochberg method. Differentially expressed genes (DEGs) were defined as those with false discovery rate (FDR) < 0.01 and |log₂ fold change| ≥ 1.

Protein sample preparation and TMT tag labeling

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

Functional annotation of DEGs and 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.

Validation of important candidate genes

To confirm the DEGs identified in RNA-Seq and the DAPs identified in TMT-based proteomics, qRT-PCR was performed on 10 genes with the same samples as those used for sequencing to measure their expression levels in different groups. The primer sequences used are listed in Supplementary Table S8. 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.

Results

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 (Supplementary Fig. S1A) and 82.3% of those in the LD muscle (Supplementary Fig. S1B) were represented by 1–10 peptides, and the molecular weight of the proteins mainly ranged from 10 to 80 kDa (Supplementary Fig. S1C, Supplementary 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 (Supplementary Fig. S2).

Screening And Functional Classification Of Daps

In total, 62 DAPs were identified between the HBF and LBF groups, of which 22 were upregulated and 40 were downregulated (Supplementary Fig. S3, Supplementary Table S6a). The heatmap of the DAPs showed that the expression levels had good repeatability within each group (Supplementary Fig. S4). The 62 DAPs were enriched in 69 GO terms. 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, Supplementary Table S6b). 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, Supplementary Table S6c). 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).

Twelve DAPs were found in the HIMF vs. LIMF groups, of which 4 were upregulated and 8 were downregulated (Supplementary Table S7a). 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 (Supplementary Table S7b). In total, 14 signal pathways enriched by KEGG analysis (Supplementary Table S7c), 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.

Integrated Analysis Of Transcriptomic And Proteomic Data

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 (Supplementary Fig. S4). 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 (Supplementary Fig. S5) 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.

Verification Of Degs And Daps

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 Supplementary Fig. S6. 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 qPCR and RNA-Seq showed the same trends (Fig. 6). The results indicated that the DEGs identified by RNA-Seq and the DAPs identified by TMT-based proteomics were reliable and efficient.

Discussion

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.

Comparison of enrichment function characterization

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 ²⁴.

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 ²⁵. As an animal model for studying human obesity, pigs exhibit obesity complications similar to those of humans ²⁶. 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 ²⁷. AMPK activity is associated with increased glucose uptake and fatty acid oxidation, as well as inhibition of glycogen synthase activity and fatty acid synthesis ²⁸. 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 ²⁹, and the AMPKγ3R200Q mutation results in a lack of AMP dependence and elevated basal activity of AMPK ²⁸. 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 ³². Acetoacetyl-CoA synthetase (AACS) is an enzyme that provides acetyl units for the biosynthesis of cholesterol and fatty acids ³³. 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 ³⁴. 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 ³⁵. Moreover, ALDH1L1 is an enzyme that metabolizes folate ³⁶. 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 ³⁷. 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 ³⁸. ACP5, which is associated with excessive fat deposition and the development of atherosclerosis in pigs ³⁹. 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 ⁴². Studies in C2C12 myoblasts have shown that overexpression of PPARD inhibits myotube formation and enhances adipocyte differentiation ⁴³. In mice, knockout of the UBD gene prevents the development of age-related obesity while prolonging lifespan and vitality ⁴⁴. 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 ⁵⁰. PRKAG2 is an important regulator of cardiac metabolism ⁵¹; an SNP of PRKAG2 is associated with insulin sensitivity and is essential for obesity, triglycerides, and HDL cholesterol ⁵². SMPDL3A—a recently identified phosphodiesterase—plays a role in cholesterol and fat metabolism by regulating the liver X-receptor ⁵³. 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 ⁵⁴. In addition, the proteins LEP ⁵⁵, GK ⁵⁶, and PPARA ⁵⁷—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. 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 ⁵⁸.

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.

Abbreviations

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.

Declarations

Data Availability Statement

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.

Author Contributions: Conceptualization, P.Z. and H.Z.; methodology, P.Z., Y,Z.,B.Z., and H.Z.; software, P.Z., Y.Z. and B.Z.; validation, P.Z., B.Z. and H.Z.; formal analysis, P.Z.; investigation, Y.W. and Q.L.; resources, Q.L.; data curation, P.Z, B.Z. and H.Z.; writing—original draft preparation, P.Z.; writing—review and editing, P.Z., B.Z. and H.Z.; visualization, P.Z.; supervision, B.Z. and H.Z.; funding acquisition, H.Z. All authors have read and agreed to the published version of the manuscript

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

Acknowledgments: We thanks Annoroad (Beijing, China) and Beijing Bangfei Bioscience &Technology Co., Ltd (Beijing, China) for the sequencing services and bioinformatic analysis.

Conflicts of Interest: The authors declare no conflict of interest.

References

1. Yang, S. L. et al. Genetic variation and relationships of eighteen Chinese indigenous pig breeds. Genet. Sel. Evol, 35, 657–671 (2003).
2. Dai, F. W. et al. Developmental differences in carcass, meat quality and muscle fibre characteristics between the Landrace and a Chinese native pig. S. Afr. J. Anim. Sci, 39, 267–273 (2009).
3. Cameron, N. D., Warriss, P. D., Porter, S. J. & Enser, M. B. Comparison of duroc and british landrace pigs for meat and eating quality. Meat Sci, 27, 227–247 (1990).
4. Wood, J. D., Kempster, A. J., David, P. J. & Bovey, M. Observations on carcass and meat quality in Duroc, Landrace and Duroc X Landrace pigs. Anim. Prod, 44, 488–488 (1987).
5. Suzuki, A. et al. Carcass composition and meat quality of chinese purebred and European X Chinese crossbred pigs. Meat Sci, 29, 31–41 (1991).
6. Zambonelli, P., Gaffo, E., Zappaterra, M., Bortoluzzi, S. & Davoli, R. Transcriptional profiling of subcutaneous adipose tissue in Italian Large White pigs divergent for backfat thickness. Anim. Genet, 47, 306–323 (2016).
7. Wood, J. D. et al. Fat deposition, fatty acid composition and meat quality: A review. Meat Sci, 78, 343–358 (2008).
8. Newcom, D. W., Baas, T. J., Schwab, C. R. & Stalder, K. J. Genetic and phenotypic relationships between individual subcutaneous backfat layers and percentage of longissimus intramuscular fat in Duroc swine. J. Anim. Sci, 83, 316–323 (2005).
9. Miao, Z. G. et al. Transcriptome analysis reveals differential gene expression in intramuscular adipose tissues of Jinhua and Landrace pigs. J. Vet. Med. Sci, 80, 953–959 (2018).
10. Zhang, Y. F. et al. Subcutaneous and intramuscular fat transcriptomes show large differences in network organization and associations with adipose traits in pigs. Sci. China Life Sci, 64, 1732–1746 (2021).
11. Canovas, A., Quintanilla, R., Amills, M. & Pena, R. N. Muscle transcriptomic profiles in pigs with divergent phenotypes for fatness traits. BMC Genom, 11, 372 (2010).
12. Thompson, A. et al. Tandem mass tags: A novel quantification strategy for comparative analysis of complex protein mixtures by MS/MS (vol 75, pg 1895, 2003). Anal Chem, 78, 4235–4235 (2006).
13. Xing, K. et al. Comparative adipose transcriptome analysis digs out genes related to fat deposition in two pig breeds. Sci. Rep, 9, 12925 (2019).
14. Liu, Y. Z. et al. Transcriptomics Analysis on Excellent Meat Quality Traits of Skeletal Muscles of the Chinese Indigenous Min Pig Compared with the Large White Breed. Int J Mol Sci, 19, 21 (2018).
15. Wu, W. J. et al. Transcriptome analysis reveals the genetic basis of skeletal muscle glycolytic potential based on a pig model., 766, 145157 (2021).
16. Ma, C. et al. TMT-labeled quantitative proteomic analyses on the longissimus dorsi to identify the proteins underlying intramuscular fat content in pigs. J. Proteomics, 213, 103630 (2020).
17. Ran, L. Y. et al. Integrated transcriptomic and proteomic analysis indicated that neurotoxicity of rats with chronic fluorosis may be in mechanism involved in the changed cholinergic pathway and oxidative stress. J. Trace. Elem. Med. Bio, 64, 126688 (2021).
18. Zhang, Y. et al. Comparative transcriptomic and proteomic analyses provide insights into functional genes for hypoxic adaptation in embryos of Tibetan chickens. Sci Rep-Uk, 10, 11213 (2020).
19. Dai, W. T. et al. Complementary transcriptomic and proteomic analyses reveal regulatory mechanisms of milk protein production in dairy cows consuming different forages. Sci. Rep, 7, 44234 (2017).
20. La Merrill, M. et al. Toxicological Function of Adipose Tissue: Focus on Persistent Organic Pollutants. Environ. Health Persp, 121, 162–169 (2013).
21. Pellegrinelli, V., Carobbio, S. & Vidal-Puig, A. Adipose tissue plasticity: how fat depots respond differently to pathophysiological cues., 59, 1075–1088 (2016).
22. Gardan, D. et al. Lipid metabolism and cellular features of skeletal muscle and subcutaneous adipose tissue in pigs differing in IGF-II genotype. Domest. Anim. Endocrin, 34, 45–53 (2008).
23. Wang, S. B. et al. Glucose Utilization, Lipid Metabolism and BMP-Smad Signaling Pathway of Porcine Intramuscular Preadipocytes Compared with Subcutaneous Preadipocytes. Cell Physiol. Biochem, 31, 981–996 (2013).
24. Zhou, G. X. et al. Global comparison of gene expression profiles between intramuscular and subcutaneous adipocytes of neonatal landrace pig using microarray. Meat Sci, 86, 440–450 (2010).
25. Sethi, J. K. Activatin' Human Adipose Progenitors in Obesity., 59, 2354–2357 (2010).
26. Renner, S. et al. Porcine models for studying complications and organ crosstalk in diabetes mellitus. Cell Tissue Res, 380, 341–378 (2020).
27. Roux, M. et al. Characterization of the bovine PRKAG3 gene: structure, polymorphism, and alternative transcripts. Mamm. Genome, 17, 83–92 (2006).
28. Barnes, B. R. et al. The 5 '-AMP-activated protein kinase gamma 3 isoform has a key role in carbohydrate and lipid metabolism in glycolytic skeletal muscle. J. Biol. Chem, 279, 38441–38447 (2004).
29. Chen, J. F. et al. New evidence of alleles (V199I and G52S) at the PRKAG3 (RN) locus affecting pork meat quality. Asian. Austral. J. Anim, 21, 471–477 (2008).
30. Kersten, S. Physiological regulation of lipoprotein lipase. Bba-Mol. Cell Biol. L, 1841, 919–933 (2014).
31. Ullrich, N. F. E., Purnell, J. Q. & Brunzell, J. D. Adipose tissue fatty acid composition in humans with lipoprotein lipase deficiency. J. Invest. Med, 49, 273–275 (2001).
32. Menendez, J. A. & Lupu, R. Fatty acid synthase and the lipogenic phenotype in cancer pathogenesis. Nat. Rev. Cancer, 7, 763–777 (2007).
33. Yamasaki, M. et al. Acetoacetyl-CoA synthetase gene is abundant in rat adipose, and related with fatty acid synthesis in mature adipocytes. Biochem. Bioph. Res. Co, 335, 215–219 (2005).
34. Yamasaki, M., Hasegawa, S., Kitani, T., Hidai, K. & Fukui, T. Differential effects of obesity on acetoacetyl-CoA synthetase gene in rat adipose tissues. Eur. J. Lipid Sci. Tech, 109, 617–622 (2007).
35. Zhao, X. et al. Weighted gene co-expression network analysis reveals potential candidate genes affecting drip loss in pork. Anim. Genet, 51, 855–865 (2020).
36. Krupenko, N. I. et al. Cytosolic 10-formyltetrahydrofolate dehydrogenase regulates glycine metabolism in mouse liver. Sci. Rep, 9, 14937 (2019).
37. Kim, J. Y., Tillison, K. S., Zhou, S. L., Lee, J. H. & Smas, C. M. Differentiation-dependent expression of Adhfe1 in adipogenesis. Arch. Biochem. Biophys, 464, 100–111 (2007).
38. Wu, J. H., Srinivasan, S. V., Neumann, J. C. & Lingrel, J. B. The KLF2 transcription factor does not affect the formation of preadipocytes but inhibits their differentiation into adipocytes., 44, 11098–11105 (2005).
39. Padilla, J. et al. Vascular transcriptional alterations produced by juvenile obesity in Ossabaw swine. Physiol. Genomics, 45, 434–446 (2013).
40. Ma, L. J. et al. Prevention of obesity and insulin resistance in mice lacking plasminogen activator inhibitor 1., 53, 336–346 (2004).
41. Liu, S. J. et al. Transplantation of adipose tissue lacking PAI-1 improves glucose tolerance and attenuates cardiac metabolic abnormalities in high-fat diet-induced obesity. Adipocyte, 9, 170–178 (2020).
42. Meidtner, K. et al. Haplotypes of the porcine peroxisome proliferator-activated receptor delta gene are associated with backfat thickness. Bmc Genet, 10, 76 (2009).
43. Yu, Y. H., Wu, S. C., Cheng, W. T. K., Mersmann, H. J. & Ding, S. T. Ectopic expression of porcine peroxisome proliferator-activated receptor delta regulates adipogenesis in mouse myoblasts. J. Anim Sci, 86, 64–72 (2008).
44. Canaan, A. et al. Extended lifespan and reduced adiposity in mice lacking the FAT10 gene. P. Natl. Acad Sci. U. S. A, 111, 5313–5318 (2014).
45. van Dierendonck, X. A. M. H., Sancerni, T., Alves-Guerra, M. C. & Stienstra, R. The role of uncoupling protein 2 in macrophages and its impact on obesity-induced adipose tissue inflammation and insulin resistance. J. Biol. Chem, 295, 17535–17548 (2020).
46. Vozza, A. et al. UCP2 transports C4 metabolites out of mitochondria, regulating glucose and glutamine oxidation. P. Natl. Acad. Sci. U. S. A, 111, 960–965 (2014).
47. Oktavianthi, S. et al. Uncoupling protein 2 gene polymorphisms are associated with obesity. Cardiovasc Diabetol, 11, 41 (2012).
48. Csernus, K., Pauler, G., Erhardt, E., Lanyi, E. & Molnar, D. Uncoupling protein-2 gene polymorphisms are associated with obesity in Hungarian children. Acta. Paediatr, 102, e200–e204 (2013).
49. Oguzkan-Balci, S. et al. Mitochondrial uncoupling protein 2 (UCP2) gene polymorphisms are associated with childhood obesity and related metabolic disorders. J. Pediatr. Endocr. Met, 26, 277–283 (2013).
50. Yuan, Z. Q., Song, D. G. & Wang, Y. Z. The novel gene pFAM134B positively regulates fat deposition in the subcutaneous fat of Sus scrofa. Biochem Bioph. Res. Co, 454, 554–559 (2014).
51. Torok, R. D. et al. PRKAG2 mutations presenting in infancy. J Inherit Metab Dis, 40, 823–830 (2017).
52. Randrianarisoa, E. et al. AMPK Subunits Harbor Largely Nonoverlapping Genetic Determinants for Body Fat Mass, Glucose Metabolism, and Cholesterol Metabolism. J. Clin. Endocr. Metab, 105, 14–25 (2020).
53. Noto, P. B. et al. Regulation of Sphingomyelin Phosphodiesterase Acid-Like 3A Gene (SMPDL3A) by Liver X Receptors. Mol. Pharmacol, 82, 719–727 (2012).
54. Rached, M. T. et al. Deletion of myeloid IRS2 enhances adipose tissue sympathetic nerve function and limits obesity. Mol. Metab, 20, 38–50 (2019).
55. Martinez-Sanchez, N., There & Again, B. Leptin Actions in White Adipose Tissue. Int. J. Mol. Sci, 21, 6039 (2020).
56. Rahib, L., MacLennan, N. K., Horvath, S., Liao, J. C. & Dipple, K. M. Glycerol kinase deficiency alters expression of genes involved in lipid metabolism, carbohydrate metabolism, and insulin signaling. Eur. J. Hum. Genet, 15, 646–657 (2007).
57. Brocker, C. N. et al. Extrahepatic PPAR modulates fatty acid oxidation and attenuates fasting-induced hepatosteatosis in mice. J. Lipid Res, 59, 2140–2152 (2018).
58. Zhao, C. et al. Knockdown of ubiquitin D inhibits adipogenesis during the differentiation of porcine intramuscular and subcutaneous preadipocytes. Cell Proliferat, 51, e12401 (2018).
59. Warr, A. et al. An improved pig reference genome sequence to enable pig genetics and genomics research. Gigascience, 9, giaa051 (2020).