Development and validation of a CTNNB1-associated metabolic prognostic model for hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a heterogeneous malignancy closely related to metabolic reprogramming. We investigated how CTNNB1 mutation regulates the HCC metabolic phenotype and thus affects the prognosis of HCC. We obtained the mRNA expression profiles and clinicopathological data from The Cancer Genome Atlas (TCGA), the International Cancer Genomics Consortium (ICGC) and the Gene Expression Omnibus database ( GSE14520 and GSE116174 ). We conducted gene set enrichment analysis on HCC patients with and without mutant CTNNB1 through TCGA dataset. The Kaplan-Meier analysis and univariate Cox regression analysis assisted in screening metabolic genes related to prognosis, and the prognosis model was constructed using the Lasso and multivariate Cox regression analysis. The prognostic model showed good prediction performance in both the training cohort (TCGA) and the validation cohorts (ICGC, GSE14520 , GSE116174 ), and the high-risk group presented obviously poorer overall survival compared with low-risk group. Cox regression analysis indicated that the risk score can be used as an independent predictor for the overall survival of HCC. The immune infiltration in different risk groups was also evaluated in this study to explore underlying mechanisms. This study is also the first to describe an metabolic prognostic model associated with CTNNB1 mutations and could be implemented for determining the prognoses of individual patients in clinical practice.     

Integrated analysis of gene expression and DNA methylation datasets identified key genes and a 6-gene prognostic signature for primary lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the main subtype of non-small cell lung cancer with a poor survival prognosis. In our study, gene expression, DNA methylation, and clinicopathological data of primary LUAD were utilized to identify potential prognostic markers for LUAD, which were recruited from The Cancer Genome Atlas (TCGA) database. Univariate regression analysis showed that there were 21 methylation-associated DEGs related to overall survival (OS), including 9 down- and 12 up-regulated genes. The 12 up-regulated genes with hypomethylation may be risky genes, whereas the other 9 down-regulated genes with hypermethylation might be protective genes. By using the Step-wise multivariate Cox analysis, a methylation-associated 6-gene (consisting of CCL20, F2, GNPNAT1, NT5E, B3GALT2, and VSIG2) prognostic signature was constructed and the risk score based on this gene signature classified patients into high- or low-risk groups. Patients of the high-risk group had shorter OS than those of the low-risk group in both the training and validation cohort. Multivariate Cox analysis and the stratified analysis revealed that the risk score was an independent prognostic factor for LUAD patients. The methylation-associated gene signature may serve as a prognostic factor for LUAD patients and the represent hypermethylated or hypomethylated genes might be potential targets for LUAD therapy.     

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.     

Prognostic and predictive role of a metabolic rate‐limiting enzyme signature in hepatocellular carcinoma

ObjectivesAbnormal expression of metabolic rate‐limiting enzymes drives the occurrence and progression of hepatocellular carcinoma (HCC). This study aimed to elucidate the comprehensive model of metabolic rate‐limiting enzymes associated with the prognosis of HCC.Materials and MethodsHCC animal model and TCGA project were used to screen out differentially expressed metabolic rate‐limiting enzyme. Cox regression, least absolute shrinkage and selection operation (LASSO) and experimentally verification were performed to identify metabolic rate‐limiting enzyme signature. The area under the receiver operating characteristic curve (AUC) and prognostic nomogram were used to assess the efficacy of the signature in the three HCC cohorts (TCGA training cohort, internal cohort and an independent validation cohort).ResultsA classifier based on three rate‐limiting enzymes (RRM1, UCK2 and G6PD) was conducted and serves as independent prognostic factor. This effect was further confirmed in an independent cohort, which indicated that the AUC at year 5 was 0.715 (95% CI: 0.653‐0.777) for clinical risk score, whereas it was significantly increased to 0.852 (95% CI: 0.798‐0.906) when combination of the clinical with signature risk score. Moreover, a comprehensive nomogram including the signature and clinicopathological aspects resulted in significantly predict the individual outcomes.ConclusionsOur results highlighted the prognostic value of rate‐limiting enzymes in HCC, which may be useful for accurate risk assessment in guiding clinical management and treatment decisions.     

Prognostic Value of E-cadherin-, CD44-, and MSH2-associated Nomograms in Patients With Stage II and III Colorectal Cancer

BACKGROUND: To evaluate the prognostic value of E-cadherin, CD44, and MSH2 expression for colorectal cancer (CRC) and construct nomograms that can predict prognosis. METHODS: We retrospectively analyzed the expression of E-cadherin, CD44, and MSH2 in 223 paraffin-embedded stage II and III CRC specimens using immunohistochemistry in the training cohort. Their prognostic values were assessed using Kaplan–Meier curves and univariate and multivariate COX regression models. Moreover, a number of risk factors were used to form nomograms to evaluate survival, and Harrell's concordance index (C-index) was used to evaluate the predictive accuracy. Further validation of the nomograms was performed in an independent cohort of 115 cases. RESULTS: Low E-cadherin expression and low CD44 expression were significantly associated with diminished overall survival (OS) and disease-free survival (DFS) in stage II and III CRC patients and patients with negative MSH2 expression had better clinical outcomes. Moreover, the multivariate COX analysis identified E-cadherin, CD44 and MSH2 expression as independent prognostic factors for DFS and OS. Using these three markers and three clinicopathological risk variables, two nomograms were constructed and externally validated for predicting OS and DFS (C-index: training cohort, 0.779 (95% CI 0.722–0.835) and 0.771 (0.720–0.822), respectively; validation cohort, 0.773 (0.709–0.837) and 0.670 (0.594–0.747), respectively). CONCLUSION: The expression levels of E-cadherin, CD44 and MSH2 were independent prognostic factors for stage II and III CRC patients. By incorporating clinicopathological features and these biomarkers, we have established two nomograms that could be used to make individualized predictions for OS and DFS.     

Identification of a three-miRNA signature as a novel potential prognostic biomarker in patients with clear cell renal cell carcinoma

Current studies suggest that some microRNAs (miRNAs) are associated with prognosis in clear cell renal cell carcinoma (ccRCC). In this paper, we aimed to identify a miRNAs signature to improve prognostic prediction for ccRCC patients. Using ccRCC RNA-Seq data of The Cancer Genome Atlas (TCGA) database, we identified 177 differentially expressed miRNAs between ccRCC and paracancerous tissue. Then all the ccRCC tumor samples were divided into training set and validation set randomly. Three-miRNA signature including miR130b, miR-18a, and miR-223 were constructed by the least absolute shrinkage and selection operator (LASSO) Cox regression model in training set. According to optimal cut-off value of three-miRNA signature risk score, all the patients could be classified into high-risk group and low-risk group significantly. Survival of patients was significantly different between two groups (hazard ratio, 5.58, 95% confidence interval, 3.17-9.80; P < 0.0001), and three-miRNA signature performed favorably prognostic and predictive accuracy. The results were further validated in the validation set and total set. Multivariate Cox regression analyses and subgroup analyses showed that three-miRNA signature was an independent prognostic factor. Two nomograms that integrated three-miRNA signature and three clinicopathological risk factors were constructed to predict overall survival and disease-free survival after surgery for ccRCC patients. Functional enrichment analysis showed the possible roles of three-miRNA signature in some cancer-associated biological processes and pathways. In conclusion, we developed a novel three-miRNA signature that performed reliable prognostic for patient survival with ccRCC, it might facilitate ccRCC patients counseling and individualize management.     

A 16-mRNA signature optimizes recurrence-free survival prediction of Stages II and III gastric cancer

High-throughput messenger RNA (mRNA) analysis has become a powerful tool for exploring tumor recurrence or metastasis mechanisms. Here, we constructed a signature to predict the recurrence risk of Stages II and III gastric cancer (GC) patients. A least absolute shrinkage and selection operator method Cox regression model was utilized to construct the signature. Using this method, a 16-mRNA signature was identified to be associated with the relapse-free survival of Stages II and III GCs in training dataset GSE62254 (n = 194). Then this signature was validated in an independent Gene Expression Omnibus cohort GSE26253 (n = 297) and a dataset of The Cancer Genome Atlas (TCGA; n = 235). This classifier could successfully screen out the high-risk Stages II and III GCs in the training cohort (hazard ratio [HR] = 40.91; 95% confidence interval [CI] = 5.58–299.7; p < .0001). Analysis in two independent validation cohorts yielded consistent results (GSE26253: HR = 1.69, 95% CI = 1.17–2.43,; p = .0045; TCGA: HR = 2.01, 95% CI = 1.13–3.56, p = .0146). Cox regression analyses revealed that the risk score derived from this signature was an independent risk factor in Stages II and III GCs. Besides, a nomogram was constructed to serve clinical practice. Through gene set variation analysis, we found several gene sets associated with chemotherapeutic drug resistance and tumor metastasis significantly enriched in high-risk patients. In summary, this 16-mRNA signature can be used as a powerful tool for prognostic evaluation and help clinicians identify high-risk patients.     

Combined analysis of RNA‐sequence and microarray data reveals effective metabolism‐based prognostic signature for neuroblastoma

The relationship between metabolism reprogramming and neuroblastoma (NB) is largely unknown. In this study, one RNA‐sequence data set (n = 153) was used as discovery cohort and two microarray data sets (n = 498 and n = 223) were used as validation cohorts. Differentially expressed metabolic genes were identified by comparing stage 4s and stage 4 NBs. Twelve metabolic genes were selected by LASSO regression analysis and integrated into the prognostic signature. The metabolic gene signature successfully stratifies NB patients into two risk groups and performs well in predicting survival of NB patients. The prognostic value of the metabolic gene signature is also independent with other clinical risk factors. Nine metabolism‐related long non‐coding RNAs (lncRNAs) were also identified and integrated into the metabolism‐related lncRNA signature. The lncRNA signature also performs well in predicting survival of NB patients. These results suggest that the metabolic signatures have the potential to be used for risk stratification of NB. Gene set enrichment analysis (GSEA) reveals that multiple metabolic processes (including oxidative phosphorylation and tricarboxylic acid cycle, both of which are emerging targets for cancer therapy) are enriched in the high‐risk NB group, and no metabolic process is enriched in the low‐risk NB group. This result indicates that metabolism reprogramming is associated with the progression of NB and targeting certain metabolic pathways might be a promising therapy for NB.     

Single-cell sequencing analysis and transcriptome analysis constructed the macrophage related gene-related signature in lung adenocarcinoma and verified by an independent cohort

BackgroundRecent studies have emphasized the close relationship between macrophages and tumor immunity, and the prognosis of lung adenocarcinoma (LUAD) patients is intimately linked to this. Nonetheless, the prognostic signature and classification of different immune patterns in LUAD patients based on the macrophages is largely unexplored.MethodsTwo sc-RNAseq datasets of LUAD patients were collected and reprocessed. The differentially expressed genes (DEGs) related to macrophages between LUAD tissues and normal lung tissues were then identified. Based upon the above genes, three distinct immune patterns in the TCGA-LUAD cohort were identified. The ssGSEA and CIBERSORT were applied for immune profiling and characterization of different subtypes. A four-gene prognostic signature for LUAD patients was established based on the DEGs between the subtypes using stepwise multi-Cox regression. TCGA-LUAD cohort was used as training set. Five GEO-LUAD datasets and an independent cohort containing 112 LUAD samples were used for validation. TIDE (tumor immune dysfunction and exclusion) and drug sensitivity analyses were also performed.ResultsMacrophage-related differentially expressed genes were found out using the publicly available scRNA-seq data of LUAD. Three different immune patterns which were proved to have distinct immune infiltration characteristics in the TCGA-LUAD cohort were recognized based on the above macrophage-related genes. Thereafter, 174 DEGs among the above three different immune patterns were figured out; on the basis of this, a four-gene prognostic signature was constructed. This signature distinguished the prognosis of LUAD patients well in various GSE datasets as well as our independent cohort. Further analyses revealed that patients which had a higher risk score also accompanied with a lower immune infiltration level and a worse response to several immunotherapy biomarkers.ConclusionThis study highlighted that macrophage were significantly associated with TME diversity and complexity. The four-gene prognostic signature could be used for predicting outcomes and immune landscapes for patients with LUAD.     

A panel of differentially methylated regions enable prognosis prediction for colorectal cancer

We aim to identify a panel of differentially methylated regions (DMRs) for predicting survival outcomes for patients with CRC from the TCGA (n = 393). Four DMRs (MUC12, TBX20, CHN2, and B3GNT7) were selected as candidate prognostic markers for CRC. The prediction potential of selected DMRs was validated by the targeted bisulfite sequencing method in an independent cohort with 251 Chinese CRC patients. DMR methylation scores (DMSs) were constructed to evaluate the prognosis of CRC. Results of the validation cohort confirmed that higher DMSs were associated with poor overall survival (OS) of CRC, with hazard ratio (HR) value ranged from 1.445 to 2.698 in multivariable Cox models. Patients in the high prognostic index (high-PI) group showed a markedly unfavorable prognosis compared to the low-PI group in both TCGA discovery cohort (HR = 3.508, 95%CI: 2.196–5.604, P < 0.001) and independent validation cohort (HR = 1.912, 95%CI: 1.258–2.907, P = 0.002).     

Integrative analysis from multi‐centre studies identifies a function‐derived personalized multi‐gene signature of outcome in colorectal cancer

Colorectal cancer (CRC) is highly heterogeneous leading to variable prognosis and treatment responses. Therefore, it is necessary to explore novel personalized and reproducible prognostic signatures to aid clinical decision‐making. The present study combined large‐scale gene expression profiles and clinical data of 1828 patients with CRC from multi‐centre studies and identified a personalized gene prognostic signature consisting of 46 unique genes (called function‐derived personalized gene signature [FunPGS]) from an integrated statistics and function‐derived perspective. In the meta‐training and multiple independent validation cohorts, the FunPGS effectively discriminated patients with CRC with significantly different prognosis at the individual level and remained as an independent factor upon adjusting for clinical covariates in multivariate analysis. Furthermore, the FunPGS demonstrated superior performance for risk stratification with respect to other recently reported signatures and clinical factors. The complementary value of the molecular signature and clinical factors was further explored, and it was observed that the composite signature called IMCPS greatly improved the predictive performance of survival estimation relative to molecular signatures or clinical factors alone. With further prospective validation in clinical trials, the FunPGS may become a promising and powerful personalized prognostic tool for stratifying patients with CRC in order to achieve an optimal systemic therapy.     

A 4-microRNA signature for survival prognosis in pediatric and adolescent acute myeloid leukemia

Acute myeloid leukemia (AML) is a hematologic malignancy with significant molecular heterogeneity. MicroRNAs (miRNAs) play a critical role in AML diagnosis, pathogenesis, and prognosis of AML. Little has been done to identify a miRNA signature in pediatric and adolescent patients for predicting overall survival. This study aims to identify a panel of miRNA signature that could predict the prognosis of all younger AML patients with all subtypes of AML by analyzing data from The Cancer Genome Atlas (TCGA). A total of 229 patients under 23 years with miRNA data and corresponding clinical data from TCGA database were enrolled in this study. Through conducting multivariate analysis in the training test, it was identified that the high expression of hsa-miR-509 and hsa-miR-542 were independent poor prognostic factors, whereas that of hsa-miR-146a and hsa-miR-3667 had a trend to be favorable factors. A 4-miRNA signature was constructed by these miRNAs considering the weight of each. In testing group and all 229 patients’ cohort as well as 59 cytogenetically normal AML (CN-AML) patients’ cohort, higher risk score was associated with shorter overall survival (OS). All results were confidential by using powerful statistical analysis. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis were carried out to further develop leukemia-relevant mechanisms supporting the model. The results indicate that the 4-miRNA-based signature is a reliable prognostic biomarker for pediatric and adolescent AML patients.     

Systematic profiling of invasion-related gene signature predicts prognostic features of lung adenocarcinoma

Due to the high heterogeneity of lung adenocarcinoma (LUAD), molecular subtype based on gene expression profiles is of great significance for diagnosis and prognosis prediction in patients with LUAD. Invasion-related genes were obtained from the CancerSEA database, and LUAD expression profiles were downloaded from The Cancer Genome Atlas. The ConsensusClusterPlus was used to obtain molecular subtypes based on invasion-related genes. The limma software package was used to identify differentially expressed genes (DEGs). A multi-gene risk model was constructed by Lasso-Cox analysis. A nomogram was also constructed based on risk scores and meaningful clinical features. 3 subtypes (C1, C2 and C3) based on the expression of 97 invasion-related genes were obtained. C3 had the worst prognosis. A total of 669 DEGs were identified among the subtypes. Pathway enrichment analysis results showed that the DEGs were mainly enriched in the cell cycle, DNA replication, the p53 signalling pathway and other tumour-related pathways. A 5-gene signature (KRT6A, MELTF, IRX5, MS4A1 and CRTAC1) was identified by using Lasso-Cox analysis. The training, validation and external independent cohorts proved that the model was robust and had better prediction ability than other lung cancer models. The gene expression results showed that the expression levels of MS4A1 and KRT6A in tumour tissues were higher than in normal tissues, while CRTAC1 expression in tumour tissues was lower than in normal tissues. The 5-gene signature prognostic stratification system based on invasion-related genes could be used to assess prognostic risk in patients with LUAD.     

Prognostic Significance of CSN2, CD8, and MMR Status-Associated Nomograms in Patients with Colorectal Cancer

BACKGROUND: COP9 signalosome subunit 2 (CSN2) is believed to be involved in human cancer, but its prognostic significance in colorectal cancer (CRC) has not been elucidated. PATIENTS AND METHODS: We retrospectively analyzed the expression of CSN2 andCD8+ tumor-infiltrating lymphocytes (TILs), and mismatch repair (MMR) status in 267 paraffin-embedded specimens using immunohistochemistry in a training cohort. A number of risk factors were used to form nomograms to evaluate survival, and Harrell's concordance index (C-index) was used to evaluate the predictive accuracy. Further validation was performed in an independent cohort of 238cases. RESULTS: Low CSN2 expression and a low number of CD8 + TILs were significantly associated with diminished disease-free survival (DFS) and overall survival (OS) in CRC patients, and patients with MMR-deficient CRC had enhanced DFS and OS. Moreover, the multivariate Cox analysis identified CSN2, CD8 + TILs, and MMR status as independent prognostic factors for DFS and OS. Using these three markers and four clinicopathological risk variables, two nomograms were constructed and validated for predicting DFS and OS (C-index: training cohort, 0.836 (95% CI:0.804–0.868) and 0.841 (0.808–0.874), respectively; validation cohort, 0.801 (0.760–843) and 0.843 (0.806–0.881), respectively). CONCLUSIONS: CSN2, CD8+ TILs, and MMR status were independent prognostic factors. The nomograms could be used to generate individualized predictions for DFS and OS.     

Development and verification of a nomogram for prediction of recurrence‐free survival in clear cell renal cell carcinoma

Nowadays, gene expression profiling has been widely used in screening out prognostic biomarkers in numerous kinds of carcinoma. Our studies attempt to construct a clinical nomogram which combines risk gene signature and clinical features for individual recurrent risk assessment and offer personalized managements for clear cell renal cell carcinoma. A total of 580 differentially expressed genes (DEGs) were identified via microarray. Functional analysis revealed that DEGs are of fundamental importance in ccRCC progression and metastasis. In our study, 338 ccRCC patients were retrospectively analysed and a risk gene signature which composed of 5 genes was obtained from a LASSO Cox regression model. Further analysis revealed that identified risk gene signature could usefully distinguish the patients with poor prognosis in training cohort (hazard ratio [HR] = 3.554, 95% confidence interval [CI] 2.261‐7.472, P < .0001, n = 107). Moreover, the prognostic value of this gene‐signature was independent of clinical features (P = .002). The efficacy of risk gene signature was verified in both internal and external cohorts. The area under receiver operating characteristic curve of this signature was 0.770, 0.765 and 0.774 in the training, testing and external validation cohorts, respectively. Finally, a nomogram was developed for clinicians and did well in the calibration plots. This nomogram based on risk gene signature and clinical features might provide a practical way for recurrence prediction and facilitating personalized managements of ccRCC patients after surgery.     

Ferroptosis-related gene signature correlates with the tumor immune features and predicts the prognosis of glioma patients

Background: Glioma is a malignant intracranial tumor and the most fatal cancer. The role of ferroptosis in the clinical progression of gliomas is unclear.Method: Univariate and least absolute shrinkage and selection operator (Lasso) Cox regression methods were used to develop a ferroptosis-related signature (FRSig) using a cohort of glioma patients from the Chinese Glioma Genome Atlas (CGGA), and was validated using an independent cohort of glioma patients from The Cancer Genome Atlas (TCGA). A single-sample gene set enrichment analysis (ssGSEA) was used to calculate levels of the immune infiltration. Multivariate Cox regression was used to determine the independent prognostic role of clinicopathological factors and to establish a nomogram model for clinical application.Results: We analyzed the correlations between the clinicopathological features and ferroptosis-related gene (FRG) expression and established an FRSig to calculate the risk score for individual glioma patients. Patients were stratified into two subgroups with distinct clinical outcomes. Immune cell infiltration in the glioma microenvironment and immune-related indexes were identified that significantly correlated with the FRSig, the tumor mutation burden (TMB), copy number alteration (CNA), and immune checkpoint expression was also significantly positively correlated with the FRSig score. Ultimately, an FRSig-based nomogram model was constructed using the independent prognostic factors age, World Health Organization (WHO) grade, and FRSig score.Conclusion: We established the FRSig to assess the prognosis of glioma patients. The FRSig also represented the glioma microenvironment status. Our FRSig will contribute to improve patient management and individualized therapy by offering a molecular biomarker signature for precise treatment.