Integrated genomic analysis identifies driver genes and cisplatin-resistant progenitor phenotype in pediatric liver cancer

Indolent cancer and pattern of progression: Two missing parameters in trial design for hepatology

The indolent and aggressive behaviors of HCC might have a role in clinical trial (CT) results; however, the indolent HCC is less analyzed compared to others cancer. Indolent profile could be characterized as follows: (1) patients with low risk of progression itself due to the HCC molecular profile and/or due to the interaction between cancer cell their microenvironment; (2) patients who achieve objective response or present spontaneous regression; and (3) patients who develop radiological progression with no consequence on either the liver function or general status, and without trigger a change in the tumor stage. Patients with "indolent HCC" generally never develop cancer-related symptoms neither die for HCC-related causes. Thus, we hypothesize that the imbalance in the proportion of "indolent" versus "aggressive HCC" between arms or the underestimation/overestimation of HCC behavior at baseline in single-arm CT could be associated with CT failure or under-overestimation of trial results. The "indolent progression" may also explain the discrepancy between radiological progression-based end points and survival. Moreover, we discuss the related causes that explain the indolent profile of HCC and propose (1) refining the progression-related end point by the pattern of progression to minimize the limitations of the current end points; (2) considering alternative statistical tools for survival analysis such as milestone survival, or restricted mean survival time to capture the value of indolent HCC. According to these considerations, we propose incorporating novel end points into the single arm of phase I/II CT as exploratory analysis or as a secondary end point in phase III CT.

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.

Molecular Profiling of a Hepatocellular Neoplasm Not Otherwise Specified (HCN-NOS) Demonstrates Distinct Molecular Features in Hepatoblastoma and HCC-Like Components

Hepatoblastomas (HB) are embryonal tumors with quiet genomes diagnosed mostly in children under 3 years of age and often cured by surgical resection and chemotherapy. However, a subset of HBs behave aggressively, displaying characteristic histologic features and higher genomic instability. Hepatocellular neoplasm-not otherwise specified (HCN-NOS) is a provisional diagnostic category for tumors exhibiting either intermediate or a combination of both HB and hepatocellular carcinoma (HCC) histological features. In this study, we characterized an HCN-NOS diagnosed in a 3-year-old patient presenting with a liver mass, in which both HB and HCC histological components were amendable to macro-dissection and molecular profiling. The spectrum of mutations, copy number changes, mRNA, and protein expression profiles within these 2 histologically distinct tumor areas demonstrate molecular heterogeneity and suggest intratumoral clonal evolution of this hepatocellular CTNNB1-mutant lesion.

A TGF-β-dominant chemoresistant phenotype of hepatoblastoma associated with aflatoxin exposure in children

BACKGROUND AND AIMS

Hepatoblastoma (HB) is the most common liver cancer in children, posing a serious threat to children's health. Chemoresistance is the leading cause of mortality in patients with HB. A more explicit definition of the features of chemotherapy resistance in HB represents a fundamental urgent need.

APPROACH AND RESULTS

We performed an integrative analysis including single-cell RNA sequencing, whole-exome sequencing, and bulk RNA sequencing in 180 HB samples, to reveal genomic features, transcriptomic profiles, and the immune microenvironment of HB. Multicolor immunohistochemistry staining and in vitro experiments were performed for validation. Here, we reported four HB transcriptional subtypes primarily defined by differential expression of transcription factors. Among them, the S2A subtype, characterized by strong expression of progenitor ( MYCN , MIXL1 ) and mesenchymal transcription factors ( TWIST1 , TBX5 ), was defined as a new chemoresistant subtype. The S2A subtype showed increased TGF-β cancer-associated fibroblast and an immunosuppressive microenvironment induced by the upregulated TGF-β of HB. Interestingly, the S2A subtype enriched SBS24 signature and significantly higher serum aflatoxin B1-albumin (AFB1-ALB) level in comparison with other subtypes. Functional assays indicated that aflatoxin promotes HB to upregulate TGF-β. Furthermore, clinical prognostic analysis showed that serum AFB1-ALB is a potential indicator of HB chemoresistance and prognosis.

CONCLUSIONS

Our studies offer new insights into the relationship between aflatoxin and HB chemoresistance and provide important implications for its diagnosis and treatment.

Mutational signature, cancer driver genes mutations and transcriptomic subgroups predict hepatoblastoma survival

BACKGROUND

Hepatoblastoma is the most frequent pediatric liver cancer. The current treatments lead to 80% of survival rate at 5 years. In this study, we evaluated the clinical relevance of molecular features to identify patients at risk of chemoresistance, relapse and death of disease.

METHODS

All the clinical data of 86 children with hepatoblastoma were retrospectively collected. Pathological slides were reviewed, tumor DNA sequencing (by whole exome, whole genome or target) and transcriptomic profiling with RNAseq or 300-genes panel were performed. Associations between the clinical, pathological, mutational and transcriptomic data were investigated.

RESULTS

High-risk patients represented 44% of our series and the median age at diagnosis was 21.9 months (range: 0-208). Alterations of the WNT/ß-catenin pathway and of the 11p15.5 imprinted locus were identified in 98% and 74% of the tumors, respectively. Other cancer driver genes mutations were only found in less than 11% of tumors. After neoadjuvant chemotherapy, disease-specific survival and poor response to neoadjuvant chemotherapy were associated with 'Liver Progenitor' (p = 0.00049, p < 0.0001) and 'Immune Cold' (p = 0.0011, p < 0.0001) transcriptomic tumor subtypes, SBS35 cisplatin mutational signature (p = 0.018, p = 0.001), mutations in rare cancer driver genes (p = 0.0039, p = 0.0017) and embryonal predominant histological type (p = 0.0013, p = 0.0077), respectively. Integration of the clinical and molecular features revealed a cluster of molecular markers associated with resistance to chemotherapy and survival, enlightening transcriptomic 'Immune Cold' and Liver Progenitor' as a predictor of survival independent of the clinical features.

CONCLUSIONS

Response to neoadjuvant chemotherapy and survival in children treated for hepatoblastoma are associated with genomic and pathological features independently of the clinical features.

Summary of biological research on hepatoblastoma: a scoping review

Background

Hepatoblastoma is the most prevalent primary hepatic malignancy in children, comprising 80% of pediatric hepatic malignancies and 1% of all pediatric malignancies. However, traditional treatments have proven inadequate in effectively curing hepatoblastoma, leading to a poor prognosis.

Methods

A literature search was conducted on multiple electronic databases (PubMed and Google Scholar). A total of 86 articles were eligible for inclusion in this review.

Result

This review aims to consolidate recent developments in hepatoblastoma research, focusing on the latest advances in cancer-associated genomics, epigenetic studies, transcriptional programs and molecular subtypes. We also discuss the current treatment approaches and forthcoming strategies to address cancer-associated biological challenges.

Conclusion

To provide a comprehensive summary of the molecular mechanisms associated with hepatoblastoma occurrence, this review highlights three key aspects: genomics, epigenetics, and transcriptomics. Our review aims to facilitate the exploration of novel molecular mechanisms and the development of innovative clinical treatment strategies for hepatoblastoma.

Targeting G9a/DNMT1 methyltransferase activity impedes IGF2-mediated survival in hepatoblastoma

BACKGROUND

As the variable clinical outcome of patients with hepatoblastoma (HB) cannot be explained by genetics alone, the identification of drugs with the potential to effectively reverse epigenetic alterations is a promising approach to overcome poor therapy response. The gene ubiquitin like with PHD and ring finger domains 1 (UHRF1) represents an encouraging epigenetic target due to its regulatory function in both DNA methylation and histone modifications and its clinical relevance in HB.

METHODS

Patient-derived xenograft in vitro and in vivo models were used to study drug response. The mechanistic basis of CM-272 treatment was elucidated using RNA sequencing and western blot experiments.

RESULTS

We validated in comprehensive data sets that UHRF1 is highly expressed in HB and associated with poor outcomes. The simultaneous pharmacological targeting of UHRF1-dependent DNA methylation and histone H3 methylation by the dual inhibitor CM-272 identified a selective impact on HB patient-derived xenograft cell viability while leaving healthy fibroblasts unaffected. RNA sequencing revealed downregulation of the IGF2-activated survival pathway as the main mode of action of CM-272 treatment, subsequently leading to loss of proliferation, hindered colony formation capability, reduced spheroid growth, decreased migration potential, and ultimately, induction of apoptosis in HB cells. Importantly, drug response depended on the level of IGF2 expression, and combination assays showed a strong synergistic effect of CM-272 with cisplatin. Preclinical testing of CM-272 in a transplanted patient-derived xenograft model proved its efficacy but also uncovered side effects presumably caused by its strong antitumor effect in IGF2-driven tumors.

CONCLUSIONS

The inhibition of UHRF1-associated epigenetic traces, such as IGF2-mediated survival, is an attractive approach to treat high-risk HB, especially when combined with the standard-of-care therapeutic cisplatin.

Computed tomography imaging phenotypes of hepatoblastoma identified from radiomics signatures are associated with the efficacy of neoadjuvant chemotherapy

BACKGROUND

Though neoadjuvant chemotherapy has been widely used in the treatment of hepatoblastoma, there still lacks an effective way to predict its effect.

OBJECTIVE

To characterize hepatoblastoma based on radiomics image features and identify radiomics-based lesion phenotypes by unsupervised machine learning, intended to build a classifier to predict the response to neoadjuvant chemotherapy.

MATERIALS AND METHODS

In this retrospective study, we segmented the arterial phase images of 137 cases of pediatric hepatoblastoma and extracted the radiomics features using PyRadiomics. Then unsupervised k-means clustering was applied to cluster the tumors, whose result was verified by t-distributed stochastic neighbor embedding (t-SNE). The least absolute shrinkage and selection operator (LASSO) regression was used for feature selection, and the clusters were visually analyzed by radiologists. The correlations between the clusters, clinical and pathological parameters, and qualitative radiological features were analyzed.

RESULTS

Hepatoblastoma was clustered into three phenotypes (homogenous type, heterogenous type, and nodulated type) based on radiomics features. The clustering results had a high correlation with response to neoadjuvant chemotherapy (P=0.02). The epithelial ratio and cystic components in radiological features were also associated with the clusters (P=0.029 and 0.008, respectively).

CONCLUSIONS

This radiomics-based cluster system may have the potential to facilitate the precise treatment of hepatoblastoma. In addition, this study further demonstrated the feasibility of using unsupervised machine learning in a disease without a proper imaging classification system.

Alpha-Momorcharin, a type I ribosome inactivating protein, induced apoptosis of hepatocellular carcinoma SK-HEP-1 cells through mitochondrial pathway

Alpha-Momorcharin (α-MMC), as one of the most important type I RIPs, has been reported to exert inhibitory effects against various tumour cells through its N-glycosidase activity. The present study was designed to propose an efficient purification strategy and explored its mechanism of apoptosis signalling pathway against human liver cancer cells SK-Hep-1. α-MMC can be successfully obtained by our purification strategy combining ion-exchange and gel-filtration chromatography. The functional studies revealed that α-MMC obviously increased the level of ROS and apoptosis rate, induced cell cycle arrest in the G1 phase, and depolarised MMP of SK-Hep-1 cells. To further confirm whether α-MMC could induce mitochondria involved apoptosis, we found that PARP-1, Caspase-3, Caspase-9, and BCL-2 were downregulated upon α-MMC. Taken together, these results suggested that this natural purified α-MMC can induce apoptosis involved mitochondria and may serve as a potential novel therapeutic drug in the treatment of human liver cancer in the future.

Genomic profiling of subcutaneous patient-derived xenografts reveals immune constraints on tumor evolution in childhood solid cancer

Subcutaneous patient-derived xenografts (PDXs) are an important tool for childhood cancer research. Here, we describe a resource of 68 early passage PDXs established from 65 pediatric solid tumor patients. Through genomic profiling of paired PDXs and patient tumors (PTs), we observe low mutational similarity in about 30% of the PT/PDX pairs. Clonal analysis in these pairs show an aggressive PT minor subclone seeds the major clone in the PDX. We show evidence that this subclone is more immunogenic and is likely suppressed by immune responses in the PT. These results suggest interplay between intratumoral heterogeneity and antitumor immunity may underlie the genetic disparity between PTs and PDXs. We further show that PDXs generally recapitulate PTs in copy number and transcriptomic profiles. Finally, we report a gene fusion LRPAP1-PDGFRA. In summary, we report a childhood cancer PDX resource and our study highlights the role of immune constraints on tumor evolution.

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.

Identification and experimental validation of druggable epigenetic targets in hepatoblastoma

BACKGROUND & AIMS

Hepatoblastoma (HB) is the most frequent childhood liver cancer. Patients with aggressive tumors have limited therapeutic options; therefore, a better understanding of HB pathogenesis is needed to improve treatment. HBs have a very low mutational burden; however, epigenetic alterations are increasingly recognized. We aimed to identify epigenetic regulators consistently dysregulated in HB and to evaluate the therapeutic efficacy of their targeting in clinically relevant models.

METHODS

We performed a comprehensive transcriptomic analysis of 180 epigenetic genes. Data from fetal, pediatric, adult, peritumoral (n = 72) and tumoral (n = 91) tissues were integrated. Selected epigenetic drugs were tested in HB cells. The most relevant epigenetic target identified was validated in primary HB cells, HB organoids, a patient-derived xenograft model, and a genetic mouse model. Transcriptomic, proteomic and metabolomic mechanistic analyses were performed.

RESULTS

Altered expression of genes regulating DNA methylation and histone modifications was consistently observed in association with molecular and clinical features of poor prognosis. The histone methyltransferase G9a was markedly upregulated in tumors with epigenetic and transcriptomic traits of increased malignancy. Pharmacological targeting of G9a significantly inhibited growth of HB cells, organoids and patient-derived xenografts. Development of HB induced by oncogenic forms of β-catenin and YAP1 was ablated in mice with hepatocyte-specific deletion of G9a. We observed that HBs undergo significant transcriptional rewiring in genes involved in amino acid metabolism and ribosomal biogenesis. G9a inhibition counteracted these pro-tumorigenic adaptations. Mechanistically, G9a targeting potently repressed the expression of c-MYC and ATF4, master regulators of HB metabolic reprogramming.

CONCLUSIONS

HBs display a profound dysregulation of the epigenetic machinery. Pharmacological targeting of key epigenetic effectors exposes metabolic vulnerabilities that can be leveraged to improve the treatment of these patients.

IMPACT AND IMPLICATIONS

In spite of recent advances in the management of hepatoblastoma (HB), treatment resistance and drug toxicity are still major concerns. This systematic study reveals the remarkable dysregulation in the expression of epigenetic genes in HB tissues. Through pharmacological and genetic experimental approaches, we demonstrate that the histone-lysine-methyltransferase G9a is an excellent drug target in HB, which can also be harnessed to enhance the efficacy of chemotherapy. Furthermore, our study highlights the profound pro-tumorigenic metabolic rewiring of HB cells orchestrated by G9a in coordination with the c-MYC oncogene. From a broader perspective, our findings suggest that anti-G9a therapies may also be effective in other c-MYC-dependent tumors.

RSK2 inactivation cooperates with AXIN1 inactivation or β-catenin activation to promote hepatocarcinogenesis

BACKGROUND & AIMS

Recurrent somatic mutations of the RPS6KA3 gene encoding for the serine/threonine kinase RSK2 were identified in hepatocellular carcinomas (HCCs), suggesting its tumour-suppressive function. Our goal was to demonstrate the tumour suppressor role of RSK2 in the liver and investigate the functional consequences of its inactivation.

METHODS

We analysed a series of 1,151 human HCCs for RSK2 mutations and 20 other driver genetic alterations. We then modelled RSK2 inactivation in mice in various mutational contexts recapitulating or not those naturally found in human HCC, using transgenic mice and liver-specific carcinogens. These models were monitored for liver tumour appearance and subjected to phenotypic and transcriptomic analyses. Functional consequences of RSK2 rescue were also investigated in a human RSK2-deficient HCC cell line.

RESULTS

RSK2-inactivating mutations are specific to human HCC and frequently co-occur with AXIN1-inactivating or β-catenin-activating mutations. Modelling of these co-occurrences in mice showed a cooperative effect in promoting liver tumours with transcriptomic profiles recapitulating those of human HCCs. By contrast, there was no cooperation in liver tumour induction between RSK2 loss and BRAF-activating mutations chemically induced by diethylnitrosamine. In human liver cancer cells, we also showed that RSK2 inactivation confers some dependency to the activation of RAS/MAPK signalling that can be targeted by MEK inhibitors.

CONCLUSIONS

Our study demonstrates the tumour suppressor role of RSK2 and its specific synergistic effect in hepatocarcinogenesis when its loss of function is specifically combined with AXIN1 inactivation or β-catenin activation. Furthermore, we identified the RAS/MAPK pathway as a potential therapeutic target for RSK2-inactivated liver tumours.

IMPACT AND IMPLICATIONS

This study demonstrated the tumour suppressor role of RSK2 in the liver and showed that its inactivation specifically synergises with AXIN1 inactivation or β-catenin activation to promote the development of HCC with similar transcriptomic profiles as found in humans. Furthermore, this study highlights that activation of the RAS/MAPK pathway is one of the key signalling pathways mediating the oncogenic effect of RSK2 inactivation that can be targeted with already available anti-MEK therapies.

Pediatric Tumors and Developmental Anomalies: A French Nationwide Cohort Study

OBJECTIVE

To assess the associations between congenital abnormalities and pediatric malignancies and evaluate the potential underlying molecular basis by collecting information on pediatric patients with cancer and congenital abnormalities.

STUDY DESIGN

Tumeur Et Développement is a national, prospective, and retrospective multicenter study recording data of children with cancer and congenital abnormalities. When feasible, blood and tumoral samples are collected for virtual biobanking.

RESULTS

From June 2013 to December 2019, 679 associations between pediatric cancers and congenital abnormalities were recorded. The most represented cancers were central nervous system tumors (n = 139; 20%), leukemia and myelodysplastic syndromes (n = 123; 18.1%), and renal tumors (n = 101; 15%). Congenital abnormalities were not related to any known genetic disorder in 66.5% of cases. In this group, the most common anomaly was intellectual disability (22.3%), followed by musculoskeletal (14.2%) and genitourinary anomalies (12.4%). Intellectual disability was mostly associated with hematologic malignancies. Embryonic tumors (neuroblastoma, Wilms tumor, and rhabdomyosarcoma) were associated with consistent abnormalities, sometimes with a close anatomical neighborhood between the abnormality and the neoplasm.

CONCLUSIONS

In the first Tumeur Et Développement analysis, 3 major themes have been identified: (1) germline mutations with or without known cancer predisposition, (2) postzygotic events responsible for genomic mosaicism, (3) coincidental associations. New pathways involved in cancer development need to be investigated to improve our understanding of childhood cancers.

Genome-wide mapping of cancer dependency genes and genetic modifiers of chemotherapy in high-risk hepatoblastoma

A lack of relevant genetic models and cell lines hampers our understanding of hepatoblastoma pathogenesis and the development of new therapies for this neoplasm. Here, we report an improved MYC-driven hepatoblastoma-like murine model that recapitulates the pathological features of embryonal type of hepatoblastoma, with transcriptomics resembling the high-risk gene signatures of the human disease. Single-cell RNA-sequencing and spatial transcriptomics identify distinct subpopulations of hepatoblastoma cells. After deriving cell lines from the mouse model, we map cancer dependency genes using CRISPR-Cas9 screening and identify druggable targets shared with human hepatoblastoma (e.g., CDK7, CDK9, PRMT1, PRMT5). Our screen also reveals oncogenes and tumor suppressor genes in hepatoblastoma that engage multiple, druggable cancer signaling pathways. Chemotherapy is critical for human hepatoblastoma treatment. A genetic mapping of doxorubicin response by CRISPR-Cas9 screening identifies modifiers whose loss-of-function synergizes with (e.g., PRKDC) or antagonizes (e.g., apoptosis genes) the effect of chemotherapy. The combination of PRKDC inhibition and doxorubicin-based chemotherapy greatly enhances therapeutic efficacy. These studies provide a set of resources including disease models suitable for identifying and validating potential therapeutic targets in human high-risk hepatoblastoma.

Single-cell transcriptomics uncovers cellular architecture and developmental trajectories in hepatoblastoma

BACKGROUND AND AIMS

Hepatoblastoma (HB) is the predominant type of childhood liver cancer. Treatment options for the clinically advanced HB remain limited. We aimed to dissect the cellular and molecular basis underlying HB oncogenesis and heterogeneity at the single-cell level, which could facilitate a better understanding of HB at both the biological and clinical levels.

APPROACH AND RESULTS

Single-cell transcriptome profiling of tumor and paired distal liver tissue samples from five patients with HB was performed. Deconvolution analysis was used for integrating the single-cell transcriptomic profiles with the bulk transcriptomes of our HB cohort of post-neoadjuvant chemotherapy tumor samples. A single-cell transcriptomic landscape of early human liver parenchymal development was established for exploring the cellular root and hierarchy of HB oncogenesis. As a result, seven distinct tumor cell subpopulations were annotated, and an effective HB subtyping method was established based on their compositions. A HB tumor cell hierarchy was further revealed to not only fit with the classical cancer stem cell (CSC) model but also mirror the early human liver parenchymal development. Moreover, FACT inhibition, which could disrupt the oncogenic positive feedback loop between MYC and SSRP1 in HB, was identified as a promising epigenetic-targeted therapeutic strategy against the CSC-like HB1-Pro-like1 subpopulation and its related high-risk "Pro-like1" subtype of HB.

CONCLUSIONS

Our findings illustrate the cellular architecture and developmental trajectories of HB via integrative bulk and single-cell transcriptome analyses, thus establishing a resourceful framework for the development of targeted diagnostics and therapeutics in the future.

Discovering cryptic splice mutations in cancers via a deep neural network framework

Somatic mutations can disrupt splicing regulatory elements and have dramatic effects on cancer genes, yet the functional consequences of mutations located in extended splice regions is difficult to predict. Here, we use a deep neural network (SpliceAI) to characterize the landscape of splice-altering mutations in cancer. In our in-house series of 401 liver cancers, SpliceAI uncovers 1244 cryptic splice mutations, located outside essential splice sites, that validate at a high rate (66%) in matched RNA-seq data. We then extend the analysis to a large pan-cancer cohort of 17 714 tumors, revealing >100 000 cryptic splice mutations. Taking into account these mutations increases the power of driver gene discovery, revealing 126 new candidate driver genes. It also reveals new driver mutations in known cancer genes, doubling the frequency of splice alterations in tumor suppressor genes. Mutational signature analysis suggests mutational processes that could give rise preferentially to splice mutations in each cancer type, with an enrichment of signatures related to clock-like processes and DNA repair deficiency. Altogether, this work sheds light on the causes and impact of cryptic splice mutations in cancer, and highlights the power of deep learning approaches to better annotate the functional consequences of mutations in oncology.

Intratumoral erythroblastic islands restrain anti-tumor immunity in hepatoblastoma

Erythroblastic islands (EBIs) are the specialized structures for erythropoiesis, but they have never been found functional in tumors. As the most common pediatric liver malignancy, hepatoblastoma (HB) requires more effective and safer therapies to prevent progression and the lifelong impact of complications on young children. However, developing such therapies is impeded by a lack of comprehensive understanding of the tumor microenvironment. By single-cell RNA sequencing of 13 treatment-naive HB patients, we discover an immune landscape characterized by aberrant accumulation of EBIs, formed by VCAM1⁺ macrophages and erythroid cells, which is inversely correlated with survival of HB. Erythroid cells inhibit the function of dendritic cells (DCs) via the LGALS9/TIM3 axis, leading to impaired anti-tumor T cell immune responses. Encouragingly, TIM3 blockades relieve the inhibitory effect of erythroid cells on DCs. Our study provides an immune evasion mechanism mediated by intratumoral EBIs and proposes TIM3 as a promising therapeutic target for HB.

Identification of IGF2 as Genomic Driver and Actionable Therapeutic Target in Hepatoblastoma

Management of hepatoblastoma (HB), the most frequent pediatric liver cancer, is based on surgical resection and perioperative chemotherapy regimens. In this study, we aimed to identify actionable targets in HB and assess the efficacy of molecular therapies in preclinical models of HB. Paired tumor and adjacent tissues from 31 HBs and a validation set of 50 HBs were analyzed using RNA-seq, SNP, and methylation arrays. IGF2 overexpression was identified as the top targetable HB driver, present in 71% of HBs (22/31). IGF2high tumors displayed progenitor cell features and shorter recurrence-free survival. IGF2 overexpression was associated in 91% of cases with fetal promoter hypomethylation, ICR1 deregulation, 11p15.5 loss of heterozygosity or miR483-5p overexpression. The antitumor effect of xentuzumab (a monoclonal antibody targeting IGF1/2) alone or in combination with the conventional therapeutic agent cisplatin was assessed in HB cell lines, in PDX-derived HB organoids and in a xenograft HB murine model. The combination of xentuzumab with cisplatin showed strong synergistic antitumor effects in organoids and in IGF2high cell lines. In mice (n = 55), the combination induced a significant decrease in tumor volume and improved survival compared with cisplatin alone. These results suggest that IGF2 is an HB actionable driver and that, in preclinical models of HB, the combination of IGF1/2 inhibition with cisplatin induces superior antitumor effects than cisplatin monotherapy. Overall, our study provides a rationale for testing IGF2 inhibitors in combination with cisplatin in HB patients with IGF2 overexpression.

Ribonucleotide reductase subunit switching in hepatoblastoma drug response and relapse

Prognosis of children with high-risk hepatoblastoma (HB), the most common pediatric liver cancer, remains poor. In this study, we found ribonucleotide reductase (RNR) subunit M2 (RRM2) was one of the key genes supporting cell proliferation in high-risk HB. While standard chemotherapies could effectively suppress RRM2 in HB cells, they induced a significant upregulation of the other RNR M2 subunit, RRM2B. Computational analysis revealed distinct signaling networks RRM2 and RRM2B were involved in HB patient tumors, with RRM2 supporting cell proliferation and RRM2B participating heavily in stress response pathways. Indeed, RRM2B upregulation in chemotherapy-treated HB cells promoted cell survival and subsequent relapse, during which RRM2B was gradually replaced back by RRM2. Combining an RRM2 inhibitor with chemotherapy showed an effective delaying of HB tumor relapse in vivo. Overall, our study revealed the distinct roles of the two RNR M2 subunits and their dynamic switching during HB cell proliferation and stress response.

Ribonucleotide Reductase Subunit Switching in Hepatoblastoma Drug Response and Relapse

Prognosis of children with high-risk hepatoblastoma (HB), the most common pediatric liver cancer, remains poor. In this study, we found ribonucleotide reductase (RNR) subunit M2 ( RRM2 ) was one of the key genes supporting cell proliferation in high-risk HB. While standard chemotherapies could effectively suppress RRM2 in HB cells, they induced a significant upregulation of the other RNR M2 subunit, RRM2B . Computational analysis revealed distinct signaling networks RRM2 and RRM2B were involved in HB patient tumors, with RRM2 supporting cell proliferation and RRM2B participating heavily in stress response pathways. Indeed, RRM2B upregulation in chemotherapy-treated HB cells promoted cell survival and subsequent relapse, during which RRM2B was gradually replaced back by RRM2. Combining an RRM2 inhibitor with chemotherapy showed an effective delaying of HB tumor relapse in vivo. Overall, our study revealed the distinct roles of the two RNR M2 subunits and their dynamic switching during HB cell proliferation and stress response.

Benign liver tumours: understanding molecular physiology to adapt clinical management

Improvements in understanding the pathophysiology of the different benign liver nodules have refined their nosological classification. New criteria have been identified using imaging, histology and molecular analyses for a precise diagnosis of these tumours. Improvement in the classification of liver tumours provides a more accurate prediction of disease progression and has modified patient management. Haemangioma and focal nodular hyperplasia, the most common benign liver tumours that develop in the absence of chronic liver disease, are usually easy to diagnose on imaging and do not require specific treatment. However, hepatocellular adenomas and cirrhotic macronodules can be difficult to discriminate from hepatocellular carcinoma. The molecular subtyping of hepatocellular adenomas in five major subgroups defined by HNF1A inactivation, β-catenin mutation in exon 3 or exon 7/8, and activation of inflammatory or Hedgehog pathways helps to identify the tumours at risk of malignant transformation or bleeding. New clinical, biological and molecular tools have gradually been included in diagnostic and treatment algorithms to classify benign liver tumours and improve patient management. This Review aims to explain the main pathogenic mechanisms of benign liver tumours and how this knowledge could influence clinical practice.

Hepatoblastomas with carcinoma features represent a biological spectrum of aggressive neoplasms in children and young adults

BACKGROUND & AIMS

Hepatoblastoma (HB) and hepatocellular carcinoma (HCC) are the predominant liver cancers in children, though their respective treatment options and associated outcomes differ dramatically. Risk stratification using a combination of clinical, histological, and molecular parameters can improve treatment selection, but it is particularly challenging for tumors with mixed histological features, including those in the recently created hepatocellular neoplasm not otherwise specified (HCN NOS) provisional category. We aimed to perform the first molecular characterization of clinically annotated cases of HCN NOS.

METHODS

We tested whether these histological features are associated with genetic alterations, cancer gene dysregulation, and outcomes. Namely, we compared the molecular features of HCN NOS, including copy number alterations, mutations, and gene expression profiles, with those in other pediatric hepatocellular neoplasms, including HBs and HCCs, as well as HBs demonstrating focal atypia or pleomorphism (HB FPAs), and HBs diagnosed in older children (>8).

RESULTS

Molecular profiles of HCN NOS and HB FPAs revealed common underlying biological features that were previously observed in HCCs. Consequently, we designated these tumor types collectively as HBs with HCC features (HBCs). These tumors were associated with high mutation rates (∼3 somatic mutations/Mb) and were enriched with mutations and alterations in key cancer genes and pathways. In addition, recurrent large-scale chromosomal gains, including gains of chromosomal arms 2q (80%), 6p (70%), and 20p (70%), were observed. Overall, HBCs were associated with poor clinical outcomes.

CONCLUSIONS

Our study indicates that histological features seen in HBCs are associated with combined molecular features of HB and HCC, that HBCs are associated with poor outcomes irrespective of patient age, and that transplanted patients are more likely to have good outcomes than those treated with chemotherapy and surgery alone. These findings highlight the importance of molecular testing and early therapeutic intervention for aggressive childhood hepatocellular neoplasms.

LAY SUMMARY

We molecularly characterized a class of histologically aggressive childhood liver cancers and showed that these tumors are clinically aggressive and that their observed histological features are associated with underlying recurrent molecular features. We proposed a diagnostic algorithm to identify these cancers using a combination of histological and molecular features, and our analysis suggested that these cancers may benefit from specialized treatment strategies that may differ from treatment guidelines for other childhood liver cancers.

Somatic variation in normal tissues: friend or foe of cancer early detection?

BACKGROUND

Seemingly normal tissues progressively become populated by mutant clones over time. Most of these clones bear mutations in well-known cancer genes but only rarely do they transform into cancer. This poses questions on what triggers cancer initiation and what implications somatic variation has for cancer early detection.

DESIGN

We analysed recent mutational screens of healthy and cancer-free diseased tissues to compare somatic drivers and the causes of somatic variation across tissues. We then reviewed the mechanisms of clonal expansion and their relationships with age and diseases other than cancer. We finally discussed the relevance of somatic variation for cancer initiation and how it can help or hinder cancer detection and prevention.

RESULTS

The extent of somatic variation is highly variable across tissues and depends on intrinsic features, such as tissue architecture and turnover, as well as the exposure to endogenous and exogenous insults. Most somatic mutations driving clonal expansion are tissue-specific and inactivate tumor suppressor genes involved in chromatin modification and cell growth signaling. Some of these genes are more frequently mutated in normal tissues than cancer, indicating a context-dependent cancer promoting or protective role. Mutant clones can persist over a long time or disappear rapidly, suggesting that their fitness depends on the dynamic equilibrium with the environment. The disruption of this equilibrium is likely responsible for their transformation into malignant clones and knowing what triggers this process is key for cancer prevention and early detection. Somatic variation should be considered in liquid biopsy, where it may contribute cancer-independent mutations, and in the identification of cancer drivers, since not all mutated genes favoring clonal expansion also drive tumorigenesis.

CONCLUSIONS

Somatic variation and the factors governing homeostasis of normal tissues should be taken into account when devising strategies for cancer prevention and early detection.

Targeting the Unwindosome by Mebendazole Is a Vulnerability of Chemoresistant Hepatoblastoma

Resistance to conventional chemotherapy remains a huge challenge in the clinical management of hepatoblastoma, the most common liver tumor in childhood. By integrating the gene expression data of hepatoblastoma patients into the perturbation prediction tool Connectivity Map, we identified the clinical widely used anthelmintic mebendazole as a drug to circumvent chemoresistance in permanent and patient-derived xenograft cell lines that are resistant to cisplatin, the therapeutic backbone of hepatoblastoma treatment. Viability assays clearly indicated a potent reduction of tumor cell growth upon mebendazole treatment in a dose-dependent manner. The combination of mebendazole and cisplatin revealed a strong synergistic effect, which was comparable to the one seen with cisplatin and doxorubicin, the current treatment for high-risk hepatoblastoma patients. Moreover, mebendazole treatment resulted in reduced colony and tumor spheroid formation capabilities, cell cycle arrest, and induction of apoptosis of hepatoblastoma cells. Mechanistically, mebendazole causes blockage of microtubule formation and transcriptional downregulation of genes encoding the unwindosome, which are highly expressed in chemoresistant tumors. Most importantly, mebendazole significantly reduced tumor growth in a subcutaneous xenograft transplantation mouse model without side effects. In conclusion, our results strongly support the clinical use of mebendazole in the treatment of chemoresistant hepatoblastoma and highlight the potential theranostic value of unwindosome-associated genes.

Single-cell analysis of hepatoblastoma identifies tumor signatures that predict chemotherapy susceptibility using patient-specific tumor spheroids

Pediatric hepatoblastoma is the most common primary liver cancer in infants and children. Studies of hepatoblastoma that focus exclusively on tumor cells demonstrate sparse somatic mutations and a common cell of origin, the hepatoblast, across patients. In contrast to the homogeneity these studies would suggest, hepatoblastoma tumors have a high degree of heterogeneity that can portend poor prognosis. In this study, we use single-cell transcriptomic techniques to analyze resected human pediatric hepatoblastoma specimens, and identify five hepatoblastoma tumor signatures that may account for the tumor heterogeneity observed in this disease. Notably, patient-derived hepatoblastoma spheroid cultures predict differential responses to treatment based on the transcriptomic signature of each tumor, suggesting a path forward for precision oncology for these tumors. In this work, we define hepatoblastoma tumor heterogeneity with single-cell resolution and demonstrate that patient-derived spheroids can be used to evaluate responses to chemotherapy.

Deleting the β-catenin degradation domain in mouse hepatocytes drives hepatocellular carcinoma or hepatoblastoma-like tumor growth

BACKGROUND & AIMS

One-third of hepatocellular carcinomas (HCCs) harbor mutations activating the β-catenin pathway, predominantly via mutations in the CTNNB1 gene itself. Mouse models of Apc loss-of-function are widely used to mimic β-catenin-dependent tumorigenesis. Given the low prevalence of APC mutations in human HCCs, we aimed to generate liver tumors through CTNNB1 exon 3 deletion (βcat^Δex3). We then compared βcat^Δex3 liver tumors with liver tumors generated via frameshift in exon 15 of Apc (Apc^fs-ex15).

METHODS

We used hepatocyte-specific and inducible mouse models generated through either a Cre-Lox or a CRISPR/Cas9 approach using adeno-associated virus vectors. Tumors generated by the Cre-Lox models were phenotypically analyzed using immunohistochemistry and were selected for transcriptomic analysis by RNA-sequencing (RNAseq). Mouse RNAseq data were compared to human RNAseq data (8 normal tissues, 48 HCCs, 9 hepatoblastomas) in an integrative analysis. Tumors generated via CRISPR were analyzed using DNA sequencing and immuno-histochemistry.

RESULTS

Mice with CTNNB1 exon 3 deletion in hepatocytes developed liver tumors indistinguishable from Apc^fs-ex15 liver tumors. Both Apc^fs-ex15 and βcat^Δex3 mouse models induced growth of phenotypically distinct tumors (differentiated or undifferentiated). Integrative analysis of human and mouse tumors showed that differentiated mouse tumors cluster with well-differentiated human CTNNB1-mutated tumors. Conversely, undifferentiated mouse tumors cluster with human mesenchymal hepatoblastomas and harbor activated YAP signaling.

CONCLUSION

Apc^fs-ex15 and βcat^Δex3 mouse models both induce growth of tumors that are transcriptionally similar to either well-differentiated and β-catenin-activated human HCCs or mesenchymal hepatoblastomas.

LAY SUMMARY

New and easy-to-use transgenic mouse models of primary liver cancers have been generated, with mutations in the gene encoding beta-catenin, which are frequent in both adult and pediatric primary liver cancers. The mice develop both types of cancer, constituting a strong preclinical model.

Treatment optimization for recurrent hepatoblastoma: retrospective study from a hepatoblastoma cohort in Southern China

PURPOSE

The present study aimed to explore the clinical characteristics and optimal treatments of RHB patients.

METHODS

We retrospectively collected 42 RHB cases and 161 primary HB (PHB) cases. Clinical characteristics were compared between RHB and PHB patients. The risk factors related to overall survival (OS) and progression-free survival (PFS) in RHB patients were explored by COX regression analysis. Patients were further divided into curable and refractory subgroup by treatments. Propensity score match (PSM) analysis was performed to match recurrent curable patients from 145 curable PHB patients from the same cohort. PFS was further compared between 34 pairs of primary and recurrent curative HB patients.

RESULTS

Recurrence treatment and number of relapsed tumors were significantly related with both OS and PFS of RHB patients (p < 0.05). Chemotherapy regimen alteration was also risk factor of PFS for RHB (HR = 4.26; 95% CI = 1.54-11.78; p = 0.005). RHB patients underwent curable treatment had better prognosis, compared with recurrent refractory subgroup (p < 0.001). Matched curable PHB patients demonstrated no significant difference of 3-year PFS with curable RHB patients (p = 0.540).

CONCLUSION

Curable RHB patients might get benefit from surgery or ablation with similar prognosis with primary curable HB patients.

Current Approaches in Hepatoblastoma-New Biological Insights to Inform Therapy

PURPOSE OF REVIEW

As the most common pediatric primary liver cancer with rising incidence, hepatoblastoma remains challenging to treat. Here, we review the current understanding of the biology of hepatoblastoma and discuss how recent advances may lead to new treatment modalities.

RECENT FINDINGS

Standard chemotherapy regimens including cisplatin, in addition to surgery, have led to high cure rates among patients with low stage hepatoblastoma; however, metastatic and relapsed disease continue to have poor outcomes. Recent genomics and functional studies in cell lines and mouse models have established a central role for the Wnt/β-catenin pathway in tumorigenesis. Targeted agents and immunotherapy approaches are emerging as potential treatment avenues. With recent gains in knowledge of the genomic and transcriptomic landscape of hepatoblastoma, new therapeutic mechanisms can now be explored to improve outcomes for metastatic and relapsed hepatoblastoma and to reduce the toxicity of current treatments.

miR-126 in Extracellular Vesicles Derived from Hepatoblastoma Cells Promotes the Tumorigenesis of Hepatoblastoma through Inducing the Differentiation of BMSCs into Cancer Stem Cells

Background

Extracellular vesicles (EVs) can deliver miRNAs between cells and play a crucial role in hepatoblastoma progression. In this study, we explored the differentially expressed miRNAs related to tumor cell-derived EVs and the mechanism by which EVs regulate hepatoblastoma progression.

Methods

Bioinformatics analysis was performed to explore the differentially expressed miRNAs between the hepatoblastoma and adjacent normal tissues. TEM, NTA, and western blotting were conducted to identify EVs. The expression of miR-126-3p, miR-126-5p, miR-30b-3p, miR-30b-3p, SRY, IL-1α, IL-6, and TGF-β was detected by RT-qPCR. Immunofluorescence (IF) was used to analyze the expression of PKH67, and flow cytometry was applied to assess the ratio of CD44+ CD90+ CD133+ cells. ELISA was used to evaluate the levels of IL-6 and TGF-β. A xenograft mouse model was constructed to detect the function of EVs with downregulated miR-126. IHC was performed to calculate β-catenin levels in tumor tissues.

Results

miR-126 was upregulated in hepatoblastoma. EVs derived from hepatoblastoma cells significantly increased the ratio of CD44+ CD90+ CD133+ cells and increased the expression of IL-6, Oct4, SRY, and TGF-β in bone marrow mesenchymal stem cells (BMSCs), while EVs with downregulated miR-126 reversed these phenomena. miR-126 downregulation notably attenuated hepatoblastoma tumor growth and decreased the ratio of CD44+ CD90+ CD133+ cells and increased the expression of IL-6, Oct4, SRY, TGF-β, and β-catenin in tumor tissues of mice. Furthermore, EVs with downregulated miR-126 inhibited the differentiation of BMSCs into cancer stem cells.

Conclusions

Exosomal miR-126 derived from hepatoblastoma cells promoted the tumorigenesis of liver cancer through inducing the differentiation of BMSCs into cancer stem cells.

Identification of distinct tumor cell populations and key genetic mechanisms through single cell sequencing in hepatoblastoma

Hepatoblastoma (HB) is the most common primary liver malignancy of childhood, and molecular investigations are limited and effective treatment options for chemoresistant disease are lacking. There is a knowledge gap in the investigation of key driver cells of HB in tumor. Here we show single cell ribonucleic acid sequencing (scRNAseq) analysis of human tumor, background liver, and patient derived xenograft (PDX) to demonstrate gene expression patterns within tumor and to identify intratumor cell subtype heterogeneity to define differing roles in pathogenesis based on intracellular signaling in pediatric HB. We have identified a driver tumor cell cluster in HB by genetic expression which can be examined to define disease mechanism and treatments. Identification of both critical mechanistic pathways combined with unique cell populations provide the basis for discovery and investigation of novel treatment strategies in vitro and in vivo.