A multimodal atlas of hepatocellular carcinoma reveals convergent evolutionary paths and 'bad apple' effect on clinical trajectory.

BACKGROUND & AIMS: Hepatocellular Carcinoma (HCC) is a highly fatal cancer characterized by high intra-tumor heterogeneity (ITH). A panoramic understanding of its tumor evolution, in relation to its clinical trajectory, may provide novel prognostic and treatment strategies. METHODS: Through the Asia-Pacific Hepatocellular Carcinoma (AHCC) trials group (NCT03267641), we recruited one of the largest prospective cohorts of HCC with over 600 whole genome and transcriptome samples from 123 treatment-naïve patients. RESULTS: Using a multi-region sampling approach, we revealed seven convergent genetic evolutionary paths governed by the early driver mutations, late copy number variations and viral integrations, which stratify patient clinical trajectories after surgical resection. Furthermore, such evolutionary paths shaped the molecular profiles, leading to distinct transcriptomic subtypes. Most significantly, although we found the coexistence of multiple transcriptomic subtypes within certain tumors, patient prognosis was best predicted by the most aggressive cell fraction of the tumor, rather than by overall degree of transcriptomic ITH level - a phenomenon we termed the 'bad apple' effect. Finally, we found that characteristics throughout early and late tumor evolution provide significant and complementary prognostic power in predicting patient survival. CONCLUSIONS: Taken together, our study generated a comprehensive landscape of evolutionary history for HCC and provided a rich multi-omics resource for understanding tumor heterogeneity and clinical trajectories. CLINICAL TRIAL NUMBER: NCT03267641 (Observational cohort) IMPACT AND IMPLICATIONS: This prospective study, utilizing comprehensive multi-sector, multi-omics sequencing and clinical data from surgically resected HCC, reveals critical insights into the role of tumor evolution and intra-tumor heterogeneity (ITH) in determining the prognosis of Hepatocellular Carcinoma (HCC). These findings are invaluable for oncology researchers and clinicians, as they underscore the influence of distinct evolutionary paths and the 'bad apple' effect, where the most aggressive tumor fraction dictates disease progression. These insights not only enhance prognostic accuracy post-surgical resection but also pave the way for developing personalized therapies tailored to specific tumor evolutionary and transcriptomic profiles. The co-existence of multiple sub-types within the same tumor prompts a re-appraisal of the utilities of depending on single samples to represent the entire tumor and suggests the need for clinical molecular imaging. This research thus marks a significant step forward in the clinical understanding and management of HCC, underscoring the importance of integrating tumor evolutionary dynamics and multi-omics biomarkers into therapeutic decision-making.

Hypoxia-driven immunosuppression by Treg and type-2 conventional dendritic cells in HCC.

BACKGROUND AND AIMS: Hypoxia is one of the central players in shaping the immune context of the tumor microenvironment (TME). However, the complex interplay between immune cell infiltrates within the hypoxic TME of HCC remains to be elucidated. APPROACH AND RESULTS: We analyzed the immune landscapes of hypoxia-low and hypoxia-high tumor regions using cytometry by time of light, immunohistochemistry, and transcriptomic analyses. The mechanisms of immunosuppression in immune subsets of interest were further explored using in vitro hypoxia assays. Regulatory T cells (Tregs) and a number of immunosuppressive myeloid subsets, including M2 macrophages and human leukocyte antigen-DR isotype (HLA-DRlo ) type 2 conventional dendritic cell (cDC2), were found to be significantly enriched in hypoxia-high tumor regions. On the other hand, the abundance of active granzyme Bhi PD-1lo CD8+ T cells in hypoxia-low tumor regions implied a relatively active immune landscape compared with hypoxia-high regions. The up-regulation of cancer-associated genes in the tumor tissues and immunosuppressive genes in the tumor-infiltrating leukocytes supported a highly pro-tumorigenic network in hypoxic HCC. Chemokine genes such as CCL20 (C-C motif chemokine ligand 20) and CXCL5 (C-X-C motif chemokine ligand 5) were associated with recruitment of both Tregs and HLA-DRlo cDC2 to hypoxia-high microenvironments. The interaction between Tregs and cDC2 under a hypoxic TME resulted in a loss of antigen-presenting HLA-DR on cDC2. CONCLUSIONS: We uncovered the unique immunosuppressive landscapes and identified key immune subsets enriched in hypoxic HCC. In particular, we identified a potential Treg-mediated immunosuppression through interaction with a cDC2 subset in HCC that could be exploited for immunotherapies.

Multi-region sampling with paired sample sequencing analyses reveals sub-groups of patients with novel patient-specific dysregulation in Hepatocellular Carcinoma.

BACKGROUND: Conventional differential expression (DE) testing compares the grouped mean value of tumour samples to the grouped mean value of the normal samples, and may miss out dysregulated genes in small subgroup of patients. This is especially so for highly heterogeneous cancer like Hepatocellular Carcinoma (HCC). METHODS: Using multi-region sampled RNA-seq data of 90 patients, we performed patient-specific differential expression testing, together with the patients' matched adjacent normal samples. RESULTS: Comparing the results from conventional DE analysis and patient-specific DE analyses, we show that the conventional DE analysis omits some genes due to high inter-individual variability present in both tumour and normal tissues. Dysregulated genes shared in small subgroup of patients were useful in stratifying patients, and presented differential prognosis. We also showed that the target genes of some of the current targeted agents used in HCC exhibited highly individualistic dysregulation pattern, which may explain the poor response rate. DISCUSSION/CONCLUSION: Our results highlight the importance of identifying patient-specific DE genes, with its potential to provide clinically valuable insights into patient subgroups for applications in precision medicine.

IFNγ-IL-17+ CD8 T cells contribute to immunosuppression and tumor progression in human hepatocellular carcinoma.

IL-17-producing CD8 (Tc17) T cells have been shown to play an important role in infection and chronic inflammation, however their implications in hepatocellular carcinoma (HCC) remain elusive. In this study, we performed cytometry by time-of-flight (CyTOF) and revealed the distinctive immunological phenotypes of two IFNγ+ and IFNγ- Tc17 subsets that were preferentially enriched in human HCC. Single-cell RNA-sequencing analysis further revealed regulatory circuits governing the different phenotypes of these Tc17 subsets. In particular, we discovered that IFNγ- Tc17 subset demonstrated pro-tumoral characteristics and expressed higher levels of CCL20. This corresponded to increased tumor infiltration of T regulatory cells (Treg) validated by immunohistochemistry in another independent HCC cohort, demonstrating the immunosuppressive functions of IFNγ- Tc17 subset. Most importantly, higher intra-tumoral proportions of IFNγ- Tc17 were associated with poorer prognosis in patients with HCC and this was further validated in The Cancer Genome Atlas (TCGA) HCC cohort. Taken together, this compendium of transcriptomic and proteomic data of Tc17 subsets sheds light on the immunosuppressive phenotypes of IFNγ- Tc17 and its implications in HCC progression.

A genomic enhancer signature associates with hepatocellular carcinoma prognosis.

Background & Aims: Lifestyle and environmental-related exposures are important risk factors for hepatocellular carcinoma (HCC), suggesting that epigenetic dysregulation significantly underpins HCC. We profiled 30 surgically resected tumours and the matched adjacent normal tissues to understand the aberrant epigenetic events associated with HCC. Methods: We identified tumour differential enhancers and the associated genes by analysing H3K27 acetylation (H3K27ac) chromatin immunoprecipitation sequencing (ChIP-seq) and Hi-C/HiChIP data from the resected tumour samples of 30 patients with early-stage HCC. This epigenome dataset was analysed with previously reported genome and transcriptome data of the overlapping group of patients from the same cohort. We performed patient-specific differential expression testing using multiregion sequencing data to identify genes that undergo both enhancer and gene expression changes. Based on the genes selected, we identified two patient groups and performed a recurrence-free survival analysis. Results: We observed large-scale changes in the enhancer distribution between HCC tumours and the adjacent normal samples. Many of the gain-in-tumour enhancers showed corresponding upregulation of the associated genes and vice versa, but much of the enhancer and gene expression changes were patient-specific. A subset of the upregulated genes was activated in a subgroup of patients' tumours. Recurrence-free survival analysis revealed that the patients with a more robust upregulation of those genes showed a worse prognosis. Conclusions: We report the genomic enhancer signature associated with differential prognosis in HCC. Findings that cohere with oncofoetal reprogramming in HCC were underpinned by genome-wide enhancer rewiring. Our results present the epigenetic changes in HCC that offer the rational selection of epigenetic-driven gene targets for therapeutic intervention or disease prognostication in HCC. Impact and Implications: Lifestyle and environmental-related exposures are the important risk factors of hepatocellular carcinoma (HCC), suggesting that tumour-associated epigenetic dysregulations may significantly underpin HCC. We profiled tumour tissues and their matched normal from 30 patients with early-stage HCC to study the dysregulated epigenetic changes associated with HCC. By also analysing the patients' RNA-seq and clinical data, we found the signature genes - with epigenetic and transcriptomic dysregulation - associated with worse prognosis. Our findings suggest that systemic approaches are needed to consider the surrounding cellular environmental and epigenetic changes in HCC tumours.

Genome instability is associated with ethnic differences between Asians and Europeans in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is one of the deadliest cancer types with diverse etiological factors across the world. Although large scale genomic studies have been conducted in different countries, integrative analysis of HCC genomes and ethnic comparison across cohorts are lacking. Methods: We first integrated genomes of 1,349 HCC patients from five large cohorts across the world and applied multiple statistical methods in identifying driver genes. Subsequently, we systematically compared HCC genomes and transcriptomes between Asians and Europeans using the TCGA cohort. Results: We identified 29 novel candidate driver genes, many of which are infrequent tumor suppressors driving late-stage tumor progression. When we systematically compared ethnic differences in the genomic landscape between Asian and European HCCs using the TCGA cohort (n = 348), we found little differences in driver frequencies. Through multi-modal integrative analysis, we found higher genomic instability in Asians together with a collection of molecular events ranging from tumor mutation burden (TMB), copy number alterations as well as transcriptomic subtypes segregating distinctively between two ethnic backgrounds. Strikingly, we identified an Asian specific transcriptomic subtype with multiple ethnically enriched genomic alterations, in particular chromosome 16 deletion, leading to a clinically aggressive RNA subgroup unique to Asians. Integrating multi-modal information, we found that survival models predict patient prognosis much better in Asians than in Europeans, demonstrating a higher potential for precision medicine applications in Asia. Conclusion: For the first time, we have uncovered an unprecedented amount of genomic differences segregating distinctively across ethnicities in HCC and highlighted the importance of differential disease biology and management in HCC across ethnic backgrounds.

Dynamic phenotypic heterogeneity and the evolution of multiple RNA subtypes in hepatocellular carcinoma: the PLANET study.

Intra-tumor heterogeneity (ITH) is a key challenge in cancer treatment, but previous studies have focused mainly on the genomic alterations without exploring phenotypic (transcriptomic and immune) heterogeneity. Using one of the largest prospective surgical cohorts for hepatocellular carcinoma (HCC) with multi-region sampling, we sequenced whole genomes and paired transcriptomes from 67 HCC patients (331 samples). We found that while genomic ITH was rather constant across stages, phenotypic ITH had a very different trajectory and quickly diversified in stage II patients. Most strikingly, 30% of patients were found to contain more than one transcriptomic subtype within a single tumor. Such phenotypic ITH was found to be much more informative in predicting patient survival than genomic ITH and explains the poor efficacy of single-target systemic therapies in HCC. Taken together, we not only revealed an unprecedentedly dynamic landscape of phenotypic heterogeneity in HCC, but also highlighted the importance of studying phenotypic evolution across cancer types.

Trajectory of immune evasion and cancer progression in hepatocellular carcinoma.

Immune evasion is key to cancer initiation and later at metastasis, but its dynamics at intermediate stages, where potential therapeutic interventions could be applied, is undefined. Here we show, using multi-dimensional analyses of resected tumours, their adjacent non-tumour tissues and peripheral blood, that extensive immune remodelling takes place in patients with stage I to III hepatocellular carcinoma (HCC). We demonstrate the depletion of anti-tumoural immune subsets and accumulation of immunosuppressive or exhausted subsets along with reduced tumour infiltration of CD8 T cells peaking at stage II tumours. Corresponding transcriptomic modification occur in the genes related to antigen presentation, immune responses, and chemotaxis. The progressive immune evasion is validated in a murine model of HCC. Our results show evidence of ongoing tumour-immune co-evolution during HCC progression and offer insights into potential interventions to reverse, prevent or limit the progression of the disease.

Intratumoural immune heterogeneity as a hallmark of tumour evolution and progression in hepatocellular carcinoma.

The clinical relevance of immune landscape intratumoural heterogeneity (immune-ITH) and its role in tumour evolution remain largely unexplored. Here, we uncover significant spatial and phenotypic immune-ITH from multiple tumour sectors and decipher its relationship with tumour evolution and disease progression in hepatocellular carcinomas (HCC). Immune-ITH is associated with tumour transcriptomic-ITH, mutational burden and distinct immune microenvironments. Tumours with low immune-ITH experience higher immunoselective pressure and escape via loss of heterozygosity in human leukocyte antigens and immunoediting. Instead, the tumours with high immune-ITH evolve to a more immunosuppressive/exhausted microenvironment. This gradient of immune pressure along with immune-ITH represents a hallmark of tumour evolution, which is closely linked to the transcriptome-immune networks contributing to disease progression and immune inactivation. Remarkably, high immune-ITH and its transcriptomic signature are predictive for worse clinical outcome in HCC patients. This in-depth investigation of ITH provides evidence on tumour-immune co-evolution along HCC progression.