SUMMARY Hepatoblastoma, the most prevalent pediatric liver cancer, almost always carries a WNT-activating CTNNB1 mutation, yet exhibits notable molecular heterogeneity. To characterize this heterogeneity and identify novel targeted therapies, we performed comprehensive analysis of hepatoblastomas and tumor-derived organoids using single-cell RNA-seq, spatial transcriptomics, single-cell ATAC-seq and high throughput drug profiling. We identified two distinct tumor epithelial signatures: hepatic ‘fetal-like’ and WNT-high ‘embryonal-like’ signatures, displaying divergent WNT signaling patterns. The liver-specific WNT targets were enriched in the fetal-like group, while the embryonal-like group was enriched in canonical WNT target genes. Gene regulatory network analysis revealed enrichment of regulons related to hepatic function such as bile acid, lipid and xenobiotic metabolism in the fetal-like subgroup but not in the embryonal-like subgroup. In addition, the dichotomous expression pattern of the transcription factors HNF4A and LEF1 allowed for a clear distinction between the fetal- and embryonal-like tumors. We also performed high-throughput drug screening using patient-derived tumor organoids and identified sensitivity to multiple inhibitor classes, most notably HDAC inhibitors. Intriguingly, embryonal-like tumor organoids, but not fetal-like tumor organoids, were sensitive to FGFR inhibitor treatments, suggesting a dependency on FGFR signaling. In summary, our data uncover the molecular and drug sensitivity landscapes of hepatoblastoma and pave the way for the development of targeted therapies.

Multi-dimensional profiling of hepatoblastomas and patient-derived tumor organoids uncovers tumor subpopulations with divergent WNT activation profiles and identifies pan-hepatoblastoma drug sensitivities