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

Construction of co-expression modules related to survival by WGCNA and identification of potential prognostic biomarkers in glioblastoma

Glioblastoma (GBM) is a malignant brain tumour with poor prognosis. The potential pathogenesis and therapeutic target are still need to be explored. Herein, TCGA expression profile data and clinical information were downloaded, and the WGCNA was conducted. Hub genes which closely related to poor prognosis of GBM were obtained. Further, the relationship between the genes of interest and prognosis of GBM, and immune microenvironment were analysed. Patients from TCGA were divided into high- and low-risk group. WGCNA was applied to the high- and low-risk group and the black module with the lowest preservation was identified which could distinguish the prognosis level of these two groups. The top 10 hub genes which were closely related to poor prognosis of patients were obtained. GO analysis showed the biological process of these genes mainly enriched in: Cell cycle, Progesterone-mediated oocyte maturation and Oocyte meiosis. CDCA5 and CDCA8 were screened out as the genes of interest. We found that their expression levels were closely related to overall survival. The difference analysis resulted from the TCGA database proved both CDCA5 and CDCA8 were highly expressed in GBM. After transfection of U87-MG cells with small interfering RNA, it revealed that knockdown of the CDCA5 and CDCA8 could influence the biological behaviours of proliferation, clonogenicity and apoptosis of GBM cells. Then, single-gene analysis was performed. CDCA5 and CDCA8 both had good correlations with genes that regulate cell cycle in the p53 signalling pathway. Moreover, it revealed that high amplification of CDCA5 was correlated with CD8⁺ T cells while CDCA8 with CD4⁺ T cells in GBM. These results might provide new molecular targets and intervention strategy for GBM.     

LncRNA and mRNA signatures associated with neoadjuvant chemoradiotherapy downstaging effects in rectal cancer

Radiotherapy plays a crucial role in combined treatment modality in local advanced rectal cancer (LARC). While neoadjuvant chemoradiotherapy responses were variable in LARC patients, so, it is important to identify genes that closely associated with short-term and long-term responses to radiotherapy. In this study, we profiled long noncoding RNAs (lncRNAs) and messenger RNAs (mRNAs) expression values of LARC patients with different neoadjuvant chemoradiotherapy downstaging depth score based on Agilent Arraystar Human LncRNA V3.0 Array(Agilent, CA). LncRNAs and mRNAs with aberrant expression values between the two groups of LARC patients were identified and lncRNA-miRNA-mRNA regulation network was also obtained through the combination of miRcode and miRTarBase database. Gene interaction network and module analysis of differential expression mRNAs contained in the lncRNA-miRNA-mRNA network identified five hub genes, including KRAS, PDPK1, PPP2R5C, PPP2R1B, and YES1, that should be closely associated with LARC’s response to chemoradiotherapy. Besides, Kaplan-Meier analysis based on the Cyber Research Center (CRC) data set from The Cancer Genome Atlas indicated that aberrant expression of the five hub genes is significantly associated with CRC overall survival. In conclusion, we obtained several biomarkers that should be associated with neoadjuvant chemoradiotherapy response in LARC, which should be helpful for individual treatment and prognosis improvement.     

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 and development of long non‐coding RNA‐associated regulatory network in colorectal cancer

Colorectal cancer (CRC) is one of the leading causes of cancer‐associated death globally. Long non‐coding RNAs (lncRNAs) have been identified as micro RNA (miRNA) sponges in a competing endogenous RNA (ceRNA) network and are involved in the regulation of mRNA expression. This study aims to construct a lncRNA‐associated ceRNA network and investigate the prognostic biomarkers in CRC. A total of 38 differentially expressed (DE) lncRNAs, 23 DEmiRNAs and 27 DEmRNAs were identified by analysing the expression profiles of CRC obtained from The Cancer Genome Atlas (TCGA). These RNAs were chosen to develop a ceRNA regulatory network of CRC, which comprised 125 edges. Survival analysis showed that four lncRNAs, six miRNAs and five mRNAs were significantly associated with overall survival. A potential regulatory axis of ADAMTS9‐AS2/miR‐32/PHLPP2 was identified from the network. Experimental validation was performed using clinical samples by quantitative real‐time PCR (qRT‐PCR), which showed that expression of the genes in the axis was associated with clinicopathological features and the correlation among them perfectly conformed to the ‘ceRNA theory’. Overexpression of ADAMTS9‐AS2 in colon cancer cell lines significantly inhibited the miR‐32 expression and promoted PHLPP2 expression, while ADAMTS9‐AS2 knockdown had the opposite effects. The constructed novel ceRNA network may provide a comprehensive understanding of the mechanisms of CRC carcinogenesis. The ADAMTS9‐AS2/miR‐32/PHLPP2 regulatory axis may serve as a potential therapeutic target for CRC.     

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.     

An 11-lncRNA expression could be potential prognostic biomarkers in head and neck squamous cell carcinoma

The aim of our study is to construct the competing endogenous RNA (ceRNA) network of head and neck squamous cell carcinoma (HNSCC) and identify key long noncoding RNAs (lncRNAs) to predict prognosis. The genes whose expression were differentially in HNSCC and normal tissues were explored by the Cancer Genome Atlas database. The ceRNA network was constructed by the Cytoscape software. The lncRNAs which could estimate the overall survival were explored from Cox proportional hazards regression. There are 1997, 589, and 82 mRNAs, lncRNAs, and miRNAs whose expression were statistically significant different, respectively. Then, the network between miRNA and mRNA or miRNA and lncRNA was constructed by miRcode, miRDB, TargetScan, and miRanda. Five mRNAs, 10 lncRNAs, and 3 miRNAs were associated with overall survival. Then, 11-lncRNAs were found to be prognostic factors. Therefore, our research analyzed the potential signature of novel 11-lncRNA as candidate prognostic biomarker from the ceRNA network for patients with HNSCC.     

The AP2M1 gene expression is a promising biomarker for predicting survival of patients with hepatocellular carcinoma

There is a growing need for the discovery of new prognostic factors for cases where the scoring and staging system of hepatocellular carcinoma (HCC) does not result in a clear definition. We analyzed whether AP-2 complex subunit mu (AP2M1) expression could be a new prognostic marker for HCC based on the roles of AP2M1 in influencing hepatocyte growth factor (HGF) promoter regulation and hepatitis C virus (HCV) assembly. Patient data were extracted from cohorts of the Gene Expression Omnibus (GSE10186), International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). Differential expression value between matched cancer and normal liver was identified using ICGC cohort. Subsequently, we compared AP2M1 expression as a prognostic gene with other well-known prognostic genes for HCC, using the time-dependent area under the curve (AUC) of the Uno's C-index, the AUC value of the receiver operating characteristics at 5 years, Kaplan-Meier survival curve, and multivariate analysis. Particularly, TCGA and GSE10186 patients were divided into subgroups based on alcohol intake, hepatitis B, and C viral infections, and analyzed in the same methods. The AP2M1 expression values in patients with cancer were much higher than matched normal liver. The AP2M1 level showed excellent prognosis predictions in comparison with existing markers in the three independent cohorts (n = 647). In particular, it was more predictive of prognosis than other markers in alcohol intake and HCV infections. In conclusion, we were confident that AP2M1 provides sufficient value as a new prognostic marker for HCC especially patients with HCV infection and/or alcohol intake.     

Pinin acts as a poor prognostic indicator for renal cell carcinoma by reducing apoptosis and promoting cell migration and invasion

Pinin (PNN) was originally characterized as a desmosome-associated molecule. Its function and the mechanism of its regulation in renal cell carcinoma (RCC) are still undefined. Data on PNN expression, clinicopathological features, and prognosis of patients with RCC were obtained from The Cancer Genome Atlas (TCGA) database. Immunohistochemistry revealed high PNN expression in tumour cells. PNN expression showed negative correlation with survival in patients with RCC, acting as an independent prognostic factor in RCC. PNN up-regulation might be attributed to epigenetic alterations in RCC. Immunofluorescence revealed PNN expression mainly in the nucleus of RCC cells. The transfection of siRNA targeting the PNN gene resulted in enhanced apoptosis, which was detected by flow cytometry, and reduced cell migration and invasion, which were assessed using wound healing and transwell migration assay. Gene set enrichment analysis revealed associations between PNN expression and several signalling pathways involved in cancer progression, as a potential mechanism underlying the carcinogenicity of PNN. The analyses of the Tumor Immune Estimation Resource platform showed significant positive associations between high PNN expression and tumour immune infiltrating cells. PNN may function as an oncogenic factor by reducing apoptosis and promoting cell migration and invasion in RCC.     

Identification of two molecular subtypes of dysregulated immune lncRNAs in ovarian cancer

Long non-coding RNA (lncRNA) has increasingly been identified as a key regulator in pathologies such as cancer. Multiple platforms were used for comprehensive analysis of ovarian cancer to identify molecular subgroups. However, lncRNA and its role in mapping the ovarian cancer subpopulation are still largely unknown. RNA-sequencing and clinical characteristics of ovarian cancer were acquired from The Cancer Genome Atlas database (TCGA). A total of 52 lncRNAs were identified as aberrant immune lncRNAs specific to ovarian cancer. We redefined two different molecular subtypes, C1(188) and C2(184 samples), in “iClusterPlus” R package, among which C2 grouped ovarian cancer samples have higher survival probability and longer median survival time (P <0.05) with activated IFN-gamma response, Wound Healing and Cytotoxic lymphocytes signal; 456 differentially expressed genes were acquired in C1 and C2 subtypes using limma (3.40.6) package, among which 419 were up-regulated and 37 were down-regulated, in TCGA dataset. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis revealed that these genes were actively involved in ECM-receptor interaction, PI3K-Akt signaling pathway interaction KEGG pathway. Compared with the existing immune subtype, the Cluster2 sample showed a substantial increase in the proportion of the existing C2 immune subtype, accounting for 81.37%, which was associated with good prognosis. Our C1 subtype contains only 56.49% of the existing immune C1 and C4, which also explains the poor prognosis of C1. Furthermore, 52 immune-related lncRNAs were used to divide the TCGA-endometrial cancer and cervical cancer samples into two categories, and C2 had a good prognosis. The differentially expressed genes were highly correlated with immune-cell-related pathways. Based on lncRNA, two molecular subtypes of ovarian cancer were identified and had significant prognostic differences and immunological characteristics.     

Identification of five long noncoding RNAs signature and risk score for prognosis of bladder urothelial carcinoma

Nowadays, an increasing number of studies illustrated that bladder urothelial cancer (BLCA) may act as the most common subtype of urological malignancies with a high rate of recurrence and metastasis. In this study, we attempted to establish a prognostic model and identify the possible pathway crosstalk. Long noncoding RNAs (lncRNAs) and mRNA expression and corresponding clinical information of patients with BLCA were downloaded from The Cancer Genome Atlas (TCGA). The differentially expressed genes analysis, univariate Cox analysis, the least absolute shrinkage, and selection operator Cox (LASSO Cox) regression model were then applied to identify five crucial lncRNAs (AC092725.1, AC104071.1, AL023584.1, AL132642.1, and AL137804.1). The multivariate cox analysis was utilized to calculate the regression coefficients (β_i). The risk-score model was subsequently constructed as follows: (0.13541AC092725.1) + (0.20968AC104071.1) + (0.1525AL023584.1) − (0.14768AL132642.1) + (0.14387AL137804.1). Nomogram and assessment of overall survival (OS) prediction were verificated by the receiver operating characteristic curve in the testing group. As to 3-, 5-year OS prediction, the area under curve (AUC) for the nomogram of training data set was 0.83 and 0.86. Besides, the AUC (0.883 and 0.879) presented excellent predictive power in the testing group. In addition, the calibration plots validated the predictive performance of the nomogram. Weighted correlation network analysis (WGCNA) coupled with functional enrichment analysis contributed to explore the potential pathways, including PI3K-Akt, HIF-1, and Jak-STAT signaling pathways. Construction of the risk-score model and data analysis were both derived from multiple packages on the basis of the R platform chiefly.     

Genome-wide identification of a competing endogenous RNA network in cholangiocarcinoma

Cholangiocarcinoma (CCA) is the second widespread liver tumor with relatively poor survival. Increasing evidence in recent studies showed long noncoding RNAs (lncRNAs) exert a crucial impact on the development and progression of CCA based on the mechanism of competing endogenous RNAs (ceRNAs). However, functional roles and regulatory mechanisms of lncRNA-regulated ceRNA in CCA, are only partially understood. The expression profile of messenger RNAs (mRNAs), lncRNAs, and microRNAs (miRNAs) downloaded from The Cancer Genome Atlas were comprehensively investigated. Differential expression of these three types of RNA between CCA and corresponding precancerous tissues were screened out for further analysis. On the basis of interactive information generated from miRDB, miRTarBase, TargetScan, and miRcode public databases, we then constructed an mRNA-miRNA-lncRNA regulatory network. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology analyses were conducted to identify the biological function of the ceRNA network involved in CCA. As a result, 2883 mRNAs, 136 miRNAs, and 993 lncRNAs were screened out as differentially expressed RNAs in CCA. In addition, a ceRNA network in CCA was constructed, composing of 50 up and 27 downregulated lncRNAs, 14 up and 7 downregulated miRNAs, 29 up and 25 downregulated mRNAs. Finally, gene set enrichment and pathway analysis indicated our CCA-specific ceRNA network was related with cancer-related pathway and molecular function. In conclusion, our research identified a novel lncRNA-related ceRNA network in CCA, which might act as a potential therapeutic target for patients with CCA.     

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.     

Different expression profile of mRNA and long noncoding RNA in autoimmune thyroid diseases patients

Increasing evidence has found that long noncoding RNAs (lncRNAs) and message RNAs (mRNAs) play an important role in the progress of autoimmune thyroid diseases (AITD). So, in this study, the different expressed of lncRNA and mRNA was screened by microarray analysis and quantitative real-time polymerase chain reaction (PCR). To further investigate the relationship among the differentially expressed genes, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and Gene ontology (GO) were used for biofunctions and signaling pathways analysis, respectively. Finally, the interaction relationship between lncRNA and mRNAs was analysis with lncRNA-mRNA coexpression network. The result found that the abnormal expression of lncRNAs and mRNAs were 1615 and 1913, respectively. The altered genes included CD40LG, IFNG, CTLA4, FAS, STAT1, STAT3, and STAT4. These were enriched in presentation and antigen processing, Th1 and Th2 cell differentiation, natural killer cell-mediated cytotoxicity, B cell receptor signaling pathway, Th17 cell differentiation, and cell adhesion molecules (CAMs), all of which had been suggested to be associated with immunopathogenic mechanisms and AITD-induced pathophysiologic changes. A coexpression network profile was contained with 126 network nodes and 477 connections which were based on seven mRNAs and 119 interacted lncRNAs. The outcomes of differentially expressed lncRNAs and their coreralated mRNAs in our study revealed that lncRNAs involved in immunopathogenic mechanisms may play a crital role in the pathogenesis of AITD.