Rongrong Liu1,2,3,4†
Weidong Liang1,5†
Qian Hua1†
Longqiu Wu1†
Xiangcai Wang1
Qiang Li1,6Fangjun Zhong1,4
Bin Li2*
Zhengang Qiu1*- 1Department of Oncology, The First Affiliated Hospital of Gannan Medical University, Ganzhou, China
- 2Department of Neurology, Ganzhou People’s Hospital, Ganzhou, China
- 3Department of Neurology, The Second Affiliated Hospital of Nanchang University, Nanchang, China
- 4Department of Neurology, The First Affiliated Hospital of Gannan Medical University, Ganzhou, China
- 5Department of Anesthesiology, The First Affiliated Hospital of Gannan Medical University, Ganzhou, China
- 6Department of Emergency, The First Affiliated Hospital of Gannan Medical University, Ganzhou, China
Introduction: Glioblastoma(GBM) is a highly malignant primary brain tumor. Even after undergoing surgery and chemotherapy, patients with this affliction still have little to no chance of survival. Current research on immunotherapy treatment for GBM shows that immune-checkpoint inhibitors (ICIs) may be a promising new treatment method. However, at present, the relationship between the fatty acid metabolic process and the prognosis of GBM patients who are receiving immunotherapy is not clear.
Methods: First, we downloaded a GBM cohort that had been treated with immunotherapy, which included the mutation and prognosis data, and the TCGA-GBM and Jonsson-GBM queues. CIBERSORT and single sample gene set enrichment analysis(ssGSEA) were used to evaluate immune cell scores. Gene set enrichment analysis (GSEA) was used to evaluate the patient’s accessment score. The pRRophetic algorithm was used to evaluate the drug sensitivity of each patient. Univariable and multivariate cox regression analyses, as well as the Kaplan-Meier (KM) method, were used to evaluate the relationship between the fatty acid metabolic process and the prognosis of GBM patients.
Results: The univariate and multivariate cox regression models showed that the fatty acid metabolic process mutant-type (MT) can be used as an independent predictor of the efficacy of immunotherapy for GBM patients. In addition, fatty acid metabolic process MT is related with significantly longer overall survival (OS) time than the wild-type(WT) variant. However, the mutation status of the fatty acid metabolic process has nothing to do with the prognosis of GBM patients who are receiving conventional treatment. Our analysis showed that fatty acid metabolic process MT correlated with significantly increased natural killer T (NKT) cells and significantly decreased CD8+T cells. At the same time, GSEA analysis revealed that fatty acid metabolic process MT was associated with significantly increased immune activation pathways and an enriched fraction of cytokine secretion compared with WT.
Conclusions: We found that fatty acid metabolic process MT may be used as an independent predictor of the efficacy of ICI treatment in GBM patients. Use of the fatty acid metabolic process MT will result in higher immunogenicity rates, a significant increase in the proportion of activated immune cells, and improvement of the immune microenvironment.
Introduction
Even if a glioblastoma (GBM), a primary brain tumor with the highest degree of intracranial malignancy, can be safely and fully resected, the median survival time for patients is still only 14 month, and there is recurrence rate of 90% (1, 2). That being said, the main treatment of early GBM is still surgical resection. Radiotherapy and chemotherapy are also important methods for treating GBM, and some clinical studies have shown that these treatments can play an active role in improving the survival time of patients. However, these methods are not without their drawbacks (3, 4). Normal brain tissue damage is caused by the low selectivity of radiotherapy and chemotherapy, while the use of chemotherapy drugs can cause systemic immunosuppression. In addition, because of the complexity and dynamic changes of tumor microenvironments, chemotherapy drugs are prone to drug resistance (4–7).
In recent years, there has been a succession of breakthroughs in the field of tumor immunotherapy. During this time, the Food and Drug Administration (FDA) has approved a number of new immunotherapy drugs for various types of tumor therapies, such as immune-checkpoint inhibitors (ICIs) and chimeric antigen receptor (CAR) T-cells (8). At present, there are many studies exploring the application effect of GBM immunotherapy. These studies have shown that in the GBM mouse model, the therapeutic inhibition of IDO, cytotoxic T lymphocyte antigen-4 (CTLA-4), or programmed death-ligand 1 (PD-L1) can significantly reduce the number of Treg cells infiltrated by tumors while significantly increasing the long-term survival rate of patients (9). Research thus far suggests that ICIs may be a promising new immunotherapy treatment method for GBM (10–12). However, the complex inhibitory immune microenvironment and the gradual formation of cancer cells during the growth of GBM severely limit the treatment’s curative effect (13). Therefore, clarifying the mechanisms behind the various aspects of the immune microenvironment would provide an important theoretical basis for better understanding the internal causes of GBM’s malignant progression while also improving GBM immunotherapy strategies (14).
In recent years, studies have shown that treatment using fatty acid metabolism can change the inflammatory microenvironment and promote the progress and metastasis of tumors (15). This interaction constitutes an intricate metabolic network in which the lipid metabolism pattern of M2 macrophages is activated by tumor cells. Zhang et al. found that M2 macrophages depend on fatty acid oxidation (FAO) pathways to obtain energy. At the same time, they also promote the secretion of IL-1β by activating reactive oxygen species and NLRP3 inflammatory corpuscles, while also accelerating the proliferation, migration, and invasion of tumor cells (16). Wu et al. found that unsaturated fatty acid derived from lipid droplets in tumor cells induced myeloid cells to polarize to M2 macrophages by promoting mitochondrial respiration. This process ultimately leads to the failure of immune cell monitoring and increased tumor growth and metastasis (17). However, despite our increasing knowledge on the subject, the role of fatty acid metabolism in the prognosis of immunotherapy with GBM is still unclear. Therefore, the purpose of this study is to explore the relationship between fatty acid metabolic pathway and the prognosis of GBM immunotherapy and condition of the immune microenvironment, in order to screen out the dominant GBM population receiving immunotherapy and further improve the clinical benefits of GBM patients.
Methods
GBM Queue
Samstein et al. published a GBM queue (ICI-GBM) that accepts ICIs, including mutation data analyzed via targeted sequencing and survival data (18). Due to the limited number of GBM queues for immune checkpoint treatment, we also collected TCGA-GBM queues (19) from the GDC database and a Jonsson-GBM queue (20) from cBioPortal. The TGCA-GBM queues included expression data, mutation data, and survival data (20). The Jonsson-GBM queue included mutation data and survival data. See Table S1 for details on the clinical features of patients in the ICI-GBM cohort.
Mutation State of the Fatty Acid Metabolic Process Pathway
In order to explore the relationship between the mutation status of the fatty acid metabolic process pathway and the clinical prognosis of GBM patients after ICI treatment, we downloaded the gene set related to the fatty acid metabolic process from the Molecular Signatures Database (MsigDB) (https://www.gsea-msigdb.org/gsea/msigdb/genesets.jsp) (21) database in the format of a GMT file. First, we screened the mutation data in the ICI-GBM, TCGA-GBMs, and Jonsson-GBM queues based on the definition of non-synonymous mutation types in the maftools R package (22). After this initial screening, only non-synonymous mutation data was retained for subsequent analysis. From here, we calculated the mutation rate of each patient’s fatty acid metabolic process and defined the mutation state of each patient’s pathway according to whether any mutation occurred. If the number of mutations in this pathway was zero, it was considered a WT. If the number of mutations in this pathway was greater than 0, it was considered an MT.
Immune Microenvironment Analysis and Drug Sensitivity Analysis
In order to explore the difference in the immune cell scores between the WT and MT groups, we used CIBERSORT and IOBR to evaluate the immune cell scores of each patient (23). In addition, we used ssGSEA and GSEA algorithms to evaluate the activation degree of patients in each pathway (24, 25). The emapplot function was used in our enrichment analysis to analyze the relationship between the pathways (26). We used the pRRophetic algorithm in our estimation of drug sensitivity to analyze the expression data of each patient (27). From this, we obtained the estimated IC50 value of each patient for 138 drugs from the Genomics of Drug Sensitivity in Cancer (GDSC) database (https://www.cancerrxgene.org/) (27).
Statistical Analysis
To explore the relationship between the fatty acid metabolic process pathway and the immunotherapy prognosis of GBM patients, we used Kaplan Meier (KM) survival analysis, univariate and multivariate cox regression analysis, and log-rank P to show the degree of statistical difference. In the difference analysis of continuous variables, the Wilcoxon-Mann-Whitney test was used to compare the amount of tumor mutation burdens (TMB), DNA damage repair (DDR) pathway mutations, and the immune cell scores. To estimate the difference of IC50 between the WT and MT groups, Fisher’s exact test was used to count the difference in the mutation frequency between each group. The difference was analyzed using the Limma Package (28). GSEA and enrichment analysis were performed using Clusterprofiler (26). A heat map was drawn using the Complexheatmap package (29). Box drawings were created using ggplot2 (30). The judging range of differential genes in the analysis was |log2FC|>1 and P < 0.05. The analysis flow of this study is shown in Figure 1.
Figure 1 An overview of this study.
Results
Prognostic Role of the Mutation State of the Fatty Acid Metabolic Process in GBM Treatment
In order to explore the relationship between the fatty acid metabolic process and the immunotherapy prognosis of GBM patients, univariate and multivariate cox regression analysis were applied to the ICI-GBM cohort. Univariate results showed that fatty acid metabolic process MT was related to significantly prolonged OS time in GBM patients (Figure 2A; P < 0.05; HR = 0.39). Through our analysis, we found that fatty acid metabolic process MT can be used as an independent predictor of ICIs in GBM patients (Figure 2B; P < 0.05; HR = 0.34). In ICI-GBM queue, survival analysis showed that fatty acid metabolic process MT had significantly longer OS time than fatty acid metabolic process WT [Figure 2C; P = 0.033; Hazard Ratio (HR) = 0.4; 95% CI: 0.22 - 0.75]. We also found that the mutation status of the fatty acid metabolic process can predict other immune queues. For example, in ICI-Bladder cancer, survival analysis showed that fatty acid metabolic process MT resulted in significantly longer OS time than fatty acid metabolic process WT (Figure 2D; P = 0.019; HR = 0.5; 95% CI: 0.31 – 0.81). In TCGA-GBM, the mutation state of fatty acid metabolic process has nothing to do with the prognosis state of GBM (Figure 2E; P > 0.05). Similarly, in Jonsson-GBM, the OS time and progression-free survival (PFS) time of patients had nothing to do with the mutation state of the fatty acid metabolic process (Figures 2F, G).
Figure 2 (A) The univariable Cox regression analysis model of ICI-GBM. (B) The multivariable Cox regression analysis model of ICI-GBM. (C) The Kaplan Meier curve of the fatty acid metabolic processes of the MT and WT groups in the ICI-GBM queues. (D) The Kaplan Meier curve of the fatty acid metabolic processes of the MT and WT groups in the ICI-BLCA queues. (E) The Kaplan Meier curve of the fatty acid metabolic processes of the MT and WT groups in the TCGA-GBM queues. (F) The Kaplan Meier curve of the fatty acid metabolic processes of the MT and WT groups in the Jonsson-GBM queues in relation to OS. (G) The Kaplan Meier curve of the fatty acid metabolic processes of the MT and WT groups in the Jonsson-GBM queues in relation to PFS. MT, mutant-type; WT, wild-type; OS, Overall Survival; PFS, Progression-free Survival.
The Mutation State of the Fatty Acid Metabolic Process and Immune Cells
Firstly, we analyzed the differences in expression data between fatty acid metabolic process MT and WT. In Figure 3A, the genes that were significantly up-regulated and significantly down-regulated in the WT and MT groups are displayed in red and blue, respectively (P < 0.05; |log2FC| > 1). Based on immune-related genes and the ssGSEA algorithm, we found that the fraction of NK cells in the MT group was significantly higher than that in WT group. On the other hand, we also found that the level of depletion among the CD8+T cells, the scores of co-inhibition antigen-presenting cells (APC), and the number of co-inhibition T cells in the MT group were significantly lower than those in WT group (Figure 3B; P < 0.05). In addition, using the CIBEROSRT algorithm to evaluate the expression data of patients, we found that the proportion of memory B cells in the MT group was significantly higher than that in the WT group (Figure 3C; P < 0.05). On the other hand, the relative ratio of resting memory CD4+ T cells in the MT group was significantly higher than that in WT group (Figure 3D; P < 0.05).
Figure 3 (A) The volcano plot of the differentially expressed genes in the TCGA-GBM queues. The red dot represents the significantly up-regulated genes (P < 0.05; log2FC > 1), while the blue dot represents the significantly down-regulated genes (P < 0.05; log2FC < -1). (B) The bar plot representing the up-regulated and down-regulated immune cells scores estimated by the ssGSEA method. The red bar represents the significantly up-regulated immune cells in the fatty acid metabolic process of the MT group, while the blue bar represents the significantly down-regulated immune cells in the fatty acid metabolic process of the MT group. The grey bar represents that there are no significantly different immune cells or pathways. (C) A comparison of the scores of the memory B cells and the fatty acid metabolic process of the MT and WT groups. (D) A comparison of the scores of the resting memory CD4+ T cells CD4 and the fatty acid metabolic processes of the MT and WT groups. MT, mutant-type; WT, wild-type; ssGSEA, Single sample Gene Set Enrichment analysis. *P < 0.05.
The Mutation State and Immunogenicity of the Fatty Acid Metabolic Process
We began by analyzing the relationship between TMB and the mutation state of the fatty acid metabolic process. We found that the level of TMB in the MT group was significantly higher than that in the WT group (Figure 4A; P < 0.05). In addition, compared with the WT group, the MT group had a high number of gene mutations in the DDR pathway (Figure 4B). We also compared the difference between MT and WT in driving mutations with the top 30 mutation frequencies in the GBM cohort (Figure 4C). We found that most of the driving mutations with significant difference belonged to oncogenes, including NF1 (25.49% vs 9.60%), ATRX (23.53% vs 7.62%), PIK3CA (17.65% vs 7.62%), KMT2C (23.53% vs 3.97%), FAT4 (25.49% vs 2.65%), TRRAP (19.61% vs 3.64%), ZFHX3 (21.57% vs 2.32%), GRIN2A (17.65% vs 2.649%), PDGFRA (13.73% vs 3.31%), SALL4 (19.61% vs 1.99%), SPEN (19.61% vs 1.99%), APC (25.49% vs 0.66%), GRM3 (11.76% vs 2.98%), PTPRT (9.80% vs 3.31%), KMT2A (17.65% vs 1.66%), SETD2 (15.69% vs 1.99%), and RANBP2 (21.57% vs 0.66%). We also compared the mutual exclusion and co-occurrence of the driver mutations with the top 30 mutation frequencies in both the MT group and WT group (Figures S1A, B).
Figure 4 (A) A comparison of the TMB levels of the fatty acid metabolic processes of the MT and WT groups. (B) A comparison of the mutation counts of DDR signaling of the fatty acid metabolic processes of the MT and WT groups. (C) The mutational landscape of the top 30 mutated driver genes in the fatty acid metabolic processes of the MT and WT groups. MT, mutant-type; WT, wild-type; TMB, Tumor mutation burdens; DDR, DNA damage repair. *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001.
The Mutation Status of the Fatty Acid Metabolic Process and Immune-Related Pathways
We used the GSEA algorithm to analyze the enrichment scores and p-values of the MT group and the WT group in combination with the expression data of GBM patients. We found that compared with the WT group, MT-GBM patients had significantly increased activity in their immune cell activation pathways (Figure 5A), such as B cell activation, lymphocyte activation, T cell activation, T cell proliferation, and lymphocyte adhesion. In addition, the immune microenvironment of the MT group was significantly enriched in the inflammatory cytokine secretion pathway (Figure 5B), as exhibited by increased IFN-γ and IL-1 secretion and production. The network diagram of enrichment analysis showed that the interaction between the immune activation pathways was very similar (Figure 5C).
Figure 5 (A) The gene sets enrichment analysis showed that the enrichment scores of the immune cells were significantly up-regulated in the fatty acid metabolic process of the MT group. (B) The gene sets enrichment analysis showed that the enrichment scores of cytokine production were significantly up-regulated in the fatty acid metabolic process of the MT group. (C) The emapplot representing the connection between the immune-related signaling pathways. MT, mutant-type.
The Mutation Status and Drug Sensitivity of the Fatty Acid Metabolic Process
In order to explore the differences in drug sensitivity between the MT group and WT group, the pRRophetic algorithm was used to estimate the drug sensitivity based on each patient’s expression data. We found that the IC50 value of methotrexate in relation to the DNA replication mechanism in the MT group was significantly higher than that in the WT group (Figure 6A). GSEA showed that the enrichment of DNA binding related pathways in the MT group was significantly lower than that in the WT group (Figure 6B). In addition, we also found that the sensitivity of MT-GBM patients to pictilisib, a PI3K/Akt inhibitor, was significantly lower than that of the WT-GBM patients (Figure 6C). GSEA also showed that the enrichment fraction of the PI3K/Akt pathway in the MT group was significantly lower than that in the WT group (Figure 6D).
Figure 6 (A) A comparison of the estimated IC50 values of methotrexate in the fatty acid metabolic processes of the MT and WT groups. (B) The results of the DNA binding signaling GSEA. (C) A comparison of the estimated of IC50 values of pictilisib in the fatty acid metabolic processes of the MT and WT groups. (D) The GSEA results of PI3K/AKT signaling. MT, mutant-type; WT, wild-type; IC50, half maximal inhibitory concentration; GSEA, Gene Set Enrichment analysis. *P < 0.05; **P < 0.01.
Discussion
In this study, univariate and multivariate Cox regression analyses showed that fatty acid metabolic process MT may be an independent prognostic factor for GBM patients receiving ICI treatment. In addition, Kaplan Meier survival analysis showed that the MT group had significantly longer OS time than the WT group. The results of immune microenvironment analysis showed that the MT group had an inflammatory immune microenvironment that manifested as significantly up-regulated NKT cells and significantly down-regulated exhausted CD8+ T cells. The results showed that the activation of immune cells, as well as the secretion and production pathways of inflammatory molecules, were significantly activated in the MT group. At the same time, the MT group displayed a significant increase in immunogenicity compared to the WT group.
The higher immunogenicity of patients with fatty acid metabolic process MT is potentially the mechanism behind patients’ better prognoses post-ICI treatment. Patients with fatty acid metabolic process MT not only had higher TMB, but also a higher number of gene mutations in their DDR pathway, which improved their immunogenicity. Due to gene fusion, deletion mutation, and point mutation of tumor genes, new abnormal proteins, which are new antigens encoded by the mutant genes of tumor cells, are produced (31). These new tumor antigens can not only enhance the immune response of T cells, but also increase the content of new antigen-specific T cells, thus improving the immune system’s ability to identify tumors (32). Studies have shown that the mutation of the DDR pathway will result in the production of higher TMB and neoantigen loads, thus significantly improving the effectiveness of ICI treatment (33–35).
An inflammatory immune microenvironment also results in greater effectiveness for ICIs treatment in patients with fatty acid metabolic process. The MT group showed a significant increase in activated immune cells such as NK cells and a significant decreased in depleted CD8+ T cells. In these groups, we found that NKT cells killed tumor cells in three ways and thus played an important anti-tumor role (1): in the Fas/FasL, perforin, and granzyme B pathways, they played a cytotoxic role and directly kill tumor cells (36) (2); they regulate the recruitment and function of other immune cells by secreting cytokines, thus indirectly playing an anti-tumor role (37) (3); NKT cells changed the immunosuppression level in the immune microenvironment, which resulted in anti-tumor activities (38). In addition, when CD8+T cells bind to antigens from dendritic cells, cytotoxic CD8+T cells can be produced. Under the induction of chemokines secreted by the dendritic cells CXCL9 and CXCL10, activated cytotoxic CD8+T cells were able to migrate to the inflammatory environment through CXCR3 expression (39). In the immune microenvironment, IL-1 helped to establish a pro-inflammatory environment by inducing pro-inflammatory cytokines and chemokines (40). IFN-γ played a central role in immunity by inducing MHCI molecules, Fas/FasL, and inducing immune proteasome-related factors (41). In the IFN-γ environment, the main cytokines secreted by activated CD8+T cells significantly increased the expression of PD-L1 (42).
This study has some limitations. Firstly, the number and sample size of clinical studies of glioma patients treated with immune checkpoint inhibitors in the database are currently small, so we hope to recruit more GBM patients who have received ICI therapy in the future to further verify the conclusions of this study. Second, few animal or cell experiments have been done on the subject. Third, the ICI-GBM cohort contained only the mutation data and clinical prognosis data of targeted sequencing, so there may have been potential bias when calculating the number of gene mutations in the fatty acid metabolic process pathway.
Conclusions
In this study, we found that fatty acid metabolic process MT may be used as an independent predictor of the efficacy of ICI treatment in GBM patients, which will allow doctors to better screen out the dominant population receiving immunotherapy in the future. Fatty acid metabolic process MT is related to higher immunogenicity, a significant increase in the quantity of activated immune cells, and improvement to the immune microenvironment.
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.
Ethics Statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author Contributions
Conceptualization, BL and ZQ. Formal analysis, RL, WL, QH, LW. Visualization, RL, QL. Writing–original draft, RL, WL, QH, LW, FZ, BL and ZQ. Writing–review & editing, RL, WL, QH, LW, FZ, BL and ZQ. All authors read and approved the final manuscript.
Funding
This work was supported by the Science and Technology Project of Education Department of Jiangxi Province (180805), the Science and Technology Project of Health and Family Planning Commission of Jiangxi Province (20195363), the CSCO-Haosen Oncology Research Fund (Y-HS2019/2-015), and the Science and Technology Project of Ganzhou City (GZ2021SF002).
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/fimmu.2022.819515/full#supplementary-material
Supplementary Figure 1 | The co-exclusivity and co-occurrence of top 30 mutated driver genes in the fatty acid metabolic processes of the WT (A) and MT (B) groups. MT, mutant-type; WT, wild-type.
Supplementary Table 1 | The clinical characteristics of the ICI-GBM cohort.
References
1. Savage N. Searching for the Roots of Brain Cancer. Nature (2018) 561:S50–1. doi: 10.1038/d41586-018-06709-2
2. Simon J-M, Toubiana T, Lang P, Taillibert S, Mazeron J-J. Radiotherapy for Glioblastomas: From Radiobiology to Concomitant Chemotherapy. Cancer Radiother (2005) 9:322–31. doi: 10.1016/j.canrad.2005.07.001
3. Lim M, Xia Y, Bettegowda C, Weller M. Current State of Immunotherapy for Glioblastoma. Nat Rev Clin Oncol (2018) 15:422–42. doi: 10.1038/s41571-018-0003-5
4. Liu F, Huang J, Liu X, Cheng Q, Luo C, Liu Z. CTLA-4 Correlates With Immune and Clinical Characteristics of Glioma. Cancer Cell International (2020) 20(1):1–0. doi: 10.1186/s12935-019-1085-6
5. Lin A, Zhang J, Luo P. Crosstalk Between the MSI Status and Tumor Microenvironment in Colorectal Cancer. Front Immunol (2020) 11:2039. doi: 10.3389/fimmu.2020.02039
6. Lin A, Wei T, Meng H, Luo P, Zhang J. Role of the Dynamic Tumor Microenvironment in Controversies Regarding Immune Checkpoint Inhibitors for the Treatment of Non-Small Cell Lung Cancer (NSCLC) With EGFR Mutations. Mol Cancer (2019) 18:139. doi: 10.1186/s12943-019-1062-7
7. Zhang J, Zhou N, Lin A, Luo P, Chen X, Deng H, et al. ZFHX3 Mutation as a Protective Biomarker for Immune Checkpoint Blockade in Non-Small Cell Lung Cancer. Cancer Immunol Immunother (2021) 70:137–51. doi: 10.1007/s00262-020-02668-8
8. June CH, O’Connor RS, Kawalekar OU, Ghassemi S, Milone MC. CAR T Cell Immunotherapy for Human Cancer. Science (2018) 359:1361–5. doi: 10.1126/science.aar6711
9. Wainwright DA, Chang AL, Dey M, Balyasnikova IV, Kim CK, Tobias A, et al. Durable Therapeutic Efficacy Utilizing Combinatorial Blockade Against IDO, CTLA-4, and PD-L1 in Mice With Brain Tumors. Clin Cancer Res an Off J Am Assoc Cancer Res (2014) 20:5290–301. doi: 10.1158/1078-0432.CCR-14-0514
10. Roth P, Valavanis A, Weller M. Long-Term Control and Partial Remission After Initial Pseudoprogression of Glioblastoma by Anti-PD-1 Treatment With Nivolumab. Neuro Oncol (2017) 19:454–6. doi: 10.1093/neuonc/now265
11. Bouffet E, Larouche V, Campbell BB, Merico D, de Borja R, Aronson M, et al. Immune Checkpoint Inhibition for Hypermutant Glioblastoma Multiforme Resulting From Germline Biallelic Mismatch Repair Deficiency. J Clin Oncol Off J Am Soc Clin Oncol (2016) 34:2206–11. doi: 10.1200/JCO.2016.66.6552
12. Johanns TM, Miller CA, Dorward IG, Tsien C, Chang E, Perry A, et al. Immunogenomics of Hypermutated Glioblastoma: A Patient With Germline POLE Deficiency Treated With Checkpoint Blockade Immunotherapy. Cancer Discov (2016) 6:1230–6. doi: 10.1158/2159-8290.CD-16-0575
13. Garner H, de Visser KE. Immune Crosstalk in Cancer Progression and Metastatic Spread: A Complex Conversation. Nat Rev Immunol (2020) 20:483–97. doi: 10.1038/s41577-019-0271-z
14. Binnewies M, Roberts EW, Kersten K, Chan V, Fearon DF, Merad M, et al. Understanding the Tumor Immune Microenvironment (TIME) for Effective Therapy. Nat Med (2018) 24:541–50. doi: 10.1038/s41591-018-0014-x
15. Long J, Zhang C-J, Zhu N, Du K, Yin Y-F, Tan X, et al. Lipid Metabolism and Carcinogenesis, Cancer Development. Am J Cancer Res (2018) 8:778–91.
16. Zhang Q, Wang H, Mao C, Sun M, Dominah G, Chen L, et al. Fatty Acid Oxidation Contributes to IL-1β Secretion in M2 Macrophages and Promotes Macrophage-Mediated Tumor Cell Migration. Mol Immunol (2018) 94:27–35. doi: 10.1016/j.molimm.2017.12.011
17. Wu H, Han Y, Rodriguez Sillke Y, Deng H, Siddiqui S, Treese C, et al. Lipid Droplet-Dependent Fatty Acid Metabolism Controls the Immune Suppressive Phenotype of Tumor-Associated Macrophages. EMBO Mol Med (2019) 11:e10698. doi: 10.15252/emmm.201910698
18. Samstein RM, Lee C-H, Shoushtari AN, Hellmann MD, Shen R, Janjigian YY, et al. Tumor Mutational Load Predicts Survival After Immunotherapy Across Multiple Cancer Types. Nat Genet (2019) 51:202–6. doi: 10.1038/s41588-018-0312-8
19. Tomczak K, Czerwińska P, Wiznerowicz M. The Cancer Genome Atlas (TCGA): An Immeasurable Source of Knowledge. Contemp Oncol (Poznan Poland) (2015) 19:A68–77. doi: 10.5114/wo.2014.47136
20. Jonsson P, Lin AL, Young RJ, DiStefano NM, Hyman DM, Li BT, et al. Genomic Correlates of Disease Progression and Treatment Response in Prospectively Characterized Gliomas. Clin Cancer Res an Off J Am Assoc Cancer Res (2019) 25:5537–47. doi: 10.1158/1078-0432.CCR-19-0032
21. Liberzon A, Subramanian A, Pinchback R, Thorvaldsdóttir H, Tamayo P, Mesirov JP. Molecular Signatures Database (MSigDB) 3.0. Bioinformatics (2011) 27:1739–40. doi: 10.1093/bioinformatics/btr260
22. Mayakonda A, Lin D-C, Assenov Y, Plass C, Koeffler HP. Maftools: Efficient and Comprehensive Analysis of Somatic Variants in Cancer. Genome Res (2018) 28:1747–56. doi: 10.1101/gr.239244.118
23. Chen B, Khodadoust MS, Liu CL, Newman AM, Alizadeh AA. Profiling Tumor Infiltrating Immune Cells With CIBERSORT. Methods Mol Biol (2018) 1711:243–59. doi: 10.1007/978-1-4939-7493-1_12
24. Hänzelmann S, Castelo R, Guinney J. GSVA: Gene Set Variation Analysis for Microarray and RNA-seq Sata. BMC Bioinformatics (2013) 14(1):1–5. doi: 10.1186/1471-2105-14-7
25. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene Set Enrichment Analysis: A Knowledge-Based Approach for Interpreting Genome-Wide Expression Profiles. Proc Natl Acad Sci USA (2005) 102:15545–50. doi: 10.1073/pnas.0506580102
26. Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, et al. Clusterprofiler 4.0: A Universal Enrichment Tool for Interpreting Omics Data. Innov (New York NY) (2021) 2:100141. doi: 10.1016/j.xinn.2021.100141
27. Geeleher P, Cox N, Huang RS. Prrophetic: An R Package for Prediction of Clinical Chemotherapeutic Response From Tumor Gene Expression Levels. PloS One (2014) 9:e107468. doi: 10.1371/journal.pone.0107468
28. Ritchie ME, Phipson B, Wu DI, Hu Y, Law CW, Shi W, et al. Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Research (2015) 43(7):e47. doi: 10.1093/nar/gkv007
29. Gu Z, Eils R, Schlesner M. Complex Heatmaps Reveal Patterns and Correlations in Multidimensional Genomic Data. Bioinformatics (2016) 32:2847–9. doi: 10.1093/bioinformatics/btw313
30. Wickham H. ggplot2: Elegant Graphics For Data Analysis [Internet]. Springer New York (2009). Available at: http://had.co.nz/ggplot2/book
31. Zhu Y, Meng X, Ruan X, Lu X, Yan F, Wang F. Characterization of Neoantigen Load Subgroups in Gynecologic and Breast Cancers. Front Bioeng Biotechnol (2020) 8:702. doi: 10.3389/fbioe.2020.00702
32. Schumacher TN, Schreiber RD. Neoantigens in Cancer Immunotherapy. Science (2015) 348:69–74. doi: 10.1126/science.aaa4971
33. Lin A, Zhang H, Meng H, Deng Z, Gu T, Luo P, et al. TNF-Alpha Pathway Alternation Predicts Survival of Immune Checkpoint Inhibitors in Non-Small Cell Lung Cancer. Front Immunol (2021) 12:667875. doi: 10.3389/fimmu.2021.667875
34. Yang Q, Zhang H, Wei T, Lin A, Sun Y, Luo P, et al. Single-Cell RNA Sequencing Reveals the Heterogeneity of Tumor-Associated Macrophage in Non-Small Cell Lung Cancer and Differences Between Sexes. Front Immunol (2021) 12:756722. doi: 10.3389/fimmu.2021.756722
35. Wang Z, Zhao J, Wang G, Zhang F, Zhang Z, Zhang F, et al. Comutations in DNA Damage Response Pathways Serve as Potential Biomarkers for Immune Checkpoint Blockade. Cancer Res (2018) 78:6486–96. doi: 10.1158/0008-5472.CAN-18-1814
36. Dashtsoodol N, Shigeura T, Tashiro T, Aihara M, Chikanishi T, Okada H, et al. Natural Killer T Cell-Targeted Immunotherapy Mediating Long-Term Memory Responses and Strong Antitumor Activity. Front Immunol (2017) 8:1206. doi: 10.3389/fimmu.2017.01206
37. Nair S, Dhodapkar MV. Natural Killer T Cells in Cancer Immunotherapy. Front Immunol (2017) 8:1178. doi: 10.3389/fimmu.2017.01178
38. Lindau D, Gielen P, Kroesen M, Wesseling P, Adema GJ. The Immunosuppressive Tumour Network: Myeloid-Derived Suppressor Cells, Regulatory T Cells and Natural Killer T Cells. Immunology (2013) 138:105–15. doi: 10.1111/imm.12036
39. Shen M, Ren X. New Insights Into the Biological Impacts of Immune Cell-Derived Exosomes Within the Tumor Environment. Cancer Lett (2018) 431:115–22. doi: 10.1016/j.canlet.2018.05.040
40. Mantovani A, Dinarello CA, Molgora M, Garlanda C. Interleukin-1 and Related Cytokines in the Regulation of Inflammation and Immunity. Immunity (2019) 50:778–95. doi: 10.1016/j.immuni.2019.03.012
41. Murata S, Takahama Y, Kasahara M, Tanaka K. The Immunoproteasome and Thymoproteasome: Functions, Evolution and Human Disease. Nat Immunol (2018) 19:923–31. doi: 10.1038/s41590-018-0186-z
Keywords: glioblastoma, immune-checkpoint inhibitors, immunotherapy, prognosis, fatty acid metabolic process
Citation: Liu R, Liang W, Hua Q, Wu L, Wang X, Li Q, Zhong F, Li B and Qiu Z (2022) Fatty Acid Metabolic Signaling Pathway Alternation Predict Prognosis of Immune Checkpoint Inhibitors in Glioblastoma. Front. Immunol. 13:819515. doi: 10.3389/fimmu.2022.819515
Received: 21 November 2021; Accepted: 10 January 2022;
Published: 18 February 2022.
Edited by:
Clare Y. Slaney, Peter MacCallum Cancer Centre, AustraliaReviewed by:
Anqi Lin, The University of Hong Kong, Hong Kong SAR, ChinaYifeng Bai, Sichuan Provincial People’s Hospital Jinniu Hospital, China
Copyright © 2022 Liu, Liang, Hua, Wu, Wang, Li, Zhong, Li and Qiu. 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: Bin Li, NTY4OTQ0MTg1QHFxLmNvbQ==; Zhengang Qiu, cWl1emhlbmdhbmdAZ211LmVkdS5jbg==
†These authors have contributed equally to this work and share first authorship