- Open Access
TMT-based quantitative proteomics analysis reveals the key proteins related with the differentiation process of goat intramuscular adipocytes
BMC Genomics volume 22, Article number: 417 (2021)
Intramuscular adipocytes differentiation is a complex process, which is regulated by various transcription factor, protein factor regulators and signal transduction pathways. However, the proteins and signal pathways that regulates goat intramuscular adipocytes differentiation remains unclear.
In this study, based on nanoscale liquid chromatography mass spectrometry analysis (LC-MS/MS), the tandem mass tag (TMT) labeling analysis was used to investigate the differentially abundant proteins (DAPs) related with the differentiation process of goat intramuscular adipocytes. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes enrichment and protein-protein interaction network analyses were performed for the characterization of the identified DAPs. The candidate proteins were verified by parallel reaction monitoring analysis. As a result, a total of 123 proteins, 70 upregulation proteins and 53 downregulation proteins, were identified as DAPs which may be related with the differentiation process of goat intramuscular adipocytes. Furthermore, the cholesterol metabolism pathway, glucagon signaling pathway and glycolysis / gluconeogenesis pathway were noticed that may be the important signal pathways for goat Intramuscular adipocytes differentiation.
By proteomic comparison between goat intramuscular preadipocytes (P_IMA) and intramuscular adipocytes (IMA), we identified a series protein that might play important role in the goat intramuscular fat differentiation, such as SRSF10, CSRP3, APOH, PPP3R1, CRTC2, FOS, SERPINE1 and AIF1L, could serve as candidates for further elucidate the molecular mechanism of IMF differentiation in goats.
Intramuscular fat (IMF) is the adipose tissue between muscle fibers, also known as marbling. IMF content has a positive effect on meat tenderness, moisture content and palatability [1,2,3,4]. As an important economic trait of lamb production, reasonable IMF content can create greater economic benefits and further improve the taste and quality of meat [2, 5].
Adipose tissue is mainly composed of a large number of adipocytes, which are generated through the proliferation and differentiation of preadipocytes. Study showed that the preadipocytes are present throughout adult life, exhibiting adipose depot specificity . In addition, the differentiation of preadipocytes into adipocytes is a complex process, which is regulated by various transcription factor, protein factor regulators and signal transduction pathway [6,7,8]. In porcine preadipocytes, miR-429 inhibits subcutaneous and intramuscular preadipocytes differentiation while promotes proliferation by directly binding to the 3′-UTRs of KLF9 and p27 . Emerging research on preadipocyte differentiation using proteome and transcriptome analysis better revealed the breadth of adipocyte differentiation. For instance, using HPLC-tandem mass spectrometry and methylated RNA immunoprecipitation (meRIP), found that MTCH2 promotes adipogenesis in pig P_IMA via an m6A-YTHDF1-dependent mechanism .
The differential expression of cellular proteins affects the cell state. Studies have showed that the relationship between protein and mRNA expression levels is a comprehensive result of translation and protein degradation, while the genome-wide correlation between mRNA expression levels and proteins is very low, and study showed that only about 40 % of protein expression differences can use the change of transcription to explain [11,12,13]. However, the significant difference in mRNA expression between conditions is often used for biological discovery. Moreover, many measurements on mRNA can be traced back to comparison with protein data, that is, attributing functional differences between conditions to protein action [14,15,16]. Adipocyte differentiation are complex quantitative traits that is regulated by multiple genes and proteins. Transcriptome and proteomics analysis of intramuscular adipocyte differentiation and fat deposition are critical to goat meat production and its quality. However, the key regulatory proteins in differentiation of goat intramuscular adipocytes still unclear.
Here, we used Tandem Mass Tags (TMT) quantitative proteomics method to analyze and screen the different-expressed proteins between goat P_IMA and IMA. Also, Parallel Reaction Monitoring (PRM) technology was used for validating the quantitative analysis of differentially expressed proteins. The results are of great significance for elucidating the molecular mechanism of IMF differentiation in goats.
Identification of goat P_IMA and IMA differentiation model
Goat IMA differentiation model in vitro was constructed in this study. As shown in Fig. 1A, cells were filled with lipid droplets after 5 days induction. Furthermore, the results of quantitative real time PCR (qRT-PCR) showed that expression of key adipogenic differentiation genes PPARg and LPL were significantly up-regulated (Fig. 1B).
Protein identification and quantitative proteome analysis based on TMT
Three samples of goat P_IMA, induced 0d, and IMA, induced 5d, respectively, were analyzed using TMT proteomics technique to identify the differentially abundant proteins (DAPs). And the experimental strategy is shown in Fig. 2.
Protein we selected containing at least one unique peptide, and the false discovery rate (FDR) of the peptide and protein were < 1 %. In this study, a total of 5929 proteins were quantified that were jointly found in two samples (Supplementary material S1). The quality control results of the proteins were showed in Fig. 3. For comparison between the intramuscular adiposes samples between 0d and 5 day, the proteins’ difference multiple (log2 value) and p-value (-log10 value) were used as the abscissa and ordinate to generate the difference protein volcano map (Fig. 4A). As shown in volcano map, we can divide the proteins into three types clearly, that are up-regulated, down-regulated and unchanged groups. According to the criterion, 123 DAPs were identified in this study, among them, 70 proteins (56.9 %) were upregulated, and 53 proteins (43.1 %) were downregulated (Fig. 4B). To provide visualization of the overall protein change effect, the DAPs of each groups were analyzed in the form of a heatmap with the hierarchical cluster analysis (Fig. 4C). Moreover, the top 20 proteins of upregulated or downregulated were listed in Table 1.
GO function and KEGG Pathway enrichment analysis
For further exploring the biological functions of DAPs, we performed GO function annotation analysis on the obtained DAPs, and the top 20 of biological processes, cellular components, and molecular functions were shown in Fig. 5. Using pie chart, we presented the enrichment results, and we found that these DAPs were involved in a variety of biological processes, such as metabolic process (42 %), transport regulation (16 %), signaling pathway (9 %) and muscle contraction (9 %). (Fig. 5A, Supplementary material S2). Cell component analysis showed that these proteins were mainly derived from organelle (36 %), cell (25 %), complex (18 %) and membrane (9 %) (Fig. 5B, Supplementary material S3). Molecular function analysis revealed that 86 % of these DAPs participated in binding and 8 % was involved in enzyme activity (Fig. 5C, Supplementary material S4). Subsequently, KEGG enrichment analysis were performed to investigate the main biochemical metabolic and signal transduction pathways in which DAPs participated. Top 20 pathways of the DAPs mapped were shown in Fig. 6A, (Supplementary material S5). We found that there were 4 proteins (XP_017910927.1, XP_017909416.1, XP_005686153.1 and XP_017908084.1) enriched in kaposi′s sarcoma-associated herpesvirus infection with the lowest P-value among identified pathways. Furthermore, three proteins were found enriched in the cholesterol metabolism signaling pathway including beta-2-glycoprotein 1 isoform X1(APOH), apolipoprotein A-II(APOA2) and sterol O-acyltransferase 1(SOAT1) (Fig. 6B). Four proteins enriched in the Glucagon signaling pathway were phosphoglycerate mutase 2(PGAM2), calcineurin subunit B type 1 isoform X1(PPP3R1), interferon regulatory factor 2-binding protein-like (IRF2BPL) and CREB-regulated transcription coactivator 2 isoform X1(CRTC2) (Fig. 6C).
Network analysis of protein-protein interactions
Using StringDB database and Cytoscape software, we constructed a protein-protein interactions (PPI) network for the DAPs. Our results showed that glyceraldehyde-3-phosphate ehydrogenase-like (LOC102181016) exhibited the highest connectivity degree. There are three proteins Fructo-oligosaccharides (FOS / c-Fos), Plasminogen activator inhibitor 1 (PAI-1/SERPINE1) and Phosphoglycerate mutase 2 (PGAM2) exhibiting the highest degree of connectivity with LOC102181016 (Fig. 7, Supplementary material S6).
PRM Validation of TMT-Based Results
To validate the differential proteins between P_IMA and IMA, eight DAPs that may be related to fat formation were selected for verifying by PRM quantitative analysis, which are Cysteine and glycine rich protein 3(CSRP3), Fibulin 1(FBLN1), Nexilin F-actin binding protein(NEXN), Serine and arginine rich splicing factor 10(SRSF10), LIM domain binding 3(LDB3), Alpha-2-macroglobulin isoform X1(ALPHA2), Allograft inflammatory factor 1 like (AIF1L) and Glutathione peroxidase 4(GPX4). According to the quantitative value of the relative expression of the target proteins in the sample, using T-test to calculate the expression difference. After normalized, the results of the relative expression of the proteins were shown in Fig. 8. As shown, there is a significant difference between P_IMA and IMA in the relative enrichment of proteins in CSRP3, FBLN1, SRSF10, LDB3 (P < 0.01) and AIF1L(P < 0.05). Comparing P_IMA to IMA in this study, CSRP3, FBLN1, SRSF10, LDB3, AIF1L and GPX4 have significant differences in both the TMT and PRM analysis. On the whole, these results are in agreement with the findings in TMT-based quantitative analysis.
Cellular differentiation is a complex process, in which, cells are affected by single-cell communication and the extracellular environment and thus differ from each other in function and morphology. However, one of the final results of cell differentiation is to produce cells with specific functions. For example, skeletal muscle cells produce large amounts of actin and myosin, and red blood cells produce hemoglobin [17,18,19]. Because of the production of these proteins the cells can perform unique functions. In adipocytes, proteins that promote adipocyte proliferation and lipid deposition are regarded as functional proteins in the process of adipocyte differentiation . Therefore, searching for DAPs during the differentiation process of adipocytes is essential to explore the mechanism of fat deposition.
In this study, a total of 5929 proteins were identified, including 123 DAPs that are 70 upregulated proteins and 53 downregulated proteins. The up-regulated and down-regulated proteins with the smallest P-value were DTD1 and ACTN2. The d-tyrosyl-tRNA deacylase 1 (DTD1) plays an important role in metabolic pathways and activation of cellular immune responses [20,21,22]. The actinin alpha 2 (ACTN2) is one of four encoding isoforms of α-actinin, and genome-wide association and multi-omic analyses reveal ACTN2 as a gene being linked to heart failure . Also, microarray analyses showed that ACTN2 is an important differentially expressed heart-related gene for pig heart steatosis and hypertrophy induced by high-energy diet . Furthermore, the most prominent upregulated protein in our result was serine and arginine rich splicing factor 10 (SRSF10). Previous study declared that SRSF10 as an essential regulator for adipocyte differentiation could controls the production of lipin1α and thus promotes adipocyte differentiation in mice . The most prominent downregulated protein was cysteine and glycine rich protein 3 (CSRP3). CSRP3, highly expressed under insulin-sensitive conditions, was highly inducible protein that plays a key role in regulating glucose homeostasis in skeletal muscle. In addition, knockdown of CSRP3 suppressed chicken satellite cell differentiation by regulating Smad3 phosphorylation in the TGF-β signaling pathway [26, 27]. Other research showed that it is related to muscle fiber hypertrophy . These finding suggested that SRSF10 and CSRP3 positively influence the goat IMF differentiation.
GO analysis on the DAPs showed that, most of the DAPs coming from organelles, may perform molecular functions in binding ways and may participate in metabolic processes. Two pathways, cholesterol metabolism and glucagon signaling pathway, were noticed in our study by KEGG pathway enrichment. And both of the two pathways may play essential roles in lipid metabolism and cell differentiation [29,30,31,32]. From the KEGG enrichment analysis, we found two pathways may involve in regulating adipocyte differentiation, that were cholesterol metabolism signaling pathway and glucagon signaling pathway. Among them, APOH is involved in lipid metabolism and synthesis . APOA2 is related to obesity, dyslipidemia and lipid metabolism [34, 35]. SOAT1 is involved in atherosclerosis, cholesterol content, glucose and lipid metabolism, study have found that overexpression of ACAT1/2 encoded by SOAT1 can significantly inhibit the differentiation of 3T3-L1 preadipocytes [36,37,38,39]. PGAM2 is a housekeeping enzyme, involved in the process of sugar metabolism, and plays an important role in muscle growth and development . Studies found that PPP3r is a ubiquitously expressed calcium-sensitive serine-threonine phosphatase, and PPP3r KO mice increase energy expenditure. In addition, skeletal muscle specific ablation Ppp3r1 promotes overall number of fat cells per fat pad . IRF2BPL as a transcriptional cofactor, is a new participant in the regulation of cell homeostasis, and also is a new genetic causes for disorders in dystonia [42, 43]. CRTC2, as a critical mediator of mTOR, can induce SREBP-1 processing and enhancement of de novo lipogenesis. mTORC1 regulates the differentiation of beige adipocytes via regulated transcriptional coactivator 2 (CRTC2) . Also, CRTC2 could modulate triglyceride synthesis through regulating of SREBP‐1 maturation . These results suggested the possible signal pathways and proteins involved in goat IMF differentiation.
According to the PPI network, we found LOC102181016 that enriched in glycolysis/gluconeogenesis pathway showed the highest degree of connectivity in DAPs, and glycolysis/gluconeogenesis pathway was a key pathway of IMF deposition [46,47,48]. The other three strongest interactions with LOC102181016 were FOS, SERPINE1 and PGAM2. FOS belongs to the activator protein 1 (AP-1) superfamily of transcription factors, activation AP-1protein transcription factor Fra-2 in adipocytes in vivo increased differentiation and apoptosis of adipocytes [49, 50]. Moreover, RARγ-C-Fos-PPARγ2 signaling was critical for ATRA-inhibited adipocyte differentiation . SERPINE1 has been reported to be related to adipocyte differentiation, and inhibition of the SERPINE1 in 3T3-L1 adipocytes can increase the expression of PPARγ, promote adipocyte differentiation, and decrease insulin resistance [52, 53]. PGAM2, insulin/IGF1-PI3K-dependent of a glycolytic enzyme, is necessary for nucleolus structure and RNA synthesis . Totally, the exact role of these interacting proteins in goat IMF differentiation still needs to be verified by further experiments.
Using PRM technique, we analyzed the potential proteins for adipogenic differentiation, which exist in top 20 up/down regulated DAPs. According to the results, we found that 5 proteins were differently enriched between goat P_IMA and IMA, that were CSRP3, SRSF10, FBLN1, LDB3 and AIF1L. And the role of CSRP3 and SRSF10 were shown in above discussion. FBLN1 as an extracellular matrix protein, is necessary for osteoblast differentiation in mice. Moreover, it was also secreted in mesenchymal stem cells which derived from human fat [55, 56]. LDB3 was well known for its role in myofibrillar myopathies. Recent study found that the LDB3 promoter responds to lipid uptake in human adipose tissue . AIF1L protein seemed to exist in all adipose tissue. A downregulated expression of AIF1L in intermediate stage of 3T3L1 cell differentiation might limit adipogenesis and/or lipogenesis .
By proteomic comparison between goat P_IMA and IMA, we identified a series protein that might play important role in the goat IMF differentiation, including SRSF10, CSRP3, APOH, PPP3R1, CRTC2, FOS, SERPINE1 and AIF1L. These proteins are of substantial significance for in-depth study of the molecular mechanism of goat IMF differentiation.
Materials and methods
Establish the model of goat intramuscular adipocyte differentiation
The samples were collected from the longissimus dorsi of three healthy Jianzhou goats whose age is 7 days old. The experimental animals were injected barbiturate injection into the intraperitoneal cavity at a dose of 100 mg/kg, then bled to death. The carcasses are temporarily stored in freezer and then handed over to a professional solid waste disposal company for unified disposal. The experiment was approved by the Academic Committee of Southwest Minzu University (Chengdu, Sichuan, China), and all the experiment complied with the requirements of the directory of the Ethical Treatment of Experimental Animals of China. Detailed procedure for the collection of intramuscular preadipocytes have previously been published [59, 60].
The goat intramuscular preadipocytes were divided into two groups, each with three biological replicates, and cultured with DMEM/F12 (Hyclone, USA) containing 10 % (v/v) fetal bovine serum (FBS, Hyclone, USA) in 5 % CO2 and 37 °C incubator. The first group preadipocytes, P_IMA, were collected for total protein extracted and RNA extraction when the preadipocytes reached 90 % or more in 6-well plates. The second group, IMA, its medium was replaced with adipogenic inducer of DMEM containing 10 % FBS and 50 µmol•L-1 oleic acid medium, when the preadipocytes reached 90 % in 6-well plates. After 5 days of induction, the cell pellet was collected for total protein extraction and RNA extraction. The medium was refreshed every two days.
Oil Red O Staining and qRT-PCR
Before staining, the P_IMA and IMA were fixed with 10 % formaldehyde for 30 min then were stained by Oil red O working solution for 20 min. Then the cells were washed by PBS for three times and photographed under microscope. Total RNA was extracted using TRIzol (TaKaRa, Otsu, Japan) and stored at -80℃. RevertAid First Strand cDNA Synthesis Kit (TaKaRa, Otsu, Japan) was used to reverse transcription of mRNA according to manufacturer instructions. Ubiquitously expressed transcript (UXT) was used as endogenous control for mRNA. qRT-PCR primers were designed using Primer Premier 5, and listed in Table 2. Reaction volume for qRT-PCR was 20 µL, including cDNA 1µL, F/R primers 1 µL, double-distilled water 7 µL and 10 µL TB Green™ Premix Ex Taq™ II (TaKaRa, Otsu, Japan).
TMT sample preparation
The cells of PIMA and IMA we collected and minced individually in liquid nitrogen and lysed with lysis buffer containing 0.2 % SDS, 6 M Urea and 100 mM NH4HCO3 (pH 8), followed by ultrasonication on ice for 5 min. Subsequently, the supernatant was transferred into a clean tube, and mixed with 2 mM DTT for 1 h at 56 °C and subsequently alkylated with enough Iodoacetamide for 1 h at room temperature in the dark. Then 4 times volume of pre-cooled acetone was added to the samples and incubated at -20 °C for 2 h. The sample was centrifuged at 12 000 rpm for 15 min at 4 °C and the precipitation was collected. After washing pellet twice with cold acetone, it was dissolved by dissolution buffer solution containing 6 M urea and 0.1 M triethylammonium bicarbonate (TEAB, pH 8.5). The protein concentration was determined by the Bradford protein assay (Bio-Rad, USA). The 12 % SDS-PAGE electrophoresis was performed for 20 µg of proteins extracted from each sample with loading buffer, and the gel was stained by Coomassie Blue R-250 for showing the protein bands.
The 120 µg protein of each sample, whose final volume of 100 µL with dissolution buffer solution, was digested with 3µL of Trypsin Gold (1 µg/µL, Promega, USA) and 500 µL of TEAB buffer (50 mM) at 37 °C during one night. Adding the same volume of 1 % formic acid, the digestion mixture was centrifuged at 12 000 rpm for 5 min. The supernatant was desalted with C18 cartridge and washed three times by 1 mL of cleaning solution containing 0.1 % formic acid and 4 % acetonitrile, followed by washing twice with 0.4 mL of eluent containing 0.1 % formic acid and 75 % acetonitrile for removing the high urea and finally was dried by vacuum centrifugation.
The desalted peptides were labeled with TMT 6-plex reagents (Thermo Fisher Scientific, USA), following the manufacturer’s instructions. In brief, 1 unit of labeling reagent was used for 0.1 mg of peptide. The Peptides, dissolved in 100 µL of 0.1 M TEAB, were incubated with the labeling reagent, which was dissolved in 41 µL of acetonitrile, for 2 h at room temperature. Then, three pooled fractions of goat intramuscular preadipocytes group have been labeled with 126, 127, and 128 tags, while those of goat intramuscular adipocytes induced at Day 5 group have been labeled with 129, 130, and 131 tags, individually. The reaction was stopped with 8 % ammonium hydroxide. The differently labeled peptides were mixed equally and desalted by peptide desalting spin columns (Thermo Fisher, USA).
Mobile phases A containing 2 % acetonitrile and the ammonium hydroxide (pH 10.0) and B containing 98 % acetonitrile and the ammonium hydroxide (pH 10.0) were used for developing a gradient elution. The TMT-labeled peptides were dissolved in 1 mL of phases A. The 1 mL supernatant was fractionated using a C18 column (4.6 × 250 mm, 5 μm) on a L3000 HPLC (Rigol, China) and the column oven was set as 50 °C. The solvent gradient was listed as follows: 97 % A and 3 % B for 0 min; 95–97 % A and 3–5 % B for 10 min; 80–95 % A and 5–20 % B for 20 min; 60–80 % A and 20–40 % B for 18 min; 50–60 % A and 40–50 % B for 2 min; 30–50 % A and 50–70 % B for 3 min; 0–30 % A and 70–100 % B for 1 min; 0 % A and 0-100 % B for 4 min, and 0 % A and B for 12 min. The eluates were collected once in 1 min and merged into 10 fractions finally, which were monitored at UV 214 nm. All fractions were dried under vacuum and re-dissolved in 0.1 % formic acid.
The shotgun proteomics analysis was performed by the EASY-nLCTM 1200 UHPLC system (Thermo Fisher Scientific, USA) coupled with the Orbitrap Q Exactive HF-X mass spectrometer (Thermo Fisher Scientific, USA). The sample with 2 µg of the total peptides was injected onto the homemade pre-column (2 cm×100 μm, 5 μm). Then, the peptides were separated using the homemade analytical column (15 cm×150 μm, 1.9 μm) in 90 min at a flow rate of 600 nL/min, with the eluent B (0.1 % formic acid in 80 % acetonitrile) in eluent A (0.1 % formic acid in H2O) for the TMT labeling 6-plex. The solvent gradient was set as follows: 95 % A and 5 % B for 0 min; 90–95 % A and 5–10 % B for 2 min; 60–90 % A and 10–40 % B for 80 min; 45–60 % A and 40–55 % B for 2 min; 10–45 % A and 55–90 % B for 1 min, and 0–10 % A and 90–100 % B for 5 min. The Orbitrap Q Exactive HF-X mass spectrometer was operated in the data-dependent acquisition mode with the spray voltage of 2.3 kV and capillary temperature of 320 °C, when the Nanospray Flex™ as the ion source. The full-scan range from 350 to 1500 m/z was acquired with a mass resolution of 6 × 104 at 200 m/z. The automatic gain control (AGC) target was set as 3 × 106, while the maximum ion injection time was 20 ms. The 40 most abundant precursor ions selected from the full-scan were fragmented using the higher energy collisional dissociation (HCD) for the MS/MS scans. The MS/MS scans were set as follows: a mass resolution of 1.5 × 104 at 200 m/z, an AGC target value of 1 × 105, a maximum ion injection time of 45 ms, and a normalized collision energy of 32 %.
The obtained spectrums were searched against NCBI nr database (X101SC19051573-Z01-Capra-hircus) by the search engines: Proteome Discoverer 2.2 (PD 2.2, Fisher Scientific, USA). A mass tolerance as10 ppm for precursor ion scan and a mass tolerance as 0.02 Da for the product ion was set as the searched parameter. For reducing the false positive rate and increasing the quality of analysis results, the confident peptide spectrums were matched that confidence more than 99 % while false discovery rate (FDR) less than 1 %. Also, the distribution of the peptide length, precursor ion tolerance, unique peptides number, protein coverage, and protein mass were regarded as the key indicators for the verification of total protein quality. The differentially expressed proteins in this study were defined based on the following criteria: the TMT ratio being > 1.5 or < 0.67 with a P-value < 0.05. While T-Test was used to analyze the proteins quantitation results statistically [61, 62].
Bioinformatics analysis of DAPs
The GO analysis was performed by interproscan-5 program against the non-redundant protein databases . According the biological processes, cellular components, and molecular functions, DAPs are classified. Moreover, the DAPs pathways were annotated using the KEGG database, a collection of databases dealing with genomes, diseases, biological pathways, drugs, and chemical materials (http://www.genome.jp/kegg/) [64, 65]. Based on the goat (Capra hircus) species, the potential interacting partners of DAPs were predicted by the StringDB server (http://string-db.org/) . Using Cytoscape software to visualize the prediction results obtained from the StringDB database.
Parallel reaction monitoring validation
To verify the results of TMT analysis coupled with LC-MS/MS, the proteins were extracted and enzymatically hydrolyzed according to the above method. Same amount of trypsin treated-peptide was taken of each sample, and labeled peptide DSPSAPVNVTVR (red bold V for heavy isotope labeling) as an internal standard of each sample and the peptides were dissolved in 0.1 % formic acid (solvent A) and solvent B (0.1 % formic acid in 80 % acetonitrile). PRM mass spectrometric analysis was performed using EASY-nLCTM 1200 UHPLC system (Thermo Fisher Scientific, USA). Liquid chromatography elution gradient was showed in Table 3. The full scan mass spectrum resolution was set to 60,000 (200 m/z), the maximum C-trap and IT were 3×106 and 20 ms respectively. The PRM resolution was set to 30,000(200 m/z), the maximum C-trap and IT were 5×104 and 80 ms respectively. Normalized collision energy is 27. The data was analyzed by Skyline software, and the peak area was corrected using the internal standard peptide.
Availability of data and materials
The datasets generated during the current study are deposited in Capra hircus -- NCBI (https://www.ncbi.nlm.nih.gov/protein).
Differentially abundant proteins
The tandem mass tag
Liquid chromatography mass spectrometry analysis
Parallel Reaction Monitoring
Kyoto Encyclopedia of Genes and Genomes
Quantitative real time PCR
- PPARg :
Preoxisome proliferator-activated receptor gamma
- LPL :
Baik M, Kang HJ, Park SJ, Na SW, Piao M, Kim SY, et al. TRIENNIAL GROWTH AND DEVELOPMENT SYMPOSIUM: Molecular mechanisms related to bovine intramuscular fat deposition in the longissimus muscle. J Anim Sci. 2017;95:2284–303.
Liu K, Yu W, Wei W, Zhang X, Tian Y, Sherif M, et al. Melatonin reduces intramuscular fat deposition by promoting lipolysis and increasing mitochondrial function. J Lipid Res. 2019;60:767–82.
Nishimura T, Hattori A, Takahashi K. Structural changes in intramuscular connective tissue during the fattening of Japanese black cattle: effect of marbling on beef tenderization. J Anim Sci. 1999;77:93–104.
Zhao X-H, Yang Z-Q, Bao L-B, Wang C-Y, -Zhou S, Gong J-M, et al. Daidzein enhances intramuscular fat deposition and improves meat quality in finishing steers. Exp Biol Med (Maywood). 2015;240:1152–7.
Hocquette JF, Gondret F, Baéza E, Médale F, Jurie C, Pethick DW. Intramuscular fat content in meat-producing animals: development, genetic and nutritional control, and identification of putative markers. Animal. 2010;4:303–19.
Sarantopoulos CN, Banyard DA, Ziegler ME, Sun B, Shaterian A, Widgerow AD. Elucidating the Preadipocyte and Its Role in Adipocyte Formation: a Comprehensive Review. Stem Cell Rev Rep. 2018;14:27–42.
Cristancho AG, Lazar MA. Forming functional fat: a growing understanding of adipocyte differentiation. Nat Rev Mol Cell Biol. 2011;12:722–34.
Wei S, Zhang M, Zheng Y, Yan P. ZBTB16 Overexpression Enhances White Adipogenesis and Induces Brown-Like Adipocyte Formation of Bovine White Intramuscular Preadipocytes. Cell Physiol Biochem. 2018;48:2528–38.
Peng Y, Chen F-F, Ge J, Zhu J-Y, Shi X-E, Li X, et al. miR-429 Inhibits Differentiation and Promotes Proliferation in Porcine Preadipocytes. Int J Mol Sci. 2016;17(12):2047.
Jiang Q, Sun B, Liu Q, Cai M, Wu R, Wang F, et al. MTCH2 promotes adipogenesis in intramuscular preadipocytes via an m6A-YTHDF1-dependent mechanism. FASEB J. 2019;33:2971–81.
de Sousa Abreu R, Penalva LO, Marcotte EM, Vogel C. Global signatures of protein and mRNA expression levels. Mol Biosyst. 2009;5:1512–26.
Vogel C, Marcotte EM. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat Rev Genet. 2012;13:227–32.
Koussounadis A, Langdon SP, Um IH, Harrison DJ, Smith VA. Relationship between differentially expressed mRNA and mRNA-protein correlations in a xenograft model system. Sci Rep. 2015;5:10775.
Wang H, Wang Q, Pape UJ, Shen B, Huang J, Wu B, et al. Systematic investigation of global coordination among mRNA and protein in cellular society. BMC Genomics. 2010;11:364.
Vogel C, Abreu R de S, Ko D, Le S-Y, Shapiro BA, Burns SC, et al. Sequence signatures and mRNA concentration can explain two-thirds of protein abundance variation in a human cell line. Mol Syst Biol. 2010;6:400.
Selbach M, Schwanhäusser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N. Widespread changes in protein synthesis induced by microRNAs. Nature. 2008;455:58–63.
Dzierzak E, Philipsen S. Erythropoiesis: development and differentiation. Cold Spring Harb Perspect Med. 2013;3:a011601.
Ali AT, Hochfeld WE, Myburgh R, Pepper MS. Adipocyte and adipogenesis. Eur J Cell Biol. 2013;92:229–36. https://doi.org/10.1016/j.ejcb.2013.06.001.
Pérez-Pérez R, López JA, García-Santos E, Camafeita E, Gómez-Serrano M, Ortega-Delgado FJ, et al. Uncovering suitable reference proteins for expression studies in human adipose tissue with relevance to obesity. PLoS One. 2012;7:e30326.
Wirleitner B, Neurauter G, Schröcksnadel K, Frick B, Fuchs D. Interferon-gamma-induced conversion of tryptophan: immunologic and neuropsychiatric aspects. Curr Med Chem. 2003;10:1581–91.
Pasaje CFA, Bae JS, Park B-L, Jang A-S, Uh S-T, Kim M-K, et al. Association analysis of DTD1 gene variations with aspirin-intolerance in asthmatics. Int J Mol Med. 2011;28:129–37.
Seroogy CM, Gern JE. The role of T regulatory cells in asthma. J Allergy Clin Immunol. 2005;116:996–9.
Arvanitis M, Tampakakis E, Zhang Y, Wang W, Auton A, 23andMe Research Team, et al. Genome-wide association and multi-omic analyses reveal ACTN2 as a gene linked to heart failure. Nat Commun. 2020;11:1122.
Xia J, Zhang Y, Xin L, Kong S, Chen Y, Yang S, et al. Global Transcriptomic Profiling of Cardiac Hypertrophy and Fatty Heart Induced by Long-Term High-Energy Diet in Bama Miniature Pigs. PLoS One. 2015;10:e0132420.
Li H, Cheng Y, Wu W, Liu Y, Wei N, Feng X, et al. SRSF10 regulates alternative splicing and is required for adipocyte differentiation. Mol Cell Biol. 2014;34:2198–207.
Hernandez-Carretero A, Weber N, LaBarge SA, Peterka V, Doan NYT, Schenk S, et al. Cysteine- and glycine-rich protein 3 regulates glucose homeostasis in skeletal muscle. Am J Physiol Endocrinol Metab. 2018;315:E267-78.
Han S, Cui C, Wang Y, He H, Liu Z, Shen X, et al. Knockdown of CSRP3 inhibits differentiation of chicken satellite cells by promoting TGF-β/Smad3 signaling. Gene. 2019;707:36–43.
Li D, Pan Z, Zhang K, Yu M, Yu D, Lu Y, et al. Identification of the Differentially Expressed Genes of Muscle Growth and Intramuscular Fat Metabolism in the Development Stage of Yellow Broilers. Genes (Basel). 2020;11.
Getz GS, Reardon CA. Apoprotein E and Reverse Cholesterol Transport. Int J Mol Sci. 2018;19(11):3479.
Lee HM, Joo BS, Lee CH, Kim HY, Ock JH, Lee YS. Effect of Glucagon-like Peptide-1 on the Differentiation of Adipose-derived Stem Cells into Osteoblasts and Adipocytes. J Menopausal Med. 2015;21:93–103.
Liu R, Li N, Lin Y, Wang M, Peng Y, Lewi K, et al. Glucagon Like Peptide-1 Promotes Adipocyte Differentiation via the Wnt4 Mediated Sequestering of Beta-Catenin. PLoS One. 2016;11:e0160212.
Zhu W, Deng X, Peng J, Zou B, Li C. A-type ECG and EGCG dimers inhibit 3T3-L1 differentiation by binding to cholesterol in lipid rafts. J Nutr Biochem. 2017;48:62–73.
Romao JM, He ML, McAllister TA, Guan LL. Effect of age on bovine subcutaneous fat proteome: molecular mechanisms of physiological variations during beef cattle growth. J Anim Sci. 2014;92:3316–27.
Domínguez-Reyes T, Astudillo-López CC, Salgado-Goytia L, Muñoz-Valle JF, Salgado-Bernabé AB, Guzmán-Guzmán IP, et al. Interaction of dietary fat intake with APOA2, APOA5 and LEPR polymorphisms and its relationship with obesity and dyslipidemia in young subjects. Lipids Health Dis. 2015;14:106.
Zaki ME, Amr KS, Abdel-Hamid M. APOA2 Polymorphism in Relation to Obesity and Lipid Metabolism. Cholesterol. 2013;2013:289481.
Wu N, Li R-Q, Li L. SOAT1 deficiency attenuates atherosclerosis by regulating inflammation and cholesterol transportation via HO-1 pathway. Biochem Biophys Res Commun. 2018;501:343–50.
Pioche T, Skiba F, Bernadet M-D, Seiliez I, Massimino W, Houssier M, et al. Kinetic study of the expression of genes related to hepatic steatosis, glucose and lipid metabolism, and cellular stress during overfeeding in mule ducks. Am J Physiol Regul Integr Comp Physiol. 2020;318:R453-67.
Do DN, Schenkel F, Miglior F, Zhao X, Ibeagha-Awemu EM. Targeted genotyping to identify potential functional variants associated with cholesterol content in bovine milk. Anim Genet. 2020;51:200–9.
Xu Y, Du X, Turner N, Brown AJ, Yang H. Enhanced acyl-CoA:cholesterol acyltransferase activity increases cholesterol levels on the lipid droplet surface and impairs adipocyte function. J Biol Chem. 2019;294:19306–21.
Qiu H, Zhao S, Xu X, Yerle M, Liu B. Assignment and expression patterns of porcine muscle-specific isoform of phosphoglycerate mutase gene. J Genet Genomics. 2008;35:257–60.
Pfluger PT, Kabra DG, Aichler M, Schriever SC, Pfuhlmann K, García VC, et al. Calcineurin Links Mitochondrial Elongation with Energy Metabolism. Cell Metab. 2015;22:838–50.
Steel D, Kurian MA. Recent genetic advances in early-onset dystonia. Curr Opin Neurol. 2020;33:500–7.
Ramalho-Oliveira R, Oliveira-Vieira B, Viola JPB. IRF2BP2: A new player in the regulation of cell homeostasis. J Leukoc Biol. 2019;106:717–23.
Zhang X, Luo Y, Wang C, Ding X, Yang X, Wu D, et al. Adipose mTORC1 Suppresses Prostaglandin Signaling and Beige Adipogenesis via the CRTC2-COX-2 Pathway. Cell Rep. 2018;24:3180–93.
Matsuzaka T, Shimano H. Novel role for the CRTC2 in lipid homeostasis. J Diabetes Investig. 2016;7:677–9.
Qi K, Men X, Wu J, Xu Z. Rearing pattern alters porcine myofiber type, fat deposition, associated microbial communities and functional capacity. BMC Microbiol. 2019;19:181.
Huang J, Feng X, Zhu R, Guo D, Wei Y, Cao X, et al. Comparative transcriptome analysis reveals that PCK1 is a potential gene affecting IMF deposition in buffalo. BMC Genomics. 2020;21:710.
Kwasiborski A, Sayd T, Chambon C, Santé-Lhoutellier V, Rocha D, Terlouw C. Pig Longissimus lumborum proteome: Part II: Relationships between protein content and meat quality. Meat Sci. 2008;80:982–96.
Choi Y, Rosewell KL, Brännström M, Akin JW, Curry TE, Jo M. FOS, a Critical Downstream Mediator of PGR and EGF Signaling Necessary for Ovulatory Prostaglandins in the Human Ovary. J Clin Endocrinol Metab. 2018;103:4241–52.
Luther J, Ubieta K, Hannemann N, Jimenez M, Garcia M, Zech C, et al. Fra-2/AP-1 controls adipocyte differentiation and survival by regulating PPARγ and hypoxia. Cell Death Differ. 2014;21:655–64.
Wang X, Yang P, Liu J, Wu H, Yu W, Zhang T, et al. RARγ-C-Fos-PPARγ2 signaling rather than ROS generation is critical for all-trans retinoic acid-inhibited adipocyte differentiation. Biochimie. 2014;106:121–30.
Liang X, Kanjanabuch T, Mao S-L, Hao C-M, Tang Y-W, Declerck PJ, et al. Plasminogen activator inhibitor-1 modulates adipocyte differentiation. Am J Physiol Endocrinol Metab. 2006;290:E103-13.
Huang W, Guo Y, Du W, Zhang X, Li A, Miao X. Global transcriptome analysis identifies differentially expressed genes related to lipid metabolism in Wagyu and Holstein cattle. Sci Rep. 2017;7:5278.
Gizak A, Grenda M, Mamczur P, Wisniewski J, Sucharski F, Silberring J, et al. Insulin/IGF1-PI3K-dependent nucleolar localization of a glycolytic enzyme–phosphoglycerate mutase 2, is necessary for proper structure of nucleolus and RNA synthesis. Oncotarget. 2015;6:17237–50.
Cooley MA, Harikrishnan K, Oppel JA, Miler SF, Barth JL, Haycraft CJ, et al. Fibulin-1 is required for bone formation and Bmp-2-mediated induction of Osterix. Bone. 2014;69:30–8.
Nakashima Y, Nahar S, Miyagi-Shiohira C, Kinjo T, Toyoda Z, Kobayashi N, et al. A Liquid Chromatography with Tandem Mass Spectrometry-Based Proteomic Analysis of the Proteins Secreted by Human Adipose-Derived Mesenchymal Stem Cells. Cell Transplant. 2018;27:1469–94.
Garske KM, Pan DZ, Miao Z, Bhagat YV, Comenho C, Robles CR, et al. Reverse gene-environment interaction approach to identify variants influencing body-mass index in humans. Nat Metab. 2019;1:630–42.
Parikh D, Riascos-Bernal DF, Egaña-Gorroño L, Jayakumar S, Almonte V, Chinnasamy P, et al. Allograft inflammatory factor-1-like is not essential for age dependent weight gain or HFD-induced obesity and glucose insensitivity. Sci Rep. 2020;10:3594.
Xu Q, Lin S, Wang Y, Zhu J, Lin Y. Fibroblast growth factor 10 (FGF10) promotes the adipogenesis of intramuscular preadipocytes in goat. Mol Biol Rep. 2018;45:1881–8.
Xu Q, Lin S, Li Q, Lin Y, Xiong Y, Zhu J, et al. Fibroblast growth factor 21 regulates lipid accumulation and adipogenesis in goat intramuscular adipocyte. Anim Biotechnol. 2019;12:1–9.
Ma C, Wang W, Wang Y, Sun Y, Kang L, Zhang Q, et al. TMT-labeled quantitative proteomic analyses on the longissimus dorsi to identify the proteins underlying intramuscular fat content in pigs. J Proteomics. 2020;213:103630.
Li L, Huang Q, Barbero M, Liu L, Nguyen T, Xu A, et al. Proteins and Signaling Pathways Response to Dry Needling Combined with Static Stretching Treatment for Chronic Myofascial Pain in a RAT Model: An Explorative Proteomic Study. Int J Mol Sci. 2019;20(3):564.
Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28:1947–51.
Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9. 1: protein-protein interaction networks, with increased coverage and integration. Nucleic acids research. 2012;41:D808-15.
We would like to thank the Beijing Novogene Technology Company for the assistance in original data analysis.
This work was supported by National Natural Sciences Foundation of China (No. 32072723, 31672395).
Ethics approval and consent to participate
All experimental procedures involving animals were performed in accordance with the guidelines and regulations approved by The Animal Care and Use Committee of the Southwest Minzu University (Chengdu, Sichuan, China). In addition, the goats’ owner was agreed us and signed the informed consent form that the experiment would be performed under the requirement of Animal Care and Use Committee of Southwest Minzu University (No. 18032). All goats were humanely sacrificed and all efforts were made to minimize the suffering.
Consent for publication
No potential conflict of interest relevant to this article was reported.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Totally proteins found in two samples.
GO functional enrichment analysis of biological processes involved in DAPs.
GO functional enrichment analysis of cell component involved in DAPs.
GO functional enrichment analysis of molecular function involved in DAPs.
KEGG enrichment analysis the main biochemical metabolic and signal transduction pathways of the DAPs.
The data of protein-protein interactions.
About this article
Cite this article
Du, Y., Wang, Y., Xu, Q. et al. TMT-based quantitative proteomics analysis reveals the key proteins related with the differentiation process of goat intramuscular adipocytes. BMC Genomics 22, 417 (2021). https://doi.org/10.1186/s12864-021-07730-y
- Differentially abundant proteins
- Intramuscular preadipocytes
- Tandem mass tag
- Parallel Reaction Monitoring