The discovery of a novel eight-mRNA-lncRNA signature predicting survival of hepatocellular carcinoma patients

Increasing evidence indicates that the expressions of messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs) undergo a frequent and aberrant change in carcinogenesis and cancer development. But some research was carried out on mRNA-lncRNA signatures for prediction of hepatocellular carcinoma (HCC) prognosis. We aimed to establish an mRNA-lncRNA signature to improve the ability to predict HCC patients’ survival. The subjects from the cancer genome atlas (TCGA) data set were randomly divided into two parts: training data set (n = 246) and testing data set (n = 124). Using computational methods, we selected eight gene signatures (five mRNAs and three lncRNAs) to generate the risk score model, which were significantly correlated with overall survival of patients with HCC in both training and testing data set. The signature had the ability to classify the patients in training data set into a high-risk group and low-risk group with significantly different overall survival (hazard ratio = 4.157, 95% confidence interval = 2.648-6.526, P < 0.001). The prognostic value was further validated in testing data set and the entire data set. Further analysis revealed that this signature was independent of tumor stage. In addition, Gene Set Enrichment Analysis suggested that high risk score group was associated with cell proliferation and division related pathways. Finally, we developed a well-performed nomogram integrating the prognostic signature and other clinical information to predict 3- and 5-year overall survival. In conclusion, the prognostic mRNAs and lncRNAs identified in our study indicate their potential role in HCC biogenesis. The risk score model based on the mRNA-lncRNA may be an efficient classification tool to evaluate the prognosis of patients’ with HCC.     

A novel prognostic signature of immune‐related genes for patients with colorectal cancer

Colorectal cancer (CRC) is one of the most commonly diagnosed cancers with an estimated 1.8 million new cases worldwide and associated with high mortality rates of 881 000 CRC‐related deaths in 2018. Screening programs and new therapies have only marginally improved the survival of CRC patients. Immune‐related genes (IRGs) have attracted attention in recent years as therapeutic targets. The aim of this study was to identify an immune‐related prognostic signature for CRC. To this end, we combined gene expression and clinical data from the CRC data sets of The Cancer Genome Atlas (TCGA) into an integrated immune landscape profile. We identified a total of 476 IRGs that were differentially expressed in CRC vs normal tissues, of which 18 were survival related according to univariate Cox analysis. Stepwise multivariate Cox proportional hazards analysis established an immune‐related prognostic signature consisting of SLC10A2, FGF2, CCL28, NDRG1, ESM1, UCN, UTS2 and TRDC. The predictive ability of this signature for 3‐ and 5‐year overall survival was determined using receiver operating characteristics (ROC), and the respective areas under the curve (AUC) were 79.2% and 76.6%. The signature showed moderate predictive accuracy in the validation and GSE38832 data sets as well. Furthermore, the 8‐IRG signature correlated significantly with tumour stage, invasion, lymph node metastasis and distant metastasis by univariate Cox analysis, and was established an independent prognostic factor by multivariate Cox regression analysis for CRC. Gene set enrichment analysis (GSEA) revealed a relationship between the IRG prognostic signature and various biological pathways. Focal adhesions and ECM‐receptor interactions were positively correlated with the risk scores, while cytosolic DNA sensing and metabolism‐related pathways were negatively correlated. Finally, the bioinformatics results were validated by real‐time RT?qPCR. In conclusion, we identified and validated a novel, immune‐related prognostic signature for patients with CRC, and this signature reflects the dysregulated tumour immune microenvironment and has a potential for better CRC patient management.     

An EMT‐related gene signature for the prognosis of human bladder cancer

The transition from non–muscle‐invasive bladder cancer (NMIBC) to muscle‐invasive bladder cancer (MIBC) is detrimental to bladder cancer (BLCA) patients. Here, we aimed to study the underlying mechanism of the subtype transition. Gene set variation analysis (GSVA) revealed the epithelial‐mesenchymal transition (EMT) signalling pathway with the most positive correlation in this transition. Then, we built a LASSO Cox regression model of an EMT‐related gene signature in BLCA. The patients with high risk scores had significantly worse overall survival (OS) and disease‐free survival (DFS) than those with low risk scores. The EMT‐related gene signature also performed favourably in the accuracy of prognosis and in the subtype survival analysis. Univariate and multivariate Cox regression analyses demonstrated that the EMT‐related gene signature, pathological N stage and age were independent prognostic factors for predicting survival in BLCA patients. Furthermore, the predictive nomogram model was able to effectively predict the outcome of BLCA patients by appropriately stratifying the risk score. In conclusion, we developed a novel EMT‐related gene signature that has tumour‐promoting effects, acts as a negative independent prognostic factor and might facilitate personalized counselling and treatment in BLCA.     

An integrated mRNA-lncRNA signature for relapse prediction in laryngeal cancer

Patients with laryngeal cancer with early relapse usually have a poor prognosis. In this study, we aimed to identify a multi-gene signature to improve the relapse prediction in laryngeal cancer. One microarray data set GSE27020 (training set, N = 109) and one RNA-sequencing data set (validation set, N = 85) were included into the analysis. In the training set, the microarray expression profile was re-annotated into an mRNA-long noncoding RNA (lncRNA) biphasic profile. Then, LASSO Cox regression model identified nine relapse-related RNA (eight mRNA and one lncRNA), and a risk score was calculated for each sample according to the model coefficients. Patients with high-risk showed poorer relapse-free survival than patients with low risk (hazard ratios (HR): 6.189, 95% confidence interval (CI): 3.075-12.460, P < 0.0001). The risk score demonstrated good accuracy in predicting the relapse (area under time-dependent receiver-operating characteristic (AUC): 0.859 at 1 year, 0.822 at 3 years, and 0.815 at 5 years). The results were validated in the validation set (HR: 3.762, 95% CI: 1.594-8.877, P = 0.011; AUC: 0.770 at 1 year, 0.769 at 3 years, and 0.728 at 5 years). The multivariate analysis reached consistent results after adjustment by multiple confounders. When compared with a 27-gene signature, a 2-lncRNA signature, and Tumor-Node-Metastasis stage, the risk score also showed better performance (P < 0.05). In conclusion, we successfully developed a robust mRNA-lncRNA signature that can accurately predict the relapse in laryngeal cancer.     

Development and validation of a novel epigenetic signature for predicting prognosis in colon cancer

Epigenetic factors play a critical role in carcinogenesis by imparting a distinct feature to the chromatin architecture. The present study aimed to develop a novel epigenetic signature for evaluating the relapse-free survival of colon cancer patients. Public microarray datasets were acquired from the Gene Expression Omnibus databases: GSE39582, GSE17538, GSE33113, and GSE37892 set. Patients from GSE39582 set were randomized 1:1 into training and internal validation series. Patients were divided into high-risk and low-risk groups in training series based on a set of 11 epigenetic factors (p < .001). The good reproducibility for the prognostic value of the epigenetic signature was confirmed in the internal validation series (p < .001), external validation series (a combination of GSE17538 set, GSE33113 set, and GSE37892 set; p = .018), and entire series (p < .001). Furthermore, a nomogram, which integrated the epigenetic signature, pathological stage, and postoperative chemotherapy, was developed based on the GSE39582 set. The time-dependent receiver operating characteristic curve at 1 year demonstrated that the comprehensive signature presented superior prognostic value than the pathological stage. In conclusion, an epigenetic signature, which could be utilized to divide colon cancer patients into two groups with significantly different risk of relapse, was established. This biomarker would aid in identifying patients who require an intensive follow-up and aggressive therapeutic intervention.     

A prognostic eight-gene expression signature for patients with breast cancer receiving adjuvant chemotherapy

Breast cancer is a popularly diagnosed malignant tumor. Genomic profiling studies suggest that breast cancer is a disease with heterogeneity. Chemotherapy is one of the chief means to treat breast cancer, while its responses and clinical outcomes vary largely due to the conventional clinicopathological factors and inherent chemosensitivity of breast cancer. Using the least absolute shrinkage and selection operator (LASSO) Cox regression model, our study established a multi-mRNA-based signature model and constructed a relative nomogram in predicting distant-recurrence-free survival for patients receiving surgery and following chemotherapy. We constructed a signature of eight mRNAs (IPCEF1, SYNDIG1, TIGIT, SPESP1, C2CD4A, CLCA2, RLN2, and CCL19) with the LASSO model, which was employed to separate subjects into groups with high- and low-risk scores. Obvious differences of distant-recurrence-free survival were found between these two groups. This eight-mRNA-based signature was independently associated with the prognosis and had better prognostic value than classical clinicopathologic factors according to multivariate Cox regression results. Receiver operating characteristic results demonstrated excellent performance in diagnosing 3-year distant-recurrence by the eight-mRNA signature. A nomogram that combined both the eight-mRNA-based signature and clinicopathological risk factors was constructed. Comparing with an ideal model, the nomograms worked well both in the training and validation sets. Through the results that the eight-mRNA signature effectively classified patients into low- and high-risk of distant recurrence, we concluded that this eight-mRNA-based signature played a promising predictive role in prognosis and could be clinically applied in breast cancer patients receiving adjuvant chemotherapy.     

Identification of a novel cell cycle-related gene signature predicting survival in patients with gastric cancer

Gastric cancer (GC) is one of the most fatal cancers in the world. Thousands of biomarkers have been explored that might be related to survival and prognosis via database mining. However, the prediction effect of single gene biomarkers is not specific enough. Increasing evidence suggests that gene signatures are emerging as a possible better alternative. We aimed to develop a novel gene signature to improve the prognosis prediction of GC. Using the messenger RNA (mRNA)-mining approach, we performed mRNA expression profiling in a large GC cohort (n = 375) from The Cancer Genome Atlas (TCGA) database. Gene Set Enrichment Analysis (GSEA) was performed, and we recovered genes related to the G2/M checkpoint, which we identified with a Cox proportional regression model. We identified a set of five genes (MARCKS, CCNF, MAPK14, INCENP, and CHAF1A), which were significantly associated with overall survival (OS) in the test series. Based on this five-gene signature, the test series patients could be classified into high-risk or low-risk subgroups. Multivariate Cox regression analysis indicated that the prognostic power of this five-gene signature was independent of clinical features. In conclusion, we developed a five-gene signature related to the cell cycle that can predict survival for GC. Our findings provide novel insight that is useful for understanding cell cycle mechanisms and for identifying patients with GC with poor prognoses.     

Identification of a four-gene metabolic signature predicting overall survival for hepatocellular carcinoma

While hundreds of consistently altered metabolic genes had been identified in hepatocellular carcinoma (HCC), the prognostic role of them remains to be further elucidated. Messenger RNA expression profiles and clinicopathological data were downloaded from The Cancer Genome Atlas—Liver Hepatocellular Carcinoma and GSE14520 data set from the Gene Expression Omnibus database. Univariate Cox regression analysis and lasso Cox regression model established a novel four-gene metabolic signature (including acetyl-CoA acetyltransferase 1, glutamic-oxaloacetic transaminase 2, phosphatidylserine synthase 2, and uridine-cytidine kinase 2) for HCC prognosis prediction. Patients in the high-risk group shown significantly poorer survival than patients in the low-risk group. The signature was significantly correlated with other negative prognostic factors such as higher α-fetoprotein. The signature was found to be an independent prognostic factor for HCC survival. Nomogram including the signature shown some clinical net benefit for overall survival prediction. Furthermore, gene set enrichment analyses revealed several significantly enriched pathways, which might help explain the underlying mechanisms. Our study identified a novel robust four-gene metabolic signature for HCC prognosis prediction. The signature might reflect the dysregulated metabolic microenvironment and provided potential biomarkers for metabolic therapy and treatment response prediction in HCC.     

Development and verification of an immune-related gene pairs prognostic signature in ovarian cancer

Ovarian cancer (OV) is the most common gynaecological cancer worldwide. Immunotherapy has recently been proven to be an effective treatment strategy. The work here attempts to produce a prognostic immune-related gene pair (IRGP) signature to estimate OV patient survival. The Gene Expression Omnibus (GEO) and Cancer Genome Atlas (TCGA) databases provided the genetic expression profiles and clinical data of OV patients. Based on the InnateDB database and the least absolute shrinkage and selection operator (LASSO) regression model, we first identified a 17-IRGP signature associated with survival. The average area under the curve (AUC) values of the training, validation, and all TCGA sets were 0.869, 0.712, and 0.778, respectively. The 17-IRGP signature noticeably split patients into high- and low-risk groups with different prognostic outcomes. As suggested by a functional study, some biological pathways, including the Toll-like receptor and chemokine signalling pathways, were significantly negatively correlated with risk scores; however, pathways such as the p53 and apoptosis signalling pathways had a positive correlation. Moreover, tumour stage III, IV, grade G1/G2, and G3/G4 samples had significant differences in risk scores. In conclusion, an effective 17-IRGP signature was produced to predict prognostic outcomes in OV, providing new insights into immunological biomarkers.     

A novel multidimensional signature predicts prognosis in hepatocellular carcinoma patients

The abnormal expression of microRNAs (miRNAs) or protein-coding genes (PCGs) have been found to be associated with the prognosis of hepatocellular carcinoma (HCC) patients. Using bioinformatics analysis methods including Cox’s proportional hazards regression analysis, the random survival forest algorithm, Kaplan–Meier, and receiver operating characteristic (ROC) curve analysis, we mined the gene expression profiles of 469 HCC patients from The Cancer Genome Atlas (n = 379) and Gene Expression Omnibus (GSE14520; n = 90) public database. We selected a signature comprising one protein-coding gene (PCG; DNA polymerase μ) and three miRNAs (hsa-miR-149-5p, hsa-miR-424-5p, hsa-miR-579-5p) with highest accurate prediction (area under the ROC curve [AUC] = 0.72; n = 189) from the training data set. The signature stratified patients into high- and low-risk groups with significantly different survival (median 27.9 vs. 55.2 months, log-rank test, p < 0.001) in the training data set, and its risk stratification ability were validated in the test data set (median 47.4 vs. 84.4 months, log-rank test, p = 0.03) and an independent data set (median 31.0 vs. 46.0 months, log-rank test, p = 0.01). Multivariable Cox regression analysis showed that the signature was an independent prognostic factor. And the signature was proved to have a better survival prediction power than tumor–node–metastasis (TNM) stage (AUC _signature = 0.72/0.64/0.62 vs. AUC _TNM = 0.65/0.61/0.61; p < 0.05). Moreover, we validated the expression of these prognostic genes from the PCG-miRNA signature in Huh-7 cell by real-time polymerase chain reaction. In conclusion, we found a signature that can predict survival of HCC patients and serve as a prognostic marker for HCC.     

A novel microRNA-based signature predicts prognosis among nasopharyngeal cancer patients

Nasopharyngeal cancer is one of the most common malignant tumors in the head and neck. Identification of promising miRNA biomarkers might benefit a lot to the detection of nasopharyngeal carcinoma. miRNA expression profile and clinical information were obtained from two microarray profiling data sets from the Gene Expression Omnibus (GEO) database. miRNA signature model was constructed via univariate Cox survival analysis, multivariate Cox survival analysis, and least absolute shrinkage and selection operator Cox regression analysis. Kaplan–Meier curve, area under the curve (AUC), decision curve analysis, Box plot, and nomogram were used to evaluate the prognosis of the model to patients. 67 up-regulated and 93 down-regulated miRNAs were identified from GEO microarray data sets (P < 0.05). A three-miRNA signature (has-miR-142-3p, has-miR-29c, and has-miR-30e) was obviously associated with the overall survival of nasopharyngeal carcinoma patients (P  < 0.001). The AUCs for the signature were 0.74, 0.7 for the training set and external validation set. The AUC of disease free survival and distant metastasis-free survival were also high. The model has better clinical independence and has better clinical prediction effect when combined with clinical characteristics (P < 0.0001). Compared with the published models, our model had a higher AUC. Our results revealed that a three-miRNA signature was a potential novel prognostic biomarker for nasopharyngeal carcinoma.Impact statementNasopharyngeal cancer is one of the most common malignant tumors in the head and neck. Identification of promising miRNA biomarkers might benefit a lot to the detection of nasopharyngeal carcinoma. A three-miRNA signature (has-miR-142-3p, has-miR-29c, and has-miR-30e) was obviously associated with the overall survival of nasopharyngeal carcinoma patients. The model has better clinical independence and has better clinical prediction effect when combined with clinical characteristics. Our results revealed that a three-miRNA signature was a potential novel prognostic biomarker for nasopharyngeal 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.     

Development and validation of a novel 15‐CpG‐based signature for predicting prognosis in triple‐negative breast cancer

DNA methylation is an important biological regulatory mechanism that changes gene expression without altering the DNA sequence. Increasing studies have revealed that DNA methylation data play a vital role in the field of oncology. However, the methylation site signature in triple‐negative breast cancer (TNBC) remains unknown. In our research, we analysed 158 TNBC samples and 98 noncancerous samples from The Cancer Genome Atlas (TCGA) in three phases. In the discovery phase, 86 CpGs were identified by univariate Cox proportional hazards regression (CPHR) analyses to be significantly correlated with overall survival (P < 0.01). In the training phase, these candidate CpGs were further narrowed down to a 15‐CpG‐based signature by conducting least absolute shrinkage and selector operator (LASSO) Cox regression in the training set. In the validation phase, the 15‐CpG‐based signature was verified using two different internal sets and one external validation set. Furthermore, a nomogram comprising the CpG‐based signature and TNM stage was generated to predict the 1‐, 3‐ and 5‐year overall survival in the primary set, and it showed excellent performance in the three validation sets (concordance indexes: 0.924, 0.974 and 0.637). This study showed that our nomogram has a precise predictive effect on the prognosis of TNBC and can potentially be implemented for clinical treatment and diagnosis.     

Inference of patient‐specific subpathway activities reveals a functional signature associated with the prognosis of patients with breast cancer

Breast cancer is one of the most deadly forms of cancer in women worldwide. Better prediction of breast cancer prognosis is essential for more personalized treatment. In this study, we aimed to infer patient‐specific subpathway activities to reveal a functional signature associated with the prognosis of patients with breast cancer. We integrated pathway structure with gene expression data to construct patient‐specific subpathway activity profiles using a greedy search algorithm. A four‐subpathway prognostic signature was developed in the training set using a random forest supervised classification algorithm and a prognostic score model with the activity profiles. According to the signature, patients were classified into high‐risk and low‐risk groups with significantly different overall survival in the training set (median survival of 65 vs 106 months, P = 1.82e‐13) and test set (median survival of 75 vs 101 months, P = 4.17e‐5). Our signature was then applied to five independent breast cancer data sets and showed similar prognostic values, confirming the accuracy and robustness of the subpathway signature. Stratified analysis suggested that the four‐subpathway signature had prognostic value within subtypes of breast cancer. Our results suggest that the four‐subpathway signature may be a useful biomarker for breast cancer prognosis.     

Analysis of methylation-driven genes for predicting the prognosis of patients with head and neck squamous cell carcinoma

To help provide evidence for prognosis prediction and personalized targeted therapy for patients with head and neck squamous cell carcinoma (HNSCC), we investigated prognosis-specific methylation-driven genes in HNSCC. Survival time data, RNA sequencing data, and methylation data for HNSCC patients were downloaded from The Cancer Genome Atlas. The MethylMix R package based on the β mixture model was utilized to screen genes with different methylation statuses in tumor tissues and adjacent normal tissues, and a total of 182 HNSCC-related methylation-driven genes were then identified. A survival prediction scoring model based on multivariate Cox analysis was developed to screen the genes related to the prognosis of HNSCC, and a linear risk model of the methylation status of six genes (INA, LINC01354, TSPYL4, MAGEB2, EPHX3, and ZNF134) was constructed. The prognostic values of the six genes were further independently explored by survival analysis combined with methylation and gene expression analyses. The 5-year survival rate in the high-risk group of patients in the test set was 30.4% (95% CI: 22.7%-40.8%) and that in the low-risk group of patients was 65.5% (95% CI: 56.1%-76.5%). The area under the receiver operating characteristic curve for the model was 0.723, which further verified the specificity and sensitivity of the model. In addition, subsequent combined survival analysis revealed that all six genes could be used as independent prognostic markers and thus might be potential drug targets. The innovative method provides new insight into the molecular mechanism and prognosis of HNSCC.