Five Long Non-Coding RNAs Establish a Prognostic Nomogram and Construct a Competing Endogenous RNA Network in the Progression of Non-Small Cell Lung Cancer

doi:10.21203/rs.3.rs-121099/v1

Download PDF

Research article

Five Long Non-Coding RNAs Establish a Prognostic Nomogram and Construct a Competing Endogenous RNA Network in the Progression of Non-Small Cell Lung Cancer

https://doi.org/10.21203/rs.3.rs-121099/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 23 Apr, 2021

Read the published version in BMC Cancer →

You are reading this latest preprint version

Background: Accumulating evidence has revealed that long non-coding RNAs (lncRNAs) play vital roles in the progression of non-small cell lung cancer (NSCLC). But the relationship between lncRNAs and survival outcome of NSCLC remains to be explored. Therefore, we attempt to figure out their survival roles and molecular connection in NSCLC.

Methods: By analyzing the transcriptome profiling of NSCLC from TCGA databases, we divided patients into three groups, and identified differentially expressed lncRNAs (DELs) of each group. Next, we explored the prognostic roles of common DELs by univariate and multivariate Cox analysis, LASSON, and Kaplan-Meier analysis. Additionally, we assessed and compared the prognostic accuracy of 5 lncRNAs through ROC curves and AUC values. Ultimately, we detected their potential function by enrichment analysis and molecular connection through establishing a ceRNA network.

Results: 197 common DELs were spotted. And we successfully screened out 5 lncRNAs related to the patient's survival, including LINC01833, AC112206.2, FAM83A-AS1, BANCR, and HOTAIR. Combing with age and AJCC stage, we constructed a nomogram that prognostic prediction was superior to the traditional parameters. Furthermore, 275 qualified mRNAs related to 5 lncRNAs were spotted. Functional analysis indicates that these lncRNAs act key roles in the progression of NSCLC, such as P53 and cell cycle signaling pathway. And ceRNA network also suggests that these lncRNAs are tightly connected with tumor progression.

Conclusions: A nomogram and ceRNA network based on 5 lncRNAs indicate that there can effectively predict the overall survival of NSCLC and potentially serve as a therapeutic guide for NSCLC.

Cancer Biology

Oncology

NSCLC

long non-coding RNA

nomogram

ceRNA

overall survival

Non-small cell lung cancer (NSCLC) is one of the most common and deadly cancers in the world. Despite advances in treatments, only 19% of patients with NSCLC survival for more than 5 years [1, 2]. What's worse, due to lacking specific symptoms in the early stage, most patients seek treatment at an advanced stage, which misses the best timing for a radical operation [3]. Therefore, it is urgent to establish a prognostic risk-score model for NSCLC patients, thereby providing a therapeutic guide for NSCLC.

Long non-coding RNAs (lncRNAs) refers to non-protein-coding transcripts over 200 nucleotides in length [4]. Although lncRNAs do not directly encode RNA, it can regulate protein expression at various stages of transcription [5, 6]. Based on the ceRNA hypothesis, lncRNAs, messenger RNAs, and pseudogenes can “talk” to each other using miRNA response elements (MREs) and assemble as a ceRNA network [7]. In this network, lncRNA act as "sponges" to absorb and bind miRNA, thereby weakening their binding ability to mRNA and regulating gene expression. Accumulating evidence has backed that lncRNAs were involved in the ceRNA network of many types of cancers, including pancreatic cancer, gastric cancer, as well as NSCLC [8–10]. Notably, LncRNAs have great advantages as biomarkers because they are stable, highly tissue-specific, and easy to detect in body fluids. Besides, some lncRNAs have been recognized as a novel biomarker, for instance, BANCR in gastric carcinoma, HOTAIR in colorectal carcinoma, and MALAT1 in lung cancer [11–13]. Therefore, it is needful to detect the prognostic relationship between lncRNAs and NSCLC.

In this study, we performed a large sample analysis to find out survival-related lncRNAs and validated it using univariate, LASSO, and multivariate Cox proportional hazards regression (CPHR). Moreover, combing with age and AJCC stage, we constructed a nomogram based on lncRNAs, which performed better prognostic prediction than clinical factors, and we successfully assessed its efficiency in LUAD and LUSC groups. Functional analysis indicated that these lncRNAs also act important roles in the progression of NSCLC, for example, the P53 signaling pathway and cell cycle pathway. Next, we successfully constructed a ceRNA network related to 5 lncRNAs. Those results indicated that those lncRNAs not only effectively predict the prognosis of NSCLC patients but also take part in the progression of NSCLC and potentially serve as a therapeutic target.

Data selection and process

As is well know that the majority of NSCLC is composed of lung squamous carcinoma (LUSC) and lung adenocarcinoma (LUAD). So we obtained the raw counts of transcriptome profiling and clinical data from TCGA-LUAD including 535 cancer samples and 59 non-tumor tissues, and data from TCGA-LUSC including 502 cancer samples and 49 non-tumor tissues. And we divided those data into three groups including the NSCLC group, LUAD group, and LUSC group. All the transcriptome data of lncRNAs, miRNAs, and mRNAs were obtained from The Cancer Genome Atlas (TCGA) database and were annotated through the Ensemble database [14, 15]. Stepwise, all the raw count data were log2 (x+1) transformed and normalized using the "LIMMA" package [16]. Then, we screened out differentially expressed miRNAs (DEMIs) in three groups by the “LIMMA” package with the criteria of |log2FC| > 1, average expression > 1, and adjust P-value < 0.05, and identified differentially expressed lncRNAs (DELs) and differentially expressed mRNAs (DEMs) in three groups with the criteria of |log2FC| > 2, average expression > 2, and adjust P-value < 0.05, respectively.

Survival analysis

Additionally, we removed samples without survival time or survival time less than 7 days to improve the reliability of our study. We first estimated the association between overall survival (OS) and clinical parameters through univariate and multivariate CPHR analysis. Furthermore, we evaluated the relationship between survival time and common DELs expression through Kaplan Meier analysis and univariate CPHR method. Only DELs that their P-value was lower than 0.05 and their expression consistent with prognosis were regarded as candidate survival-related lncRNAs. Combing with clinical risk factors, we performed LASSON Cox regression analysis to obtain the best fitting variables. After selecting the best-fit of OS-related variables by the calculation mentioned above, we further verified their prognostic value through multivariate CPHR analysis. And only variables that P-value was lower than 0.05 in univariate and multivariate CPHR calculations were selected as qualified lncRNAs and were chose for the next analysis.

Constructing a risk score formula and nomogram model

Combining with clinical risk variables and qualified lncRNAs, we performed univariate and multivariate CPHR analysis to identify OS-related biomarkers in 970 patients with NSCLC. Also, we calculated the prognostic risk score of each patient through multivariate CPHR analysis according to the formula as follows: risk score = X1α1 + X2α2 + X3α3 +...+ Xnαn. And patients in the NSCLC group were divided into high-risk groups and low-risk groups based on the median value of risk score. C-indexes were performed to assess the predictive performances of our risk score formula. Next, we evaluated the prognostic differences in high-risk and low-risk groups by Kaplan Meier analysis and T-test. We assessed the prediction performance of the risk formula through ROC curves at 3 and 5 years and computed their AUC values in three groups. Moreover, we established a nomogram to vividly depict the predictive relationship among clinical factors, lncRNAs, and OS. Calibration curves of 3 and 5 years were calculated to assess the reliability of OS prediction between predicted performance and actual ability. All the analyses mentioned above were conducted in NSCLC, LUAD, and LUSC groups, respectively.

Function enrichment analysis

We performed the Gene Ontology (GO) terms and Kyoto Encyclopedia of Gene and Genomes (KEGG) pathways enrichment analysis to elucidate the potential functions of lncRNAs in the nomogram. The common DEMs in the three groups were first identified. Next, the correlation coefficient of each lncRNA with common DEMs was calculated, respectively. To obtain an accurate result, we only selected DEMs with a correlation coefficient greater than 0.2 for further enrichment analysis. Then, we performed the GO and KEGG analysis of lncRNA-related DEMs via the DAVID database, and the Enrichr database [17, 18]. And we depicted the top 10 enriched GO terms and KEGG pathways through a bar plot with the criteria of adjusting P-value < 0.05.

Establishing a ceRNA network

To explore the potential interaction between 5 lincRNAs and miRNAs, the LncBase database that provided miRNAs and lncRNAs interactions according to MREs sites was applied to predict the downstream miRNA of 5 lncRNAs with the criteria of Prediction score > 0.8 [19]. Only miRNAs expressed on the NSCLC group, LUAD group, LUSC group, and LncBase database were chosen as qualified miRNAs. Additionally, the miRDB database and the miRTarBase database are a widely-used tool for miRNA target prediction that were employed to find out potential mRNAs binding to qualified miRNAs [20, 21]. Only miRNAs that validated in miRDB database, miRTarBase database, and DEMs were considered as candidate mRNAs. Finally, to vividly display the interaction of 5 lncRNAs with qualified miRNAs and candidate mRNAs, we constructed a ceRNA network by Cytoscape software (Version 3.7.2) [22].

Statistical analysis

We performed all the statistical analyses mentioned above using R software (version 3.6.1). Briefly, OS was analyzed using the Kaplan–Meier test, and the log-rank T-test was applied to calculate the statistical significance. Univariate and multivariate CPHR analyses were conducted through "Survival" packages. LASSON CPHR was performed using "Survival" and "glmnet" packages. A nomogram was constructed by "Survival" and "rms" packages. And a time-dependent ROC analysis was performed by the ‘‘survivalROC’’ package, C-index by ‘‘survival’’ package, and calibration curve by ‘‘rms’’ package.

The ceRNA network was constructed by Cytoscape software. And we set a p-value lower than 0.05 as s statistical significance.

Screening differentially expressed RNAs

A total of 1145 NSCLC patients were enrolled from the TCGA database. We divided those data into three groups including the NSCLC group, LUAD group, and LUSC group. Then, we identified the DELs and DEMs in three groups by standards of |log2FC| >2, average expression > 2, and adjust P value < 0.05, respectively. There are a total of 426 DELs in the NSCLC group, 312 DELs in the LUAD group, 687 DELs in the LUSC group, and 197 common DELs in three groups (Fig. 1A–D). There is a total of 1905 DEMs in the NSCLC group, 1434 DEMs in the LUAD group, 2641 DEMs in the LUSC group, and 1131 common DEMs in three groups (Supplementary Fig. 1A–D). And there is a total of 130 DEMIs in the NSCLC group, 133 DEMIs in the LUAD group, 163 DEMIs in the LUSC group, and 89 common DEMIs in three groups (Supplementary Fig. 5A–D) with the criteria of |log2FC| > 1, average expression > 1, and adjust P-value < 0.05. Also, we present the top 10 up and down-regulated DELs spotted in three groups (Fig. 1E). Common DELs in three groups were selected for the next survival analysis.

Developing A Risk Score Formula And Prognostic Nomogram

Using multivariable CPHR analysis, we divided NSCLC patients into high-risk and low-risk groups based on the value of the median risk score. And we developed a risk score formula to understand the relationship between overall survival and lncRNAs with clinical variables. The formula was as shown as follows: Risk Score = (0.152 × age) + (0.370 × AJCC) + (0.127 × Expression LINC01833) + (0.250 × Expression FAM83A-AS1) + (0.127 × Expression HOTAIR) - (0.124 × Expression BANCR) - (0.283 × Expression AC112206.2). We depicted the distributions of risk score and the status of overall survival in the NSCLC group (Fig. 4A). Also, we validate the distributions of risk score and the status of overall survival in the LUAD group (Supplementary Fig. 2A) and LUSC group (Supplementary Fig. 2B). Furthermore, we found that the high-risk group of NSCLC patients had a worse prognostic outcome than patients with lower risk scores (Fig. 4B). And this tendency also appeared in the LUAD group (Supplementary Fig. 2C) and LUSC group (Supplementary Fig. 2D). We performed the time-related ROC curves to compare the sensitivities and specificities of predictive formula in NSCLC. The result suggested that the AUC value of 3 and 5 years was 0.703, 0.667, respectively (Fig. 4C). But the AUC value of clinical variables was only 0.632 and 0.618, which indicate that the prognostic prediction of our nomogram was better than the traditional age and AJCC stage (Fig. 4D). Also, we assessed the predictive ability of formula in LUAD and LUSC group. In the LUAD group, the AUC value of 3 and 5 years was 0.749, 0.735, respectively (Supplementary Fig. 2E). The AUC value of 3 and 5 years was 0.65, 0.619 in the LUSC group, respectively (Supplementary Fig. 2F). To vividly displayed the prognostic performance of OS-related variables, a nomogram was established. As shown in Fig. 5A, the nomogram could usefully predict the prognosis of 3 years and 5 years in NSCLC patients (Fig. 5A). Additionally, calibration curves suggested that the nomogram had a superior agreement between the predicted and actual OS of 3-year and 5-year in the NSCLC group (Fig. 5B-C) as well as in the LUAD group(Supplementary Fig. 3A-B) and LUSC group(Supplementary Fig. 3C-D).

Construction Of A Cerna Network

Based on the miRNA prediction from 5 lncRNAs, we identified 294 downstream miRNAs by the LncBase database (Supplementary Table S3). And we evaluated the expression of 294 miRNAs in the NSCLC group, LUAD group, and LUSC group. There are a total of 20 common miRNAs in four groups (Supplementary Fig. 5D). Next, miRDB and miRTarBase databases were employed for screening miRNA-linked mRNAs. And we validated their expression in the DEMs of three groups. There are a total of 91 common mRNAs in four groups. According to prediction in Supplementary Table S4, we have found out 22 pairs of lncRNA-miRNA interactions, 145 pairs of miRNA-mRNA interactions. Finally, we constructed a ceRNA network in Fig. 6 to vividly display the interactions of 5 lncRNAs with 20 miRNAs and 91 mRNAs. Overall, this evidence has revealed that those lncRNAs not only effectively predict the survival outcome of NSCLC patients but also take part in the progression of NSCLC and potentially serve as a therapeutic target.

It is well acknowledged that the AJCC stage has been extensively used to estimate the survival outcome of tumor patients [23]. However, some limitations of the AJCC stage can be found in our clinical practice. For example, patients with similar anatomic sites and AJCC staging can exhibit variable responses to treatment and different survival outcomes. This difference may result from tumor heterogeneity, which partly arises from genetic mutations [24, 25]. Furthermore, recent studies have indicated that age and gender are also effective predictors for OS [26, 27]. So, we attempted to establish a new staging system that combines clinical variables with genetic mutations. Also, some evidence has indicated that lncRNAs not only act a regulatory role in the progression of NSCLC, but also have great potentials as biomarkers because they are stable, highly tissue-specific, and easy to detect in body fluids. For example, serum exosomal MALAT-1 was identified as a diagnostic predictor for NSCLC patients when they are in early-phase or metastasis [13]. High expression of HOTAIR was closely related to progressive disease, worse survival, and more potential in tumor recurrence after radical operation [28]. Low expression of GAS6-AS1 was connected with the occurrence of lymph node metastasis and an independent biomarker for the prognostic outcome of NSCLC [29]. Therefore, it is urgent to set up a reliable prognostic model for NSCLC.

In this research, we successfully identified a survival nomogram based on lncRNAs and clinical variables of NSCLC. We combined age and AJCC staging with 5 lncRNAs into a risk formula and weighted each parameter to detect their relationship with overall survival. And a nomogram was established to vividly quantify the OS probability of each variable. Notably, we weighted each lncRNA into a nomogram instead of integrating 5 lncRNAs as a whole [30, 31]. Because it is difficult to test all variables at a time in clinical practice. In our nomogram, we weighted a risk point to every variable and calculated patient survival at 3 or 5 years based on the total risk point. And the predictive performance of our nomogram was superior to the traditional age and AJCC stage. Moreover, our nomogram is easy to understand. Its simplicity will allow clinicians to quickly evaluate survival outcomes and make decisions about individual NSCLC patients. Even individuals without a medical background can easily understand the meaning of our nomogram. Those features will make our nomogram an accurate and effective biomarker for clinical applications.

Additionally, the functional analysis indicated that 5 lncRNAs were also involved in several cancer pathways, for example, p53 and MicroRNAs in cancer pathways. So, the lncRNAs in the nomogram can not only serve as a prognostic biomarker but also function as a regulator in the occurrence and progression of tumors. Moreover, some studies have found their role in the progression of NSCLC. The high expression of HOTAIR is involved in many types of biological processes in NSCLC. For instance, HOTAIR contributes to the down-regulating expression of p21WAF1/CIP1, thereby inducing the cisplatin resistance of A549 cells [32]. Wang et al. demonstrated that HOTAIR was reported to facilitate the proliferation and migration of A549 and H838 cells through sponging with miR-326, thus control the expression of phox2a [33]. The high expression of HOTAIR can be controlled by Col-1, thereby promoting the formation of microenvironment and progression of NSCLC [34]. Several studies have recently revealed that the lower expression of BANCR was related to the initiation and progression of NSCLC. Sun et al. found that lower expression of BANCR can encourage the epithelial-mesenchymal transition of A549 and SPC-A1 cells and improve their ability in invasion and metastasis [35]. Up-regulating BANCR was tightly linked with radiotherapy for lung cancer [36]. Jiang et al. discovered that BANCR was able to moderate the ability of invasion and metastasis in lung cancer through the p38 MAPK and JNK pathway rather than the ERK MAPK pathway [37]. Notably, those lncRNAs can act as a ceRNA network, thereby participating in cancer progression. For instance, HOTAIR/miR-149-5p/HNRNPA1 axis promotes the cell growth, migration, and invasion in NSCLC [38]; FAM83A-AS1/miR-150-5p/MMP14 regulates LUAD progression and invasion [39]; BANCR/miR-338-3p/IGF1R network regulated Raf/MEK/ERK pathway, thereby encouraging the proliferation, migration, invasion and epithelial-mesenchymal transition (EMT) of esophageal cancer [40]. Those studies indicated that the lncRNAs in the nomogram can function as a therapeutic target for NSCLC.

Although our lncRNAs-related nomogram showed good performance in survival prediction of NSCLC, some limitations can be detected in our research. First, our study is a retrospective study. But 970 patients is a large sample. So our results are acceptable. Additionally, our data lacked information such as chemo-radiotherapy history, smoking history, and patients' disease history. This may result from the limitations of our data. So, further clinical studies are needed to verify our results when applied to clinical practice. Last, we identified survival-related lncRNAs to construct the nomogram, which might overlook some valuable information. All in all, despite these limitations in our study, we believe that our persistent efforts will eventually establish an ideal prognostic model in clinical practice.

In summary, by a large sample analysis, we successfully constructed a nomogram based on lncRNAs and clinical variables that predicts the survival of NSCLC patients. And the predictive performance of our prognostic nomogram was better than the traditional AJCC stage and age. In addition to the survival prediction of our nomogram, functional analysis and ceRNA network also indicate that they might involve in cancer progression and potentially serve as a therapeutic target for NSCLC.

NSCLC: non-small cell lung cancer; LUAD: lung adenocarcinoma; LUSC: lung squamous carcinoma;

lncRNAs: long non-coding RNAs; CPHR: Cox proportional hazards regression; TCGA: The Cancer Genome Atlas; GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Gene and Genomes; AJCC: American Joint Committee on Cancer; OS: overall survival; HR: hazard ratio; CI: confidence interval; ROC: receiver operating curves; AUC: Area Under Curve; ceRNA: competing endogenous RNA.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Availability of data and material

All data generated or analyzed during this study are included in this published article and its supplementary information files.

Competing interests

The authors declare no competing interests.

Funding

Not applicable.

Authors' contributions

YY, RKM designed this project, performed data analysis and wrote the manuscript. All authors have read and approved the final manuscript.

Acknowledgements

Not applicable.

Siegel R, Miller K, Jemal A. Cancer statistics, 2019. Cancer J Clin. 2019;69(1):7–34.
New M, Keith R: Early Detection and Chemoprevention of Lung Cancer. F1000Research 2018, 7:61.
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 therapy. 2017;10:4527–34.
Morris K, Mattick J. The rise of regulatory RNA. Nature reviews Genetics. 2014;15(6):423–37.
Engreitz J, Haines J, Perez E, Munson G, Chen J, Kane M, McDonel P, Guttman M, Lander E. Local regulation of gene expression by lncRNA promoters, transcription and splicing. Nature. 2016;539(7629):452–5.
Bhan A, Soleimani M, Mandal S. Long Noncoding RNA and Cancer: A New Paradigm. Cancer research. 2017;77(15):3965–81.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP: A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell 2011, 146(3):353–358.
Xiong G, Liu C, Yang G, Feng M, Xu J, Zhao F, You L, Zhou L, Zheng L, Hu Y, et al. Long noncoding RNA GSTM3TV2 upregulates LAT2 and OLR1 by competitively sponging let-7 to promote gemcitabine resistance in pancreatic cancer. J Hematol Oncol. 2019;12(1):97.
Zu F, Han H, Sheng W, Sun J, Zang H, Liang Y, Liu Q. Identification of a competing endogenous RNA axis related to gastric cancer. Aging (Albany NY) 2020, 12.
Chen X, Wang Z, Tong F, Dong X, Wu G, Zhang R. lncRNA UCA1 Promotes Gefitinib Resistance as a ceRNA to Target FOSL2 by Sponging miR-143 in Non-small Cell Lung Cancer. Mol Ther Nucleic Acids. 2020;19:643–53.
Li L, Zhang L, Zhang Y, Zhou F. Increased expression of LncRNA BANCR is associated with clinical progression and poor prognosis in gastric cancer. Biomedicine pharmacotherapy = Biomedecine pharmacotherapie. 2015;72:109–12.
Shen X, Bai Y, Luo B, Zhou X. Upregulation of lncRNA BANCR associated with the lymph node metastasis and poor prognosis in colorectal cancer. Biological research. 2017;50(1):32.
Zhang R, Xia Y, Wang Z, Zheng J, Chen Y, Li X, Wang Y, Ming H. Serum long non coding RNA MALAT-1 protected by exosomes is up-regulated and promotes cell proliferation and migration in non-small cell lung cancer. Biochem Biophys Res Commun. 2017;490(2):406–14.
Weinstein J, Collisson E, Mills G, Shaw K, Ozenberger B, Ellrott K, Shmulevich I, Sander C, Stuart J. The Cancer Genome Atlas Pan-Cancer analysis project. Nat Genet. 2013;45(10):1113–20.
Yates A, Achuthan P, Akanni W, Allen J, Allen J, Alvarez-Jarreta J, Amode M, Armean I, Azov A, Bennett R, et al: Ensembl 2020. Nucleic acids research 2020, 48:D682-D688.
Ritchie M, Phipson B, Wu D, Hu Y, Law C, Shi W, Smyth G. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic acids research. 2015;43(7):e47.
Huang dW, Sherman B, Lempicki R. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature protocols. 2009;4(1):44–57.
Kuleshov M, Jones M, Rouillard A, Fernandez N, Duan Q, Wang Z, Koplev S, Jenkins S, Jagodnik K, Lachmann A, et al: Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic acids research 2016, 44:W90-97.
Karagkouni D, Paraskevopoulou MD, Tastsoglou S, Skoufos G, Karavangeli A, Pierros V, Zacharopoulou E, Hatzigeorgiou AG. DIANA-LncBase v3: indexing experimentally supported miRNA targets on non-coding transcripts. Nucleic acids research. 2020;48(D1):D101–10.
Chou C-H, Shrestha S, Yang C-D, Chang N-W, Lin Y-L, Liao K-W, Huang W-C, Sun T-H, Tu S-J, Lee W-H, et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic acids research. 2018;46(D1):D296–302.
Chen Y, Wang X. miRDB: an online database for prediction of functional microRNA targets. Nucleic acids research. 2020;48(D1):D127–31.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Rami-Porta R, Asamura H, Travis W, Rusch V. Lung cancer - major changes in the American Joint Committee on Cancer eighth edition cancer staging manual. Cancer J Clin. 2017;67(2):138–55.
Jain D, Roy-Chowdhuri S. Molecular Pathology of Lung Cancer Cytology Specimens: A Concise Review. Arch Pathol Lab Med. 2018;142(9):1127–33.
Rooney C, Sethi T. Advances in molecular biology of lung disease: aiming for precision therapy in non-small cell lung cancer. Chest. 2015;148(4):1063–72.
Radkiewicz C, Dickman P, Johansson A, Wagenius G, Edgren G, Lambe M. Sex and survival in non-small cell lung cancer: A nationwide cohort study. PloS one. 2019;14(6):e0219206.
Ahmed T, Lycan T, Dothard A, Ehrlichman P, Ruiz J, Farris M, Topaloglu U, Levine B, Grant S, Klepin H, et al. Performance Status and Age as Predictors of Immunotherapy Outcomes in Advanced Non-Small-Cell Lung Cancer. Clin Lung Cancer. 2020;21(4):e286–93.
Liu X, Liu Z, Sun M, Liu J, Wang Z, De W. The long non-coding RNA HOTAIR indicates a poor prognosis and promotes metastasis in non-small cell lung cancer. BMC Cancer. 2013;13:464.
Han L, Kong R, Yin D, Zhang E, Xu T, De W, Shu Y: Low expression of long noncoding RNA GAS6-AS1 predicts a poor prognosis in patients with NSCLC. Medical oncology (Northwood, London, England) 2013, 30(4):694.
Wang Y, Du L, Yang X, Li J, Li P, Zhao Y, Duan W, Chen Y, Wang Y, Mao H, et al. A nomogram combining long non-coding RNA expression profiles and clinical factors predicts survival in patients with bladder cancer. Aging. 2020;12(3):2857–79.
Liu Z, Sun M, Lu K, Liu J, Zhang M, Wu W, De W, Wang Z, Wang R: The long noncoding RNA HOTAIR contributes to cisplatin resistance of human lung adenocarcinoma cells via downregualtion of p21(WAF1/CIP1) expression. PloS one 2013, 8(10):e77293.
Wang R, Chen X, Xu T, Xia R, Han L, Chen W, De W, Shu Y. MiR-326 regulates cell proliferation and migration in lung cancer by targeting phox2a and is regulated by HOTAIR. American journal of cancer research. 2016;6(2):173–86.
Zhuang Y, Wang X, Nguyen H, Zhuo Y, Cui X, Fewell C, Flemington E, Shan B. Induction of long intergenic non-coding RNA HOTAIR in lung cancer cells by type I collagen. J Hematol Oncol. 2013;6:35.
Lin T, Fu Y, Zhang X, Gu J, Ma X, Miao R, Xiang X, Niu W, Qu K, Liu C, et al. A seven-long noncoding RNA signature predicts overall survival for patients with early stage non-small cell lung cancer. Aging. 2018;10(9):2356–66.
Sun M, Liu X, Wang K, Nie F, Kong R, Yang J, Xia R, Xu T, Jin F, Liu Z, et al. Downregulation of BRAF activated non-coding RNA is associated with poor prognosis for non-small cell lung cancer and promotes metastasis by affecting epithelial-mesenchymal transition. Mol Cancer. 2014;13:68.
Chen J, Chen M, Zheng Y, Wang S, Shen Z. Up-regulation of BRAF activated non-coding RNA is associated with radiation therapy for lung cancer. Biomedicine pharmacotherapy = Biomedecine pharmacotherapie. 2015;71:79–83.
Jiang W, Zhang D, Xu B, Wu Z, Liu S, Zhang L, Tian Y, Han X, Tian D. Long non-coding RNA BANCR promotes proliferation and migration of lung carcinoma via MAPK pathways. Biomedicine pharmacotherapy = Biomedecine pharmacotherapie. 2015;69:90–5.
Li H, Cui Z, Lv X, Li J, Gao M, Yang Z, Bi Y, Zhang Z, Wang S, Li S, et al. Long Non-coding RNA HOTAIR Function as a Competing Endogenous RNA for miR-149-5p to Promote the Cell Growth, Migration, and Invasion in Non-small Cell Lung Cancer. Front Oncol. 2020;10:528520.
Xiao G, Wang P, Zheng X, Liu D, Sun X. FAM83A-AS1 promotes lung adenocarcinoma cell migration and invasion by targeting miR-150-5p and modifying MMP14. Cell Cycle. 2019;18(21):2972–85.
Song W, Wang K, Yang X, Dai W, Fan Z. Long non–coding RNA BANCR mediates esophageal squamous cell carcinoma progression by regulating the IGF1R/Raf/MEK/ERK pathway via miR–338–3p. Int J Mol Med. 2020;46(4):1377–88.

SupplementaryFigure1.tif
Supplementary Figure 1. Screening differentially expressed mRNAs (DEMs) in three groups. (A–C) The volcano plots of DEMs in the TCGA_NSCLC group, TCGA_LUAD group, and TCGA_LUSC group with thresholds of |log2FC| > 2, average expression > 2, and adjust P-value < 0.05, respectively. The red dots and blue dots represent the up-regulated and down-regulated DEMs, separately. (D) The intersection of DELs in three groups.
SupplementaryFigure2.tif
Supplementary Figure 2. Assessing the prognostic performance of the risk score formula in LUAD and LUSC group. (A-B) The risk score distribution and OS status of the formula in LUAD and LUSC group, respectively. (C-D) Kaplan-Meier curves for OS based on the formula in LUAD and LUSC group, separately. The tick-marks on the curve represent the censored patients. (E-F) ROC curve analysis of the formula for predicting OS in LUAD and LUSC group, respectively.
SupplementaryFigure3.tif
Supplementary Figure 3. Evaluating the prediction performance of the nomogram. (A-B) Calibration curves of a nomogram to evaluate the prediction performance of 3-years and 5-years in the LUAD group. (C-D) Calibration curves of a nomogram to evaluate the prediction performance of 3-years and 5-years in the LUSC group.
SupplementaryFigure4.tif
Supplementary Figure 4. Functional enrichment analysis for the differentially expressed mRNAs (DEMs) of 5 lncRNAs. (A-C) The top ten enriched GO terms of qualified DEMs with a correlation coefficient greater than 0.2 in biological processes (BP), cellular components (CC), and molecular functions (MF), respectively. (D) The top ten enriched KEGG pathways of qualified DEMs in the NSCLC database.
SupplementaryFigure5.tif
Supplementary Figure 5. Screening differentially expressed miRNA (DEMIs) in three groups. (A–C) The volcano plots of DEMIs in the TCGA_NSCLC group, TCGA_LUAD group, and TCGA_LUSC group with thresholds of |log2FC| > 1, average expression > 1, and adjust P-value < 0.05, respectively. The red dots and blue dots represent the up-regulated and down-regulated DEMIs, separately. (D) The intersection of DEMIs in three groups and the LncBase database.
Supplementarytable1TableS1.docx
Supplementary Table S1. Clinical characteristics of NSCLC patients in this study.
Supplementarytable2TableS2.docx
Supplementary Table S2. Univariate and multivariate Cox regression analysis in NSCLC.
Supplementarytable3TableS3.docx
Supplementary Table S3. The lncRNA-miRNA pairs were predicted by the LncBase database.
Supplementarytable4TableS4.docx
Supplementary Table S4. The qualified lncRNA-miRNA and miRNA-mRNA pairs were constructed in the ceRNA network.

Download PDF

Journal Publication

published 23 Apr, 2021

Read the published version in BMC Cancer →

Editor assigned by journal
10 Dec, 2020
Reviewers invited by journal
10 Dec, 2020
Submission checks completed at journal
03 Dec, 2020
Editor invited by journal
24 Nov, 2020
First submitted to journal
02 Nov, 2020

You are reading this latest preprint version

Five Long Non-Coding RNAs Establish a Prognostic Nomogram and Construct a Competing Endogenous RNA Network in the Progression of Non-Small Cell Lung Cancer

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Materials And Methods

Result

Screening differentially expressed RNAs

Developing A Risk Score Formula And Prognostic Nomogram

Construction Of A Cerna Network

Discussion

Conclusions

Abbreviations

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1