An eight-long noncoding RNA expression signature for colorectal cancer patients’ prognosis

Long noncoding RNAs (lncRNAs) have recently emerged as important biomarkers of cancer progression. Here, we proposed to develop a lncRNA-based signature with a prognostic value for colorectal cancer (CRC) overall survival (OS). Through mining microarray datasets, we analyzed the lncRNA expression profiles of 122 patients with CRC from Gene Expression Omnibus. Associations between lncRNA and CRC OS were firstly evaluated through univariate Cox regression analysis. A random survival forest method was applied for further screening of the lncRNA signature, which resulted in eight lncRNAs, including PEG3-AS1, LOC100505715, MINCR, DBH-AS1, LINC00664, FAM224A, LOC642852, and LINC00662. Combination of the eight lncRNAs weighted by their multivariate Cox regression coefficients formed a prognostic signature, through which, we could divide the 122 patients with CRC into two subgroups with significantly different OS. Good robustness of the lncRNA signature's prognostic value was verified through an independent data set consisting of 55 patients with CRC. In addition, gene set enrichment analysis indicated the potential association between high prognostic value and oxygen metabolism-related processes. This result should indicate that lncRNAs could be a useful signature for CRC prognosis.     

Genome-wide screening and cohorts validation identifying novel lncRNAs as prognostic biomarkers for clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the main subtype of renal cell carcinoma with varied prognosis. We aimed to identify and assess the possible prognostic long noncoding RNA (lncRNA) biomarkers. LncRNAs expression data and corresponding clinical information of 619 ccRCC patients were downloaded from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. Differentially expressed genes analysis, univariate Cox regression, the least absolute shrinkage and selection operator Cox regression model were utilized to identify hub lncRNAs. Multivariate Cox regression was used to establish the risk model. Statistical analysis was performed using R 3.5.3. The expression value of five lncRNAs and the risk-score levels were significantly associated with a survival prognosis of ccRCC patients (all P < .001). In the TCGA validation cohort, the area under the curve (AUC) for the integrated nomogram was 0.905 and 0.91 for 3-, 5-year prediction separately. The AUC reached up to 0.757 in an independent ICGC cohort. Besides, the calibration plots also illustrated well curve-fitting between observation values and predictive values. Weighted gene co-expression network analysis and subsequent pathway analysis revealed that the PI3K-Akt-mTOR and hypoxia-inducible factor signaling crosstalk might function as the most essential mechanisms related to the five-lncRNAs signature. Our study suggested that lncRNA AC009654.1, AC092490.2, LINC00524, LINC01234, and LINC01885 were significantly associated with ccRCC prognosis. The prognostic model based on this five lncRNA may predict the overall survival of ccRCC.     

An Immune-Related Six-lncRNA Signature to Improve Prognosis Prediction of Glioblastoma Multiforme

Recent studies have demonstrated the utility and superiority of long non-coding RNAs (lncRNAs) as novel biomarkers for cancer diagnosis, prognosis, and therapy. In the present study, the prognostic value of lncRNAs in glioblastoma multiforme was systematically investigated by performing a genome-wide analysis of lncRNA expression profiles in 419 glioblastoma patients from The Cancer Genome Atlas (TCGA) project. Using survival analysis and Cox regression model, we identified a set of six lncRNAs (AC005013.5, UBE2R2-AS1, ENTPD1-AS1, RP11-89C21.2, AC073115.6, and XLOC_004803) demonstrating an ability to stratify patients into high- and low-risk groups with significantly different survival (median 0.899 vs. 1.611 years, p = 3.87e?09, log-rank test) in the training cohort. The six-lncRNA signature was successfully validated on independent test cohort of 219 patients with glioblastoma, and it revealed superior performance for risk stratification with respect to existing lncRNA-related signatures. Multivariate Cox and stratification analysis indicated that the six-lncRNA signature was an independent prognostic factor after adjusting for other clinical covariates. Further in silico functional analysis suggested that the six-lncRNA signature may be involved in the immune-related biological processes and pathways which are very well known in the context of glioblastoma tumorigenesis. The identified lncRNA signature had important clinical implication for improving outcome prediction and guiding the tailored therapy for glioblastoma patients with further prospective validation.     

Identification of pathological grade and prognosis-associated lncRNA for ovarian cancer

Ovarian carcinoma (OC) is one of the most common malignant tumors in female genitals. In recent years, the therapeutic effect of OC has been significantly improved through the application of effective chemotherapy regimen. However, the 5-year survival rate is also lower than 30% due to high rate of relapse. So, it is needed to screen reliable predictive and prognostic markers of OC. Ovarian cancer gene expression data and corresponding clinical data used were downloaded from Gene Expression Omnibus database. Weighted gene expression network analysis (WGCNA) and Cox proportional hazards regression (PHR) were used to screen Pathological Grade and Prognosis-associated long noncoding RNA (lncRNA). Kaplan-Meier analysis and receiver operating characteristic curves analysis were performed to evaluate the predictive ability of the selected lncRNA. Gene Ontology (GO) enrichment and Gene Set Enrichment Analysis (GSEA) enrichment analysis methods were used to explore the possible mechanisms of the selected lncRNA affecting the development of OC. Five reliably lncRNAs (LINC00664, LINC00667, LINC01139, LINC01419, and LOC286437) was identified through a series of bioinformatics methods. In testing cohorts, we found that the five lncRNAs in predicting the risk of OC recurrence is robustness, and multivariate Cox PHR analysis indicate that the five lncRNAs is an independent risk factor for OC recurrence. Moreover, GO and GSEA enrichment analysis showed that the five lncRNAs are involved in multiple ovarian cancer occurrence mechanism. In summary, all these findings indicated that the five lncRNAs can effectively predict the risk of recurrence of ovarian cancer.     

A prognostic 11 long noncoding RNA expression signature for breast invasive carcinoma

Breast cancer, the most common cancer in women worldwide, is associated with high mortality. The long non-coding RNAs (lncRNAs) with a little capacity of coding proteins is playing an increasingly important role in the cancer paradigm. Accumulating evidences demonstrate that lncRNAs have crucial connections with breast cancer prognosis while the studies of lncRNAs in breast cancer are still in its primary stage. In this study, we collected 1052 clinical patient samples, a comparatively large sample size, including 13 159 lncRNA expression profiles of breast invasive carcinoma (BRCA) from The Cancer Genome Atlas database to identify prognosis-related lncRNAs. We randomly separated all of these clinical patient samples into training and testing sets. In the training set, we performed univariable Cox regression analysis for primary screening and played the model for Robust likelihood-based survival for 1000 times. Then 11 lncRNAs with a frequency more than 600 were selected for prediction of the prognosis of BRCA. Using the analysis of multivariate Cox regression, we established a signature risk-score formula for 11 lncRNA to identify the relationship between lncRNA signatures and overall survival. The 11 lncRNA signature was validated both in the testing and the complete set and could effectively classify the high-/low-risk group with different OS. We also verified our results in different stages. Moreover, we analyzed the connection between the 11 lncRNAs and the genes of ESR1, PGR, and Her2, of which protein products (ESR, PGR, and HER2) were used to classify the breast cancer subtypes widely. The results indicated correlations between 11 lncRNAs and the gene of PGR and ESR1. Thus, a prognostic model for 11 lncRNA expression was developed to classify the BRAC clinical patient samples, providing new avenues in understanding the potential therapeutic methods of breast cancer.     

Identification of LINC01615 as potential metastasis-related long noncoding RNA in hepatocellular carcinoma

Hepatocellular carcinoma is one of the most prevalent and fatal cancers. Studying the long noncoding RNA (lncRNA) alterations in hepatocellular carcinoma may lead to new therapeutic strategies. We checked whether there were correlations between The Cancer Genome Atlas expression profiles of the differentially expressed lncRNAs and their DNA methylation status or the copy number variations for hepatocellular carcinoma. We obtained 41 lncRNAs that were differentially expressed between tumor and normal samples, and their DNA methylation status was negatively correlated with the expression levels. We identified five lncRNAs that were recurrently amplified or deleted in tumor samples, but none of them were associated with the messenger RNA (mRNA) expression levels. To obtain the biological function of these lncRNAs, the coexpressed mRNAs in the hepatocellular carcinoma were figured out. A total of 10 lncRNAs were highly correlated with at least one gene. Six out of the ten lncRNAs were already known to be related with cancer previously. LINC01615 had 72 coexpressed genes, and we carried out the gene ontology (GO) term enrichment for these protein-coding genes. The results suggested that these lncRNAs were associated with extracellular matrix organization. To summarize, we identified 41 potentially cancer-related lncRNAs. In particular, we proposed that LINC01615 potentially affected the extracellular matrix and had further impacts on the metastasis of hepatocellular carcinoma.     

Identification of autophagy-related long non-coding RNA prognostic signature for breast cancer

Autophagy-related long non-coding RNAs (lncRNAs) disorders are related to the occurrence and development of breast cancer. The purpose of this study is to explore whether autophagy-related lncRNA can predict the prognosis of breast cancer patients. The autophagy-related lncRNAs prognostic signature was constructed by Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression. We identified five autophagy-related lncRNAs (MAPT-AS1, LINC01871, AL122010.1, AC090912.1, AC061992.1) associated with prognostic value, and they were used to construct an autophagy-related lncRNA prognostic signature (ALPS) model. ALPS model offered an independent prognostic value (HR = 1.664, 1.381-2.006), where this risk score of the model was significantly related to the TNM stage, ER, PR and HER2 status in breast cancer patients. Nomogram could be utilized to predict survival for patients with breast cancer. Principal component analysis and Sankey Diagram results indicated that the distribution of five lncRNAs from the ALPS model tends to be low-risk. Gene set enrichment analysis showed that the high-risk group was enriched in autophagy and cancer-related pathways, and the low-risk group was enriched in regulatory immune-related pathways. These results indicated that the ALPS model composed of five autophagy-related lncRNAs could predict the prognosis of breast cancer patients.     

Comprehensive analysis of novel three-long noncoding RNA signatures as a diagnostic and prognostic biomarkers of human triple-negative breast cancer

Currently, traditional predictors of prognosis (tumor size, nodal status, progesterone receptor [PR], estrogen receptor [ER], or human epidermal growth factor receptor-2 [HER2]) are insufficient for precise survival prediction for triple-negative breast cancer (TNBC). Long noncoding RNAs (lncRNAs) have been observed to exert critical functions in cancer, including in TNBC. Nevertheless, systematically tracking expression-based lncRNA biomarkers based on the sequence data for the prediction of prognosis in TNBC has not yet been investigated. To ascertain whether biomarkers exist that can distinguish TNBC from adjacent normal tissue or nTNBC, we implemented a comprehensive analysis of lncRNA expression profiles and clinical data of 1097 BC samples from The Cancer Genome Atlas database. A total of 1510 differentially expressed lncRNAs in normal and TNBC samples were extracted. Similarly, 672 differentially expressed lncRNAs between nTNBC and TNBC samples were detected. The receiver operating characteristic curve analysis indicated that three upregulated lncRNAs (AC091043.1, AP000924.1, and FOXCUT) may be of strong diagnostic value for predicting the existence of TNBC in the training and validation sets (area under the curve (AUC > 0.85). Kaplan-Meier analysis demonstrated that the other three lncRNAs (AC010343.3, AL354793.1, and FGF10-AS1) were associated with the prognosis of TNBC patients (P < 0.05). We used the three overall survival (OS)-related lncRNAs to establish a three-lncRNA signature. Multivariate Cox regression analysis suggested that the three-lncRNA signature was a prognostic factor independent of other clinical variables ( P < 0.01) for predicting OS in TNBC patients that could be utilized to classify patients into high- or low-risk subgroups. Our results might provide efficient signatures for clinical diagnosis and prognostic evaluation of TNBC.     

Downregulation of long noncoding RNA LINC00683 associated with unfavorable prognosis in prostate cancer based on TCGA

Prostate cancer (PCa) is the third most common reason of cancer-related deaths in men. Accumulating evidence has shown that dysregulation of long noncoding RNAs (lncRNAs) is closely related to cancer initiation and development. Although large numbers of lncRNAs have been discovered, knowledge regarding their function and physiological/pathological significance remains limited. In this study, we aimed to reveal functional lncRNAs and identify prognosis-related RNAs in PCa by analyzing data from The Cancer Genome Atlas (TCGA). To achieve this, an lncRNA-mRNA coexpression network was constructed by weighted correlation network analysis. Additionally, a subnetwork was extracted from this weighted correlation network, and seven lncRNAs were identified as core nodes. Further Kaplan-Meier survival analysis showed that three lncRNAs (LINC00683, LINC00857, and FENDRR) were significantly downregulated in PCa samples, and there was a strong positive correlation with patient survival. Importantly, LINC00683 has not been fully reported as related with PCa. Additionally, gene set enrichment analysis indicated that LINC00683 might be involved in cancer-related pathways such as the Wnt pathway. Based on the findings of this study, lncRNA LINC00683 is likely to provide a new diagnostic biomarker and therapeutic target for future PCa treatments.     

基于TCGA和GEO数据库的胃癌lncRNA筛选与ceRNA网络构建

胃癌（gastric cancer,GC）是我国最常见的恶性肿瘤之一,严重危害人类健康。胃癌发病机制复杂,缺乏特异性预后生物标志物。长链非编码RNA（long non?coding RNA,lncRNA）可作为竞争性内源RNA（competing endogenous RNA,ceRNA）,影响microRNA（miRNA）与mRNA的结合,从而影响胃癌的发生、发展。基于TCGA和GEO数据库的转录组数据,筛选GC中差异表达的lncRNAs,并构建基于6条lncRNAs（HAGLROS、TMEM92?AS1、LINC01745、HOXC?AS3、SEMA3B?AS1、FEZF1?AS1）的lncRNA?miRNA?mRNA网络。网络核心基因的KEGG/GO富集和蛋白质互作分析结果显示,lncRNA可能通过miRNA海绵吸附作用,调控胃癌的发生、发展与转移。AC011352.1、AC087636.1、AC093627.1、GAS1RR与胃癌患者的预后相关性具有统计学意义（P<0.05）,并可能成为胃癌患者潜在的预后生物标志物。     

Integrated analysis of a competing endogenous RNA network reveals key long noncoding RNAs as potential prognostic biomarkers for hepatocellular carcinoma

Growing evidence has revealed that long noncoding RNAs (lncRNAs) have an important impact on tumorigenesis and tumor progression via a mechanism involving competing endogenous RNAs (ceRNAs). However, their use in predicting the survival of a patient with hepatocellular carcinoma (HCC) remains unclear. The aim of this study was to develop a novel lncRNA expression–based risk score system to accurately predict the survival of patients with HCC. In our study, using expression profiles downloaded from The Cancer Genome Atlas database, the differentially expressed messenger RNAs (mRNAs), lncRNAs, and microRNAs (miRNAs) were explored in patients with HCC and normal liver tissues, and then a ceRNA network constructed. A risk score system was established between lncRNA expression of the ceRNA network and overall survival (OS) or recurrence-free survival (RFS); it was further analyzed for associations with the clinical features of patients with HCC. In HCC, 473 differentially expressed lncRNAs, 63 differentially expressed miRNAs, and 1417 differentially expressed mRNAs were detected. The ceRNA network comprised 41 lncRNA nodes, 12 miRNA nodes, 24 mRNA nodes, and 172 edges. The lncRNA expression–based risk score system for OS was constructed based on six lncRNAs (MYLK-AS1, AL359878.1, PART1, TSPEAR-AS1, C10orf91, and LINC00501), while the risk score system for RFS was based on four lncRNAs (WARS2-IT1, AL359878.1, AL357060.1, and PART1). Univariate and multivariate Cox analyses showed the risk score systems for OS or RFS were significant independent factors adjusted for clinical factors. Receiver operating characteristic curve analysis showed the area under the curve for the risk score system was 0.704 for OS, and 0.71 for RFS. Our result revealed a lncRNA expression–based risk score system for OS or RFS can effectively predict the survival of patients with HCC and aid in good clinical decision-making.     

Construction of an immune-related LncRNA signature with prognostic significance for bladder cancer

Bladder cancer (BLCA) is one of the most common urological cancer with increasing cases and deaths every year. In the present study, we aim to construct an immune-related prognostic lncRNA signature (IRPLS) in bladder cancer (BLCA) patients and explore its immunogenomic implications in pan-cancers. First, the immune-related differentially expressed lncRNAs (IRDELs) were identified by ‘limma’ R package and the score of IRPLS in every patient were evaluated by Cox regression. The dysregulation of IRDELs expression between cancer and para-cancer normal tissues was validated through RT-qPCR. Then, we further explore the biological functions of a novel lncRNA from IRPLS, RP11-89 in BLCA using CCK8 assay, Transwell assay and Apoptosis analysis, which indicated that RP11-89 was able to promote cell proliferation and invasive capacity while inhibits cell apoptosis in BLCA. In addition, we performed bioinformatic methods and RIP to investigate and validate the RP11-89/miR-27a-3p/PPARγ pathway in order to explore the mechanism. Next, CIBERSORT and ESTIMATE algorithm were used to evaluate abundance of tumour-infiltrating immune cells and scores of tumour environment elements in BLCA with different level of IRPLS risk scores. Finally, multiple bioinformatic methods were performed to show us the immune landscape of these four lncRNAs for pan-cancers. In conclusion, this study first constructed an immune-related prognostic lncRNA signature, which consists of RP11-89, PSORS1C3, LINC02672 and MIR100HG and might shed lights on novel targets for individualized immunotherapy for BLCA patients.     

A novel seven-lncRNA signature for prognosis prediction in hepatocellular carcinoma

Liver cancer is still one of the leading causes of cancer-related death worldwide. This study is dedicated to developing a multi–long noncoding RNA (lncRNA) model for risk stratification and prognosis prediction on patients with hepatocellular carcinoma (HCC). We first downloaded lncRNA expression profiles and corresponding clinical information of patients with liver cancer from The Cancer Genome Atlas database. Differentially expressed (DE) lncRNAs between HCC samples and normal samples were identified. In total, 308 patients with HCC were randomly divided into a training group (n = 154) and a testing group (n = 154). Univariate Cox regression and least absolute shrinkage and selection operator Cox regression analyses were performed to select the best survival-related candidates from these DE lncRNAs in the training set. Seven lncRNAs (AC009005.2, RP11-363N22.3, RP11-932O9.10, RP11-572O6.1, RP11-190C22.8, RP11-388C12.8, and ZFPM2-AS1) were finally identified and used to construct a seven-lncRNA signature. The signature could classify patients into high-risk and low-risk groups with significantly different overall survival. The area under the curve of receiver operating characteristic curve for the signature to predict 5-year survival reached more than 0.75. Besides, the prognostic value of the seven-lncRNA signature was independent of conventional clinical factors. The predictive performance of the signature was further validated in the testing set and the whole set. Functional enrichment analysis indicated that the seven prognostic lncRNAs may be involved in several essential biological processes and pathways. The current study demonstrated the potential clinical implications of the seven-lncRNA signature for survival prediction of patients with HCC.     

An angiogenesis‐related long noncoding RNA signature correlates with prognosis in patients with hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most prevalent and lethal cancers worldwide. Neovascularization is closely related to the malignancy of tumors. We constructed a signature of angiogenesis-related long noncoding RNA (lncRNA) to predict the prognosis of patients with HCC. The lncRNA expression matrix of 424 HCC patients was downloaded from The Cancer Genome Atlas (TCGA). First, gene set enrichment analysis (GSEA) was used to distinguish the differentially expressed genes of the angiogenesis genes in liver cancer and adjacent tissues. Next, a signature of angiogenesis-related lncRNAs was constructed using univariate and multivariate analyses, and receiver operating characteristic (ROC) curves were used to assess the accuracy. The signature and relevant clinical information were used to construct the nomogram. A 5-lncRNA signature was highly correlated with overall survival (OS) in HCC patients and performed well in evaluations using the C-index, areas under the curve, and calibration curves. In summary, the 5-lncRNA model can serve as an accurate signature to predict the prognosis of patients with liver cancer, but its mechanism of action must be further elucidated by experiments.     

A 7-lncRNA signature predict prognosis of Uterine corpus endometrial carcinoma

Current research indicate that long noncoding RNAs (lncRNAs) are associated with the progression of various cancers and can be used as prognostic biomarkers. This study aims to construct a prognostic lncRNA signature for the risk assessment of Uterine corpus endometrial carcinoma (UCEC). The RNA-Seq expression profile and corresponding clinical data of UCEC patients obtained from The Cancer Genome Atlas database. First, some prognosis-related lncRNAs were obtained by univariate Cox analysis. The minimum absolute contraction and selection operator (LASSO) regression and the Cox proportional hazard regression method were used to further identify the lncRNA prognostic model. Finally, seven lncRNAs (AC110491.1, AL451137.1, AC005381.1, AC103563.2, AC007422.2, AC108025.2, and MIR7-3HG) were identified as potential prognostic factors. According to the model constructed by the above analysis, the risk score of each UCEC patient was calculated, and the patients were classified into high and low-risk groups. The low-risk group had significant survival benefits. Moreover, we constructed a nomogram that incorporated independent prognostic factors (age, tumor stage, tumor grade, and risk score). The c-index value for evaluating the predictive nomogram model was 0.801. The area under the curve was 0.797 (3-year survival). The calibration curve also showed that there was a satisfactory agreement between the predicted and observed values in the probability of 1-, 3-, and 5-year overall survival. On the basis of the coexpression relationship, we established a coexpression network of lncRNA-messenger RNA (mRNA) of the 7-lncRNA. The Kyoto Encyclopedia of Genes and Genomes analysis of the coexpressing mRNAs showed that the main pathways related to the 7-lncRNA signature were neuroactive ligand-receptor interaction, serotonergic synapse, and gastric cancer pathway. Therefore, our study revealed that the 7-lncRNA could be used to predict the prognosis of UCEC and for postoperative treatment and follow-up.     

A novel autophagy-related lncRNA survival model for lung adenocarcinoma

Long non-coding RNA (lncRNA) is an important regulatory factor in the development of lung adenocarcinoma, which is related to the control of autophagy. LncRNA can also be used as a biomarker of prognosis in patients with lung adenocarcinoma. Therefore, it is important to determine the prognostic value of autophagy-related lncRNA in lung adenocarcinoma. In this study, autophagy-related mRNAs-lncRNAs were screened from lung adenocarcinoma and a co-expression network of autophagy-related mRNAs-lncRNAs was constructed by using The Cancer Genome Atlas (TCGA). The univariate and multivariate Cox proportional hazard analyses were used to evaluate the prognostic value of the autophagy-related lncRNAs and finally obtained a survival model composed of 11 autophagy-related lncRNAs. Through Kaplan-Meier analysis, univariate and multivariate Cox regression analysis and time-dependent receiver operating characteristic (ROC) curve analysis, it was further verified that the survival model was a new independent prognostic factor for patients with lung adenocarcinoma. In addition, based on the survival model, gene set enrichment analysis (GSEA) was used to illustrate the function of genes in low-risk and high-risk groups. These 11 lncRNAs were GAS6-AS1, AC106047.1, AC010980.2, AL034397.3, NKILA, AL606489.1, HLA-DQB1-AS1, LINC01116, LINC01806, FAM83A-AS1 and AC090559.1. The hazard ratio (HR) of the risk score was 1.256 (1.196-1.320) (P < .001) in univariate Cox regression analysis and 1.215 (1.149-1.286) (P < .001) in multivariate Cox regression analysis. And the AUC value of the risk score was 0.809. The 11 autophagy-related lncRNA survival models had important predictive value for the prognosis of lung adenocarcinoma and may become clinical autophagy-related therapeutic targets.