- 1Department of Clinical Oncology, Taihe Hospital, Jinzhou Medical University Union Training Base, Shiyan, China
- 2Department of Clinical Medicine, The First Clinical College of Hubei University of Medicine, Shiyan, China
- 3Department of Clinical Oncology, Taihe Hospital, Hubei University of Medicine, Shiyan, China
- 4Hubei Key Laboratory of Stem Cell Research, Taihe Hospital, Hubei University of Medicine, Shiyan, China
- 5Department of Oncology, Danjiangkou First Hospital, Danjiangkou, China
- 6Biomedical Research Institute, Hubei University of Medicine, Shiyan, China
Background: Glucose-6-phosphate isomerase (GPI) plays an important role in glycolysis and gluconeogenesis. However, the role of GPI in lung adenocarcinoma (LUAD) remains unclear.
Methods: All original data were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases and integrated via R 3.2.2. GPI expression was explored with TCGA, GEO, and Oncomine databases. Immunohistochemistry staining was used to analyze GPI expression in clinical specimens. The correlations between GPI and cancer immune characteristics were analyzed via the TIMER and TISIDB databases. GPI-specific siRNAs were used to verify the role of GPI expression on cell proliferation and cell cycle distribution.
Results: In general, GPI is predominantly overexpressed and has reference value in the diagnosis and prognostic estimation of LUAD. Upregulated GPI was associated with poorer overall survival, clinical stage, N stage, and primary therapy outcome in LUAD. Mechanistically, we identified a hub gene that included a total of 56 GPI-related genes, which were tightly associated with the cell cycle pathway in LUAD patients. Knockdown of GPI induced cell proliferation inhibition and cell cycle arrest. GPI expression was positively correlated with infiltrating levels of Th2 cells and regulatory T cells (Tregs); in contrast, GPI expression was negatively correlated with infiltrating levels of CD8+ T cells, central memory T cells, dendritic cells, macrophages, mast cells, and eosinophils. GPI was negatively correlated with the expression of immunostimulators, such as CD40L, IL6R, and TMEM173, in LUAD.
Conclusion: GPI may play an important role in the cell cycle and can be used as a prognostic biomarker for determining the prognosis and immune infiltration in LUAD.
Introduction
Lung cancer is one of the malignant cancers with the highest incidence and the worst prognosis worldwide. Lung cancer is divided into small cell lung cancer and non-small-cell lung cancer (NSCLC) (1, 2). Lung adenocarcinoma (LUAD) is the largest subgroup of NSCLC (3). Recently, new therapies have been created for NSCLC, including surgery, chemotherapy, radiotherapy, targeted therapy, and immunotherapy (4). However, the recurrence rate and mortality rate of LUAD are still high, and the prognosis of patients is poor. Therefore, it is imperative to find more effective biomarkers to facilitate novel therapeutic methods.
Immunotherapy has achieved significant results in a variety of tumors and is expected to improve the prognosis of NSCLC (5). However, whether the interaction between NSCLC and the immune system could affect the prognosis of NSCLC patients is largely unclear. Tumor-infiltrating lymphocytes (TILs) have been shown to play an important role in tumor-associated immune responses (6). A number of studies have shown that TILs can effectively improve the prognosis of various tumors (7, 8). Therefore, it is of great clinical significance to investigate the infiltration of immune cells in NSCLC.
Glucose-6-phosphate isomerase (GPI) is a member of the glucose phosphate isomerase protein family. GPI plays an important role in glycolysis and gluconeogenesis, interconverting D-glucose-6-phosphate and D-fructose-6-phosphate (9, 10). GPI can be secreted to the outside of cells, functioning as a cytokine or growth factor (11). GPI is an important autoantigen in human rheumatoid arthritis (RA), and GPI levels, including soluble GPI and GPI immune complexes, are significantly increased in serum and synovial fluid in RA patients (11). In the past decade, aberrant expression of GPI has also been involved in the development and progression of several cancers (12–14). An increasing number of studies have shown that GPI is not secreted by normal cells, but GPI levels are elevated in the serum or urine of patients with malignant tumors. A high expression of GPI in glioblastoma and renal cell carcinoma is associated with poor prognosis. GPI promotes the migration of melanoma cells and the development of invasive breast cancer stem cells and may be involved in cancer metastasis and invasion (9, 15). However, the role of GPI in LUAD has not been systematically explored in multiple aspects.
In recent years, a growing number of platforms and databases have enabled cancer researchers to use multiple sets of data for bioinformatics analysis of cancer. To better understand the role of the GPI gene in NSCLC, we comprehensively evaluated the relationship between GPI expression and the prognosis of LUAD patients through The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), Oncomine, and Kaplan–Meier plotter databases. In addition, we also studied the correlation between GPI and TILs using the Tumor Immunity Estimation Resource (TIMER) and TISIDB databases. Our results demonstrate that GPI can be used as a biomarker to predict the prognosis and immune infiltration of LUAD patients.
Materials and Methods
Patient Datasets
The messenger RNA (mRNA) expression data (including 535 LUAD samples and 59 adjacent nontumor samples) and clinical information were downloaded from TCGA database (https://cancergenome.nih.gov). We also downloaded the following gene expression profiles from the GEO: GSE10072 (including 33 LUAD samples and 33 paired adjacent nontumor samples), GSE32863 (including 58 LUAD samples and 58 paired adjacent nontumor samples), GSE7503 (including 83 LUAD samples and 83 paired adjacent nontumor samples), GSE30219 (including 84 LUAD samples and 14 adjacent nontumor samples), GSE31210 (including 226 LUAD samples), GSE37745 (including 106 LUAD samples), and GSE3141 (including 58 LUAD samples). These were extracted from the GEO database (https://www.ncbi.nlm.nih.gov/geo/) to further validate our results. The Human Protein Atlas (HPA) database (http://www.proteinatlas.org/) was used to verify the expression of GPI in LUAD at the protein level. We obtained the gene amplification and mutation status of GPI through the cBioPortal for Cancer Genomics (http://www.cbioportal.org/).
Oncomine Database Analysis
The expression level of the GPI gene in various types of cancers was identified in the Oncomine database (https://www.oncomine.org/resource/login.html) (16). Student’s t-test was used to compare the transcription levels of the GPI gene in cancer samples with those in normal controls. The cutoff p-value and fold change were defined as 0.05 and 1.5, respectively.
Survival Analysis
The Kaplan–Meier plotter (http://kmplot.com/analysis/) (17) and GEPIA were used to estimate the correlation between GPI expression and the survival rate of different clinical features in LUAD patients, and the hazard ratio (HR) and log-rank p-value of the 95% confidence interval were calculated.
Protein–Protein Interaction Network Analysis
Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) interaction (18), we constructed a pathway and gene regulatory network. Protein–protein interaction (PPI) data were extracted from the STRING database based on protein interactions and signaling pathways (https://string-db.org), and the network was constructed using Cytoscape 3.7.2 applications.
Establishment and Evaluation of the Nomograms for LUAD Survival Prediction
In this study, we selected all independent clinicopathological prognostic factors from Cox regression analysis and constructed a contingency table to evaluate the 1-, 3-, and 5-year overall survival (OS) probability of LUAD patients. By comparing the predicted probability of the line chart with the observed actual probability through a calibration curve, the accuracy of the line chart was verified. The overlapping reference lines show that the model is accurate. Receiver operating characteristic (ROC) analysis was used to compare the prediction accuracy of the line chart of the combined model and the line chart of various clinicopathological prognostic factors.
Functional Enrichment Analysis
To understand the biological processes and pathways that GPI may participate in, we conducted the following analysis. The first 600 genes related to GPI in LUAD and 56 genes related to LUAD survival were obtained from GEPIA (http://gepia.cancer-pku.cn/) (19). These genes were enriched by Gene Ontology (GO) [including biological processes (BP), cellular components (CC), and molecular function (MF)] and KEGG pathway analyses using the Database for Annotation, Visualization and Integrated Discover (DAVID) and visualized by the R package ggplot2. A corrected p < 0.05 was determined to be statistically significant.
LinkedOmics Database Analysis
The LinkedOmics database (http://www.linkedomics.org) (20) contains multiomics data and clinical data for 32 cancer types and a total of 11,158 patients from TCGA project (20). In this study, the LinkedOmics database was used to explore the expression profile of GPI. We explored the GO biological process (GO_BP) and KEGG pathways of GPI and its co-expressed genes by using gene set enrichment analysis (GSEA) in the Link Interpreter module.
Timer Database Analysis
TIMER (http://www.cistrome.shinyapps.io) (21) is a convenient and accurate online analysis tool that can explore the expression levels of genes in normal and tumor tissues in multiple tumors in TCGA datasets, infer the abundance of tumor-infiltrating immune cells from the gene expression profiles, and evaluate their clinical impact. In this study, the TIMER database was used to analyze the expression levels of genes in normal and tumor tissues in multiple tumors. The deconvolution statistical method in TIMER was used to infer the abundance of tumor-infiltrating immune cells (TIICs) from the gene expression profiles of LUAD samples in TCGA datasets.
TISIDB Database Analysis
TISIDB is an integrated knowledge base portal that plays a vital role in detecting the interaction between tumors and the immune system (http://cis.hku.hk/TISIDB/) (22). To further clarify the immune correlation of GPI in cancer, we used the “Immunomodulator” module of the TISIDB database to analyze and evaluate the correlation between GPI expression and the levels of immune checkpoint genes. To further study the association between GPI and chemokine/chemokine receptor expression, we evaluated the expression levels of chemokine/chemokine receptors of TIICs through the “chemokine” module.
Cell Culture and Transfection
The lung adenocarcinoma cell lines H1975 and A549 were purchased from ATCC and cultured in Dulbecco’s modified Eagle’s medium (DMEM) (Gibco, Waltham, MA, USA) containing 10% fetal bovine serum (FBS) (Biological Industries, Israel) in a CO2 incubator at 37°C. The cells were grown to ~60% confluency in six-well plates (NEST Biotechnology Co., Ltd., Wuxi, China) 1 day before transfection. We used Lipofectamine 3000 (Invitrogen, San Diego, CA, USA) to transfect the mall interfering RNA (siRNA) sequences according to the manufacturer’s instructions (23). The cells were cultured with serum-free medium supplemented with 10% FBS 6 h post-transfection and then incubated for 24–72 h. The siRNA sequences of GPI in H1975 and A549 cells were used as follows: siRNA negative control (NC) (5′-TTCTCCGAACGTCACGT-3′), siRNA-GPI #1 (5′-CCAAGATGATACCCTGTGA-3′), siRNA-GPI #2 (5′-GTGCTCA AGTGACCTCTCA-3′), and siRNA-GPI #3 (5′-GGCATATCCTGGTGGATTA-3′).
Cell Proliferative Assay
H1975 and A549 cells transfected with GPI siRNA were placed in a 96-well plate at 3000/well and allowed to adhere overnight. Ten microliters of Cell Counting Kit-8 (CCK-8 solution) was added to each well and incubated for 2 h. A microplate reader (Thermo Fisher Scientific, Shanghai, China) was used to measure the absorbance at 450 nm every day, and the experiment was repeated three times independently (24).
Western Blot
Lung cancer cells were lysed with protease inhibitor, phosphatase inhibitor, and RIPA cleavage buffer (Thermo Fisher Scientific, Shanghai, China). The protein lysates were separated by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) (Servicebio, Wuhan, China) and imprinted on polyvinylidene fluoride (PVDF) membranes (Millipore, Billerica, MA, USA) for analysis (25). Anti-GPI (1:1,000 dilution, #57893; Cell Signaling Technology, Danvers, MA, USA) and anti-GAPDH (1:5,000 dilution, 60004-1-Ig; Proteintech, Wuhan, China) were incubated overnight at 4°C. Horseradish peroxidase (HRP)-conjugated secondary antibodies (1:2,000 dilution; A0181, Beyotime, Shanghai, China) were added for 2 h at room temperature. Western blot analysis was performed with ImageJ software.
Cell Cycle Analysis
A549 and H1975 cells were cultured in six-well serum-free culture with 1 × 105 cells for 12 h and were transfected with GPI siRNA. The cells were collected after 48 h, washed twice with phosphate-buffered saline (PBS), and then fixed with 70% pre-cooled ethanol at 4°C for 12 h. The cells were stained with RNase A-containing PI buffer (KGA511; Keygen, Jiangsu, China) for 30 min. Flow cytometry (Beckman Coulter Quanta SC System) was used to detect the cell cycle distribution.
Immunohistochemistry Staining
This study was approved by the Institutional Research Ethics Committee of Taihe Hospital. Forty-six paraffin-embedded lung adenocarcinoma tissues and paracarcinoma tissues were used for immunohistochemistry (IHC) staining. Each tissue was incubated overnight with a primary antibody for GPI (1:1,000 dilution, #94068; Cell Signaling Technology, Danvers, MA, USA) at 4°C. After washing with PBS, each section was incubated with an HRP-labeled goat anti-rabbit secondary antibody (1:2,000 dilution, A0181; Beyotime, Shanghai, China) for 2 h at room temperature. IHC staining was achieved using 3,3′-diaminobenzidine and counterstained with hematoxylin (Beyotime, Shanghai, China). The IHC staining results were analyzed and scored by two pathologists who were blinded to the sources of the clinical samples. ImageJ software was used to analyze the intensity of staining using a semiquantitative integration method.
Statistical Analysis
Box plot and scatter plot methods were used to detect the expression level of the GPI gene in LUAD patients. The cutoff value of GPI expression was selected as the median method of gene expression. The Wilcoxon signed-rank test and logistic regression were used to analyze the relationship between the clinical features of LUAD and the expression of GPI. Log-rank test was used to test the p-value. Single-factor and multivariate Cox analyses were used to screen potential prognostic factors. In all analyses, *, **, and *** indicate p < 0.05, p < 0.01 and p < 0.001, respectively.
Results
GPI Is Highly Expressed in LUAD
We developed a flow diagram to show our process (Figure 1). To explore the expression level of GPI in normal and tumor tissues, we used the Oncomine database to analyze the expression levels of GPI mRNA in different tumors and normal tissues. The results showed that, compared with that in normal tissues, GPI expression was significantly higher in tumor tissues such as bladder cancer, colorectal cancer, gastric cancer, kidney cancer, lung cancer, ovarian cancer, and lymphoma. In addition, a lower expression was observed in brain tumors, central nervous system (CNS) tumors, leukemia, and sarcoma (Figure 2A). Similar results were also obtained from the TIMER database of TCGA clinical patients (Figure 2B). In addition, to understand the mutation level of GPI in LUAD, we analyzed its genome and copy number. The OncoPrint map of the GPI gene of LUAD patients in TCGA dataset was analyzed using a cBioPortal map (Figure 2C), and the results showed that GPI had less than 6% gene amplification and missense mutations.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-g002.jpg)
Figure 2 The expression level of glucose-6-phosphate isomerase (GPI) in different human cancers. (A) Increased or decreased GPI expression in datasets of different cancers compared with normal tissues in the Oncomine database. (B) TIMER was used to detect the expression levels of GPI in different tumors in The Cancer Genome Atlas (TCGA) database. (C) The cBioPortal OncoPrint map shows the distribution of GPI genome changes in lung adenocarcinoma (LUAD) patients. (D–J) Expression level of GPI in normal tissues and paired adjacent tissues [unmatched tissues (D) and matched tissues (E)] and the tumor tissues from patients with different clinical characteristics in TCGA [clinical stage (F), N stage (G), primary therapy outcome (H), smoking pack-years (I), and OS event (J)]. **p < 0.01, ***p < 0.001.
We used TCGA database to evaluate the mRNA expression levels of GPI in LUAD patients and compared these with the expression levels in adjacent tissues. The results showed that the expression level of GPI in LUAD was significantly higher than that in paracarcinoma tissues (p < 0.01) (Figure 2D). The results were verified in LUAD and paired paracarcinoma tissues (Figure 2E).
To analyze the correlation between GPI expression and clinical characteristics in LUAD patients, we analyzed the mRNA expression levels of GPI in different clinical categories in TCGA database. The association identified between GPI expression and clinical features in patients with LUAD is summarized in Table 1. The results showed that the high expression of GPI was significantly related to the clinical stage (p < 0.05), N stage (p < 0.05), initial treatment result (p < 0.05), and OS event (p < 0.01) of these patients (Figures 2F–J and Table 1). To further elucidate the abnormal upregulated mechanisms of GPI in LUAD tissues, we also explored the correlation of the expression levels of GPI and their methylation status using public databases (DiseaseMeth, version 2.0). The data showed that the difference in the methylation levels between tumor tissues and paracarcinoma tissues was not significant (Supplementary Table S1). These data suggest that GPI is significantly upregulated in LUAD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-t001.jpg)
Table 1 Correlation between glucose-6-phosphate isomerase (GPI) expression and the clinicopathological features of the lung adenocarcinoma (LUAD) cases from The Cancer Genome Atlas (TCGA).
Validation Using Independent External Databases and Clinical Specimens
To further verify the expression level of GPI in LUAD, we selected another four independent external GEO datasets (validation cohort) to analyze the GPI transcription levels of cancer tissues and adjacent tissues in LUAD, including GSE10072, GSE32863, GSE7503, and GSE30219. The results showed that the transcription level of GPI in LUAD was significantly higher than that in noncancerous adjacent tissues from four LUAD datasets (p < 0.001) (Figures 3A–D). The results of this comparison were also verified in LUAD tissues and accurately matched noncancerous tissues (GSE30219 does not contain accurately matched noncancerous tissues) and lung cancer cell lines (Figures 3E–H).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-g003.jpg)
Figure 3 Validation of the high expression of glucose-6-phosphate isomerase (GPI) in lung adenocarcinoma (LUAD) using independent external databases and cell lines. (A–D) Expression of GPI in tumor and unpaired paracarcinoma tissues of the GSE10072, GSE32863, GSE7503, and GSE30219 datasets in the Gene Expression Omnibus (GEO) database. (E–G) Expression level of GPI in tumor and paired adjacent tissues in the GSE10072, GSE32863, and GSE7503 datasets of the GEO database. (H) Transcriptional level of GPI in different non-small-cell lung cancer (NSCLC) cell lines. (I) Validation of the expression level of GPI in LUAD using the Human Protein Atlas database (immunohistochemistry). (J) Western blot detecting the protein expression level of GPI in different NSCLC cell lines. ***p < 0.001.
At the same time, we analyzed the protein expression level of GPI in LUAD based on the HPA database. The data showed that the expression level of GPI in tumors was higher than that in adjacent tissues (Figure 3I). In addition, the expression level of GPI in NSCLC cell lines was detected by Western blotting. The results showed that the expression level of GPI in A549, H1975, H1299, SK-MES-1, and H460 cells was significantly higher than that in the normal lung epithelial cell line BEAS-2B (Figure 3J). Finally, clinical specimens of LUAD were collected to investigate GPI expression, and the results indicated that the expression level of GPI in tumors was higher than that in paracarcinoma (Figures 4A, B). These data verify that GPI is highly expressed in LUAD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-g004.jpg)
Figure 4 Protein expression of glucose-6-phosphate isomerase (GPI) in clinical specimens of lung adenocarcinoma (LUAD). (A) Immunohistochemical staining of GPI was performed in tumor tissues (n = 46) and paracarcinoma tissues (n = 46). Representative images are shown. Score bars, 50 μm. (B) Staining was quantified, as shown. The dot plot depicts the means and standard deviation of 46 images of tumor tissues and adjacent normal tissues. *p < 0.05.
High Expression of GPI Is an Independent Prognostic Factor for the Overall Survival of LUAD
To identify whether GPI expression affects patient survival, we classified LUAD patients in TCGA database into a high GPI expression group (the top 50% of samples with the highest expression) and a low GPI expression group (the remaining 50% of the samples) in order to perform survival analysis according to the mean expression value of GPI. The Kaplan–Meier survival analysis showed that the high expression of GPI was related to the poor prognosis of LUAD patients (HR = 1.44, p = 0.014) (Figure 5A). Subgroup analysis showed that a high GPI expression was significantly correlated with poor prognosis in LUAD in the following cases: patients over 65 years old (HR = 1.94, p = 0.002), female patients (HR = 1.60, p = 0.021), T1 (HR = 2.05, p = 0.25), T3 (HR = 4.22, p = 0.005), pathological stage II (HR = 1.90, p = 0.024), CR/PR (complete/partial response; HR = 1.59, p = 0.035), smokers (HR = 1.48, p = 0.018), and smokers over 40 years old (HR = 2.10, p = 0.006). These data are shown in Figures 5B–I.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-g005.jpg)
Figure 5 Kaplan–Meier survival curve analysis of the prognostic significance of a high and a low expression of glucose-6-phosphate isomerase (GPI) in lung adenocarcinoma (LUAD) using The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. (A) Kaplan–Meier estimates of the overall survival probability of TCGA patients in all LUAD patients. (B–L) Subgroup analysis for age greater than 65 years (B), female patients (C), T1 (D), T3 (E), pathological staging stage II (F), complete/partial response (CR/PR) (G), smokers (H), and smokers over 40 years old (I) and tumor patients in the GEO datasets: GSE31210 (J), GSE37745 (K), and GSE3141 (L).
In addition, we used three other independent external GEO datasets—GSE31210, GSE37745, and GSE3141—as the validation cohort to verify the reproducibility and transferability of GPI expression data in the prognosis of LUAD patients. The Kaplan–Meier survival analysis showed that patients with increased GPI expression in three independent external datasets had shorter OS (Figures 5J–L). Univariate Cox analysis demonstrated that a high GPI expression was significantly correlated with poor OS (HR = 1.435, 95% CI = 1.075–1.916, p = 0.014) (Table 2). These data suggest that a high expression of GPI is an independent prognostic factor for the OS of LUAD patients.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-t002.jpg)
Table 2 Univariate and multivariate Cox regression analyses of the clinical characteristics associated with overall survival in lung adenocarcinoma (LUAD) in The Cancer Genome Atlas (TCGA).
Diagnostic Value of GPI Expression in LUAD
To analyze the diagnostic value of GPI expression in LUAD, we performed ROC curve and nomogram analysis on the GPI gene expression data of TCGA database to evaluate the diagnostic value of the gene. The area under the ROC curve (AUC) was 0.940, suggesting a higher diagnostic value, as shown in Figure 6A. We combined the expression level of GPI with the clinical variables to construct a nomogram in order to predict the survival probability of patients at 1, 3, and 5 years. The nomogram indicated that the prognostic prediction of the expression level of GPI was better than that of the traditional clinical features of age and sex (Figure 6B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-g006.jpg)
Figure 6 Diagnostic value of glucose-6-phosphate isomerase (GPI) expression in lung adenocarcinoma (LUAD). (A) Receiver operating characteristic (ROC) curve analysis for GPI expression in LUAD and adjacent tissue. (B) Nomogram survival prediction chart for predicting the 1-, 3-, and 5-year overall survival rates.
GPI Is Closely Related to LUAD Cell Cycle Regulation
To understand the biological function of GPI in LUAD, we used the LinkFinder module of the LnkedOmics website to detect the co-expression pattern of GPI in LUAD of TCGA database. The red dot indicates the top 25 genes that were positively correlated with GPI, and the green dot represents the bottom 25 genes that were negatively correlated with GPI (Figure 7A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-g007.jpg)
Figure 7 Glucose-6-phosphate isomerase (GPI) functional clustering and interaction network analysis of GPI-related genes. (A) Heatmap showing the top 50 genes in lung adenocarcinoma (LUAD) that were positively and negatively related to GPI. Red represents positively related genes and blue represents negatively related genes. (B) Gene Ontology (GO) term and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of GPI-related genes in LUAD. (C) Venn diagram of the GPI-related genes and the survival-related and upregulated genes in LUAD. (D) GO term and KEGG pathway analyses of the GPI-related genes and LUAD survival-related genes in LUAD.
We used DAVID Functional Annotation Bioinformatics Microarray Analysis to identify the enriched GO functional enrichment and KEGG pathways among the GPI-related genes (top 600) and found that cell cycle signaling pathways, RNA replication, and mismatch repair were enriched among these genes (Figure 7B).
To identify genes with the same regulatory direction in high-GPI and non-survival patients, we crossed the 600 genes with the highest GPI correlation with the 148 survival-related upregulated genes in LUAD and detected 56 genes at the intersection that were related to GPI and LUAD survival (Figure 7C). These 56 protein-coding genes may be potential genetic biomarkers for LUAD patients. GO functional enrichment and KEGG pathway analysis were performed for these 56 genes, and the results showed that differentially expressed genes (DEGs) were significantly enriched in the cell cycle from the GO terms and KEGG pathways (Figure 7D).
After discovering significantly different pathways, PPIs and correlated analysis were used to identify the interactions between these 56 proteins. We found that there was a stronger enrichment network among these proteins than in random proteins, and these genes were particularly related to cell cycle pathways (Figure 8A). Gene co-expression correlation analysis showed that most of the proteins in the network had a strong positive correlation with each other (Figure 8B). Therefore, these GPI-associated established genes, especially cell cycle-related genes, have a strong intertwined interaction and can be used as multigene biomarkers to predict the survival of LUAD patients.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-g008.jpg)
Figure 8 Glucose-6-phosphate isomerase (GPI)-related gene interaction network (A) and gene co-expression matrix (B).
GPI Affects LUAD Cell Proliferation by Regulating the Cell Cycle
To explore the regulatory effect of GPI on the cell cycle, GPI-targeted siRNA sequences were used to interfere with GPI expression, and the knockout efficiency was detected by Western blotting. The results showed that the expression of GPI was significantly decreased after transfection of GPI siRNA (Figure 9A). CCK-8 assays were applied to investigate the effect of GPI on cell proliferation. The results showed that silencing GPI significantly inhibited the proliferation of A549 and H1975 cells (Figures 9B, C). The effect of silencing GPI on the cell cycle of LUAD cells was further determined by flow cytometry. The results showed that A549 and H1975 cells revealed G2/M cell cycle arrest after GPI knockdown (Figure 9D). These data indicated that GPI affects the proliferation of LUAD cells by regulating cell cycle progression.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-g009.jpg)
Figure 9 Glucose-6-phosphate isomerase (GPI) regulates the cycle progression of lung adenocarcinoma (LUAD) cells. (A) Western blot detection of the interfering ability of GPI siRNA in A549 cells. A549 and H1975 cells were transfected with GPI siRNA #SI1 and #SI2. (B–D) CCK-8 was used to detect the proliferation ability of A549 cells (B) and H1975 cells (C), and flow cytometry was used to detect the cell cycle distribution of A549 and H1975 cells (D). *p < 0.05, **p < 0.01.
GPI Is Closely Related to Cell Cycle Regulatory Genes in LUAD
The above bioinformatics analysis indicated that GPI is significantly enriched in the cell cycle pathway. Therefore, we further analyzed the correlation between GPI and the cell cycle regulation genes in TCGA-LUAD, and we found that the expressions of the cell cycle regulatory genes CCNA2, CCNB1, CCNB2, CCNE1, CHEK1, BUB1B, ESPL1, PTTG1, PCNA, PKMYT1, CDC45, PLK1, MCM2, MCM4, MCM6, E2F1, CDC6, CDC20, CDC25A, and CDC25C were positively correlated with GPI (r > 0.3, p < 0.001) (Figures 10A–T). These data verify that GPI is closely related to LUAD cell cycle regulation genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-g010.jpg)
Figure 10 Correlation analysis between glucose-6-phosphate isomerase (GPI) and the cell cycle regulatory genes in lung adenocarcinoma (LUAD) in The Cancer Genome Atlas (TCGA). (A) CCNA2, (B) CCNB1, (C) CCNB2, (D) CCNE1, (E) CHEK1, (F) BUB1B, (G) ESPL1, (H) PTTG1, (I) PCNA, (J) PKMYT1, (K) CDC45, (L) PLK1, (M) MCM2, (N) MCM4, (O) MCM6, (P) E2F1, (Q) CDC6, (R) CDC20, (S) CDC25A, and (T) CDC25C.
Correlation of GPI Expression With Immune Characteristics
Research has shown that the levels of GPI are elevated in autoimmune diseases (26). In this study, we found that GPI is abnormally highly expressed in LUAD, so we speculate that GPI may be involved in regulating the tumor immune response. To explore the correlation between the expression level of GPI and tumor immune response, we used the TIMER database to investigate immune infiltration in LUAD with different GPI expression levels. The results showed that the infiltration levels of CD8+ T cells, dendritic cells (DCs), eosinophils, interdigitating cells, macrophages, mast cells, NK CD56 dim cells, central memory T cells, effective memory T cells, and follicular helper T cells in LUAD patients with high GPI expression were significantly lower than those in patients with low GPI expression. In contrast, the infiltration of Th2 cells, regulatory T cells (Tregs), and gamma delta T cells in LUAD patients with high GPI expression was significantly higher than that in patients with low GPI expression. In addition, there was no significant difference in the infiltration levels of Th1 cells, Th17 cells, B cells, natural killer (NK) cells, and neutrophils between patients with high and low GPI expression (Figures 11A, B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-g011.jpg)
Figure 11 Correlation analysis of glucose-6-phosphate isomerase (GPI) expression and immune infiltration in lung adenocarcinoma (LUAD). (A, B) Differential distribution of immune cells in patients with high GPI expression and low GPI expression. (C–N) Correlation between the expression level of GPI and immune infiltration in lung adenocarcinoma (LUAD): (C) Th2 cells, (D) regulatory T cells, (E) gamma delta T cells, (F) NK CD56bright cells, (G) dendritic cells (DC), (H) immature DC, (I) Mast cells, (J) central memory T cells, (K) Macrophages, (L) Eosinophils, (M) effective memory T cell, (N) CD8+ T cells. *p < 0.05, **p < 0.01, ***p < 0.001, ns, no significance.
We further analyzed the correlation between the expression level of GPI and immune infiltration in LUAD, and the results showed that the expression level of GPI was positively correlated with the infiltrating levels of Th2 cells (r = 0.47, p < 0.001) (Figure 11C), Tregs (r = 0.13, p < 0.002) (Figure 11D), gamma delta T cells (r = 0.110, p =0.008) (Figure 11E) and NK CD56bright cells (r = 0.230, p <0.001) (Figure 11F). In contrast, the expression level of GPI was negatively correlated with the infiltrating levels of DC (r = -0.200, p < 0.001) (Figure 11G), immature DC (r = -0.200, p < 0.001) (Figure 11H), Mast cells (r = -0.260, p < 0.001) (Figure 11I), central memory T cells (r = -0.240, p < 0.001) (Figure 11J), (K) Macrophages (r = -0.140, p = 0.001) (Figure 11K), Eosinophils (r = -0.190, p < 0.001) (Figure 11L), effective memory T cell (r = -0.150, p < 0.001) (Figure 11M), CD8+ T cells (r = -0.0960, p = 0.026) (Figure 11N). A Th1/Th2 balance drift exists in many tumors, such as lung cancer, cervical cancer, breast cancer, gastric cancer, and other cancers (27). Th2 cells are often dominant, which may be related to tumor immune escape (28). Tregs also play an important role in tumor immune escape (29). Therefore, these data suggest that GPI may inhibit the tumor immune response of LUAD by positively regulating Th2 and Treg infiltration into tumors and negatively regulating CD8+ T-cell, central memory T-cell, DC, macrophage, mast cell, and eosinophil infiltration into tumors.
Chemokines and chemokine receptors are essential for the infiltration of immune cells into tumors (30). Therefore, we used TISIDB to analyze the correlation between the expression level of GPI and immune cell chemokines and chemokine receptors in LUAD. The heatmap results showed that several chemokines and chemokine receptors were significantly correlated with the expression of GPI in LUAD (Figures 12A, B). These results demonstrate that the GPI gene may play an important role in tumor immunity. To further clarify the relationship between GPI expression and immune cell migration, we comprehensively analyzed the correlation between GPI expression and chemokines/chemokine receptors. The results showed that GPI expression was negatively correlated with CCL14 (r = −0.373, p < 2.2e−16), CCR4 (r = −0.321, p < 9.17e−14), CCR6 (r = −0.369, p < 2.2e−16), and CX3CR1 (r = −0.376, p < 2.2e−16) (Figures 12C–F), whereas GPI was not highly correlated with other chemokines/chemokine receptors (−0.3 < r < 0.3). These data indicated that GPI is negatively correlated with the expressions of chemokines/chemokine receptors in LUAD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-g012.jpg)
Figure 12 Correlation analysis between glucose-6-phosphate isomerase (GPI) expression and chemokines and/or chemokine receptors. (A) Heatmap analysis of the correlation between GPI and chemokines in tumors. (B) Heatmap analysis of the correlation between GPI and chemokine receptors in tumors. (C–F) GPI expression in lung adenocarcinoma (LUAD) is negatively correlated with CCL14, CCR4, CCR6, and CX3CR1.
Immune checkpoint inhibitors (ICIs) are a significantly novel strategy for tumor immunotherapy that has gradually improved the prognosis of patients with many types of cancers (31). Subsequently, we analyzed the correlation between GPI and the expressions of immunoinhibitors and immunostimulators in different types of human cancers using the TISIDB database (Figures 13A, B). Interestingly, there was no significant correlation between GPI and immunoinhibitor expression. However, GPI was negatively correlated with the expression of several immunostimulators, such as CD40L (r = −0.366, p < 6.85e−19), IL6R (r = −0.365, p < 2.02e−18), and TMEM173 (r = −0.303, p < 2.37e−12), in LUAD (Figures 13C–E). Therefore, these results suggest that GPI may play a role in regulating tumor immunity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/752642/fonc-11-752642-HTML/image_m/fonc-11-752642-g013.jpg)
Figure 13 Correlation analysis between glucose-6-phosphate isomerase (GPI) expression and immunoinhibitors and immunostimulators. (A) Correlation between GPI and immunoinhibitors in tumors by heatmap analysis. (B) Correlation between GPI and immunostimulators in tumors by heatmap analysis. (C–E) GPI expression is negatively correlated with the immunostimulators CD40L, IL6R, and TMEM173 in lung adenocarcinoma (LUAD).
Discussion
GPI encodes a member of the protein family of glucose phosphate isomerases. The encoded protein has been identified as a part-time protein based on its ability to perform mechanically different functions (10). In recent years, the role of GPI in tumors has gradually attracted attention. Studies have found that GPI is related to metastasis and poor prognosis in colorectal cancer, renal cell carcinoma, breast cancer, and endometrial cancer (13, 14, 32). At present, there are few studies on GPI in tumors. Considering the limited research on the GPI gene in cancer, we decided to conduct a fully integrated bioinformatics analysis to identify its biological function and potential regulatory pathway in LUAD. In this study, we first tried to determine the expression and prognostic value of the GPI gene in cancer and found that GPI expression was upregulated in LUAD. A high GPI expression is correlated with poor OS in LUAD. At the same time, a high GPI expression is related to poor clinicopathological characteristics. In addition, GPI has certain reference value for the diagnosis and prognosis of LUAD. These findings strongly suggest that GPI may be used as a biomarker for oncogenes and prognosis. To explain the underlying molecular mechanism of GPI affecting the prognosis of LUAD, we defined 56 genes from the GPI gene network as our hub genes, which can predict the survival rate and pathological stage of patients with LUAD. We also analyzed the expression and prognostic value of the GPI gene in lung squamous cell carcinoma (LUSC) and found that GPI expression was upregulated in LUSC, but not significantly associated with prognosis (data not shown).
Tumors are a type of disease characterized by uncontrolled cell growth and proliferation. Tumor cells are different from normal cells in terms of proliferation, differentiation, and apoptosis (33). Among these cells, cell cycle disorder is one of the most important features (34). In this study, to analyze the molecular mechanism of GPI-mediated LUAD progression and poor prognosis, we performed GO functional enrichment and KEGG pathway analysis and found that the cell cycle gene set was most significantly enriched in the GPI high-expression group. This conclusion was also verified by GSEA. To identify the underlying molecular mechanisms of GPI affecting the prognosis of LUAD, we defined 56 overlapping genes as our hub genes derived from crossing the 600 genes with the highest GPI correlation with the 148 survival-related upregulated genes in LUAD. GO functional enrichment and KEGG pathway analysis were used to analyze the biological processes and signaling pathways of the hub genes. Interestingly, cell cycle BP were the most enriched among the hub genes. The expression levels of cyclin genes in cancer cells are often different from those in normal cells (34). In this study, we found that GPI was positively correlated with the expressions of cell cycle regulation genes, such as CCNA2, CCNB1, CCNB1, and CCNE1. Our cell cycle experiments found that silencing GPI prevented LUAD cells from transitioning from the G2 phase to the M phase. Our data suggest that GPI may regulate the progression of LUAD by regulating the cell cycle transition from the G2 to the M phase.
The development of a cancer is related to its surrounding environment and tumor microenvironment (TME). The TME is composed of immune cells, the extracellular matrix, mesenchymal cells, and inflammatory mediators, which have an impact on tumor growth, metastasis, and clinical survival results (35). Although an effective immune response can have an antitumor effect, cancer cells have evolved a variety of mechanisms, including antigen presentation dysfunction and the recruitment of immunosuppressive cells, which promote tumors to evade the attack of immune cells. Previous studies have reported that immune infiltration can affect the prognosis of patients, and tumor-infiltrating lymphocyte grade is an independent predictor of sentinel lymph node status in patients with tumors (36). This study explored the correlation between the expression of GPI and the level of immune infiltration of LUAD. Our results showed that a high expression of GPI is positively correlated with the infiltration level of Th2 cells and Tregs; however, it is negatively correlated with the infiltration level of CD8+ T cells, DCs, macrophages, and eosinophils. These findings indicate that GPI may play an essential regulatory role in the tumor immune microenvironment and the development of LUAD.
Chemokines and their receptors play an important role in the directional migration of immune cells (37–40). In this study, the TISIDB database was used to analyze the correlation between the expression level of GPI and the expressions of immune cell chemokines and chemokine receptors in LUAD. The results showed that the expression level of GPI was negatively correlated with the expressions of CCL14, CCR4, CCR6, and CX3CR1, suggesting that a high expression of GPI may inhibit the migration of immune cells into the TME. Studies show that the active form of CCL14 binds to CCR1, CCR3, and CCR5 to promote the chemotaxis of monocytes, eosinophils, and T lymphocytes (41, 42). The inflammatory receptor CCR4 promotes the respiratory burst and phagocytosis of macrophages (43) and controls the differentiation of M1/M2 macrophages (44). CCR6 controls integrin-mediated adhesion in B cells (45, 46). CX3CR1 controls the survival of white blood cells and the activation of NK cells (47), which may explain how GPI regulates immune infiltration in LUAD.
In summary, our study found that the expression of GPI is significantly upregulated and is closely related to the poor prognosis of LUAD patients. GPI has certain reference value for the diagnosis and prognosis of LUAD. GPI may affect the progression of LUAD by regulating the cell cycle and immune infiltration. It may serve as a novel prognostic biomarker for LUAD patients.
However, even though we conducted a systematic analysis on GPI and used different databases for cross-validation, there were some limitations to this study. Firstly, the comparability and reproducibility of the microarray data generated in different laboratories remain controversial, which would lead to systemic bias. Secondly, in vivo/in vitro experiments are required to demonstrate the effect of GPI on the tumor cell cycle in order to improve the reliability of our results. Thirdly, even though we concluded that GPI expression is strongly related to immune infiltration and the prognosis of LUAD, we lack direct evidence that GPI influences prognosis by participating in immune infiltration. These issues are worth further exploration in the future.
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 authors.
Ethics Statement
The studies involving human participants were reviewed and approved by the Institutional Research Ethics Committee of Taihe Hospital. The patients/participants provided written informed consent to participate in this study.
Author Contributions
ZL and YY conceived and designed the work. JH and XD acquired, analyzed, and interpreted the data. RS, MLu, and MLi drafted and revised the article. BG, ZP, and YL helped with the cell culture and transfection. TZ and CT performed the cell cycle analysis and Western blot experiments. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Nature Science Foundation of China (nos. 81802997 and 81602391), the Foundation for Free Exploration of Hubei University of Medicine (no. FDFR201802), the Natural Science Foundation of Hubei Province of China (nos. 2019CFA034, 2017CFB167, 2018CFB405, and 2017CFB456), the Natural Science Foundation of Hubei Provincial Department of Education (nos. Q20202106, D20172102, and Q20162109), and Hubei Provincial Health Commission (no. WJ2021M048).
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/fonc.2021.752642/full#supplementary-material
References
1. Siegel F, Miller J, Jemal I, Siegel RL, Torre LA, Jemal A. Erratum: Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2020) 70(4):313. doi: 10.3322/caac.21609
2. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2020. CA Cancer J Clin (2020) 70(1):7–30. doi: 10.3322/caac.21590
3. AbdulJabbar K, Raza SEA, Rosenthal R, Jamal-Hanjani M, Veeriah S, Akarca A, et al. Geospatial Immune Variability Illuminates Differential Evolution of Lung Adenocarcinoma. Nat Med (2020) 26(7):1054–62. doi: 10.1038/s41591-020-0900-x
4. Herbst RS, Morgensztern D, Boshoff C. The Biology and Management of Non-Small Cell Lung Cancer. Nature (2018) 553(7689):446–54. doi: 10.1038/nature25183
5. Thomas A, Hassan R. Immunotherapies for Non-Small-Cell Lung Cancer and Mesothelioma. Lancet Oncol (2012) 13(7):e301–10. doi: 10.1016/S1470-2045(12)70126-2
6. Caushi JX, Zhang J, Ji Z, Vaghasia A, Zhang B, Hsiue EH, et al. Transcriptional Programs of Neoantigen-Specific TIL in Anti-PD-1-Treated Lung Cancers. Nature (2021) 596(7870):126–32. doi: 10.1038/s41586-021-03752-4
7. Sinicrope FA, Graham RP. Tumor-Infiltrating Lymphocytes for Prognostic Stratification in Nonmetastatic Colon Cancer-Are We There Yet? JAMA Oncol (2021) 7(7):969–70. doi: 10.1001/jamaoncol.2021.0177
8. Chapman PB. Targeting Tumor-Rejection Antigens in Melanoma With Tumor-Infiltrating Lymphocytes. J Clin Oncol (2021) 39(24):2640–2. doi: 10.1200/JCO.21.01012
9. Knight AL, Yan X, Hamamichi S, Ajjuri RR, Mazzulli JR, Zhang MW, et al. The Glycolytic Enzyme, GPI, Is a Functionally Conserved Modifier of Dopaminergic Neurodegeneration in Parkinson’s Models. Cell Metab (2014) 20(1):145–57. doi: 10.1016/j.cmet.2014.04.017
10. Kroemer G, Lopez-Otin C, Madeo F, de Cabo R. Carbotoxicity-Noxious Effects of Carbohydrates. Cell (2018) 175(3):605–14. doi: 10.1016/j.cell.2018.07.044
11. Kassahn D, Kolb C, Solomon S, Bochtler P, Illges H. Few Human Autoimmune Sera Detect GPI. Nat Immunol (2002) 3(5):411–2. doi: 10.1038/ni0502-411b
12. Lucarelli G, Rutigliano M, Sanguedolce F, Galleggiante V, Giglio A, Cagiano S, et al. Increased Expression of the Autocrine Motility Factor Is Associated With Poor Prognosis in Patients With Clear Cell-Renal Cell Carcinoma. Medicine (Baltimore) (2015) 94(46):e2117. doi: 10.1097/MD.0000000000002117
13. van Veen M, Matas-Rico E, van de Wetering K, Leyton-Puig D, Kedziora KM, De Lorenzi V, et al. Negative Regulation of Urokinase Receptor Activity by a GPI-Specific Phospholipase C in Breast Cancer Cells. Elife (2017) 6:e23649. doi: 10.7554/eLife.23649
14. Wu ST, Liu B, Ai ZZ, Hong ZC, You PT, Wu HZ, et al. Esculetin Inhibits Cancer Cell Glycolysis by Binding Tumor PGK2, GPD2, and GPI. Front Pharmacol (2020) 11:379. doi: 10.3389/fphar.2020.00379
15. Gallardo-Pérez JC, Adán-Ladrón de Guevara A, Marín-Hernández A, Moreno-Sánchez R, Rodríguez-Enríquez S. HPI/AMF Inhibition Halts the Development of the Aggressive Phenotype of Breast Cancer Stem Cells. Biochim Biophys Acta Mol Cell Res (2017) 1864(10):1679–90. doi: 10.1016/j.bbamcr.2017.06.015
16. Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, et al. Oncomine 3.0: Genes, Pathways, and Networks in a Collection of 18,000 Cancer Gene Expression Profiles. Neoplasia (New York NY) (2007) 9(2):166–80. doi: 10.1593/neo.07112
17. Lánczky A, Nagy Á, Bottai G, Munkácsy G, Szabó A, Santarpia L, et al. Mirpower: A Web-Tool to Validate Survival-Associated miRNAs Utilizing Expression Data From 2178 Breast Cancer Patients. Breast Cancer Res Treat (2016) 160(3):439–46. doi: 10.1007/s10549-016-4013-7
18. Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: New Perspectives on Genomes, Pathways, Diseases and Drugs. Nucleic Acids Res (2017) 45(D1):D353–61. doi: 10.1093/nar/gkw1092
19. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: A Web Server for Cancer and Normal Gene Expression Profiling and Interactive Analyses. Nucleic Acids Res (2017) 45(W1):W98–102. doi: 10.1093/nar/gkx247
20. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene Set Enrichment Analysis: A Knowledge-Based Approach for Interpreting Genome-Wide Expression Profiles. Proc Natl Acad Sci USA (2005) 102(43):15545–50. doi: 10.1073/pnas.0506580102
21. Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, et al. PGC-1alpha-Responsive Genes Involved in Oxidative Phosphorylation Are Coordinately Downregulated in Human Diabetes. Nat Genet (2003) 34(3):267–73. doi: 10.1038/ng1180
22. Yu G, Wang LG, Han Y, He QY. Clusterprofiler: An R Package for Comparing Biological Themes Among Gene Clusters. Omics J Integr Biol (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118
23. Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: Analyzing Multi-Omics Data Within and Across 32 Cancer Types. Nucleic Acids Res (2018) 46(D1):D956–D63. doi: 10.1093/nar/gkx1090
24. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res (2017) 77(21):e108–e10. doi: 10.1158/0008-5472.CAN-17-0307
25. Ru B, Wong CN, Tong Y, Zhong JY, Zhong SSW, Wu WC, et al. TISIDB: An Integrated Repository Portal for Tumor-Immune System Interactions. Bioinformatics (Oxford England) (2019) 35(20):4200–2. doi: 10.1093/bioinformatics/btz210
26. Deng XZ, Geng SS, Luo M, Chai JJ, Xu Y, Chen CL, et al. Curcumin Potentiates Laryngeal Squamous Carcinoma Radiosensitivity via NF-KappaB Inhibition by Suppressing IKKgamma Expression. J Recept Signal Transduct Res (2020) 40(6):541–9. doi: 10.1080/10799893.2020.1767649
27. Deng X, Chen C, Wu F, Qiu L, Ke Q, Sun R, et al. Curcumin Inhibits the Migration and Invasion of Non-Small-Cell Lung Cancer Cells Through Radiation-Induced Suppression of Epithelial-Mesenchymal Transition and Soluble E-Cadherin Expression. Technol Cancer Res Treat (2020) 19:1533033820947485. doi: 10.1177/1533033820947485
28. Peng S, Wang J, Lu C, Xu Z, Chai JJ, Ke Q, et al. Emodin Enhances Cisplatin Sensitivity in Non-Small Cell Lung Cancer Through Pgp Downregulation. Oncol Lett (2021) 21(3):230. doi: 10.3892/ol.2021.12491
29. Schaller M, Stohl W, Tan SM, Benoit VM, Hilbert DM, Ditzel HJ. Raised Levels of Anti-Glucose-6-Phosphate Isomerase IgG in Serum and Synovial Fluid From Patients With Inflammatory Arthritis. Ann Rheum Dis (2005) 64(5):743–9. doi: 10.1136/ard.2004.025502
30. Li J, Jie HB, Lei Y, Gildener-Leapman N, Trivedi S, Green T, et al. PD-1/SHP-2 Inhibits Tc1/Th1 Phenotypic Responses and the Activation of T Cells in the Tumor Microenvironment. Cancer Res (2015) 75(3):508–18. doi: 10.1158/0008-5472.CAN-14-1215
31. Zhao X, Liu J, Ge S, Chen C, Li S, Wu X, et al. Saikosaponin A Inhibits Breast Cancer by Regulating Th1/Th2 Balance. Front Pharmacol (2019) 10:624. doi: 10.3389/fphar.2019.00624
32. Luo CT, Liao W, Dadi S, Toure A, Li MO. Graded Foxo1 Activity in Treg Cells Differentiates Tumour Immunity From Spontaneous Autoimmunity. Nature (2016) 529(7587):532–6. doi: 10.1038/nature16486
33. Dyer DP, Medina-Ruiz L, Bartolini R, Schuette F, Hughes CE, Pallas K, et al. Chemokine Receptor Redundancy and Specificity Are Context Dependent. Immunity (2019) 50(2):378–89.e5. doi: 10.1016/j.immuni.2019.01.009
34. Sharma P, Hu-Lieskovan S, Wargo JA, Ribas A. Primary, Adaptive, and Acquired Resistance to Cancer Immunotherapy. Cell (2017) 168(4):707–23. doi: 10.1016/j.cell.2017.01.017
35. Wang M, Chang M, Li C, Chen Q, Hou Z, Xing B, et al. Tumor Microenvironment-activated Reactive Oxygen Species Amplifier for Enzymatic Cascade Cancer Starvation/Chemodynamic /Immunotherapy. Adv Mater (2021) e2106010. doi: 10.1002/adma.202106010
36. Slack FJ, Chinnaiyan AM. The Role of Non-Coding RNAs in Oncology. Cell (2019) 179(5):1033–55. doi: 10.1016/j.cell.2019.10.017
37. Phua SC, Chiba S, Suzuki M, Su E, Roberson EC, Pusapati GV, et al. Dynamic Remodeling of Membrane Composition Drives Cell Cycle Through Primary Cilia Excision. Cell (2019) 178(1):261. doi: 10.1016/j.cell.2019.06.015
38. Ley K. Arrest Chemokines. Microcirculation (New York NY 1994) (2003) 10(3-4):289–95. doi: 10.1038/sj.mn.7800194
39. Kanemitsu N, Ebisuno Y, Tanaka T, Otani K, Hayasaka H, Kaisho T, et al. CXCL13 Is an Arrest Chemokine for B Cells in High Endothelial Venules. Blood (2005) 106(8):2613–8. doi: 10.1182/blood-2005-01-0133
40. Pallandre JR, Krzewski K, Bedel R, Ryffel B, Caignard A, Rohrlich PS, et al. Dendritic Cell and Natural Killer Cell Cross-Talk: A Pivotal Role of CX3CL1 in NK Cytoskeleton Organization and Activation. Blood (2008) 112(12):4420–4. doi: 10.1182/blood-2007-12-126888
41. Lee SS, Cheah YK. The Interplay Between MicroRNAs and Cellular Components of Tumour Microenvironment (TME) on Non-Small-Cell Lung Cancer (NSCLC) Progression. J Immunol Res (2019) 2019:3046379. doi: 10.1155/2019/3046379
42. Azimi F, Scolyer RA, Rumcheva P, Moncrieff M, Murali R, McCarthy SW, et al. Tumor-Infiltrating Lymphocyte Grade Is an Independent Predictor of Sentinel Lymph Node Status and Survival in Patients With Cutaneous Melanoma. J Clin Oncol (2012) 30(21):2678–83. doi: 10.1200/jco.2011.37.8539
43. Peng D, Kryczek I, Nagarsheth N, Zhao L, Wei S, Wang W, et al. Epigenetic Silencing of TH1-Type Chemokines Shapes Tumour Immunity and Immunotherapy. Nature (2015) 527(7577):249–53. doi: 10.1038/nature15520
44. Nagarsheth N, Wicha MS, Zou W. Chemokines in the Cancer Microenvironment and Their Relevance in Cancer Immunotherapy. Nat Rev Immunol (2017) 17(9):559–72. doi: 10.1038/nri.2017.49
45. Zlotnik A, Yoshie O. The Chemokine Superfamily Revisited. Immunity (2012) 36(5):705–16. doi: 10.1016/j.immuni.2012.05.008
46. Matsukawa A, Hogaboam CM, Lukacs NW, Lincoln PM, Evanoff HL, Kunkel SL. Pivotal Role of the CC Chemokine, Macrophage-Derived Chemokine, in the Innate Immune Response. J Immunol (Baltimore Md 1950) (2000) 164(10):5362–8. doi: 10.4049/jimmunol.164.10.5362
Keywords: glucose-6-phosphate isomerase (GPI), prognosis (carcinoma), immune infiltrate, cell cycle, lung adenocarcinoma
Citation: Han J, Deng X, Sun R, Luo M, Liang M, Gu B, Zhang T, Peng Z, Lu Y, Tian C, Yan Y and Luo Z (2021) GPI Is a Prognostic Biomarker and Correlates With Immune Infiltrates in Lung Adenocarcinoma. Front. Oncol. 11:752642. doi: 10.3389/fonc.2021.752642
Received: 03 August 2021; Accepted: 03 November 2021;
Published: 29 November 2021.
Edited by:
Jawed Akhtar Siddiqui, University of Nebraska Medical Center, United StatesReviewed by:
Musarat Ishaq, University of Melbourne, AustraliaYoulang Zhou, Affiliated Hospital of Nantong University, China
Copyright © 2021 Han, Deng, Sun, Luo, Liang, Gu, Zhang, Peng, Lu, Tian, Yan and Luo. 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: Zhiguo Luo, Luozhiguo@hbmu.edu.cn; Yutao Yan, yanyutao_111@163.com
†These authors have contributed equally to this work