Non-terminally exhausted tumor-resident memory HBV-specific T cell responses correlate with relapse-free survival in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) often develops following chronic hepatitis B virus (HBV) infection and responds poorly to immune checkpoint blockade. Here, we examined the antigen specificities of HCC-infiltrating T cells and their relevance to tumor control. Using highly multiplexed peptide-MHC tetramer staining of unexpanded cells from blood, liver, and tumor tissues from 46 HCC patients, we detected 91 different antigen-specific CD8+ T cell populations targeting HBV, neoantigen, tumor-associated, and disease-unrelated antigens. Parallel high-dimensional analysis delineated five distinct antigen-specific tissue-resident memory T (Trm) cell populations. Intratumoral and intrahepatic HBV-specific T cells were enriched for two Trm cell subsets that were PD-1loTOXlo, despite being clonally expanded. High frequencies of intratumoral terminally exhausted T cells were uncommon. Patients with tumor-infiltrating HBV-specific CD8+ Trm cells exhibited longer-term relapse-free survival. Thus, non-terminally exhausted HBV-specific CD8+ Trm cells show hallmarks of active involvement and effective antitumor response, implying that these cells could be harnessed for therapeutic purposes.

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 trials group (NCT03267641), we recruited one of the largest prospective cohorts of patients with 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 provides a rich multi-omics resource for understanding tumor heterogeneity and clinical trajectories. IMPACT AND IMPLICATIONS: This prospective study, utilizing comprehensive multi-sector, multi-omics sequencing and clinical data from surgically resected hepatocellular carcinoma (HCC), reveals critical insights into the role of tumor evolution and intra-tumor heterogeneity (ITH) in determining the prognosis of 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 personalized treatment strategies tailored to specific tumor evolutionary and transcriptomic profiles. The coexistence of multiple subtypes 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. CLINICAL TRIAL NUMBER: NCT03267641 (Observational cohort).

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.

Epigenetic promoter alterations in GI tumour immune-editing and resistance to immune checkpoint inhibition.

OBJECTIVES: Epigenomic alterations in cancer interact with the immune microenvironment to dictate tumour evolution and therapeutic response. We aimed to study the regulation of the tumour immune microenvironment through epigenetic alternate promoter use in gastric cancer and to expand our findings to other gastrointestinal tumours. DESIGN: Alternate promoter burden (APB) was quantified using a novel bioinformatic algorithm (proActiv) to infer promoter activity from short-read RNA sequencing and samples categorised into APBhigh, APBint and APBlow. Single-cell RNA sequencing was performed to analyse the intratumour immune microenvironment. A humanised mouse cancer in vivo model was used to explore dynamic temporal interactions between tumour kinetics, alternate promoter usage and the human immune system. Multiple cohorts of gastrointestinal tumours treated with immunotherapy were assessed for correlation between APB and treatment outcomes. RESULTS: APBhigh gastric cancer tumours expressed decreased levels of T-cell cytolytic activity and exhibited signatures of immune depletion. Single-cell RNAsequencing analysis confirmed distinct immunological populations and lower T-cell proportions in APBhigh tumours. Functional in vivo studies using 'humanised mice' harbouring an active human immune system revealed distinct temporal relationships between APB and tumour growth, with APBhigh tumours having almost no human T-cell infiltration. Analysis of immunotherapy-treated patients with GI cancer confirmed resistance of APBhigh tumours to immune checkpoint inhibition. APBhigh gastric cancer exhibited significantly poorer progression-free survival compared with APBlow (median 55 days vs 121 days, HR 0.40, 95% CI 0.18 to 0.93, p=0.032). CONCLUSION: These findings demonstrate an association between alternate promoter use and the tumour microenvironment, leading to immune evasion and immunotherapy resistance.

Uncoupling immune trajectories of response and adverse events from anti-PD-1 immunotherapy in hepatocellular carcinoma.

BACKGROUND & AIMS: While immune checkpoint blockade (ICB) has shown promise in patients with hepatocellular carcinoma (HCC), it is associated with modest response rates and immune-related adverse events (irAEs) are common. In this study, we aimed to decipher immune trajectories and mechanisms of response and/or irAEs in patients with HCC receiving anti-programmed cell death 1 (anti-PD-1) therapy. METHODS: Pre- and on-treatment peripheral blood samples (n = 60) obtained from 32 patients with HCC (Singapore cohort) were analysed by cytometry by time-of-flight and single-cell RNA sequencing, with flow cytometric validation in an independent Korean cohort (n = 29). Mechanistic validation was conducted by bulk RNA sequencing of 20 pre- and on-treatment tumour biopsies and using a murine HCC model treated with different immunotherapeutic combinations. RESULTS: Single-cell analyses identified CXCR3+CD8+ effector memory T (TEM) cells and CD11c+ antigen-presenting cells (APC) as associated with response (p = 0.0004 and 0.0255, respectively), progression-free survival (p = 0.00079 and 0.0015, respectively), and irAEs (p = 0.0034 and 0.0125, respectively) in anti-PD-1-treated patients with HCC. Type-1 conventional dendritic cells were identified as the specific APC associated with response, while 2 immunosuppressive CD14+ myeloid clusters were linked to reduced irAEs. Further analyses of CXCR3+CD8+ TEM cells showed cell-cell interactions specific to response vs. irAEs, from which the anti-PD-1 and anti-TNFR2 combination was harnessed to uncouple these effects, resulting in enhanced response without increased irAEs in a murine HCC model. CONCLUSIONS: This study identifies early predictors of clinical response to anti-PD-1 ICB in patients with HCC and offers mechanistic insights into the immune trajectories of these immune subsets at the interface between response and toxicity. We also propose a new combination immunotherapy for HCC to enhance response without exacerbating irAEs. CLINICAL TRIAL NUMBER: NCT03695952. LAY SUMMARY: Response rates to immune checkpoint blockade (ICB) treatment in hepatocellular carcinoma (HCC) remain modest and adverse events are common. Herein, we identified early predictors of response and gained an in-depth understanding of the immunological mechanisms behind response and adverse events in patients with HCC treated with ICB. We also proposed a new combination immunotherapy for HCC that enhances response without exacerbating adverse events.

A metagenomics-based workflow for the detection and genomic characterization of GBS in raw freshwater fish.

The unexpected foodborne outbreak in Singapore in 2015 has accentuated Group B Streptococcus (GBS, Streptococcus agalactiae) sequence type 283 as an emerging foodborne pathogen transmitted via the consumption of contaminated raw freshwater fish. Isolation-based workflows utilizing conventional microbiological and whole-genome sequencing methods are commonly used to support biosurveillance efforts critical for the control management of this emerging foodborne pathogen. However, these isolation-based workflows tend to have relatively long turnaround times that hamper a timely response for implementing risk mitigation. To address this gap, we have developed a metagenomics-based workflow for the simultaneous detection and genomic characterization of GBS in raw freshwater fish. Notably, our validation results showed that this metagenomics-based workflow could achieve comparable accuracy and potentially better detection limits while halving the turnaround time (from 2 weeks to 5 days) relative to an isolation-based workflow. The metagenomics-based workflow was also successfully adapted for use on a portable long-read nanopore sequencer, demonstrating its potential applicability for real-time point-of-need testing. Using GBS in freshwater fish as an example, this work represents a proof-of-concept study that supports the feasibility and validity of metagenomics as a rapid and accurate test methodology for the detection and genomic characterization of foodborne pathogens in complex food matrices. IMPORTANCE: The need for a rapid and accurate food microbiological testing method is apparent for a timely and effective foodborne outbreak response. This is particularly relevant for emerging foodborne pathogens such as Group B Streptococcus (GBS) whose associated food safety risk might be undercharacterized. By using GBS in raw freshwater fish as a case example, this study describes the development of a metagenomics-based workflow for rapid food microbiological safety testing and surveillance. This study can inform as a working model for various foodborne pathogens in other complex food matrices, paving the way for future methodological development of metagenomics for food microbiological safety testing.

Sequence Types and Antimicrobial Resistance Profiles of Salmonella Typhimurium in the Food Chain in Singapore.

Salmonella remains a significant foodborne pathogen globally with S. Typhimurium presenting as a frequently occurring serovar. This study aimed to characterize 67 S. Typhimurium isolates from humans, food, farms, and slaughterhouses collected in Singapore from 2016 to 2017. Using whole-genome sequencing analysis, the isolates were found to belong to either ST19 (n = 33) or ST36 (n = 34). ST36 predominated in human intestinal and chicken isolates, while human extra-intestinal and non-chicken food isolates belonged to ST19. Plasmids were predicted in 88.1% (n = 59) of the isolates with the most common incompatibility group profiles being IncFIB(S), IncFII(S) and IncQ1. IncFIB(S) (adjusted p-value < 0.05) and IncFII(S) (adjusted p-value < 0.05) were significantly more prevalent in ST19 isolates, while Col156 (adjusted p-value < 0.05) was more significantly found in ST36 isolates. ST36 isolates exhibited higher resistance to multiple antibiotic classes such as penicillins, phenicols, folate pathway inhibitors, aminoglycosides, ß-lactam/ß-lactamase inhibitor combinations, tetracyclines, and fluoroquinolones. Phylogenetics analysis suggested potential shared routes of transmission among human, chicken, farm and slaughterhouse environments. Taken together, this study offers a cross-sectional epidemiological insight into the genomic epidemiology and antimicrobial landscape of S. Typhimurium isolates in Singapore, informing strategies for future public health and food safety surveillance.

Multimodal molecular landscape of response to Y90-resin microsphere radioembolization followed by nivolumab for advanced hepatocellular carcinoma.

BACKGROUND: Combination therapy with radioembolization (yttrium-90)-resin microspheres) followed by nivolumab has shown a promising response rate of 30.6% in a Phase II trial (CA209-678) for advanced hepatocellular carcinoma (HCC); however, the response mechanisms and relevant biomarkers remain unknown. METHODS: By collecting both pretreatment and on-treatment samples, we performed multimodal profiling of tissue and blood samples and investigated molecular changes associated with favorable responses in 33 patients from the trial. RESULTS: We found that higher tumor mutation burden, NCOR1 mutations and higher expression of interferon gamma pathways occurred more frequently in responders. Meanwhile, non-responders tended to be enriched for a novel Asian-specific transcriptomic subtype (Kaya_P2) with a high frequency of chromosome 16 deletions and upregulated cell cycle pathways. Strikingly, unlike other cancer types, we did not observe any association between T-cell populations and treatment response, but tumors from responders had a higher proportion of CXCL9+/CXCR3+ macrophages. Moreover, biomarkers discovered in previous immunotherapy trials were not predictive in the current cohort, suggesting a distinctive molecular landscape associated with differential responses to the combination therapy. CONCLUSIONS: This study unraveled extensive molecular changes underlying distinctive responses to the novel treatment and pinpointed new directions for harnessing combination therapy in patients with advanced HCC.

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.

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.

Radioembolisation with Y90-resin microspheres followed by nivolumab for advanced hepatocellular carcinoma (CA 209-678): a single arm, single centre, phase 2 trial.

BACKGROUND: Therapeutic synergism between radiotherapy and immune checkpoint blockade has been observed in preclinical models of hepatocellular carcinoma. We aimed to study the safety and efficacy of sequential radioembolisation with yttrium-90-resin microspheres (Y90-radioembolisation) followed by nivolumab in patients with advanced hepatocellular carcinoma. METHODS: Patients with Child-Pugh A cirrhosis and advanced hepatocellular carcinoma not suitable for curative surgery were treated with Y90-radioembolisation followed by intravenous nivolumab 240 mg 21 days after Y90-radioembolisation and every 2 weeks thereafter. The primary endpoint, assessed in the per-protocol population, was the objective response rate, determined by RECIST version 1.1, defined as the proportion of patients with a confirmed complete or partial response observed for lesions both within and outside the Y90-radioembolisation field. This study is registered with ClinicalTrials.gov, NCT03033446 and has been completed. FINDINGS: 40 patients were enrolled, of whom 36 received Y90-radioembolisation followed by nivolumab. One (3%) patient had a complete response and ten (28%) had a partial response; the objective response rate was 30·6% (95% CI 16·4-48·1). The most common treatment-related adverse events of any grade were pruritus (18 [50%] of 36 patients) and maculopapular rash (13 [36%]). Two (6%) patients experienced grade 3-4 treatment-related adverse events: one patient had a grade 3 increase in alanine aminotransferase levels, grade 3 bilirubin increase, and grade 4 increase in aspartate aminotransferase levels, while the other had a grade 3 maculopapular rash. Five (14%) patients had a treatment-related serious adverse event (Steven-Johnson syndrome, hepatitis E infection, fever, liver abscesses, and ascites). INTERPRETATION: Y90-radioembolisation followed by nivolumab resulted in an encouraging objective response rate in patients with advanced hepatocellular carcinoma, although the activity observed was not as high as the study was powered for. This strategy should be further evaluated in patients with Barcelona Clinic Liver Clinic (BCLC) stage B hepatocellular carcinoma that is ineligible or refractory to transarterial chemoembolisation and patients with BCLC C disease without extrahepatic spread. FUNDING: National Medical Research Council Singapore, Bristol-Myers Squibb, Sirtex.

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.

Genomic landscape of lung adenocarcinoma in East Asians.

Lung cancer is the world's leading cause of cancer death and shows strong ancestry disparities. By sequencing and assembling a large genomic and transcriptomic dataset of lung adenocarcinoma (LUAD) in individuals of East Asian ancestry (EAS; n = 305), we found that East Asian LUADs had more stable genomes characterized by fewer mutations and fewer copy number alterations than LUADs from individuals of European ancestry. This difference is much stronger in smokers as compared to nonsmokers. Transcriptomic clustering identified a new EAS-specific LUAD subgroup with a less complex genomic profile and upregulated immune-related genes, allowing the possibility of immunotherapy-based approaches. Integrative analysis across clinical and molecular features showed the importance of molecular phenotypes in patient prognostic stratification. EAS LUADs had better prediction accuracy than those of European ancestry, potentially due to their less complex genomic architecture. This study elucidated a comprehensive genomic landscape of EAS LUADs and highlighted important ancestry differences between the two cohorts.

The spatial organization of intra-tumour heterogeneity and evolutionary trajectories of metastases in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) has one of the poorest survival rates among cancers. Using multi-regional sampling of nine resected HCC with different aetiologies, here we construct phylogenetic relationships of these sectors, showing diverse levels of genetic sharing, spanning early to late diversification. Unlike the variegated pattern found in colorectal cancers, a large proportion of HCC display a clear isolation-by-distance pattern where spatially closer sectors are genetically more similar. Two resected intra-hepatic metastases showed genetic divergence occurring before and after primary tumour diversification, respectively. Metastatic tumours had much higher variability than their primary tumours, suggesting that intra-hepatic metastasis is accompanied by rapid diversification at the distant location. The presence of co-existing mutations offers the possibility of drug repositioning for HCC treatment. Taken together, these insights into intra-tumour heterogeneity allow for a comprehensive understanding of the evolutionary trajectories of HCC and suggest novel avenues for personalized therapy.