Introduction
Globally, tumors continue to pose a threat to human life and health. In 2021, statistics from the American Cancer Society showed that there were 338,090 new cases and 169,280 deaths were estimated for digestive tract tumors in the United States (Siegel et al.
2021). DSC mainly include cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), esophageal carcinoma (ESCA), liver hepatocellular carcinoma (LIHC), pancreatic adenocarcinoma (PAAD), rectum adenocarcinoma (READ), and stomach adenocarcinoma (STAD), continue to be major types of cancer threatening human health and survival. DSC are affected by several risk factors, including lifestyle habits, tobacco use, alcoholism, physical activity, and gender (Xie et al.
2021). Therefore, it is of great importance to identify a biological target with ideal diagnostic and therapeutic importance for patients with tumors.
NFE2L3 was discovered in the 1990s, belonging to the cap 'n' collar basic-region leucine zipper family, mainly located in the endoplasmic reticulum and the nuclear envelope Chevillard and Blank (
2011). Currently, researches related to ap 'n' collar basic-region leucine zipper family primarily focus on nuclear factor erythroid-derived 2-like 2 (NFE2L2). NFE2L2 plays the role of the primary transcriptional regulator for the cell's antioxidant response controls, which are associated with cancer and drug resistance. The binding of kelch-like ECH-associated protein 1 (KEAP1) to NFE2L2 causes polyubiquitylation, leading to swift degradation of NFE2L2 and significantly influencing the initiation, progression, and evolution of tumors. The mechanism of NFE2L3 in cancers is entirely different from NRF2, but limited research has evaluated the significance of NFE2L3 as an innovative indicator in malignancies.
NFE2L3, as a crucial regulator of the cellular stress response, expressed in diverse tissues of the human, different levels of expression highlight the unique functionality of NFE2L3 in different tissues (Qian et al.
2022; Chen et al.
2019; Zhang et al.
2022). In recent years, a growing number of evidence has indicated the participation of NFE2L3 in the processes of differentiation, inflammation, and carcinogenesis. James Saliba demonstrated that the deletion of NFE2L3 results in decreased inflammation and a shift in the tumor microenvironment via the IL33 and RAB pathways (Saliba et al.
2022). Knockdown of NFE2L3 inhibited bladder cancer cell proliferation by inducing the cell cycle and apoptosis (Qian et al.
2022). In vitro and in vivo experiments support NFE2L3 as an oncogene involved in the process of tumor formation. Although findings have suggested that NFE2L3 may be a crucial regulator of cancer progression, the underlying molecular mechanisms remain poorly understood.
The prevalent approach to treating cancer involves a combination of chemotherapy, radiotherapy, and surgical excision of the tumor. Tumor immunotherapy as a novel therapeutic modality has revolutionized cancer treatment (Galluzzi et al.
2020; Wang et al.
2020; Feng et al.
2021). The investigators discovered that NFE2 family is a potential prognostic biomarker and is associated with immune infiltration in many cancers (Ju et al.
2020).
The goal of the study was to investigate the possibility of NFE2L3 as novel biomarker for DSC. We import data from different database to provide a systematic and comprehensive insight into NFEE2L3 expression level, diagnostic and prognostic value, and underling mechanism in DSC. We also used patients’ specimens and in vivo models test to confirm the biological functions of NFE2L3 in DSC.
Materials and methods
Gene expression analysis
All raw data were downloaded from The Cancer Genome Atlas (TCGA), University of California Santa Cruz (UCSC) Xena browser (
https://xena.ucsc.edu/), and the Human Protein Atlas (HPA) (
https://www.proteinatlas.org/). Raw data were processed using Strawberry Perl 5.30.0.1 and R 4.1.1. The Wilcoxon test was used to analyze the data.
p < 0.05 (*) was considered statistically significant.
Twelve gastric cancer and normal tissues from patients admitted to Second Affiliated Hospital of Nanchang University. Written informed consent was obtained from the patients. This study was conducted in accordance with the guidelines of the Declaration of Helsinki, and approved by the Medical Ethics Committee of the Second Affiliated Hospital of Nanchang University (No. 2022096 Approve data: 2022.9.7). The recruitment period for this study from October 1, 2022 to January 30, 2023.
Tumor mutational burden and microsatellite instability of NFE2L3 in pan-cancer
We extracted relevant data of Tumor mutation burden (TMB) and Microsatellite Instability (MSI) from raw data, and the relationship between NFE2L3 and pan-cancer was analyzed using Spearman’s rank correlation coefficient.
Clinical relevance of NEF2L3 in DSC
We created receiver operating characteristic (ROC) curves using the R package “pROC” to examine the diagnostic value of NEF2L3 in DSC. The ROC curve area (AUC) between 0.7 and 0.9 indicated sure precision, and the AUC over 0.9 indicated great accuracy. The prognostic effect of NFE2L3 in DSC was investigated by performing Kaplan–Meier survival analysis. We selected three markers for the Kaplan–Meier analysis: overall survival (OS), disease-specific survival (DSS), progression-free interval (PFI). Both the “survival” and “survminer” packages were used to plot the survival curves. The relationship between NFE2L3 expression and clinical features is studied using “limma” and “ggpubr” packages.
NFE2L3 expression in subtypes and immune infiltration
NFE2L3 expression in immunological or molecular subtypes was measured in TISIDB (
http://cis.hku.hk/TISIDB/), which was a website for tumor immune interaction. The immune cell composition of DSC was analyzed using CIBERSORT (
https://cibersort.stanford.edu/). R packages “ggplot2”, “ggpubr”, and “ggExtra” were used to process the data and calculated the correlation coefficient. |R|≥ 0.3 and
p < 0.05 (*) were considered statistically significant.
Co-expression genes of NFE2L3 in DSC
We used the packages “limma” and “BiocManager” to investigate the co-expressed genes and chemokine receptors that were both positively and negatively related to NFE2L3 expression in DSC.
Protein–protein interaction network building
A protein–protein interaction (PPI) network was constructed using Search Tool for the Retrieval of Interacting Genes (STRING) (
https://string-db.org/), with the following fixed principal parameters: a minimum required interaction score of 0.4, a maximum of 50 inter-actors to be displayed, and the sources of active interaction from experiments, text mining, and databases. The PPI network was visualized using Cytoscape (version 3.9.0).
Gene ontology and kyoto encyclopedia of genes and genomes pathway enrichment analysis
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses and Genomes Enrichment Analyses (GSEA) were performed to assess the potential biological functions and pathways of NFE2L3 in DSC, using the packages “clusterProfiler”, “ggplot2” and “limma”.
qRT-PCR
Total RNA extraction was extracted using the TRIzol Up Plus RNA Kit (Trnsgen). The cDNA synthesis of cDNA is done using M-MLV reverse transcription (Trnsgen). Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed by SYBR Green PCR Master Mix (Takara).
The primer sequences of NFE2L3 for qRT-PCR were as follows:
forward, 5′-TTCAGCCAGGCTATAAGTCAGG-3′,
reverse, 5′-GCTCAGGATTGGTGGTATGAGA-3′.
The primer sequences of β-actin for qRT-PCR were as follows:
forward, 5′-CTCCATCCTGGCCTCGCTGT-3′,
reverse, 5′-GCTGTCACCGTTCC-3′.
Cell line and culture
HGC-27 and AGS were purchased from the American Type Culture Collection (ATCC). Cells were cultured in RPMI-1640 containing 10% fetal bovine serum and cultured at 37℃ containing 5% CO2.
siRNAs and plasmids transfection
Plasmids (pcDNA3.1) were purchased from Shanghai Genechem. Plasmids were transfected into AGS via Lipofectamine 2000 (ThermoFisher). siRNA (siNFE2L3#1, siNFE2L3#2, siNFE2L3#3) purchased from RIBOBIO and transfected into HGC-27 using Lipofectamine 2000. qRT-PCR was applied to detect the efficiency of NFE2L3 overexpression or knockdown.
The sequences of siRNA were as follows:
siNFE2L3#1 GCATGACCTTTATCATGAT
siNFE2L3#2 GTATGACCTTGACATAAAT
siNFE2L3#3 GTCAGGATGTGAATCTTCA
Proliferation and migration assays
For proliferation experiments, approximately 1000 cells per well for each cell line were seeded in 96-well plates. The 450 nm absorbance was tested for 6 days (0, 1, 2, 3, 4, 5 days). The CCK-8 reagent was purchased from APExBIO (USA). We seeded 3 × 105 cells per well in 6-well plate and collect images at 0 h, 12 h, and 24 h to calculate migration rates over time.
Statistical analysis
The statistical analysis in this study was automatically calculated using the aforementioned online databases. p value < 0.05 was checked as the statistical significance. *p < 0.05, **p < 0.01, ***p < 0.001.
Discussion
Our work in this study provides a comprehensive insight into the importance of NFE2L3 in DSC and revealed that NFE2L3 has great potential to become a widely utilized biomarker. NFE2L3 is found to be abnormally expressed in cancers and regulates DSC through multiple mechanisms. NFE2L3 has ideal clinical relevance for its good prognostic and diagnostic value.
Accumulating evidence suggests that NFE2L3 is crucial for the development of cancer. NFE2L3 belongs to the Cap ‘
n’ Collar family, which plays a vital role in the cancer-promoting function by various mechanisms (Zhang et al.
2009; Sekine and Motohashi
2021). For example, NFE2L3 functions as an oncogene by triggering the Wnt/-catenin signaling pathway (Ren et al.
2020), absence of the NFE2L3 affects the formation and growth of tumors (Bury et al.
2019a,
b). NFE2L3 regulates the immune microenvironment and has potential prognostic value in patients with kidney cancer (Feng et al.
2021). The absence of NFE2L3 results in significantly less inflammation in the colon, which inhibits an increase in the number and volume of tumors (Ju et al.
2020). Although the close relationship between NFE2L3 and cancers has been demonstrated, the clinical value and the underlying mechanism of NFE2L3 in cancers remain unclear.
A pan-cancer analysis showed that NFE2L3 is abnormally expressed in most malignancies, predominantly in the DSC. These findings have also been confirmed in cell lines and in cancer tissue data obtained from the HPA database. We further confirm this speculation in gastric cancer tissues and cell lines, the qRT-PCR results showed significant levels of NFE2L3. The above results lead to an important conclusion that NFE2L3 may involve in promoting tumorigenesis.
While prognostic clinical markers are useful, existing molecular makers are still of limited value in cancers. Recent studies suggested that NFE2L3 has an ideal diagnostic and prognostic value in bladder, colorectal, and renal cancers (Feng et al.
2021; Waku et al.
2020). Consistent with our findings, NFE2L3 functioned as an oncogene in PAAD and LIHC based on the fact that patients in the high-expression group had a worse prognosis than those in the low-expression group. Our results also confirmed that NFE2L3 showed ideal diagnostic and prognostic value for DSC. Therefore, NFE2L3 is expected to emerge as a novel indicator for the diagnosis and prognosis of DSC.
It is well established that cancer is a complex and heterogeneous genomic disease, and gene functions cannot be guaranteed based on the basis of standard definitions. One explanation for the aberrant function of genes in cancer is that different cellular contexts result in a wide range of pathological outcomes (Datta et al.
2020; Lin et al.
2020). A detailed mechanistic explanation for this problem has not yet been fully elucidated, and it is worth investigating in future.
ICIs as a mainstay of cancer treatment have provided novel strategies and improved clinical outcomes in multiple cancer types. However, response rates to immunotherapy vary due to several critical factors, including mutation burden, tumor-infiltrating lymphocytes, and regulatory checkpoint receptors (Carlino et al.
2021)
. Researchers found that NFE2L3 played an important role in various immune-related pathways and immune cells. In Nfe2l3-/- mice, the levels of FOXP3 and the immune checkpoint markers cytotoxic T lymphocyte-associated protein 4 (CTLA4), TIM3, and lymphocyte activating 3 (LAG3) were increased (Bury et al.
2019a,
b). The enrichment analyses revealed that NFE2L3 was associated with various immune-relevant pathways in kidney cancer and related to the infiltration ratios of 17 types of immune cells in patients with kidney cancer (Zhang et al.
2022). Nevertheless, the molecular mechanisms behind it remain mysterious.
The results of our study indicated a strong correlation between NFE2L3 expression and MIS in DSC. MSI is a certain phenomenon that often arises in tumors that have an underlying defect in DNA mismatch repair (Tomlinson et al.
2002). Patients with high MSI rates indicating poor median PFS or OS (De' Angelis et al.
2018; Ganesh et al.
2019; Asaoka et al.
2015). Therefore, it can be inferred that MSI could play a crucial role in shaping the progression of DSC. Relevant studies have indicated that cancers harbor high TMB means the potential to generate many neoantigens, making tumors more immunogenic for ICIs (Bravaccini et al.
2021; Chan et al.
2019). Expect MSI, the TMB also demonstrated the significance of the DSC. A higher TMB is normally related to an adverse outcome in advanced tumors.
According to our results, we revealed that NFE2L3 was closely correlated with immunity. It is expected to become novel immunotherapeutic target for DSC in future. Further studies should continue to investigate the mechanisms of NFE2L3 in immunotherapy.
Our findings indicated a significant positive correlation between NFE2L3 and immune or molecular subtypes. Cellular and immune subtypes are the main factors responsible for tumor heterogeneity. Treatment elicits varied responses from tumors, owing to the genetic diversity within the tumor. The combination of immune and molecular subtypes for precisely targeted tumor therapy will also greatly enhance the current status of immunotherapeutic efficacy. Therefore, we believe that NFE2L3 plays a vital role in DSC immunotherapy. Therapeutic cancer vaccines remain a valid immunotherapy option and hold significant implications for the clinical advancement of cancer patients (Pe'er et al.
2021; Baharom et al.
2022). Over the past few years, there has been a significant amount of research conducted on vaccines and therapeutics specifically designed for the innate and adaptive immune systems in the field of cancer treatment (Peng et al.
2022; Ratnam et al.
2019). Thus, the important functions of immune cells in tumor-specific immunity against immune cells can be understood. Our findings showed that compared with non-responders, patients who generated antibody responses to therapy exhibited prolonged DFS. We have reason to believe that NFE2L3 to be a potential target for immunotherapy in future.
A multitude of studies have indicated that patients with immune cell infiltration tumors exhibit a more favorable response to ICI than patients with non-immunological tumors. Recent research has revealed that NFE2L3 plays a pivotal role in cancer development (Wang et al.
2018; Zhang et al.
2019), but the precise mechanisms have yet to be elucidated. The findings of this research indicated that DSCs experienced immune cell infiltration. The immune cell infiltration encompasses a wide range of cell types, such as Tregs, T cells, dendritic cells, NK cells, T follicular helper cells, macrophages, mast cells, and monocytes. In STAD, we observed a strong correlation between NFE2L3 expression and immune cells abundance. The research revealing the correlation of immune cell infiltration with NFE2L3 may open unique avenues for precise cancer therapeutics against cancers.
Co-expression analysis can be employed as a means to identify disease-associated genes and gene functions in tumors. Chemokines combine with the cell surface chemokine receptors to perform biological functions, such as the chemotaxis, leukocyte migration, and inflammatory activities (Charo and Ransohoff
2006; Arimont et al.
2017). The causative role of chemokines is variable, and they play an important role in many diseases, such as cancer, viral infections, inflammatory, and autoimmune diseases. Researchers have found that chemokines affect anti-tumorigenic activity by regulating tumor angiogenesis and infiltration of immune cells (Lacalle et al.
2017). The results of our study showed a strong correlation between NFE2L3 and immune-related genes or chemokine receptors. In COAD, LIHC, and STAD, there was a close correlation between NFE2L3 expression and the chemokine receptors, such as C-X3-C motif chemokine receptor 1 (CX3CR1), plexin B2 (PLXNB2), atypical chemokine receptor 1 (ACKR1), and formyl peptide receptor 2 (FPR2). Anja Schmall found that the tumor-associated crosstalk between macrophages and cancer cells via the CCR2 and CX3CR1 signaling pathways directed the lung cancer growth and metastasis (Schmall et al.
2015). Circ_0013958 plays an oncogenic role in ovarian cancer by regulating the miR-637/PLXNB2 axis (Liang et al.
2021). Immune-related genes, such as PD-1, PD-L1, and CTLA4, have been successfully used in tumor immunotherapy and have been shown to have a significant effect on B cell lymphomas, HCC, and non-small cell lung cancer (Xu-Monette and Zhou
2018; Kamada et al.
2019), (Thommen et al.
2018). This study revealed novel immune-related genes that are co-expressed with NFE2L3 in DSC, including TIGIT, IDO1, ADORA2A, and CD70, which can be used as immune therapy, checkpoints, and biomarkers (Zhai et al.
2018; Chauvin and Zarour
2020). Therefore, combining immune-related genes or chemokine receptors with NFE2L3 may enhance the efficacy of DSC diagnosis and immunotherapy.
We conducted GO enrichment and KEGG analysis on the NFE2L3 binding proteins, we found that the Hedgehog signaling pathway is the primary target of protein catabolic processes and cell differentiation. Hedgehog signaling pathway activated tumors by driving EMT, and the inhibitor of Hedgehog signaling pathway exhibited remarkable clinical outcomes in various types of cancer Zhang et al. (
2009). Nevertheless, the investigation into the differential regulation of target genes by NFE2L3 remains to be investigated.
NFE2L3 is involved in the regulation of various biological and cellular processes, such as cell cycle, cell differentiation, or inflammatory processes. We performed GSEA and found that NFE2L3 was significantly associated with epigenetics, methylation modification, cell cycle regulation, calcium signaling pathway and steroid hormone biosynthesis. Consequently, the distinct manifestation of NFE2L3 assumes a crucial function in DSC through various biological functions and signaling pathways, ultimately impacting the process of tumorigenesis. Nevertheless, a comprehensive understanding of the findings would require further experimental validation.
To validate the bioinformatics findings, in vivo experiments conducted on gastric cancer cells demonstrated the role of NFE2L3 as an oncogene, facilitating the proliferation and movement of gastric cancer cells. NFE2L3 enhances tumor cell migration ability by affecting the EMT through Wnt/β-catenin signaling pathway (Ren et al.
2020; Zeng and Ju
2018). These findings and our hypothesis are validated experimental results. More in vitro research is needed to gain a better understanding of how things work in future.
The outcomes of these bioinformatics analyses lay the groundwork for future comprehensive investigations into the mechanisms underlying tumorigenesis and evolution. NFE2L3 was found to be aberrantly expressed in DSC, and significantly correlated with the diagnosis and prognosis of DSC patients. In addition, in vivo results suggest that NFE2L3 may act as an oncogene roles in gastric cancer through promoting the gastric cancer cells growth and migration.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.