Immune-related long noncoding RNA signature for predicting survival and immune checkpoint blockade in hepatocellular carcinoma

Long noncoding RNAs (lncRNAs) show multiple functions, including immune response. Recently, the immune-related lncRNAs have been reported in some cancers. We first investigated the immune-related lncRNA signature as a potential target in hepatocellular carcinoma (HCC) survival. The training set (n = 368) and the independent external validation cohort (n = 115) were used. Immune genes and lncRNAs coexpression were constructed to identify immune-related lncRNAs. Cox regression analyses were perfumed to establish the immune-related lncRNA signature. Regulatory roles of this signature on cancer pathways and the immunologic features were investigated. The correlation between immune checkpoint inhibitors and this signature was examined. In this study, the immune-related lncRNA signature was identified in HCC, which could stratify patients into high- and low-risk groups. This immune-related lncRNA signature was correlated with disease progression and worse survival and was an independent prognostic biomarker. Our immune-related lncRNA signature was still a powerful tool in predicting survival in each stratum of age, gender, and tumor stage. This signature mediated cell cycle, glycolysis, DNA repair, mammalian target of rapamycin signaling, and immunologic characteristics (i.e., natural killer cells vs. Th1 cells down, etc). This signature was associated with immune cell infiltration (i.e., macrophages M0, Tregs, CD4 memory T cells, and macrophages M1, etc.,) and immune checkpoint blockade (ICB) immunotherapy-related molecules (i.e., PD-L1, PD-L2, and IDO1). Our findings suggested that the immune-related lncRNA signature had an important value for survival prediction and may have the potential to measure the response to ICB immunotherapy. This signature may guide the selection of the immunotherapy for HCC.     

Identification of 4-lncRNA prognostic signature in head and neck squamous cell carcinoma

Deregulated long noncoding RNAs (lncRNA) have been critically implicated in tumorigenesis and serve as novel diagnostic and prognostic biomarkers. Here we sought to develop a prognostic lncRNA signature in patients with head and neck squamous cell carcinoma (HNSCC). Original RNA-seq data of 499 HNSCC samples were retrieved from The Cancer Genome Atlas database, which was randomly divided into training and testing set. Univariate Cox regression survival analysis, robust likelihood-based survival model and random sampling iterations were applied to identify prognostic lncRNA candidates in the training cohort. A prognostic risk score was developed based on the Cox coefficient of four individual lncRNA imputed as follows: (0.14546 × expression level of RP11-366H4.1) + (0.27106 × expression level of LINC01123) + (0.54316 × expression level of RP11-110I1.14) + (−0.48794 × expression level of CTD-2506J14.1). Kaplan-Meier analysis revealed that patients with high-risk score had significantly reduced overall survival as compared with those with low-risk score when patients in training, testing, and validation cohorts were stratified into high- or low-risk subgroups. Multivariate survival analysis further revealed that this 4-lncRNA signature was a novel and important prognostic factor independent of multiple clinicopathological parameters. Importantly, ROC analyses indicated that predictive accuracy and sensitivity of this 4-lncRNA signature outperformed those previously well-established prognostic factors. Noticeably, prognostic score based on quantification of these 4-lncRNA via qRT-PCR in another independent HNSCC cohort robustly stratified patients into subgroups with high or low survival. Taken together, we developed a robust 4-lncRNA prognostic signature for HNSCC that might provide a novel powerful prognostic biomarker for precision oncology.     

A panel of eight autophagy-related long non-coding RNAs is a good predictive parameter for clear cell renal cell carcinoma

Clear-cell renal cell carcinoma (ccRCC) carries a variable prognosis. Prognostic biomarkers can stratify patients according to risk, and can provide crucial information for clinical decision-making. We screened for an autophagy-related long non-coding lncRNA (lncRNA) signature to improve postoperative risk stratification in The Cancer Genome Atlas (TCGA) database. We confirmed this model in ICGC and SYSU cohorts as a significant and independent prognostic signature. Western blotting, autophagic-flux assay and transmission electron microscopy were used to verify that regulation of expression of 8 lncRNAs related to autophagy affected changes in autophagic flow in vitro. Our data suggest that 8-lncRNA signature related to autophagy is a promising prognostic tool in predicting the survival of patients with ccRCC. Combination of this signature with clinical and pathologic parameters could aid accurate risk assessment to guide clinical management, and this 8-lncRNAs signature related to autophagy may serve as a therapeutic target.     

Identification of m6A methyltransferase-related lncRNA signature for predicting immunotherapy and prognosis in patients with hepatocellular carcinoma

N6-methyladenosine (m6A) methyltransferase has been shown to be an oncogene in a variety of cancers. Nevertheless, the relationship between the long non-coding RNAs (lncRNAs) and hepatocellular carcinoma (HCC) remains elusive. We integrated the gene expression data of 371 HCC and 50 normal tissues from The Cancer Genome Atlas (TCGA) database. Differentially expressed protein-coding genes (DE-PCGs)/lncRNAs (DE-lncRs) analysis and univariate regression and Kaplan–Meier (K–M) analysis were performed to identify m6A methyltransferase-related lncRNAs. Three prognostic lncRNAs were selected by univariate and LASSO Cox regression analyses to construct the m6A methyltransferase-related lncRNA signature. Multivariate Cox regression analyses illustrated that this signature was an independent prognostic factor for overall survival (OS) prediction. The Gene Set Enrichment Analysis (GSEA) suggested that the m6A methyltransferase-related lncRNAs were involved in the immune-related biological processes (BPs) and pathways. Besides, we discovered that the lncRNAs signature was correlated with the tumor microenvironment (TME) and the expression of critical immune checkpoints. Tumor Immune Dysfunction and Exclusion (TIDE) analysis revealed that the lncRNAs could predict the clinical response to immunotherapy. Our study had originated a prognostic signature for HCC based on the potential prognostic m6A methyltransferase-related lncRNAs. The present study had deepened the understanding of the TME status of HCC patients and laid a theoretical foundation for the choice of immunotherapy.     

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 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.     

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.     

Identification of an immune-related RNA-binding protein signature to predict survival and targeted therapy responses in liver cancer

RNA-binding proteins (RBPs) play crucial roles in multiple cancers. However, very few RBPs and their association with immune genes have been systematically studied in liver cancer (LC). We aimed to identify an immune-related RBP signature to predict the survival of LC patients. Bioinformatics methods were used to identify differentially expressed, immune-related, and prognostic RBPs and to develop an immune-related RBP signature based on data from the Cancer Genome Atlas (TCGA) cohort. We obtained eight differentially expressed, immune-related, and prognostic RBPs to construct a risk signature. The signature could effectively distinguish between high- and low-risk patients, and its predictive capacity was validated in the International Cancer Genomics Consortium (ICGC) cohort. We speculated that the high-risk group was more sensitive to targeted therapy. The immune-related RBP signature is an independent prognostic biomarker for LC patients and can expand the application of targeted therapy through patient stratification.     

Identification of a glycolysis-related lncRNA prognostic signature for clear cell renal cell carcinoma

Background: The present study investigated the independent prognostic value of glycolysis-related long noncoding (lnc)RNAs in clear cell renal cell carcinoma (ccRCC).Methods: A coexpression analysis of glycolysis-related mRNAs–long noncoding RNAs (lncRNAs) in ccRCC from The Cancer Genome Atlas (TCGA) was carried out. Clinical samples were randomly divided into training and validation sets. Univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses were performed to establish a glycolysis risk model with prognostic value for ccRCC, which was validated in the training and validation sets and in the whole cohort by Kaplan–Meier, univariate and multivariate Cox regression, and receiver operating characteristic (ROC) curve analyses. Principal component analysis (PCA) and functional annotation by gene set enrichment analysis (GSEA) were performed to evaluate the risk model.Results: We identified 297 glycolysis-associated lncRNAs in ccRCC; of these, 7 were found to have prognostic value in ccRCC patients by Kaplan–Meier, univariate and multivariate Cox regression, and ROC curve analyses. The results of the GSEA suggested a close association between the 7-lncRNA signature and glycolysis-related biological processes and pathways.Conclusion: The seven identified glycolysis-related lncRNAs constitute an lncRNA signature with prognostic value for ccRCC and provide potential therapeutic targets for the treatment of ccRCC patients.     

Identification of a universal 6-lncRNA prognostic signature for three pathologic subtypes of renal cell carcinoma

Renal cell carcinoma (RCC) is the most common adult renal epithelial cancer susceptible to metastasis and patients with irresectable RCC always have a poor prognosis. Long noncoding RNAs (lncRNAs) have recently been documented as having critical roles in the etiology of RCC. Nevertheless, the prognostic significance of lncRNA-based signature for outcome prediction in patients with RCC has not been well investigated. Therefore, it is essential to identify a lncRNA-based signature for predicting RCC prognosis. In the current study, we comprehensively analyzed the RNA sequencing data of the three main pathological subtypes of RCC (kidney renal clear cell carcinoma [KIRC], kidney renal papillary cell carcinoma [KIRP], and kidney chromophobe carcinoma [KICH]) from The Cancer Genome Atlas (TCGA) database, and identified a 6-lncRNA prognostic signature with the help of a step-wise multivariate Cox regression model. The 6-lncRNA signature stratified the patients into low- and high-risk groups with significantly different prognosis. Multivariate Cox regression analysis showed that predictive value of the 6-lncRNA signature was independent of other clinical or pathological factors in the entire cohort and in each cohort of RCC subtypes. In addition, the three independent prognostic clinical factors (including age, pathologic stage III, and stage IV) was also stratified into low- and high-risk groups basis on the risk score, and the stratification analyses demonstrated that the high-risk score was a poor prognostic factor. In conclusion, these findings indicate that the 6-lncRNA signature is a novel prognostic biomarker for all three subtypes of RCC, and can increase the accuracy of predicting overall survival.     

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.     

Survival prediction and response to immune checkpoint inhibitors: A prognostic immune signature for hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most common cancers all over the world. Several studies have explored if immune-related genes and tumor immune microenvironment could play roles in HCC prognoses. This study is aimed at developing a prognostic signature of HCC based on immune-related genes or tumor immune microenvironment to predict survival and response to immune checkpoint inhibitors (ICIs). We constructed a prognostic signature using bioinformatics method and validated its predictive capability. The mechanisms of the signature prediction were explored with The Cancer Immunome Atlas (TCIA) and mutation analysis. We also explored the association between the signature and immunophenoscore (IPS), which is the marker of ICIs response. A 6 immune-related-gene (6-IRG) signature was developed. It was revealed in a multivariate analysis that the 6-IRG signature was an independent prognostic factor of overall survival and progression-free interval among HCC patients. In the high-risk group of 6-IRG signature score, macrophage M0 cells and regulatory T cells, which are observed associated with poor overall survival in our study, were higher. The low-risk group had a higher IPS, which meant a better response to ICIs. Taken together, we constructed a reliable 6-IRG signature for prediction of survival and response to ICIs. The signature needs further testing for clinical application.     

Exploration of a novel prognostic risk signatures and immune checkpoint molecules in endometrial carcinoma microenvironment

In this study, we devoted to investigate immune-related genes and tumor microenvironment (TME) in EC based on The Cancer Genome Atlas (TCGA) database. In total 799 up-regulated and 139 down-regulated immune-related and differentially expressed genes in EC were investigated for functional annotations and prognosis. By a conjoint Cox regression analysis, we built two risk models for OS and DFS, as well as the consistent nomograms. Immune-related pathways were revealed mostly enriched in the low-risk group. By further analyzing TME based on the risk signatures, the higher immune cell infiltration and activation, lower tumor purity and higher tumor mutational burden were found in low-risk group, which presented a better prognosis. Both the expression and immunophenoscore of immune checkpoints PD-1, CTLA4, PD-L1 and PD-L2 increased significantly in low-risk group. These findings may provide new ideas for novel biomarkers and immunotherapy targets in EC.     

Development and validation of a five-lncRNA signature with prognostic value in colon cancer

Dysregulation of long noncoding RNAs (lncRNAs) has been found in a large number of human cancers, including colon cancer. Therefore, the implementation of potential lncRNAs biomarkers with prognostic prediction value are very much essential. GSE39582 data set was downloaded from database of Gene Expression Omnibus. Re-annotation analysis of lncRNA expression profiles was performed by NetAffx annotation files. Univariate and multivariate Cox proportional analyses helped select prognostic lncRNAs. Algorithm of random survival forest-variable hunting (RSF-VH) together with stepwise multivariate Cox proportional analysis were performed to establish lncRNA signature. The log-rank test was carried out to analyze and compare the Kaplan-Meier survival curves of patients’ overall survival (OS). Receiver operating characteristic (ROC) analysis was used for comparing the survival prediction regarding its specificity and sensitivity based on lncRNA risk score, followed by calculating the values of area under the curve (AUC). The single-sample GSEA (ssGSEA) analysis was used to describe biological functions associated with this signature. Finally, to determine the robustness of this model, we used the validation sets including GSE17536 and The Cancer Genome Atlas data set. After re-annotation analysis of lncRNAs, a total of 14 lncRNA probes were obtained by univariate and multivariate Cox proportional analysis. Then, the RSF-VH algorithm and stepwise multivariate Cox analysis helped to build a five-lncRNA prognostic signature for colon cancer. The patients in group with high risk showed an obviously shorter survival time compared with patients in group with low risk with AUC of 0.75. In addition, the five-lncRNA signature can be used to independently predict the survival of patients with colon cancer. The ssGSEA analysis revealed that pathways such as extracellular matrix-receptor interaction was activated with an increase in risk score. These findings determined the strong power of prognostic prediction value of this five-lncRNA signature for colon cancer.     

Developing a risk scoring system based on immune-related lncRNAs for patients with gastric cancer

The immune system and the tumor interact closely during tumor development. Aberrantly expressed long non-coding RNAs (lncRNAs) may be potentially applied as diagnostic and prognostic markers for gastric cancer (GC). At present, the diagnosis and treatment of GC patients remain a formidable clinical challenge. The present study aimed to build a risk scoring system to improve the prognosis of GC patients. In the present study, ssGSEA was used to evaluate the infiltration of immune cells in GC tumor tissue samples, and the samples were split into a high immune cell infiltration group and a low immune cell infiltration group. About 1262 differentially expressed lncRNAs between the high immune cell infiltration group and the low immune cell infiltration group. About 3204 differentially expressed lncRNAs between GC tumor tissues and paracancerous tissues were identified. Then, 621 immune-related lncRNAs were screened using a Venn analysis based on the above results, and 85 prognostic lncRNAs were identified using a univariate Cox analysis. We constructed a prognostic signature using LASSO analysis and evaluated the predictive performance of the signature using ROC analysis. GO and KEGG enrichment analyses were performed on the lncRNAs using the R package, ‘clusterProfiler’. The TIMER online database was used to analyze correlations between the risk score and the abundances of the six types of immune cells. In conclusion, our study found that specific immune-related lncRNAs were clinically significant. These lncRNAs were used to construct a reliable prognostic signature and analyzed immune infiltrates, which may assist clinicians in developing individualized treatment strategies for GC patients.     

A prognostic eight-lncRNA expression signature in predicting recurrence of ER-positive breast cancer receiving endocrine therapy

Long noncoding RNAs (lncRNAs) have the main role in the tumorigenesis of breast cancer. In the present study, lncRNA expression profiling was collected to identify a lncRNA expression signature from the Gene Expression Omnibus database. An eight-lncRNA signature was established to predict the survival of patients with estrogen receptor (ER)-positive breast cancer receiving endocrine therapy. Patients were separated into a low-risk group and a high-risk group based on this signature. Patients in high-risk group have worse survival compared to those in low-risk group using Kaplan–Meier curve analysis with log-rank test. Receiver operating characteristic analysis suggested good diagnostic efficiency of the eight-lncRNA signature. When adjusting the clinical features, including age, grade, lymph node status, and tumor size, this signature was independently associated with the relapse-free survival. The prognostic value of the lncRNA prognostic model was then validated in validation sets. When validated in a cohort of patients treated with neoadjuvant chemotherapy and endocrine therapy, this signature demonstrated good performance as well. Besides, we have built a nomogram that integrated the conventional clinicopathological features and the eight-lncRNA-based signature. To sum up, our results indicated that the eight-lncRNA prognostic model was a reliable tool to group patients at high and low risk of disease relapse. This signature may have possible implication in prognostic evaluations of patients with ER-positive breast cancer receiving endocrine therapy.     

Identification of a prognostic signature based on immune-related genes in bladder cancer

Bladder cancer (BLCA) has a high incidence and recurrence rate, and the effect of immunotherapy varies from person to person. Immune-related genes (IRGs) have been shown to be associated with immunotherapy and prognosis in many other cancers, but their role in immunogenic BLCA is less well defined. In this study, we constructed an eight-IRG risk model, which demonstrated strong prognostic and immunotherapeutic predictive power. The signature was significantly related to tumor clinicopathological characteristics, tumor class, immune cell infiltration and mutation status. Additionally, a nomogram containing the risk score and other potential risk factors could effectively predict the long-term overall survival probability of BLCA patients. The enriched mechanisms identified by gene set enrichment analysis suggested that the reason why this signature can accurately distinguish high- and low-risk populations may be closely related to the different degrees of innate immune response and T cell activation in different patients.