- Department of Thoracic Oncology, Tianjin Medical University Cancer Institute and Hospital, National Clinical Research Center for Cancer, Tianjin Key Laboratory of Cancer Prevention and Therapy, Tianjin’s Clinical Research Center for Cancer, Tianjin, China
Background: Glucosamine 6-phosphate N-acetyltransferase (GNPNAT1) is a key enzyme in the hexosamine biosynthetic pathway (HBP), which functions as promoting proliferation in some tumors, yet its potential biological function and mechanism in lung adenocarcinoma (LUAD) have not been explored.
Methods: The mRNA differential expression of GNPNAT1 in LUAD and normal tissues was analyzed using the Cancer Genome Atlas (TCGA) database and validated by real-time PCR. The clinical value of GNPNAT1 in LUAD was investigated based on the data from the TCGA database. Then, immunohistochemistry (IHC) of GNPNAT1 was applied to verify the expression and clinical significance in LUAD from the protein level. The relationship between GNPNAT1 and epigenetics was explored using the cBioPortal database, and the miRNAs regulating GNPNAT1 were found using the miRNA database. The association between GNPNAT1 expression and tumor-infiltrating immune cells in LUAD was observed through the Tumor IMmune Estimation Resource (TIMER). Finally, Gene set enrichment analysis (GSEA) was used to explore the biological signaling pathways involved in GNPNAT1 in LUAD.
Results: GNPNAT1 was upregulated in LUAD compared with normal tissues, which was verified through qRT-PCR in different cell lines (P < 0.05), and associated with patients’ clinical stage, tumor size, and lymphatic metastasis status (all P < 0.01). Kaplan–Meier (KM) analysis suggested that patients with upregulated GNPNAT1 had a relatively poor prognosis (P < 0.0001). Furthermore, multivariate Cox regression analysis indicated that GNPNAT1 was an independent prognostic factor for LUAD (OS, TCGA dataset: HR = 1.028, 95% CI: 1.013–1.044, P < 0.001; OS, validation set: HR = 1.313, 95% CI: 1.130–1.526, P < 0.001). GNPNAT1 overexpression was correlated with DNA copy amplification (P < 0.0001), low DNA methylation (R = −0.52, P < 0.0001), and downregulation of hsa-miR-30d-3p (R = −0.17, P < 0.001). GNPNAT1 expression was linked to B cells (R = −0.304, P < 0.0001), CD4+T cells (R = −0.218, P < 0.0001), and dendritic cells (R = −0.137, P = 0.002). Eventually, GSEA showed that the signaling pathways of the cell cycle, ubiquitin-mediated proteolysis, mismatch repair and p53 were enriched in the GNPNAT1 overexpression group.
Conclusion: GNPNAT1 may be a potential prognostic biomarker and novel target for intervention in LUAD.
Introduction
Lung cancer is the leading cause of cancer-related death around the world. In recent years, adenocarcinoma has become the predominant pathological type in lung cancer, accounting for about 40% of all lung cancer patients (Molina et al., 2008). Although the surgery and conventional chemoradiotherapy have improved patient survival, the overall 5-year survival rate was around 15%. Nowadays, the rapid development of molecular detection technology has made cancer treatment more precise and the clinical application of targeted therapy significantly improved the survival rate of patients with positive driver genes. However, there are still many undetected genetic changes that may be functionally important in lung cancer. Therefore, it is urgent to find potential genetic changes to prompt drug development and improve patients’ prognosis and survival.
There is obvious heterogeneity of different cancer patients in the incentive and etiology, but metabolic abnormalities, especially glycometabolism disorder, occur in nearly all tumors (Hanahan and Weinberg, 2011). GNPNAT1 is a crucial enzyme in hexosamine biosynthetic pathway (HBP), one of the essential glucose metabolism pathways branching off from glycolysis. Active metabolism was often occurred in rapidly proliferating tumor cells to promote the progression and metastasis in cancers. Under this circumstances, cancer cells primarily choose to raise the product conversion rate of HBP by increasing the glucose, glutamine, and other nutrients intake, and change cancer-related signaling pathways, such as the Ras (Ying et al., 2012), mechanistic target of rapamycin 2 (mTORC2) (Moloughney et al., 2018), and transforming growth factor-beta 1 (TGF-β1) (Lucena et al., 2016). Studies found that embryonic cells lacking GNPNAT1 exhibited defects in proliferation and adhesiveness, but increased apoptosis ability (Boehmelt et al., 2000). Furthermore, researchers have shown that GNPNAT1 was upregulated in breast cancer and prostate cancer, which were related to tumor proliferation and metastasis (Kaushik et al., 2016; Chokchaitaweesuk et al., 2019). However, its expression and clinical value in lung adenocarcinoma (LUAD) are still unclear.
In this study, we analyzed the differential expression of GNPNAT1 in LUAD and normal tissues and investigated its clinical implication based on the Cancer Genome Atlas (TCGA) database, which was confirmed by immunohistochemistry (IHC) staining. Then, we preliminarily explored the potential mechanism of GNPNAT1 in the development and progression in LUAD. Our research provided a new perspective on the diagnosis and even treatment of LUAD in the future.
Materials and Methods
Samples and Cell Lines
The training set samples were obtained from the TCGA database, including 535 LUAD and 59 adjacent normal tissue samples. The transcriptome expression data and clinical information (type: HTSeq–FPKM; time: May 20, 2020) were both downloaded from the TCGA GDC1. In the validation set, 116 samples of LUAD and 18 adjacent non-tumor samples were acquired from LUAD patients who underwent surgery between December 2012 and February 2014 at the Tianjin Cancer Institute and Hospital, which was approved by the Ethics Committee of the Tianjin Cancer Institute and Hospital and was consistent with the ethical guidelines of the Helsinki Declaration. All patients signed informed consent, none of who underwent chemotherapy or radiotherapy before surgery. In addition, GSE19188 (Hou et al., 2010), GSE19804 (Lu et al., 2010), GSE31210 (Okayama et al., 2012), and GSE32863 (Selamat et al., 2012) from the GEO database2 were used to verify the differential expression at the mRNA level. The baseline information of LUAD patients from the TCGA database and our hospital are shown in Tables 1, 2. The expression levels of GNPNAT1 in different tumor cells and tissues were acquired from the Cancer Cell Line Encyclopedia (CCLE)3 (Cerami et al., 2012), TCGA, and Genotype-Tissue Expression (GTEx) databases. Normal lung cell lines (BEAS-2B) and LUAD cell lines (including NCI-H1975, NCI-H358, PC-9, HCC827, and NCI-H1299) in this study were purchased from American Type Culture Collection (ATCC) cell bank.
Quantitative Real-Time PCR
Trizol reagent was used to extract the total RNA according to the product protocol (Invitrogen, No.15596026). Reverse transcription reaction was carried out to acquire cDNA to prepare for the quantitative real-time PCR with the PrimeScriptTM RT Master Mix (TaKaRa, RR036A). qPCR was cycled with the CFX96 quantitative real time gene amplification instrument (Bio-Rad) using 2X SG Fast qPCR Master Mix (Sangon Biotech, No. B639271). Primers for GNPNAT1 are followed from 5′ to 3′: AGGGCCTCTACGGTTCCTGT (F), GTGTTGGGGAAATGGCTGGA (R). HS-ACTB was used as reference primers (Sangon Biotech, B662102-0001). The amplification efficiency was assessed by the standard curve. The experiment was repeated three times.
Immunohistochemistry Staining in the Validation Set
Fresh LUAD and normal lung specimens after surgery were fixed with formalin solution and embedded in paraffin blocks. Slides were baked to dry in an oven at 70°C overnight. Then, deparaffinized in three changes of xylene, and rehydrated in sequential incubation with 100, 95, 85, and 75% ethanol. Antigen retrieval was performed with pressure cooking in citrate buffer (pH = 6.0) and allowed to cool at room temperature. Quench endogenous peroxidase activity through incubating slides in 3% H2O2 for 20 min. Subsequently, slides were incubated with the rabbit anti-human GNPNAT1 antibody (Proteintech Cat# 16282-1-AP, RRID: AB_2110243, 1:200) at 4°C overnight. Employing envision secondary antibody for 30 min at room temporary. Slides were stained with 3,3′-diaminobenzidine (DAB) and watched under the microscope. Next, counterstain sections in mayor hematoxylin and rinse slides gently in a distilled water. Transfer slides to a solution containing 1% HCl and 99% ethanol for 6–10 s, and then rinse by water immediately. Eventually, dehydrate by graded ethanol series, incubated in xylene and mount the slides with neutral balsam. GNPNAT1 was localized primarily in the cytoplasm of LUAD cells.
The GNPNAT1 staining intensity was scored across four grades (0, no immune response; 1, weak immune response; 2, moderate immune response; and 3, strong immune response). The percentage of tumor cells positive was evaluated across a range from one to four (1, <25% positive cells; 2, 25–50% positive cells; 3, 50–75% positive cells; and 4, 75–100% positive cells). Each individual intensity score was multiplied by the percentage score to obtain the final IHC score, ranging from 0 to 12.
Effects of DNA Copy Number Alternations and Methylation Status on GNPNAT1
To observe the potential role of GNPNAT1 in LUAD from the epigenetic perspective, we analyzed the different DNA copy number alterations (CNAs) and methylation of GNPNAT1 expression in LUAD by the cBioPortal database4 (Cerami et al., 2012). Moreover, previous studies have shown that abnormal methylation of CpG island was closely related to the occurrence of tumors (Costello et al., 2000; Ehrlich, 2002), so we used MEXPRESS5 (Koch et al., 2019), which can visualize DNA methylation, gene expression and clinical data of TCGA, to study the correlation between methylation status of CpG island and its expression. The correlation coefficient |R| > 0.2 and P < 0.00001 were considered as having a significant difference in this analysis.
The Prediction for Upstream miRNAs of Regulating GNPNAT1
DIANA- microT6 (Paraskevopoulou et al., 2013) and miRWalk7 (Sticht et al., 2018) databases were applied to find out the potential miRNAs regulating GNPNAT1. Then, differential expression miRNAs in LUAD were screened through the OncomiR database8 (Wong et al., 2018). The common miRNAs that must be downregulated in LUAD from the TCGA database in the above databases were seen as potential regulatory miRNAs of GNPNAT1 in LUAD. The TargetScanHuman 7.2 database9 (Agarwal et al., 2015) was used to predict the possible regions combined with GNPNAT1 3′ UTR binding site.
The Correlation Between GNPNAT1 and Tumor-Infiltrating Immune Cells in LUAD
To explore the potential immunomodulatory effect of GNPNAT1 in LUAD, the Tumor IMmune Estimation Resource(TIMER) database10 (Li et al., 2017) was performed to probe the relevance of GNPNAT1 expression with immune cells (B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells). As we all know, tumor cells overexpressed MKI67 and its protein ki67 was the most commonly used as an indicator of tumor proliferation in pathology. Importantly, one recent study has suggested that ki67 overexpression was associated with decreased immune cells and tumor invasion (Mitchell et al., 2019). So, the MKI67 expression in LUAD and normal tissues was observed from the TCGA database and the correlation between GNPNAT1 and MKI67 was evaluated to investigate the role of GNPNAT1 in the immune microenvironment and tumor proliferation.
Protein-Protein Interaction Network and Gene Set Enrichment Analysis
Gene co-expression analysis has already been proven to be useful in exploring specific genes’ functions. Studies have noticed that co-expressed genes participated in similar biological processes in function (Stuart et al., 2003; Pan et al., 2019). Hence, the co-expressed genes of GNPNAT1 were obtained from the Tissue and Cancer Specific Biological Networks (TCSBN)11 (Lee et al., 2018). Then, the acquired genes were imported into Cytoscape12 (Shannon et al., 2003) to construct the protein-protein interaction (PPI) network. Gene set enrichment analysis (GSEA) software (version 4.0.3) (Subramanian et al., 2005) was applied to get the relevant signaling pathways involved in GNPNAT1 enrichment genes in LUAD. The annotated gene set of c2.cp.kegg.v6.2.symbols.gmt was chosen to conduct the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The screening criteria were P < 0.05 and false discovery rate (FDR) q-value < 0.25, which were deemed as significantly enriched.
Statistical Analyses
The correlation between GNPNAT1 and LUAD clinical features was evaluated by non-parametric tests. We performed the Mann–Whitney U test for two groups and the Kruskal–Wallis test for no less than three groups, respectively. In the TCGA database, LUAD patients were divided into expression-high and expression-low groups by the median expression of GNPNAT1 (median value = 9.3). In the validation set, the IHC score was seen as GNPNAT1 expression at the protein level. The score no more than seven was regarded as low expression and more than seven scores were the high expression. Kaplan–Meier (KM) analysis was conducted based on the survival time in the low and high group of GNPNAT1 expression. Moreover, univariate and multivariate Cox regression analyses were used to test the correlation between survival time and clinical prognostic indicators and GAPNAT1 expression. Nomograms were constructed based on the independent factors of Cox multivariate analyses in the TCGA set. The concordance index (C-index) and calibration were also assessed to effectively measure the performance of constructed nomograms. The correlation between GNPNAT1 expression and DNA methylation level or immune cell infiltration levels was analyzed by Pearson correlation. R version 3.6.0, GraphPad Prism 7 and SPSS 26.0 were employed for statistical analyses. P < 0.05 was considered as statistically significant and all statistical analyses applied bilateral detection (“∗”: P < 0.05, “∗∗”: P < 0.01, “∗∗∗”: P < 0.001, and “****”: P < 0.0001, ns: insignificant).
Results
The Upregulation of GNPNAT1 in LUAD
At the cell level, the CCLE database revealed that GNPNAT1 was up-regulated in NSCLC and other 39 different tumor cells (Figure 1A). At the tissue level, TCGA combined with GTEx database analysis indicated that GNPNAT1 was highly expressed in LUAD and other 27 tumors (P < 0.05) (Figure 1B). At the mRNA level, we analyzed 515 patients, including 59 normal tissues and 535 LUAD samples from TCGA. The results unveiled that GNPNAT1 was overexpressed in LUAD (P < 0.0001, Figure 2A), which was further certified by GSE19188 (P = 5.4e-10), GSE19804 (P = 5.4e-9), GSE31210 (P = 3.5e-8), and GSE32683 (P = 2.1e-10) (Figures 2B–E). To verify this finding, we performed the quantitative real-time PCR assay to compare the relative mRNA level between the normal lung epithelial cell (BEAS-2B) and various LUAD cell lines, including H1975, H358, PC-9, HCC827, and H1299. The results demonstrated that the mRNA level of GNPNAT1 was much higher in LUAD cells than that of BEAS-2B (H1975 vs. EBAS-2B, P = 0.0009; H358 vs. EBAS-2B, P = 0.0083; PC-9 vs. EBAS-2B, P = 0.0009; HCC827 vs. BEAS-2B, P = 0.2366; H1299 vs. EBAS-2B, P = 0.0267) (Figure 2F). At the protein level, the IHC results of 116 LUAD and 18 normal lung samples showed that GNPNAT1 had an elevated expression level compared with normal lung tissues (P = 5.2e-12) (Figures 2G,H). The positive control of GNPNAT1 for IHC is human liver tissue (Supplementary Figure 1A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/605754/fmolb-08-605754-HTML/image_m/fmolb-08-605754-g001.jpg)
Figure 1. The expression of GNPNAT1 in cancers. (A) The mRNA expression level of GNPNAT1 in cancer cells. (B) The expression level of GNPNAT1 in different tumor and normal tissues.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/605754/fmolb-08-605754-HTML/image_m/fmolb-08-605754-g002.jpg)
Figure 2. The differential expression of GNPNAT1 in normal tissues and LUAD. (A) The differential expression of GNPNAT1 from the TCGA database. The differential expression of GNPNAT1 in GSE19188 (B), GSE19804 (C), GSE31210 (D), and GSE32863 (E). The expression level in different LUAD and normal lung cell lines (F). The differential expression of GNPNAT1 from clinical samples (G). The typical IHC staining results of GNPNAT1 in normal lung tissue and LUAD tissues (H). *P < 0.05, **P < 0.01, ***P < 0.001, and ****P < 0.0001, ns: insignificant.
The GNPNAT1 Overexpression Predicts Poor Prognosis in LUAD
Survival analysis was conducted based on survival time (OS, overall survival; PFS, progression-free survival) and GNPNAT1 expression. The Kaplan–Meier analysis implied that the high expression group of GNPNAT1 had a shortened OS (P < 0.0001; Figure 3A), but there was no significant difference on PFS between low and high groups (P = 0.084; Supplementary Figure 1B). Correlation analysis displayed that GNPNAT1 was clearly correlated with gender (P = 0.00042), clinical stage (stage II vs. stage I, P = 0.0067; stage III vs. stage I, P = 1.8e-05; stage IV vs. stage I, P > 0.05), tumor size (T2 vs. T1, P = 3.3e-05; T3 vs. T1, P = 0.018; T4 vs. T1, P = 0.0064), lymph node metastasis (P = 2.9e-5) and patients’ status at the end of follow-up (P = 0.023) (Figures 3B–F). There were no significant differences in race (P = 0.075), smoking history (P = 0.190), age (P = 0.140), the Eastern Cooperative Oncology Group (ECOG) score (P = 0.600) and distant metastasis (P = 0.640) (Supplementary Figures 1C–G). Multivariate Cox regression analysis proved that GNPNAT1 expression level and TNM stage were both independent prognostic factors in LUAD (GNPNAT1: HR = 1.028, 95% CI: 1.013–1.044, P = 0.00021; TNM: HR = 1.542, 95% CI: 1.041–2.284, P = 0.0307; Table 3). These results were verified again through the IHC staining in clinical variables, including clinical stage (stage III vs. Early stage, P = 2.0e-6), tumor size (No-T1 vs. T1: P = 0.006) and lymph node metastasis (positive vs. negative, P = 3.4e-8) (Figures 4A–C). Survival analysis and Cox regression analysis further validated the outcomes found from the TCGA database (Figures 4D,E, Table 4, and Supplementary Table 1). To further validate the prognostic value of GNPNAT1 in LUAD, we incorporated GNPNAT1 expression and/or TNM staging to construct predictive nomograms. The calibration plots noted that the nomogram, including the expression of GNPNAT1 (C-index = 0.713), was more precise than another one (C-index = 0.677) (Figure 5 and Supplementary Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/605754/fmolb-08-605754-HTML/image_m/fmolb-08-605754-g003.jpg)
Figure 3. Association of GNPNAT1 expression with clinicopathologic characteristics. (A) The survival curve between low and high expression of GNPNAT1 in LUAD. The correlation between GNPNAT1 and clinicopathologic variables involved in gender (B), cancer status (C), clinical stage (D), lymph node metastasis status (E), and tumor size (F). *P < 0.05, **P < 0.01, ***P < 0.001, and ****P < 0.0001, ns: insignificant.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/605754/fmolb-08-605754-HTML/image_m/fmolb-08-605754-g004.jpg)
Figure 4. Association of GNPNAT1 expression with clinical features in validation set, including TNM stage (A), tumor size (B), and lymph node metastasis status (C). The Kaplan–Meier curves about the correlation between GNPNAT1 expression and overall survival (D) and disease-free survival (E). **P < 0.01 and ****P < 0.0001, ns: insignificant).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/605754/fmolb-08-605754-HTML/image_m/fmolb-08-605754-g005.jpg)
Figure 5. The prognostic nomograms based on GNPNAT1 and TNM stage (A) and TNM stage alone (C). The calibration curve of (A) nomogram for predicting OS at 2 year (B). The calibration curve of (C) nomogram for predicting OS at 2 year (D).
DNA Copy Number Amplification and Hypomethylation Led to GNPNAT1 Upregulation in LUAD
Based on the LUAD data from the cBioPortal database containing mRNA, CNAs, and DNA methylation (n = 450), we discovered that GNPNAT1 amplification was significantly associated with its mRNA overexpression (Total: P = 7.7e-10; Amplification vs. Diploid, P = 1.9e-05; Amplification vs. Gain, P = 0.002; Gain vs. Shallow Deletion, P = 1.5e-06; Amplification vs. Shallow Deletion, P = 6.2e-06) (Figure 6A). Besides, GNPNAT1 expression also had a closely inverse correlation with DNA methylation (Pearson’s R = −0.52, P < 0.01) (Figure 6B). We also noticed a negative relationship between GNPNAT1 expression and 10 CpG islands through MEXPRESS (Supplementary Table 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/605754/fmolb-08-605754-HTML/image_m/fmolb-08-605754-g006.jpg)
Figure 6. The correlation between GNPNAT1 and some epigenetic factors. The association of GNPNAT1 with its DNA copy alternations (A), methylation (B), and has-miR-30d-3p expression (C). (D) The KM curves of the correlation between has-miR-30d-3p expression and overall survival in LUAD based on TCGA data (Cutoff value of high/low group was the median of miRNA expression). (E) The predicted target region of GNPNAT1 and has-miR-30d-3p.
The miR-30d-3p Downregulation Correlated to GNPNAT1 Overexpression in LUAD
We screened sixteen potential target miRNAs by three databases (Supplementary Figure 3 and Table 5). Subsequently, the above miRNAs’ differential expressions in normal and LUAD samples were assessed by the OncomiR database. We discerned that four miRNAs met our research requirements, including hsa-miR-1323, hsa-miR-30a-3p, hsa-miR-30d-3p, and hsa-miR-30d-3p. Multivariate analysis suggested that only hsa-miR-30d-3p was related to pathologic stage (P = 0.031) and tumor size (P = 0.018) (Table 6). Moreover, we discovered that GNPNAT1 expression had certain negative correlation with hsa-miR-30d-3p expression (R = −0.170, P = 0.0002) (Figure 6C). KM analysis showed that patients with low expression of hsa-miR-30d-3p had a worse prognosis compared with the high expression group (P = 0.044) (Figure 6D). The binding site of hsa-miR-30d-3p and GNPNAT1 at 3′UTR was seen in Figure 6E.
GNPNAT1 Expression Associated With Immune Cell Infiltration
Tumor IMmune Estimation Resource results noted that GNPNAT1 expression was negatively correlated with B cells (R = −0.304, P = 8.04e-12), CD4+ T cells (R = −0.218, P = 1.24e-06), and dendritic cells (R = −0.137, P = 0.002) (Figure 7A). Additionally, MKI67 was upregulated in LUAD (P < 0.01), and we also observed a closely positive relation between GNPNAT1 and MKI67 (Pearson R = 0.548, P < 0.01) (Figures 7B,C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/605754/fmolb-08-605754-HTML/image_m/fmolb-08-605754-g007.jpg)
Figure 7. The role of GNPNAT1 in immunomodulation for LUAD. (A) The association between GNPNAT1 expression and tumor infiltration immune cells. (B) The expression level of MKI67 in LUAD and normal tissues from TCGA. (C) The correlation of GNPNAT1 and MKI67 based on TIMER database.
PPI and GSEA
The PPI network showed that GNPNAT1 had a positive correlation with 25 genes (P < 0.05) (Figure 8A), among which CXCL5 and EIF2S1 had a relatively stronger correlation with GNPNAT1 (RCXCL5 = 0.62, P < 2.2e-16; REIF2S1 = 0.52, P < 2.2e-16) (Figures 8B,C and Supplementary Table 3). Besides, GSEA revealed that upregulated GNPNAT1 was mainly enriched in the cell cycle, ubiquitin-mediated proteolysis, mismatch repair, and p53 signaling pathways (Figures 9A–D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/605754/fmolb-08-605754-HTML/image_m/fmolb-08-605754-g008.jpg)
Figure 8. The PPI network of GNPNAT1 in LUAD. (A) The PPI network was constructed by database. The correlation between GNPNAT1 and CXCL5 (B)/EIF2SI (C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/605754/fmolb-08-605754-HTML/image_m/fmolb-08-605754-g009.jpg)
Figure 9. The GSEA results about GNPNAT1 in LUAD from TCGA. The signaling pathways of the cell cycle (A), ubiquitin mediated proteolysis (B), mismatch repair (C), and p53 signaling pathway (D).
Discussion
Glucose metabolism disorder is a common feature in many cancers, which has become the focused area of anti-tumor drug research (Abdel-Wahab et al., 2019). The principle that tumors have increased demand for glucose has already been widely used by FDG-PET/CT to focus on lesions (Gambhir, 2002). HBP is a branch of glycometabolism in which only about 3–5% of glucose in tumor cells can enter (Akella et al., 2019), but this pathway is the hub of energy metabolism that links glucose, lipid, and protein metabolism. Except for increasing glucose consumption, tumor cells also enhanced the demand for glutamine, an essential substrate for HBP. Studies have found that HBP was involved in the energy metabolism of lung cancer, prostate cancer, and other tumors (Itkonen et al., 2015; Taparra et al., 2018). Furthermore, the final product of HBP, uracil diphosphate n-acetylglucosamine (UDP-GlcNAc), is an essential cell signal regulator and contributes to tumor growth. Also, researchers observed that GNPNAT1 was significantly up-regulated in prostate cancer (Kaushik et al., 2016), and GNPNAT1 low expression led to reduced proliferation of tumor cells in lung cancer receiving chemotherapy (Zhao et al., 2017). However, there has no study regarding the expression and clinical role of GNPNAT1 in LUAD.
In this study, our results indicated that GNPNAT1 was increased in various tumors, including LUAD, and its overexpression was related to advanced staging, lymph node metastasis, and poor prognosis. Additionally, multivariate Cox analysis revealed that GNPNAT1 was an independent prognostic factor in LUAD. The nomograms suggested that GNPNAT1 may play a potential role in clinical diagnosis and prognosis assessment.
Genomic instability is another critical trait of tumors (Hanahan and Weinberg, 2011), so we probed the possible mechanism of GNPNAT1 overexpression in LUAD from epigenetics. GNPNAT1 was located at chromosome 14q22.1, and our study illustrated that the DNA amplification of GNPNAT1 had a close correlation with its overexpression (P < 0.05), which was consistent with one research that the region of 14q22-q24 was significantly amplified in prostate cancer (Bernardino et al., 1997). Furthermore, DNA hypomethylation can also induce genomic instability and play a pro-oncogenic role (Esteller and Herman, 2002), which has been regarded as a prognostic marker in cancers, such as lung cancer and bladder cancer (Moore et al., 2008; Pfeifer and Rauch, 2009). We also noticed that upregulated GNPNAT1 had a positive relation with DNA hypomethylation. Further analysis found that there was a negative correlation between GNPNAT1 overexpression and multiple sites on CpG island, which implied GNPNAT1 might protect tumors from damage through DNA copy number amplification and CpG island hypomethylation.
It was generally reckoned that miRNAs involved multiple tumorigenesis processes and the primary mechanism was that miRNAs could cause the abnormal expression of target genes (Esquela-Kerscher and Slack, 2006). In recent years, the function of miRNA, LncRNA and other non-coding RNAs (ncRNAs) have been gradually established in tumor metabolism. Increasing evidence showed that ncRNAs could directly target glucose transport or consumption or indirectly regulate cancer-related signaling pathways to change glucose metabolism and affect tumor occurrence and development (Shankaraiah et al., 2018). Therefore, we investigated the possible mechanism of GNPNAT1 over-expression in LUAD from its upstream, and results showed that hsa-miR-30d-3p was significantly down-regulated in LUAD and associated with some clinical features. Researchers have reported that hsa-miR-30d-3p was down-regulated in lung cancer and its down-regulation promoted tolerance to EGFR-targeted drug in lung cancer patients (Pan et al., 2019). Additionally, correlation analysis demonstrated that hsa-miR-30d-3p was conversely associated with GNPNAT1 expression and patients with low hsa-miR-30d-3p expression had a poorer prognosis. The above results manifested the down-regulation of hsa-miR-30d-3p may promote GNPNAT1 overexpression by certain pathways, which led to a worse prognosis among LUAD patients.
Tumor microenvironment, the basis of tumor growth and development, was infiltrated by immune cells which, to some extent, could determine the effect of immunotherapy on patients (Li et al., 2016). According to some studies, abnormal glucose metabolism resulted in the accumulation of lactic acid in the tumor microenvironment, promoting tumor growth and invasion by weakening T cells’ activation and migration of monocytes (Gottfried et al., 2006; Fischer et al., 2007). A report at the 20th World Lung Cancer Congress (WCLC) also pointed out that the glucose transporter (GLUT3/GLUT1 ratio) could predict the immunotherapy response of lung cancer (Na et al., 2019). As we know, GNPNAT1 and GLUT1 were both involved in the glucose metabolism pathway, and they had a positive correlation on expression level. Thus, we explored the correlation between GNPNAT1 expression and immune cells through the TIMER database. The analysis revealed that GNPNAT1 had a converse correlation with the infiltration of B cells, CD4+ T cells, and dendritic cells, which all had an anti-tumor effect in NSCLC (Inoshima et al., 2002; Germain et al., 2014; Bruno et al., 2017). GNPNAT1 is one of the key enzymes of HBP process, the end product of which is an essential substrate of glycosylation. Studies suggested that immune cells could express different types of glycosylation-dependent lectin receptors, which affect the function of antigen-presenting cell (APC) and inhibit T-cell activity to promote immune escape by binding to tumor cell surface glycoproteins or glycolipids (RodrÍguez et al., 2018; Sasawatari et al., 2020). Hence, we speculated that the overexpression of GNPNAT1 might affect the function of tumor-infiltrating immune cells by regulating glycosylation modifications. We also noticed a positive correlation between GNPNAT1 and MKI67. Importantly, a previous study has reported that ki67 was associated with decreased immune cells in NSCLC (Mitchell et al., 2019), suggesting that GNPNAT1 may indirectly affect immune function through interacting with MKI67.
Since there were no studies to explore the biological functions of GNPNAT1 overexpression in tumors, we constructed the PPI network which found CXCL5 and EIF2S1 both had a relatively stronger correlation with GNPNAT1. Previous studies have pointed out that CXCL5 was closely related to the invasion and progression of lung cancer (Saintigny et al., 2013), and EIF2S1 could activate autophagy and promote the occurrence of tumors (Dey et al., 2013), so it was reasonable to conjecture that GNPNAT1 might play a synergistic role with these co-expression genes in promoting tumor proliferation.
O-GlcNAcylation, one type of post-translational modification in HBP, was prevalent in tumors, and its high expression level was found to cause p53 instability and promote ubiquitin-mediated proteasomal degradation, thereby leading to resistance to cisplatin-induced apoptosis in lung cancer (Luanpitpong et al., 2017). Recent evidence showed that O-GlcNAcylation promoted mutant KRAS-induced lung tumorigenesis (Taparra et al., 2018). GSEA results also reminded us that GNPNAT1-high group was mainly enriched in the cell cycle, ubiquitin-mediated proteolysis, mismatch repair and p53 signaling pathway, in line with the above study, which proved our GSEA enrichment results were meaningful.
Meanwhile, considering that the partial analysis results of this study are obtained through mining sequencing data, there are several limitations in our study. Firstly, the association between GNPNAT1 and tumor-infiltrating immune cells in lung cancer was acquired via cancer database and bioinformatics analysis methods, it is necessary for researchers to further explore the immune regulatory function of GNPNAT1 in vivo and in vitro experiments. Secondly, the signaling pathways analyzed in this study were discovered based on data mining, experiments were also needed to verify its causal relationship in lung cancer. Finally, clinical samples in the validation set were relatively small, and we will expand the sample size for analysis. In the future, we will investigate the effect of GNPNAT1 on lung cancer cells through invasion and migration experiments in vitro, verify the regulatory relationship between GNPANT1 and miRNA-30d-3p and finally propose to explore the effect of GNPNAT1 on lung cancer through animal models.
In conclusion, GNPNAT1 may serve as a prognostic biomarker and indirectly involve in immune regulation for LUAD. Its upregulation might be regulated by DNA copy number amplification, hypomethylation, and miR-30d-3p down-regulation.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article.
Ethics Statement
The studies involving human participants were reviewed and approved by Tianjin Medical University Cancer Institute and Hospital. The patients/participants provided their written informed consent to participate in this study. 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
WL designed the study, performed the experiment, analyzed the data, and wrote the first draft. KJ participated in the experiment and revised the manuscript. JW assisted in analyzing and revising the manuscript. TM and MZ performed some data analyses. DH directed and supervised the whole work. All participants have read and approved for this manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (No. 81572321).
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.605754/full#supplementary-material
Footnotes
- ^ https://portal.gdc.cancer.gov/
- ^ https://www.ncbi.nlm.nih.gov/geo/
- ^ https://portals.broadinstitute.org/ccle
- ^ https://www.cbioportal.org/
- ^ https://mexpress.be/
- ^ http://www.microrna.gr/webServer
- ^ http://mirwalk.umm.uni-heidelberg.de/
- ^ http://oncomir.org/oncomir/index.html
- ^ http://www.targetscan.org/
- ^ https://cistrome.shinyapps.io/timer/
- ^ http://inetmodels.com
- ^ https://cytoscape.org/
References
Abdel-Wahab, A. F., Mahmoud, W., and Al-Harizy, R. M. (2019). Targeting glucose metabolism to suppress cancer progression: prospective of anti-glycolytic cancer therapy. Pharmacol. Res. 150:104511. doi: 10.1016/j.phrs.2019.104511
Agarwal, V., Bell, G. W., Nam, J. W., and Bartel, D. P. (2015). Predicting effective microRNA target sites in mammalian mRNAs. Elife 4:e05005. doi: 10.7554/eLife.05005
Akella, N. M., Ciraku, L., and Reginato, M. J. (2019). Fueling the fire: emerging role of the hexosamine biosynthetic pathway in cancer. BMC Biol. 17:52. doi: 10.1186/s12915-019-0671-3
Bernardino, J., Bourgeois, C. A., Muleris, M., Dutrillaux, A. M., Malfoy, B., and Dutrillaux, B. (1997). Characterization of chromosome changes in two human prostatic carcinoma cell lines (PC-3 and DU145) using chromosome painting and comparative genomic hybridization. Cancer Genet. Cytogenet. 96, 123–128. doi: 10.1016/s0165-4608(96)00258-0
Boehmelt, G., Wakeham, A., Elia, A., Sasaki, T., Plyte, S., Potter, J., et al. (2000). Decreased UDP-GlcNAc levels abrogate proliferation control in EMeg32-deficient cells. Embo. J. 19, 5092–5104. doi: 10.1093/emboj/19.19.5092
Bruno, T. C., Ebner, P. J., Moore, B. L., Squalls, O. G., Waugh, K. A., Eruslanov, E. B., et al. (2017). Antigen-presenting intratumoral B cells affect CD4(+) TIL phenotypes in non-small cell lung cancer patients. Cancer Immunol. Res. 5, 898–907. doi: 10.1158/2326-6066.Cir-17-0075
Cerami, E., Gao, J., Dogrusoz, U., Gross, B. E., Sumer, S. O., Aksoy, B. A., et al. (2012). The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2, 401–404. doi: 10.1158/2159-8290.Cd-12-0095
Chokchaitaweesuk, C., Kobayashi, T., Izumikawa, T., and Itano, N. (2019). Enhanced hexosamine metabolism drives metabolic and signaling networks involving hyaluronan production and O-GlcNAcylation to exacerbate breast cancer. Cell. Death Dis. 10:803. doi: 10.1038/s41419-019-2034-y
Costello, J. F., Frühwald, M. C., Smiraglia, D. J., Rush, L. J., Robertson, G. P., Gao, X., et al. (2000). Aberrant CpG-island methylation has non-random and tumour-type-specific patterns. Nat. Genet. 24, 132–138. doi: 10.1038/72785
Dey, S., Tameire, F., and Koumenis, C. (2013). PERK-ing up autophagy during MYC-induced tumorigenesis. Autophagy. 9, 612–614. doi: 10.4161/auto.23486
Ehrlich, M. (2002). DNA methylation in cancer: too much, but also too little. Oncogene 21, 5400–5413. doi: 10.1038/sj.onc.1205651
Esquela-Kerscher, A., and Slack, F. J. (2006). Oncomirs - microRNAs with a role in cancer. Nat. Rev. Cancer. 6, 259–269. doi: 10.1038/nrc1840
Esteller, M., and Herman, J. G. (2002). Cancer as an epigenetic disease: DNA methylation and chromatin alterations in human tumours. J. Pathol. 196, 1–7. doi: 10.1002/path.1024
Fischer, K., Hoffmann, P., Voelkl, S., Meidenbauer, N., Ammer, J., Edinger, M., et al. (2007). Inhibitory effect of tumor cell-derived lactic acid on human T cells. Blood 109, 3812–3819. doi: 10.1182/blood-2006-07-035972
Gambhir, S. S. (2002). Molecular imaging of cancer with positron emission tomography. Nat. Rev. Cancer 2, 683–693. doi: 10.1038/nrc882
Germain, C., Gnjatic, S., Tamzalit, F., Knockaert, S., Remark, R., Goc, J., et al. (2014). Presence of B cells in tertiary lymphoid structures is associated with a protective immunity in patients with lung cancer. Am. J. Respir. Crit. Care Med. 189, 832–844. doi: 10.1164/rccm.201309-1611OC
Gottfried, E., Kunz-Schughart, L. A., Ebner, S., Mueller-Klieser, W., Hoves, S., Andreesen, R., et al. (2006). Tumor-derived lactic acid modulates dendritic cell activation and antigen expression. Blood 107, 2013–2021. doi: 10.1182/blood-2005-05-1795
Hanahan, D., and Weinberg, R. A. (2011). Hallmarks of cancer: the next generation. Cell 144, 646–674. doi: 10.1016/j.cell.2011.02.013
Hou, J., Aerts, J., Den Hamer, B., Van Ijcken, W., Den Bakker, M., Riegman, P., et al. (2010). Gene expression-based classification of non-small cell lung carcinomas and survival prediction. PLoS One 5:e10312. doi: 10.1371/journal.pone.0010312
Inoshima, N., Nakanishi, Y., Minami, T., Izumi, M., Takayama, K., Yoshino, I., et al. (2002). The influence of dendritic cell infiltration and vascular endothelial growth factor expression on the prognosis of non-small cell lung cancer. Clin. Cancer Res. 8, 3480–3486.
Itkonen, H. M., Engedal, N., Babaie, E., Luhr, M., Guldvik, I. J., Minner, S., et al. (2015). UAP1 is overexpressed in prostate cancer and is protective against inhibitors of N-linked glycosylation. Oncogene 34, 3744–3750. doi: 10.1038/onc.2014.307
Kaushik, A. K., Shojaie, A., Panzitt, K., Sonavane, R., Venghatakrishnan, H., Manikkam, M., et al. (2016). Inhibition of the hexosamine biosynthetic pathway promotes castration-resistant prostate cancer. Nat. Commun. 7:11612. doi: 10.1038/ncomms11612
Koch, A., Jeschke, J., Van Criekinge, W., Van Engeland, M., and De Meyer, T. (2019). MEXPRESS update 2019. Nucleic Acids Res. 47, W561–W565. doi: 10.1093/nar/gkz445
Lee, S., Zhang, C., Arif, M., Liu, Z., Benfeitas, R., Bidkhori, G., et al. (2018). TCSBN: a database of tissue and cancer specific biological networks. Nucleic Acids Res. 46, D595–D600. doi: 10.1093/nar/gkx994
Li, B., Severson, E., Pignon, J. C., Zhao, H., Li, T., Novak, J., et al. (2016). Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol. 17:174. doi: 10.1186/s13059-016-1028-7
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 77, e108–e110. doi: 10.1158/0008-5472.Can-17-0307
Lu, T. P., Tsai, M. H., Lee, J. M., Hsu, C. P., Chen, P. C., Lin, C. W., et al. (2010). Identification of a novel biomarker, SEMA5A, for non-small cell lung carcinoma in nonsmoking women. Cancer Epidemiol Biomarkers Prev. 19, 2590–2597. doi: 10.1158/1055-9965.Epi-10-0332
Luanpitpong, S., Angsutararux, P., Samart, P., Chanthra, N., Chanvorachote, P., and Issaragrisil, S. (2017). Hyper-O-GlcNAcylation induces cisplatin resistance via regulation of p53 and c-Myc in human lung carcinoma. Sci. Rep. 7:10607. doi: 10.1038/s41598-017-10886-x
Lucena, M. C., Carvalho-Cruz, P., Donadio, J. L., Oliveira, I. A., De Queiroz, R. M., Marinho-Carvalho, M. M., et al. (2016). Epithelial mesenchymal transition induces aberrant glycosylation through hexosamine biosynthetic pathway activation. J. Biol. Chem. 291, 12917–12929. doi: 10.1074/jbc.M116.729236
Mitchell, K. G., Parra, E. R., Nelson, D. B., Zhang, J., Wistuba, I. I., Fujimoto, J., et al. (2019). Tumor cellular proliferation is associated with enhanced immune checkpoint expression in stage I non-small cell lung cancer. J. Thorac. Cardiovasc. Surg. 158, 911–919. doi: 10.1016/j.jtcvs.2019.04.084
Molina, J. R., Yang, P., Cassivi, S. D., Schild, S. E., and Adjei, A. A. (2008). Non-small cell lung cancer: epidemiology, risk factors, treatment, and survivorship. Mayo Clin. Proc. 83, 584–594. doi: 10.4065/83.5.584
Moloughney, J. G., Vega-Cotto, N. M., Liu, S., Patel, C., Kim, P. K., Wu, C. C., et al. (2018). Jacinto E. mTORC2 modulates the amplitude and duration of GFAT1 Ser-243 phosphorylation to maintain flux through the hexosamine pathway during starvation. J. Biol. Chem. 293, 16464–16478. doi: 10.1074/jbc.RA118.003991
Moore, L. E., Pfeiffer, R. M., Poscablo, C., Real, F. X., Kogevinas, M., Silverman, D., et al. (2008). Genomic DNA hypomethylation as a biomarker for bladder cancer susceptibility in the spanish bladder cancer study: a case-control study. Lancet Oncol. 9, 359–366. doi: 10.1016/s1470-2045(08)70038-x
Na, K. J., Choi, H., and Kim, Y. T. (2019). OA08.06 reciprocal change in glucose metabolism of cancer and immune cells mediated by different GLUT predicts immunotherapy response. J. Thoracic Oncol. 14, S227–S228. doi: 10.1016/j.jtho.2019.08.451
Okayama, H., Kohno, T., Ishii, Y., Shimada, Y., Shiraishi, K., Iwakawa, R., et al. (2012). Identification of genes upregulated in ALK-positive and EGFR/KRAS/ALK-negative lung adenocarcinomas. Cancer Res. 72, 100–11. doi: 10.1158/0008-5472.Can-11-1403
Pan, L., Wang, H., Jiang, C., Li, W., Chen, Y., and Ying, G. (2019). Multiple MicroRNAs synergistically promote tolerance to epidermal growth factor receptor-targeted drugs in smoked lung cancer therapies. J. Cancer Res. Ther. 15, 876–881. doi: 10.4103/jcrt.JCRT_208_18
Paraskevopoulou, M. D., Georgakilas, G., Kostoulas, N., Vlachos, I. S., Vergoulis, T., Reczko, M., et al. (2013). DIANA-microT web server v5.0: service integration into miRNA functional analysis workflows. Nucleic Acids Res. 41, W169–W173. doi: 10.1093/nar/gkt393
Pfeifer, G. P., and Rauch, T. A. (2009). DNA methylation patterns in lung carcinomas. Semin. Cancer Biol. 19, 181–187. doi: 10.1016/j.semcancer.2009.02.008
RodrÍguez, E., Schetters, S. T. T., and Van Kooyk, Y. (2018). The tumour glyco-code as a novel immune checkpoint for immunotherapy. Nat. Rev. Immunol. 18, 204–211. doi: 10.1038/nri.2018.3
Saintigny, P., Massarelli, E., Lin, S., Ahn, Y. H., Chen, Y., Goswami, S., et al. (2013). Koo JS. CXCR2 expression in tumor cells is a poor prognostic factor and promotes invasion and metastasis in lung adenocarcinoma. Cancer Res. 73, 571–582. doi: 10.1158/0008-5472.Can-12-0263
Sasawatari, S., Okamoto, Y., Kumanogoh, A., and Toyofuku, T. (2020). Blockade of N-glycosylation promotes antitumor immune response of T cells. J. Immunol. 204, 1373–1385. doi: 10.4049/jimmunol.1900937
Selamat, S. A., Chung, B. S., Girard, L., Zhang, W., Zhang, Y., Campan, M., et al. (2012). Genome-scale analysis of DNA methylation in lung adenocarcinoma and integration with mRNA expression. Genome Res. 22, 1197–211. doi: 10.1101/gr.132662.111
Shankaraiah, R. C., Veronese, A., Sabbioni, S., and Negrini, M. (2018). Non-coding RNAs in the reprogramming of glucose metabolism in cancer. Cancer Lett. 419, 167–174. doi: 10.1016/j.canlet.2018.01.048
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Sticht, C., De La Torre, C., Parveen, A., and Gretztz, N. (2018). miRWalk: An online resource for prediction of microRNA binding sites. PLoS One. 13:e0206239. doi: 10.1371/journal.pone.0206239
Stuart, J. M., Segal, E., Koller, D., and Kim, S. K. (2003). A gene-coexpression network for global discovery of conserved genetic modules. Science 302, 249–255. doi: 10.1126/science.1087447
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U S A. 102, 15545–15550. doi: 10.1073/pnas.0506580102
Taparra, K., Wang, H., Malek, R., Lafargue, A., Barbhuiya, M. A., Wang, X., et al. (2018). O-GlcNAcylation is required for mutant KRAS-induced lung tumorigenesis. J. Clin. Invest. 128, 4924–4937. doi: 10.1172/jci94844
Wong, N. W., Chen, Y., Chen, S., and Wang, X. (2018). OncomiR: an online resource for exploring pan-cancer microRNA dysregulation. Bioinformatics 34, 713–715. doi: 10.1093/bioinformatics/btx627
Ying, H., Kimmelman, A. C., Lyssiotis, C. A., Hua, S., Chu, G. C., Fletcher-Sananikone, E., et al. (2012). Oncogenic Kras maintains pancreatic tumors through regulation of anabolic glucose metabolism. Cell 149, 656–670. doi: 10.1016/j.cell.2012.01.058
Keywords: lung adenocarcinoma, GNPNAT1, prognosis, biomarker, immune infiltration
Citation: Liu W, Jiang K, Wang J, Mei T, Zhao M and Huang D (2021) Upregulation of GNPNAT1 Predicts Poor Prognosis and Correlates With Immune Infiltration in Lung Adenocarcinoma. Front. Mol. Biosci. 8:605754. doi: 10.3389/fmolb.2021.605754
Received: 20 December 2020; Accepted: 23 February 2021;
Published: 25 March 2021.
Edited by:
Zhijie Wang, National Cancer Center of China, ChinaCopyright © 2021 Liu, Jiang, Wang, Mei, Zhao and Huang. 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: Dingzhi Huang, dingzhih72@163.com