Identification of a five-gene signature with prognostic value in colorectal cancer

Colorectal cancer (CRC) ranks as one of the most commonly diagnosed malignancies worldwide. Although mortality rates have been decreasing, the prognosis of CRC patients is still highly dependent on the individual. Therefore, identifying and understanding novel biomarkers for CRC prognosis remains crucial. The gene expression profiles of five-gene expression omnibus (GEO) data sets of CRC were first downloaded. A total of 352 consistent differentially expressed genes (DEGs) were identified for CRC and paired with normal tissues. Functional analysis including gene ontology and Kyoto encyclopedia of genes and genomes pathway enrichment revealed that these DEGs were related to metabolic pathways, tight junctions, and the cell cycle. Ten hub DEGs were identified based on the search tool for the retrieval of interacting genes database and protein–protein interaction networks. By using univariate Cox proportional hazard regression analysis, we found 11 survival-related genes among these DEGs. We finally established a five-gene signature (kinesin family member 15, N-acetyltransferase 2, glutathione peroxidase 3, secretogranin II, and chloride channel accessory 1) with prognostic value in CRC by step multivariate Cox regression analysis. Based on this risk scoring system, patients in the high-risk group had significantly poorer survival results compared with those in the low-risk group (log-rank test, p < 0.0001). Finally, we validated our gene signature scoring system in two independent GEO cohorts (GSE17536 and GSE33113). We found all five of the signature genes to be DEGs in The Cancer Genome Atlas database. In conclusion, our findings suggest that our five DEG-based signature can provide a novel biomarker with useful applications in CRC prognosis.     

Development of a novel gene signature in patients without Helicobacter pylori infection gastric cancer

Gastric cancer (GC) is one of the most fatal common cancers in worldwide. Helicobacter pylori (H. pylori) infection is closely related to the development of GC, although the mechanism is still unclear. In our study, we aim to develop a robust messenger RNA (mRNA) signature associated with H. pylori (-) GC that can sensitively and efficiently predict the prognostic. The RNA-seq expression profile and corresponding clinical data of 598 gastric cancer samples and 63 normal samples obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus database. Using gene set enrichment analysis H. pylori (+) GC and H. pylori (-) GC patients and normal samples to select certain genes for further analysis. Using univariate and multivariate Cox regression model to establish a gene signature for predicting the overall survival (OS). Finally, we identified G2/M related seven-mRNA signature (TGFB1, EGF, MKI67, ILF3, INCENP, TNPO2, and CHAF1A) closely related to the prognosis of patients with H. pylori (-) GC. The seven-mRNA signature was identified to act as an independent prognostic biomarker by stratified analysis and multivariate Cox regression analysis. It was also validated on two test groups from TCGA and GSE15460 and shown that patients with high-risk scores based on the expression of the seven mRNAs had significantly shorter survival times compared to patients with low-risk scores (P < .0001). In this study, we developed a seven-mRNA signature related to G2/M checkpoint from H. pylori (-) GCs that as an independent biomarker potentially with a good performance in predicting OS and might be valuable for the clinical management for patients with GC.     

Five genes as a novel signature for predicting the prognosis of patients with laryngeal cancer

In this study, we purpose to investigate a novel five-gene signature for predicting the prognosis of patients with laryngeal cancer. The laryngeal cancer datasets were obtained from The Cancer Genome Atlas (TCGA). Both univariate and multivariate Cox regression analysis was applied to screening for prognostic differential expressed genes (DEGs), and a novel gene signature was obtained. The performance of this Cox regression model was tested by receiver operating characteristic (ROC) curves and area under the curve (AUC). Further survival analysis for each of the five genes was carried out through the Kaplan-Meier curve and Log-rank test. Totally, 622 DEGs were screened from the TCGA datasets in this study. We construct a five-gene signature through Cox survival analysis. Patients were divided into low- and high-risk groups depending on the median risk score, and a significant difference of the 5-year overall survival was found between these two groups (P < .05). ROC curves verified that this five-gene signature had good performance to predict the prognosis of laryngeal cancer (AUC = 0.862, P < .05). In conclusion, the five-gene signature consist of EMP1, HOXB9, DPY19L2P1, MMP1, and KLHDC7B might be applied as an independent prognosis predictor of laryngeal cancer.     

Identification of a novel glycolysis-related gene signature that can predict the survival of patients with lung adenocarcinoma

Lung cancer is one of the most malignant cancers worldwide, and lung adenocarcinoma (LUAD) is the most common histologic subtype. Thousands of biomarkers related to the survival and prognosis of patients with this cancer type have been investigated through database mining; however, the prediction effect of a single gene biomarker is not satisfactorily specific or sensitive. Thus, the present study aimed to develop a novel gene signature of prognostic values for patients with LUAD. Using a data-mining method, we performed expression profiling of 1145 mRNAs in large cohorts with LUAD (n = 511) from The Cancer Genome Atlas database. Using the Gene Set Enrichment Analysis, we selected 198 genes related to GLYCOLYSIS, which is the most important enrichment gene set. Moreover, these genes were identified using Cox proportional regression modeling. We established a risk score staging system to predict the outcome of patients with LUAD and subsequently identified four genes (AGRN, AKR1A1, DDIT4, and HMMR) that were closely related to the prognosis of patients with LUAD. The identified genes allowed us to classify patients into the high-risk group (with poor outcome) and low-risk group (with better outcome). Compared with other clinical factors, the risk score has a better performance in predicting the outcome of patients with LUAD, particularly in the early stage of LUAD. In conclusion, we developed a four-gene signature related to glycolysis by utilizing the Cox regression model and a risk staging model for LUAD, which might prove valuable for the clinical management of patients with LUAD.     

Identification of DNA methylation signature to predict prognosis in gastric adenocarcinoma

Gastric adenocarcinoma is an important death-related cancer. To find factors related to survival and prognosis, and thus improve recovery prospects, a powerful signature is needed. DNA methylation plays an important role in gastric adenocarcinoma processes and development, and here we report on the search for a significant DNA methylation gene to aid with the earlier diagnosis of gastric adenocarcinoma patients. A Cox proportional risk regression analysis and random survival forest algorithm were used to analyze gastric adenocarcinoma patients’ DNA methylation data from The Cancer Genome Atlas, a public database. DNA methylation gene signature consisting of five genes (SERPINA3, AP000357.4, GZMA, AC004702.2, and GREB1L) were selected. As the most accurate predictor, the area under the curve in the training and test group were 0.72 and 0.61, respectively. The signature was able to sort patients into high- and low-risk groups with meaningful overall survival rates (median: 18.36 vs 72.23 months, log-rank test, P < 0.001) in the training group, which predictive ability was validated in a test data set (median: 25.56 vs 58.80 months, log-rank test, P < 0.016). A multivariate Cox regression analysis showed the significant DNA methylation was an independent prediction prognostic factor for gastric adenocarcinoma patients. Functional analysis suggests that these signature genes may be related to pathways and biological processes associated with tumorigenesis. The significant DNA methylation gene could be a novel prediction and prognostic biomarker that both aids in the treatment and predicts the overall survival likelihoods of gastric adenocarcinoma patients.     

Identification of a five-gene signature of the RGS gene family with prognostic value in ovarian cancer

The RGS (regulator of G protein signaling) gene family, which includes negative regulators of G protein-coupled receptors, comprises important drug targets for malignant tumors. It is thus of great significance to explore the value of RGS family genes for diagnostic and prognostic prediction in ovarian cancer. The RNA-seq, immunophenotype, and stem cell index data of pan-cancer, The Cancer Genome Atlas (TCGA) data, and GTEx data of ovarian cancer were downloaded from the UCSC Xena database. In the pan-cancer database, the expression level of RGS1, RGS18, RGS19, and RGS13 was positively correlated with stromal and immune cell scores. Cancer patients with high RGS18 expression were more sensitive to cyclophosphamide and nelarabine, whereas those with high RGS19 expression were more sensitive to cladribine and nelarabine. The relationship between RGS family gene expression and overall survival (OS) and progression-free survival (PFS) of ovarian cancer patients was analyzed using the KM-plotter database, RGS17, RGS16, RGS1, and RGS8 could be used as diagnostic biomarkers of the immune subtype of ovarian cancer, and RGS10 and RGS16 could be used as biomarkers to predict the clinical stage of this disease. Further, Lasso cox analysis identified a five-gene risk score (RGS11, RGS10, RGS13, RGS4, and RGS3). Multivariate COX analysis showed that the risk score was an independent prognostic factor for patients with ovarian cancer. Immunohistochemistry and the HPA protein database confirmed that the five-gene signature is overexpressed in ovarian cancer. GSEA showed that it is mainly involved in the ECM-receptor interaction, TGF-beta signaling pathway, Wnt signaling pathway, and chemokine signaling pathway, which promote the occurrence and development of ovarian cancer. The prediction model of ovarian cancer constructed using RGS family genes is of great significance for clinical decision making and the personalized treatment of patients with ovarian cancer.     

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.     

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.     

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 validation of autophagy-related prognostic signature for head and neck squamous cell carcinoma

BackgroundMany studies have demonstrated that autophagy plays a significant role in regulating tumor growth and progression. However, the effect of autophagy-related genes (ARGs) on the prognosis have rarely been analyzed in head and neck squamous cell carcinoma (HNSCC).MethodsWe obtained differentially expressed ARGs from HNSCC mRNA data in The Cancer Genome Atlas (TCGA) database. And then we performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses to explore the autophagy-related biological functions. The overall survival (OS)-related and disease specific survival (DSS)-related ARGs were identified by univariate Cox regression analyses. With these genes, we established OS-related and DSS-related risk signature by LASSO regression method, respectively. We validated the reliability of the risk signature with receiver operating characteristic (ROC) analysis, Kaplan-Meier survival curves, clinical correlation analysis, and nomogram. Then we analyzed relationships between risk signature and immune cell infiltration.ResultsWe established the prognostic signatures based on 14 ARGs for OS and 12 ARGs for DSS. The ROC curves, survival analysis, and nomogram validated the predictive accuracy of the models. Clinic correlation analysis showed that the risk group was closely related to Stage, pathological T stage, pathological N stage and human papilloma virus (HPV) subtype. Cox regression demonstrated that the risk score was an independent predictor for the prognosis of HNSCC patients. Furthermore, patients in low-risk score group exhibited higher immunescore and distinct immune cell infiltration than high-risk score group. And we further analysis revealed that the copy number alterations (CNAs) of ARGs-based signature affected the abundance of tumor-infiltrating immune cells.ConclusionIn this study, we identified novel autophagy-related signature for the prediction of OS and DSS in patients with HNSCC. Meanwhile, our study provides a novel sight to understand the role of autophagy and elucidate the important role of autophagy in tumor immune microenvironment (TIME) of HNSCC.     

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.     

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.     

Screening of significant biomarkers related with prognosis of liver cancer by lncRNA-associated ceRNAs analysis

This study aimed to identify significant biomarkers related to the prognosis of liver cancer using long noncoding RNA (lncRNA)-associated competing endogenous RNAs (ceRNAs) analysis. Differentially expressed mRNA and lncRNAs between liver cancer and paracancerous tissues were screened, and the functions of these mRNAs were predicted by gene ontology and pathway enrichment analyses. A ceRNA network consisting of differentially expressed mRNAs and lncRNAs was constructed. LncRNA FENDRR and lncRNA HAND2-AS1 were hub nodes in the ceRNA network. A risk score assessment model consisting of eight genes (PDE2A, ESR1, FBLN5, ALDH8A1, AKR1D1, EHHADH, ADRA1A, and GNE) associated with prognosis were developed. Multivariate Cox regression suggested that both pathologic_T and risk group could be regarded as independent prognostic factors. Furthermore, a nomogram model consisting of pathologic_T and risk group showed a good prediction ability for predicting the survival rate of liver cancer patients. The nomogram model consisting of pathologic_T and a risk score assessment model could be regarded as an independent factor for predicting prognosis of liver cancer.     

Novel gene signatures for prognosis prediction in ovarian cancer

Ovarian cancer (OV) is one of the leading causes of cancer deaths in women worldwide. Late diagnosis and heterogeneous treatment result to poor survival outcomes for patients with OV. Therefore, we aimed to develop novel biomarkers for prognosis prediction from the potential molecular mechanism of tumorigenesis. Eight eligible data sets related to OV in GEO database were integrated to identify differential expression genes (DEGs) between tumour tissues and normal. Enrichment analyses discovered DEGs were most significantly enriched in G2/M checkpoint signalling pathway. Subsequently, we constructed a multi‐gene signature based on the LASSO Cox regression model in the TCGA database and time‐dependent ROC curves showed good predictive accuracy for 1‐, 3‐ and 5‐year overall survival. Utility in various types of OV was validated through subgroup survival analysis. Risk scores formulated by the multi‐gene signature stratified patients into high‐risk and low‐risk, and the former inclined worse overall survival than the latter. By incorporating this signature with age and pathological tumour stage, a visual predictive nomogram was established, which was useful for clinicians to predict survival outcome of patients. Furthermore, SNRPD1 and EFNA5 were selected from the multi‐gene signature as simplified prognostic indicators. Higher EFNA5 expression or lower SNRPD1 indicated poorer outcome. The correlation between signature gene expression and clinical characteristics was observed through WGCNA. Drug‐gene interaction was used to identify 16 potentially targeted drugs for OV treatment. In conclusion, we established novel gene signatures as independent prognostic factors to stratify the risk of OV patients and facilitate the implementation of personalized therapies.     

Identification of a 4‐mRNA metastasis‐related prognostic signature for patients with breast cancer

Metastasis‐related mRNAs have showed great promise as prognostic biomarkers in various types of cancers. Therefore, we attempted to develop a metastasis‐associated gene signature to enhance prognostic prediction of breast cancer (BC) based on gene expression profiling. We firstly screened and identified 56 differentially expressed mRNAs by analysing BC tumour tissues with and without metastasis in the discovery cohort (GSE102484, n = 683). We then found 26 of these differentially expressed genes were associated with metastasis‐free survival (MFS) in the training set (GSE20685, n = 319). A metastasis‐associated gene signature built using a LASSO Cox regression model, which consisted of four mRNAs, can classify patients into high‐ and low‐risk groups in the training cohort. Patients with high‐risk scores in the training cohort had shorter MFS (hazard ratio [HR] 3.89, 95% CI 2.53‐5.98; P < 0.001), disease‐free survival (DFS) (HR 4.69, 2.93‐7.50; P < 0.001) and overall survival (HR 4.06, 2.56‐6.45; P < 0.001) than patients with low‐risk scores. The prognostic accuracy of mRNAs signature was validated in the two independent validation cohorts (GSE21653, n = 248; GSE31448, n = 246). We then developed a nomogram based on the mRNAs signature and clinical‐related risk factors (T stage and N stage) that predicted an individual's risk of disease, which can be assessed by calibration curves. Our study demonstrated that this 4‐mRNA signature might be a reliable and useful prognostic tool for DFS evaluation and will facilitate tailored therapy for BC patients at different risk of disease.     

DNA methylation profiling to predict recurrence risk in stage Ι lung adenocarcinoma: Development and validation of a nomogram to clinical management

Increasing evidence suggested DNA methylation may serve as potential prognostic biomarkers; however, few related DNA methylation signatures have been established for prediction of lung cancer prognosis. We aimed at developing DNA methylation signature to improve prognosis prediction of stage I lung adenocarcinoma (LUAD). A total of 268 stage I LUAD patients from the Cancer Genome Atlas (TCGA) database were included. These patients were separated into training and internal validation datasets. GSE39279 was used as an external validation set. A 13‐DNA methylation signature was identified to be crucially relevant to the relapse‐free survival (RFS) of patients with stage I LUAD by the univariate Cox proportional hazard analysis and the least absolute shrinkage and selection operator (LASSO) Cox regression analysis and multivariate Cox proportional hazard analysis in the training dataset. The Kaplan‐Meier analysis indicated that the 13‐DNA methylation signature could significantly distinguish the high‐ and low‐risk patients in entire TCGA dataset, internal validation and external validation datasets. The receiver operating characteristic (ROC) analysis further verified that the 13‐DNA methylation signature had a better value to predict the RFS of stage I LUAD patients in internal validation, external validation and entire TCGA datasets. In addition, a nomogram combining methylomic risk scores with other clinicopathological factors was performed and the result suggested the good predictive value of the nomogram. In conclusion, we successfully built a DNA methylation‐associated nomogram, enabling prediction of the RFS of patients with stage I LUAD.     

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.