Next Article in Journal
Effect on Broiler Production Performance and Meat Quality of Feeding Ulva lactuca Supplemented with Carbohydrases
Previous Article in Journal
Coexistence between Humans and ‘Misunderstood’ Domestic Cats in the Anthropocene: Exploring Behavioural Plasticity as a Gatekeeper of Evolution
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 12
Issue 13
10.3390/ani12131718
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proteomics Insights into the Gene Network of cis9, trans11-Conjugated Linoleic Acid Biosynthesis in Bovine Mammary Gland Epithelial Cells

by
Liying Peng
1,†,
Ge Bai
1,2,†,
Chunzheng Wang
1,
Jianan Dong
1,
Yongjun Liu
1,3,
Zhe Sun
1,4,5,
Yuguo Zhen
1,4,
Guixin Qin
1,
Xuefeng Zhang
1,4,
Natnael Demelash
1,6 and
Tao Wang
1,4,*
1
JLAU-Borui Dairy Science and Technology R&D Center, Key Laboratory of Animal Nutrition and Feed Science of Jilin Province, Key Laboratory of Animal Production Product Quality and Security Ministry of Education, College of Animal Science and Technology, Jilin Agricultural University, Changchun 130118, China
2
Changchun Institute of Biological Products Co., Ltd., Changchun 130012, China
3
Ningxia Agricultural Reclamation Helan Mountain Dairy Co., Ltd., Yinchuan 750028, China
4
Postdoctoral Scientific Research Workstation, Feed Engineering Technology Research Center of Jilin Province, Changchun Borui Science & Technology Co., Ltd., Changchun 130118, China
5
College of Life Science, Jilin Agricultural University, Changchun 130118, China
6
College of Agriculture and Environmental Science, Dilla University, Dilla P.O. Box 419, Ethiopia
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Animals 2022, 12(13), 1718; https://0-doi-org.brum.beds.ac.uk/10.3390/ani12131718
Submission received: 9 April 2022 / Revised: 28 June 2022 / Accepted: 29 June 2022 / Published: 2 July 2022
(This article belongs to the Special Issue The Mammary Gland Identifies Genes and Networks in Dairy Livestock)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Cis9, trans11-conjugated linoleic acid (c9, t11-CLA) is a bioactive compound that is mainly found in ruminant products and has a variety of beneficial biological functions in human health. However, the gene network of c9, t11-CLA biosynthesis in the mammary gland is still not well documented. Therefore, this study aimed to elucidate the stearoyl-coenzyme A desaturase (SCD1)-dependent gene network of c9, t11-CLA biosynthesis in MAC-T cells from the perspective of energy metabolism. The results showed that the deficiency of stearoyl-coenzyme A desaturase (SCD1), led by CAY10566, blocked the biosynthesis of c9, t11-CLA. In total, 60 SCD1-related proteins that were mainly involved in energy metabolism pathways were primary screened by Tandem mass tag-based quantitative proteomics analysis. Moreover, 17 proteins were validated under the parallel reaction monitoring analysis. Then, 11 genes involved in energy metabolism pathways were verified to have negative relationships with SCD1 after small RNA interference analysis. Based on the above results, we concluded that the 11 genes involved in energy metabolism pathways had an impact on the SCD1-dependent molecular mechanism of c9, t11-CLA biosynthesis. This study provided a fresh understanding of the gene network of c9, t11-CLA biosynthesis in the mammary gland from a metabolism perspective.

Abstract

The objective of the study was to elucidate the stearoyl-coenzyme A desaturase (SCD1)-dependent gene network of c9, t11-CLA biosynthesis in MAC-T cells from an energy metabolism perspective. The cells were divided into the CAY group (firstly incubated with CAY10566, a chemical inhibitor of SCD1, then incubated with trans-11-octadecenoic acid, (TVA)), the TVA group (only TVA), and the control group (without CAY, TVA). The c9, t11-CLA, and TVA contents were determined by gas chromatography. The mRNA levels of SCD1 and candidate genes were analyzed via real-time PCR. Tandem mass tag (TMT)-based quantitative proteomics, bioinformatic analysis, parallel reaction monitoring (PRM), and small RNA interference were used to explore genes involved in the SCD1-dependent c9, t11-CLA biosynthesis. The results showed that the SCD1 deficiency led by CAY10566 blocked the biosynthesis of c9, t11-CLA. In total, 60 SCD1-related proteins mainly involved in energy metabolism pathways were primarily screened by TMT-based quantitative proteomics analysis. Moreover, 17 proteins were validated using PRM analysis. Then, 11 genes were verified to have negative relationships with SCD1 after the small RNA interference analysis. Based on the above results, we concluded that genes involved in energy metabolism pathways have an impact on the SCD1-dependent molecular mechanism of c9, t11-CLA biosynthesis.

1. Introduction

Cis9, trans11-conjugated linoleic acid (c9, t11-CLA) is a bioactive compound with a variety of beneficial biological functions, including an important role in improving the nutritional value of milk [1,2,3]. For dairy cows, the main source of c9, t11-CLA in milk is de novo biosynthesis in the mammary glands and the remaining portion is synthesized by rumen microorganisms [4,5,6]. The de novo biosynthesis of c9, t11-CLA in the mammary glands refers to the c9, t11-CLA dehydrogenated from trans-11 vaccenic acid (TVA) by stearoyl-coenzyme A desaturase (SCD1) [7,8,9]. TVA is synthesized by rumen microorganisms and is subsequently absorbed and transported to the mammary tissue for c9, t11-CLA biosynthesis. The rumen microbial synthesis of c9, t11-CLA refers to the transformation of dietary linoleic acid into c9, t11-CLA under the isomerization by microorganisms in the rumen [10]. TVA and SCD1 are two essential components of the biosynthesis process of c9, t11-CLA. To date, five SCD isoforms have been discovered [11]. SCD1 is located in the endoplasmic reticulum and mainly expressed in the mammary glands, liver, and adipose tissue [12,13]. SCD1 is crucial for the c9, t11-CLA content in milk [14,15,16,17]. Factors that regulate the SCD1 activity or the synthesis of TVA in rumen can alter the milk content of cis9, trans11-CLA [4]. However, the SCD1-dependent gene network of c9, t11-CLA biosynthesis in mammary gland is still not well elucidated. It is important for more studies to examine the molecular mechanism of c9, t11-CLA synthesis in order to reveal the individual difference in c9, t11-CLA synthesis ability of dairy cows. This seems to be a potential option that would improve the nutritive value of milk.
It is well known that pentose phosphate pathway (PPP) is one of the ubiquitous catabolic pathways of glucose, and the reductase II (NADPH) produced in PPP is important for the desaturation of fatty acids [18]. Therefore, we proposed a hypothesis that the genes involved in energy metabolism pathways may have an impact on the SCD1-dependent molecular mechanism of c9, t11-CLA biosynthesis. The objective of the study was to elucidate the SCD1-dependent gene network of c9, t11-CLA biosynthesis from the perspective of energy metabolism. MAC-T cells, an immortalized bovine mammary gland epithelial cell line, were used in this study [19,20,21]. The cells were divided into the CAY group (firstly incubated with CAY10566, a chemical inhibitor of SCD1 then incubated with trans-11-octadecenoic acid (TVA)), the TVA group (only TVA), and control group (without CAY, TVA). The c9, t11-CLA and TVA contents were determined by gas chromatography. The mRNA levels of SCD1 and candidate genes were analyzed using real-time PCR. Tandem mass tag (TMT)-based quantitative proteomics, bioinformatic analysis, parallel reaction monitoring (PRM), and small RNA interference were used to explore the genes involved in the c9, t11-CLA biosynthesis. This study may provide a fresh understanding of the SCD1-dependent gene network of c9, t11-CLA biosynthesis in the mammary gland.

2. Materials and Methods

2.1. Materials

Total RNA was extracted from MAC-T cells and complementary DNA (cDNA) was generated using commercial kits (TIANGEN Biotech Co., Ltd., Beijing, China). The kit for real-time polymerase chain reaction (real-time PCR) was from Foregene Co., Ltd., Beijing, China. CAY10566 (Cayman Chemical Company, Ann Arbor, MI, USA), a chemical inhibitor of SCD1, was used for the construction of SCD1-inhibited MAC-T cells model [20]. Trans-11-octadecenoic acid (Sigma-Aldrich, St. Louis, MO, USA), fatty acid standard 37 component FAME mix (FAME37), and trans-11-octadecenoic methyl ester (NU-CHEKPREP, Elysian, MN, USA), and methyl 9(Z), 11(E)-octadecadienoate (Matreya LLC, State College, PA, USA) were purchased for the fatty acids’ analysis.

2.2. Cell Treatment and Real-Time PCR Analysis

MAC-T cells were maintained in DMEM/High glucose medium, and the culture conditions were similar to those described previously [22]. MAC-T cells were divided into three groups: CAY group (incubated with 50 µM TVA for 2.5 h after the treatment of CAY10566 at concentrations of 10 nM for 12 h), which was considered as the SCD1-inhibited MAC-T cells model, TVA group (only incubated with 50 µM TVA for 2.5 h), and control group (without CAY, TVA). CAY10566 was proved to inhibit SCD1 at the transcriptional level in our previous works [20,23]. RNA was reverse transcribed after measuring the concentration using a NanoDrop 2000 ultra-micro spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). After reverse transcription, the final assay was carried out in the real-time fluorescence quantification system (Bio-Rad, Hercules, CA, USA). The primer sequences are shown in Supplementary Table S1 [23]. The β-actin was used for a threshold comparison with the target gene for final relative quantification [24].

2.3. Fatty Acids Analysis

The total lipids of the cells in those three groups (CAY group, TVA group, control group; three replicates in each group) were extracted after adding chloroform/methanol (2:1, v/v), and the specific steps were as follows [25,26,27]: 10 mL of chloroform/methanol was added to the cells, shaking and ultrasound were applied to break the cells for lipid extraction, and excess pure water was added to wash the cells. Subsequently, a centrifugation at 2500× g for 5 min was performed to take the organic layer; 10 mol/L KOH solution and methanol solution were added for methyl esterification, after which 12 mol/L H2SO4 was added to samples for methyl esterification. Finally, hexane (3 mL) was added, followed by centrifugation for 2000× g for 10 min. The supernatant (fatty-acid methyl esters (FAMEs)) was subjected to gas chromatography (GC) analysis. FAMEs were analyzed using an Agilent 7890B GC system (Agilent Technologies, Santa Clara, CA, USA) equipped with a DB-Fast FAME fused silica capillary column (90 m × 0.32 mm × 0.2 µm). Front inlet and FID temperatures were maintained at 240 °C and 260 °C, respectively. The hydrogen flow rate to the detector was 30 mL/min, and the airflow rate was 400 mL/min. The split ratio was 10:1 [28,29]. Percentages of individual fatty acids were calculated as the ratio of individual GC areas to the total area of identified fatty acids [27].

2.4. Identification of Differentially Expressed Proteins by Tandem Mass Tag (TMT)-Based Quantitative Proteomics

The protein samples of those three groups (CAY group, TVA group, control group; three replicates in each group) were prepared as follows: a 4-fold volume of lysis buffer (8 mol/L urea, 1% protease inhibitor) was used to crack the cells [30]. To obtain the protein sample, the supernatant was collected through centrifuging at 12,000× g for 10 min at 4 °C then the protein concentration was measured using a BCA kit, taking an equal amount of protein for enzyme digestion. TCA was slowly added to the sample, then mixed to precipitate for 2 h at 4 °C. The sample was centrifuged at 4500× g for 5 min, and the precipitate was washed with pre-cooled acetone after discarding the supernatant. Dried, TEAB was then added after the precipitate was dried and broken up by sonication and then trypsin was added overnight at a ratio of 1:50 (protease:protein). Dithiothreitol (DTT) and iodoacetamide (IAA) were individually added to a final concentration of 5 mM and 11 mM, respectively [31,32]. The trypsinized peptides were desalted using a Strata X C18 column and freeze-dried [33].
Then, the peptides were dissolved in 0.5M TEAB, and the peptide was labeled according to the instructions of the tandem mass tags (TMT) kit. The peptides were graded by high-pH, reversed-phase HPLC using an Agilent 300 Extend column C18 (5 μm × 4.6 mm × 250 mm), separated combined, and then freeze-dried under vacuum for subsequent usage [31,34]. The peptides were dissolved in liquid chromatography mobile phase A (aqueous solution of 0.1% HCOOH + 2% acetonitrile) and separated using an ultra-performance liquid chromatography system. Mobile phase B contained 0.1% HCOOH and 90% acetonitrile [35,36]. The liquid phase gradient was set as follows: 0–38 min, 6–22% B; 38–52 min, 22–32% B; 52–56 min, 32–80% B; 56–60 min, 80% B. Flow rate was maintained at 500 nL/min; ion source voltage: 2.1 kV, primary mass spectrometry scan range set to 350–1600 m/z and the scan resolution was set to 120,000; the secondary mass spectrometry scan range was fixed at 100 m/z and the secondary scan resolution was set to 30,000 [37]. The mass spectrometry proteomics data were deposited to the ProteomeXchange Consortium via the PRIDE [38] partner repository with the dataset identifier PXD028620.

2.5. Protein Functional Annotation and Bioinformatic Analysis

Gene ontology (GO) was used to describe the properties of genes and gene products in living organisms, and was divided into biological process, cellular component, and molecular function [39,40]. GO annotation analysis was performed using the UniProt-GOA database (http://www.ebi.ac.uk/GOA/, accessed on 28 December 2020) and the protein sequence-based algorithm software InterProScan. The protein structural domain annotation was performed using the InterProScan software and the corresponding InterPro structural domain database (http://www.ebi.ac.uk/interpro/, accessed on 28 December 2020) for the protein domain annotation of the identified proteins. KEGG pathway annotation was performed on the submitted proteins using the online service tool KAAS.

2.6. Parallel Reaction Monitoring (PRM) Validation

To assess the reliability of proteomic differential protein data, parallel reaction monitoring (PRM) was used to verify the differential proteins related to SCD1. PRM is a method used for the selective detection and quantification of target proteins and peptides using ion monitoring technology [41,42].

2.7. Small RNA Interference

MAC-T cells were transfected with specific siRNA for bovine 6-phosphogluconolactonase (PGLS), phosphofructokinase, liver type (PFKL), aldolase, fructose-bisphosphate A (ALDOA), triosephosphate isomerase 1 (TPI1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), phosphoglycerate kinase 1 (PGK1), phosphoglycerate mutase 1 (PGAM1), enolase 1 (ENO1), pyruvate kinase M1/2 (PKM), lactate dehydrogenase A (LDHA), and lactate dehydrogenase B (LDHB), or a scrambled siRNA using Lipofectamine® RNAIMAZ Reagent (Invitrogen, CA, USA). The specific siRNAs are shown in Supplementary Table S2 [20]. After 24 h, the cells were incubated with 50 µM trans-11 C18:1 for 2.5 h. Then, cells were harvested for total mRNA extraction and real-time PCR analysis [20]. The reaction system of real-time PCR and the fluorescence quantitative reaction procedure are listed in Supplementary Tables S3 and S4.

2.8. Statistical Analysis

All data are presented as the mean ± standard error (S.E.). The data of mRNA relative expression, TVA and c9, t11-CLA level were analyzed with the one-way analysis of variance (one-way ANOVA, SPSS 24.0, SPSS Inc., Chicago, IL, USA). Differences were considered significant at p < 0.05, highly significant at p < 0.01. Differential proteins were screened from quantifiable proteins with p-value < 0.05, and a differential expression change more than 1.3/1 was defined as the threshold of change for significant up-regulation, and a change of less than 1/1.3 as the threshold of change for significant down-regulation. TMT labeling and PRM verification results were compared using the GraphPad Prism 8.0.2 (GraphPad Software, San Diego, CA, USA).

3. Results

3.1. Effects of CAY and TVA on SCD1 mRNA Expression

The results (Figure 1A) indicated that compared with the control and TVA group, the mRNA relative expression of SCD1 was significantly decreased after CAY treatment (p < 0.01), and the inhibitory effect was about 90%. After adding TVA, SCD1 activity was promoted, and its mRNA expression was significantly higher than that in the control group (p < 0.05).

3.2. Effects of CAY and TVA on TVA Accumulation and c9, t11-CLA Synthesis

After the addition of TVA and CAY, TVA accumulation in the CAY group was significantly higher than that in the TVA group (Figure 1B; p < 0.01), while the c9, t11-CLA biosynthesis was significantly lower than that in the TVA group (Figure 1C; p < 0.01).

3.3. Identification of Differentially Expressed Proteins

A total of 312,032 secondary spectra were generated using mass spectrometry analysis of the obtained protein samples and 6455 quantifiable proteins were obtained (Figure 2A). The two-way comparison between the control, TVA, and CAY groups regarding protein expression is summarized in Figure 2B. A volcanic map of these differentially expressed proteins (red dots indicate up-regulated proteins and blue dots indicate down-regulated proteins) in MAC-T cells under different incubation conditions is shown in Figure 3.

3.4. Functional Annotation of Differentially Expressed Proteins

The differential proteins identified in each group were mainly located in the cytoplasm, nucleus, extracellular space, mitochondria, and cell membrane, with the highest number of proteins being localized in the cytoplasm and nucleus (Figure 4A). Differential proteins were found to be involved in a total of 28 terms of biological functions, mainly focusing on biological regulation, metabolic processes, binding, and catalytic activity (Figure 4B).

3.5. GO and KEGG Analysis of Differentially Expressed Proteins

Following GO classification with biological process (BP), cellular component (CC), and molecular function (MF), a cluster analysis was conducted to compare the functional correlations between DE proteins in the experimental groups; the results are displayed as a heat map (Figure 5).
The up-regulated and down-regulated differential proteins of the two comparison groups were enriched by about 30 terms under the three categories, and we mainly focused on the pathways related to energy metabolism. Among the secondary pathways, the glycolytic pathway and nucleotide catabolic pathway of the up-regulated proteins were found in the CAY/control group, and the functional pathways of lipid acyl coenzyme A biological processes, fatty acid elongase activity, and fatty acid synthase activity of the down-regulated proteins were associated with fatty acids; in the CAY/TVA group, the glycolytic pathway and nucleotide catabolic pathway were equally enriched in the up-regulated differential proteins, as well as in the glycolytic pathway and nucleotide catabolic pathway. Regarding the down-regulated differential protein enrichment pathways, lipoprotein catabolism and sphingolipid catabolism were associated with fatty acids. The enrichment process of the TVA/control group did not have the desired related pathways. KEGG pathway enrichment statistical analysis was performed for up-regulated and down-regulated differential proteins (Figure 6), similar to the GO analysis. Glycolysis, pentose phosphate, and amino acid and nucleotide energy metabolism pathways, fatty acid extension pathways, and unsaturated fatty acid biosynthesis pathways were enriched. In addition, we determined the common pathways from GO and KEGG analyses, and the results are shown in Supplementary Tables S5–S8.

3.6. PRM Validation

Among the differential proteins screened from the proteomics analysis, 17 were selected for PRM validation, and a comparison between the verification results and the TMT labeling results is shown in Figure 7. The ploidy changes are based on 1; a score greater than 1 represents up-regulated proteins and a score of less than 1 represents down-regulated proteins. PRM validation results and proteomics trends were consistent among the groups, except for PKM and PPT1 in the TVA/control group, indicating the reliability of the results of this experiment.

3.7. The Responses of SCD1 after the Knockdown of PGLS, PFKL, ALDOA, TPI1, GAPDH, PGK1, PGAM1, ENO1, PKM, LDHA, and LDHB

The results (Figure 8) indicated that PGLS, PFKL, ALDOA, TPI1, GAPDH, PGK1, PGAM1, ENO1, PKM, LDHA, and LDHB were successfully and individually reduced by their specific siRNA sequences (p < 0.01).
After the knockdown of PGLS (siRNA-1, siRNA-2, siRNA-3), PFKL (siRNA-1, siRNA-2, siRNA-3), ALDOA (siRNA-1, siRNA-2, siRNA-3), TPI1 (siRNA-1, siRNA-2, siRNA-3), GAPDH (siRNA-1, siRNA-2, siRNA-3), PGK (siRNA-1, siRNA-2, siRNA-3), PGAM1 (siRNA-1), ENO1 (siRNA-1, siRNA-2, siRNA-3), PKM (siRNA-2, siRNA-3), LDHA (siRNA-1, siRNA-2, siRNA-3), or LDHB (siRNA-2, siRNA-3), the mRNA expression of SCD1 was significantly elevated (p < 0.05; Figure 9).
Therefore, based on the above data the SCD1-dependent gene network of c9, t11-CLA biosynthesis is outlined for the mammary gland in Figure 10, from the energy metabolism perspective.

4. Discussion

As a key enzyme in fatty acid biosynthesis, SCD1 is highly expressed in adipose tissue, liver, and especially in lactating mammary tissue, which plays an important role in the biosynthesis of monounsaturated fatty acids [43,44]. The inhibition of SCD1 expression causes a series of physiological changes in the body related to fatty acid oxidation, lipid biosynthesis, and triacylglycerol (TAG) deficiency [45,46]. In this study, the contents of c9, t11-CLA, and TVA were determined by gas chromatography after the incubation with CAY10566, while the mRNA levels of SCD1 were analyzed by real-time PCR. Our results showed that CAY10566 successfully inhibited SCD1, which could not catalyze the conversion of TVA to c9, t11-CLA and is consistent with a previous study by our group [23]. The MAC-T cells model with loss of SCD1 function was successfully constructed, and also showed that SCD1 plays a critical role in regulating the c9, t11-CLA content. Proteins are essential to cell activity, and their characterization is extremely important for revealing life phenomena and understanding life activities; therefore, research at the protein level is crucial [47]. SCD1 has been directly or indirectly confirmed in several studies to be a desaturase affecting the c9, t11-CLA content in dairy cows and goats, acting together with other proteins for fatty acid biosynthesis [5,15,25,48].
After verifying the function of SCD1 in regulating c9, 11-CLA, using more accurate and specific omics techniques to search SCD1-related proteins is necessary for the better understanding of the gene network of c9, t11-CLA biosynthesis in the mammary glands of dairy cows. Therefore, in this study, tandem mass tag (TMT)-based quantitative proteomics, bioinformatic analysis, parallel reaction monitoring (PRM), and small RNA interference were used to explore the genes involved in the SCD1-dependent c9, t11-CLA biosynthesis. Through the tandem mass tag (TMT)-based quantitative proteomics and bioinformatic analysis, 60 differentially expressed proteins that were mainly enriched in the pathways of energy metabolism were preliminary screened in this study. Energy metabolism and lipid metabolism are inextricably linked, glycolysis-producing dihydroxyacetone phosphate reduction forms of glycerol, pyruvate oxidative decarboxylation forms acetyl coenzyme A, a raw material for fatty acid biosynthesis [49]. In addition, the pentose phosphate pathway is also closely related to fatty acid formation, and the reductase II (NADPH) produced in the pentose phosphate pathway is used in the biosynthesis of fatty acids [18,50,51,52]. Nucleotides, such as NAD+, NADP+, AD, and CoA, can regulate many biological processes and participate in the biosynthesis of coenzymes. Several important coenzymes of nucleotide coenzymes are adenosine derivatives. NAD+ and FAD can participate in many redox reactions through the transfer of hydrogen atoms, while CoA activates fatty acids, which are related to fatty acid oxidation and steroid biosynthesis [53].
To check the reliability of the proteomics analysis results, a comparison was made between the PRM validation results and the TMT labeling results. In total, 17 proteins were validated using this comparison. The consistency of the PRM validation results and proteomics trends indicates the reliability of this experiment’s results. Then, these 17 proteins underwent individual knockdown by their specific siRNA sequences in MAC-T cells. After the siRNA interference, PGLS, PFKL, ALDOA, TPI1, GAPDH, PGK1, PGAM1, ENO1, PKM, LDHA, and LDHB were verified to have a negative relationship with SCD1. PGLS, a hydrolase, specifically catalyzes the hydrolysis from 6-phosphogluconolactone to 6-phosphogluconic acid [54,55]. PGLS plays an important role in the production of metabolites 5-phosphoribose (R5P) and NADPH in the pentose phosphate pathway [55]. PFKL is the most important rate-limiting enzyme in glycolysis [56,57,58]. ALDOA is a glycolytic enzyme that catalyzes the reversible conversion of fructose-1,6-bisphosphate to glyceraldehyde 3-phosphate and dihydroxyacetone phosphate [59]. TPI1 catalyzes the mutual transformation of dihydroxyacetone phosphate (DHAP) and d-glyceraldehyde-3-phosphate (G3P) during glycolysis and gluconeogenesis. GAPDH catalyzes the phosphorylation and oxidation of glyceraldehyde-3-phosphate at the same time, and NAD+ is used as the electron acceptor to generate 1,3-diphosphoglyceric acid, thus generating NADH [60]. PGK1 plays a role in the glycolytic pathway, catalyzing the transfer of high-energy phosphate from high-energy phosphate from the 1-position of 1,3-biphosphoglycerate (1,3-BPG) to ADP, thus generating of 3-phosphoglycerate (3-PG) and ATP [61]. PGAM1 is a glycolytic enzyme, which can catalyze the mutual transformation of 3-phosphoglycerate (3PG) and 2-phosphoglycerate (2PG) and 2,3-bisphosphoglycerate (2,3-BPG) as a cofactor in glycolysis [62]. As a moonlight protein, ENO1 performs multiple biochemical functions, such as catalyzing the conversion of 2-phosphoglycerate (2-PGA) to phosphoenolpyruvic acid (PEP) in the glycolytic pathway [63]. PKM is a subtype of pyruvate kinase (PK) which is known as the glycolytic enzyme that catalyzes the final step in glycolysis [64]. The enzyme LDH is a major player in the glucose metabolism, which can catalyze pyruvate to lactic acid [65].
These proteins play a key role in the SCD1-dependent gene network of c9, t11-CLA biosynthesis. Through the above analysis and a literature review, the 11 genes that are mainly involved in the energy metabolism pathways found in this study have close links with the SCD1-dependent molecular mechanism of c9, t11-CLA biosynthesis. This study not only provides a theoretical basis and the technical means for further elucidating the molecular mechanism of c9, t11-CLA synthesis in mammary gland of dairy cows, it also has important clinical implications for revealing the individual difference in c9, t11-CLA synthesis ability of dairy cows.

5. Conclusions

In total, 60 SCD1-related proteins mainly involved in energy metabolism pathways were screened in MAC-T cells using the proteomics approaches in this study. Additionally, 17 proteins were validated under PRM analysis. Then, PGLS, PFKL, ALDOA, TPI1, GAPDH, PGK1, PGAM1, ENO1, PKM, LDHA, and LDHB were verified to have a negative relationship with SCD1 after the small RNA interference analysis. Subsequently, the SCD1-dependent gene network of c9, t11-CLA biosynthesis was outlined, which provided a fresh understanding of this process from the perspective of energy metabolism.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ani12131718/s1, Table S1: Primer sequences for real-time polymerase chain reaction; Table S2: The sequences of scrambled and specific siRNA; Table S3: Reaction system of real-time PCR; Table S4: Fluorescence quantitative reaction procedure; Table S5: The proteins involved in glycolysis pathway screened from the tandem mass tag-based quantitative proteomics analysis; Table S6: The proteins involved in nucleotide catabolism pathway screened from the tandem mass tag-based quantitative proteomics analysis; Table S7: The proteins involved in other energy metabolism pathway screened from the tandem mass tag-based quantitative proteomics analysis; Table S8: The proteins involved in fat metabolism pathway screened from the tandem mass tag-based quantitative proteomics analysis.

Author Contributions

Conceptualization, T.W.; data curation, L.P. and G.B.; formal analysis, L.P., G.B., C.W. and J.D.; investigation, L.P. and G.B.; methodology, Z.S.; project administration, Y.Z.; resources, Y.L. and X.Z.; supervision, T.W.; writing—original draft, L.P. and G.B.; writing—review and editing, G.Q., N.D. and T.W. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by the National Natural Science Foundation of China (grant number 31902185).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The MAC-T cells used in this study were donated by Jianxin Liu and Hongyun Liu, Zhejiang University. We thanks a lot for their helps.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lavillonniere, F.; Chajes, V.; Martin, J.C.; Sébédio, J.L.; Lhuillery, C.; Bougnoux, P. Dietary purified cis9, trans11 conjugated linoleic acid isomer has anticarcinogenic properties in chemically induced mammary tumors in rats. Nutr. Cancer 2003, 45, 190–194. [Google Scholar] [CrossRef] [PubMed]
  2. Huebners, S.M.; Olson, J.M.; Campbell, J.P.; Bishop, J.W.; Crump, P.M.; Cook, M.E. Low dietary c9t11-conjugated linoleic acid intake from dairy fat or supplements reduces inflammation in collagen-induced arthritis. Lipids 2016, 51, 807–819. [Google Scholar] [CrossRef] [PubMed]
  3. Schaftenaa, F.; Frodermann, V.; Kuiper, J.; Lutgens, E. Atherosclerosis: The interplay between lipids and immune cells. Curr. Opin. Lipidol. 2016, 27, 209–215. [Google Scholar] [CrossRef]
  4. Wang, T.; Lee, H.G. Advances in Research on cis9, trans11 conjugated linoleic acid: A major functional conjugated linoleic acid isomer. Crit. Rev. Food Sci. 2015, 55, 720–731. [Google Scholar] [CrossRef] [PubMed]
  5. Yamaguchi, M.; Naito, T.; Nagao, Y.; Kabuyama, Y.; Hashimoto, K.; Azuma, N.; Maeda, I. Effect of increased feeding of dietary α-linolenic acid by grazing on formation of the cis9, trans11-18:2 isoform of conjugated linoleic acid in bovine milk. Anim. Sci. J. 2017, 88, 1006–1011. [Google Scholar] [CrossRef] [PubMed]
  6. Robyn, B.; Robyn, F.; Orina, B. Atheroprotective effects of conjugated linoleic acid. Br. J. Clin. Pharmacol. 2017, 83, 46–53. [Google Scholar] [CrossRef] [Green Version]
  7. Mosley, E.E.; Shafii, B.; Moate, P.J.; McGuire, M.A. cis9, trans11 Conjugated linoleic acid is synthesized directly from vaccenic acid in lactating dairy cattle. J. Nutr. 2006, 136, 570–575. [Google Scholar] [CrossRef] [Green Version]
  8. Griinari, J.M.; Corl, B.A.; Lacy, S.H.; Chouinard, P.Y.; Nurmela, K.V.; Bauman, D.E. Conjugated linoleic acid is synthesized endogenously in lactating dairy cows by D9-desaturase. J. Nutr. 2000, 130, 2285–2291. [Google Scholar] [CrossRef]
  9. Kay, J.K.; Mackle, T.R.; Auldist, M.J.; Thomson, N.A.; Bauman, D.E. Endogenous synthesis of cis9, trans11 conjugated linoleic acid in dairy cows fed fresh pasture. J. Dairy Sci. 2004, 87, 369–378. [Google Scholar] [CrossRef] [Green Version]
  10. Sun, J.H.; Hyeong, S.K.; Young, Y.B.; Yeonhwa, P. Overview of conjugated linoleic acid formation and accumulation in animal products. Livest. Sci. 2016, 195, 105–111. [Google Scholar] [CrossRef]
  11. Iwai, T.; Kume, S.; Chin-Kanasaki, M.; Kuwagata, S.; Araki, H.; Takeda, N.; Sugaya, T.; Uzu, T.; Maegawa, H.; Araki, S.I. Stearoyl-CoA desaturase-1 protects cells against lipotoxicity-mediated apoptosis in proximal tubular cells. Int. J. Mol. Sci. 2016, 17, 1868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Lengi, A.J.; Corl, B.A. Identification and characterization of a novel bovine stearoyl-CoA desaturase isoform with homology to human SCD5. Lipids 2007, 42, 499–508. [Google Scholar] [CrossRef] [PubMed]
  13. Rezamand, P.; Watts, J.S.; Yavah, K.M.; Mosley, E.E.; Ma, L.Y.; Corl, B.A.; McGuire, M.A. Relationship between stearoyl-CoA desaturase 1 gene expression, relative protein abundance, and its fatty acid products in bovine tissues. J. Dairy Res. 2014, 81, 333–339. [Google Scholar] [CrossRef] [PubMed]
  14. Ntambi, J.M.; Buhrow, S.A.; Kaestner, K.H.; Christy, R.J.; Sibley, E.; KellyJr, T.J.; Lane, M.D. Differentiation-induced gene expression in 3T3-L1 preadipocytes. Characterization of a differentially expressed gene encoding stearoyl-CoA desaturase. J. Biol. Chem. 1988, 263, 17291–17300. [Google Scholar] [CrossRef]
  15. Wang, L.Z.; You, J.H.; Zhong, B.S.; Ren, C.F.; Zhang, Y.L.; Li, M.; Zhang, G.M.; Jia, R.X.; Ying, S.J.; Wang, F. Scd1 mammary-specific vector constructed and overexpressed in goat fibroblast cells resulting in an increase of palmitoleic acid and oleic acid. Biochem. Biophys. Res. Commun. 2013, 443, 389–394. [Google Scholar] [CrossRef]
  16. Yao, D.W.; Luo, J.; He, Q.Y.; Shi, H.B.; Li, J.; Wang, H.; Xu, H.F.; Chen, Z.; Yi, Y.Q.; Loor, J.J. SCD1 alters long-chain fatty acid (LCFA) composition and its expression is directly regulated by SREBP-1 and PPARgamma 1 in dairy goat mammary cells. J. Cell Physiol. 2017, 232, 635–649. [Google Scholar] [CrossRef]
  17. Lounis, M.A.; Escoula, Q.; Veillette, C.; Bergeron, K.F.; Ntambi, J.M.; Mounier, C. SCD1 deficiency protects mice against ethanol-induced liver injury. BBA-Mol. Cell Biol. Lipids 2016, 1861, 1662–1670. [Google Scholar] [CrossRef] [Green Version]
  18. Paton, C.M.; Ntambi, J.M. Biochemical and physiological function of stearoyl-CoA desaturase. Am. J. Physiol.-Endoc. Metab. 2009, 297, E28–E37. [Google Scholar] [CrossRef] [Green Version]
  19. Jedrzejczak, M.; Szatkowska, I. Bovine mammary epithelial cell cultures for the study of mammary gland functions. In Vitro Cell. Dev. Biol. Anim. 2014, 50, 389–398. [Google Scholar] [CrossRef] [Green Version]
  20. Wang, T.; Lee, S.B.; Hwang, J.H.; Lim, J.N.; Jung, U.S.; Kim, M.J.; Kim, H.S.; Kang, S.H.; Choi, J.S.; Lee, S.G.; et al. Proteomic analysis reveals PGAM1 altering cis9, trans11 conjugated linoleic acid synthesis in bovine mammary gland. Lipids 2015, 50, 469–481. [Google Scholar] [CrossRef]
  21. Mitz, C.A.; Viloria-Petit, A.M. Contrasting effects of transforming growth factor β1 on programmed cell death of bovine mammary epithelial cell lines MAC-T and BME-UV1. J. Dairy Sci. 2020, 103, 5532–5549. [Google Scholar] [CrossRef] [PubMed]
  22. Wang, T.; Lim, J.N.; Bok, J.D.; Kim, J.H.; Kang, S.K.; Lee, S.B.; Hwang, J.H.; Lee, K.H.; Kang, H.S.; Choi, Y.J.; et al. Association of protein expression in isolated milk epithelial cells and cis9, trans11 conjugated linoleic acid proportions in milk from dairy cow. J. Sci. Food Agric. 2013, 94, 1835–1843. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, T.; Lee, H.G. Responses of MAC-T cells to inhibited stearoyl-CoA desaturase 1 during cis9, trans11 conjugated linoleic acid synthesis. Lipids 2018, 53, 647–652. [Google Scholar] [CrossRef] [PubMed]
  24. Johnson, T.L.; Fujimoto, B.A.; Jiménez-Flores, R.; Peterson, D.G. Growth hormone alters lipid composition and increases the abundance of casein and lactalbumin mRNA in the MAC-T cell line. J. Dairy Res. 2010, 77, 199–204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Bernard, L.; Leroux, C.; Bonnet, M.; Rouel, J.; Martin, P.; Chilliard, Y. Expression and nutritional regulation of lipogenic genes in mammary gland and adipose tissues of lactating goats. J. Dairy Res. 2005, 72, 250–255. [Google Scholar] [CrossRef]
  26. O’Fallon, J.V.; Busboom, J.R.; Nelson, M.L.; Gaskins, C.T. A direct method for fatty acid methyl ester synthesis: Application to wet meat tissues, oils, and feedstuffs. J. Anim. Sci. 2007, 85, 1511–1521. [Google Scholar] [CrossRef] [Green Version]
  27. Wang, T.; Lim, J.N.; Choi, S.H.; Kang, H.S.; Lee, H.G. Comparative studies on derivatization methods of single or mixed fatty acids. Food Sci. Biotechnol. 2013, 22, 1573–1579. [Google Scholar] [CrossRef]
  28. Wang, T.; Oh, J.J.; Lim, J.N.; Hong, J.E.; Kim, J.H.; Kim, J.H.; Kang, H.S.; Choi, Y.J.; Lee, H.G. Effects of lactation stage and individual performance on milk cis9, trans11 conjugated linoleic acids content in dairy cows. Asian-Australas. J. Anim. Sci. 2013, 26, 189–194. [Google Scholar] [CrossRef] [Green Version]
  29. Min, C.W.; Park, J.; Bae, J.W.; Agrawal, G.K.; Rakwal, R.; Kim, Y.; Yang, P.F.; Kim, S.T.; Gupta, R. In-depth investigation of low-abundance proteins in matured and filling stages seeds of glycine max employing a combination of protamine sulfate precipitation and TMT-based quantitative proteomic analysis. Cells 2020, 9, 1517. [Google Scholar] [CrossRef]
  30. Duan, Z.Q.; Yuan, C.; Han, Y.F.; Zhou, L.; Zhao, J.F.; Ruan, Y.; Chen, J.Q.; Ni, M.M.; Ji, X.Q. TMT-based quantitative proteomics analysis reveals the attenuated replication mechanism of Newcastle disease virus caused by nuclear localization signal mutation in viral matrix protein. Virulence 2020, 11, 607–635. [Google Scholar] [CrossRef]
  31. Zhao, L.; Geng, T.T.; Sun, K.C.; Su, S.; Zhao, Y.; Bao, N.; Pan, L.; Sun, H.; Li, M.Y. Proteomic analysis reveals the molecular mechanism of Hippophae rhamnoides polysaccharide intervention in LPS-induced inflammation of IPEC-J2 cells in piglets. Int. J. Biol. Macromol. 2020, 164, 3294–3304. [Google Scholar] [CrossRef] [PubMed]
  32. Mirzaei, M.; Pascovici, D.; Wu, J.X.; Chick, J.; Wu, Y.; Cooke, B.; Haynes, P.; Molloy, M.P. TMT one-stop shop: From reliable sample preparation to computational analysis platform. In Proteome Bioinformatics; Keerthikumar, S., Mathivanan, S., Eds.; Humana Press: New York, NY, USA, 2010; Volume 1549. [Google Scholar] [CrossRef]
  33. Ma, C.; Wang, W.W.; Wang, Y.D.; Sun, Y.; Kang, L.; Zhang, Q.; Jiang, Y.L. TMT-labeled quantitative proteomic analyses on the longissimus dorsi to identify the proteins underlying intramuscular fat content in pigs. J. Proteom. 2019, 213, 103630. [Google Scholar] [CrossRef] [PubMed]
  34. Hou, C.Y.; Guo, D.Q.; Yu, X.; Wang, S.Y.; Liu, T.H. TMT-based proteomics analysis of the anti-hepatocellular carcinoma effect of combined dihydroartemisinin and sorafenib. Biomed. Pharmacother. 2020, 126, 109862. [Google Scholar] [CrossRef]
  35. Sun, J.; Jiang, T.C.; Gu, F.; Ma, D.H.; Liang, J.M. TMT-based proteomic analysis of plasma from children with rolandic epilepsy. Dis. Markers 2020, 2020, 8840482. [Google Scholar] [CrossRef]
  36. Jia, B.Y.; Xiang, D.C.; Liu, S.N.; Zhang, B.; Shao, Q.Y.; Hong, Q.H.; Quan, G.B.; Wu, G.Q. TMT-based quantitative proteomic analysis of cumulus cells derived from vitrified porcine immature oocytes following in vitro maturation. Theriogenology 2020, 152, 8–17. [Google Scholar] [CrossRef]
  37. Wang, Z.K.; Liu, F.J.; Ye, S.L.; Jiang, P.; Yu, X.C.; Xu, J.; Du, X.; Ma, L.; Cao, H.G.; Yuan, C.; et al. Plasma proteome profiling of high-altitude polycythemia using TMT-based quantitative proteomics approach. J. Proteom. 2018, 194, 60–69. [Google Scholar] [CrossRef] [PubMed]
  38. Deutsch, E.W.; Bandeira, N.; Sharma, V.; Perez-Riverol, Y.; Carver, J.J.; Kundu, D.J.; García-Seisdedos, D.; Jarnuczak, A.F.; Hewapathirana, S.; Pullman, B.S.; et al. The ProteomeXchange consortium in 2020: Enabling ‘big data’ approaches in proteomics. Nucleic Acids Res. 2020, 48, D1145–D1152. [Google Scholar] [CrossRef] [Green Version]
  39. Kan, L.B.; Cui, D.Z.; Chai, Y.Y.; Ma, L.; Li, X.Y.; Zhao, M. TMT-based quantitative proteomic analysis of antitumor mechanism of Sporisorium reilianum polysaccharide WM-NP-60 against HCT116 cells. Int. J. Biol. Macromol. 2020, 165, 1755–1764. [Google Scholar] [CrossRef]
  40. Chen, C.; Hou, J.; Tanner, J.J.; Cheng, J.L. Bioinformatics methods for mass spectrometry-based proteomics data analysis. Int. J. Mol. Sci. 2020, 21, 2873. [Google Scholar] [CrossRef] [Green Version]
  41. Xue, T.; Liu, P.; Zhou, Y.; Liu, K.; Yang, L.; Moritz, R.L.; Yan, W.; Xu, L.X. Interleukin-6 induced “Acute” phenotypic microenvironment promotes Th1 anti-tumor immunity in cryo-thermal therapy revealed by shotgun and parallel peaction monitoring proteomics. Theranostics 2016, 6, 773–794. [Google Scholar] [CrossRef]
  42. Yan, X.Y.; Wu, Y.J.; Zhong, F.J.; Jiang, Q.Y.; Zhou, T.T.; Guo, Y.F.; Yang, X.R.; Liang, J.; Liao, D.J.; Lan, G.Q. iTRAQ and PRM-based quantitative proteomics in T2DM-susceptible and-tolerant models of Bama mini-pig. Gene 2018, 675, 119–127. [Google Scholar] [CrossRef] [PubMed]
  43. Hwang, J.; Singh, N.; Long, C.; Smith, S.B. The lentiviral system construction for high expressed porcine stearoyl CoA desaturase-l and functional characterization in stably transduced porcine swine kidney cells. Lipids 2018, 53, 933–945. [Google Scholar] [CrossRef] [PubMed]
  44. Mahesh, M.; Bharathi, M.; Reddy, M.R.G.; Kumar, M.S.; Putcha, U.K.; Vajreswari, A.; Jeyakumar, S.M. Carrot juice administration decreases liver stearoyl-CoA desaturase I and improves docosahexaenoic acid levels but not steatosis in high fructose diet-fed weanling wistar rats. Prev. Nutr. Food Sci. 2016, 21, 171–180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Cohen, P.; Miyazaki, M.; Socci, N.D.; Hagge-Greenberg, A.; Liedtke, W.; Soukas, A.A.; Sharma, R.; Hudgins, L.C.; Ntambi, J.M.; Friedman, J.M. Role for stearoyl-coa desaturase eptin-mediated weight loss. Science 2002, 297, 240–243. [Google Scholar] [CrossRef] [PubMed]
  46. Zou, Y.; Wang, Y.N.; Ma, H.; He, Z.H.; Tang, Y.; Guo, L.; Liu, Y.; Ding, M.; Qian, S.W.; Tang, Q.Q. SCD1 promotes lipid mobilization in subcutaneous white adipose tissue. J. Lipid Res. 2020, 61, 1589–1604. [Google Scholar] [CrossRef]
  47. Teng, Z.X.; Guo, M.Z.; Liu, X.Y.; Tian, Z.; Che, K. Revealing protein functions based on relationships of interacting proteins and GO terms. J. Biomed. Semant. 2017, 8, 27. [Google Scholar] [CrossRef] [Green Version]
  48. Conte, G.; Mele, M.; Chessa, S.; Castiglioni, B.; Serra, A.; Pagnacco, G.; Secchiari, P. Diacylglycerol acyltransferase 1, stearoyl-CoA desaturase 1, and sterol regulatory element binding protein 1 gene polymorphisms and milk fatty acid composition in Italian Brown cattle. J. Dairy Sci. 2010, 93, 753–763. [Google Scholar] [CrossRef] [Green Version]
  49. Yan, P.; Eng, O.C.; Yu, C.J. A review on the expression and metabolic features of orphan human cytochrome P450 2S1 (CYP2S1). Curr. Drug Metab. 2018, 19, 917–929. [Google Scholar] [CrossRef]
  50. Bray, G.A.; Champagne, C.M. Obesity and the metabolic syndrome: Implications for dietetics practitioners. J. Am. Diet. Assoc. 2004, 104, 86–89. [Google Scholar] [CrossRef]
  51. Ebert, K.; Witt, H. Fructose malabsorption. Mol. Cell. Pediatr. 2016, 3, 10. [Google Scholar] [CrossRef]
  52. Roy, K.; Wu, Y.; Meitzler, J.L.; Juhasz, A.; Liu, H.; Jiang, G.J.; Lu, J.M.; Antony, S.; Doroshow, J.H. NADPH oxidases and cancer. Clin. Sci. 2015, 128, 863–875. [Google Scholar] [CrossRef] [PubMed]
  53. Opielka, M.; Sobocki, B.; Mierzejewska, P.; Smolenski, R.T. The effect of trehalose on intracellular and extracellular nucleotide metabolism. A pilot study. Nucleosides Nucleotides Nucleic Acids 2020, 39, 1400–1409. [Google Scholar] [CrossRef] [PubMed]
  54. Beutler, E.; Kuhl, W.; Gelbart, T. 6-Phosphogluconolactonase deficiency, a hereditary erythrocyte enzyme deficiency: Possible interaction with glucose-6-phosphate dehydrogenase deficiency. Proc. Natl. Acad. Sci. USA 1985, 82, 3876–3878. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Li, C.Z.; Chen, J.; Li, Y.S.; Ye, Z.T.; Tian, X.B.; Wei, Y.; Hao, Z.C.; Pan, Y.; Zhou, H.L. 6-Phosphogluconolactonase promotes hepatocellular carcinogenesis by activating pentose phosphate pathway. Front. Cell Dev. Biol. 2021, 9, 753196. [Google Scholar] [CrossRef] [PubMed]
  56. Feng, Y.L.; Zhang, Y.; Cai, Y.; Liu, R.J.; Lu, M.L.; Li, T.Z.M.; Fu, Y.; Guo, M.; Huang, H.C.; Ou, Y.F.; et al. A20 targets PFKL and glycolysis to inhibit the progression of hepatocellular carcinoma. Cell Death Dis. 2020, 11, 89. [Google Scholar] [CrossRef]
  57. Mor, I.; Cheung, E.C.; Vousden, K.H. Control of glycolysis through regulation of PFK1: Old friends and recent additions. Cold Spring Harb. Symp. Quant. Biol. 2011, 76, 211–216. [Google Scholar] [CrossRef]
  58. Moreno-Sanchez, R.; Rodriguez-Enriquez, S.; Marin-Hernandez, A.; Saavedra, E. Energy metabolism in tumor cells. FEBS J. 2007, 274, 1393–1418. [Google Scholar] [CrossRef]
  59. Huang, Z.K.; Hua, Y.B.; Tian, Y.; Qin, C.; Qian, J.; Bao, M.L.; Liu, Y.Y.; Wang, S.Q.; Cao, Q.; Ju, X.B.; et al. High expression of fructose-bisphosphate aldolase A induces progression of renal cell carcinoma. Oncol. Rep. 2018, 39, 2996–3006. [Google Scholar] [CrossRef] [Green Version]
  60. Vladimir, F.L.; Irina, V.G.; Boris, A.M. Glyceraldehyde-3-phosphate dehydrogenase is a multifaceted therapeutic target. Pharmaceutics 2020, 12, 416. [Google Scholar] [CrossRef]
  61. Qian, X.; Li, X.J.; Cai, Q.S.; Zhang, C.B.; Yu, Q.J.; Jiang, Y.H.; Lee, J.H.; Hawke, D.; Wang, Y.G.; Xia, Y.; et al. Phosphoglycerate kinase 1 phosphorylates beclin1 to induce autophagy. Mol. Cell 2017, 65, 917–931. [Google Scholar] [CrossRef] [Green Version]
  62. Wang, P.H.; Jiang, L.L.; Cao, Y.; Ye, D.Y.; Zhou, L. The design and synthesis of N-xanthone benzenesulfonamides as novel phosphoglycerate mutase 1 (PGAM1) inhibitors. Molecules 2018, 23, 1396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Qiao, G.; Wu, A.G.; Chen, X.L.; Tian, Y.; Lin, X.K. Enolase 1, a moonlighting protein, as a potential target for cancer treatment. Int. J. Biol. Sci. 2021, 17, 3981–3992. [Google Scholar] [CrossRef] [PubMed]
  64. Taniguchi, K.; Sugito, N.; Shinohara, H.; Kuranaga, Y.; Inomata, Y.; Komura, K.; Uchiyama, K.; Akao, Y. Organ-specific microRNAs (MIR122, 137, and 206) contribute to tissue characteristics and carcinogenesis by regulating pyruvate kinase M1/2 (PKM) expression. Int. J. Mol. Sci. 2018, 19, 1276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Wilpe, S.V.; Koornstra, R.; Brok, M.D.; Groot, J.W.D.; Blank, C.; Vries, J.D.; Gerritsen, W.; Mehra, N. Lactate dehydrogenase: Actate dehydrogenase: A marker of diminished antitumor immunity. OncoImmunology 2020, 9, 1731942. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Animals 12 01718 g001 550
Figure 1. The SCD1 mRNA relative expression (A); TVA accumulation (B); % of total identified fatty acids) and c9, t11-CLA synthesis (C); % of total identified fatty acids in MAC-T cells under different incubation conditions. CAY group (SCD1-inhibited MAC-T cells model): incubated with 50 μM TVA for 2.5 h after CAY10566 treatment at concentrations of 10 nM for 12 h, TVA group: only incubated with 50 μM TVA for 2.5 h, and control group: without CAY, TVA. ** p < 0.01.
Figure 1. The SCD1 mRNA relative expression (A); TVA accumulation (B); % of total identified fatty acids) and c9, t11-CLA synthesis (C); % of total identified fatty acids in MAC-T cells under different incubation conditions. CAY group (SCD1-inhibited MAC-T cells model): incubated with 50 μM TVA for 2.5 h after CAY10566 treatment at concentrations of 10 nM for 12 h, TVA group: only incubated with 50 μM TVA for 2.5 h, and control group: without CAY, TVA. ** p < 0.01.
Animals 12 01718 g001
Animals 12 01718 g002 550
Figure 2. Statistics of mass spectrum data (A) and the number of differentially expressed proteins difference protein between groups of MAC-T cells under different incubation conditions (B). CAY group (SCD1-inhibited MAC-T cells model): incubated with 50 μM TVA for 2.5 h after CAY10566 treatment at concentrations of 10 nM for 12 h, TVA group: only incubated with 50 μM TVA for 2.5 h, and control group: without CAY, TVA. Differential proteins were screened from quantifiable proteins with p-value < 0.05, and a differential expression change more than 1.3/1 was defined as the threshold of change for significant up-regulation, while a change of less than 1/1.3 was the threshold of change for significant down-regulation.
Figure 2. Statistics of mass spectrum data (A) and the number of differentially expressed proteins difference protein between groups of MAC-T cells under different incubation conditions (B). CAY group (SCD1-inhibited MAC-T cells model): incubated with 50 μM TVA for 2.5 h after CAY10566 treatment at concentrations of 10 nM for 12 h, TVA group: only incubated with 50 μM TVA for 2.5 h, and control group: without CAY, TVA. Differential proteins were screened from quantifiable proteins with p-value < 0.05, and a differential expression change more than 1.3/1 was defined as the threshold of change for significant up-regulation, while a change of less than 1/1.3 was the threshold of change for significant down-regulation.
Animals 12 01718 g002
Animals 12 01718 g003 550
Figure 3. Volcanic map of differentially expressed proteins in MAC-T cells under different incubation conditions. The horizontal coordinate indicates the fold change (log2 transformation), while the vertical coordinate indicates the p-values (log10 transformation). Red dots indicate upregulation and blue dots indicate downregulation.
Figure 3. Volcanic map of differentially expressed proteins in MAC-T cells under different incubation conditions. The horizontal coordinate indicates the fold change (log2 transformation), while the vertical coordinate indicates the p-values (log10 transformation). Red dots indicate upregulation and blue dots indicate downregulation.
Animals 12 01718 g003
Animals 12 01718 g004 550
Figure 4. The subcellular structure localization of differentially expressed proteins (A) and GO secondary annotation (B) in MAC-T cells under different incubation conditions. CAY group (SCD1-inhibited MAC-T cells model): incubated with 50 μM TVA for 2.5 h after the treatment of CAY10566 at concentrations of 10 nM for 12 h, TVA group: only incubated with 50 μM TVA for 2.5 h, and control group: without CAY, TVA.
Figure 4. The subcellular structure localization of differentially expressed proteins (A) and GO secondary annotation (B) in MAC-T cells under different incubation conditions. CAY group (SCD1-inhibited MAC-T cells model): incubated with 50 μM TVA for 2.5 h after the treatment of CAY10566 at concentrations of 10 nM for 12 h, TVA group: only incubated with 50 μM TVA for 2.5 h, and control group: without CAY, TVA.
Animals 12 01718 g004
Animals 12 01718 g005 550
Figure 5. Cluster analysis heat map based on GO enrichment classifications: biological process (A); cellular component (B); molecular function (C). GO enrichment was performed for differential proteins in MAC-T cells under different incubation conditions, and cluster analysis was performed to determine the correlation between the functions of differential proteins. The horizontal direction indicates the enrichment test results of different groups, and the vertical direction shows a description of differential proteins’ enrichment-related functions. Different sets of differential proteins and color blocks correspond with the functional descriptions and indicate the degree of enrichment, where red indicates stronger enrichment (the deeper the red, the stronger the enrichment) and blue indicates weaker enrichment (the lighter the blue, the weaker the enrichment).
Figure 5. Cluster analysis heat map based on GO enrichment classifications: biological process (A); cellular component (B); molecular function (C). GO enrichment was performed for differential proteins in MAC-T cells under different incubation conditions, and cluster analysis was performed to determine the correlation between the functions of differential proteins. The horizontal direction indicates the enrichment test results of different groups, and the vertical direction shows a description of differential proteins’ enrichment-related functions. Different sets of differential proteins and color blocks correspond with the functional descriptions and indicate the degree of enrichment, where red indicates stronger enrichment (the deeper the red, the stronger the enrichment) and blue indicates weaker enrichment (the lighter the blue, the weaker the enrichment).
Animals 12 01718 g005
Animals 12 01718 g006 550
Figure 6. KEGG enrichment pathway map of the up-regulated proteins (A) and down-regulated proteins (B) in MAC-T cells under different incubation conditions. CAY group (SCD1-inhibited MAC-T cells model): incubated with 50 μM TVA for 2.5 h after the treatment of CAY10566 at concentrations of 10 nM for 12 h, TVA group: only incubated with 50 μM TVA for 2.5 h, and control group: without CAY, TVA.
Figure 6. KEGG enrichment pathway map of the up-regulated proteins (A) and down-regulated proteins (B) in MAC-T cells under different incubation conditions. CAY group (SCD1-inhibited MAC-T cells model): incubated with 50 μM TVA for 2.5 h after the treatment of CAY10566 at concentrations of 10 nM for 12 h, TVA group: only incubated with 50 μM TVA for 2.5 h, and control group: without CAY, TVA.
Animals 12 01718 g006
Animals 12 01718 g007 550
Figure 7. Tandem mass tag labeling and parallel reaction monitoring analysis among the differentially expressed proteins found in MAC-T cells under different incubation conditions. CAY group (SCD1-inhibited MAC-T cells model): incubated with 50 μM TVA for 2.5 h after the treatment of CAY10566 at concentrations of 10 nM for 12 h, TVA group: only incubated with 50 μM TVA for 2.5 h, and control group: without CAY, TVA.
Figure 7. Tandem mass tag labeling and parallel reaction monitoring analysis among the differentially expressed proteins found in MAC-T cells under different incubation conditions. CAY group (SCD1-inhibited MAC-T cells model): incubated with 50 μM TVA for 2.5 h after the treatment of CAY10566 at concentrations of 10 nM for 12 h, TVA group: only incubated with 50 μM TVA for 2.5 h, and control group: without CAY, TVA.
Animals 12 01718 g007
Animals 12 01718 g008 550
Figure 8. The mRNA relative expression levels of 6-phosphogluconolactonase (PGLS), phosphofructokinase, liver type (PFKL), aldolase, fructose-bisphosphate A (ALDOA), triosephosphate isomerase 1 (TPI1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), phosphoglycerate kinase 1 (PGK1), phosphoglycerate mutase 1 (PGAM1), enolase 1 (ENO1), pyruvate kinase M1/2 (PKM), lactate dehydrogenase A (LDHA), and lactate dehydrogenase B (LDHB) were successfully and individually reduced by their specific siRNA sequences in MAC-T cells. ** p < 0.01.
Figure 8. The mRNA relative expression levels of 6-phosphogluconolactonase (PGLS), phosphofructokinase, liver type (PFKL), aldolase, fructose-bisphosphate A (ALDOA), triosephosphate isomerase 1 (TPI1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), phosphoglycerate kinase 1 (PGK1), phosphoglycerate mutase 1 (PGAM1), enolase 1 (ENO1), pyruvate kinase M1/2 (PKM), lactate dehydrogenase A (LDHA), and lactate dehydrogenase B (LDHB) were successfully and individually reduced by their specific siRNA sequences in MAC-T cells. ** p < 0.01.
Animals 12 01718 g008
Animals 12 01718 g009 550
Figure 9. After individual PGLS, PFKL, ALDOA, TPI1, GAPDH, PGK1, PGAM1, ENO1, PKM, LDHA, and LDHB knockdown, the mRNA relative expression levels of SCD1 (stearoyl CoA desaturase 1) were altered in MAC-T cells. * p < 0.05, ** p < 0.01.
Figure 9. After individual PGLS, PFKL, ALDOA, TPI1, GAPDH, PGK1, PGAM1, ENO1, PKM, LDHA, and LDHB knockdown, the mRNA relative expression levels of SCD1 (stearoyl CoA desaturase 1) were altered in MAC-T cells. * p < 0.05, ** p < 0.01.
Animals 12 01718 g009
Animals 12 01718 g010 550
Figure 10. The SCD1-dependent gene network of c9, t11-CLA biosynthesis from the perspective of energy metabolism. The green color protein indicates that they were verified to have negative relationships with SCD1.
Figure 10. The SCD1-dependent gene network of c9, t11-CLA biosynthesis from the perspective of energy metabolism. The green color protein indicates that they were verified to have negative relationships with SCD1.
Animals 12 01718 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peng, L.; Bai, G.; Wang, C.; Dong, J.; Liu, Y.; Sun, Z.; Zhen, Y.; Qin, G.; Zhang, X.; Demelash, N.; et al. Proteomics Insights into the Gene Network of cis9, trans11-Conjugated Linoleic Acid Biosynthesis in Bovine Mammary Gland Epithelial Cells. Animals 2022, 12, 1718. https://0-doi-org.brum.beds.ac.uk/10.3390/ani12131718

AMA Style

Peng L, Bai G, Wang C, Dong J, Liu Y, Sun Z, Zhen Y, Qin G, Zhang X, Demelash N, et al. Proteomics Insights into the Gene Network of cis9, trans11-Conjugated Linoleic Acid Biosynthesis in Bovine Mammary Gland Epithelial Cells. Animals. 2022; 12(13):1718. https://0-doi-org.brum.beds.ac.uk/10.3390/ani12131718

Chicago/Turabian Style

Peng, Liying, Ge Bai, Chunzheng Wang, Jianan Dong, Yongjun Liu, Zhe Sun, Yuguo Zhen, Guixin Qin, Xuefeng Zhang, Natnael Demelash, and et al. 2022. "Proteomics Insights into the Gene Network of cis9, trans11-Conjugated Linoleic Acid Biosynthesis in Bovine Mammary Gland Epithelial Cells" Animals 12, no. 13: 1718. https://0-doi-org.brum.beds.ac.uk/10.3390/ani12131718

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop