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.     

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.     

Five key lncRNAs considered as prognostic targets for predicting pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) has a poor prognosis, and the 5‐year survival rate was only 7.7%. To improve prognosis, a screening biomarker for early diagnosis of pancreatic cancer is in urgent need. Long non‐coding RNA (lncRNA) expression profiles as potential cancer prognostic biomarkers play critical roles in development of tumorigenesis and metastasis of cancer. However, lncRNA signatures in predicting the survival of a patient with PDAC remain unknown. In the current study, we try to identify potential lncRNA biomarkers and their prognostic values in PDAC. LncRNAs expression profiles and corresponding clinical information for 182 cases with PDAC were acquired from The Cancer Genome Atlas (TCGA). A total of 14 470 lncRNA were identified in the cohort, and 175 PDAC patients had clinical variables. We obtained 108 differential expressed lncRNA via R packages. Univariate and multivariate Cox proportional hazards regression, lasso regression was performed to screen the potential prognostic lncRNA. Five lncRNAs have been recognized to significantly correlate with OS. We established a linear prognostic model of five lncRNA (C9orf139, MIR600HG, RP5‐965G21.4, RP11‐436K8.1, and CTC‐327F10.4) and divided patients into high‐ and low‐risk group according to the prognostic index. The five lncRNAs played independent prognostic biomarkers of OS of PDAC patients and the AUC of the ROC curve for the five lncRNAs signatures prediction 5‐year survival was 0.742. In addition, targeted genes of MIR600HG, C9orf139, and CTC‐327F10.4 were explored and functional enrichment was also conducted. These results suggested that this five‐lncRNAs signature could act as potential prognostic biomarkers in the prediction of PDAC patient's survival.     

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

Genome‐wide methylomic analyses identify prognostic epigenetic signature in lower grade glioma

Glioma is the most malignant and aggressive type of brain tumour with high heterogeneity and mortality. Although some clinicopathological factors have been identified as prognostic biomarkers, the individual variants and risk stratification in patients with lower grade glioma (LGG) have not been fully elucidated. The primary aim of this study was to identify an efficient DNA methylation combination biomarker for risk stratification and prognosis in LGG. We conducted a retrospective cohort study by analysing whole genome DNA methylation data of 646 patients with LGG from the TCGA and GEO database. Cox proportional hazard analysis was carried out to screen and construct biomarker model that predicted overall survival (OS). The Kaplan‐Meier survival curves and time‐dependent ROC were constructed to prove the efficiency of the signature. Then, another independent cohort was used to further validate the finding. A two‐CpG site DNA methylation signature was identified by multivariate Cox proportional hazard analysis. Further analysis indicated that the signature was an independent survival predictor from other clinical factors and exhibited higher predictive accuracy compared with known biomarkers. This signature was significantly correlated with immune‐checkpoint blockade, immunotherapy‐related signatures and ferroptosis regulator genes. The expression pattern and functional analysis showed that these two genes corresponding with two methylation sites contained in the model were correlated with immune infiltration level, and involved in MAPK and Rap1 signalling pathway. The signature may contribute to improve the risk stratification of patients and provide a more accurate assessment for precision medicine in the clinic.     

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.     

Genome-wide screening and cohorts validation identifying novel lncRNAs as prognostic biomarkers for clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the main subtype of renal cell carcinoma with varied prognosis. We aimed to identify and assess the possible prognostic long noncoding RNA (lncRNA) biomarkers. LncRNAs expression data and corresponding clinical information of 619 ccRCC patients were downloaded from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. Differentially expressed genes analysis, univariate Cox regression, the least absolute shrinkage and selection operator Cox regression model were utilized to identify hub lncRNAs. Multivariate Cox regression was used to establish the risk model. Statistical analysis was performed using R 3.5.3. The expression value of five lncRNAs and the risk-score levels were significantly associated with a survival prognosis of ccRCC patients (all P < .001). In the TCGA validation cohort, the area under the curve (AUC) for the integrated nomogram was 0.905 and 0.91 for 3-, 5-year prediction separately. The AUC reached up to 0.757 in an independent ICGC cohort. Besides, the calibration plots also illustrated well curve-fitting between observation values and predictive values. Weighted gene co-expression network analysis and subsequent pathway analysis revealed that the PI3K-Akt-mTOR and hypoxia-inducible factor signaling crosstalk might function as the most essential mechanisms related to the five-lncRNAs signature. Our study suggested that lncRNA AC009654.1, AC092490.2, LINC00524, LINC01234, and LINC01885 were significantly associated with ccRCC prognosis. The prognostic model based on this five lncRNA may predict the overall survival of ccRCC.     

Integrative analysis identifies potential DNA methylation biomarkers for pan-cancer diagnosis and prognosis

DNA methylation status is closely associated with diverse diseases, and is generally more stable than gene expression, thus abnormal DNA methylation could be important biomarkers for tumor diagnosis, treatment and prognosis. However, the signatures regarding DNA methylation changes for pan-cancer diagnosis and prognosis are less explored. Here we systematically analyzed the genome-wide DNA methylation patterns in diverse TCGA cancers with machine learning. We identified seven CpG sites that could effectively discriminate tumor samples from adjacent normal tissue samples for 12 main cancers of TCGA (1216 samples, AUC > 0.99). Those seven potential diagnostic biomarkers were further validated in the other 9 different TCGA cancers and 4 independent datasets (AUC > 0.92). Three out of the seven CpG sites were correlated with cell division, DNA replication and cell cycle. We also identified 12 CpG sites that can effectively distinguish 26 different cancers (7605 samples), and the result was repeatable in independent datasets as well as two disparate tumors with metastases (micro-average AUC > 0.89). Furthermore, a series of potential signatures that could significantly predict the prognosis of tumor patients for 7 different cancer were identified via survival analysis (p-value < 1e-4). Collectively, DNA methylation patterns vary greatly between tumor and adjacent normal tissues, as well as among different types of cancers. Our identified signatures may aid the decision of clinical diagnosis and prognosis for pan-cancer and the potential cancer-specific biomarkers could be used to predict the primary site of metastatic breast and prostate cancers.     

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.     

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.     

The impact on high‐grade serous ovarian cancer of obesity and lipid metabolism‐related gene expression patterns: the underestimated driving force affecting prognosis

To investigate whether specific obesity/metabolism‐related gene expression patterns affect the survival of patients with ovarian cancer. Clinical and genomic data of 590 samples from the high‐grade ovarian serous carcinoma (HGOSC) study of The Cancer Genome Atlas (TCGA) and 91 samples from the Australian Ovarian Cancer Study were downloaded from the International Cancer Genome Consortium (ICGC) portal. Clustering of mRNA microarray and reverse‐phase protein array (RPPA) data was performed with 83 consensus driver genes and 144 obesity and lipid metabolism‐related genes. Association between different clusters and survival was analyzed with the Kaplan–Meier method and a Cox regression. Mutually exclusive, co‐occurrence and network analyses were also carried out. Using RNA and RPPA data, it was possible to identify two subsets of HGOSCs with similar clinical characteristics and cancer driver mutation profiles (e.g. TP53), but with different outcome. These differences depend more on up‐regulation of specific obesity and lipid metabolism‐related genes than on the number of gene mutations or copy number alterations. It was also found that CD36 and TGF‐ß are highly up‐regulated at the protein levels in the cluster with the poorer outcome. In contrast, BSCL2 is highly up‐regulated in the cluster with better progression‐free and overall survival. Different obesity/metabolism‐related gene expression patterns constitute a risk factor for prognosis independent of the therapy results in the Cox regression. Prognoses were conditioned by the differential expression of obesity and lipid metabolism‐related genes in HGOSCs with similar cancer driver mutation profiles, independent of the initial therapeutic response.     

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.     

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

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

Identification of a novel 5-methylcytosine-related signature for prognostic prediction of kidney renal papillary cell carcinoma and a Putative target for drug repurposing

BackgroundMany studies have demonstrated the crucial roles of 5-methylcytosine (m5C) RNA methylation in cancer pathogenesis.MethodsTwo datasets, including TCGA-KIRP and ICGC, and related clinical information were downloaded, where the expression of 13 m5C regulators was examined. We applied LASSO regression to construct a multi-m5C-regulator-based signature in the TCGA cohort, which was further validated using the ICGC cohort. Univariate and multivariate Cox regressions were applied to evaluate the independent prognostic value of our model. The differences in biological functions and immune characterizations between high and low-risk groups divided based on the risk scores were also investigated via multiple approaches, such as enrichment analyses, mutation mining, and immune scoring. Finally, the sensitivities of commonly used targeted drugs were tested, and the connectivity MAP (cMAP) was utilized to screen potentially effective molecules for patients in the high-risk group. Experimental validation was done following qPCR tests in Caki-2 and HK-2 cell lines.Results3 m5C regulators, including ALYREF, DNMT3B and YBX1, were involved in our model. Survival analysis revealed a worse prognosis for patients in the high-risk group. Cox regression results indicated our model's superior predictive performance compared to single-factor prognostic evaluation. Functional enrichment analyses indicated a higher mutation frequency and poorer tumor microenvironment of patients in the high-risk group. qPCR-based results revealed that ALYREF, DNMT3B, and YBX1 were significantly up-regulated in Caki-2 cell lines compared with HK-2 cell lines. Molecules like BRD-K72451865, Levosimendan, and BRD-K03515135 were advised by cMAP for patients in the high-risk group.ConclusionOur study presented a novel predictive model for KIRP prognosis. Furthermore, the results of our analysis provide new insights for investigating m5C events in KIRP pathogenesis.     

A large-scale transcriptome analysis identified ELANE and PRTN3 as novel methylation prognostic signatures for clear cell renal cell carcinoma

DNA methylation was involved in the progress of many types of cancer including clear cell renal cell carcinomas (ccRCCs). This study aimed to identify the prognostic DNA methylation biomarkers for the ccRCCs by a large-scale RNA-seq analysis. The DNA methylation data and the corresponding clinical information of the patients with ccRCCs were extracted from TCGA database and randomly divided into the training group and the validation group. The differentially expressed CpG sites and the survival-related CpG sites were further identified, which was combined into CpG sites pair and followed by screening the survival-related pairs. The C-index and the forward search algorithms were constructed to identify the prognostic signatures for the patients with ccRCCs. The prognostic signatures were verified by the validation dataset and the protein–protein interactions (PPI) network analysis was performed on the CPG sites of the signature. A total of 9,861 differentially expressed CPG sites were identified and 567 CpG sites were found to relate to the overall survival (OS) of the patients with ccRCCs. Besides, 1,146 CPG sites pairs were found to be related to the OS of the ccRCCs samples and the signature composed of seven CpG sites pairs were obtained to predict the prognosis of patients with ccRCCs and the results were verified in the validation dataset. Besides, the PPI network analysis showed that ELANE and PRTN3 gene may be associated with the invasion and metastasis of ccRCCs and could function as potential prognostic and therapeutic signatures for ccRCCs.     

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.     

A robust six-miRNA prognostic signature for head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) remains a major health problem worldwide. We aimed to identify a robust microRNA (miRNA)-based signature for predicting HNSCC prognosis. The miRNA expression profiles of HNSCC were obtained from The Cancer Genome Atlas (TCGA) database. The TCGA HNSCC cohort was randomly divided into the discovery and validation cohort. A miRNA-based prognostic signature was built up based on TGCA discovery cohort, and then further validated. The downstream targets of prognostic miRNAs were subjected to functional enrichment analyses. The role of miR-1229-3p, a prognosis-related miRNA, in tumorigenesis of HNSCC was further evaluated. A total of 305 significantly differentially expressed miRNAs were found between HNSCC samples and normal tissues. A six-miRNA prognostic signature was constructed, which exhibited a strong association with overall survival (OS) in the TCGA discovery cohort. In addition, these findings were successfully confirmed in TCGA validation cohort and our own independent cohort. The miRNA-based signature was demonstrated as an independent prognostic indicator for HNSCC. A risk signature-based nomogram model was constructed and showed good performance for predicting the OS for HNSCC. The functional analyses revealed that the downstream targets of these prognostic miRNAs were closely linked to cancer progression. Mechanistically, in vitro analysis revealed that miR-1229-3p played a tumor promoting role in HNSCC. In conclusion, our study has developed a robust miRNA-based signature for predicting the prognosis of HNSCC with high accuracy, which will contribute to improve the therapeutic outcome.     

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.