Mutation-Derived Long Noncoding RNA Signature Predicts Survival in Lung Adenocarcinoma.

Genomic instability is one of the representative features of cancer evolution. Recent research has revealed that long noncoding RNAs (lncRNAs) play a critical role in maintaining genomic instability. Our work proposed a gene signature (GILncSig) based on genomic instability-derived lncRNAs to probe the possibility of lncRNA signatures as an index of genomic instability, providing a potential new approach to identify genomic instability-related cancer biomarkers. Lung adenocarcinoma (LUAD) gene expression data from an RNA-seq FPKM dataset, somatic mutation information and relevant clinical materials were downloaded from The Cancer Genome Atlas (TCGA). A prognostic model consisting of genomic instability-related lncRNAs was constructed, termed GILncSig, to calculate the risk score. We validated GILncSig using data from the Gene Expression Omnibus (GEO) database. In this study, we used R software for data analysis. Through univariate and multivariate Cox regression analyses, five genomic instability-associated lncRNAs (LINC01671, LINC01116, LINC01214, lncRNA PTCSC3, and LINC02555) were identified. We constructed a lncRNA signature (GILncSig) related to genomic instability. LUAD patients were classified into two risk groups by GILncSig. The results showed that the survival rate of LUAD patients in the low-risk group was higher than that of those in the high-risk group. Then, we verified GILncSig in the GEO database. GILncSig was associated with the genomic mutation rate of LUAD. We also used GILncSig to divide TP53 mutant-type patients and TP53 wild-type patients into two groups and performed prognostic analysis. The results suggested that compared with TP53 mutation status, GILncSig may have better prognostic significance. By combining the lncRNA expression profiles associated with somatic mutations and the corresponding clinical characteristics of LUAD, a lncRNA signature (GILncSig) related to genomic instability was established.

Yang L ，Guo G ，Yu X ，Wen Y ，Lin Y ，Zhang R ，Zhao D ，Huang Z ，Wang G ，Yan Y ，Zhang X ，Chen D ，Xing W ，Wang W ，Zeng W ，Zhang L ... -  《Frontiers in Oncology》

Identification of a Novel Prognostic Signature of Genome Instability-Related LncRNAs in Early Stage Lung Adenocarcinoma.

Increasing evidence has demonstrated that long non-coding RNAs (lncRNAs) play a crucial part in maintaining genomic instability. We therefore identified genome instability-related lncRNAs and constructed a prediction signature for early stage lung adenocarcinoma (LUAD) as well in order for classification of high-risk group of patients and improvement of individualized therapies. Early stage LUAD RNA-seq and clinical data from The Cancer Genome Atlas (TCGA) were randomly divided into training set (n = 177) and testing set (n = 176). A total of 146 genomic instability-associated lncRNAs were identified based on somatic mutation profiles combining lncRNA expression profiles from TCGA by the "limma R" package. We performed Cox regression analysis to develop this predictive indicator. We validated the prognostic signature by an external independent LUAD cohort with microarray platform acquired from the Gene Expression Omnibus (GEO). A genome instability-related six-lncRNA-based gene signature (GILncSig) was established to divide subjects into high-risk and low-risk groups with different outcomes at statistically significant levels. According to the multivariate Cox regression and stratification analysis, the GILncSig was an independent predictive factor. Furthermore, the six-lncRNA signature achieved AUC values of 0.745, 0.659, and 0.708 in the training set, testing set, and TCGA set, respectively. When compared with other prognostic lncRNA signatures, the GILncSig also exhibited better prediction performance. The prognostic lncRNA signature is a potent tool for risk stratification of early stage LUAD patients. Our study also provided new insights for identifying genome instability-related cancer biomarkers.

Peng B ，Li H ，Na R ，Lu T ，Li Y ，Zhao J ，Zhang H ，Zhang L ... -  《-》

Prognostic Value of Genomic Instability of m(6)A-Related lncRNAs in Lung Adenocarcinoma.

Background: Genomic instability of N6-methyladenosine (m6A)-related long noncoding RNAs (lncRNAs) plays a pivotal role in the tumorigenesis of lung adenocarcinoma (LUAD). Our study identified a signature of genomic instability of m6A-associated lncRNA signature and revealed its prognostic role in LUAD. Methods: We downloaded RNA-sequencing data and somatic mutation data for LUAD from The Cancer Genome Atlas (TCGA) and the GSE102287 dataset from the Gene Expression Omnibus (GEO) database. The "Limma" R package was used to identify a network of regulatory m6A-related lncRNAs. We used the Wilcoxon test method to identify genomic-instability-derived m6A-related lncRNAs. A competing endogenous RNA (ceRNA) network was constructed to identify the mechanism of the genomic instability of m6A-related lncRNAs. Univariate and multivariate Cox regression analyses were performed to construct a prognostic model for internal testing and validation of the prognostic m6A-related lncRNAs using the GEO dataset. Performance analysis was conducted to compare our prognostic model with the previously published lncRNA models. The CIBERSORT algorithm was used to explore the relationship of m6A-related lncRNAs and the immune microenvironment. Prognostic m6A-related lncRNAs in prognosis, the tumor microenvironment, stemness scores, and anticancer drug sensitivity were analyzed to explore the role of prognostic m6A-related lncRNAs in LUAD. Results: A total of 42 genomic instability-derived m6A-related lncRNAs were differentially expressed between the GS (genomic stable) and GU (genomic unstable) groups of LUAD patients. Four differentially expressed lncRNAs, 17 differentially expressed microRNAs, and 75 differentially expressed mRNAs were involved in the genomic-instability-derived m6A-related lncRNA-mediated ceRNA network. A prediction model based on 17 prognostic m6A-associated lncRNAs was constructed based on three TCGA datasets (all, training, and testing) and validated in the GSE102287 dataset. Performance comparison analysis showed that our prediction model (area under the curve [AUC] = 0.746) could better predict the survival of LUAD patients than the previously published lncRNA models (AUC = 0.577, AUC = 0.681). Prognostic m6A-related-lncRNAs have pivotal roles in the tumor microenvironment, stemness scores, and anticancer drug sensitivity of LUAD. Conclusion: A signature of genomic instability of m6A-associated lncRNAs to predict the survival of LUAD patients was validated. The prognostic, immune microenvironment and anticancer drug sensitivity analysis shed new light on the potential novel therapeutic targets in LUAD.

Li R ，Li JP ，Liu TT ，Huo C ，Yao J ，Ji XL ，Qu YQ ... -  《Frontiers in Cell and Developmental Biology》

Mutator-Derived lncRNA Landscape: A Novel Insight Into the Genomic Instability of Prostate Cancer.

Increasing evidence has emerged to reveal the correlation between genomic instability and long non-coding RNAs (lncRNAs). The genomic instability-derived lncRNA landscape of prostate cancer (PCa) and its critical clinical implications remain to be understood. Patients diagnosed with PCa were recruited from The Cancer Genome Atlas (TCGA) program. Genomic instability-associated lncRNAs were identified by a mutator hypothesis-originated calculative approach. A signature (GILncSig) was derived from genomic instability-associated lncRNAs to classify PCa patients into high-risk and low-risk groups. The biochemical recurrence (BCR) model of a genomic instability-derived lncRNA signature (GILncSig) was established by Cox regression and stratified analysis in the train set. Then its prognostic value and association with clinical features were verified by Kaplan-Meier (K-M) analysis and receiver operating characteristic (ROC) curve in the test set and the total patient set. The regulatory network of transcription factors (TFs) and lncRNAs was established to evaluate TF-lncRNA interactions. A total of 95 genomic instability-associated lncRNAs of PCa were identified. We constructed the GILncSig based on 10 lncRNAs with independent prognostic value. GILncSig separated patients into the high-risk (n = 121) group and the low-risk (n = 121) group in the train set. Patients with high GILncSig score suffered from more frequent BCR than those with low GILncSig score. The results were further validated in the test set, the whole TCGA cohort, and different subgroups stratified by age and Gleason score (GS). A high GILncSig risk score was significantly associated with a high mutation burden and a low critical gene expression (PTEN and CDK12) in PCa. The predictive performance of our BCR model based on GILncSig outperformed other existing BCR models of PCa based on lncRNAs. The GILncSig also showed a remarkable ability to predict BCR in the subgroup of patients with TP53 mutation or wild type. Transcription factors, such as FOXA1, JUND, and SRF, were found to participate in the regulation of lncRNAs with prognostic value. In summary, we developed a prognostic signature of BCR based on genomic instability-associated lncRNAs for PCa, which may provide new insights into the epigenetic mechanism of BCR.

Tang L ，Li W ，Xu H ，Zheng X ，Qiu S ，He W ，Wei Q ，Ai J ，Yang L ，Liu J ... -  《-》

Comprehensive Analyses of Mutation-Derived Long-Chain Noncoding RNA Signatures of Genome Instability in Kidney Renal Papillary Cell Carcinoma.

Background: The role of long-chain noncoding RNA (lncRNA) in genomic instability has been demonstrated to be increasingly importance. Therefore, in this study, lncRNAs associated with genomic instability were identified and kidney renal papillary cell carcinoma (KIRP)-associated predictive features were analysed to classify high-risk patients and improve individualised treatment. Methods: The training (n = 142) and test (n = 144) sets were created using raw RNA-seq and patient's clinical data of KIRP obtained from The Cancer Genome Atlas (TCGA).There are 27 long-chain noncoding RNAs (lncRNAs) that are connected with genomic instability, these lncRNAs were identified using the 'limma' R package based on the numbers of somatic mutations and lncRNA expression profiles acquired from KIRP TCGA cohort. Furthermore, Cox regression analysis was carried out to develop a genome instability-derived lncRNA-based gene signature (GILncSig), whose prognostic value was confirmed in the test cohort as well as across the entire KIRP TCGA dataset. Results: A GILncSig derived from three lncRNAs (BOLA3-AS1, AC004870, and LINC00839), which were related with poor KIRP survival, was identified, which was split up into high- and low-risk groups. Additionally, the GILncSig was found to be an independent prognostic predictive index in KIRP using univariate and multivariate Cox analysis. Furthermore, the prognostic significance and characteristics of GilncSig were confirmed in the training test and TCGA sets. GilncSig also showed better predictive performance than other prognostic lncRNA features. Conclusion: The function of lncRNAs in genomic instability and the genetic diversity of KIRP were elucidated in this work. Moreover, three lncRNAs were screened for prediction of the outcome of KIRP survival and novel insights into identifying cancer biomarkers related to genomic instability were discussed.

Li J ，Wei S ，Zhang Y ，Lu S ，Zhang X ，Wang Q ，Yan J ，Yang S ，Chen L ，Liu Y ，Huang Z ... -  《-》