Lineage-specific detection of residual disease predicts relapse in patients with chronic myeloid leukemia stopping therapy.

ABSTRACT: Patients with chronic myeloid leukemia who are eligible for treatment-free remission (TFR) may still relapse after tyrosine kinase inhibitor (TKI) cessation. There is a need for accurate predictors of outcome to enable patients with a favorable profile to proceed while avoiding futile attempts. Sensitive detection of residual disease in total leukocytes at treatment cessation is associated with relapse but is not highly discriminatory, likely because it is a composite measure of residual leukemia derived from different cell lineages, whereas only some lineages are relevant for relapse. We prospectively measured BCR::ABL1 DNA as a predictive yes/no binary test in 5 cellular fractions from 48 patients meeting conventional criteria for TKI discontinuation. The median BCR::ABL1 DNA level was higher in granulocytes and T cells, but not in other lineages, in patients who relapsed. Among the 40 patients undergoing their first TFR attempt, we defined 3 groups with differing relapse risk: granulocyte-positive group (100%), granulocyte-negative/T-cell-positive group (67%), and granulocyte-negative /T-cell-negative group (25%). These data show the critical importance of lineage-specific assessment of residual disease in the selection of patients who can attempt to achieve TFR with a high expectation of success and, concurrently, defer patients who have a high probability of relapse.

Immune modulation in chronic myeloid leukaemia patients treated with nilotinib and interferon-alpha.

The addition of interferon to tyrosine kinase inhibitors (TKIs), to improve deep molecular response (DMR) and potentially treatment-free remission (TFR) rates in chronic-phase chronic myeloid leukaemia (CP-CML) patients is under active investigation. However, the immunobiology of this combination is poorly understood. We performed a comprehensive longitudinal assessment of immunological changes in CML patients treated with nilotinib and interferon-alpha (IFN-α) within the ALLG CML11 trial (n = 12) or nilotinib alone (n = 17). We demonstrate that nilotinib+IFN transiently reduced absolute counts of natural killer (NK) cells, compared with nilotinib alone. Furthermore, CD16+ -cytolytic and CD57+ CD62L- -mature NK cells were transiently reduced during IFN therapy, without affecting NK-cell function. IFN transiently increased cytotoxic T-lymphocyte (CTL) responses to leukaemia-associated antigens (LAAs) proteinase-3, BMI-1 and PRAME; and had no effect on regulatory T cells, or myeloid-derived suppressor cells. Patients on nilotinib+IFN who achieved MR4.5 by 12 months had a significantly higher proportion of NK cells expressing NKp46, NKp30 and NKG2D compared with patients not achieving this milestone. This difference was not observed in the nilotinib-alone group. The addition of IFN to nilotinib drives an increase in NK-activating receptors, CTLs responding to LAAs and results in transient immune modulation, which may influence earlier DMR, and its effect on long-term outcomes warrants further investigation.

Solute Carrier Family 29A1 Mediates In Vitro Resistance to Azacitidine in Acute Myeloid Leukemia Cell Lines.

Azacitidine (AZA) is commonly used hypomethylating agent for higher risk myelodysplastic syndromes and acute myeloid leukemia (AML). Although some patients achieve remission, eventually most patients fail AZA therapy. Comprehensive analysis of intracellular uptake and retention (IUR) of carbon-labeled AZA (14C-AZA), gene expression, transporter pump activity with or without inhibitors, and cytotoxicity in naïve and resistant cell lines provided insight into the mechanism of AZA resistance. AML cell lines were exposed to increasing concentrations of AZA to create resistant clones. 14C-AZA IUR was significantly lower in MOLM-13- (1.65 ± 0.08 ng vs. 5.79 ± 0.18 ng; p < 0.0001) and SKM-1- (1.10 ± 0.08 vs. 5.08 ± 0.26 ng; p < 0.0001) resistant cells compared to respective parental cells. Importantly, 14C-AZA IUR progressively reduced with downregulation of SLC29A1 expression in MOLM-13- and SKM-1-resistant cells. Furthermore, nitrobenzyl mercaptopurine riboside, an SLC29A inhibitor, reduced 14C-AZA IUR in MOLM-13 (5.79 ± 0.18 vs. 2.07 ± 0.23, p < 0.0001) and SKM-1-naive cells (5.08 ± 2.59 vs. 1.39 ± 0.19, p = 0.0002) and reduced efficacy of AZA. As the expression of cellular efflux pumps such as ABCB1 and ABCG2 did not change in AZA-resistant cells, they are unlikely contribute to AZA resistance. Therefore, the current study provides a causal link between in vitro AZA resistance and downregulation of cellular influx transporter SLC29A1.

Outcomes following venetoclax-based treatment in therapy-related myeloid neoplasms.

Therapy-related myeloid neoplasms (t-MN) are aggressive malignancies in need of effective therapies. The BCL-2 inhibitor venetoclax represents a paradigm shift in the treatment of acute myeloid leukemia. However, the effectiveness of venetoclax has not been studied in a large cohort of t-MN. We retrospectively analyzed 378 t-MN patients, of which 96 (25.4%, 47 therapy-related acute myeloid leukemia, 1 therapy-related chronic myelomonocytic leukemia, 48 therapy-related myelodysplastic syndrome) received venetoclax. Median interval from t-MN to venetoclax initiation was 2.9 (Interquartile range [IQR] 0.7-12) months, and patients received a median of 3 (IQR 1-4) cycles. The composite complete remission (CRc) rate, median progression-free survival (PFS), and overall survival (OS) were 39.1%, 4.9 months, and 7 months, respectively. The upfront use of venetoclax and achieving CRc were associated with improved survival, whereas the presence of Chromosome 7 abnormalities was associated with an inferior survival. Neither the TP53-status nor the percent bone marrow blast predicted the likelihood of CRc or survival. Paired genetic analysis performed at venetoclax initiation and failure did not show the evidence of the selection of the TP53-mutated clone. In a propensity-matched analysis, the use of venetoclax-based regimen as the first-line therapy was associated with a superior survival compared to hypomethylating agent (HMA)-based first-line therapy (9.4 vs. 6.1 months, p = .01). We conclude that the upfront use of venetoclax with HMA improved survival, though PFS and OS remain poor. As the phenotype at diagnosis or the percent blasts did not predict outcomes, venetoclax should be studied in all t-MN phenotypes.

Successful treatment-free remission in chronic myeloid leukaemia and its association with reduced immune suppressors and increased natural killer cells.

There is currently no biomarker that reliably predicts treatment-free remission (TFR) in chronic myeloid leukaemia (CML). We characterised effector and suppressor immune responses at the time of tyrosine kinase inhibitor (TKI) cessation in patients from the CML8 and CML10 clinical studies. Natural killer (NK) cells with increased expression of activating NK receptors were higher in patients who achieved TFR. There was no difference in the proportion of CD4+ or CD8+ T cells. Furthermore, we found that FoxP3+ regulatory T cells (T reg) and monocytic myeloid-derived suppressor cells (Mo-MDSCs) were concomitantly decreased in TFR patients, suggesting that the effector and suppressor arms of the immune system work in concert to mediate TFR. A discovery cohort (CML10) was used to generate a predictive model, using logistic regression. Patients classified into the high-risk group were more likely to relapse when compared with the low-risk group (HR 7·4, 95% CI 2·9-19·1). The model was successfully validated on the independent CML8 cohort (HR 8·3, 95% CI 2·2-31·3). Effective prediction of TFR success may be obtained with an effector-suppressor score, calculated using absolute NK cell, T reg, and Mo-MDSC counts, at TKI cessation, reflecting the contribution of both immune suppressors and effectors in the immunobiology underlying successful TFR.

Pre-B acute lymphoblastic leukaemia recurrent fusion, EP300-ZNF384, is associated with a distinct gene expression.

BACKGROUND: Zinc-finger protein 384 (ZNF384) fusions are an emerging subtype of precursor B-cell acute lymphoblastic leukaemia (pre-B-ALL) and here we further characterised their prevalence, survival outcomes and transcriptome. METHODS: Bone marrow mononuclear cells from 274 BCR-ABL1-negative pre-B-ALL patients were immunophenotyped and transcriptome molecularly characterised. Transcriptomic data was analysed by principal component analysis and gene-set enrichment analysis to identify gene and pathway expression changes. RESULTS: We exclusively detect E1A-associated protein p300 (EP300)-ZNF384 in 5.7% of BCR-ABL1-negative adolescent/young adult (AYA)/adult pre-B-ALL patients. EP300-ZNF384 patients do not appear to be a high-risk subgroup. Transcriptomic analysis revealed that EP300-ZNF384 samples have a distinct gene expression profile that results in the up-regulation of Janus kinase/signal transducers and activators of transcription (JAK/STAT) and cell adhesion pathways and down-regulation of cell cycle and DNA repair pathways. CONCLUSIONS: Importantly, this report contributes to a better overview of the incidence of EP300-ZNF384 patients and show that they have a distinct gene signature with concurrent up-regulation of JAK-STAT pathway, reduced expression of B-cell regulators and reduced DNA repair capacity.

Increased peroxisome proliferator-activated receptor γ activity reduces imatinib uptake and efficacy in chronic myeloid leukemia mononuclear cells.

Imatinib is actively transported by organic cation transporter-1 (OCT-1) influx transporter, and low OCT-1 activity in diagnostic chronic myeloid leukemia blood mononuclear cells is significantly associated with poor molecular response to imatinib. Herein we report that, in diagnostic chronic myeloid leukemia mononuclear cells and BCR-ABL1+ cell lines, peroxisome proliferator-activated receptor γ agonists (GW1929, rosiglitazone, pioglitazone) significantly decrease OCT-1 activity; conversely, peroxisome proliferator-activated receptor γ antagonists (GW9662, T0070907) increase OCT-1 activity. Importantly, these effects can lead to corresponding changes in sensitivity to BCR-ABL kinase inhibition. Results were confirmed in peroxisome proliferator-activated receptor γ-transduced K562 cells. Furthermore, we identified a strong negative correlation between OCT-1 activity and peroxisome proliferator-activated receptor γ transcriptional activity in diagnostic chronic myeloid leukemia patients (n=84; P<0.0001), suggesting that peroxisome proliferator-activated receptor γ activation has a negative impact on the intracellular uptake of imatinib and consequent BCR-ABL kinase inhibition. The inter-patient variability of peroxisome proliferator-activated receptor γ activation likely accounts for the heterogeneity observed in patient OCT-1 activity at diagnosis. Recently, the peroxisome proliferator-activated receptor γ agonist pioglitazone was reported to act synergistically with imatinib, targeting the residual chronic myeloid leukemia stem cell pool. Our findings suggest that peroxisome proliferator-activated receptor γ ligands have differential effects on circulating mononuclear cells compared to stem cells. Since the effect of peroxisome proliferator-activated receptor γ activation on imatinib uptake in mononuclear cells may counteract the clinical benefit of this activation in stem cells, caution should be applied when combining these therapies, especially in patients with high peroxisome proliferator-activated receptor γ transcriptional activity.

Ecotopic viral integration site 1 (EVI1) regulates multiple cellular processes important for cancer and is a synergistic partner for FOS protein in invasive tumors.

Ecotropic viral integration site 1 (EVI1) is an oncogenic dual domain zinc finger transcription factor that plays an essential role in the regulation of hematopoietic stem cell renewal, and its overexpression in myeloid leukemia and epithelial cancers is associated with poor patient survival. Despite the discovery of EVI1 in 1988 and its emerging role as a dominant oncogene in various types of cancer, few EVI1 target genes are known. This lack of knowledge has precluded a clear understanding of exactly how EVI1 contributes to cancer. Using a combination of ChIP-Seq and microarray studies in human ovarian carcinoma cells, we show that the two zinc finger domains of EVI1 bind to DNA independently and regulate different sets of target genes. Strikingly, an enriched fraction of EVI1 target genes are cancer genes or genes associated with cancer. We also show that more than 25% of EVI1-occupied genes contain linked EVI1 and activator protein (AP)1 DNA binding sites, and this finding provides evidence for a synergistic cooperative interaction between EVI1 and the AP1 family member FOS in the regulation of cell adhesion, proliferation, and colony formation. An increased number of dual EVI1/AP1 target genes are also differentially regulated in late-stage ovarian carcinomas, further confirming the importance of the functional cooperation between EVI1 and FOS. Collectively, our data indicate that EVI1 is a multipurpose transcription factor that synergizes with FOS in invasive tumors.

The GM-CSF receptor utilizes ß-catenin and Tcf4 to specify macrophage lineage differentiation.

Granulocyte-macrophage colony stimulating factor (GM-CSF) promotes the growth, survival, differentiation and activation of normal myeloid cells and is essential for fully functional macrophage differentiation in vivo. To better understand the mechanisms by which growth factors control the balance between proliferation and self-renewal versus growth-suppression and differentiation we have used the bi-potent FDB1 myeloid cell line, which proliferates in IL-3 and differentiates to granulocytes and macrophages in response to GM-CSF. This provides a manipulable model in which to dissect the switch between growth and differentiation. We show that, in the context of signaling from an activating mutant of the GM-CSF receptor ß subunit, a single intracellular tyrosine residue (Y577) mediates the granulocyte fate decision. Loss of granulocyte differentiation in a Y577F second-site mutant is accompanied by enhanced macrophage differentiation and accumulation of ß-catenin together with activation of Tcf4 and other Wnt target genes. These include the known macrophage lineage inducer, Egr1. We show that forced expression of Tcf4 or a stabilised ß-catenin mutant is sufficient to promote macrophage differentiation in response to GM-CSF and that GM-CSF can regulate ß-catenin stability, most likely via GSK3ß. Consistent with this pathway being active in primary cells we show that inhibition of GSK3ß activity promotes the formation of macrophage colonies at the expense of granulocyte colonies in response to GM-CSF. This study therefore identifies a novel pathway through which growth factor receptor signaling can interact with transcriptional regulators to influence lineage choice during myeloid differentiation.

The acquisition order of leukemic drug resistance mutations is directed by the selective fitness associated with each resistance mechanism.

In Chronic Myeloid Leukemia, the transition from drug sensitive to drug resistant disease is poorly understood. Here, we used exploratory sequencing of gene transcripts to determine the mechanisms of drug resistance in a dasatinib resistant cell line model. Importantly, cell samples were collected sequentially during drug exposure and dose escalation, revealing several resistance mechanisms which fluctuated over time. BCR::ABL1 overexpression, BCR::ABL1 kinase domain mutation, and overexpression of the small molecule transporter ABCG2, were identified as dasatinib resistance mechanisms. The acquisition of mutations followed an order corresponding with the increase in selective fitness associated with each resistance mechanism. Additionally, it was demonstrated that ABCG2 overexpression confers partial ponatinib resistance. The results of this study have broad applicability and help direct effective therapeutic drug usage and dosing regimens and may be useful for clinicians to select the most efficacious therapy at the most beneficial time.

Adverse outcomes for chronic myeloid leukemia patients with splenomegaly and low in vivo kinase inhibition on imatinib.

Variability in the molecular response to frontline tyrosine kinase inhibitor (TKI) therapy in chronic myeloid leukemia may be partially driven by differences in the level of kinase inhibition induced. We measured in vivo BCR::ABL1 kinase inhibition (IVKI) in circulating mononuclear cells after 7 days of therapy. In 173 patients on imatinib 600 mg/day, 23% had low IVKI (<11% reduction in kinase activity from baseline); this was associated with higher rates of early molecular response (EMR) failure; lower rates of major molecular response (MMR), and MR4.5 by 36 months, compared to high IVKI patients. Low IVKI was more common (39%) in patients with large spleens (≥10 cm by palpation). Notably 55% of patients with large spleens and low IVKI experienced EMR failure whereas the EMR failure rate in patients with large spleens and high IVKI was only 12% (p = 0.014). Furthermore, patients with large spleen and low IVKI had a higher incidence of blast crisis, inferior MMR, MR4.5, and event-free survival compared to patients with large spleen and high IVKI and remaining patients. In nilotinib-treated patients (n = 73), only 4% had low IVKI. The combination of low IVKI and large spleen is associated with markedly inferior outcomes and interventions in this setting warrant further studies.

A Personalized Risk Model for Azacitidine Outcome in Myelodysplastic Syndrome and Other Myeloid Neoplasms Identified by Machine Learning Model Utilizing Real-World Data.

Azacitidine is an approved therapy for higher-risk myelodysplastic syndrome (MDS). However, only 30-40% patients respond to azacitidine, and the responses may take up to six cycles to become evident. Delayed responses and the myelosuppressive effects of azacitidine make it challenging to predict which patients will benefit. This is further compounded by a lack of uniform prognostic tools to identify patients at risk of early treatment failure. Hence, we performed a retrospective analysis of 273 consecutive azacytidine-treated patients. The median overall survival was 16.25 months with only 9% alive at 5 years. By using pre-treatment variables incorporated into a random forest machine learning model, we successfully identified those patients unlikely to benefit from azacytidine upfront (7.99 vs. 22.8 months, p < 0.0001). This model also identified those who required significantly more hospitalizations and transfusion support. Notably, it accurately predicted survival outcomes, outperforming the existing prognostic scoring system. By integrating somatic mutations, we further refined the model and identified three distinct risk groups with significant differences in survival (5.6 vs. 10.5 vs. 43.5 months, p < 0.0001). These real-world findings emphasize the urgent need for personalized prediction tools tailored to hypomethylating agents, reducing unnecessary complications and resource utilization in MDS treatment.

Association of TIM-3 checkpoint receptor expression on T cells with treatment-free remission in chronic myeloid leukemia.

Dysregulation of immune-checkpoint receptors has been reported at diagnosis of chronic myeloid leukemia (CML), however, their role in the maintenance of remission after tyrosine kinase inhibitor (TKI) cessation is unclear. We assessed programmed cell death-1 (PD-1), T-cell immunoglobulin, and mucin-domain containing protein-3 (TIM-3), cytotoxic T-lymphocyte-associated protein-4 (CTLA-4), lymphocyte-activation gene-3 (LAG-3), and T-cell immunoreceptor with immunoglobulin and immunoreceptor tyrosine-based inhibitory motif (ITIM) domains (TIGIT) expression on T-cell subsets, regulatory T cells (T-regs), and natural killer (NK) cells at the time of TKI cessation in 44 patients (22 patients sustained treatment-free remission [TFR] and 22 experienced molecular relapse [MolR]). We confirmed our previous finding that absolute numbers of T-regs are increased in patients who experienced MolR compared with those who sustained TFR. The immune-checkpoint receptors PD-1, CTLA-4, LAG-3, and TIGIT on T or NK cells were not differentially expressed between the MolR and TFR groups. However, TIM-3 was consistently upregulated on bulk T cells (CD3+) and T-cell subsets including, CD4+ T cells, CD8+ T cells, and T-regs, in patients who relapsed in comparison with those who maintained TFR after discontinuation. Furthermore, gene expression analysis from publicly available data sets showed increased TIM-3 expression on CML stem cells compared with normal hematopoietic stem cells. These findings suggest that among the targetable immune-checkpoint molecules, TIM-3 blockade may potentially improve effector immune response in patients with CML stopping TKI, while concomitantly targeting leukemic stem cells and could be a promising therapeutic strategy for preventing relapse after cessation of TKI in patients with CML.

Factors predicting survival following alloSCT in patients with therapy-related AML and MDS: a multicenter study.

Therapy-related myeloid neoplasms (t-MN) are aggressive myeloid neoplasms. Factors predicting post-allogeneic stem cell transplant (alloSCT) survival are not well-known. We studied the prognostic utility of factors at: t-MN diagnosis, pre-alloSCT, and post-alloSCT. Primary endpoints were 3-year overall survival (OS), relapse incidence (RI), and non-relapse mortality (NRM). Post-alloSCT OS did not differ between t-MDS and t-AML (20.1 vs. 19.6 months, P = 1), though t-MDS had a significantly higher 3-year RI compared to t-AML (45.1% vs. 26.9%, P = 0.03). In t-MDS, the presence of monosomy 5 (HR 3.63, P = 0.006) or monosomy 17 (HR 11.81, P = 0.01) pre-alloSCT were associated with higher RI. Complex karyotype was the only factor adversely influencing survival at all the timepoints. The inclusion of genetic information yielded 2 risk-categories: high-risk defined by the presence of pathogenic variants (PV) in (TP53/BCOR/IDH1/GATA2/BCORL1) and standard-risk (remainder of the patients) with 3-year post-alloSCT OS of 0% and 64.6%, respectively (P = 0.001). We concluded that while alloSCT was curative in a subset of t-MN patients, outcomes remained poor, specifically in the high-risk category. t-MDS patients, especially those with persistent disease pre-alloSCT were at increased risk of relapse. Disease-related factors at t-MN diagnosis were the most prognostic of post-alloSCT survival; utility of factors available later in the course, was incremental.

RNA-Based Targeted Gene Sequencing Improves the Diagnostic Yield of Mutant Detection in Chronic Myeloid Leukemia.

Mutation detection is increasingly used for the management of hematological malignancies. Prior whole transcriptome and whole exome sequencing studies using total RNA and DNA identified diverse mutation types in cancer-related genes associated with treatment failure in patients with chronic myeloid leukemia. Variants included single-nucleotide variants and small insertions/deletions, plus fusion transcripts and partial or whole gene deletions. The hypothesis that all of these mutation types could be detected by a single cost-effective hybridization capture next-generation sequencing method using total RNA was assessed. A method was developed that targeted 130 genes relevant for myeloid and lymphoid leukemia. Retrospective samples with 121 precharacterized variants were tested using total RNA and/or DNA. Concordance of detection of precharacterized variants using RNA or DNA was 96%, whereas the enhanced sensitivity identified additional variants. Comparison between 24 matched DNA and RNA samples demonstrated 95.3% of 170 variants detectable using DNA were detected using RNA, including all but one variant predicted to activate nonsense-mediated decay. RNA identified an additional 10 variants, including fusion transcripts. Furthermore, the true effect of splice variants on RNA splicing was only evident using RNA. In conclusion, capture sequencing using total RNA alone is suitable for detecting a range of variants relevant in chronic myeloid leukemia and may be more broadly applied to other hematological malignancies where diverse variant types define risk groups.

Expression profiling of a hemopoietic cell survival transcriptome implicates osteopontin as a functional prognostic factor in AML.

Deregulated cell survival programs are a classic hallmark of cancer. We have previously identified a serine residue (Ser585) in the betac subunit of the granulocyte-macrophage colony-stimulating factor receptor that selectively and independently promotes cell survival. We now show that Ser585 phosphorylation is constitutive in 20 (87%) of 23 acute myeloid leukemia (AML) patient samples, indicating that this survival-only pathway is frequently deregulated in leukemia. We performed a global expression screen to identify gene targets of this survival pathway and report a 138-gene betac Ser585-regulated transcriptome. Pathway analysis defines a gene network enriched for PI3-kinase target genes and a cluster of genes involved in cancer and cell survival. We show that one such gene, osteopontin (OPN), is a functionally relevant target of the Ser585-survival pathway as shown by siRNA-mediated knockdown of OPN expression that induces cell death in both AML blasts and CD34(+)CD38(-)CD123(+) leukemic progenitors. Increased expression of OPN at diagnosis is associated with poor prognosis with multivariate analysis indicating that it is an independent predictor of overall patient survival in normal karyotype AML (n = 60; HR = 2.2; P = .01). These results delineate a novel cytokine-regulated Ser585/PI3-kinase signaling network that is deregulated in AML and identify OPN as a potential prognostic and therapeutic target.

Polycomb Factor PHF19 Controls Cell Growth and Differentiation Toward Erythroid Pathway in Chronic Myeloid Leukemia Cells.

Polycomb group (PcG) of proteins are a group of highly conserved epigenetic regulators involved in many biological functions, such as embryonic development, cell proliferation, and adult stem cell determination. PHD finger protein 19 (PHF19) is an associated factor of Polycomb repressor complex 2 (PRC2), often upregulated in human cancers. In particular, myeloid leukemia cell lines show increased levels of PHF19, yet little is known about its function. Here, we have characterized the role of PHF19 in myeloid leukemia cells. We demonstrated that PHF19 depletion decreases cell proliferation and promotes chronic myeloid leukemia (CML) differentiation. Mechanistically, we have shown how PHF19 regulates the proliferation of CML through a direct regulation of the cell cycle inhibitor p21. Furthermore, we observed that MTF2, a PHF19 homolog, partially compensates for PHF19 depletion in a subset of target genes, instructing specific erythroid differentiation. Taken together, our results show that PHF19 is a key transcriptional regulator for cell fate determination and could be a potential therapeutic target for myeloid leukemia treatment.

Integrated Bioinformatics Analysis Reveals Key Candidate Genes and Pathways Associated With Clinical Outcome in Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) accounts for approximately 85-90% of all liver cancer cases and has poor relapse-free survival. There are many gene expression studies that have been performed to elucidate the genetic landscape and driver pathways leading to HCC. However, existing studies have been limited by the sample size and thus the pathogenesis of HCC is still unclear. In this study, we performed an integrated characterization using four independent datasets including 320 HCC samples and 270 normal liver tissues to identify the candidate genes and pathways in the progression of HCC. A total of 89 consistent differentially expression genes (DEGs) were identified. Gene-set enrichment analysis revealed that these genes were significantly enriched for cellular response to zinc ion in biological process group, collagen trimer in the cellular component group, extracellular matrix (ECM) structural constituent conferring tensile strength in the molecular function group, protein digestion and absorption, mineral absorption and ECM-receptor interaction. Network system biology based on the protein-protein interaction (PPI) network was also performed to identify the most connected and important genes based on our DEGs. The top five hub genes including osteopontin (SPP1), Collagen alpha-2(I) chain (COL1A2), Insulin-like growth factor I (IGF1), lipoprotein A (LPA), and Galectin-3 (LGALS3) were identified. Western blot and immunohistochemistry analysis were employed to verify the differential protein expression of hub genes in HCC patients. More importantly, we identified that these five hub genes were significantly associated with poor disease-free survival and overall survival. In summary, we have identified a potential clinical significance of these genes as prognostic biomarkers for HCC patients who would benefit from experimental approaches to obtain optimal outcome.