Single-cell multiomics reveals the interplay of clonal evolution and cellular plasticity in hepatoblastoma.

Hepatoblastomas (HB) display heterogeneous cellular phenotypes that influence the clinical outcome, but the underlying mechanisms are poorly understood. Here, we use a single-cell multiomic strategy to unravel the molecular determinants of this plasticity. We identify a continuum of HB cell states between hepatocytic (scH), liver progenitor (scLP) and mesenchymal (scM) differentiation poles, with an intermediate scH/LP population bordering scLP and scH areas in spatial transcriptomics. Chromatin accessibility landscapes reveal the gene regulatory networks of each differentiation pole, and the sequence of transcription factor activations underlying cell state transitions. Single-cell mapping of somatic alterations reveals the clonal architecture of each tumor, showing that each genetic subclone displays its own range of cellular plasticity across differentiation states. The most scLP subclones, overexpressing stem cell and DNA repair genes, proliferate faster after neo-adjuvant chemotherapy. These results highlight how the interplay of clonal evolution and epigenetic plasticity shapes the potential of HB subclones to respond to chemotherapy.

DLK1/DIO3 locus upregulation by a ß-catenin-dependent enhancer drives cell proliferation and liver tumorigenesis.

The CTNNB1 gene, encoding ß-catenin, is frequently mutated in hepatocellular carcinoma (HCC, ∼30%) and in hepatoblastoma (HB, >80%), in which DLK1/DIO3 locus induction is correlated with CTNNB1 mutations. Here, we aim to decipher how sustained ß-catenin activation regulates DLK1/DIO3 locus expression and the role this locus plays in HB and HCC development in mouse models deleted for Apc (ApcΔhep) or Ctnnb1-exon 3 (ß-cateninΔExon3) and in human CTNNB1-mutated hepatic cancer cells. We identified an enhancer site bound by TCF-4/ß-catenin complexes in an open conformation upon sustained ß-catenin activation (DLK1-Wnt responsive element [WRE]) and increasing DLK1/DIO3 locus transcription in ß-catenin-mutated human HB and mouse models. DLK1-WRE editing by CRISPR-Cas9 approach impaired DLK1/DIO3 locus expression and slowed tumor growth in subcutaneous CTNNB1-mutated tumor cell grafts, ApcΔhep HB and ß-cateninΔExon3 HCC. Tumor growth inhibition resulted either from increased FADD expression and subsequent caspase-3 cleavage in the first case or from decreased expression of cell cycle actors regulated by FoxM1 in the others. Therefore, the DLK1/DIO3 locus is an essential determinant of FoxM1-dependent cell proliferation during ß-catenin-driven liver tumorigenesis. Targeting the DLK1-WRE enhancer to silence the DLK1/DIO3 locus might thus represent an interesting therapeutic strategy to restrict tumor growth in primary liver cancers with CTNNB1 mutations.

Preneoplastic liver colonization by 11p15.5 altered mosaic cells in young children with hepatoblastoma.

Pediatric liver tumors are very rare tumors with the most common diagnosis being hepatoblastoma. While hepatoblastomas are predominantly sporadic, around 15% of cases develop as part of predisposition syndromes such as Beckwith-Wiedemann (11p15.5 locus altered). Here, we identify mosaic genetic alterations of 11p15.5 locus in the liver of hepatoblastoma patients without a clinical diagnosis of Beckwith-Wiedemann syndrome. We do not retrieve these alterations in children with other types of pediatric liver tumors. We show that mosaic 11p15.5 alterations in liver FFPE sections of hepatoblastoma patients display IGF2 overexpression and H19 downregulation together with an alteration of the liver zonation. Moreover, mosaic livers' microenvironment is enriched in extracellular matrix and angiogenesis. Spatial transcriptomics and single-nucleus RNAseq analyses identify a 60-gene signature in 11p15.5 altered hepatocytes. These data provide insights for 11p15.5 mosaicism detection and its functional consequences during the early steps of carcinogenesis.

Engineering of ultraID, a compact and hyperactive enzyme for proximity-dependent biotinylation in living cells.

Proximity-dependent biotinylation (PDB) combined with mass spectrometry analysis has established itself as a key technology to study protein-protein interactions in living cells. A widespread approach, BioID, uses an abortive variant of the E. coli BirA biotin protein ligase, a quite bulky enzyme with slow labeling kinetics. To improve PDB versatility and speed, various enzymes have been developed by different approaches. Here we present a small-size engineered enzyme: ultraID. We show its practical use to probe the interactome of Argonaute-2 after a 10 min labeling pulse and expression at physiological levels. Moreover, using ultraID, we provide a membrane-associated interactome of coatomer, the coat protein complex of COPI vesicles. To date, ultraID is the smallest and most efficient biotin ligase available for PDB and offers the possibility of investigating interactomes at a high temporal resolution.

Integrated Genomic Analysis Identifies Driver Genes and Cisplatin-Resistant Progenitor Phenotype in Pediatric Liver Cancer.

Pediatric liver cancers (PLC) comprise diverse diseases affecting infants, children, and adolescents. Despite overall good prognosis, PLCs display heterogeneous response to chemotherapy. Integrated genomic analysis of 126 pediatric liver tumors showed a continuum of driver mechanisms associated with patient age, including new targetable oncogenes. In 10% of patients with hepatoblastoma, all before three years old, we identified a mosaic premalignant clonal expansion of cells altered at the 11p15.5 locus. Analysis of spatial and longitudinal heterogeneity revealed an important plasticity between "hepatocytic," "liver progenitor," and "mesenchymal" molecular subgroups of hepatoblastoma. We showed that during chemotherapy, "liver progenitor" cells accumulated massive loads of cisplatin-induced mutations with a specific mutational signature, leading to the development of heavily mutated relapses and metastases. Drug screening in PLC cell lines identified promising targets for cisplatin-resistant progenitor cells, validated in mouse xenograft experiments. These data provide new insights into cisplatin resistance mechanisms in PLC and suggest alternative therapies. SIGNIFICANCE: PLCs are deadly when they resist chemotherapy, with limited alternative treatment options. Using a multiomics approach, we identified PLC driver genes and the cellular phenotype at the origin of cisplatin resistance. We validated new treatments targeting these molecular features in cell lines and xenografts.This article is highlighted in the In This Issue feature, p. 2355.

DNA Methylation Signatures Reveal the Diversity of Processes Remodeling Hepatocellular Carcinoma Methylomes.

BACKGROUND AND AIMS: DNA methylation patterns are highly rearranged in HCCs. However, diverse sources of variation are intermingled in cancer methylomes, precluding the precise characterization of underlying molecular mechanisms. We developed a computational framework (methylation signature analysis with independent component analysis [MethICA]) leveraging independent component analysis to disentangle the diverse processes contributing to DNA methylation changes in tumors. APPROACH AND RESULTS: Applied to a collection of 738 HCCs, MethICA unraveled 13 stable methylation components preferentially active in specific chromatin states, sequence contexts, and replication timings. These included signatures of general processes associated with sex and age but also signatures related to specific driver events and molecular subgroups. Catenin beta 1 mutations were major modulators of methylation patterns in HCC, characterized by a targeted hypomethylation of transcription factor 7-bound enhancers in the vicinity of Wnt target genes as well as a widespread hypomethylation of late-replicated partially methylated domains. By contrast, demethylation of early replicated highly methylated domains was a signature of replication stress, leading to an extensive hypomethylator phenotype in cyclin-activated HCC. Inactivating mutations of the chromatin remodeler AT-rich interactive domain-containing protein 1A were associated with epigenetic silencing of differentiation-promoting transcriptional networks, also detectable in cirrhotic liver. Finally, a hypermethylation signature targeting polycomb-repressed chromatin domains was identified in the G1 molecular subgroup with progenitor features. CONCLUSIONS: This study elucidates the diversity of processes remodeling HCC methylomes and reveals the epigenetic and transcriptional impact of driver alterations.