Next Article in Journal
The Rise of Retinal Organoids for Vision Research
Previous Article in Journal
Seroconversion and Skin Mucosal Parameters during Koi Herpesvirus Shedding in Common Carp, Cyprinus carpio
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 22
10.3390/ijms21228479
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Establishment of a Gene Signature to Predict Prognosis for Patients with Lung Adenocarcinoma

by
Zhaodong Li
1,
Fangyuan Qi
1 and
Fan Li
1,2,3,4,5,*
1
The Key Laboratory of Zoonosis, Department of Pathogenobiology, College of Basic Medicine, Jilin University, Changchun 130021, China
2
The Key Laboratory for Bionics Engineering, Jilin University, Changchun 130021, China
3
Engineering Research Center for Medical Biomaterials of Jilin Province, Jilin University, Changchun 130021, China
4
Key Laboratory for Biomedical Materials of Jilin Province, Jilin University, Changchun 130021, China
5
State Key Laboratory of Pathogenesis, Prevention and Treatment of HighIncidence Diseases in Central Asia, Urumqi 830000, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(22), 8479; https://doi.org/10.3390/ijms21228479
Submission received: 28 August 2020 / Revised: 31 October 2020 / Accepted: 4 November 2020 / Published: 11 November 2020
(This article belongs to the Section Molecular Genetics and Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:
Accumulating evidence indicates that the reliable gene signature may serve as an independent prognosis factor for lung adenocarcinoma (LUAD) diagnosis. Here, we sought to identify a risk score signature for survival prediction of LUAD patients. In the Gene Expression Omnibus (GEO) database, GSE18842, GSE75037, GSE101929, and GSE19188 mRNA expression profiles were downloaded to screen differentially expressed genes (DEGs), which were used to establish a protein-protein interaction network and perform clustering module analysis. Univariate and multivariate proportional hazards regression analyses were applied to develop and validate the gene signature based on the TCGA dataset. The signature genes were then verified on GEPIA, Oncomine, and HPA platforms. Expression levels of corresponding genes were also measured by qRT-PCR and Western blotting in HBE, A549, and PC-9 cell lines. The prognostic signature based on eight genes (TTK, HMMR, ASPM, CDCA8, KIF2C, CCNA2, CCNB2, and MKI67) was established, which was independent of other clinical factors. The risk model offered better discrimination between risk groups, and patients with high-risk scores tended to have poor survival rate at 1-, 3- and 5-year follow-up. The model also presented better survival prediction in cancer-specific cohorts of age, gender, clinical stage III/IV, primary tumor 1/2, and lymph node metastasis 1/2. The signature genes, moreover, were highly expressed in A549 and PC-9 cells. In conclusion, the risk score signature could be used for prognostic estimation and as an independent risk factor for survival prediction in patients with LUAD.

Graphical Abstract

1. Introduction

Lung cancer is still the deadliest cancer type worldwide. More than 1.6 million new patients are diagnosed each year with lung cancer, and the disease is associated with poor quality of life and high expenditure [1,2]. Lung adenocarcinoma (LUAD) is considered to be the most frequent type of lung cancer, and approximately 40% of lung cancer patients have LUAD [3]. Despite advances in cancer therapy in the past years, the prognosis of LUAD patients is still unsatisfactory, with poor overall survival (OS) rate. It is well known that the genetic alterations have critical functions on critical biological pathways during LUAD development and progression, and further promote the recurrence of cancer as well as lower the survival rate of patients [4,5]. Thus, it is necessary to identify reliable predictors for prognostic estimation, which could bring tremendously guiding value in the management of the LUAD. Especially, the construction of multiple-gene signature would benefit LUAD patients from prognostic prediction.
The widespread application of gene expression profiles and the availability of public genomic datasets enable researchers to study and elucidate the underlying molecular mechanisms of diseases at genome and transcriptome levels [6,7,8]. With the integration of genomics technology and bioinformatics analysis, gene microarrays have been widely applied to collect and study gene expression profiles in several cancer types [9,10]. While these tools allowed researchers to screen for tumor-associated genes and identify the core prognosis factors, a single gene biomarker cannot offer efficient survival prediction. A risk model based on multiple genes may offer a better prognostic factor for predicting the survival of patients [11,12]. In this study, we aimed to establish a multiple-gene risk score signature via systematic bioinformatics analysis.
Here, we obtained four gene expression profiles from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/) and used these profiles to collect differentially expressed genes (DEGs) and to identify the signature genes. We then constructed protein-protein interaction (PPI) and subsequently performed a clustering module and biological pathways enrichment analyses. We adopted univariate and multivariate Cox proportional hazards regression analyses to construct the risk score model, and then used survival analyses as well as the receiver operating characteristic (ROC) curve to verify the sensitivity and specificity of this model. We also illustrate the prognostic value of the risk score signature. The detailed workflow of this work is provided in Figure 1.

2. Results

2.1. Identification of 573 DEGs Shared by Four GEO Profiles

To screen the DEGs between LUAD and normal lung tissues, we utilized the GSE18842, GSE75037, GSE101929, and GSE19188 gene expression profiles from the GEO database. The clinical characteristics of lung cancer patients contained in the Cancer Genome Atlas (TCGA) and these four GEO datasets are shown in Table 1. It is noteworthy that data on some clinical characteristics (particularly T classification, N classification, and M classification) were not provided for the four GEO datasets. The DEGs were screened in each of the four GEO datasets, respectively. First, we found 1840 DEGs in GSE18842, 2856 DEGs in GSE19188, 1959 DEGs in GSE75037, and 1458 DEGs in GSE101929, respectively. Of them, 812, 1313, 792, and 553 genes were upregulated and 1028, 1543, 1167, and 905 genes were downregulated in GSE18842, GSE19188, GSE 75037, and GSE101929 datasets, respectively (Figure 2A–D). Then, a total of 573 overlapping DEGs (197 upregulated and 376 downregulated)were screened among these four profiles (Figure 2E,F), which were used for the following analyses.These common DEGs significantly forcused on mitotic nuclear division (biological process, BP),extracellular matrix (cellular component, CC),and extracellular matrix structural constituent (molecular function, MF), as well as other functional pathways (Supplementary Figure S1A). Moreover, cytokine-cytokine receptor interaction, IL-17 signaling pathway, p53 signaling pathway and cell cycle were also main signaling pathways enriched by the 573 overlapping DEGs (Supplementary Figure S1B).

2.2. Candidate Core Genes Identification

STRING and Cytoscape 3.6.0 software programs were used to construct PPI network for the 573 overlapping DEGs and the network included 502 nodes and 4320 edges (Figure 3A). Moreover, we obtained four clustering modules from Cytoscape’s MCODE plug-in. Supplementary Table S1 has shown that Cluster 1 included 72 nodes and 2326 edges with the highest MCODE score 65.521. Furthermore, it significantly focused on cell cycle, DNA replication, Mitotic M-M/G1 phases, M Phase, and other signaling pathways (adjusted p < 0.05) (Figure 3B). Meanwhile, other cluster modules were mainly enriched in cell surface interactions at the vascular wall, Hemostasis, Beta3 integrin cell surface interactions, and Peptide ligand-binding receptors (adjusted p < 0.05) (Figure 3C–E). More details are provided in Supplementary Table S2. In this study, clustering module 1 was considered as the key module. Then, the top 30 genes with largest degrees of connectivity in cluster 1 were shown in Table 2. We selected the first 20 genes with the highest degrees of connectivity as the candidate core genes for further study.

2.3. The Risk Model Based on the Eight Genes is Verified as an Independent Prognosis Factor

To establish the gene expression signature, we performed univariate and multivariate Cox regression analyses for the 20 candidate core genes and identified eight DEGs (TTK, HMMR, ASPM, CDCA8, KIF2C, CCNB2, MKI67, and CCNA2), as shown in Table 3. The risk score for OS was calculated as follows:
Risk score = ((−0.462) × TTK expression level) + (0.633 × HMMT expression level)
+ ((−0.550) × ASPM expression level) + ((−1.309) × CDCA8 expression level)
+ (1.188 × KIF2C expression level) + ((−0.474) × CCNB2 expression level)
+ (0.570 × MKI67 expression level) + (0.484 × CCNA2 expression level).
In the training group, LUAD patients were divided into high-risk and low-risk groups according to median risk score, which was used as cut-off value. The survival analysis revealed that the survival rate was distinctively lower in the high-risk group than low-risk group (p < 0.01), as shown in Figure 4A. Additionally, the risk model offered a survival prediction with AUCs: 0.706, 0.748, and 0.698 at 1-, 3- and 5-yearfollow-up, respectively (Figure 4B). The heatmap indicated that eight signature genes tended to have higher expression in high-risk patients (Figure 4C). The distribution of LUAD patients was plotted based on median cut-off value (Figure 4D). The survival status of patients in the training group was illustrated in Figure 4E. Meanwhile, we observed similar results in the testing group (Figure 5). Moreover, our risk model could be used as an independent prognosis predictor for LUAD patients in the training (p < 0.01) and testing (p < 0.001) groups (Table 4) in multivariate Cox analysis.
We further analyzed the TCGA dataset to evaluate the prognostic value of the risk score signature in LUAD patients (Figure 6). We found that the high-risk group was distinctively associated with worse OS in >65 age (p = 0.008), male (p = 0.023), stage III/IV (p = 0.008), N1-3 (p = 0.024), and T1-2 (p = 0.027) subgroups. Then, the ROC analysis verified the accuracy of our gene signature in OS prediction in subgroups with different clinical characteristics (Supplementary Figure S2).The values of AUCs at 1, 3, and 5 years were between 0.965 and 0.610, which indicated that the risk model had a robust performance to predict OS probabilities for LUAD patients.

2.4. Signature Genes Show High Expression in LUAD Samples

The results from our correlation analysis showed that the risk score was positively associated with mRNA expression levels of all signature genes (p < 0.0001; Figure 7A), which was significantly higher in LUAD tissues compared with normal lung tissues (Figure 7B). We also found that TTK, HMMR, CDCA8, CCNA2, CCNB2, and MKI67 proteins were over-expressed in LUAD tissues versus normal tissues (Supplementary Figure S3A). Compared to normal tissues, ASPM and KIF2C expressions were significantly elevated in lung cancer tissues within three datasets (Hou Lung, Selamat Lung, and Okayama Lung) [13,14,15] and five datasets (Hou Lung, Landi Lung, Okayama Lung, Stearman Lung, and Su Lung) [13,15,16,17,18], respectively (Supplementary Figure S3B).

2.5. Signature Genes with High Expression in LUAD Cells

Compared with human bronchial epithelial (HBE) cells, the mRNA levels of TTK, HMMR, ASPM, CDCA8, KIF2C, CCNB2, MKI67, and CCNA2 were significantly upregulated in A549 and PC-9 cell lines(p < 0.0001, Figure 8A). Meanwhile, TTK, HMMR, ASPM, CDCA8, KIF2C, CCNB2, MKI67, and CCNA2 proteins were also significantly over-expressed in A549 and PC-9 compared to HBE cells(p < 0.05, Figure 8B, Supplementary Figure S4).These results were consistent with expression levels of riskscore signature genes in LUAD tissues.

3. Discussion

Currently, LUAD is regarded as the most common subtype of lung cancer with a five-year OS rate between 4% and 17% [3]. The pathogenesis of LUAD involves a diverse and complex set of events, such as gene expression, autophagy activation, unexpected tumor microenvironment, immune cell infiltration, abnormal cell cycle, DNA methylation, epigenetic interactions, and other molecular and cellular events [19,20,21,22,23]. Moreover, accumulating evidence suggests that the corresponding signaling pathways, such as the MAPK, AKT-mTOR, and Wnt/β-catenin signaling pathways, are more frequently activated in LUAD [24,25,26]. While research into the tumor mechanism and treatment has made significant progress, the prognosis still remains poor in LUAD. Hence, it is urgent to identify the precise and effective prognostic signature to predict the survival of patients with LUAD.
Several risk score models based on gene expression for LUAD prognosis have been identified in the previous studies. The 9-gene signature (HMMR, B4GALT1, SLC16A3, ANGPTL4, EXT1, GPC1, RBCK1, SOD1, and AGRN) has been identified as an independent prognostic factor for OS in LUAD patients and is significantly associated with metastasis in the test series [2]. A 4-gene signature (CTNNB1 or β-catenin, SOX9, DVL3, and Wnt2b) involved in Wnt/β-catenin pathway can significantly divide LUAD patients into high- and low-risk groups with different OS rates [27]. An immune signature consisting of 40 genes can effectively differentiate high- and low-risk groups among patients with stage I or II lung adenocarcinoma and independently predict OS [28]. Meanwhile, the clinical characteristics such as age and stage can improve the prognostic accuracy of the immune signature. A 6-gene signature (RRAGB, RSPH9, RPS6KL1, RXFP1, RRM2, and RTL1) can significantly stratify LUAD patients into high- and low-risk groups, and remain as an independent prognostic factor to estimate OS in a multivariate Cox proportional hazards model analysis [11]. An 8-gene signature (DLGAP5, KIF11, RAD51AP1, CCNB1, AURKA, CDC6, OIP5, and NCAPG) has been constructed by a multivariate Cox regression model that can clearly separate LUAD patients into groups with significantly different OS rates [29]. Our 8-gene prognostic signature (TTK, HMMR, ASPM, CDCA8, KIF2C, MKI67, CCNA2, and CCNB2) was different from these previously identified gene signatures, indicating it might be a novel tool for LUAD prognosis.
TTK is a mitotic checkpoint kinase that is highly expressed in several human cancers [30]. A high TTK expression, moreover, is positively associated with higher grade aggressiveness and therapeutic resistance in breast cancer, implicating the TTK might be involved in cancer cell proliferation and poor patient survival rate [31]. HMMR is also highly expressed in human multiple malignancies, such as gastric cancer and lung cancer, which may promote cancer progression and lower patient survival rate [2,32]. Ye et al. suggested that TGF-β signaling and Hippo pathway could contribute to sarcomagenesis and metastasis by upregulating the HMMR expression [33]. ASPM, the novel Wnt co-activator, augments the Wnt-β-catenin signaling to maintain the subpopulation of cancer stem cells in prostate cancer [34]. The expression of ASPM is incrementally upregulated in primary and metastatic lung cancer, indicating its potential roles in the occurrence and progression of lung cancer [35]. CDCA8, on the other hand, is involved in cell cycle regulation [36]. For example, CDCA8 suppression leads to cell cycle G1 phase arrest, upregulations of p21 and p27 expressions, and the downregulations of CCND1 and Bcl-2 expressions [37]. Additionally, high CDCA8 expression showed a significant correlation with poor survival in patients with cutaneous melanoma [38]. Results from integrated analyses show that KIF2C is overexpressed in several solid cancers and has already been identified as an important prognostic factor for cancers [39,40]. KIF2C, a mitotic centromere-associated kinase, acts on cellular senescence mainly via p53 signaling; it frequently induced T cell responses in certain cancers [39,41]. CCNA2 has critical functions in controlling cell cycle at the G1/S and the G2/M transitions and is necessary for embryonic cell and hematopoietic lineage [42]. It has been reported that CCNA2 might be associated with the progression of epithelial-mesenchymal transitions (EMT) and metastasis [43]. CCNA2 expression, moreover, is significantly upregulated in many cancer types according to the Human Protein Atlas database, implying its potential role in cancer development and progression. Similarly, CCNB2 is also a key regulator of the cell cycle and may play a role in the development and progression of cancers in humans [44,45,46]. In this study, we established the risk score signature based on these eight genes and elucidated the prognostic value of the gene signature for LUAD.
Despite the encouraging predictive value in LUAD prognosis, our risk model still has some limitations. First, since our model is based on retrospective datasets and these datasets may have unbalanced clinical features with treatment heterogeneity, its efficacy needs to be validated in clinical trials with enough LUAD samples. Second, our risk model could not accurately predict survival possibilities for LUAD patients in <=65 age, female, Stage I and II, N0, and T3-4subgroups of LUAD patients, which might be due to limited samples. Meanwhile, the values of AUC at 1, 3, and 5 years were lower than 0.8, indicating that we should try to establish a new and powerful risk score model in the further analyses. Finally, the signaling pathways involved in the progression of LUAD need to be examined more thoroughly, which is beyond the scope of our current study.

4. Materials and Methods

4.1. Data Collection

Four gene expression profiles including GSE18842, GSE75037, GSE101929, and GSE19188 were obtained from the GEO online public database. GSE18842 contained 45 normal lung samples and 46 LUAD samples; GSE75037 contained 83 normal lung samples and 83 LUAD samples; GSE101929 contained 34 normal lung samples and 32 LUAD samples; and GSE19188 contained 65 normal lung samples and 91 LUAD samples.

4.2. Data Processing and Identification of DEGs

The normalization and log2 conversion for each GEO matrix data were performed by the limma package in R software (version 3.5.3). The limma package was applied to screen the DEGs in GSE18842, GSE75037, GSE19188, and GSE101929 profiles. The adjusted p-value < 0.05 and |log2FC| ≥ 1 were set as the screening criteria for DEGs.

4.3. PPI Network Construction and Clustering Module Analysis

The information regarding interactions of proteins was retrieved from the Search Tool for the Retrieval of Interacting Genes database (STRING) (https://string-db.org/), and a combined score of ≥0.4 was used as the cut-off criterion. Next, a protein-protein interaction (PPI) network was constructed using Cytoscape software (version 3.6.0) (https://cytoscape.org/). To glean specific biological information, molecularcomplex detection (MCODE), a Cytoscape plug-in, was used to cluster modules from the PPI network with the following parameters: degree cutoff = 2, Node Score cutoff = 0.2, and K-Core = 2. The enrichment analysis for the clustering modules was investigated with FunRich software (version 3.1.3).

4.4. Establishment and Validation of a Prognostic Signature

Patients in the Cancer Genome Atlas (TCGA) dataset were randomly divided into either a training group (228 LUAD cases) or a testing group (230 LUAD cases). After univariate Cox regression analysis, the DEGs that were significantly associated with OS were subjected to the least absolute shrinkage and selection operator (LASSO) regression analysis. The analysis was performed with the R “glmnet” package to screen out candidate prognostic genes in the training group. Subsequently, the candidate prognostic genes were included in the multivariate Cox regression analysis to establish the risk score formula, which could separate the high- and low-risk subgroups according to median risk score in patients with LUAD. Kaplan–Meier survival analysis and ROC curve were performed to determine the predictive value of the risk model. Univariate and multivariate Cox regression analyses were performed on the clinical data and risk score to evaluate whether the risk model was an independent prognosis factor from clinical parameters; the clinical relevance of the gene expression signature was also validated.

4.5. Validation of Corresponding Genes in Risk Score Signature

The Gene Expression Profiling Interactive Analysis (GEPIA) database (http://gepia.cancer-pku.cn) was adopted to illustrate the mRNA expression levels of signature genes in LUAD and normal lung tissues. The immunohistochemical results of corresponding genes from the Human Protein Atlas (HPA) (https://www.proteinatlas.org) were obtained to elucidate the protein expression levels of these genes. Furthermore, a meta-analysis was performed on the Oncomine database (https://www.oncomine.org) to identify the expression patterns of certain genes.

4.6. Cell Culture

HBE, A549, and PC-9 cells were grown in high glucose Dulbecco’s Modified Eagle’s media (DMEM; Hyclone, Logan, UT, USA) containing 10 % (v/v) fetal bovine serum (FBS; Gibco, Grand Island, NY, USA), and 1 % penicillin/streptomycin (MRC, Jintan, China) at 37 °C and 5% CO2. All cell lines were obtained from the Shanghai Cell Bank of Chinese Academy of Medical Sciences (Shanghai, China).

4.7. Quantitative Real-Time Polymerase Chain Reaction

Total RNA was extracted by using Total RNA Extraction Kit (Solarbo, Beijing, China), and reverse transcription was performed using the first-strand cDNA synthesis kit (Invitrogen, Carlsbad, CA, USA) following the manufacturers’ instructions. Real-time PCR was performed using Premix Ex Taq SYBR Green PCR (TaKaRa, Dalian, China) on an ABI PRISM 7300 (Applied Biosystems, Foster City, CA, USA) following the manufacturer’s guidelines. The primer sequences used are given in Supplementary Table S3.

4.8. Western Blot Analysis

After extraction of cell lysate, protein concentrations were quantified by the Bradford method (Beyotime, Shanghai, China) and proteins were separated on5–12% SDS-PAGE gels, which were then transferred to PVDF membranes (ThermoFisher, Waltham, MA, USA). The membranes were blocked with 5 % non-fat milk, incubated with primary antibody at 4 °C overnight, treated with horseradish peroxidase-conjugated secondary antibody (Bioss, Beijing, China) at room temperature for 1 h, and visualized on a Tanon 5200 (Tanon, Shanghai, China). Primary antibodies were purchased from Abcam (Cambridge, UK), and the following antibodies were used: TTK, HMMR, ASPM, CDCA8, KIF2C, CCNB2, MKI67, CCNA2, and β-actin. More information about primary antibodies was provided in Supplementary Table S4.

4.9. Statistical Analysis

All statistical analysis in this work was conducted using R software (version 3.6.0) and GraphPad Prism 7.0. The correlations between risk scores and expression levels of prognostic signature genes were analyzed via Pearson’s correlation test. Kaplan–Meier analysis was performed to estimate survival, and a log–rank test was used to compare between-group survival distributions. The predictive accuracy of our risk model was assessed by ROC analysis with the R “SurvivalROC” package. The univariate and multivariate Cox proportional hazards regression models were used to calculate regression coefficients and hazard ratios and to establish the risk score model. Two-tailed p < 0.05 was considered to be statistically significant.

5. Conclusions

We established a risk score signature for survival prediction in LUAD patients where patients with high-risk scores exhibited a poor survival rate. To some extent, our model may improve the prognostic accuracy in clinical application. Further studies are needed to improve the precision and reliability of the signature before our model can be implemented in clinical settings.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/22/8479/s1.

Author Contributions

Conceptualization, F.L. and Z.L.; methodology, Z.L. and F.Q.; software, Z.L.; validation, Z.L. and F.Q.; formal analysis, Z.L.; investigation, Z.L. and F.Q.; resources, F.Q.; data curation, F.L.; writing—original draft preparation, Z.L.; writing—review and editing, F.L. and F.Q.; visualization, Z.L. and F.Q.; supervision, F.L.; project administration, F.L.; funding acquisition, F.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We are extremely grateful for reviewers’ input in helping this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

LUADLung adenocarcinoma
DEGsDifferentially expressed genes
TNMTumor, node, and metastases
GEOGene Expression Omnibus
TCGAThe Cancer Genome Atlas
PPIProtein-protein interaction
ROCReceiver operating characteristic curve
STRINGSearch Tool for the Retrieval of Interacting Genes database
MCODEMolecularcomplex detection
OSOverall survival
LASSOLeast absolute shrinkage and selection operator
GEPIAGene Expression Profiling Interactive Analysis
HPAHuman Protein Atlas
GOGene Ontology
KEGGKyoto Encyclopedia of Genes and Genomes

References

  1. Ling, B.; Liao, X.; Huang, Y.; Liang, L.; Jiang, Y.; Pang, Y.; Qi, G. Identification of prognostic markers of lung cancer through bioinformatics analysis and in vitro experiments. Int. J. Oncol. 2020, 56, 193–205. [Google Scholar] [CrossRef] [PubMed]
  2. Zhang, L.; Zhang, Z.; Yu, Z. Identification of a novel glycolysis-related gene signature for predicting metastasis and survival in patients with lung adenocarcinoma. J. Transl. Med. 2019, 17, 423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Song, Q.; Shang, J.; Yang, Z.; Zhang, L.; Zhang, C.; Chen, J.; Wu, X. Identification of an immune signature predicting prognosis risk of patients in lung adenocarcinoma. J. Transl. Med. 2019, 17, 70. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Denisenko, T.V.; Budkevich, I.N.; Zhivotovsky, B. Cell death-based treatment of lung adenocarcinoma. Cell Death Dis. 2018, 9, 117. [Google Scholar] [CrossRef] [PubMed]
  5. Niemira, M.; Collin, F.; Szalkowska, A.; Bielska, A.; Chwialkowska, K.; Reszec, J.; Niklinski, J.; Kwasniewski, M.; Kretowski, A. Molecular Signature of Subtypes of Non-Small-Cell Lung Cancer by Large-Scale Transcriptional Profiling: Identification of Key Modules and Genes by Weighted Gene Co-Expression Network Analysis (WGCNA). Cancers 2019, 12, 37. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Barresi, V.; Cinnirella, G.; Valenti, G.; Spampinato, G.; Musso, N.; Castorina, S.; Condorelli, D.F. Gene expression profiles in genome instability-based classes of colorectal cancer. BMC Cancer 2018, 18, 1265. [Google Scholar] [CrossRef]
  7. Yoshimi, A.; Lin, K.T.; Wiseman, D.H.; Rahman, M.A.; Pastore, A.; Wang, B.; Lee, S.C.; Micol, J.B.; Zhang, X.J.; de Botton, S.; et al. Coordinated alterations in RNA splicing and epigenetic regulation drive leukaemogenesis. Nature 2019, 574, 273–277. [Google Scholar] [CrossRef]
  8. Sung, J.Y.; Shin, H.T.; Sohn, K.A.; Shin, S.Y.; Park, W.Y.; Joung, J.G. Assessment of intratumoral heterogeneity with mutations and gene expression profiles. PLoS ONE 2019, 14, e0219682. [Google Scholar] [CrossRef] [Green Version]
  9. Chen, H.; Sheng, L.; Gong, Z.; Ru, S.; Bian, H. Investigation of the molecular mechanisms of hepatic injury upon naphthalene exposure in zebrafish (Danio rerio). Ecotoxicology 2018, 27, 650–660. [Google Scholar] [CrossRef]
  10. Zhou, C.; Ji, J.; Cai, Q.; Shi, M.; Chen, X.; Yu, Y.; Zhu, Z.; Zhang, J. MTA2 enhances colony formation and tumor growth of gastric cancer cells through IL-11. BMC Cancer 2015, 15, 343. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Xie, H.; Xie, C. A Six-Gene Signature Predicts Survival of Adenocarcinoma Type of Non-Small-Cell Lung Cancer Patients: A Comprehensive Study Based on Integrated Analysis and Weighted Gene Coexpression Network. BioMed Res. Int. 2019, 2019, 4250613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Zhu, J.; Wang, M.; Hu, D. Development of an autophagy-related gene prognostic signature in lung adenocarcinoma and lung squamous cell carcinoma. PeerJ 2020, 8, e8288. [Google Scholar] [CrossRef] [PubMed]
  13. Hou, J.; Aerts, J.; den-Hamer, B.; van-Ijcken, W.; den-Bakker, M.; Riegman, P.; van-der-Leest, C.; van-der-Spek, P.; Foekens, J.A.; Hoogsteden, H.C.; et al. Gene expression-based classification of non-small cell lung carcinomas and survival prediction. PLoS ONE 2010, 5, e10312. [Google Scholar] [CrossRef]
  14. Selamat, S.A.; Chung, B.S.; Girard, L.; Zhang, W.; Zhang, Y.; Campan, M.; Siegmund, K.D.; Koss, M.N.; Hagen, J.A.; Lam, W.L.; et al. Genome-scale analysis of DNA methylation in lung adenocarcinoma and integration with mRNA expression. Genome Res. 2012, 22, 1197–1211. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Okayama, H.; Kohno, T.; Ishii, Y.; Shimada, Y.; Shiraishi, K.; Iwakawa, R.; Furuta, K.; Tsuta, K.; Shibata, T.; Yamamoto, S.; et al. Identification of genes upregulated in ALK-positive and EGFR/KRAS/ALK-negative lung adenocarcinomas. Cancer Res. 2012, 72, 100–111. [Google Scholar] [CrossRef] [Green Version]
  16. Landi, M.T.; Dracheva, T.; Rotunno, M.; Figueroa, J.D.; Liu, H.; Dasgupta, A.; Mann, F.E.; Fukuoka, J.; Hames, M.; Bergen, A.W.; et al. Gene expression signature of cigarette smoking and its role in lung adenocarcinoma development and survival. PLoS ONE 2008, 3, e1651. [Google Scholar] [CrossRef]
  17. Stearman, R.S.; Dwyer-Nield, L.; Zerbe, L.; Blaine, S.A.; Chan, Z.; Bunn, P.A.J.; Johnson, G.L.; Hirsch, F.R.; Merrick, D.T.; Franklin, W.A.; et al. Analysis of orthologous gene expression between human pulmonary adenocarcinoma and a carcinogen-induced murine model. Am. J. Pathol. 2005, 167, 1763–1775. [Google Scholar] [CrossRef] [Green Version]
  18. Su, L.J.; Chang, C.W.; Wu, Y.C.; Chen, K.C.; Lin, C.J.; Liang, S.C.; Lin, C.H.; Whang-Peng, J.; Hsu, S.L.; Chen, C.H.; et al. Selection of DDX5 as a novel internal control for Q-RT-PCR from microarray data using a block bootstrap re-sampling scheme. BMC Genom. 2007, 8, 140. [Google Scholar] [CrossRef] [Green Version]
  19. Wu, T.; Wang, M.C.; Jing, L.; Liu, Z.Y.; Guo, H.; Liu, Y.; Bai, Y.Y.; Cheng, Y.Z.; Nan, K.J.; Liang, X. Autophagy facilitates lung adenocarcinoma resistance to cisplatin treatment by activation of AMPK/mTOR signaling pathway. Drug Des. Dev. Ther. 2015, 9, 6421–6431. [Google Scholar] [CrossRef] [Green Version]
  20. Zhang, C.; Zhang, J.; Xu, F.P.; Wang, Y.G.; Xie, Z.; Su, J.; Dong, S.; Nie, Q.; Shao, Y.; Zhou, Q.; et al. Genomic Landscape and Immune Microenvironment Features of Preinvasive and Early Invasive Lung Adenocarcinoma. J. Thorac. Oncol. 2019, 14, 1912–1923. [Google Scholar] [CrossRef]
  21. Guo, D.; Wang, M.; Shen, Z.; Zhu, J. A new immune signature for survival prediction and immune checkpoint molecules in lung adenocarcinoma. J. Transl. Med. 2020, 18, 123. [Google Scholar] [CrossRef] [Green Version]
  22. Li, R.; Yang, Y.E.; Yin, Y.H.; Zhang, M.Y.; Li, H.; Qu, Y.Q. Methylation and transcriptome analysis reveal lung adenocarcinoma-specific diagnostic biomarkers. J. Transl. Med. 2019, 17, 324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Zhang, Y.; Yuan, Y.; Li, Y.; Zhang, P.; Chen, P.; Sun, S. An inverse interaction between HOXA11 and HOXA11-AS is associated with cisplatin resistance in lung adenocarcinoma. Epigenetics 2019, 14, 949–960. [Google Scholar] [CrossRef]
  24. Stutvoet, T.S.; Kol, A.; de Vries, E.G.; de Bruyn, M.; Fehrmann, R.S.; Terwisscha van Scheltinga, A.G.; de Jong, S. MAPK pathway activity plays a key role in PD-L1 expression of lung adenocarcinoma cells. J. Pathol. 2019, 249, 52–64. [Google Scholar] [CrossRef] [Green Version]
  25. Cao, P.; Liu, B.; Du, F.; Li, D.; Wang, Y.; Yan, X.; Li, X.; Li, Y. Scutellarin suppresses proliferation and promotes apoptosis in A549 lung adenocarcinoma cells via AKT/mTOR/4EBP1 and STAT3 pathways. Thorac. Cancer 2019, 10, 492–500. [Google Scholar] [CrossRef]
  26. Zhao, C.; Zhang, J.; Ma, L.; Wu, H.; Zhang, H.; Su, J.; Geng, B.; Yao, Q.; Zheng, J. GOLPH3 Promotes Angiogenesis of Lung Adenocarcinoma by Regulating the Wnt/beta-Catenin Signaling Pathway. OncoTargets Ther. 2020, 13, 6265–6277. [Google Scholar] [CrossRef]
  27. Cui, Y.; Fang, W.; Li, C.; Tang, K.; Zhang, J.; Lei, Y.; He, W.; Peng, S.; Kuang, M.; Zhang, H.; et al. Development and Validation of a Novel Signature to Predict Overall Survival in “Driver Gene-negative” Lung Adenocarcinoma (LUAD): Results of a Multicenter Study. Clin. Cancer Res. 2019, 25, 1546–1556. [Google Scholar] [CrossRef] [Green Version]
  28. Li, B.; Cui, Y.; Diehn, M.; Li, R. Development and Validation of an Individualized Immune Prognostic Signature in Early-Stage Nonsquamous Non-Small Cell Lung Cancer. JAMA Oncol. 2017, 3, 1529–1537. [Google Scholar] [CrossRef]
  29. Li, S.; Xuan, Y.; Gao, B.; Sun, X.; Miao, S.; Lu, T.; Wang, Y.; Jiao, W. Identification of an eight-gene prognostic signature for lung adenocarcinoma. Cancer Manag. Res. 2018, 10, 3383–3392. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, L.; Jiang, B.; Zhu, N.; Tao, M.; Jun, Y.; Chen, X.; Wang, Q.; Luo, C. Mitotic checkpoint kinase Mps1/TTK predicts prognosis of colon cancer patients and regulates tumor proliferation and differentiation via PKCalpha/ERK1/2 and PI3K/Akt pathway. Med. Oncol. 2019, 37, 5. [Google Scholar] [CrossRef] [PubMed]
  31. King, J.L.; Zhang, B.; Li, Y.; Li, K.P.; Ni, J.J.; Saavedra, H.I.; Dong, J.T. TTK promotes mesenchymal signaling via multiple mechanisms in triple negative breast cancer. Oncogenesis 2018, 7, 69. [Google Scholar] [CrossRef] [PubMed]
  32. Ishigami, S.; Ueno, S.; Nishizono, Y.; Matsumoto, M.; Kurahara, H.; Arigami, T.; Uchikado, Y.; Setoyama, T.; Arima, H.; Yoshiaki, K.; et al. Prognostic impact of CD168 expression in gastric cancer. BMC Cancer 2011, 11, 106. [Google Scholar] [CrossRef] [Green Version]
  33. Ye, S.; Liu, Y.; Fuller, A.M.; Katti, R.; Ciotti, G.E.; Chor, S.; Alam, M.Z.; Devalaraja, S.; Lorent, K.; Weber, K.; et al. TGFbeta and Hippo Pathways Cooperate to Enhance Sarcomagenesis and Metastasis through the Hyaluronan-Mediated Motility Receptor (HMMR). Mol. Cancer Res. 2020, 18, 560–573. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Pai, V.C.; Hsu, C.C.; Chan, T.S.; Liao, W.Y.; Chuu, C.P.; Chen, W.Y.; Li, C.R.; Lin, C.Y.; Huang, S.P.; Chen, L.T.; et al. ASPM promotes prostate cancer stemness and progression by augmenting Wnt-Dvl-3-beta-catenin signaling. Oncogene 2019, 38, 1340–1353. [Google Scholar] [CrossRef]
  35. Liu, X.; Liu, X.; Li, J.; Ren, F. Identification and Integrated Analysis of Key Biomarkers for Diagnosis and Prognosis of Non-Small Cell Lung Cancer. Med. Sci. Monit. 2019, 25, 9280–9289. [Google Scholar] [CrossRef]
  36. Yu, D.; Shi, L.; Bu, Y.; Li, W. Cell Division Cycle Associated 8 Is a Key Regulator of Tamoxifen Resistance in Breast Cancer. J. Breast Cancer 2019, 22, 237–247. [Google Scholar] [CrossRef]
  37. Bu, Y.; Shi, L.; Yu, D.; Liang, Z.; Li, W. CDCA8 is a key mediator of estrogen-stimulated cell proliferation in breast cancer cells. Gene 2019, 703, 1–6. [Google Scholar] [CrossRef]
  38. Ci, C.; Tang, B.; Lyu, D.; Liu, W.; Qiang, D.; Ji, X.; Qiu, X.; Chen, L.; Ding, W. Overexpression of CDCA8 promotes the malignant progression of cutaneous melanoma and leads to poor prognosis. Int. J. Mol. Med. 2019, 43, 404–412. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Gnjatic, S.; Cao, Y.; Reichelt, U.; Yekebas, E.F.; Nolker, C.; Marx, A.H.; Erbersdobler, A.; Nishikawa, H.; Hildebrandt, Y.; Bartels, K.; et al. NY-CO-58/KIF2C is overexpressed in a variety of solid tumors and induces frequent T cell responses in patients with colorectal cancer. Int. J. Cancer 2010, 127, 381–393. [Google Scholar] [CrossRef]
  40. Song, X.; Zhang, T.; Wang, X.; Liao, X.; Han, C.; Yang, C.; Su, K.; Cao, W.; Gong, Y.; Chen, Z.; et al. Distinct Diagnostic and Prognostic Values of Kinesin Family Member Genes Expression in Patients with Breast Cancer. Med. Sci. Monit. 2018, 24, 9442–9464. [Google Scholar] [CrossRef]
  41. Gwon, M.R.; Cho, J.H.; Kim, J.R. Mitotic centromere-associated kinase (MCAK/Kif2C) regulates cellular senescence in human primary cells through a p53-dependent pathway. FEBS Lett. 2012, 586, 4148–4156. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Arsic, N.; Bendris, N.; Peter, M.; Begon-Pescia, C.; Rebouissou, C.; Gadea, G.; Bouquier, N.; Bibeau, F.; Lemmers, B.; Blanchard, J.M. A novel function for Cyclin A2: Control of cell invasion via RhoA signaling. J. Cell. Biol. 2012, 196, 147–162. [Google Scholar] [CrossRef] [Green Version]
  43. Bendris, N.; Arsic, N.; Lemmers, B.; Blanchard, J.M. Cyclin A2, Rho GTPases and EMT. Small GTPases 2012, 3, 225–228. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Li, H.; Tian, X.; Wang, P.; Huang, M.; Xu, R.; Nie, T. MicroRNA-582-3p negatively regulates cell proliferation and cell cycle progression in acute myeloid leukemia by targeting cyclin B2. Cell. Mol. Biol. Lett. 2019, 24, 66. [Google Scholar] [CrossRef]
  45. Gao, C.L.; Wang, G.W.; Yang, G.Q.; Yang, H.; Zhuang, L. Karyopherin subunit-alpha 2 expression accelerates cell cycle progression by upregulating CCNB2 and CDK1 in hepatocellular carcinoma. Oncol. Lett. 2018, 15, 2815–2820. [Google Scholar]
  46. Gao, X.; Chen, Y.; Chen, M.; Wang, S.; Wen, X.; Zhang, S. Identification of key candidate genes and biological pathways in bladder cancer. PeerJ 2018, 6, e6036. [Google Scholar] [CrossRef] [PubMed]
Ijms 21 08479 g001 550
Figure 1. Flow chart for the construction of a risk score model in lung cancer.
Figure 1. Flow chart for the construction of a risk score model in lung cancer.
Ijms 21 08479 g001
Ijms 21 08479 g002 550
Figure 2. Identification of differentially expressed genes (DEGs) from four Gene Expression Omnibus (GEO) profiles. (A–D) Upregulated (redcolored spots) and downregulated (greencolored spots) DEGs in lung adenocarcinoma (LUAD) (compared to normal lung tissues) screened from the GEO profile GSE18842 (A), GSE75037 (B), GSE101929(C), and GSE19188 (D). (E–F) A total of 197 upregulated (E) and 376 downregulated (F) DEGs were shared among the four GEO expression profiles.
Figure 2. Identification of differentially expressed genes (DEGs) from four Gene Expression Omnibus (GEO) profiles. (A–D) Upregulated (redcolored spots) and downregulated (greencolored spots) DEGs in lung adenocarcinoma (LUAD) (compared to normal lung tissues) screened from the GEO profile GSE18842 (A), GSE75037 (B), GSE101929(C), and GSE19188 (D). (E–F) A total of 197 upregulated (E) and 376 downregulated (F) DEGs were shared among the four GEO expression profiles.
Ijms 21 08479 g002
Ijms 21 08479 g003 550
Figure 3. Protein-protein interaction (PPI) network construction and enrichment analysis for modules. (A) PPI network and module clustering analysis for the 573 shared DEGs. Blue, green, pink, and orange nodes represent the four clustering modules, of which cluster 1 (blue) had the highest MCODE score. (B–E) Biological pathway analysis for cluster 1 (B), cluster 2 (C), cluster 3 (D), and cluster 4 (E). Significant signaling pathways were mainly involved in cell cycle and DNA Replication.
Figure 3. Protein-protein interaction (PPI) network construction and enrichment analysis for modules. (A) PPI network and module clustering analysis for the 573 shared DEGs. Blue, green, pink, and orange nodes represent the four clustering modules, of which cluster 1 (blue) had the highest MCODE score. (B–E) Biological pathway analysis for cluster 1 (B), cluster 2 (C), cluster 3 (D), and cluster 4 (E). Significant signaling pathways were mainly involved in cell cycle and DNA Replication.
Ijms 21 08479 g003
Ijms 21 08479 g004 550
Figure 4. Development of gene signature in the training group. (A) Survival curves of high- and low-risk groups separated by gene signature. (B) Receiver operating characteristic (ROC) curves for survival risk predicted by gene signature for 1-, 3-, and 5-year follow-ups. (C–E) Expression heatmap of eight signature genes (C), risk score distribution (D), and survival status of patients (E). The risk score based on gene signature appears to be correlated with survival status for lung cancer patients in low- and high-risk groups.
Figure 4. Development of gene signature in the training group. (A) Survival curves of high- and low-risk groups separated by gene signature. (B) Receiver operating characteristic (ROC) curves for survival risk predicted by gene signature for 1-, 3-, and 5-year follow-ups. (C–E) Expression heatmap of eight signature genes (C), risk score distribution (D), and survival status of patients (E). The risk score based on gene signature appears to be correlated with survival status for lung cancer patients in low- and high-risk groups.
Ijms 21 08479 g004
Ijms 21 08479 g005 550
Figure 5. Validation of gene signature in the testing group. (A) Survival curve of high- and low-risk groups separated by gene signature. (B) ROC analysis for survival rate predicted by gene signaturefor 1-, 3-, and 5-year follow-ups. (C–E) Expression heatmap of eight signature genes (C), risk score distribution (D), and survival status of patients (E). The risk score based on gene signature in the testing group also appears to be correlated with survival status for lung cancer patients in low- and high-risk groups, which was consistent with that in training group.
Figure 5. Validation of gene signature in the testing group. (A) Survival curve of high- and low-risk groups separated by gene signature. (B) ROC analysis for survival rate predicted by gene signaturefor 1-, 3-, and 5-year follow-ups. (C–E) Expression heatmap of eight signature genes (C), risk score distribution (D), and survival status of patients (E). The risk score based on gene signature in the testing group also appears to be correlated with survival status for lung cancer patients in low- and high-risk groups, which was consistent with that in training group.
Ijms 21 08479 g005
Ijms 21 08479 g006 550
Figure 6. The prognostic value of riskscore signature in LUAD. (A–J) Kaplan–Meier curves for OS in patients with LUAD; clinical characteristics include age: ≤65 (A) and >65 (B), gender: female (C) and male (D), clinical stage I/II (E) and III/IV (F), lymph node metastasis: N0 (G) and N1-3 (H), and primary tumor: T1-2 (I) and T3-4 (J). Male, >65 age, stage III/IV, N1-2 and T1-2 subgroups were significantly associated with worse OS in high-risk group with respective p < 0.05.
Figure 6. The prognostic value of riskscore signature in LUAD. (A–J) Kaplan–Meier curves for OS in patients with LUAD; clinical characteristics include age: ≤65 (A) and >65 (B), gender: female (C) and male (D), clinical stage I/II (E) and III/IV (F), lymph node metastasis: N0 (G) and N1-3 (H), and primary tumor: T1-2 (I) and T3-4 (J). Male, >65 age, stage III/IV, N1-2 and T1-2 subgroups were significantly associated with worse OS in high-risk group with respective p < 0.05.
Ijms 21 08479 g006
Ijms 21 08479 g007 550
Figure 7. Validations of eight genes contained in the risk score model. (A) Scatter plots between mRNA expression of each gene (y-axis) and corresponding risk score (x-axis) in patients with LUAD for these eight genes, which showed significant positive Pearson correlation coefficients. (B) Comparisons of mRNA expression levels of each gene in LUAD tissues versus normal lung tissues, respectively. All genes had higher expression levels in cancer tissues than in normal lung tissues. The red and gray boxes represent cancer and normal tissues, respectively. The red-marked asterisk indicates that the differential expression of each mRNA is significant (p < 0.05).
Figure 7. Validations of eight genes contained in the risk score model. (A) Scatter plots between mRNA expression of each gene (y-axis) and corresponding risk score (x-axis) in patients with LUAD for these eight genes, which showed significant positive Pearson correlation coefficients. (B) Comparisons of mRNA expression levels of each gene in LUAD tissues versus normal lung tissues, respectively. All genes had higher expression levels in cancer tissues than in normal lung tissues. The red and gray boxes represent cancer and normal tissues, respectively. The red-marked asterisk indicates that the differential expression of each mRNA is significant (p < 0.05).
Ijms 21 08479 g007
Ijms 21 08479 g008 550
Figure 8. Expression profiles of eight signature genes in the HBE, A549, and PC-9 cell lines. (A) RT-PCR analysis of relative mRNA expression. The genes expression levels for each of the eight genes were significantly increased in A549 and PC-9 cells compared to HBE cells (p < 0.0001). (B) Western blot analysis of protein expression levels. Every protein appeared to have higher expression levels in A549 and PC-9 cells compared with HBE cells. The measurements of mRNA and protein expression levels were performed in triplicate.
Figure 8. Expression profiles of eight signature genes in the HBE, A549, and PC-9 cell lines. (A) RT-PCR analysis of relative mRNA expression. The genes expression levels for each of the eight genes were significantly increased in A549 and PC-9 cells compared to HBE cells (p < 0.0001). (B) Western blot analysis of protein expression levels. Every protein appeared to have higher expression levels in A549 and PC-9 cells compared with HBE cells. The measurements of mRNA and protein expression levels were performed in triplicate.
Ijms 21 08479 g008
Table
Table 1. Clinical characteristics of patients with lung cancer in different datasets.
Table 1. Clinical characteristics of patients with lung cancer in different datasets.
CharacteristicTCGA Data (n, %)GSE18842 (n, %)GSE75037 (n, %)GSE101929 (n, %)GSE19188 (n, %)
PlatformIllumina HiSeq2000 RNA sequencing platformAffymetrix Human Genome U133 Plus 2.0 ArrayIllumina HumanWG-6 v3.0 expression beadchipAffymetrix Human Genome U133 Plus 2.0 ArrayAffymetrix Human Genome U133 Plus 2.0 Array
Samples551 (100.0%)91 (100.0%)166 (100.0%)66 (100.0%)156 (100.0%)
Normal54 (9.8%)45 (49.5%)83 (50.0%)34 (51.5%)65 (41.7%)
Tumor497 (90.2%)46 (50.5%)83 (50.0%)32 (48.5%)91 (58.3%)
Survival Status486 (88.2%)NANA66 (100.0%)82 (52.6%)
Death162 (29.4%)NANA40 (60.6%)50 (32.1%)
Survival324 (58.8%)NANA26 (39.4%)32 (20.5%)
Age467 (84.8%)NA166 (100.0%)66 (100.0%)NA
<=65227 (41.2%)NA58 (34.9%)43 (65.2%)NA
>65240 (43.6%)NA108 (65.1%)23 (34.8%)NA
Gender486 (88.2%)NA166 (100.0%)66 (100.0%)134 (85.9%)
Female264 (47.9%)NA118 (71.1%)38 (57.6%)34 (21.8%)
Male222 (40.3%)NA48 (28.9%)28 (42.4%)100 (64.1%)
Stage478 (86.8%)NA33 (19.9%)NANA
I262 (47.5%)NA20 (12.0%)NANA
II112 (20.3%)NA8 (4.8%)NANA
III79 (14.3%)NA5 (3.0%)NANA
IV25 (4.5%)NANANANA
T classification483 (87.7%)NANANANA
T1163 (29.6%)NANANANA
T2260 (47.2%)NANANANA
T341 (7.4%)NANANANA
T419 (3.4%)NANANANA
N classification474 (86.0%)NANANANA
N0312 (56.6%)NANANANA
N190 (16.3%)NANANANA
N270 (12.7%)NANANANA
N32 (0.4%)NANANANA
M classification357 (64.8%)NANANANA
M0333 (60.4%)NANANANA
M124 (4.4%)NANANANA
For TNM classification, T, N, and M refer to primary tumor, regional lymph nodes, and distant metastasis, respectively. Abbreviations: TCGA, The Cancer Genome Atlas; NA, not available.
Table
Table 2. Top 30 genes with highest degrees of connectivity in clustering module 1.
Table 2. Top 30 genes with highest degrees of connectivity in clustering module 1.
NO.GeneDegreeNO.GeneDegreeNO.GeneDegree
1UBE2C7111ASPM7121KIAA010170
2NUSAP17112CENPF7122SPAG570
3TPX27113CDCA87123KIF1570
4PBK7114KIF2C7124CEP5570
5MELK7115AURKB7125CENPE70
6TTK7116CCNB27126CDC2070
7KIF117117KIF20A7127BIRC570
8TOP2A7118MKI677128MCM1070
9HMMR7119CCNA27129MAD2L170
10RRM27120CCNB17130AURKA70
Table
Table 3. Identification of gene expression signature by univariate and multivariate Cox regression analysis.
Table 3. Identification of gene expression signature by univariate and multivariate Cox regression analysis.
NO.GeneUnivariate Analysis *Multivariate Analysis **
HR95%CIpHR95%CICoef.
1UBE2C1.1451.033–1.2700.010 --- --- ---
2TPX21.2261.089–1.3810.001 --- --- ---
3PBK1.2641.100–1.4530.001 --- --- ---
4MELK1.2331.069–1.4220.004 --- --- ---
5TTK1.2471.053–1.4770.0100.6300.341–1.165−0.462
6KIF111.3581.148–1.608<0.001 --- --- ---
7TOP2A1.1781.043–1.3310.008 --- --- ---
8HMMR1.4721.243–1.742<0.0011.8831.153–3.0740.633
9RRM21.2981.128–1.493<0.001 --- --- ---
10ASPM1.4091.169–1.698<0.0010.5770.287–1.159−0.550
11CENPF1.2931.112–1.5030.001 --- --- ---
12CDCA81.2061.039–1.4010.0140.2700.100–0.730−1.309
13KIF2C1.2341.074–1.4170.0033.2811.232–8.7381.188
14AURKB1.1881.039–1.3580.012 --- --- ---
15CCNB21.2581.085–1.4580.0020.6220.329–1.178−0.474
16KIF20A1.3501.136–1.6050.001 --- --- ---
17MKI671.3091.129–1.518<0.0011.7681.103–2.8350.570
18CCNA21.3281.150–1.533<0.0011.6220.889–2.9590.484
19CCNB11.3211.136–1.535<0.001 --- --- ---
20NUSAP11.2931.107–1.5110.001 --- --- ---
* The 20 DEGs were significantly associated with overall survival (p < 0.05) using univariate Cox regression analysis. ** Then, a least absolute shrinkage and selection operator (LASSO) regression on these 20 DEGs was performed to identify the most informative gene set for survival prediction. Finally, eight genes marked in gray in the table were selected to perform multivariate Cox regression analysis and to generate a prognostic risk model according to their respective regression coefficients. Abbreviations: HR, hazard ratio; CI, confidence interval; p, p-value; Coef., coefficient.
Table
Table 4. Univariate and multivariate Cox regression analysis for risk score on overall survival of patients with LUAD.
Table 4. Univariate and multivariate Cox regression analysis for risk score on overall survival of patients with LUAD.
ParameterUnivariate AnalysisMultivariate Analysis
HR95%CIpHR95%CIp
Training Group
Age0.4920.169–1.4310.1930.6120.204–1.8380.381
Gender1.1280.656–1.9400.6631.0620.605–1.8650.833
Stage0.8080.421–1.5540.5230.2950.061–1.4120.126
T classification1.7530.826–3.7210.1442.7551.086–6.9880.033
M classification0.9610.233–3.970.9563.1620.328–30.5010.320
N classification0.9140.578–1.4440.6991.5340.637–3.6900.340
RiskScore3.2851.681–6.4200.0012.9311.474–5.8290.002
Testing group
Age1.0080.978–1.0380.6120.9910.962–1.0210.556
Gender0.6230.352–1.1050.1060.5950.327–1.0830.089
Stage1.0390.792–1.3630.7820.8610.371–2.0010.728
T classification1.0630.755–1.4960.7261.0830.704–1.6640.717
M classification0.9430.372–2.3930.9020.9160.142–5.9090.926
N classification1.1780.802–1.7300.4041.6280.751–3.5290.217
RiskScore1.5941.256–2.022<0.0011.6621.284–2.152<0.001
For TNM classification, T, N, and M refer to primary tumor, regional lymph nodes, and distant metastasis, respectively. Abbreviations: LUAD, lung adenocarcinoma; HR, hazard ratio; CI, confidence interval; p, p-value.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Z.; Qi, F.; Li, F. Establishment of a Gene Signature to Predict Prognosis for Patients with Lung Adenocarcinoma. Int. J. Mol. Sci. 2020, 21, 8479. https://doi.org/10.3390/ijms21228479

AMA Style

Li Z, Qi F, Li F. Establishment of a Gene Signature to Predict Prognosis for Patients with Lung Adenocarcinoma. International Journal of Molecular Sciences. 2020; 21(22):8479. https://doi.org/10.3390/ijms21228479

Chicago/Turabian Style

Li, Zhaodong, Fangyuan Qi, and Fan Li. 2020. "Establishment of a Gene Signature to Predict Prognosis for Patients with Lung Adenocarcinoma" International Journal of Molecular Sciences 21, no. 22: 8479. https://doi.org/10.3390/ijms21228479

APA Style

Li, Z., Qi, F., & Li, F. (2020). Establishment of a Gene Signature to Predict Prognosis for Patients with Lung Adenocarcinoma. International Journal of Molecular Sciences, 21(22), 8479. https://doi.org/10.3390/ijms21228479

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop