Research Paper Volume 10, Issue 9 pp 2356—2366

A seven-long noncoding RNA signature predicts overall survival for patients with early stage non-small cell lung cancer

Ting Lin1,2, *,, Yunong Fu1, *,, Xing Zhang1, *,, Jingxian Gu1, , Xiaohua Ma1, , Runchen Miao1, , Xiaohong Xiang1, , Wenquan Niu3, , Kai Qu1, , Chang Liu1, , Qifei Wu4, ,

  • 1 Department of Hepatobiliary Surgery, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’a, Shaanxi 710061, China
  • 2 Department of Surgical Intensive Care Units, The First Affiliated Hospital of Xi'an Jiaotong University, Xi’a, Shaanxi 710061, China
  • 3 Institute of Clinical Medical Sciences, China-Japan Friendship Hospital, Beijing 100029, China
  • 4 Department of Thoracic Surgery, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’a, Shaanxi 710061, China
* Equal contribution

Received: July 5, 2018       Accepted: September 6, 2018       Published: September 11, 2018      

https://doi.org/10.18632/aging.101550
How to Cite

Copyright: Lin et al. This is an open‐access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Non-small cell lung cancer (NSCLC) is the most common cancer and cause of cancer-related mortality globally. Increasing evidence suggested that the long non-coding RNAs (lncRNAs) were involved in cancer-related death. To explore the possible prognostic lncRNA biomarkers for NSCLC patients, in the present study, we conducted a comprehensive lncRNA profiling analysis based on 1902 patients from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) datasets. In the discovery phase, we employed 682 patients from the combination of four GEO datasets (GSE30219, GSE31546, GSE33745 and GSE50081) and conducted a seven-lncRNA formula to predict overall survival (OS). Next, we validated our risk-score formula in two independent datasets, TCGA (n=994) and GSE31210 (n=226). Stratified analysis revealed that the seven-lncRNA signature was significantly associated with OS in stage I patients from both discovery and validation groups (all P<0.001). Additionally, the prognostic value of the seven-lncRNA signature was also found to be favorable in patients carrying wild-type KRAS or EGFR. Bioinformatical analysis suggested that the seven-lncRNA signature affected patients’ prognosis by influencing cell cycle-related pathways. In summary, our findings revealed a seven-lncRNA signature that predicted OS of NSCLC patients, especially in those with early tumor stage and carrying wild-type KRAS or EGFR.

Introduction

Non-small cell lung cancer (NSCLC) is one of the most common and lethal malignant diseases worldwide during the past decades due to lacking of early diagnostic and predictive biomarkers [1,2]. With the development of the molecular targeted agents [3] and the therapeutic strategy [4], the 5-year overall survival (OS) rate of NSCLC has been prolonged. However, only one third of NSCLC patients are diagnosed at an early stage [5]. To date, there are limited methods for us to provide a prognostic prediction of NSCLC patients, which is important to choose appropriate therapies for those patients [6]. Therefore, our study focused on finding a prognostic risk-score model derived from the whole-genome data.

Long non-coding RNAs (lncRNAs) are defined as non-protein coding transcripts longer than 200 nucleotides [7,8]. During the past decade, genome-wide sequencing helped researchers discover a large amount of lncRNAs [8]. Increasing studies suggested that the lncRNAs were involved in multiple cellular biological process and were associated with different human diseases including cancer [9,10]. Furthermore, several lncRNAs have been identified as oncogenes, such as MALAT1 in lung cancer [11], HOTAIR in breast cancer [12], and PANDA in hepatocellular carcinoma [13]. Above lncRNAs were also found to be associated with the prognosis of cancer patients [1315].

Recently, the prognostic values of lncRNAs have attracted much attention in the field of cancer research. Researchers put great effort into investigating the lncRNA biomarkers in multiple malignancies, including gastric cancer [16], breast cancer [17], hepatocellular carcinoma [18] and ovarian cancer [19]. In lung cancer, Tu et al explored the lncRNA biomarkers in based on 739 patients derived from GEO datasets [20] and identified an eight-lncRNA signature which was associated with clinical outcomes. In the present study, we conducted a large sample size (n=1902) and cross-platform (including RNA-seq and microarray data) analysis, and performed multiple validation. Considering the above clear advantages, our findings will provide robust evidence for the exploration of lncRNA biomarkers in NSCLC.

Results

Identification of seven lncRNAs as prognostic biomarkers in discovery group

By subjecting the lncRNA profiling data of discovery group to univariable Cox proportional hazards regression analysis, we obtained 457 lncRNAs whose P-value were less than 0.05. Among them, seven lncRNAs were selected for further analysis, including APTR, DHRS4-AS1, ITGA9-AS1, LINC01137, LOC101927972, RPARP-AS1 and SH3BP5-AS1. Next, we fitted the 7 lncRNAs expression, the corresponding survival status and survival time into the multivariable Cox regression. Using the coefficients obtained from the multivariable Cox regression, a risk-score formula was constructed as following: risk score = - 0.06 x APTR - 0.20 x DHRS4-AS1 - 0.35 x ITGA9-AS1 - 0.23 x LINC01137 - 0.25 x LOC101927972 - 0.06 x RPARP-AS1 - 0.20 x SH3BP5-AS1. Notably, in the formula, every lncRNA expression value is weight by a negative coefficient (Table 1). As shown in Figure 1A and 1B, the heatmap revealed that the expression level of seven lncRNAs was decreased accompanying with the higher risk scores. Furthermore, we calculated the association between risk score and cancer-related death (Figure 1C). Our data showed that mortality rate in high risk group was significant higher than low risk group (Figure 1D), indicating the seven lncRNAs might play a protective role in NSCLC.

Table 1. Seven lncRNAs significantly associated with the OS of NSCLC patients in the training group.

Gene nameEnsemble IDChr.CoordinateZ scoreP value aCoefficient b
APTRENSG00000214293777,657,660-77,697,345-3.251.16E-03-0.06
DHRS4-AS1ENSG000002152561423,938,731-23,988,839-3.141.69E-03-0.20
ITGA9-AS1ENSG00000235257337,745,432-37,861,780-2.834.62E-03-0.35
LINC01137ENSG00000233621137,454,879-37,474,4433.261.12E-03-0.23
LOC101927972ENSG0000026998565,451,683-5,458,075-2.982.89E-03-0.25
RPARP-AS1ENSG0000026960910102,449,817-102,461,106-3.358.05E-04-0.06
SH3BP5-AS1ENSG00000224660315,254,184-15,264,498-2.923.45E-03-0.20
a Derived from the univariable Cox proportional hazards regression analysis in the discovery group.
b Derived from the multivariable Cox proportional hazards regression analysis in the discovery group.
The seven-lncRNA signature-based risk score in prognosis of overall survival in the discovery group. (A) The seven lncRNA-based risk score distribution. (B) The heatmap of the seven lncRNA expression profiles. (C) The vital status of patients in high- and low-risk groups. (D) The mortality rate in low- and high-risk score groups.

Figure 1. The seven-lncRNA signature-based risk score in prognosis of overall survival in the discovery group. (A) The seven lncRNA-based risk score distribution. (B) The heatmap of the seven lncRNA expression profiles. (C) The vital status of patients in high- and low-risk groups. (D) The mortality rate in low- and high-risk score groups.

The prognostic values of seven-lncRNA signature in discovery and validation groups

Considering the relationship between the seven-lncRNA signature and patients’ deaths, we next explored whether the seven-lncRNA signature could affect patients’ OS time. Our data in discovery group (Meta-GEO dataset, n=682) showed that the patients who had low risk-scores were supposed to have a longer OS time than higher risk group (Figure 2A). To validate above finding, we employed two validation groups, TCGA and GSE31210 dataset, containing 994 and 226 patients, respectively. As expected, patients in high risk group had a significant increased mortality risk than low risk group either in TCGA (HR=1.34, 95% CI= 1.10-1.63, P =0.004) or GSE31210 datasets (HR=3.06, 95% CI= 1.58-5.93, P =0.002) (Figure 2B and 2C). Moreover, prognostic meta-analysis based on above three groups (n=1902) confirmed that the seven-lncRNA signature was a risk factor for NSCLC patients (combined HR=1.76, 95% CI= 1.26-2.47, P=0.001) (Figure 2D).

The association between seven-lncRNA signature and overall survival in discovery and two validation groups. Kaplan-Meier survival curves were plotted to estimate the overall survival probabilities for the low- versus high-risk group in the discovery group (A), validation group-1 (B) and validation group-2 (C). (D) Meta-analysis was performed using the prognostic results of three groups.

Figure 2. The association between seven-lncRNA signature and overall survival in discovery and two validation groups. Kaplan-Meier survival curves were plotted to estimate the overall survival probabilities for the low- versus high-risk group in the discovery group (A), validation group-1 (B) and validation group-2 (C). (D) Meta-analysis was performed using the prognostic results of three groups.

The seven-lncRNA signature was associated with OS in stage I patients

To explore the impacts of clinical characteristics on the prognostic values of the seven-lncRNA signature, we performed a set of predefined stratified analysis. As shown in Table 2, we stratified the NSCLC patients by five clinical characteristics (including age, gender, smoking status, pathological subtypes and AJCC stage). According to the results of stratified analysis, the most significant result was found in stage I patients in both discovery and two validation groups (all P<0.001) (Table 2). To validate above findings, we further divided stage I patients into Ia and Ib according AJCC system, and calculated the prognostic values of the seven-lncRNA signature in each subgroups, respectively. Kaplan-Meier method was used to visualize the OS probabilities between high- and low-risk groups. As expected, the overall survival time of high risk group was all significantly shorter than low risk group in either stage Ia or Ib patients (Figure 3). The prognostic meta-analysis showed that the stage Ia patients with high risk scores had a two-fold increased mortality risk than those with low risk scores (combined HR=2.17, 95% CI= 1.15-4.09, P=0.017) (Figure 3D). Similarly, the association between the seven-lncRNA signature and OS also reached a statistical significance in stage Ib patients (combined HR=1.88, 95% CI= 1.30-2.72, P=0.001) (Figure 3H).

Table 2. The association between seven-lncRNA signature and OS of NSCLC patients in discovery and validating groups.

VariableDiscovery GroupValidation Group-1Validation Group-2
Number
(High/Low)
HR (95%CI)P valueNumber
(High/Low)
HR (95%CI)P valueNumber
(High/Low)
HR (95%CI)P value
Total341/3411.85 (1.53-2.25)<0.001497/4971.34 (1.10-1.63)0.004113/1133.06 (1.45-5.44)0.002
Age
   ≤70214/2251.73 (1.36-2.16)<0.001291/2841.30 (0.99-1.70)0.056111/1103.18 (1.45-5.79)0.003
   >70111/852.16 (1.61-3.28)<0.001189/2061.39 (1.04-1.87)0.0292/31.24 (0.12-12.79)0.855
Gender
   Male203/2521.72 (1.36-2.12)<0.001281/3151.37 (1.06-1.75)0.01547/585.04 (1.52-9.20)0.004
   Female138/891.95 (1.38-3.02)<0.001216/1821.25 (0.91-1.75)0.17066/551.69 (0.64-4.52)0.298
Smoking history
   Never smoker23/52.31 (0.41-26.27)0.27557/322.07(1.10-5.46)0.03262/532.14(0.77-5.80)0.152
   Ever smoker66/251.66 (0.80-3.83)0.163105/1471.18(0.81-1.71)0.39851/603.84(1.33-7.68)0.009
   Current smoker44/131.68 (0.75-4.54)0.193319/2991.35(1.04-1.74)0.0230/0NANA
Histology
   Squamous Carcinoma55/1361.32 (0.90-1.88)0.170185/3091.41 (1.05-1.83)0.020NANANA
   Adenocarcinoma199/1191.71 (1.30-2.44)<0.001312/1881.26 (0.94-1.73)0.120NANANA
AJCC stage
Stage I254/2031.70 (1.35-2.23)< 0.001275/2351.68 (1.25-2.31)< 0.00191/779.32 (2.27-15.35)< 0.001
   Stage II56/701.97 (1.24-3.02)0.004119/1581.09 (0.75-1.57)0.64922/361.03 (0.40-2.66)0.944
   Stage III24/591.13 (0.71-1.82)0.60874/891.05 (0.69-1.59)0.8180/0NANA
   Stage IV6/73.09 (1.25-14.41)0.03021/111.25 (0.51-3.17)0.6070/0NANA
Abbreviations: HR, Hazard ratio; 95%CI, 95% confidence interval; AJCC, the American Joint Committee on Cancer.
The association between seven-lncRNA signature and overall survival in patients with stage I. The Kaplan-Meier survival curves of discovery group (A), validation group-1 (B) and validation group-2 (C) were plotted, and meta-analysis (D) was conducted in patients with stage Ia. The similar results were obtained from the patients with stage Ib in discovery group (E), validation group-1 (F) and validation group-2 (G), and the prognostic meta-analysis (H) was also conducted.

Figure 3. The association between seven-lncRNA signature and overall survival in patients with stage I. The Kaplan-Meier survival curves of discovery group (A), validation group-1 (B) and validation group-2 (C) were plotted, and meta-analysis (D) was conducted in patients with stage Ia. The similar results were obtained from the patients with stage Ib in discovery group (E), validation group-1 (F) and validation group-2 (G), and the prognostic meta-analysis (H) was also conducted.

The seven-lncRNA signature was associated with OS in patients carrying wild-type KRAS or EGFR genes

EGFR and KRAS were the most frequent mutant genes and associated with poor prognosis in NSCLC. Previous studies suggested that the mutation frequency of above two genes was significantly decreased in patients with early stage. Bearing this in mind, we also performed stratified analysis based on EGFR or KRAS mutation status. Our data showed that the higher risk score was associated with higher mortality risk in wild-type KRAS or EGFR subgroups (Figure 4). Above results were consistent in the two validation groups, suggesting the seven-lncRNA signature acted as a risk factor for patients carrying wild-type KRAS or EGFR genes.

Kaplan-Meier estimates of the overall survival of patients carrying wild-type KRAS or EGFR genes. (A) Kaplan-Meier survival curves were plotted to estimate the overall survival for patients carrying wild-type KRAS gene in two validation groups. (B) Kaplan-Meier survival curves were plotted to estimate the overall survival for patients carrying wild-type EGFR gene in two validation groups.

Figure 4. Kaplan-Meier estimates of the overall survival of patients carrying wild-type KRAS or EGFR genes. (A) Kaplan-Meier survival curves were plotted to estimate the overall survival for patients carrying wild-type KRAS gene in two validation groups. (B) Kaplan-Meier survival curves were plotted to estimate the overall survival for patients carrying wild-type EGFR gene in two validation groups.

Cell cycle-related pathways were associated with the seven-lncRNA signature

Finally, to explore the potential mechanisms of the seven-lncRNA signature, we used the WGCNA method to cluster genes that highly correlated with the risk scores based on the profiling data of TCGA dataset [21]. We identified a total of 11 modules and found that pink module was the only positively correlated with the risk-score (P =2x10-105) (Figure 5A and 5B). Pathway enrichment analysis was then performed on an online-based web tool “Metascape” (http://metascape.org/) using the genes in pink module. As shown in Figure 5C, genes were significantly enriched in cell cycle-related pathways, including “cell cycle”, “cell cycle phase transition” and “regulation of cell cycle process” pathways, suggesting the activation of cell cycle-related pathways might contribute to higher mortality risk in patients with high risk scores.

WGCNA predicted biological pathways associated with the seven-lncRNA signature. (A) The gene clusters obtained by WGCNA method. (B) The relationship between pink module membership and gene significance for risk score. (C) Significantly enriched pathways of the co-expressed genes in pink module.

Figure 5. WGCNA predicted biological pathways associated with the seven-lncRNA signature. (A) The gene clusters obtained by WGCNA method. (B) The relationship between pink module membership and gene significance for risk score. (C) Significantly enriched pathways of the co-expressed genes in pink module.

Discussion

In the past decade, lncRNAs have been demonstrated by accumulating evidence to contribute in carcinogenesis and tumor progression [22]. In the present study, we conducted a seven-lncRNA signature, and validated their association with OS in two independent datasets. Moreover, we found the prognostic values of the seven-lncRNA signature reached a high statistical significance in patients with stage I and carrying wild-type KRAS and EGFR genes. To the best of our knowledge, this is the first study, which was based on the high-throughput lncRNA profiling data from nearly two thousand patients, to explore the lncRNAs as prognostic biomarkers for early stage NSCLC patients. Our findings might provide convincing evidence for clinical application.

In the present study, we identified a total seven lncRNAs which were associated with OS in NSCLC, including APTR, DHRS4-AS1, ITGA9-AS1, LINC01137, LOC101927972, RPARP-AS1 and SH3BP5-AS1. Most of these identified lncRNAs were reported in cancer for the first time, except APTR (Alu-mediated p21 transcriptional regulator). Previous studies reported that APTR could promote the cell proliferation and contribute to carcinogenesis by repressing p21 promoter activities in human glioblastomas [21]. By contrast, in this study, we found that the coefficient of APTR was minus, indicating that APTR might play a protective role in NSCLC. The detailed molecular mechanism warrant further investigations.

The most important finding in this study was that the seven-lncRNA signature was significantly associated with OS in stage I patients. Our data showed their association reached a high statistical significance (P<0.001) either in discovery group or in validation groups, indicating the seven-lncRNA signature can provide a powerful prognostic tool for NSCLC with early stage. Recently, the oncogenic roles of EGFR and KRAS mutations in NSCLC have attracted increasing attentions. Patients carrying mutant EGFR or KRAS always had a poorer prognosis than those carrying wild-type genes [23,24]. Thus, the frequencies of wild-type EGFR and KRAS were always high in patients with early stage. We reviewed the EGFR and KRAS mutation status in TCGA-NSCLC dataset, and found that the mutation frequencies of EGFR and KRAS in stage 1 patients were only 10.2% and 5.7%, respectively, which were lower than advanced stage according to literature reports [25,26]. Interestingly, our data also showed that the seven-lncRNA signature was significantly associated with OS in patients carrying wild-type EGFR or KRAS genes, which might partially explain the favorable prognostic values of the seven-lncRNA signature in early stage patients.

In order to investigate the potential mechanisms affected by the seven lncRNAs, we also performed bioinformatic analysis. By using the WGCNA method, we clustered 11 gene modules from more than 5000 differentially expressed genes, and found that only one of them (pink) was positively correlated with seven-lncRNA signature. Pathway enrichment analysis further suggested that the genes in pink module were mostly enriched in cell cycle-related signaling pathways, indicating the seven-lncRNA signature might affect cell cycle-related pathways and consequently contributed to tumor progression.

In summary, we constructed a risk-score model derived from seven lncRNAs to predict the OS time of NSCLC patients. The risk association of the seven-lncRNA signature with survival was more evident in patients with early stage and carrying wild-type EGFR or KRAS genes. We warrant further studies, especially the large, well-performed cohorts to confirm or refuse our findings.

Materials and Methods

Sample sources and study design

The raw data of gene expression and corresponding clinical information of NSCLC patients were downloaded from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) websites. Four GEO datasets (GSE30219, GSE31546, GSE37745 and GSE50081) were combined and used as the discovery group. After removal the five samples without survival data or enough clinical information, 682 patients from four GEO datasets were finally used for further analysis. Meanwhile, 994 patients from TCGA dataset were employed as validation group. Moreover, we also employed 226 patients from another independent GEO dataset (GSE31210) as the second validation group. All patients in GSE31210 belonged to the stage I and II of the American Joint Committee on Cancer (AJCC) system. The RNA expression profiling data of five GEO datasets were all performed on Affymetrix U133 Plus 2.0 microarray platform and the TCGA data was performed on Illumina sequencing platform. The clinical characteristics of discovery and validation groups were shown in Table 3.

Table 3. Clinical features of NSCLC patients in the training and validating groups.

FeaturesDiscovery group
(n=682)
Validation group-1
(n=994)
Validation group-2
(n=226)
Age (years), no (%)
  ≤70469 (71.7)575 (58.7)221 (97.8)
  >70196 (28.3)404 (41.3)5 (2.2)
Gender, no (%)
  Male227 (33.3)596 (60.0)105 (46.5)
  Female455 (66.7)398 (40.0)121 (53.5)
Smoking status, no (%)
  Never smoker28 (15.9)89 (9.3)115 (50.9)
  Ever smoker
  Current smoker
91 (51.7)
57 (32.4)
252 (26.3)
618 (64.4)
111 (49.1)
0 (0)
Pathological grade, no (%)
  Squamous Carcinoma191 (37.5)494 (49.7)NA
  Adenocarcinoma318 (62.5)500 (50.3)NA
AJCC stage, no (%)
  I457 (67.3)510 (51.9)168 (74.3)
  II126 (18.6)277 (28.2)58 (25.7)
  III83 (12.2)163 (16.6)0 (0)
  IV13 (1.9)32 (3.3)0 (0)
KRAS
Mutant type
Wild type
EGFR
Mutant type
Wild type
NA
NA
NA
NA
24 (31.2)
53 (68.8)
101 (17.9)
463 (82.1)
20 (8.8)
206 (91.2)
127 (56.2)
99 (43.8)
Death status (%)
  Yes418 (61.3)394 (39.6)9 (8.0)
  No264 (38.7)600 (60.4)104 (92.0)
Recurrence status (%)
  Yes222 (39.4)221 (22.2)64 (28.3)
  No341 (60.6)773 (77.8)162 (71.7)
Abbreviations: AJCC, the American Joint Committee on Cancer; KRAS, Kirsten rat sarcoma viral oncogene homolog; EGFR, epidermal growth factor receptor.

Normalization and lncRNA annotation of GEO data and TCGA data

The raw data of five GEO datasets (GSE30219, GSE31210, GSE31546, GSE37745 and GSE50081) were downloaded as probe-level CEL files. Then, all raw data were quantile normalized using Robust Multi-array Average (RMA) method. With the lncRNA-specific probes presented on Affymetrix U133 Plus 2.0 downloaded from Affymetrix website (http://www.affymetrix.com) [15], we finally identified 2986 lncRNA transcripts with RefSeq transcript IDs. As for TCGA data, we downloaded the sequencing data with the type of fragments per kilobase of exon per million fragments (FPKM) from the TCGA website (https://tcga-data.nci.nih.gov/).

Construction of risk-score formula

After normalization and annotation of GEO raw data, we subjected 2986 lncRNAs to univariable Cox regression proportional hazards regression analysis to select lncRNAs which were associated with OS of NSCLC patients. Those lncRNAs with a statistical significance in univariable Cox regression were then selected into multivariable Cox regression to obtain the coefficients. By linearly combining the expression value of selected lncRNAs weighted by their coefficients, a risk-score formula was constructed as following [20]:

risk score = i=1N(Expi*Coei)

In our formula, the N is the number of 7, the Expi are the expression value of every 7 lncRNAs and the Coei are their corresponding coefficients from the multivariable Cox regression analysis [27].

Survival analysis

According to above formula, the risk scores of NSCLC patients were calculated, by which patients were divided into high- and low-risk group with the cutoff of the median. The Kaplan-Meier method was used to assess the difference in survival time of high- and low- risk NSCLC patients. Difference of survival times between the two groups was considered as significant according to the P value of two-sided log-rank test less than 0.05.

Prognostic meta-analysis

To investigate the combined prognostic values of the seven-lncRNA signature of discovery and two validation groups, we performed a prognostic meta-analysis. All meta-analysis was performed on STATA software, version 12.0 (StataCorp LP, College Station, TX, United States). Pooled HR value was calculated using the random-effects model.

Weighted Correlation Network analysis (WGCNA)

To find the gene modules associated with our risk scores, we construct a co-expression network using the R package “WGCNA” according to our previous reports [28]. The soft thresholding power was selected to 9 to produce a weighted network. The enrolled genes were hierarchically clustered into 11 modules except the gray module.

Pathway enrichment analysis

The pink modules, the most significant modules being associated with risk score, were picked out to perform the pathway enrichment analysis. Pathway enrichment analysis of genes in pink module was performed on an online-based web tool “Metascape” (http://metascape.org/).

Conflicts of Interest

The authors have no conflict of interests to declare.

Funding

This work was supported by the National Natural Science Foundation of China (No. 81871998) and the Natural Science Foundation of Shaanxi Province (Nos. 2018JM7013).

References

  • 1. Fintelmann FJ, Gottumukkala RV, McDermott S, Gilman MD, Lennes IT, Shepard JO. Lung Cancer Screening: Why, When, and How? Radiol Clin North Am. 2017; 55:1163–81. https://doi.org/10.1016/j.rcl.2017.06.003 [PubMed]
  • 2. New M, Keith R. Early Detection and Chemoprevention of Lung Cancer. F1000 Res. 2018; 7:61. https://doi.org/10.12688/f1000research.12433.1 [PubMed]
  • 3. Reungwetwattana T, Dy GK. Targeted therapies in development for non-small cell lung cancer. J Carcinog. 2013; 12:22. https://doi.org/10.4103/1477-3163.123972 [PubMed]
  • 4. Petrella F, Spaggiari L. Therapeutic options following pneumonectomy in non-small cell lung cancer. Expert Rev Respir Med. 2016; 10:919–25. https://doi.org/10.1080/17476348.2016.1188694 [PubMed]
  • 5. Tang Y, Qiao G, Xu E, Xuan Y, Liao M, Yin G. Biomarkers for early diagnosis, prognosis, prediction, and recurrence monitoring of non-small cell lung cancer. OncoTargets Ther. 2017; 10:4527–34. https://doi.org/10.2147/OTT.S142149 [PubMed]
  • 6. i H, Cho JY. Lung Cancer Biomarkers. Adv Clin Chem. 2015; 72:107–70. https://doi.org/10.1016/bs.acc.2015.07.003 [PubMed]
  • 7. Mercer TR, Dinger ME, Mattick JS. Long non-coding RNAs: insights into functions. Nat Rev Genet. 2009; 10:155–59. https://doi.org/10.1038/nrg2521 [PubMed]
  • 8. Rinn JL, Chang HY. Genome regulation by long noncoding RNAs. Annu Rev Biochem. 2012; 81:145–66. https://doi.org/10.1146/annurev-biochem-051410-092902 [PubMed]
  • 9. Kino T, Hurt DE, Ichijo T, Nader N, Chrousos GP. Noncoding RNA gas5 is a growth arrest- and starvation-associated repressor of the glucocorticoid receptor. Sci Signal. 2010; 3:ra8. https://doi.org/10.1126/scisignal.2000568 [PubMed]
  • 10. Wang KC, Yang YW, Liu B, Sanyal A, Corces-Zimmerman R, Chen Y, Lajoie BR, Protacio A, Flynn RA, Gupta RA, Wysocka J, Lei M, Dekker J, et al. A long noncoding RNA maintains active chromatin to coordinate homeotic gene expression. Nature. 2011; 472:120–24. https://doi.org/10.1038/nature09819 [PubMed]
  • 11. Wang JZ, Xiang JJ, Wu LG, Bai YS, Chen ZW, Yin XQ, Wang Q, Guo WH, Peng Y, Guo H, Xu P. A genetic variant in long non-coding RNA MALAT1 associated with survival outcome among patients with advanced lung adenocarcinoma: a survival cohort analysis. BMC Cancer. 2017; 17:167. https://doi.org/10.1186/s12885-017-3151-6 [PubMed]
  • 12. Soudyab M, Iranpour M, Ghafouri-Fard S. The Role of Long Non-Coding RNAs in Breast Cancer. Arch Iran Med. 2016; 19:508–17. [PubMed]
  • 13. Mehra M, Chauhan R. Long Noncoding RNAs as a Key Player in Hepatocellular Carcinoma. Biomark Cancer. 2017; 9:X17737301. https://doi.org/10.1177/1179299X17737301 [PubMed]
  • 14. Tian X, Xu G. Clinical value of lncRNA MALAT1 as a prognostic marker in human cancer: systematic review and meta-analysis. BMJ Open. 2015; 5:e008653. https://doi.org/10.1136/bmjopen-2015-008653 [PubMed]
  • 15. Zhang KJ, Zhang Y, Luo ZL, Liu L, Yang J, Wu LC, Yu SS, Liu JB. [Long non-coding RNA HOTAIR in plasma as a potential biomarker for breast cancer diagnosis]. Nan Fang Yi Ke Da Xue Xue Bao. 2016; 36:488–92. [PubMed]
  • 16. Zhu X, Tian X, Yu C, Shen C, Yan T, Hong J, Wang Z, Fang JY, Chen H. A long non-coding RNA signature to improve prognosis prediction of gastric cancer. Mol Cancer. 2016; 15:60. https://doi.org/10.1186/s12943-016-0544-0 [PubMed]
  • 17. Meng J, Li P, Zhang Q, Yang Z, Fu S. A four-long non-coding RNA signature in predicting breast cancer survival. J Exp Clin Cancer Res. 2014; 33:84. https://doi.org/10.1186/s13046-014-0084-7 [PubMed]
  • 18. Wang Z, Wu Q, Feng S, Zhao Y, Tao C. Identification of four prognostic LncRNAs for survival prediction of patients with hepatocellular carcinoma. PeerJ. 2017; 5:e3575. https://doi.org/10.7717/peerj.3575 [PubMed]
  • 19. Yang K, Hou Y, Li A, Li Z, Wang W, Xie H, Rong Z, Lou G, Li K. Identification of a six-lncRNA signature associated with recurrence of ovarian cancer. Sci Rep. 2017; 7:752. https://doi.org/10.1038/s41598-017-00763-y [PubMed]
  • 20. Tu Z, He D, Deng X, Xiong M, Huang X, Li X, Hao L, Ding Q, Zhang Q. An eight-long non-coding RNA signature as a candidate prognostic biomarker for lung cancer. Oncol Rep. 2016; 36:215–22. https://doi.org/10.3892/or.2016.4817 [PubMed]
  • 21. Negishi M, Wongpalee SP, Sarkar S, Park J, Lee KY, Shibata Y, Reon BJ, Abounader R, Suzuki Y, Sugano S, Dutta A. A new lncRNA, APTR, associates with and represses the CDKN1A/p21 promoter by recruiting polycomb proteins. PLoS One. 2014; 9:e95216. https://doi.org/10.1371/journal.pone.0095216 [PubMed]
  • 22. Luo X, Qiu Y, Jiang Y, Chen F, Jiang L, Zhou Y, Dan H, Zeng X, Lei YL, Chen Q. Long non-coding RNA implicated in the invasion and metastasis of head and neck cancer: possible function and mechanisms. Mol Cancer. 2018; 17:14. https://doi.org/10.1186/s12943-018-0763-7 [PubMed]
  • 23. Tu HY, Ke EE, Yang JJ, Sun YL, Yan HH, Zheng MY, Bai XY, Wang Z, Su J, Chen ZH, Zhang XC, Dong ZY, Wu SP, et al. A comprehensive review of uncommon EGFR mutations in patients with non-small cell lung cancer. Lung Cancer. 2017; 114:96–102. https://doi.org/10.1016/j.lungcan.2017.11.005 [PubMed]
  • 24. Lee CK, Davies L, Wu YL, Mitsudomi T, Inoue A, Rosell R, Zhou C, Nakagawa K, Thongprasert S, Fukuoka M, Lord S, Marschner I, Tu YK, et al. Gefitinib or Erlotinib vs Chemotherapy for EGFR Mutation-Positive Lung Cancer: Individual Patient Data Meta-Analysis of Overall Survival. J Natl Cancer Inst. 2017; 109. https://doi.org/10.1093/jnci/djw279 [PubMed]
  • 25. Dogan S, Shen R, Ang DC, Johnson ML, D’Angelo SP, Paik PK, Brzostowski EB, Riely GJ, Kris MG, Zakowski MF, Ladanyi M. Molecular epidemiology of EGFR and KRAS mutations in 3,026 lung adenocarcinomas: higher susceptibility of women to smoking-related KRAS-mutant cancers. Clin Cancer Res. 2012; 18:6169–77. https://doi.org/10.1158/1078-0432.CCR-11-3265 [PubMed]
  • 26. Marchetti A, Milella M, Felicioni L, Cappuzzo F, Irtelli L, Del Grammastro M, Sciarrotta M, Malatesta S, Nuzzo C, Finocchiaro G, Perrucci B, Carlone D, Gelibter AJ, et al. Clinical implications of KRAS mutations in lung cancer patients treated with tyrosine kinase inhibitors: an important role for mutations in minor clones. Neoplasia. 2009; 11:1084–92. https://doi.org/10.1593/neo.09814 [PubMed]
  • 27. Zhou M, Guo M, He D, Wang X, Cui Y, Yang H, Hao D, Sun J. A potential signature of eight long non-coding RNAs predicts survival in patients with non-small cell lung cancer. J Transl Med. 2015; 13:231. https://doi.org/10.1186/s12967-015-0556-3 [PubMed]
  • 28. Lv B, Zhang L, Miao R, Xiang X, Dong S, Lin T, Li K, Qu K. Comprehensive analysis and experimental verification of LINC01314 as a tumor suppressor in hepatoblastoma. Biomed Pharmacother. 2018; 98:783–92. https://doi.org/10.1016/j.biopha.2018.01.013 [PubMed]