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

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.

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.

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.

A new perspective on the immune escape mechanism in HCC: onco-foetal reprogramming.

We found a shared immunosuppressive microenvironment between foetal liver and hepatocellular carcinoma (HCC) which includes the re-emergence of foetal-associated endothelial cells (PLVAP/VEGFR2) and foetal-like (FOLR2) tumour-associated macrophages in HCC, mediated via VEGF-NOTCH signalling. The discoveries suggest possible novel targets for therapeutic interventions in HCC.

Pharmacogenetic effects of regulatory nuclear receptors (PXR, CAR, RXRα and HNF4α) on docetaxel disposition in Chinese nasopharyngeal cancer patients.

PURPOSE: This exploratory study was aimed at elucidating the pharmacogenetics of regulatory nuclear receptors (PXR, CAR, RXRα and HNF4α) and their implications on docetaxel pharmacokinetics and pharmacodynamics in local Chinese nasopharyngeal cancer patients. METHODS: A total of 59 single nucleotide polymorphisms (SNPs), including tag-SNPs and functionally relevant SNPs of the genes encoding these regulatory nuclear receptors (PXR/NR1I2, CAR/NR1I3, RXRα/NR2B1 and HNF4α/NR2A1), were profiled in the patients enrolled in our study by direct sequencing (N = 50). The generalized linear model was employed to estimate the haplotypic effects on the pharmacokinetics and pharmacodynamics of the patients. RESULTS: The pharmacokinetic profiles of docetaxel in these patients were characterized by marked interindividual variability, with approximately four- to sixfold variations observed in Cmax, AUC0-∞ and CL. Individual SNP association tests revealed that polymorphisms in NR2B1 and NR2A1 were significantly correlated with altered docetaxel pharmacokinetics. Subsequent haplotype association analysis identified the NR2B1 LD block 2 AG haplotype [*+4458G>A(rs3132291) and *+4988A>G(rs4842198)] to be significantly associated with altered pharmacokinetics, in which patients carrying two copies of the AG haplotype had approximately a 20 % decreased Cmax and AUC0-∞ and a 21 % increased CL compared to those who carried only one copy or no copies of the haplotype. A number of SNPs in NR1I2, NR1I3, NR2B1 and NR2A1 were also associated with a significant decrease in blood counts from baseline. No haplotype was found to exert any effects on the pharmacodynamics parameters. CONCLUSIONS: The present exploratory study identified several SNPs in the genes encoding regulatory nuclear receptors which may account for the interpatient variability in docetaxel pharmacokinetics and pharmacodynamics. These findings highlight the important role of regulatory nuclear receptors on the disposition of docetaxel.

Genetic variations of NR1I3 and NR2B1 in Asian populations.

Several nuclear receptors are being increasingly recognized for their role as master xenosensors. Among them, CAR-RXRα heterodimer, as encoded by NR1I3 and NR2B1, responds to the presence of drug compounds and regulates the transcription of a wide array of genes involved in their disposition. To investigate the frequency distribution and linkage disequilibrium patterns of NR1I3 and NR2B1 genetic variations, these genes were screened in 168 healthy local Asian subjects, namely Chinese, Malays, and Indians (n=56 subjects each). A total of 38 and 88 SNPs were identified in NR1I3 and NR2B1, respectively. Among them, there were 13 and 43 novel SNPs present at low allelic frequencies (<10%) in NR1I3 and NR2B1, respectively. Notably, the genetic variations in the NR1I3 and NR2B1 genes were mainly confined to the introns whilst the exons were highly conserved across the ethnic populations. Indians harboured distinct frequency distributions from Chinese and Malays in both genes. Based on the linkage disequilibrium patterns of both genes, a number of tag-SNPs were selected for each population (n=8-13 for NR1I3; n=12-18 for NR2B1). In-silico prediction analyses revealed a number of possible functional SNPs. Our data would be valuable for future pharmacogenetic studies on the drug substrates of CAR-RXRα target genes.

Genetic variations of NR2A1 in Asian populations: implications in pharmacogenetics studies.

HNF4α (encoded by gene NR2A1) is a dominant transcriptional regulator of various drug disposition genes. It forms a circuitry of molecular cross-talk with other nuclear receptors such as PXR and CAR to synergistically initiate transcription. This study reports on the frequency, linkage disequilibrium pattern and tag-SNP selection of NR2A1 polymorphisms in three local Asian populations, namely Chinese, Malays and Indians (n = 56 subjects each). A total of 69 polymorphisms were identified in the genomic region of NR2A1, of which thirty-three were novel polymorphisms with low allelic frequencies (<0.02). The exonic region of NR2A1 was highly conserved with only 4 novel and 1 reported SNPs identified at low allelic frequencies of less than 0.02. Based on the criteria of MAF ≥ 0.05 and R(2) ≥ 0.80, there were 19, 20 and 22 tag-SNPs selected to represent the genetic polymorphisms of NR2A1 in Chinese, Malays and Indians, respectively. In-silico predictions suggested that some of these polymorphic variants may exert functional effects through affecting the binding sites of transcription and splicing factors. Our study provides valuable information on the genetic variability of NR2A1 which would be useful for pharmacogenetics studies in the local Asian populations.

Influence of SLCO1B3 haplotype-tag SNPs on docetaxel disposition in Chinese nasopharyngeal cancer patients.

WHAT IS ALREADY KNOWN ABOUT THIS SUBJECT: SLCO1B3 is an influx transporter located at the hepatocyte basolateral membrane and it is involved in the uptake of a broad range of drug substrates including docetaxel. The pharmacogenetics of SLCO1B3 is not well characterized and previous in vivo and in vitro studies reported conflicting results with regards to the functional effects of the limited number of SLCO1B3 polymorphisms that were studied. Docetaxel displays a wide interindividual variability in its pharmacokinetics and pharmacodynamics and an understanding of SLCO1B3 pharmacogenetics might provide clinical benefits in guiding docetaxel dosing. WHAT THIS STUDY ADDS: The SLCO1B3 gene was comprehensively screened in the local healthy Asian populations (n= 168). A strong linkage disequilibrium pattern was detected across a total of 88 polymorphisms and 15 haplotype-tag SNPs (htSNPs) were identified. These htSNPs were profiled in a cohort of Chinese nasopharyngeal cancer (NPC) patients (n= 50). Genotypic-phenotypic analysis showed that a haplotypic construct comprising of four variants [IVS4+76G>A, 699G>A(Met233Ile), IVS12-5676A>G, and *347_*348insA] was the critical determinant of docetaxel disposition. This study suggests that the comprehensive screening and haplotypic linkage analysis of SLCO1B3 can better elucidate its pharmacogenetic effects on interpatient variability of docetaxel and other putative drug substrates. Further studies are warranted in cancer patients belonging to other ethnic groups. AIMS To completely screen the SLCO1B3 gene in three distinct healthy Asian populations (Chinese, Malay and Indian, n= 168) and investigate the influence of haplotype-tag SNPs (htSNPs) on docetaxel disposition in 50 nasopharyngeal carcinoma patients. METHODS: Genomic DNA of individuals was screened for SLCO1B3 polymorphisms by direct sequencing. htSNPs were derived based on the sequence clustering algorithm and profiled in the patients. Population based genetic association analysis was performed using Haplostats package implemented in R and PLINK. RESULTS: A strong linkage disequilibrium pattern was detected across a total of 88 polymorphisms and 15-htSNPs were identified. The SLCO1B3 haplotypic region comprising seven htSNPs was found to be significantly associated with docetaxel clearance (P= 0.003). Conditional haplotype analyses revealed that the haplotypic constructs comprising the IVS4+76G>A, 699G>A(Met233Ile), IVS12-5676A>G, and *347_*348insA polymorphisms were critical determinants of variability in docetaxel disposition [clearance and area under the plasma concentration-time curve (AUC(0,∞)): r(2) = 29% and 22%, respectively]. Patients harbouring the GAG*347insA haplotype were significantly associated with a 30% decrease in clearance and a 40% increase in AUC(0,∞) of docetaxel compared with patients harbouring the reference haplotype, GGA*347wt (P= 0.025 and 0.018, respectively). In contrast, a 50% higher clearance was observed in patients carrying the GAG*347wt haplotype compared with those with the reference haplotype (P= 0.002). The functional SLCO1B3 haplotypic constructs included the widely studied Met233Ile variant and *347_*348insA located in the putative miR-890 binding site in the 3'-untranslated region which may influence the transport characteristics of SLCO1B3. CONCLUSIONS: This study highlights the importance of SLCO1B3 polymorphic variations in influencing docetaxel disposition in nasopharyngeal carcinoma patients.

The effects of CYP3A4, CYP3A5, ABCB1, ABCC2, ABCG2 and SLCO1B3 single nucleotide polymorphisms on the pharmacokinetics and pharmacodynamics of docetaxel in nasopharyngeal carcinoma patients.

PURPOSE: This exploratory study aimed to explain the interindividual variabilities of docetaxel pharmacokinetics and pharmacodynamics in Asian nasopharyngeal carcinoma patients (n = 54) through the genotyping of CYP3A4, CYP3A5, ABCB1, ABCC2, ABCG2 and SLCO1B3 genes. METHODS: Docetaxel was administered over 1 h on days 1, 8, and 15 every 28 days at 30 mg/m(2)/dose. Genomic DNA was isolated from peripheral blood and genotyped for the selected polymorphisms in the candidate genes. Docetaxel pharmacokinetic parameters were estimated by non-compartmental modelling. RESULTS: Patients homozygous for the variant allele (GG) of SLCO1B3 rs11045585 (IVS12-5676A > G) had significantly higher area under the plasma concentration-time curve of docetaxel (P = 0.026) and lower clearance (P = 0.036) compared to patients with AA/AG genotypes. Patients harbouring the heterozygous genotype (GA + GT + TA) for ABCB1 rs2032582 (2677G > T/A) had the highest percentage decrease in nadir haemoglobin from cycle 1 baseline compared to those with GG/TT genotypes (P = 0.006). Similar trend was observed for ABCB1 rs1045642 (3435C > T) with heterozygotes (CT) having the highest percentage decrease in nadir haemoglobin from cycle 1 baseline compared to those with CC/TT genotypes (P = 0.066). CONCLUSIONS: This study suggests that the cooperative influence of functional polymorphisms in SLCO1B3 and ABCB1 genes may be responsible for the interindividual variability in docetaxel disposition in Asian nasopharyngeal cancer patients.