- Department of Thoracic Oncology, State Key Laboratory of Oncology in South China, Collaborative Innovation Center for Cancer Medicine, Sun Yat-sen University Cancer Center, Guangzhou, China
Background: Malignant pleural mesothelioma (MPM) is a malignant tumor originating from pleural mesothelial cells and has a high mortality rate worldwide. With the advent of immunotherapy in MPM treatment, there is an urgent need to elucidate the immune-related mechanisms in this caner.
Methods: Single-sample gene set enrichment analysis (ssGSEA) was used to score the immunocytes infiltration of data from different database sources. Identification of immunocyte-related genes was performed with weighted gene co-expression network analysis (WGCNA), differentially expressed genes (DEGs) analysis, and correlation analysis. Pan-caner analysis was performed using “DiffExp” and “Correlation” modules in TIMER.
Results: T-helper 2 (Th2) cell was found to be a poor prognostic factor for patients with MPM. Then a transcription factor, NFE2L3, was identified as a biomarker that showed a strong positive correlation with Th2 cell infiltration, and was highly expressed in MPM tissues and was related to the poor prognosis of these patients. At the same time, multiple NFE2L3 methylation sites were negatively correlated with Th2 cell infiltration, and patients with a high degree of methylation enjoy a better prognosis. Pan-caner analysis indicated that NFE2L3 might promote the differentiation of Th2 cells through the IL-2/STAT5/NLRP3 signaling pathway in MPM and many other cancers.
Conclusion: We believe that NFE2L3 can serve as a potential biomarker related to the diagnosis and prognosis of patients with MPM, and speculate that NFE2L3 could promote Th2 cell differentiation via IL-2/STAT5/NLRP3 signaling pathway in MPM and many other cancers.
Introduction
Malignant pleural mesothelioma (MPM) originates from pleural mesothelial cells and is a relatively rare type of cancer, which accounts for 0.3% of all cancer cases (Fernandez-Cuesta et al., 2021). Due to its aggressiveness and difficulty in early diagnosis, MPM is difficult to be cured. Because of the limited role of surgery in the treatment of MPM, the chemotherapy regimen of pemetrexed combined with platinum has occupied a dominant position in the treatment of MPM for a long time (de Gooijer et al., 2018). However, survival benefit from chemotherapy is limited, and the 5-years survival rate of patients with MPM is still less than 10% (Kindler et al., 2018).
In recent years, immune checkpoint inhibitors (ICIs) have been proven to improve the prognosis of various solid tumors, and their anti-tumor effects in MPM have gradually become clear (Lievense et al., 2017). A phase III clinical trial (Checkmate 743) has proved for the first time that compared with chemotherapy, the first-line treatment of nivolumab combined with ipilimumab can provide significant improvements in overall survival (OS) for patients with advanced MPM (Baas et al., 2021). However, considering the complexity of the tumor immune microenvironment, we still need to explore more immune-related mechanisms and targets to increase our understanding of MPM and treatment methods for this disease.
As an essential part of the human immune system, CD4+ T cells play a pivotal role in adaptive immune responses. Under the activation of external cytokines, naive CD4+ T cells can differentiate into multiple T helper cells, including Th1, Th2, Th17 (Zhu et al., 2010). Among them, Th2 cells have been confirmed to have immunosuppressive effects in many tumors (Kusuda et al., 2005; Nevala et al., 2009; De Monte et al., 2011). In normal tissues, T helper lymphocyte subsets are in equilibrium. However, tumor cells can secrete a variety of cytokines, putting the body in a state where Th2 cells dominate, leading to immune escape and tumor progression. Therefore, in this study, we focused on the effect of Th2 cells in patients with malignant pleural mesothelioma and explored a novel regulatory network of Th2 cell differentiation using a series of bioinformatics analysis methods.
Materials and Methods
Analysis Overview
In this study, utilizing transcriptomes downloaded from TCGA and GEO, we firstly identified that Th2 cell was associated with prognosis of patients with MPM, and NFE2L3, a transcription factor (TF), was associated with infiltration of this immunocyte. Then, using DNA methylation data from TCGA, we analyzed the correlation of NFE2L3’s methylation site with the infiltration of Th2 cell and the prognosis of patients with MPM. Meanwhile, using tissue sections from patients with MPM, we verified the expression of NFE2L3 in normal pleural tissues and tumor tissues. Finally, with the help of TIMER database, we analyzed the correlation between NFE2L3 and Th2 cell regulatory pathway-related genes in mesothelioma (MESO) and many other cancers. The detailed research process is shown in the flow chart (Figure 1).
FIGURE 1. The detailed workflow of this study. MPM, malignant pleural mesothelioma; TCGA, The Cancer Genome Atlas; ssGSEA, single-sample gene set enrichment analysis; Th2 cells, type 2 T helper cells; WGCNA, weighted gene co-expression network analysis; GEO, Gene Expression Omnibus; DEGs, differentially expressed genes.
Date Source
All data used in this study are from public databases, including Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) (Barrett et al., 2013; Cancer Genome Atlas Research et al., 2013). The clinical data, transcriptome data in FPKM format and DNA methylation data of malignant pleural mesothelioma (MPM) in TCGA were downloaded from the National Cancer Institute’s (NCI’s) Genomic Data Commons (GDC) (https://portal.gdc.cancer.gov/) (Zhang et al., 2021). Search with “mesothelioma” as a keyword, we found and downloaded gene expression profiles of GSE51024 and GSE163720 from GEO database (http://www.ncbi.nlm.nih.gov/geo/). GSE51024 is composed of 55 MPM tissues and 41 normal paired lung parenchyma tissues. And GSE163720 contains 131 tumor samples from patients with MPM.
Single-Sample Gene Set Enrichment Analysis
Immunocyte-related gene sets were got from The Cancer Immunome database (TCIA) (Charoentong et al., 2017). Through the expression value of 782 immune-related genes, we scored infiltration levels of 28 types of immunocytes. Immunocytes infiltration levels of samples were quantified by the ssGSEA algorithm in R package GSVA (Hanzelmann et al., 2013). Through Kaplan-Meier survival analysis, we found out the target immunocyte associated with the prognosis of patients with MPM.
Weighted Gene Co-Expression Network Analysis
The transcriptome data of samples from TCGA were further analyzed utilizing WGCNA to find a gene set highly correlated with infiltration level of the target immunocyte. As an algorithm for transcriptome analysis, WGCNA can identify genes with highly correlated expression patterns, and calculate the correlation between the gene set and clinicopathological traits of samples (Langfelder and Horvath, 2008).
R package “WGCNA” was used to complete the calculation process in this step. First, the 5,000 genes with the highest average expression were selected for the subsequent analysis. Then, the optimal soft-thresholding power for network construction was calculated, and module eigengenes (MEs) containing a series of co-expressed genes were constructed using a dynamic tree-cutting algorithm. Finally, the co-expressed gene set with the strongest correlation with Th2 cell infiltration can be found by analyzing the correlation between MEs and clinicopathological traits.
Functional Enrichment Analysis
Metascape, a meta-analysis website, was used for functional enrichment analysis of the co-expressed gene set (Zhou et al., 2019). By using the “Express Analysis” module in Metascape, the enriched biological processes and pathways of the selected co-expressed gene set were obtained.
Differentially Expressed Genes and Correlation Analysis
GSE51024 is composed of 96 samples, including 55 MPM tissues and 41 normal paired lung parenchyma tissues. Using GSE51024 for DEGs analysis, we can screen out genes that are significantly up-regulated or down-regulated in MPM tissues compared to normal lung parenchyma tissues. DEGs were analyzed using “limma” package in R, with an adjusted p value < 0.05 and |logFC| >1 (Ritchie et al., 2015). Mann-Whitney test was performed to calculate expression differences between tumor and normal tissues (Perme and Manevski, 2019).
The intersection of DEGs and selected MEs was used to further analyze the correlation with infiltration level of the target immunocyte. Transcriptome data and methylation data from TCGA was used for correlation analysis, and GSE163720 from GEO was used for verification. Spearman’s correlation test was performed to screen the gene with the highest correlation with infiltration level of the target immunocyte (Bishara and Hittner, 2012). And the correlation coefficient greater than 0.3 is considered vital (Funder and Ozer, 2019). Then, Kaplan-Meier survival analysis was used to explore the relationship between transcriptome data and methylation data of the screened gene and prognosis of patients with MPM.
Immunohistochemistry Staining
MPM and normal pleural tissues were obtained from patients who undergone surgery in Sun Yat-Sen University Cancer Center. All specimens were diagnosed as MPM by pathologist in Sun Yat-Sen University Cancer Center.
Immunohistochemistry (IHC) staining was used to examine NFE2L3 expression in MPM tissues and paired normal pleural tissues. All paraffin-embedded specimens were cut into 5-μm sections and placed on glass slides, then baked at 60°C for 1 h. Firstly, dewaxing all specimens with xylene and rehydrating with ethanol, then immerse the specimens in sodium citrate-EDTA buffer, and using microwave heating for antigen retrieval. Secondly, using 3% hydrogen peroxide to inactivate endogenous peroxidase, and then blocking non-specific binding with 10% goat serum. Thirdly, incubating the slides with anti-NFE2L3 rabbit polyclonal antibody (1:200; NBP2-30870; Novus Biologicals) overnight at 4°C, after washing with PBS for 4 times, then adding secondary antibody polymer horseradish peroxidase to sections. Finally, the slices wer stained with DAB (3ʹ-diaminobenzidine) and hematoxylin sequentially, dehydrated with gradient ethanol and mounted with neutral resin.
Pan-Cancer Analysis
TIMER (Tumor IMmune Estimation Resource) is a website that can provide a comprehensive analysis of transcriptome data from TCGA (Li et al., 2020). First, we study the differential expression of the screened gene between tumor and normal tissues across all TCGA samples using “DiffExp” module in TIMER. Then, using “Correlation” module, we performed correlation analysis of two related genes not only in MESO but also in many other cancer types. In this step, correlation analysis was adjusted by tumor purity.
Results
Patients With High Th2 Cell Infiltration Suffer a Poor Prognosis
The heat map of the infiltration level of 28 immunocytes is shown in Figure 2A. Through Kaplan-Meier survival analysis, it can be seen that patients with high Th2 cell infiltration suffer a significantly poor prognosis (Figure 2B, p = 0.00014). A Sankey plot presents the correlation between stage, recurrence and survival status of patients stratified by Th2 cell infiltration level (Figure 2C).
FIGURE 2. (A) Immunocyte infiltration levels of 85 samples from TCGA. (B) Kaplan–Meier curve of patients from TCGA when using the median of Th2 cell infiltration levels as the cut-off value (C) Sankey plot depicting the relationship across the tumor stage, new tumor event and survival status of patients stratified by median value of Th2 cell infiltration levels.
Identification of Co-Expressed Genes Associated With Th2 Cell
WGCNA can identify co-expressed genes that are highly correlated with Th2 cell infiltration. With the help of a dynamic tree-cutting algorithm, the 5,000 genes with the highest average expression were divided into 18 module eigengenes (MEs). Then Pearson’s correlation coefficient was used to calculate the correlation between the MEs and clinicopathological traits. Figures 3A,B shows that the red module exhibits the highest correlation coefficient with Th2 cells (Cor = 0.42, P = 5e-05). The correlation across each gene in the red module and Th2 cells is plotted in Figure 3C (Cor = 0.53, p = 3.1e-20).
FIGURE 3. Weighted gene co-expression network analysis (WGCNA) of MPM samples from TCGA. (A) Correlations between module eigengenes (MEs) and Th2 cell infiltration levels. (B) Gene significance (GS) for Th2 cell infiltration levels across all modules. (C) Gene significance (GS) of the genes contained in the red module versus module membership (Th2 cell infiltration levels). (D) Functional enrichment analysis of genes contained in the red module (Functional enrichment terms colored by p-values).
Functional enrichment analysis revealed that genes contained in the red module are mainly involved in various functions of the immune system (Figure 3D), including regulation of cell activation, neutrophil degranulation, antigen processing and presentation of peptide antigen via MHC class II, cytokine signaling in immune system, lysosome, positive regulation of immune response, leukocyte differentiation, leukocyte chemotaxis, microglia pathogen phagocytosis pathway and negative regulation of immune system process.
Identification of Th2 Cell Infiltration-Related Genes
Differentially expressed genes (DEGs) analysis was performed to identify up-regulated and down-regulated genes between MPM tissues and normal lung parenchyma tissues. Through DEGs analysis of the transcriptome data of GSE51024, 1,211 differentially expressed genes were screened, including 364 genes up-regulated in MPM tissues and 847 down-regulated genes. The expression profile of DEGs is shown in a volcano map (Figure 4A, adjusted p value <0.05 and |logFC| >1).
FIGURE 4. (A) Volcano plot of differentially expressed genes (DEGs) between tumor tissues and paired normal tissues of samples from GSE51024 (adjusted p value < 0.05 and |logFC| >1). (B) Venn plot of the intersection of DEGs from GSE51024 and genes contained in red module from WGCNA (C) Correlation between NFE2L3 expression and Th2 cell infiltration levels of samples from TCGA. (D) Correlation between NFE2L3 expression and Th2 cell infiltration levels of samples from GSE163720 (E) Correlations between DNA methylation levels of NFE2L3 related sites and Th2 cell infiltration levels of samples from TCGA.
The red module obtained by WGCNA of data from TCGA contains 260 genes, of which 10 genes are differentially expressed in MPM tissues (Figure 4B). Through correlation analysis with Th2 cell infiltration level of samples from TCGA, NFE2L3 was identified that its expression level has the highest correlation with Th2 cell infiltration level (Figure 4C, Cor = 0.47). Then we utilized the transcriptome data of GSE163720 to validate and confirmed the strong correlation between NFE2L3 and Th2 cells (Figure 4D, Cor = 0.30).
By comparing MPM tissues and normal lung parenchyma tissues of GSE51024, NFE2L3 is highly expressed in MPM tissues (Figure 5B, p = 5.04e-11). And using tissue sections for immunohistochemical staining, we can observe that NFE2L3 is mainly expressed in tumor cell nuclei, but not in normal pleural tissues (Figure 5D). It is also confirmed by analyzing the transcriptome data of TCGA that patients with higher NFE2L3 expression suffer a worse prognosis (Figure 5C, p < 0.0001).
FIGURE 5. (A) Pan-cancer analysis of NFE2L3 expression of all samples in TCGA. (B) Box plot of NFE2L3 expression in tumor tissues and paired normal tissues of samples from GSE51024 (C) Kaplan–Meier curve of patients from TCGA when using the median of NFE2L3 expression as the cut-off value. (D) Immunohistochemical staining of NFE2L3 in normal pleural epithelial tissue (black arrow) and malignant pleural mesothelioma tumor tissue (E) Kaplan–Meier curves depicting the relationships between DNA methylation levels of NFE2L3 related sites and prognosis of patients with MPM.
In TCGA, DNA methylation levels of MPM tissues were determined with the Illumina Infinium Methylation 450 K array. After matching sample names, 9 of the 19 methylation sites of NFE2L3 were found to be correlated with the infiltration of Th2 cells, and all of them were negatively correlated (Figure 4E). Table 1 details the basic information of NFE2L3 methylation sites and their correlations with Th2 cell infiltration level. Kaplan-Meier survival analysis also identified that the methylation levels of cg08822075, cg10536999, cg12510708, and cg19310148, were associated with the prognosis of patients, and patients with higher methylation levels had better prognosis (Figure 5D).
TABLE 1. The basic information of DNA methylation sites of NFE2L3 and the correlations between their methylation levels and Th2 cell infiltration. Chrom, chromosome; ChromStart, starting position in the chromosome; ChromEnd, end position in the chromosome.
NFE2L3 Could Promote Th2 Cell Differentiation via IL-2/STAT5/NLRP3 Signaling Pathway
Through pan-cancer analysis of samples in TCGA, it can be seen that compared with normal tissues, NFE2L3 is highly expressed in almost all kinds of tumor tissues (Figure 5A).
According to previous research, NLRP3 is a transcriptional regulator of Th2 cell differentiation, and signal transducer IL-2R and STAT5 triggers its expression (Bruchard et al., 2015). Therefore, we further studied the correlation between NFE2L3 and the expression of related genes involved in this pathway.
As shown in Table 2, in MESO, NFE2L3 has a strong correlation with IL-2RA (Cor = 0.295), IL-2RB (Cor = 0.303) and IL-2RG (Cor = 0.228), which constitute the high-affinity IL2 receptor (IL-2R) (Wang et al., 2005). STAT5A and STAT5B are part of JAK/STAT signaling pathway, they can mediate transcriptional signals by forming homodimers or heterodimers (Maurer et al., 2019). In MESO, NFE2L3 and STAT5B show a strong correlation (Table 2, Cor = 0.35), but the correlation between NFE2L3 and STAT5A is not statistically significant. And there is also a strong correlation between NFE2L3 and NLRP3 (Table 2, Cor = 0.254). The potential interactions across NFE2L3, IL-2R, STAT5, NLRP3, Th2 cell and tumor cells are shown in Figure 6.
TABLE 2. Correlation analyses between NFE2L3 and IL- 2R/STAT5/NLRP3 related genes in TIMER, which were adjusted by tumor purity. MESO, mesothelioma; BRCA (Her2), Her2 positive breast invasive carcinoma; DLBC, diffuse large B-cell lymphoma; ESCA, esophageal carcinoma; HNSC-HPVpos, HPV positive head and neck cancer; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LGG, low grade glioma; LIHC, liver hepatocellular carcinoma; LUSC, lung squamous cell carcinoma; PRAD, prostate adenocarcinoma; SARC, sarcoma; SKCM, skin cutaneous melanoma; THCA, thyroid carcinoma; Cor, R value of Spearman’s correlation.
FIGURE 6. The potential interactions across NFE2L3, IL-2R, STAT5, NLRP3, Th2 cell and tumor cells in the tumor microenvironment (TME).
In addition to MESO, NFE2L3 also shows powerful correlations with IL-2RA, IL-2RB and IL-2RG in many other tumors, including BRCA-Her2, DLBC, ESCA, HNSC-HPVpos, KIRC, KIRP, LGG, LIHC, LUSC, PRAD, SARC, SKCM and THCA (Table 2). And among the above cancers, NFE2L3 shows a strong correlation with STAT5A or STAT5B as well (Table 2). Among the cancers analyzed, the correlations between NFE2L3 and NLRP3 are also statistically significant, among which the correlation is relatively low in LUSC (Table 2, Cor = 0.158), and the correlation is highest in DLBC (Table 2, Cor = 0.638).
Discussion
As a very aggressive malignant tumor, malignant pleural mesothelioma (MPM) is believed to be closely related to asbestos exposure, BRCA1-associated protein 1 (BAP1) mutation and ionizing radiation to chest (Carbone et al., 2019). And a variety of somatic mutations including BAP1, TP53, NF2 and LATS1/2, are closely related to the occurrence and development of MPM (Bueno et al., 2016; Hmeljak et al., 2018; Yang et al., 2020a; Yang et al., 2020b). Due to the low incidence of MPM and the difficulty in diagnosis, there is still no unified and effective model for its treatment. Surgery used to be the only treatment for MPM, but the indications, extent of surgical resection and survival benefits of surgery are still controversial (Schipper et al., 2008; Cao et al., 2014). Beginning in 2004, the US Food and Drug Administration (FDA) approved the combination of cisplatin and pemetrexed as the first-line regimen for treatment of mesothelioma (Vogelzang et al., 2003). Although it has dramatically improved the survival of patients, the median survival period of patients who received combined chemotherapy after surgery still hovered between 17 and 25 months (Tsao et al., 2018).
Of note, with the advent of the age of immunotherapy, a variety of immune checkpoint inhibitors (ICIs) have brought a new dawn to the treatment of MPM. The DREAM study (Durvalumab with First-Line Chemotherapy in Mesothelioma) investigated the combination of PD-L1 inhibitor durvalumab and first-line chemotherapy (cisplatin and pemetrexed), and has brought survival benefits to patients (Nowak et al., 2020). Then a subsequent international, randomized, phase 3 study (CheckMate743) investigated Nivolumab in combination with Ipilimumab versus Pemetrexed with Cisplatin or Carboplatin as first line treatment in unresectable MPM, and identified that dual immune checkpoint inhibitors could bring long-term survival for these patients regardless of histological type (Baas et al., 2021).
Therefore, a better understanding of the tumor immune microenvironment and the construction of a more precise immune regulatory network will bring more individualized immunotherapy and survival benefits to patients. However, the tumor immune microenvironment is shaped by tumor cells and immunocytes together, and is in dynamic change (Schreiber et al., 2021). It is well known that CD4+ T cell populations are abundant in this environment, including pro-tumor CD4+ regulatory T cells (Tregs) and anti-tumor Th1 cells, but the role of Th2 cells is not clear so far (Facciabene et al., 2012; Tay et al., 2021).
In our study, the patients with high Th2 cell infiltration levels suffer a poor prognosis, which is consistent with some previous studies that Th2 cells are associated with tumor progression and poor prognosis in many cancers such as pancreatic cancer, breast cancer and melanoma (Nevala et al., 2009; De Monte et al., 2011; Zhang et al., 2015). The promotion effect of Th2 cells on tumors is probably due to the cytokines secreted by them, including IL-4, IL-5, IL-10 and IL-13 (Lee, 2014; Mollazadeh et al., 2019). IL-4 and IL-13 are highly similar in structure and function (Shi et al., 2021). After binding to their receptors, they can promote the proliferation, adhesion and metastasis of tumor cells, and may become potential targets for tumor treatment (Suzuki et al., 2015; Ghilardi et al., 2020). IL-5 can create a tumor-promoting immune microenvironment locally by recruiting eosinophils, thereby promoting the metastasis of tumor cells (Zaynagetdinov et al., 2015; Reichman et al., 2016). IL-10 can create an immunosuppressive tumor microenvironment through multiple pathways including NF-κB, and promote the transformation of cancer stemness (Yang et al., 2019; Saraiva et al., 2020).
Through further analysis, both the mRNA expression level and DNA methylation level of NFE2L3 was found to be highly correlated with the infiltration level of Th2 cells. NFE2L3 is a family member of the Cap’n’collar (CNC) transcription factors, and this family also include NFE2L1, NFE2L2, NF-E2, Bach1 and Bach2 (Ren et al., 2020). Among them, the family member that is widely investigated is NFE2L2, which has been confirmed as a driver gene of malignant tumor (DeNicola et al., 2011). As a homolog of NFE2L2, NFE2L3 has been proven to be related to multiple phenotypes of malignant tumors as well, including proliferation and epithelial-mesenchymal transition (EMT) (Bury et al., 2019; Ren et al., 2020). In our study, NFE2L3 was detected to be highly expressed in MPM tumor tissues, and the higher expression level is associated with poor prognosis. And the hypermethylation of multiple sites of NFE2L3 was also associated with better prognosis of patients with MPM. In multiple previous studies, somatic mutations of NFE2L2, the homolog of NFE2L3, was detected in plasma cell-free DNA (cfDNA) in hepatocellular carcinoma (HCC) and lung squamous cell carcinoma (LUSC), and was regarded as a non-invasive biomarker for tumor risk prediction and overall survival (Jeong et al., 2017; Jiao et al., 2021). Therefore, through next-generation sequencing (NGS) of tumor tissues and detection of cfDNA in peripheral blood, NFE2L3 may serve as a potential marker for the diagnosis and prognosis prediction of patients with MPM.
Then, we explored the underlying mechanism between NFE2L3 expression and Th2 cell differentiation. Th2 cell differentiation can be activated and modulated by a variety of regulators, including IL-4/STAT6 signaling pathway, IL-2/STAT5 signaling pathway and transcriptional regulator NLRP3, of which NLRP3 expression is triggered via IL-2/STAT5 signaling pathway (Lee, 2014; Bruchard et al., 2015; Jones et al., 2020). NLRP3 inflammasome is a kind of cytoplasmic protein complex, which has been proven that it can recruit myeloid-derived suppressor cells (MDSCs) and tumor-associated macrophages (TAMs) to promote tumor progression and metastasis (Moossavi et al., 2018; Hamarsheh and Zeiser, 2020). In the process of Th2 cells differentiation, NLRP3 was found to be localized in the nucleus and act as a transcription factor for Th2 cells (Bruchard et al., 2015).
In our study, the expression of NFE2L3 shows strong positive correlations with the expression of IL-2 receptor-related genes, STAT5 related genes and NLRP3. These correlations can be observed in multiple cancers, which are even more significant. Taken together, we speculated that NFE2L3, a novel biomarker in malignant pleural mesothelioma, can promote Th2 cell differentiation via IL-2/STAT5/NLRP3 signaling pathway in mesothelioma and many other cancers.
We acknowledge that there are several limitations in this study. In this article, our analysis is based only on transcriptome and DNA methylation data from TCGA and GEO without biological validation. Therefore, complete biological experiments are urgently needed to verify our conclusions in the future.
Conclusion
In our study, with the help of transcriptome data from multiple databases and a variety of bioinformatics analysis methods, we found that Th2 cell is a poor prognostic factor for patients with MPM. Through further screening, we found that NFE2L3 was highly expressed in tumor tissues of patients with MPM and both its mRNA expression level and DNA methylation level was highly correlated with the infiltration level of Th2 cells. Moreover, the correlation analysis in multiple cancers indicated that NFE2L3 was strongly correlated with the expression level of IL-2RA, IL-2RB, IL-2RG, STAT5A, STAT5B and NLRP3, which constitute the IL-2/STAT5/NLRP3 signaling pathway. Therefore, we hypothesize that NFE2L3, a novel biomarker in malignant pleural mesothelioma, may promote the differentiation of Th2 cells through the IL-2/STAT5/NLRP3 signaling pathway in multiple cancers.
Data Availability Statement
The original contributions presented in the study are included in the article Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics Statement
Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
PL conceived this study; BL and PD retrieved and disposed data from GEO database; ZW and HY performed data analysis and prepared the manuscript; PL revised the manuscript. All authors approved the final version of the manuscript.
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.
Acknowledgments
We thank the contributors to the profile of GSE51024 and GSE163720 in the GEO datasets (https://www.ncbi.nlm.nih.gov/geoprofiles).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.805256/full#supplementary-material
References
Baas, P., Scherpereel, A., Nowak, A. K., Fujimoto, N., Peters, S., Tsao, A. S., et al. (2021). First-line Nivolumab Plus Ipilimumab in Unresectable Malignant Pleural Mesothelioma (CheckMate 743): a Multicentre, Randomised, Open-Label, Phase 3 Trial. The Lancet 397 (10272), 375–386. doi:10.1016/S0140-6736(20)32714-8
Barrett, T., Wilhite, S. E., Ledoux, P., Evangelista, C., Kim, I. F., Tomashevsky, M., et al. (2013). NCBI GEO: Archive for Functional Genomics Data Sets-Update. Nucleic Acids Res. 41, D991–D995. doi:10.1093/nar/gks1193
Bishara, A. J., and Hittner, J. B. (2012). Testing the Significance of a Correlation with Nonnormal Data: Comparison of Pearson, Spearman, Transformation, and Resampling Approaches. Psychol. Methods 17 (3), 399–417. doi:10.1037/a0028087
Bruchard, M., Rebé, C., Derangère, V., Togbé, D., Ryffel, B., Boidot, R., et al. (2015). The Receptor NLRP3 Is a Transcriptional Regulator of TH2 Differentiation. Nat. Immunol. 16 (8), 859–870. doi:10.1038/ni.3202
Bueno, R., Stawiski, E. W., Goldstein, L. D., Durinck, S., De Rienzo, A., Modrusan, Z., et al. (2016). Comprehensive Genomic Analysis of Malignant Pleural Mesothelioma Identifies Recurrent Mutations, Gene Fusions and Splicing Alterations. Nat. Genet. 48 (4), 407–416. doi:10.1038/ng.3520
Bury, M., Le Calvé, B., Lessard, F., Dal Maso, T., Saliba, J., Michiels, C., et al. (2019). NFE2L3 Controls Colon Cancer Cell Growth through Regulation of DUX4, a CDK1 Inhibitor. Cel Rep. 29 (6), 1469–1481. doi:10.1016/j.celrep.2019.09.087
Cancer Genome Atlas Research Weinstein, J. N., Weinstein, J. N., Collisson, E. A., Mills, G. B., Shaw, K. R. M., Ozenberger, B. A., et al. (2013). The Cancer Genome Atlas Pan-Cancer Analysis Project. Nat. Genet. 45 (10), 1113–1120. doi:10.1038/ng.2764
Cao, C., Tian, D., Park, J., Allan, J., Pataky, K. A., and Yan, T. D. (2014). A Systematic Review and Meta-Analysis of Surgical Treatments for Malignant Pleural Mesothelioma. Lung Cancer 83 (2), 240–245. doi:10.1016/j.lungcan.2013.11.026
Carbone, M., Adusumilli, P. S., Alexander, H. R., Baas, P., Bardelli, F., Bononi, A., et al. (2019). Mesothelioma: Scientific Clues for Prevention, Diagnosis, and Therapy. CA A. Cancer J. Clin. 69 (5), 402–429. doi:10.3322/caac.21572
Charoentong, P., Finotello, F., Angelova, M., Mayer, C., Efremova, M., Rieder, D., et al. (2017). Pan-cancer Immunogenomic Analyses Reveal Genotype-Immunophenotype Relationships and Predictors of Response to Checkpoint Blockade. Cel Rep. 18 (1), 248–262. doi:10.1016/j.celrep.2016.12.019
de Gooijer, C. J., Baas, P., and Burgers, J. A. (2018). Current Chemotherapy Strategies in Malignant Pleural Mesothelioma. Transl. Lung Cancer Res. 7 (5), 574–583. doi:10.21037/tlcr.2018.04.10
De Monte, L., Reni, M., Tassi, E., Clavenna, D., Papa, I., Recalde, H., et al. (2011). Intratumor T Helper Type 2 Cell Infiltrate Correlates with Cancer-Associated Fibroblast Thymic Stromal Lymphopoietin Production and Reduced Survival in Pancreatic Cancer. J. Exp. Med. 208 (3), 469–478. doi:10.1084/jem.20101876
DeNicola, G. M., Karreth, F. A., Humpton, T. J., Gopinathan, A., Wei, C., Frese, K., et al. (2011). Oncogene-induced Nrf2 Transcription Promotes ROS Detoxification and Tumorigenesis. Nature 475 (7354), 106–109. doi:10.1038/nature10189
Facciabene, A., Motz, G. T., and Coukos, G. (2012). T-regulatory Cells: Key Players in Tumor Immune Escape and Angiogenesis: Figure 1. Cancer Res. 72 (9), 2162–2171. doi:10.1158/0008-5472.CAN-11-3687
Fernandez-Cuesta, L., Mangiante, L., Alcala, N., and Foll, M. (2021). Challenges in Lung and Thoracic Pathology: Molecular Advances in the Classification of Pleural Mesotheliomas. Virchows Arch. 478 (1), 73–80. doi:10.1007/s00428-020-02980-9
Funder, D. C., and Ozer, D. J. (2019). Evaluating Effect Size in Psychological Research: Sense and Nonsense. Adv. Methods Practices Psychol. Sci. 2, 156–168. doi:10.1177/2515245919847202
Ghilardi, N., Pappu, R., Arron, J. R., and Chan, A. C. (2020). 30 Years of Biotherapeutics Development-What Have We Learned? Annu. Rev. Immunol. 38, 249–287. doi:10.1146/annurev-immunol-101619-031510
Hamarsheh, S. A., and Zeiser, R. (2020). NLRP3 Inflammasome Activation in Cancer: A Double-Edged Sword. Front. Immunol. 11, 1444. doi:10.3389/fimmu.2020.01444
Hänzelmann, S., Castelo, R., and Guinney, J. (2013). GSVA: Gene Set Variation Analysis for Microarray and RNA-Seq Data. BMC Bioinformatics 14, 7. doi:10.1186/1471-2105-14-7
Hmeljak, J., Sanchez-Vega, F., Hoadley, K. A., Shih, J., Stewart, C., Heiman, D., et al. (2018). Integrative Molecular Characterization of Malignant Pleural Mesothelioma. Cancer Discov. 8 (12), 1548–1565. doi:10.1158/2159-8290.CD-18-0804
Jeong, Y., Hoang, N. T., Lovejoy, A., Stehr, H., Newman, A. M., Gentles, A. J., et al. (2017). Role of KEAP1/NRF2 and TP53 Mutations in Lung Squamous Cell Carcinoma Development and Radiation Resistance. Cancer Discov. 7 (1), 86–101. doi:10.1158/2159-8290.CD-16-0127
Jiao, J., Sanchez, J. I., Thompson, E. J., Mao, X., McCormick, J. B., Fisher-Hoch, S. P., et al. (2021). Somatic Mutations in Circulating Cell-free DNA and Risk for Hepatocellular Carcinoma in Hispanics. Ijms 22 (14), 7411. doi:10.3390/ijms22147411
Jones, D. M., Read, K. A., and Oestreich, K. J. (2020). Dynamic Roles for IL-2-STAT5 Signaling in Effector and Regulatory CD4+ T Cell Populations. J.I. 205 (7), 1721–1730. doi:10.4049/jimmunol.2000612
Kindler, H. L., Ismaila, N., Armato, S. G., Bueno, R., Hesdorffer, M., Jahan, T., et al. (2018). Treatment of Malignant Pleural Mesothelioma: American Society of Clinical Oncology Clinical Practice Guideline. Jco 36 (13), 1343–1373. doi:10.1200/JCO.2017.76.6394
Kusuda, T., Shigemasa, K., Arihiro, K., Fujii, T., Nagai, N., and Ohama, K. (2005). Relative Expression Levels of Th1 and Th2 Cytokine mRNA Are Independent Prognostic Factors in Patients with Ovarian Cancer. Oncol. Rep. 13 (6), 1153–1158. doi:10.3892/or.13.6.1153
Langfelder, P., and Horvath, S. (2008). WGCNA: an R Package for Weighted Correlation Network Analysis. BMC Bioinformatics 9, 559. doi:10.1186/1471-2105-9-559
Lee, G. R. (2014). Transcriptional Regulation of T Helper Type 2 Differentiation. Immunology 141 (4), 498–505. doi:10.1111/imm.12216
Li, T., Fu, J., Zeng, Z., Cohen, D., Li, J., Chen, Q., et al. (2020). TIMER2.0 for Analysis of Tumor-Infiltrating Immune Cells. Nucleic Acids Res. 48 (W1), W509–W514. doi:10.1093/nar/gkaa407
Lievense, L. A., Sterman, D. H., Cornelissen, R., and Aerts, J. G. (2017). Checkpoint Blockade in Lung Cancer and Mesothelioma. Am. J. Respir. Crit. Care Med. 196 (3), 274–282. doi:10.1164/rccm.201608-1755CI
Maurer, B., Kollmann, S., Pickem, J., Hoelbl-Kovacic, A., and Sexl, V. (2019). STAT5A and STAT5B-Twins with Different Personalities in Hematopoiesis and Leukemia. Cancers 11 (11), 1726. doi:10.3390/cancers11111726
Mollazadeh, H., Cicero, A. F. G., Blesso, C. N., Pirro, M., Majeed, M., and Sahebkar, A. (2019). Immune Modulation by Curcumin: The Role of Interleukin-10. Crit. Rev. Food Sci. Nutr. 59 (1), 89–101. doi:10.1080/10408398.2017.1358139
Moossavi, M., Parsamanesh, N., Bahrami, A., Atkin, S. L., and Sahebkar, A. (2018). Role of the NLRP3 Inflammasome in Cancer. Mol. Cancer 17 (1), 158. doi:10.1186/s12943-018-0900-3
Nevala, W. K., Vachon, C. M., Leontovich, A. A., Scott, C. G., Thompson, M. A., Markovic, S. N., et al. (2009). Evidence of Systemic Th2-Driven Chronic Inflammation in Patients with Metastatic Melanoma. Clin. Cancer Res. 15 (6), 1931–1939. doi:10.1158/1078-0432.CCR-08-1980
Nowak, A. K., Lesterhuis, W. J., Kok, P.-S., Brown, C., Hughes, B. G., Karikios, D. J., et al. (2020). Durvalumab with First-Line Chemotherapy in Previously Untreated Malignant Pleural Mesothelioma (DREAM): a Multicentre, Single-Arm, Phase 2 Trial with a Safety Run-In. Lancet Oncol. 21 (9), 1213–1223. doi:10.1016/S1470-2045(20)30462-9
Perme, M. P., and Manevski, D. (2019). Confidence Intervals for the Mann-Whitney Test. Stat. Methods Med. Res. 28 (12), 3755–3768. doi:10.1177/0962280218814556
Reichman, H., Karo-Atar, D., and Munitz, A. (2016). Emerging Roles for Eosinophils in the Tumor Microenvironment. Trends Cancer 2 (11), 664–675. doi:10.1016/j.trecan.2016.10.002
Ren, Y., Wang, Y., Hao, S., Yang, Y., Xiong, W., Qiu, L., et al. (2020). NFE2L3 Promotes Malignant Behavior and EMT of Human Hepatocellular Carcinoma (HepG2) Cells via Wnt/β-catenin P-athway. J. Cancer 11 (23), 6939–6949. doi:10.7150/jca.48100
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res. 43 (7), e47. doi:10.1093/nar/gkv007
Saraiva, M., Vieira, P., and O’Garra, A. (2020). Biology and Therapeutic Potential of Interleukin-10. J. Exp. Med. 217 (1), e20190418. doi:10.1084/jem.20190418
Schipper, P. H., Nichols, F. C., Thomse, K. M., Deschamps, C., Cassivi, S. D., Allen, M. S., et al. (2008). Malignant Pleural Mesothelioma: Surgical Management in 285 Patients. Ann. Thorac. Surg. 85 (1), 257–264. doi:10.1016/j.athoracsur.2007.06.066
Schreiber, S., Hammers, C. M., Kaasch, A. J., Schraven, B., Dudeck, A., and Kahlfuss, S. (2021). Metabolic Interdependency of Th2 Cell-Mediated Type 2 Immunity and the Tumor Microenvironment. Front. Immunol. 12, 632581. doi:10.3389/fimmu.2021.632581
Shi, J., Song, X., Traub, B., Luxenhofer, M., and Kornmann, M. (2021). Involvement of IL-4, IL-13 and Their Receptors in Pancreatic Cancer. Ijms 22 (6), 2998. doi:10.3390/ijms22062998
Suzuki, A., Leland, P., Joshi, B. H., and Puri, R. K. (2015). Targeting of IL-4 and IL-13 Receptors for Cancer Therapy. Cytokine 75 (1), 79–88. doi:10.1016/j.cyto.2015.05.026
Tay, R. E., Richardson, E. K., and Toh, H. C. (2021). Revisiting the Role of CD4+ T Cells in Cancer Immunotherapy-New Insights into Old Paradigms. Cancer Gene Ther. 28 (1-2), 5–17. doi:10.1038/s41417-020-0183-x
Tsao, A. S., Lindwasser, O. W., Adjei, A. A., Adusumilli, P. S., Beyers, M. L., Blumenthal, G. M., et al. (2018). Current and Future Management of Malignant Mesothelioma: A Consensus Report from the National Cancer Institute Thoracic Malignancy Steering Committee, International Association for the Study of Lung Cancer, and Mesothelioma Applied Research Foundation. J. Thorac. Oncol. 13 (11), 1655–1667. doi:10.1016/j.jtho.2018.08.2036
Vogelzang, N. J., Rusthoven, J. J., Symanowski, J., Denham, C., Kaukel, E., Ruffie, P., et al. (2003). Phase III Study of Pemetrexed in Combination with Cisplatin versus Cisplatin Alone in Patients with Malignant Pleural Mesothelioma. Jco 21 (14), 2636–2644. doi:10.1200/JCO.2003.11.136
Wang, X., Rickert, M., and Garcia, K. C. (2005). Structure of the Quaternary Complex of Interleukin-2 with its α, SS, and γ C Receptors. Science 310 (5751), 1159–1163. doi:10.1126/science.1117893
Yang, H., Xu, D., Schmid, R. A., and Peng, R.-W. (2020a). Biomarker-guided Targeted and Immunotherapies in Malignant Pleural Mesothelioma. Ther. Adv. Med. Oncol. 12, 175883592097142. doi:10.1177/1758835920971421
Yang, H., Xu, D., Yang, Z., Yao, F., Zhao, H., Schmid, R. A., et al. (2020b). Systematic Analysis of Aberrant Biochemical Networks and Potential Drug Vulnerabilities Induced by Tumor Suppressor Loss in Malignant Pleural Mesothelioma. Cancers 12 (8), 2310. doi:10.3390/cancers12082310
Yang, L., Dong, Y., Li, Y., Wang, D., Liu, S., Wang, D., et al. (2019). IL‐10 Derived from M2 Macrophage Promotes Cancer Stemness via JAK1/STAT1/NF‐κB/Notch1 Pathway in Non‐small Cell Lung Cancer. Int. J. Cancer 145 (4), 1099–1110. doi:10.1002/ijc.32151
Zaynagetdinov, R., Sherrill, T. P., Gleaves, L. A., McLoed, A. G., Saxon, J. A., Habermann, A. C., et al. (2015). Interleukin-5 Facilitates Lung Metastasis by Modulating the Immune Microenvironment. Cancer Res. 75 (8), 1624–1634. doi:10.1158/0008-5472.CAN-14-2379
Zhang, Q., Qin, J., Zhong, L., Gong, L., Zhang, B., Zhang, Y., et al. (2015). CCL5-Mediated Th2 Immune Polarization Promotes Metastasis in Luminal Breast Cancer. Cancer Res. 75 (20), 4312–4321. doi:10.1158/0008-5472.CAN-14-3590
Zhang, Z., Hernandez, K., Savage, J., Li, S., Miller, D., Agrawal, S., et al. (2021). Uniform Genomic Data Analysis in the NCI Genomic Data Commons. Nat. Commun. 12 (1), 1226. doi:10.1038/s41467-021-21254-9
Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A. H., Tanaseichuk, O., et al. (2019). Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat. Commun. 10 (1), 1523. doi:10.1038/s41467-019-09234-6
Keywords: NFE2L3, IL-2, Stat5, NLRP3, t-helper 2 cell, malignant pleural mesothelioma
Citation: Wang Z, Yang H, Luo B, Duan P and Lin P (2022) NFE2L3 as a Novel Biomarker Associated With IL-2/STAT5/NLRP3 Signaling Pathway in Malignant Pleural Mesothelioma and Other Cancers. Front. Genet. 13:805256. doi: 10.3389/fgene.2022.805256
Received: 30 October 2021; Accepted: 19 April 2022;
Published: 18 May 2022.
Edited by:
Eszter Lakatos, Queen Mary University of London, United KingdomReviewed by:
Haitang Yang, Shanghai Jiao Tong University, ChinaAkinobu Ota, Aichi Medical University, Japan
Copyright © 2022 Wang, Yang, Luo, Duan and Lin. 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: Peng Lin, bGlucGVuZ0BzeXN1Y2Mub3JnLmNu
†These authors have contributed equally to this work