- 1Department of Pathology, Xuzhou Medical University, Xuzhou, China
- 2Department of Human Anatomy, Xuzhou Medical University, Xuzhou, China
- 3Department of Orthopedics, The Affiliated Suqian Hospital of Xuzhou Medical University, Suqian, China
- 4Department of Oncological Radiotherapy, The Second Affiliated Hospital of Xuzhou Medical University, Xuzhou, China
Background: Tumors frequently evade immune surveillance through multiple pathways to escape T cell recognition and destruction. Previous studies indicated that lipid metabolism alteration could affect the anti-tumor immunity of cancer cells. Nonetheless, the studies that investigated lipid metabolism-related gene for cancer immunotherapy are still few.
Materials and methods: By mining the TCGA database, we screened out carnitine palmitoyltransferase-2 (CPT2), a key enzyme in the fatty acid β-oxidation (FAO) process associated with anti-tumor immunity. We then analyzed the gene expression and clinicopathological features of CPT2 using open-source platforms and databases. Molecular proteins interacting with CPT2 were also identified using web interaction tools. Subsequently, the relationship between CPT2 and survival was analyzed in cancer patients.
Results: Our study revealed that CPT2 played a vital role in tumor microenvironment and immune response signaling pathways. We have also demonstrated that increased CPT2 gene expression could enhance the level of tumor immune cell infiltration. Furthermore, high CPT2 expression positively related with overall survival associated with immunotherapy. CPT2 expression was also associated with the prognosis of human cancers, suggesting that CPT2 may be a potential biomarker for predicting the efficacy of cancer immunotherapy.
Conclusion: To the best of our knowledge, the relationship between CPT2 and tumor immune microenvironment was first proposed in this study. Therefore, further studies on CPT2 may provide new insights into the development of effective cancer immunotherapy.
Introduction
It is well established that immune escape and immune infiltration are associated with tumor growth and development (Alemohammad et al., 2022). Different types of cells such as T (Treg) cells, M2-type macrophages and immune checkpoints during the immune response are involved in the regulation of the immune regulatory system (Butt and Mills, 2014). Immune checkpoint therapy has led to major breakthroughs in the treatment of tumors previously considered to be highly lethal. Specifically, the presence of tumor-infiltrating lymphocytes (TILs), T cell (CD4+, CD8+, NK, γδ), macrophages, dendritic cells (DCs), B cells, mast cells, and neutrophils, as well as their localization, formation and density, can influence both the immune response and patient survival (Nossal, 1975; Fridman et al., 2012). Improving the survival of cancer patients through tumor immunotherapy has brought a new revolution in the field of oncology (Waldman et al., 2020). Tumor immunotherapy has become a popular and advanced treatment nowadays. However, it now appears that immunotherapy has reached a bottleneck. Only by better understanding how cancer cells evade immunotherapy can cancer be better treated (Cerezo and Rocchi, 2020; Waldman et al., 2020).
Cancer cells take up nutrients in tumor immune microenvironment (TIME) to create a favorable growing environment for tumor cells. Multiple immune cells and stromal cells make up the tumor tissue and they collaborate with each other to facilitate the reprogramming of the metabolic environment (Frades et al., 2021). A key determinant of the efficacy of immunotherapy appears to be metabolic interference between immune effector cells and tumor cells (Cerezo and Rocchi, 2020). The well-known Warburg effect is that cancer cells rely on glycolysis even in the presence of sufficient oxygen (Warburg, 1956; Li and Zhang, 2016). In addition, metabolic reprogramming of cancer cells also manifests as aberrant lipid metabolism and the dysregulation of amino acid metabolism (Li and Zhang, 2016). Although most of the knowledge about metabolic dysregulation in cancer has a focus on carbohydrates, we decided to focus our interest on the lipid metabolism of cancer cells (Merino Salvador et al., 2017; Cerezo and Rocchi, 2020). And the increase in de novo lipogenesis is considered to be a new hallmark of many aggressive forms of cancer (Robey et al., 2015). The metabolism of lipids is greatly altered in the proliferating cells. Unlike normal cells, the biosynthesis of fatty acids (FAs) in tumor tissue is significantly increased to meet the demand for local nutrients, oxygen, lipid synthesis of membrane and signaling molecules (Robey et al., 2015). Lipid metabolism dysfunction in the tumor microenvironment (TME) has been shown to affect tumor-associated immune responses (Yu et al., 2021). And a number of enzymes involved in lipid metabolism have been identified to be associated with carcinogenesis (Menendez and Lupu, 2007; Zhou et al., 2013; Mashimo et al., 2014; Bian et al., 2021).
Carnitine palmitoyltransferase (CPT) is an acyltransferase localised to the inner mitochondrial membrane in the process of fatty acid oxidation (FAO). Carnitine palmitoyltransferase (CPT) deficiency is a common disease that occurs during mitochondrial fatty acid oxidation. (Liu et al., 2022). The CPT system consists of two proteins, CPT1 and CPT2, which have different functions and localizations (Bonnefont et al., 2004). CPT1 locates the external mitochondrial membrane and acts as a catalyst for the formation of acylcarnitins based on carnitine and acyl CoA. The main function of inner mitochondrial membrane CPT2 is to form acyl coenzyme A from acylcarnitine and CoA (Bonnefont et al., 2004). In recent years, the role of the CPT system in cancer metabolic mechanisms has been updated due to the increasing understanding of fatty acid oxidation in cancer (Li et al., 2021). Studies have shown that PD-1 supports lipolysis by upregulating CPT1A and FAO (Patsoukis et al., 2015). Moreover, targeted CPT1 therapy reduces the formation of invasive pseudopodia in hepatocellular cancer cells and eliminates the tumor-promoting effect of tumor-associated macrophages (TAMs) (Su et al., 2020). CPT2 deficiency or dysfunction usually results in some lipid metabolic disorders such as diabetes, cardiac hypertrophy, NAFLD and obesity (Gu et al., 2017; Pereyra et al., 2020). CPT2 is highly expressed in recurrent breast cancer and associated with poor prognosis of patients. CPT2 increase promotes the growth and radiation resistance of breast cancer cells (Han et al., 2019). However, CPT2 was shown to be downregulated in HCC tissues and associated with vascular invasion and tumor tissue differentiation (Fujiwara et al., 2018; Lin et al., 2018). Silencing of CPT2 contributes to the metastatic and tumorigenic activity of hepatocellular carcinoma cells (Lin et al., 2018). Similarly, CPT2 knockdown has been demonstrated to promote the proliferation, migration and invasion of CRC cells in vitro, and to induce cell cycle arrest and apoptosis (Li et al., 2021). In addition, upregulation of CPT2 has been shown to attenuate CRC tumor development and increase chemosensitivity via a hyperactive Wnt/β-catenin pathway in vivo (Liu et al., 2022). Unfortunately, the role of CPT2 in tumor immune microenvironment remains unknown.
In this study, we explored the expression of CPT2 in human cancer and identified its relationship with the level of immune cell infiltration. And the entire process of exploring the function of CPT2 in the pan-cancer landscape was showed in Figure 1.
Methods
Enrichment analysis of CPT2-related genes
TCGA database were used to screen lipid metabolism enzymes related to effector T cells. The screened genes were then plotted via a Venn diagram using the JVENN (http://jvenn.toulouse.inra.fr/app/example.html) website. Through the utilization of the purity-adjusted Spearman’s rank correlation test, p-values and partial cor values were obtained. STRING is a database dedicated to the study of protein interactions throughout the organism (Huang et al., 2018). We submitted CPT2 on the STRING (https://string-db.org) website. STRING generated a network of key proteins that have interactions with CPT2. BioGRID (Biological General Repository for Interaction Datasets, thebiogrid.org) is a comprehensive, open access biomedical database resource that collates protein, genetic and chemical interactions from multiple species including yeast, worms, flies, mice and humans. We searched for CPT2 on BioGRID (https://thebiogrid.org) website and then obtained a genes network associated with CPT2 from the same organism and different organism.
Linkedomics database
The open multi-omics database LinkedOmics (http://www.linkedomics.org/login.php) contains data from 10 cancer cohorts from the Clinical Proteomics Tumor Analysis Consortium (CPTAC) and 32 TCGA cancer types. The gene set enrichment analysis (GSEA) and GO analysis (biological process) of CPT2 were generated from LinkedOmics.
Differential expression analysis
The Human Protein Atlas (HPA) (https://www.proteinatlas.org) provides information on the tissue and cellular distribution of 26,000 human proteins. HPA covers the expression of each protein in 64 cell lines, 48 human normal tissues and 20 tumor tissues. We used the HPA database to evaluate immunohistochemical staining of CPT2 in liver cancer, stomach cancer, lung cancer, renal cancer, pancreatic cancer, colorectal cancer and normal tissue.
The Timer website’s DiffExp module displays the difference gene expression patterns in all TCGA tumors between the tumor and normal tissue. In a box plot, the distribution of gene expression levels is displayed, and the Wilcoxon test determines the statistical significance of differential expression. We analyzed whether CPT2 was up- or downregulated in tumors for each cancer type compared to normal tissues.
TNMplot: Differential gene expression analysis
TNM plot (https://tnmplot.com/analysis/) is a comprehensive database using existing transcriptome-level datasets, which can be mined by comparing normal, tumor and metastasis data for all genes in real time. The entire database contains 56,938 samples, including 33,520 samples from 3,180 gene chip-based studies (453 metastatic, 29,376 tumorigenic and 3,691 normal samples), 11,010 samples from TCGA (394 metastatic, 9,886 tumorigenic and 730 normal samples), 1,193 samples from TARGET (1 metastatic, 1,180 tumorigenic and 12 normal samples) and 11,215 normal samples from GTEx. The Normal and Tumor analysis page provides detailed analysis for CPT2 in a selected tissue type using gene chip based data.
Immune infiltration analysis
A thorough investigation of the relationship between cancerous cells and the host immune system is necessary in light of the recent clinical success of cancer immunotherapy in a range of tumor types. TIMER (https://cistrome.shinyapps.io/timer/) provides 6 major analytic modules that allow users to interactively explore the associations between immune infiltrates and a wide spectrum of factors, including gene expression, clinical outcomes, somatic mutations, somatic copy number alterations, differential expression and correlation analysis. TIMER applies the deconvolution algorithm to estimate the number of six different immune infiltrating cell types (B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells). TISIDB (http://cis.hku.hk/TISIDB) is a web portal for studying function of specific genes in tumor-immune interactions through literature mining and high-throughput data analysis. Correlation of CPT2 and immune effector cells in Low grade glioma (LGG) was analyzed by TISIDB website.
Survival analysis
The main objective of the Kaplan Meier plotter is to assess the correlation between expression of genes (more than 30,000 samples from 21 tumor types) and survival rates. Sources of the database include GEO, EGA and TCGA. One of the finest methods for calculating the percentage of participants who continue to live following therapy is the Kaplan-Meier estimate. The number of participants who survive or are saved following an intervention is measured over time in a clinical study or community trial to determine the intervention’s effectiveness. The study of cohort data is known as survival analysis, and survival time is the period of time between a stated point and the occurrence of an event, such as death (Goel et al., 2010). Survival curves were calculated using the Kaplan-Meier method, and differences between groups were compared using the log-rank test (Zou et al., 2021). p-values less than 0.05 were considered as a significant difference. We have not only investigated the prognosis of CPT2 in patients with liver, lung, kidney, thyroid, pancreatic and gastric cancers, but also analyzed the clinical association between immune characteristics and ultimate patient survival by Kaplan Meier.
Genetic alteration analysis
The cBioPortal for Cancer Genomics is an online tool for exploring, visualizing, and interpreting multimodal cancer genomes data (http://cbioportal.org). The portal transforms genetic, epigenetic, gene expression, and proteomic processes into plainly understandable terms using molecular profiling data from cancer tissues and cell types. We obtained the frequency, mutation type, copy number alteration (CNA) and structural variants of TCGA tumors in the “Cancer Type Summary” by a quick search of CPT2 on the cBioPortal (https://www.cbioportal.org/) website. We obtained the mutation site information and 3D structure of CPT2 protein by “mutation” model.
Proportion of immune cells analysis
Cancer Single-cell Expression Map (https://ngdc.cncb.ac.cn/cancerscem/index) is a freely accessible database for the gathering, examination, and presentation of single-cell RNA-Seq data from human cancer. With a thorough online analysis platform built into the database, a multi-level study sheds light on the tumor microenvironment of many types of human cancers.
Gene-metabolite interaction analysis
Metabolic Atlas is a web platform for exploring and analyzing open-source genome scale metabolic models (GEMs). The purpose is to collect curated GEMs and align these models with FAIR principles. The website features GEM visualizations and comparisons, as well as access to resources, algorithms, other databases, and more general software applications. Metabolic Atlas is designed for use in metabolomics, clinical chemistry, biomarker development, and general education.
Results
Identification of CPT2 as the hub gene
To explore the relationship between lipid metabolism and tumor immunity, we screen out remarkable lipid metabolism-related genes associated with effector T cells through TCGA database. In total, 1,468 differentially expressed genes (DEGs) related to effector T cells were obtained. Among them, there are five lipid metabolic enzymes (LPCAT2, PLCG2, PLCB2, CPT2 and PLA2G4A) that are significantly associated with effector T cells (Figure 2A). The two types of searches on the string website, single protein by name/identifier and geneset by pathway/process/disease/publication, allow us to obtain two different maps of the CPT2 protein interaction network (Figure 2B, C). Through BioGRID websites, we were able to get the PPI network to analyze the CPT2 protein interactions (Figure 2D). Moreover, to understand the biological function of CPT2 in LGG, GO functional enrichment analysis revealed that CPT2 was largely associated with immune regulation including neutrophil mediated immunity, adaptive immune response, lymphocyte mediated immunity, T cell activation, interleukin-1 production, regulation of immune effector process, interleukin-6 production and immune response regulating signaling pathway (Figure 2E). More specifically, gene set enrichment analysis (GSEA) demonstrated a positive correlation between CPT2 expression levels and immune regulation including neutrophil mediated immunity, adaptive immune response, lymphocyte mediated immunity, T cell activation, interleukin-1 production, regulation of immune effector process, interleukin-6 production and immune response regulating signaling pathway (Figure 2F).
FIGURE 2. Identification of CPT2 as the hub gene. (A) Venn diagram displaying the overlaps of top 75 enzymes involved in lipid metabolism and Genes with a significantly positive correlation with Effector T cell DEGs. (B, C) Two alternative maps of the CPT2 protein interaction network using the two search methods available on the string website: single protein by name/identifier and geneset by pathway/process/disease/publication. (D) Examining significant CPT2-related protein interactions in the BioGRID database. (E, F) GO biological process analysis and GESA analysis of CPT2 in LGG.
Analysis of CPT2 and clinical features in human cancer
To assess the clinical features of CPT2 expression, we evaluated the immunohistochemical staining of CPT2 via HPA database. Results indicated that CPT2 was strongly marked in normal tissues than tumor tissues (Figure 3A). In liver hepatocellular Carcinoma (LIHC), lung adenocarcinoma (LUAD), kidney renal clear cell carcinoma (KIRC), thyroid carcinoma (THCA), pancreatic adenocarcinoma (PAAD), and stomach adenocarcinoma (STAD), low levels of CPT2 expression were found to be related with a poor overall survival (OS), according to an investigation of the Kaplan-Meier Plottr database. (Figure 3B). Subsequently, we next used the TIMER website to examine the differential expression of CPT2 in various tumor tissues and normal tissues. The findings demonstrated that CPT2 expression was significantly lower in cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), head-neck squamous cell carcinoma (HNSC), kidney chromophobe (KICH), kidney renal papillary cell carcinoma (KIRP), LIHC, lung squamous cell carcinoma (LUSC) and rectal adenocarcinoma (READ) (Figure 3C). Furthermore, it is identified that CPT2 was overexpressed in normal tissues compared with tumor tissues by TNMplot analysis (Figure 3D).
FIGURE 3. Analysis of CPT2 and clinical features in human cancer. (A) CPT2 protein expression levels in normal and tumor tissues were compared through immunohistochemical images involving liver cancer, stomach cancer, lung cancer, kidney cancer, pancreatic cancer, colorectal cancer, and low grade glioma. (B) Survival curves for patients with high (red) and low (black) expression of CPT2 were displayed using the Kaplan-Meier Plotter database for liver hepatocallular carcinoma, lung adenocarcinoma, kidney renal clear cell thyroid carcinoma, pancreatic ductal adenocarcinoma, and stomach adenocarcinoma. The threshold of Cox p-value <0.05. (C) The expression level of CPT2 in different types of tumor tissues and normal tissues in TIMER database. Distributions of gene expression levels are displayed using box plots. Red boxes represent tumor tissues; blue boxes represent normal tissues. The Wilcoxon test results’ statistical significance is indicated by the number of stars (∗:
CPT2 is associated with immune cells in the tumor microenvironment of LGG
It is well established that cells of the immune system in the tumor microenvironment may have an impact on the prognosis and survival of cancer patients. Based on the TIMER 2.0 database and taking into account the strong connection between CPT2 and immune infiltration, we discovered a substantial positive correlation between CPT2 and the degree of immune cell infiltration in LGG. (Figure 4A). Furthermore, we discovered a positive correlation between CPT2 expression and the infiltration levels of effector T cells, exhausted T cells, resident memory T cells, Th1-like cells, and naive T cells (Figure 4B). The correlations of immune cell infiltrations and CPT2 in LGG were systematically investigated through TISIDB database, which revealed significant positive correlations between CPT2 expression and Tcm CD8 T-cell subsets, Imm B-cell, Act DC-cell, iDC-cell, Macrophage cell, Mast cell, Tcm CD4 T-cell, Mem B-cell, NK cell, NK T-cell, Tfh cell and Tgd cell (Figure 4C). These results indicated that CPT2 was positively correlated with immune cells in tumor microenvironment.
FIGURE 4. Correlations of CPT2 expression with immune cells. (A) Correlation between CPT2 expression and B cell, CD8+ T cell, CD4+ T cell, Macrophage, Neutrophil and Dentritic cell in LGG. Scatter plots display purity-corrected partial Spearman’s rho values as well as statistical significance. The gene expression levels against tumor purity are always displayed on the left-most panel. (B) Correlation between CPT2 expression and highly connected immune cell, including Effector T cell, Exhausted T cell, Resident memory T cell Th1-like cell and Naïve T cell in LGG. (C) Bioinformatics analysis of the correlation between CPT2 and immune effector cells using TISIDB database.
Correlation between CPT2 expression and overall survival associated with immunotherapy
Figure 5A illustrates that a reduction in CD8+ T cells and type 1 T-helper cells decreases survival in lung adenocarcinoma patients with high CPT2 expression. Moreover, Figure 5B indicates that overall survival was significantly higher when CD8+ T cells, CD4+ T cells and eosinophils were enriched in squamous cell lung cancer patients with high CPT2 expression. In order to learn more about the connection between CPT2 and immunotherapy, we discovered that patients with lung squamous cell carcinoma who received anti-PDL1 and anti-CTLA-4 treatment had higher survival rates when CPT2 was highly expressed (Figure 5C). These results confirmed the positive relationship between CPT2 expression and overall survival associated with immunotherapy.
FIGURE 5. Correlation between CPT2 expression and overall survival associated with immunotherapy. (A, B) Kaplan-Meier curves show the impact of CPT2 and immune cells (CD4+ T cells, CD8+ T cells, type 1 T-helper cells, type 2 T-helper cells) expression levels on patient survival prognosis. (C) Survival analysis of all tumor types in Kaplan-Meier include Bladder, Esophageal adeno, Glioblastoma, Hepatocellular carcinoma, HNSCC, Melanoma, NSCLC, NSLC, Urothelial treated with immune checkpoint (PD-L1 and CTLA-4) or without any immunotherapy.
CPT2 positively regulates immune cells in LGG
TIMER 2.0 web server was applied to assess the effects of CPT2 on the level of immune cell infiltration in pan-cancer. It was demonstrated that CPT2 promoted the infiltration of CD8+ T cells and monocytes in LGG (Figure 6A); CPT2 promoted the infiltration of neutrophil in LGG (Figure 6B); CPT2 promoted the infiltration of CD4+ T cell in LGG (Figure 6C); CPT2 promoted the infiltration of Macrophage in LGG (Figure 6D). These data demonstrated that CPT2 might positively regulate immune cells in LGG.
FIGURE 6. Heat maps and scatter plots of CPT2 correlated with immune cell infiltration level. The left panels are heat maps showing the correlation between all tumor types and immune cells in TCGA, with red representing positive correlation and blue representing negative correlation, the darker the color the higher the correlation. The right panels are scatter plot. (A) The relationship between CPT2 expression and CD8+ T cells, and monocytes infiltration levels was analyzed by the XCELL algorithm in LGG. (B) The relationship between CPT2 expression and neutrophil infiltration level was analyzed by the TIMER algorithm in LGG. (C) The relationship between CPT2 expression and CD4+ T cells infiltration level was analyzed by the CIBERSORT-ABS and XCELL algorithms in LGG. (D) The relationship between CPT2 expression and Macrophage infiltration level was analyzed by the EPIC, XCELL, and QUANTISEQ algorithms in LGG.
Mutation analysis of CPT2 in human cancer
Figure 7A illustrated the X113_splice altered site of 3D structure of CPT2 protein. An alteration at position 113 of the CPT2 protein has been found in cutaneous melanoma (X113_splice). According to the mutation sites and types, we found missense mutation of CPT2 as the main mutation type (Figure 7B). The mutation type, frequency, CAN, and structural variant of CPT2 in all TCGA cancers were obtained from the cBioPortal website. According to Figure 7C, ovarian epithelial tumors had the greatest modification frequency of any TCGA tumors (3.94%), with amplification accounting for the majority of the frequency (3.25%). Analysis of CPT2 in pan-cancer illustrated that CPT2 is amplified in a subset. The cBioPortal results demonstrate that the samples with amplified CPT2 have higher levels of CPT2 mRNA. (Figure 7D).
FIGURE 7. Mutation analysis of CPT2 in human cancer. (A) X113_splice mutation site of CPT2 in the 3D structure at cBioPortal website. (B) A total of 74 mutation sites were detected between amino acids 0 and 658 of CPT2 protein at cBioPortal website. (C) Alteration frequency of CPT2 in different tumor types. (D) A plot showing the relationship between CPT2 mRNA abundance and CNA in the CPT2 gene in tumors from the selected pan-cancer atlas study.
Discussion
Regulatory B cells, CD4+ Th2 cells, bone marrow myeloid-derived suppressor cells (BM-MDSCs), CD4+ T regulatory cells suppressor cells and alternatively activated macrophages (M2) are engulfed by the tumor development (Edechi et al., 2019). The findings suggested that CPT-1A expression and methylation is negatively associated with these cells, suggesting that survival and prognosis may be affected by inhibiting tumor-infiltrating lymphocytes (TILs) response (Das et al., 2022). In an earlier study, CPT-1A was shown to be a novel biomarker for the diagnosis and prediction of breast cancer (Rhodes et al., 2004). CPT-1A genomic alterations have been detected, in 9% of breast cancers. Alterations in the CPT1A genome have been confirmed to be associated with poor prognosis. Moreover, CPT-1A reduces the expression of CPT-1A through methylation and may be associated with a good prognosis. To date, numerous studies have demonstrated that CPT1 is closely associated with the occurrence and development of tumor cells in lung cancer (Zaugg et al., 2011), colorectal cancer (Wang et al., 2018), breast cancer (Linher-Melville et al., 2011), hepatocellular carcinoma (Brown et al., 2018), gastric cancer (Wang et al., 2020), ovarian cancer (Shao et al., 2016) and prostate cancer (Joshi et al., 2019). Nevertheless, the impact of CPT2 on the progress of cancers is still not well understood.
To better understand the role of CPT2 in cancer development and prognosis, we performed comparative data mining on numerous gene expression data sets. In the current study, we investigated the impact of CPT2 on pan cancer and the resulting clinical implications. We first performed differential and survival analysis based on TCGA and GEO databases. We found that CPT2 was highly expressed in most normal tissues compared to cancers, and its low expression was associated with poor prognosis in LGG. A total of 40 CPT2 mutations have been identified, distributed throughout the coding sequence of the gene. In addition to the prevalent S113L exon 3, exons 4 and 5 were hotspots for mutations. Missense mutations were the predominant mutation type, accounting for about 77%, while the remaining 23% were predicted to be truncated proteins (Bonnefont et al., 2004). Through the mutation analysis website, we found that CPT2 had a high mutation rate in ovarian epithelial tumors and sarcoma, and the mutation types were predominantly missense mutations.
Fatty acids are a large, critical, energy source that can meet the substantial fuel needs of acutely propagating cancer cells (De Rosa et al., 2012; Hoang-Minh et al., 2018; Duman et al., 2019). Recently, two groups have further revealed that FA metabolism accelerates the development of brain metastases from breast cancer (Gupta et al., 2020), while glioblastoma multiforme (GBM) metabolism is considerably less dependent on glycolysis (Reinfeld et al., 2021). Moreover, it has been reported that FAs provide essential nutrients for the growth of tumor cells, and that a deficiency of FAs can lead directly to the death of FA-dependent tumor cells (Aloia et al., 2019). Previously, many studies have demonstrated regulation of FA synthesis in immune cells is also essential for immunity against cancer and the synthesis of FA is regulated by the metabolism of lipids and lipoproteins through the PPAR activators (Nakamura et al., 2014; Qian et al., 2018). However, it is unclear whether enhanced FA metabolism in recurrent and radioresistant GBM, can generate or promote a state of immunosuppression leading to immune evasion (Jiang et al., 2022). FA deficiency also blocks the growth of MDSC and Tregs that use FAOs. Promoting the FOA process may be an effective strategy to improve antitumor immunity, as FAO has been suggested to favour the longevity of memory T cells (Chowdhury et al., 2018). Furthermore, several lipid metabolism-related genes were identified to be associated with tumor immunology. Inhibition of ZDHHC9, a palmitoyl transferase, could interrupt the palmitoylation of PD-L1, sensitized tumor cells to immune T cell attacking (Yang et al., 2019; Liao et al., 2021). The palmitoyltransferase ZDHHC3 (DHHC3) has also been shown to be involved in acetyl transfer during PD-L1 palmitoylation (Yao et al., 2019; Liao et al., 2021). PD-L1 palmitoylation can inhibit PD-L1 expression and function to promote cancer cell killing by T cells (Liao et al., 2021). Accumulating studies strengthened the role of lipid metabolism enzymes for TIME modulation in tumor expansion (Hodson and Fielding, 2013; Bai et al., 2015; Cotte et al., 2018).
Hence, our study revealed that CPT2, a key enzyme in the FAO process, was positively associated with the level of immune cell infiltration, which suggested that CPT2 could affect tumor progression by reprogramming lipid metabolism in the TIME. Based on previous reports of PD-1 inducing upregulation of CPT1A and FAO (Patsoukis et al., 2015), we further explored the relationship between CPT2 and immunotherapy. It was demonstrated that survival in squamous cell lung cancer was significantly increased when CPT2 was highly expressed and treated with anti-PDL1 or anti-CTLA-4. Therefore, further studies would be warranted to verify the value of the prognostic and predictive benefit of CPT2 agonist.
In addition, there is increasing evidence that tumor metabolism can affect immunological molecules by releasing metabolites as well as to being essential for supporting growth and survival. CPT2, a key molecule in lipid metabolism, and the crosstalk between tumor metabolites including Malonyl-CoA, 1,1-Dimethylbiguanide, L-Palmitoylcarnitine, Estradiol, Propionyl-CoA, Palmityl-CoA, Palmitaldehyde, Coenzyme A, and Palmitic acid may offer new hope for tumor immunotherapy.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
DP and LZ developed the study concept and were responsible for the study design. XJ, WD, and CS wrote the manuscript and analyzed data. MK and QS revised the manuscript. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2023.1186064/full#supplementary-material
References
Alemohammad, H., Najafzadeh, B., Asadzadeh, Z., Baghbanzadeh, A., Ghorbaninezhad, F., Najafzadeh, A., et al. (2022). The importance of immune checkpoints in immune monitoring: A future paradigm shift in the treatment of cancer. Biomed. Pharmacother. 146, 112516. doi:10.1016/j.biopha.2021.112516
Aloia, A., Mullhaupt, D., Chabbert, C. D., Eberhart, T., Fluckiger-Mangual, S., Vukolic, A., et al. (2019). A fatty acid oxidation-dependent metabolic shift regulates the adaptation of BRAF-mutated melanoma to MAPK inhibitors. Clin. Cancer Res. 25, 6852–6867. doi:10.1158/1078-0432.CCR-19-0253
Bai, Y., Mccoy, J. G., Levin, E. J., Sobrado, P., Rajashankar, K. R., Fox, B. G., et al. (2015). X-ray structure of a mammalian stearoyl-CoA desaturase. Nature 524, 252–256. doi:10.1038/nature14549
Bian, X., Liu, R., Meng, Y., Xing, D., Xu, D., and Lu, Z. (2021). Lipid metabolism and cancer. J. Exp. Med. 218, e20201606. doi:10.1084/jem.20201606
Bonnefont, J. P., Djouadi, F., Prip-Buus, C., Gobin, S., Munnich, A., and Bastin, J. (2004). Carnitine palmitoyltransferases 1 and 2: Biochemical, molecular and medical aspects. Mol. Asp. Med. 25, 495–520. doi:10.1016/j.mam.2004.06.004
Brown, Z. J., Fu, Q., Ma, C., Kruhlak, M., Zhang, H., Luo, J., et al. (2018). Carnitine palmitoyltransferase gene upregulation by linoleic acid induces CD4(+) T cell apoptosis promoting HCC development. Cell Death Dis. 9, 620. doi:10.1038/s41419-018-0687-6
Butt, A. Q., and Mills, K. H. (2014). Immunosuppressive networks and checkpoints controlling antitumor immunity and their blockade in the development of cancer immunotherapeutics and vaccines. Oncogene 33, 4623–4631. doi:10.1038/onc.2013.432
Cerezo, M., and Rocchi, S. (2020). Cancer cell metabolic reprogramming: A keystone for the response to immunotherapy. Cell Death Dis. 11, 964. doi:10.1038/s41419-020-03175-5
Chowdhury, P. S., Chamoto, K., Kumar, A., and Honjo, T. (2018). PPAR-induced fatty acid oxidation in T cells increases the number of tumor-reactive CD8(+) T cells and facilitates anti-PD-1 therapy. Cancer Immunol. Res. 6, 1375–1387. doi:10.1158/2326-6066.CIR-18-0095
Cotte, A. K., Aires, V., Fredon, M., Limagne, E., Derangere, V., Thibaudin, M., et al. (2018). Lysophosphatidylcholine acyltransferase 2-mediated lipid droplet production supports colorectal cancer chemoresistance. Nat. Commun. 9, 322. doi:10.1038/s41467-017-02732-5
Das, M., Giannoudis, A., and Sharma, V. (2022). The role of CPT1A as a biomarker of breast cancer progression: A bioinformatic approach. Sci. Rep. 12, 16441. doi:10.1038/s41598-022-20585-x
De Rosa, A., Pellegatta, S., Rossi, M., Tunici, P., Magnoni, L., Speranza, M. C., et al. (2012). A radial glia gene marker, fatty acid binding protein 7 (FABP7), is involved in proliferation and invasion of glioblastoma cells. PLoS One 7, e52113. doi:10.1371/journal.pone.0052113
Duman, C., Yaqubi, K., Hoffmann, A., Acikgoz, A. A., Korshunov, A., Bendszus, M., et al. (2019). Acyl-CoA-Binding protein drives glioblastoma tumorigenesis by sustaining fatty acid oxidation. Cell Metab. 30, 274–289. doi:10.1016/j.cmet.2019.04.004
Edechi, C. A., Ikeogu, N., Uzonna, J. E., and Myal, Y. (2019). Regulation of immunity in breast cancer. Cancers (Basel) 11, 1080. doi:10.3390/cancers11081080
Frades, I., Foguet, C., Cascante, M., and Arauzo-Bravo, M. J. (2021). Genome scale modeling to study the metabolic competition between cells in the tumor microenvironment. Cancers (Basel) 13, 4609. doi:10.3390/cancers13184609
Fridman, W. H., Pages, F., Sautes-Fridman, C., and Galon, J. (2012). The immune contexture in human tumours: Impact on clinical outcome. Nat. Rev. Cancer 12, 298–306. doi:10.1038/nrc3245
Fujiwara, N., Nakagawa, H., Enooku, K., Kudo, Y., Hayata, Y., Nakatsuka, T., et al. (2018). CPT2 downregulation adapts HCC to lipid-rich environment and promotes carcinogenesis via acylcarnitine accumulation in obesity. Gut 67, 1493–1504. doi:10.1136/gutjnl-2017-315193
Goel, M. K., Khanna, P., and Kishore, J. (2010). Understanding survival analysis: Kaplan-Meier estimate. Int. J. Ayurveda Res. 1, 274–278. doi:10.4103/0974-7788.76794
Gu, J. J., Yao, M., Yang, J., Cai, Y., Zheng, W. J., Wang, L., et al. (2017). Mitochondrial carnitine palmitoyl transferase-II inactivity aggravates lipid accumulation in rat hepatocarcinogenesis. World J. Gastroenterol. 23, 256–264. doi:10.3748/wjg.v23.i2.256
Gupta, K., Vuckovic, I., Zhang, S., Xiong, Y., Carlson, B. L., Jacobs, J., et al. (2020). Radiation induced metabolic alterations associate with tumor aggressiveness and poor outcome in glioblastoma. Front. Oncol. 10, 535. doi:10.3389/fonc.2020.00535
Han, S., Wei, R., Zhang, X., Jiang, N., Fan, M., Huang, J. H., et al. (2019). CPT1A/2-Mediated FAO enhancement-A metabolic target in radioresistant breast cancer. Front. Oncol. 9, 1201. doi:10.3389/fonc.2019.01201
Hoang-Minh, L. B., Siebzehnrubl, F. A., Yang, C., Suzuki-Hatano, S., Dajac, K., Loche, T., et al. (2018). Infiltrative and drug-resistant slow-cycling cells support metabolic heterogeneity in glioblastoma. EMBO J. 37, e98772. doi:10.15252/embj.201798772
Hodson, L., and Fielding, B. A. (2013). Stearoyl-CoA desaturase: Rogue or innocent bystander? Prog. Lipid Res. 52, 15–42. doi:10.1016/j.plipres.2012.08.002
Huang, J. K., Carlin, D. E., Yu, M. K., Zhang, W., Kreisberg, J. F., Tamayo, P., et al. (2018). Systematic evaluation of molecular networks for discovery of disease genes. Cell Syst. 6, 484–495. doi:10.1016/j.cels.2018.03.001
Jiang, N., Xie, B., Xiao, W., Fan, M., Xu, S., Duan, Y., et al. (2022). Fatty acid oxidation fuels glioblastoma radioresistance with CD47-mediated immune evasion. Nat. Commun. 13, 1511. doi:10.1038/s41467-022-29137-3
Joshi, M., Stoykova, G. E., Salzmann-Sullivan, M., Dzieciatkowska, M., Liebman, L. N., Deep, G., et al. (2019). CPT1A supports castration-resistant prostate cancer in androgen-deprived conditions. Cells 8, 1115. doi:10.3390/cells8101115
Li, H., Chen, J., Liu, J., Lai, Y., Huang, S., Zheng, L., et al. (2021). CPT2 downregulation triggers stemness and oxaliplatin resistance in colorectal cancer via activating the ROS/Wnt/β-catenin-induced glycolytic metabolism. Exp. Cell Res. 409, 112892. doi:10.1016/j.yexcr.2021.112892
Li, Z., and Zhang, H. (2016). Reprogramming of glucose, fatty acid and amino acid metabolism for cancer progression. Cell Mol. Life Sci. 73, 377–392. doi:10.1007/s00018-015-2070-4
Liao, Q., Zhou, Y., Xia, L., and Cao, D. (2021). Lipid metabolism and immune checkpoints. Adv. Exp. Med. Biol. 1316, 191–211. doi:10.1007/978-981-33-6785-2_12
Lin, M., Lv, D., Zheng, Y., Wu, M., Xu, C., Zhang, Q., et al. (2018). Downregulation of CPT2 promotes tumorigenesis and chemoresistance to cisplatin in hepatocellular carcinoma. Onco Targets Ther. 11, 3101–3110. doi:10.2147/OTT.S163266
Linher-Melville, K., Zantinge, S., Sanli, T., Gerstein, H., Tsakiridis, T., and Singh, G. (2011). Establishing a relationship between prolactin and altered fatty acid beta-oxidation via carnitine palmitoyl transferase 1 in breast cancer cells. BMC Cancer 11, 56. doi:10.1186/1471-2407-11-56
Liu, F., Li, X., Yan, H., Wu, J., Yang, Y., He, J., et al. (2022). Downregulation of CPT2 promotes proliferation and inhibits apoptosis through p53 pathway in colorectal cancer. Cell Signal 92, 110267. doi:10.1016/j.cellsig.2022.110267
Mashimo, T., Pichumani, K., Vemireddy, V., Hatanpaa, K. J., Singh, D. K., Sirasanagandla, S., et al. (2014). Acetate is a bioenergetic substrate for human glioblastoma and brain metastases. Cell 159, 1603–1614. doi:10.1016/j.cell.2014.11.025
Menendez, J. A., and Lupu, R. (2007). Fatty acid synthase and the lipogenic phenotype in cancer pathogenesis. Nat. Rev. Cancer 7, 763–777. doi:10.1038/nrc2222
Merino Salvador, M., Gomez De Cedron, M., Moreno Rubio, J., Falagan Martinez, S., Sanchez Martinez, R., Casado, E., et al. (2017). Lipid metabolism and lung cancer. Crit. Rev. Oncol. Hematol. 112, 31–40. doi:10.1016/j.critrevonc.2017.02.001
Nakamura, M. T., Yudell, B. E., and Loor, J. J. (2014). Regulation of energy metabolism by long-chain fatty acids. Prog. Lipid Res. 53, 124–144. doi:10.1016/j.plipres.2013.12.001
Nossal, G. J. (1975). Kinetics of antibody formation and regulatory aspects of immunity. Acta Endocrinol. Suppl. (Copenh) 194, 96–116. doi:10.1530/acta.0.080s096
Patsoukis, N., Bardhan, K., Chatterjee, P., Sari, D., Liu, B., Bell, L. N., et al. (2015). PD-1 alters T-cell metabolic reprogramming by inhibiting glycolysis and promoting lipolysis and fatty acid oxidation. Nat. Commun. 6, 6692. doi:10.1038/ncomms7692
Pereyra, A. S., Rajan, A., Ferreira, C. R., and Ellis, J. M. (2020). Loss of muscle carnitine palmitoyltransferase 2 prevents diet-induced obesity and insulin resistance despite long-chain acylcarnitine accumulation. Cell Rep. 33, 108374. doi:10.1016/j.celrep.2020.108374
Qian, X., Yang, Z., Mao, E., and Chen, E. (2018). Regulation of fatty acid synthesis in immune cells. Scand. J. Immunol. 88, e12713. doi:10.1111/sji.12713
Reinfeld, B. I., Madden, M. Z., Wolf, M. M., Chytil, A., Bader, J. E., Patterson, A. R., et al. (2021). Cell-programmed nutrient partitioning in the tumour microenvironment. Nature 593, 282–288. doi:10.1038/s41586-021-03442-1
Rhodes, D. R., Yu, J., Shanker, K., Deshpande, N., Varambally, R., Ghosh, D., et al. (2004). Oncomine: A cancer microarray database and integrated data-mining platform. Neoplasia 6, 1–6. doi:10.1016/s1476-5586(04)80047-2
Robey, R. B., Weisz, J., Kuemmerle, N. B., Salzberg, A. C., Berg, A., Brown, D. G., et al. (2015). Metabolic reprogramming and dysregulated metabolism: Cause, consequence and/or enabler of environmental carcinogenesis? Carcinogenesis 36 (1), S203–S231. Suppl. doi:10.1093/carcin/bgv037
Shao, H., Mohamed, E. M., Xu, G. G., Waters, M., Jing, K., Ma, Y., et al. (2016). Carnitine palmitoyltransferase 1A functions to repress FoxO transcription factors to allow cell cycle progression in ovarian cancer. Oncotarget 7, 3832–3846. doi:10.18632/oncotarget.6757
Su, P., Wang, Q., Bi, E., Ma, X., Liu, L., Yang, M., et al. (2020). Enhanced lipid accumulation and metabolism are required for the differentiation and activation of tumor-associated macrophages. Cancer Res. 80, 1438–1450. doi:10.1158/0008-5472.CAN-19-2994
Waldman, A. D., Fritz, J. M., and Lenardo, M. J. (2020). A guide to cancer immunotherapy: from T cell basic science to clinical practice. Nat. Rev. Immunol. 20, 651–668. doi:10.1038/s41577-020-0306-5
Wang, L., Li, C., Song, Y., and Yan, Z. (2020). Inhibition of carnitine palmitoyl transferase 1A-induced fatty acid oxidation suppresses cell progression in gastric cancer. Arch. Biochem. Biophys. 696, 108664. doi:10.1016/j.abb.2020.108664
Wang, Y. N., Zeng, Z. L., Lu, J., Wang, Y., Liu, Z. X., He, M. M., et al. (2018). CPT1A-mediated fatty acid oxidation promotes colorectal cancer cell metastasis by inhibiting anoikis. Oncogene 37, 6025–6040. doi:10.1038/s41388-018-0384-z
Warburg, O. (1956). On the origin of cancer cells. Science 123, 309–314. doi:10.1126/science.123.3191.309
Yang, Y., Hsu, J. M., Sun, L., Chan, L. C., Li, C. W., Hsu, J. L., et al. (2019). Palmitoylation stabilizes PD-L1 to promote breast tumor growth. Cell Res. 29, 83–86. doi:10.1038/s41422-018-0124-5
Yao, H., Lan, J., Li, C., Shi, H., Brosseau, J. P., Wang, H., et al. (2019). Author Correction: Inhibiting PD-L1 palmitoylation enhances T-cell immune responses against tumours. Nat. Biomed. Eng. 3, 414. doi:10.1038/s41551-019-0402-7
Yu, W., Lei, Q., Yang, L., Qin, G., Liu, S., Wang, D., et al. (2021). Contradictory roles of lipid metabolism in immune response within the tumor microenvironment. J. Hematol. Oncol. 14, 187. doi:10.1186/s13045-021-01200-4
Zaugg, K., Yao, Y., Reilly, P. T., Kannan, K., Kiarash, R., Mason, J., et al. (2011). Carnitine palmitoyltransferase 1C promotes cell survival and tumor growth under conditions of metabolic stress. Genes Dev. 25, 1041–1051. doi:10.1101/gad.1987211
Zhou, Y., Bollu, L. R., Tozzi, F., Ye, X., Bhattacharya, R., Gao, G., et al. (2013). ATP citrate lyase mediates resistance of colorectal cancer cells to SN38. Mol. Cancer Ther. 12, 2782–2791. doi:10.1158/1535-7163.MCT-13-0098
Keywords: lipid metabolism, CPT2, FAO, immune cell infiltration, immunotherapy
Citation: Jiang X, Du W, Shi C, Kang M, Song Q, Zhang L and Pei D (2023) Identification of a lipid metabolism-related gene for cancer immunotherapy. Front. Pharmacol. 14:1186064. doi: 10.3389/fphar.2023.1186064
Received: 14 March 2023; Accepted: 03 May 2023;
Published: 12 May 2023.
Edited by:
Karunakaran Kalesh, Teesside University, United KingdomReviewed by:
Baharan Fekry, University of Texas Health Science Center at Houston, United StatesAmrita Nandan, Independent Researcher, Prayagraj, India
Copyright © 2023 Jiang, Du, Shi, Kang, Song, Zhang and Pei. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lansheng Zhang, bGFuc2hlbmd5YW5AMTYzLmNvbQ==; Dongsheng Pei, ZHNwZWlAeHpobXUuZWR1LmNu
†These authors have contributed equally to this work