A Novel Hypoxia-Related Gene Signature with Strong Predicting Ability in Non-Small-Cell Lung Cancer Identified by Comprehensive Profiling

Yang, Huajun; Wang, Zhongan; Gong, Ling; Huang, Guichuan; Chen, Daigang; Li, Xiaoping; Du, Fei; Lin, Jiang; Yang, Xueyi

doi:https://doi.org/10.1155/2022/8594658

International Journal of Genomics

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 8594658 | https://doi.org/10.1155/2022/8594658

A Novel Hypoxia-Related Gene Signature with Strong Predicting Ability in Non-Small-Cell Lung Cancer Identified by Comprehensive Profiling

Huajun Yang,¹Zhongan Wang,¹Ling Gong,²Guichuan Huang,²Daigang Chen,²Xiaoping Li,²Fei Du,²Jiang Lin,²and Xueyi Yang³

Academic Editor: Mohsen Nikseresht

Received16 May 2021

Revised10 Dec 2021

Accepted16 Apr 2022

Published19 May 2022

Abstract

Background. Non-small-cell lung cancer (NSCLC) is the most common malignant tumor among males and females worldwide. Hypoxia is a typical feature of the tumor microenvironment, and it affects cancer development. Circular RNAs (circRNAs) have been reported to sponge miRNAs to regulate target gene expression and play an essential role in tumorigenesis and progression. This study is aimed at identifying whether circRNAs could be used as the diagnostic biomarkers for NSCLC. Methods. The heterogeneity of samples in this study was assessed by principal component analysis (PCA). Furthermore, the Gene Expression Omnibus (GEO) database was normalized by the affy R package. We further screened the differentially expressed genes (DEGs) and differentially expressed circular RNAs (DEcircRNAs) using the DEseq2 R package. Moreover, we analyzed the Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment of DEGs using the cluster profile R package. Besides, the Gene Set Enrichment Analysis (GSEA) was used to identify the biological function of DEGs. The interaction between DEGs and the competing endogenous RNAs (ceRNA) network was detected using STRING and visualized using Cytoscape. Starbase predicted the miRNAs of target hub genes, and miRanda predicted the target miRNAs of circRNAs. The RNA-seq profiler and clinical information were downloaded from The Cancer Genome Atlas (TCGA) database. Then, the variables were assessed by the univariate and multivariate Cox proportional hazard regression models. Significant variables in the univariate Cox proportional hazard regression model were included in the multivariate Cox proportional hazard regression model to analyze the association between the variables of clinical features. Furthermore, the overall survival of variables was determined by the Kaplan-Meier survival curve, and the time-dependent receiver operating characteristic (ROC) curve analysis was used to calculate and validate the risk score in NSCLC patients. Moreover, predictive nomograms were constructed and used to predict the prognostic features between the high-risk and low-risk score groups. Results. We screened a total of 2039 DEGs, including 1293 upregulated DEGs and 746 downregulated DEGs in hypoxia-treated A549 cells. A549 cells treated with hypoxia had a total of 70 DEcircRNAs, including 21 upregulated and 49 downregulated DEcircRNAs, compared to A549 cells treated with normoxia. The upregulated genes were significantly enriched in 284 GO terms and 42 KEGG pathways, while the downregulated genes were significantly enriched in 184 GO terms and 25 KEGG pathways. Moreover, the function analysis by GSEA showed enrichment in the enzyme-linked receptor protein signaling pathway, hypoxia-inducible factor- (HIF-) 1 signaling pathway, and G protein-coupled receptor (GPCR) downstream signaling. Furthermore, six hub modules and 10 hub genes, CDC45, EXO1, PLK1, RFC4, CCNB1, CDC6, MCM10, DLGAP5, AURKA, and POLE2, were identified. The ceRNA network was constructed, and it consisted of 4 circRNAs, 14 miRNAs, and 38 mRNAs. The ROC curve was constructed and calculated. The area under the curve (AUC) value was 0.62, and the optimal threshold was 0.28. Based on the optimal threshold, the patients were divided into the high-risk score and low-risk score groups. The survival rate in the high-risk score group was lower than that in the low-risk score group. The expression of SERPINE1, STC2, and LPCAT1; clinical stage; and age of the patient were significantly correlated with the high-risk score. Moreover, nomograms were established based on the risk factors in multivariate analysis, and the median survival time, 3-year survival probability, and 5-year survival were possibly predicted according to nomograms. Conclusion. The ceRNA network associated with NSCLC was identified, and the hub genes, circRNAs, might act as the potential biomarkers for NSCLC.

1. Introduction

Lung cancer is the leading cause of cancer-related mortality worldwide [1]. More than 80% of lung cancer cases are diagnosed as non-small-cell lung cancer (NSCLC) [2]. According to the histopathology, there are about 40% adenocarcinomas, 25% squamous cell carcinomas, 10% lung squamous cell carcinomas, large cell carcinomas, and sarcomatoid carcinomas in NSCLC. The main manifestations of NSCLC are cough, expectoration, hemoptysis, wheezing, dyspnea, fever, and weight loss. Intrathoracic metastasis may be accompanied by chest pain, pleural effusion, hoarseness, superior vena cava obstruction syndrome, Horner syndrome, and dysphagia [3]. At present, lung cancer is mainly screened by X-ray and chest computed tomography (CT). The final diagnosis depends on exfoliative sputum cytology, bronchoscopy, video-assisted thoracoscopic surgery, lung biopsy, and surgical procedures for pathological examination [7]. As a result, most NSCLC patients are already in an advanced stage at the time of detection, which seriously affects the treatment of patients [4]. On the other hand, although surgery is advocated in the early stage, and radiotherapy, chemotherapy, targeted therapy, or combined immunotherapy are advocated in the middle and late stages, there are often some problems, such as recurrence after surgery, gene mutation, and drug resistance after chemotherapy and targeted therapy [5, 6]. Moreover, many patients cannot tolerate the toxic effects and side effects, and the 5-year survival rate is still lower than 15% [7, 8]. Therefore, it is very important to better understand the molecular mechanisms underlying NCSLC and explore the potential biomarkers for the diagnosis and treatment of NSCLC.

Hypoxia is an inherent feature of solid tumors because of the imbalance between tumor cell proliferation rates and vascular nutrient supply [9]. Increasing evidence has revealed that hypoxia plays a critical role in tumor occurrence and development, including lung cancer [10]. For example, under the hypoxic microenvironment, tumor cells can undergo epithelial-mesenchymal transition (EMT), thus increasing the degree of malignancy of the tumor [11]. Hypoxia can enhance Wnt signal activity by stabilizing β-catenin and changing its location in the nucleus, promoting proliferation, migration, invasion, and EMT of lung cancer cells, and inhibiting apoptosis [12]. In addition, hypoxia also activates PI3K/Akt and Wnt signaling in a HIF-2 α-dependent manner, thereby enhancing the resistance of lung cancer cells to chronic hypoxia-induced stress, inducing EMT of tumor cells and increasing the malignancy of tumor cells [13]. Furthermore, the expression of mir-191 is upregulated after hypoxia in NSCLC, which leads to a decrease in NFIA expression and promotes the proliferation and migration of lung cancer cells [14]. Under hypoxic conditions, autophagy enhances the antiradiation ability of NSCLC by regulating the level of reactive oxygen species (ROS), affecting the effect of radiotherapy and leading to a poor prognosis of lung cancer [15]. Hypoxia can also upregulate the level of galectin-3, thus enhancing the activity of RhoA and significantly increasing the cell migration and invasion activity [16]. In particular, recent research studies have found that hypoxia is associated with the prognosis of cancers and can act as a target for treatment [17, 18]. However, how hypoxia regulates NSCLC has not been fully elucidated. Therefore, research studies focusing on the molecular mechanism of hypoxia in the occurrence and progression may contribute to the screening of novel biomarkers for the diagnosis and treatment of NSCLC.

Circular RNAs (circRNAs) are a group of noncoding RNAs derived from precursor mRNA back-splicing, and they form a closed covalent circular structure [19]. circRNAs are abundant in humans and have many biological functions [20], including acting as “sponges” of miRNAs that regulate gene expression at the posttranscriptional level [21]. Recent studies have suggested that numerous circRNAs have been reported to play an important role in the proliferation, migration, and invasion of NSCLC cells [7]. For instance, circ_0000284 [19] and circ_0074027 [22] have been reported to promote NSCLC progression by increasing CUL4B through sponging miR-335-5p and by upregulating PD-L1 expression through sponging miR-337-3p, respectively. Upregulated circular RNA vangl1 is involved in NSCLC progression by inhibiting miR-195 and activating Bcl-2 [23]. Besides, abnormal expression of circRNAs is significantly associated with cisplatin resistance in NSCLC, which provides a new therapeutic target for reversing cisplatin resistance in NSCLC [24]. Song et al. have reported that the host genes of circRNAs identified in cisplatin-resistant NSCLC cell lines were involved in biological processes and pathways contributing to cisplatin resistance in NSCLC, indicating the therapeutic role of circRNAs in NSCLC [24]. However, to our knowledge, circRNAs associated with hypoxia in NSCLC have not yet been fully explored.

Therefore, the current study was designed to explore circRNAs in NSCLC associated with hypoxia, followed by constructing a hypoxia-related ceRNA regulatory network and investigating the potential prognostic biomarkers in NSCLC via bioinformatic analysis. Thus, the results of this study will enhance our understanding of the mechanisms underlying hypoxia-related NSCLC and provide novel insights into the diagnosis and treatment of NSCLC.

2. Materials and Methods

2.1. Data Acquisition and Processing

The circRNA expression profiles and transcriptome expression profiles were obtained from the GSE131378 dataset of the GEO database (https://www.ncbi.nlm.nih.gov/geo/). The GSE131378 dataset included 2 normoxia- and 2 hypoxia-treated A549 cells. The RNA-sequencing (RNA-seq) data and paired clinical information from 497 NSCLC samples and 49 standard samples were obtained from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/). The raw data obtained from the Gene Expression Omnibus (GEO) database was normalized using the affy R package for subsequent analysis. To eliminate the system error and verify the reliable data for the subsequent analysis, quality control (QC) of data in this study was assessed by the principal component analysis (PCA).

2.2. Screening of Differentially Expressed Genes (DEGs) and Differentially Expressed circRNAs (DEcircRNAs)

DEGs and DEcircRNAs between normoxia-treated A549 cells and hypoxia-treated A549 cells were screened using the DEseq2 R package with a cutoff value of and a value < 0.05 (FC, fold change).

2.3. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway Enrichment Analysis

GO and KEGG enrichment analysis was performed using the clusterProfiler R package. was used as a threshold to identify significant pathways. GO biological function analysis consisted of biological process (BP), molecular function (MF), and cellular component (CC).

2.4. Gene Set Enrichment Analysis (GSEA)

GSEA performed the gene function analysis of DEGs. The function pathways were assessed between the high- and low-expression groups of DEGs. A threshold and were used as the criteria.

2.5. Protein-Protein Interaction (PPI) Network Construction and Hub Gene Identification

To further explore the interaction among the DEGs and identify the hub genes, a PPI network of DEGs was constructed using the Search Tool for the Retrieval of Interaction Genes (STRING) database (http://string-db.org/). The interaction network was constructed with a . The PPI network was analyzed and visualized using Cytoscape software. The hub modules were identified with a cutoff value , and the hub genes were identified with a cutoff value .

2.6. miRNA of Target mRNA and Target circRNA of miRNA Prediction

The miRNAs targeting hub genes were screened using the Starbase online database (http://starbase.sysu.edu.cn/) with criteria, including , , and . Moreover, the target circRNAs of miRNAs were identified using miRanda online software with a mirSVR score. The lower the score presented, the more reliable the binding site.

2.7. circRNA-miRNA-mRNA (ceRNA) Network Construction

The hypergeometric distribution test was used to assess whether target miRNAs of DEGs were enriched in the miRNA set of circRNAs. Then, the circRNA-miRNA-mRNA network was constructed and visualized using Cytoscape based on the hypergeometric distribution test; .

2.8. Patients and Tissue Samples

Cancer tissues and adjacent tissues of 5 patients with NSCLC who were surgically treated in the Third Affiliated Hospital of Zunyi Medical University from September 2019 to March 2020 were collected. All patients were diagnosed with NSCLC by histopathological examination and did not receive radiotherapy and chemotherapy before the operation. After collection, the specimens were frozen in liquid nitrogen for later use. The study was approved by the Ethics Committee of the Third Affiliated Hospital of Zunyi Medical University, and written informed consent was obtained from all patients or their guardians.

2.9. RNA Isolation and Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR)

All cancer tissues and adjacent tissues were lysed with TRIzol Reagent (Life Technologies-Invitrogen, Carlsbad, CA, USA), and total RNA was isolated following the manufacturer’s instructions. Then, the concentration and purity of the RNA solution were quantified using a NanoDrop 2000FC-3100 nucleic acid protein quantifier (Thermo Fisher Scientific, Waltham, MA, USA Life Real). The extracted RNA was reverse-transcribed to cDNA using the SureScript First-Strand cDNA Synthesis Kit (Genecopoeia, Guangzhou, China) prior to qRT-PCR. The qRT-PCR reaction consisted of 3 μl of the reverse transcription product, 5 μl of 5x BlazeTaq qPCR Mix (Genecopoeia, Guangzhou, China), and 1 μl of forward and reverse primer each. PCR was performed in a Bio-Rad CFX96 Touch™ PCR detection system (Bio-Rad Laboratories, Inc., USA) under the following conditions: initial denaturation at 95°C for 1 min, followed by 40 cycles that involved incubation at 95°C for 20 s, 55°C for 20 s, and 72°C for 30 s. The forward primer of LPCAT1 was “ACCTGCCTAATTACCTTCAAAC,” and the reverse primer of LPCAT1 was “TCCGCAATACCTATCTTCTCTC.” The forward primer of SERPINE1 was “GACTCCCTTCCCCGACTCCA,” and the reverse primer of SERPINE1 was “CGGTCATTCCCAGGTTCTCT.” The forward primer of STC2 was “GTGGGGTGTGGCGTGTTT,” and the reverse primer of STC2 was “TGGGAGGCTTCTGGATGG.” The forward primer of β-actin was “TCCCTGGAGAAGAGCTATGA,” and the reverse primer of β-actin was “AGGAAGGAAGGCTGGAAAAG.” All primers were synthesized by Servicebio (Servicebio, Wuhan, China). The β-actin gene served as an internal control, and the relative expression of 3 hub genes was determined using the 2^-ΔΔCt method [25]. The experiment was performed in triplicate on independent occasions. Statistical differences in the 3 hub genes between the control and NSCLC samples were detected by paired -tests, using GraphPad Prism V6 (GraphPad Software, La Jolla, CA, USA), and the level of statistical significance was tested and presented as .

2.10. Statistical Analysis

The different expression of hub genes simultaneously regulated by multiple circRNAs was determined between standard samples and NSCLC samples of pathologic stages using the one-way ANOVA. Survival curves were performed using the survival R package to assess the association between variables and survival. Candidate variables with were included in the multivariate model. Independent prognostic factors were determined by stepwise regression of multivariate Cox analysis based on the Akaike Information Criterion (AIC) value. The risk score was calculated and normalized using multivariate Cox models. The pROC-R package was then used to create the ROC curve and calculate the AUC and optimal threshold values based on the risk score. NSCLC samples were divided into a high-risk group and a low-risk group. Nomograms of the median survival time and 3-/5-year survival probability were constructed based on multivariate Cox models.

3. Results

3.1. Data Selected and Patients’ Characteristics

The circRNA expression profiles and transcriptome expression profiles were obtained from the GSE131378 dataset. To avoid the system error and verify the reliable data for the subsequent analysis, the data were assessed by PCA. The results showed that differences between groups were more remarkable than those within groups before and after data normalization (Figures 1(a)–1(d)).

(a)

(b)

(c)

(d)

3.2. Identification of DEGs and DEcircRNAs in NSCLC Cells

A total of 2039 DEGs and 70 DEcircRNAs were screened by the DEseq2 R package using the criterion and . The volcano plot map showed significant differences and distribution of the fold change in DEGs (Figure 2(a)) and DEcircRNAs (Figure 2(c)). The heat map illustrated that there were 1293 upregulated DEGs and 746 downregulated DEGs in hypoxia-treated A549 cells compared to normoxia-treated A549 cells (Figure 2(b), Table S1). Moreover, there were 21 upregulated DEcircRNAs and 49 downregulated DEcircRNAs in hypoxia-treated A549 cells compared to normoxia-treated A549 cells, as shown in Figure 2(d) (Table S2).

(a)

(b)

(c)

(d)

3.3. GO and KEGG Functional Enrichment Analyses of DEGs

To further explore the function of DEGs, GO enrichment, including BP term, CC term, and MF term, was analyzed. We found that 284 GO and 42 KEGG pathways were significantly enriched by the upregulated DEGs (Tables S3 and S4). Furthermore, 184 GO and 25 KEGG pathways were significantly enriched by the downregulated DEGs (Tables S5 and S6). As shown in Figure 3(a), the first 30 GO terms for upregulation of DEGs included NADH regeneration, extracellular structural tissue, cell-substrate adhesion, neuronal projection guidance, response to hypoxia, collagen-containing extracellular matrix, desmosomes, extracellular matrix components, myositis, integrin complex, monosaccharide binding, growth factor binding, collagen binding, fibronectin-binding, integrin binding, splinter cell cycle checkpoint, cell cycle DNA replication, nuclear chromosome segregation, DNA integrity checkpoint, negative regulation of cell cycle processes, nucleosome, condensed chromosome kinesis, noncore complex, ribosome, DNA helicase activity, nucleosome binding, protein binding involved in protein folding, chromatin DNA binding, oxidoreductase activity, and CH-OH group acting on the donor (Figure 3(c)). Besides, the 10 KEGG pathways of the upregulated DEGs were as follows: fructose- and mannose-rich metabolism, glycolysis/glycogenesis, biosynthesis of amniotic acid, HIF-1 signaling pathway, focal adhesion, carbon metabolism, carbon metabolism in cancer, protein digestion and absorption, ECM-receptor interaction, and bile secretion (Figure 3(b)). The top 10 KEGG pathways for the downregulated DEGs were as follows: systemic lupus erythematosus-rich, alcoholism, cell cycle, DNA replication, homologous recombination, Fanconi anemia pathway, viral carcinogenesis, necrotizing disease, progesterone-mediated oocyte maturation, and glutathione metabolism (Figure 3(d)).

(a)

(b)

(c)

(d)

3.4. Function Analysis of DEGs by GSEA

To further investigate the function of DEGs, the DEGs were assessed by GSEA with a threshold of and . The results of GSEA showed that the top 30 rich GO terms in the DEGs with the highest enrichment score included enzyme-linked receptor protein signaling pathway, regulation of anatomical structure morphogenesis, cation transport, metal ion transport, regulation of cell motility, cell surface, extracellular matrix, cell-substrate junction, cell-substrate adherens junction, focal adhesion, signaling receptor activity, receptor regulator activity, receptor-ligand activity, growth factor binding, and cytokine receptor binding (Figures 4(a)–4(c)). Ten KEGG pathways were enriched in the HIF-1 signaling pathway, fructose and mannose metabolism, homologous recombination, pentose and glucuronate interconversions, tyrosine metabolism, glutathione metabolism, DNA replication, ascorbate, and alternate metabolism, porphyrin, and chlorophyll metabolism, proteasome (Figure 4(d)). Moreover, the top 10 Reactome terms were enriched in signaling by GPCR, GPCR downstream signaling, extracellular matrix organization, degradation of the extracellular matrix, class B/2 (secretin family receptors), glucose metabolism, collagen formation, gluconeogenesis, snRNP assembly, and metabolism of noncoding RNA (Figure 4(e)).

(a)

(b)

(c)

(d)

(e)

3.5. PPI Network Construction and Hub Gene Identification

We analyzed the interaction between 2039 DEGs and identified the hub genes by constructing a PPI network. The unconnected nodes were removed, and the PPI network included 206 nodes, and 551 edges were constructed (Figure 5(a)). The following hub genes in the PPI network with were identified: CDC45, EXO1, PLK1, RFC4, CCNB1, CDC6, MCM10, DLGAP5, AURKA, and POLE2. Moreover, 6 hub connective modules in the PPI network were identified with connectivity , and 10 hub genes were located in the most extensive module of the six modules (Figure 5(b)).

(a)

(b)

3.6. ceRNA Network Construction

The ceRNA network was constructed to perform further synthetic analysis of the role of hub genes and circRNAs. Firstly, the miRNAs targeting the hub DEGs of 6 hub modules were predicted (Table 1). We then predicted the miRNAs targeting the DEcircRNAs (Table 2). After identifying target miRNAs of DEGs, the miRNA set of circRNAs was determined by the hypergeometric distribution test. The ceRNA network was constructed based on the intersection of hub DEGs and circRNAs, including 4 circRNA nodes, 14 miRNA nodes, and 38 mRNA nodes in hypoxia-treated A549 cells (Figure 6).

3.7. Correlation between Hub Genes and Clinicopathological Characteristics

We further tested the expression of hub DEGs between standard samples and NSCLC samples from TCGA database. The results showed that the expression of ADM, BIRC5, C1QL1, CCNA2, DKK1, FAM160A1, HMGA2, HNRNPA2B1, HOXCB, NECTIN1, PPIH, PRMT3, STC2, and ZWILCH was significantly increased in NSCLC samples compared to standard samples (Figure 7). Inversely, the expression of BHLHE40, C11orf86, CCDC85A, CCND3, DKK3, ETV1, HECA, ISOC1, KDM7A, LPCAT1, PAM, PEA15, PPP1R3B, PTPRE, RASGEF1B, SLC12A2, SNAP25, and WSB1 was strongly decreased in NSCLC samples compared to normal samples (Figure 7). Furthermore, other hub DEGs showed no significant difference between NSCLC samples and normal samples (Supplementary Figure 1). Furthermore, the correlation between hub DEGs and NSCLC patients was analyzed, and the results indicated that high expression of CCDC85A, LPCAT1, PTPRE, SERPINE1, SNAP25, and TCFB1 was correlated with more prolonged survival (Figures 8(c)–8(h)). Besides, clinical features, including age and pathological stages, significantly affected the survival. Patients aged less than 65 years showed more prolonged survival, and patients diagnosed in the I/II stage showed more prolonged survival than those diagnosed in the III/IV stage (Figures 8(a) and 8(b)). However, gender and other hub expression did not significantly affect NSCLC (Supplementary Figure 2).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

3.8. Nomogram Construction and Validation

Candidate variables with a on univariate analysis were included in the multivariable model (Table 3). An ROC curve was drawn using the pROC R package; the AUC value was 0.62, and the optimal threshold was 0.28 (Figure 9(a)). NSCLC samples were grouped into a high-risk group and a low-risk group based on the risk score. Moreover, the survival rate showed a significant difference between the high-risk and low-risk groups, and the high-risk group revealed shorter 5-year survival rates (Figure 9(b)). NSCLC patients were divided into high-risk score and low-risk score groups based on the median risk score, and the number of NSCLC patients with high-risk scores was less than the number of NSCLC patients with low-risk scores (Figure 9(c), top figure). The survival status map showed that the survival rate in the high-risk and low-risk groups declined with time, and most patients in the high-risk score group showed short survival of less than 5 years (Figure 9(c), central figure). Moreover, the heat map indicated that high expression levels of SERPINE1, STC2, and LPCAT1, clinical stage, and age of patients were significantly associated with high-risk scores (Figure 9(c), bottom figure). Moreover, prognostic nomograms were established based on the risk factors in multivariate analysis, and the median survival time, 3-year survival probability, and 5-year survival probability were predicted according to the nomograms (Figure 9(e)).

(a)

(b)

(c)

(d)

(e)

Figure 9

Nomogram construction and validation. (a) The ROC curve for risk prediction. (b) Survival curve for the high-risk score and low-risk score groups. (c) Top figure: the risk score curve. Middle figure: the survival status map. Bottom figure: the heat map of prognostic risk factors (SERPINE1, LPCAT1, STC2, clinical stage, age, and risk score) levels. (d) Nomogram for the median survival time prediction of NSCLC patients. (e) Nomogram for the 3-year and 5-year survival probability prediction in NSCLC patients.

3.9. Validation of the Expression of 3 Hub Genes in NSCLC Patients

We further validated the expression of the 3 above-mentioned hub genes in clinical tissues using qRT-PCR. High expression of SERPINE1 and STC2 in lung cancer tissues of NSCLC patients () and low expression of LPCAT1 in the NSCLC group were confirmed at the mRNA level compared to those in the normal lung tissues () (Figure 10).

(a)

(b)

(c)

4. Discussion

Hypoxia is a typical feature of NSCLC. The molecular mechanisms involved in NSCLC have not yet been completely clarified. In this study, we identified circRNAs using the expression profiling of hypoxia-treated NSCLC cells and constructed the ceRNA regulatory network to further reveal the role of hypoxia in NSCLC. Additionally, we explored the potential hypoxia-related biomarkers for predicting the prognosis of NSCLC patients.

On comparing the expression profiling between hypoxia-treated and nontreated NSCLC cells, we found 2039 DEGs. To investigate the biological functions of these DEGs, we performed GO, KEGG, and GSEA analyses. We found that the upregulated DEGs were significantly enriched in response to hypoxia and the HIF-1 signaling pathway, demonstrating that hypoxia did play an important role in NSCLC. Also, the metastasis-associated protein 2 promotes the metastasis of NSCLC by regulating the ERK/AKT and vascular endothelial growth factor (VEGF) signaling pathways [26]. By regulating the miR-29c/vascular endothelial growth factor (VEGF) signaling pathway, PVT1 promotes angiogenesis in non-small-cell lung cancer (NSCLC) [27]. Thus, these hypoxia-related DEGs may also regulate the progression of NSCLC via the ERK/AKT and VEGF pathways.

Furthermore, we explored the interactions of these DEGs and identified ten hub genes, including LOXL2, STC2, SERPINE1, SLC2A3, PFKFB4, EGLN3, NDRG1, MT1X, PGK1, and ADM. LOXL2 is associated with the progression of lung adenocarcinoma [28] and poor prognosis of NSCLC [29]. Mir-504 can also inhibit the proliferation and invasion of NSCLC cells by targeting LOXL2 [30]. Stc2/Jun/Axl signal activation can mediate the acquired resistance of lung cancer patients to EGFR tyrosine kinase inhibitors [31]. Serpine1 is upregulated in mesenchymal lung cancer cells and promotes cell invasion [32]. Abnormal serpine1 DNA methylation participates in EMT of ovarian cancer induced by carboplatin [33]. Also, serpine1, as an oncogene of gastric adenocarcinoma, promotes tumor cell proliferation, migration, and invasion by regulating EMT [34]. Slc2a3 promotes glycolysis and provides energy for gastric cancer cell proliferation [35]. Pfkfb4 is a biomarker for predicting the poor prognosis of gastric cancer patients [36]. It can promote lung adenocarcinoma progression through phosphorylation and activation of transcription coactivator SRC-2 [37]. Egln2 DNA methylation and expression interact with HIF1A to affect the survival of early NSCLC [38]. Increased expression of NDRG1 in NSCLC is associated with advanced T stage and inadequate angiogenesis [39]. Mtlx can promote the migration and invasion of spc-a-1sci and PC-9 lung cancer cell lines [40]. GBP1 promotes erlotinib resistance in NSCLC through the PGK1-activated EMT signaling pathway [41]. Rab11-FIP2 inhibits the growth of NSCLC by regulating the ubiquitination of PGK1 [42]. Thus, these essential genes can affect the development of NSCLC by affecting the proliferation, migration, invasion, and drug resistance of cancer cells and then affect the survival and prognosis of NSCLC patients.

Meanwhile, we identified XX DEcircRNAs between hypoxia-treated and nontreated NSCLC cells and constructed a ceRNA regulatory network, including 4 circRNAs, 14 miRNA, and 38 mRNAs. These four circRNAs may act as endogenous sponges of the corresponding miRNAs to regulate the expressions of corresponding mRNAs to further affect NSCLC. Moreover, we found that CCDC85A, LPCAT1, PTPRE, SERPINE1, SNAP25, and TCFB1 in the ceRNA network were significantly correlated with NSCLC patients’ survival on univariate Cox analysis. To further obtain more robust biomarkers for predicting the prognosis of NSCLC patients, we performed multivariate Cox analysis and found that SERPINE1, STC2, and LPCAT1 were independent risk factors for NSCLC. Serpine1 is the most closely related gene with invasion, which can be used to evaluate the importance in the invasion and metastasis of NSCLC [43]. Stc2 can be used as an independent prognostic factor to predict the overall survival rate of NSCLC patients [44]. The acquired resistance of lung cancer patients to EGFR tyrosine kinase inhibitors is mediated by reactivation of the stc2/Jun/Axl signal [31]. Lpcat1 promotes brain metastasis of lung adenocarcinoma by upregulating the PI3K/Akt/myc pathway [45].

Finally, we constructed a related risk score that divided NSCLC patients into high- and low-risk groups to accurately predict their clinical outcomes based on the features of SERPINE1, STC2, and LPCAT1.

In conclusion, for the first time, we constructed the ceRNA regulatory network in NSCLC associated with hypoxia and identified SERPINE1, STC2, and LPCAT1 as potential prognostic biomarkers for NSCLC. Further fundamental in vitro and in vivo experiments are needed to provide more solid evidence for our study. Our findings increase our knowledge of the molecular mechanisms involved in hypoxia-related NSCLC and guide the discovery of new therapies for NSCLC.

Data Availability

The data are obtained from GEO and TCGA databases.

Conflicts of Interest

The authors declare no conflict of interests, financial or otherwise.

Acknowledgments

This research was supported by grants from the Science and Technology Bureau Project of Zunyi City (Zunshi Kehe [2019]146 ,[2019]167, [2020]8, [2020]104, and [2020]123).

Supplementary Materials

Supplementary Figure S1: the expressions ofCNKSR3, DGAT2, FAMB1A, SERPINE1, TGFB1, and TMEM132B in the ceRNA network showed no significant difference between NSCLC and normal samples. Supplementary Figure S2: no significant difference of survival between groups divided by gender or expression of ADM, BHLHE40, BIRC5, C1QL1, C11orf86, CCNA2, CCND3, CNKSR3, DKK1, DKK3, DGAT2, ETV1, FAM81A, FAM160A1, HECA, HMGA2, HOXC8, ISOC1, KDM7A, NECTIN1, HNRNPA2B1, PAM, PEA15, PPIH, PPP1R3B, RASGEF1B, SLC12A2, ZWILCH, WSB1, TMEM132B, or STC2. Supplementary Table S1: identification of 1293 upregulated DEGs and 746 downregulated DEGs in hypoxia-treated A549 cells compared to normoxia-treated A549 cells displayed in the heat map. Supplementary Table S2: 21 upregulated DEcircRNAs and 49 downregulated DEcircRNAs identified in hypoxia-treated A549 cells compared to normoxia-treated A549 cells. Supplementary Table S3 and S4: upregulated DEGs were significantly enriched into 284 GO terms and 42 KEGG pathways. Supplementary Table S5 and S6: downregulated DEGs were significantly enriched into 184 GO terms and 25 KEGG pathways. (Supplementary Materials)

References

F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: globocan estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: a Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
View at: Publisher Site | Google Scholar
H. Yin, X. Wang, X. Zhang et al., “UBE2T promotes radiation resistance in non-small cell lung cancer via inducing epithelial-mesenchymal transition and the ubiquitination-mediated FOXO1 degradation,” Cancer Letters, vol. 494, pp. 121–131, 2020.
View at: Publisher Site | Google Scholar
G. C. Brady, J. W. G. Roe, M. O. Brien, A. Boaz, and C. Shaw, “An investigation of the prevalence of swallowing difficulties and impact on quality of life in patients with advanced lung cancer,” Supportive Care in Cancer, vol. 26, no. 2, pp. 515–519, 2018.
View at: Publisher Site | Google Scholar
R. L. Siegel, K. D. Miller, and A. Jemal, “Cancer statistics, 2018,” CA: a Cancer Journal for Clinicians, vol. 68, no. 1, pp. 7–30, 2018.
View at: Publisher Site | Google Scholar
Y. Zhao, H. Wang, and C. He, “Drug resistance of targeted therapy for advanced non-small cell lung cancer harbored egfr mutation: from mechanism analysis to clinical strategy,” Journal of Cancer Research and Clinical Oncology, vol. 147, no. 12, pp. 3653–3664, 2021.
View at: Publisher Site | Google Scholar
T. Lv, Z. Li, L. Xu, Y. Zhang, H. Chen, and Y. Gao, “Chloroquine in combination with aptamer-modified nanocomplexes for tumor vessel normalization and efficient erlotinib/survivin shRNA co-delivery to overcome drug resistance in EGFR-mutated non-small cell lung cancer,” Acta Biomaterialia, vol. 76, pp. 257–274, 2018.
View at: Publisher Site | Google Scholar
C. Wang, H. Huang, Y. Su, C. Tu, T. Hsia, and S. Huang, “Conventional treatment integrated with Chinese herbal medicine improves the survival rate of patients with advanced non-small cell lung cancer,” Complementary Therapies in Medicine, vol. 40, pp. 29–36, 2018.
View at: Publisher Site | Google Scholar
L. A. Torre, F. Bray, R. L. Siegel, J. Ferlay, J. Lortet-Tieulent, and A. Jemal, “Global cancer statistics, 2012,” CA: a Cancer Journal for Clinicians, vol. 65, no. 2, pp. 87–108, 2015.
View at: Publisher Site | Google Scholar
L. H. Gray, A. D. Conger, M. Ebert, S. Hornsey, and O. C. A. Scott, “The concentration of oxygen dissolved in tissues at the time of irradiation as a factor in radiotherapy,” The British Journal of Radiology, vol. 26, no. 312, pp. 638–648, 1953.
View at: Publisher Site | Google Scholar
A. R. Farina, L. Cappabianca, M. Sebastiano, V. Zelli, S. Guadagni, and A. R. Mackay, “Hypoxia-induced alternative splicing: the 11th hallmark of cancer,” Journal of Experimental & Clinical Cancer Research, vol. 39, no. 1, p. 110, 2020.
View at: Publisher Site | Google Scholar
Y. Peng, J. Syu, S. Wang, P. Pan, and H. Kung, “BSA-bounded p-cresyl sulfate potentiates the malignancy of bladder carcinoma by triggering cell migration and emt through the ros/src/fak signaling pathway,” Cell Biology and Toxicology, vol. 36, no. 4, pp. 287–300, 2020.
View at: Publisher Site | Google Scholar
S. Liu, N. Yang, L. Wang, B. Wei, J. Chen, and Y. Gao, “lncRNA SNHG11 promotes lung cancer cell proliferation and migration via activation of Wnt/β-catenin signaling pathway,” Journal of Cellular Physiology, vol. 235, no. 10, pp. 7541–7553, 2020.
View at: Publisher Site | Google Scholar
C. Hong, W. Chen, and C. Wu, “Upregulation of wnt signaling under hypoxia promotes lung cancer progression,” Oncology Reports, vol. 38, no. 3, pp. 1706–1714, 2017.
View at: Publisher Site | Google Scholar
J. Zhao, C. Qiao, Z. Ding et al., “A novel pathway in nsclc cells: mir-191, targeting nfia, is induced by chronic hypoxia, and promotes cell proliferation and migration,” Molecular Medicine Reports, vol. 15, no. 3, pp. 1319–1325, 2017.
View at: Publisher Site | Google Scholar
X. Chen, P. Wang, F. Guo et al., “Autophagy enhanced the radioresistance of non-small cell lung cancer by regulating ros level under hypoxia condition,” International Journal of Radiation Biology, vol. 93, no. 8, pp. 764–770, 2017.
View at: Publisher Site | Google Scholar
Y. Kataoka, Y. Ohshio, K. Teramoto, T. Igarashi, T. Asai, and J. Hanaoka, “Hypoxia-induced galectin-3 enhances rhoa function to activate the motility of tumor cells in non-small cell lung cancer,” Oncology Reports, vol. 41, no. 2, pp. 853–862, 2019.
View at: Google Scholar
P. Gong, Y. Shao, S. Huang et al., “Hypoxia-associated prognostic markers and competing endogenous rna co-expression networks in breast cancer,” Frontiers in Oncology, vol. 10, article 579868, 2020.
View at: Publisher Site | Google Scholar
B. Zhang, B. Tang, J. Gao, J. Li, L. Kong, and L. Qin, “A hypoxia-related signature for clinically predicting diagnosis, prognosis and immune microenvironment of hepatocellular carcinoma patients,” Journal of Translational Medicine, vol. 18, no. 1, p. 342, 2020.
View at: Publisher Site | Google Scholar
L. Li, Q. Zhang, and K. Lian, “Circular rna circ_0000284 plays an oncogenic role in the progression of non-small cell lung cancer through the mir-377-3p-mediated pd-l1 promotion,” Cancer Cell International, vol. 20, no. 1, p. 247, 2020.
View at: Publisher Site | Google Scholar
J. Salzman, R. E. Chen, M. N. Olsen, P. L. Wang, and P. O. Brown, “Cell-type specific features of circular rna expression,” PLoS Genetics, vol. 9, no. 9, article e1003777, 2013.
View at: Publisher Site | Google Scholar
S. Memczak, M. Jens, A. Elefsinioti et al., “Circular RNAs are a large class of animal RNAs with regulatory potency,” Nature, vol. 495, no. 7441, pp. 333–338, 2013.
View at: Publisher Site | Google Scholar
C. Yu, J. Ying, K. Yu, W. Shen, and M. Jiang, “Circ_0074027 contributes to nonsmall cell lung cancer progression by upregulatingcul4b expression through mir-335-5p,” Cancer Biotherapy and Radiopharmaceuticals, vol. 37, no. 2, pp. 1–11, 2020.
View at: Google Scholar
L. Wang, H. Ma, W. Kong, B. Liu, and X. Zhang, “Up-regulated circular rna vangl1 contributes to progression of non-small cell lung cancer through inhibition of mir-195 and activation of bcl-2,” Bioscience Reports, vol. 39, no. 6, 2019.
View at: Google Scholar
L. Song, Z. Cui, and X. Guo, “Comprehensive analysis of circular rna expression profiles in cisplatin-resistant non-small cell lung cancer cell lines,” Acta Biochimica et Biophysica Sinica, vol. 52, no. 9, pp. 944–953, 2020.
View at: Publisher Site | Google Scholar
J. Xu, W. Wu, C. Wu et al., “A large-scale, multicentered trial evaluating the sensitivity and specificity of digital pcr versus arms-pcr for detecting ctdna-based egfr p.t790m in non-small-cell lung cancer patients,” Translational Lung Cancer Research, vol. 10, no. 10, pp. 3888–3901, 2021.
View at: Publisher Site | Google Scholar
B. Zhang, F. Tao, and H. Zhang, “Metastasis-associated protein 2 promotes the metastasis of non-small cell lung carcinoma by regulating the erk/akt and vegf signaling pathways,” Molecular Medicine Reports, vol. 17, no. 4, pp. 4899–4908, 2018.
View at: Publisher Site | Google Scholar
Z. Mao, B. Xu, L. He, and G. Zhang, “PVT1 promotes angiogenesis by regulating mir-29c/vascular endothelial growth factor (vegf) signaling pathway in non-small-cell lung cancer (nsclc),” Medical Science Monitor, vol. 25, pp. 5418–5425, 2019.
View at: Publisher Site | Google Scholar
P. Zhan, X. Shen, Q. Qian et al., “Down-regulation of lysyl oxidase-like 2 (loxl2) is associated with disease progression in lung adenocarcinomas,” Medical Oncology, vol. 29, no. 2, pp. 648–655, 2012.
View at: Publisher Site | Google Scholar
P. Zhan, X. Lv, Y. Ji, H. Xie, and L. Yu, “Increased lysyl oxidase-like 2 associates with a poor prognosis in non-small cell lung cancer,” The Clinical Respiratory Journal, vol. 12, no. 2, pp. 712–720, 2018.
View at: Publisher Site | Google Scholar
M. Ye, J. Zhang, T. Guo, and X. Pan, “MiR-504 inhibits cell proliferation and invasion by targeting LOXL2 in non small cell lung cancer,” Biomedicine & Pharmacotherapy, vol. 97, pp. 1289–1295, 2018.
View at: Publisher Site | Google Scholar
Y. N. Liu, M. F. Tsai, S. G. Wu et al., “Acquired resistance to EGFR tyrosine kinase inhibitors is mediated by the reactivation of STC2/JUN/AXL signaling in lung cancer,” International Journal of Cancer, vol. 145, no. 6, pp. 1609–1624, 2019.
View at: Publisher Site | Google Scholar
H. J. Kong, E. J. Kwon, O. S. Kwon et al., “Crosstalk between yap and tgfβ regulates serpine1 expression in mesenchymal lung cancer cells,” International Journal of Oncology, vol. 58, no. 1, pp. 111–121, 2020.
View at: Publisher Site | Google Scholar
J. Pan, F. Qu, F. Wang et al., “Aberrant serpine1 dna methylation is involved in carboplatin induced epithelial-mesenchymal transition in epithelial ovarian cancer,” Archives of Gynecology and Obstetrics, vol. 296, no. 6, pp. 1145–1152, 2017.
View at: Publisher Site | Google Scholar
J. Yang, L. Ma, and Z. Zhu, “SERPINE1 as a cancer-promoting gene in gastric adenocarcinoma: facilitates tumour cell proliferation, migration, and invasion by regulating emt,” Journal of Chemotherapy, vol. 31, no. 7-8, pp. 408–418, 2019.
View at: Publisher Site | Google Scholar
D. Chen, H. Wang, J. Chen et al., “MicroRNA-129-5p regulates glycolysis and cell proliferation by targeting the glucose transporter slc2a3 in gastric cancer cells,” Frontiers in Pharmacology, vol. 9, p. 502, 2018.
View at: Publisher Site | Google Scholar
F. Wang, X. Wu, Y. Li, X. Cao, C. Zhang, and Y. Gao, “PFKFB4 as a promising biomarker to predict a poor prognosis in patients with gastric cancer,” Oncology Letters, vol. 21, no. 4, p. 296, 2021.
View at: Publisher Site | Google Scholar
J. Meng, X. Chen, and Z. Han, “PFKFB4 promotes lung adenocarcinoma progression via phosphorylating and activating transcriptional coactivator src-2,” BMC Pulmonary Medicine, vol. 21, no. 1, p. 60, 2021.
View at: Publisher Site | Google Scholar
R. Zhang, L. Lai, J. He et al., “EGLN2DNA methylation and expression interact withHIF1Ato affect survival of early-stage nsclc,” Epigenetics, vol. 14, no. 2, pp. 118–129, 2019.
View at: Publisher Site | Google Scholar
C. Fan, J. Yu, Y. Liu, H. Xu, and E. Wang, “Increased ndrg1 expression is associated with advanced t stages and poor vascularization in non-small cell lung cancer,” Pathology & Oncology Research, vol. 18, no. 3, pp. 549–556, 2012.
View at: Publisher Site | Google Scholar
Y. Liu, L. Liu, T. Yu et al., “Systematic analysis of mrna expression profiles in nsclc cell lines to screen metastasis-related genes,” Molecular Medicine Reports, vol. 14, no. 6, pp. 5093–5103, 2016.
View at: Publisher Site | Google Scholar
L. Cheng, L. Gou, T. Wei, and J. Zhang, “GBP1 promotes erlotinib resistance via pgk1-activated emt signaling in non-small cell lung cancer,” International Journal of Oncology, vol. 57, no. 3, pp. 858–870, 2020.
View at: Publisher Site | Google Scholar
W. Dong, H. Li, and X. Wu, “Rab11-FIP2 suppressed tumor growth via regulation of PGK1 ubiquitination in non-small cell lung cancer,” Biochemical and Biophysical Research Communications, vol. 508, no. 1, pp. 60–65, 2019.
View at: Publisher Site | Google Scholar
A. S. Lader, M. F. Ramoni, B. R. Zetter, I. S. Kohane, and D. J. Kwiatkowski, “Identification of a transcriptional profile associated with in vitro invasion in non-small cell lung cancer cell lines,” Cancer Biology & Therapy, vol. 3, no. 7, pp. 624–631, 2004.
View at: Google Scholar
R. Li, X. Liu, X. J. Zhou et al., “Identification and validation of the prognostic value of immune-related genes in non-small cell lung cancer,” American Journal of Translational Research, vol. 12, no. 9, pp. 5844–5865, 2020.
View at: Google Scholar
C. Wei, X. Dong, H. Lu et al., “LPCAT1 promotes brain metastasis of lung adenocarcinoma by up-regulating pi3k/akt/myc pathway,” Journal of Experimental & Clinical Cancer Research, vol. 38, no. 1, p. 95, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Huajun Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

551

Downloads

871

Citations