Skip to main content

ORIGINAL RESEARCH article

Front. Pharmacol., 12 May 2023
Sec. Pharmacogenetics and Pharmacogenomics
This article is part of the Research Topic Integrating Genetics and Proteomics for Drug Discovery View all articles

Identification of a lipid metabolism-related gene for cancer immunotherapy

Xin Jiang&#x;Xin Jiang1Wenqi Du,&#x;Wenqi Du1,2Ce Shi&#x;Ce Shi3Mengjie KangMengjie Kang1Qiuya SongQiuya Song1Lansheng Zhang
Lansheng Zhang4*Dongsheng Pei
Dongsheng Pei1*
  • 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.

FIGURE 1
www.frontiersin.org

FIGURE 1. The analytical process of exploring the function of CPT2 in the pan-cancer landscape.

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
www.frontiersin.org

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
www.frontiersin.org

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 (∗: pvalue<0.05; ∗∗: pvalue<0.01; ∗∗∗: pvalue<0.001). As displayed in gray columns when normal data are available. (D) TNMplot database data concerning the expression pattern of CPT2 in human cancer. Red boxes represent tumor tissues; green boxes represent normal tissues. ∗ P<0.05, ∗∗ P<0.01.

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
www.frontiersin.org

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
www.frontiersin.org

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
www.frontiersin.org

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
www.frontiersin.org

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Warburg, O. (1956). On the origin of cancer cells. Science 123, 309–314. doi:10.1126/science.123.3191.309

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, J., Du, K., Li, S., Lu, L., Mei, J., Lin, W., et al. (2021). Glutamine metabolism regulators associated with cancer development and the tumor microenvironment: A pan-cancer multi-omics analysis. GenesGenes (Basel) 12, 1305. doi:10.3390/genes12091305

CrossRef Full Text | Google Scholar

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 Kingdom

Reviewed by:

Baharan Fekry, University of Texas Health Science Center at Houston, United States
Amrita 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

Disclaimer: 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.