Notch activity characterizes a common hepatocellular carcinoma subtype with unique molecular and clinicopathologic features

Notch activity characterizes a common hepatocellular carcinoma subtype with unique molecular and clinicopathologic features Changyu Zhu1,2, Yu-Jui Ho2, Marcela A. Salomao3, Dianne H. Dapito1, Alberto Bartolome1, Robert F. Schwabe1, Ju-Seog Lee4, Scott W. Lowe2,5, Utpal B. Pajvani1,* 1Department of Medicine, Columbia University, New York, NY, USA; 2Department of Cancer Biology and Genetics, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY, USA; 3Department of Laboratory Medicine and Pathology, Mayo Clinic, Scottsdale, AZ, USA; 4Department of Systems Biology, University of Texas MD Anderson Cancer Center, Houston, TX, USA; 5Howard Hughes Medical Institute, Chevy Chase, MD, USA. Abstract Background & Aims: Hepatocellular carcinoma (HCC) is an increasing cause of cancer-related deaths worldwide, due in part through prevalent obesity-related nonalcoholic steatohepatitis (NASH). Hepatocyte Notch pathway is a pathogenic factor to NASH-associated fibrosis, but its role in HCC is less defined. We characterized the molecular and clinical features of Notch-active human HCC, and investigated the mechanisms of how Notch affects NASH-driven HCC. Methods: Using a 14-gene Notch score, we stratified human HCCs in multiple datasets with comprehensive profiling. We performed gene set enrichment analyses comparing Notch-active HCCs to published HCC subtype signatures. Next, we sorted Notch-active hepatocytes from Notch reporter mice for RNA sequencing to study how Notch activation contributes to tumorigenesis in NASH, and characterized Notch-active tumors in a HCC model combining carcinogen and NASH- *Corresponding author: Utpal B. Pajvani, Department of Medicine, Columbia University, Russ Berrie Medical Science Pavilion, 1150 St Nicholas Ave, New York, NY, 10032. [email protected]; fax: (212) 851-5493. Author Contributions: C.Z. designed, performed and interpreted experiments, analyzed data, and wrote the manuscript; Y.J.H. and J.S.L. performed data analyses; M.A.S. performed pathologic reviews of mouse and human HCCs; D.H.D. and A.B. performed and interpreted experiments; R.F.S. and S.W.L. designed and interpreted experiments; U.B.P. designed and interpreted experiments and wrote the manuscript. Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. Conflict of interest statement: The authors declare that they have no competing interest in the work described. Data availability statement: All data supporting the findings of this study are available from the first and corresponding authors upon request. Data and materials availability. RNA sequencing data (Venus mouse hepatocytes and L-NICD mouse HCC) that support the findings of this study have been deposited in the Gene Expression Omnibus under the accession code GSE154553. All other data supporting the findings of this study are available from the first and corresponding authors upon request. HHS Public Access Author manuscript J Hepatol . Author manuscript; available in PMC 2022 March 01. Published in final edited form as: J Hepatol . 2021 March ; 74(3): 613–626. doi:10.1016/j.jhep.2020.09.032. Author Manuscript Author Manuscript Author Manuscript Author Manuscript inducing diet. We used genetic mouse models to manipulate hepatocyte Notch to investigate the sufficiency and necessity of Notch in NASH-driven tumorigenesis. Results: Notch-active signatures were found in ~30% of human HCCs that transcriptionally resemble cholangiocarcinoma-like HCC with lack of activating CTNNB1 (β-catenin) mutations and overall poor prognosis. Endogenous Notch activation in hepatocytes is associated with repressed β-catenin signaling and hepatic metabolic functions, in lieu of increased interactions with extracellular matrix in NASH. Constitutive hepatocyte Notch activation is sufficient to induce β-catenin-inactive HCC in mice with NASH. Notch and β-catenin show a pattern of mutual exclusivity in carcinogen-induced HCC; in this mouse model, chronic blockade of Notch led to the tumor development reliant on β-catenin. Conclusions: Notch activity characterizes a distinct HCC molecular subtype with unique histology and prognosis. Sustained Notch signaling in chronic liver diseases is a driver of tumor formation without acquiring specific genomic mutations. Lay Summary: Notch pathway activity is increased in a subset of liver cancers with poor outcomes, and aberrant Notch activity in obese mice can drive HCC. Graphical Abstract Keywords Hepatocellular carcinoma; Nonalcoholic steatohepatitis; Fibrosis; Notch signaling; β-catenin; Mouse Model Introduction Hepatocellular carcinoma (HCC) is a malignant primary liver cancer, and the fourth leading cause of cancer-related deaths worldwide,1 with a 5-year survival rate of less than 18%.2 HCC typically arises in the setting of chronic liver diseases of various etiologies, with approximately 90% of tumors occurring in fibrotic livers.3 Of chronic liver diseases, obesityZhu et al. Page 2 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript related nonalcoholic steatohepatitis (NASH) is the fastest-growing etiology for HCC in the United States.4 Comprehensive characterizations of large human HCC datasets, including The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC), revealed that HCCs are highly heterogenous with distinct histological and molecular subtypes.5 Only a few genes are frequently mutated (> 5% of HCCs), and aside from “hotspot” mutations in CTNNB1 (β-catenin),6 TP53 ,7 and TERT (Telomerase reverse transcriptase) promoter,8 most mutations are not recurrent.9 Adding to this complexity, recent sequencing efforts reported that heterogeneity also exists at interand intra-tumoral levels in HCC.10 Elucidation of additional molecular mechanisms underlying HCC diversity is necessary for personalized medicine to cope with the growing obese population. Notch signaling is a conserved cell-to-cell communication pathway. Notch ligands (Jagged1/2, and Delta-like ligands 1/3/4) on a signal-sending cell bind Notch receptors (Notch1–4) on a physically proximal signal-receiving cell, resulting in proteolytic cleavage and release of the Notch intracellular domain (NICD). NICD is a transcriptional coactivator that binds to Rbpj and Mastermind (MAM) to activate expression of Notch target genes, including Hairy enhancer of split ( Hes ) and Hes -related Hey family genes.11 Notch is best known for its role in cell fate determination. During liver development, Notch activation in progenitor cells near the portal vein prompts commitment to the biliary lineage.12 Postdevelopment, Notch signaling remains dormant in hepatocytes, but is aberrantly activated in the setting of obesity and NASH as a maladaptive repair signal that facilitates fibrosis development, which can be ameliorated by pharmacologic Notch inhibitors13. As NOTCH1 gain-of-function mutations are frequently associated with human T-cell acute lymphoblastic leukemia (T-ALL),14 we hypothesized that increased hepatocyte Notch activity might push NASH to HCC. But Notch pathway mutations are rarely reported in HCC,15 although this has not been systematically evaluated. More intriguingly, Notch actions in solid tumors is highly context-dependent,16 and there are even reported Notch loss-of-function mutations in some lung cancers,17, 18 suggesting potential tumor suppressive functions of this pathway. Animal studies have been similarly inconclusive, with hepatoblast-specific NICD transgenic mice developing HCC,19, 20 but Notch activity also found to be tumor-suppressive in a HCC model of Retinoblastoma (Rb) pathway deficiency. 21 To clarify these issues, we created a 14-gene transcriptional Notch score to stratify human HCCs. Using this scoring system, we found that Notch-active HCCs were not enriched in somatic mutations in Notch pathways, but were notable for lack of activating CTNNB1 mutations and repressed metabolic processes. Notch-active HCCs expressed a cholangiocarcinoma-like signature with enriched intratumoral fibrosis and had higher tumor stages and worse prognosis than Notch-inactive HCCs. We observed a similar transcriptional profile – decreased hepatocyte metabolic processes and β-catenin signaling, but increased ECM/stromal interactions and growth factor signaling – in Notch-active hepatocytes in reporter mice fed a NASH provoking-diet, as well as fully-penetrant, diet-dependent HCC in mice with constitutive hepatocyte Notch activity. Consistently, in HCC induced by carcinogen and NASH diet-feeding, Notch activity defined the majority of tumors. In this model, Notch inhibition did not impact overall tumor burden, but provoked compensatory βZhu et al. Page 3 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript catenin activity in residual tumors. Overall, these data suggest that NASH-induced persistent Notch activation prompts tumorigenesis even in the absence of genomic driver events such as activating mutations in CTNNB1 or Notch pathway components. Materials and Methods Human HCC datasets. For analyses of human HCC datasets (TCGA, NCI, Korean, and Fudan cohorts), please see Supplementary Materials and Methods. Mouse Studies. Notch-Venus,22 homozygous RosaNICD ,23 and RosaDNMAM 24 male mice were crossed with female C57BL/6J (strain #000664, the Jackson Laboratory) mice to generate heterozygous transgenic mice for experiments. Mice were weaned to standard chow (PicoLab, #5053) in all experiments. For the DEN NASH model, 8-week-old male C57BL/6J mice were first treated with 100 mg/kg diethylnitrosamine (DEN, Sigma #73861), then started on NASH diet (Teklad, TD.190142) with fructose-containing drinking water (23.1 g fructose and 18.9 g glucose dissolved in 1 L water, then filter-sterilized) two weeks later. 8-week-old RosaNICD or RosaDNMAM mice were transduced by tail vein injections with 1.5×1011 genome copies/mouse of Adeno-associated virus subtype 8 (AAV8) expressing hepatocytespecific Cre recombinase (AAV-TBGCre , 107787-AAV8) or control (AAV-TBGLacZ , 105534-AAV8) from Addgene, to generate Cre− control, L-NICD, and L-DNMAM mice. After transduction, mice were maintained on chow diet or started on NASH or high-fat (Research Diets, D12492) diet for 16 weeks. Animals were housed in standard cages at 22°C in a 12-hour light and 12-hour dark cycle, and monitored for overall well-being and signs of distress with body weight measured weekly. Upon completion of each studies, all mice were weighed and euthanized. Numbers and diameters of macroscopically visible tumors were recorded, and livers were weighed, excised and split for histology and molecular characterization. Specifically, tumors (diameter ≥ 5 mm) and surrounding non-tumor tissues were carefully separated for comparison including qPCR and RNA sequencing. The Columbia University Institutional Animal Care and Use Committee approved all animal procedures. Statistical analysis. Data are shown as mean ± SEM unless indicated otherwise. Statistical analysis was performed using Prism software (version 8, GraphPad Software). Contingency comparison was done by Fisher exact tests (2 groups) or Chi-square tests (3 or more groups). Survival analyses were done by log-rank tests. Differences between 2 groups were calculated using 2tailed Student’s t tests. Analysis involving multiple groups was performed using one-way ANOVA, followed by Dunnett’s multiple comparison tests for two specific groups. P < 0.05 was considered statistically significant. For further details, please refer to CTAT table and Supplementary Materials and Methods. Zhu et al. Page 4 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript Results Notch signaling is transcriptionally active in a subset of human HCCs that lack activating CTNNB1 mutations. To define the relationship between tumor Notch activity and HCC progression, we applied a 14-gene Notch transcriptional signature score (hence, Notch score) to the TCGA HCC dataset.9 Specifically, we calculated the z-score average of Notch receptors ( NOTCH1 , 2 , 3 , 4 ) and ligands ( JAG1 , JAG2 , DLL1 , DLL4 ), core transcription coactivator ( MAML1 ), canonical targets ( HES1 , HEY1 , HEYL ) and liver-enriched downstream effectors ( SOX9 , SPP1 ) in 373 HCCs of mixed etiologies with available transcription data (Fig. S1A), and defined three groups: Notch-low (29.5%), -mid (41%), and -high (29.5%) (Fig. 1A). Of note, Notch-high HCCs were observed at a similar frequency to the percentage (31.8%) of human hepatitis C virus-related HCCs reported as having an active Notch signature.19 We next defined molecular characteristics of Notch-high HCCs. Compared to Notch-inactive groups, Notch-high HCCs are mostly (~90%) deprived of activating mutations of CTNNB1 6 (Fig. 1B, C), a frequently mutated oncogene in HCC that activates a unique set of downstream genes associated with periportal hepatocytes, including GLUL , AXIN2 , RHBG , and NOTUM 25 (Fig. S1B). To validate these data, we next applied our Notch score to an additional two cohorts with both genomic and transcriptomic data – LIRI-JP26 and LICAFR27 from ICGC (Fig. S1C–F) – both of which showed a lower number of activating CTNNB1 mutations in Notch-high tumors (Fig. 1D–G). The inverse correlation between Notch activity and CTNNB1 mutations are reminiscent of opposing actions of Notch and Wnt/β-catenin in development28 and liver injury response,29 suggesting two distinct mechanistic routes of HCC development. Notch genomic mutations are not causal to transcriptional Notch activity in HCC. Although activating NOTCH1 mutations are found in over 50% of T-ALL patients,14 the role of Notch in solid tumors is highly context-dependent with both gainand loss-offunction mutations reported,16 though none thus far in HCC.15, 19 We queried 1221 HCC samples from five different datasets on cBioPortal, and found that ~10% of tumors harbor Notch pathway somatic mutations, with the majority in Notch receptor genes (Fig. S2A). To investigate the correlation between these mutations and tumoral Notch activity, we grouped these mutations based on the Notch score in TCGA, LIRI-JP and LICA-FR datasets (Table S1–3). This analysis revealed primarily missense mutations without recurrent “hotspots” associated with known oncogenic or tumor suppressive functions, with inconsistent relationship between Notch score and mutations (Fig. S2B), suggesting likely passenger mutations in genes of relatively large size. We did find several mutations, especially in Notch receptor genes (Fig. S2B) that could lead to Notch gain or loss of function, and theoretically drive HCC. For instance, the frameshift insertion mutation E2143Lfs*5 in the last exon of NOTCH2 results in loss of the C-terminal proline-glutamic acid-serinethreonine-rich (PEST) domain responsible for ubiquitination and degradation of NICD, resulting in persistent NOTCH2-ICD activation.30 Other truncations (G1453Vfs*100 and A1917Efs*13) likely result in dominant-negative proteins that bind to ligands without transcriptional capabilities.31 But these data overall suggest that somatic mutations in the Zhu et al. Page 5 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript Notch pathway that may drive HCC are rare, whereas the 14-gene Notch score revealed far greater frequency of Notch activity in tumors. Notch-active HCCs are cholangiocarcinoma-like and enriched with intratumoral fibrosis. Next, we compared transcriptomics from Notch-high and -low HCCs. Gene set enrichment analysis (GSEA) further confirmed that Notch-high HCCs display a previously identified Notch-induced liver cancer signature19 (Fig. S3A), while Notch-low HCCs fit the CTNNB1 activation subclass32 (Fig. S3B). Canonical β-catenin target genes associated with the presence of activating CTNNB1 mutations, including GLUL , AXIN2 , NOTUM , are downregulated in Notch-high HCCs; whereas genes related to extracellular matrix (ECM) and myofibroblast markers ( PDGFRA/B and VIM ) are highly induced in these tumors (Fig. 2A). Pathway analysis confirmed that the top pathways upregulated in Notch-high HCCs are related to ECM deposition, while downregulated pathways are primarily metabolic, suggesting that Notch-high HCCs are less differentiated and more fibrotic (Fig. 2B and Fig. S3C). Consistently, Notch-high tumors show marked intratumoral fibrosis (Fig. 2C), and enrichment of stroma transcription signatures33, 34 (Fig. S3D). Notch-high HCCs also showed increased expression of cholangiocyte/progenitor markers ( KRT19 and DCLK1 ) (Fig. 2A), suggestive of cholangiocarcinoma (CC)-type features. CCs also show a highly desmoplastic stroma,3 and these data are reminiscent of the role of Notch in promoting biliary lineage during liver development and mouse data suggesting potential of Notch to transform mouse hepatocytes into intrahepatic CC (ICC).35 Thus, we applied a CC-like HCC transcriptional signature, which marks a HCC subtype with similarly poor outcomes as combined HCC-CC and CC despite lacking the overall histopathological features of CC36, and found that over 50% of Notch-high tumors fit this signature in TCGA HCC datasets (Fig. 2D), while GSEA using another signature of Cytokeratin 19 (KRT19)- positive HCCs37 demonstrate similar results (Fig. S3E). Consistently, expression of Notch pathway genes, β-catenin targets, and hepatocyte markers of Notch-high HCC fall between Notch-low HCC and ICC (Fig. S4A–C), with some notable differences in expression of HCC markers (i.e. AFP and GPC3 ) which are far higher in Notch-high HCC compared to ICC (Fig. S4D). Consistent with the RNA seq results, Notch-high HCCs have higher expression of genes related to cholangiocyte/progenitor cell identity and ECM/fibrosis than Notch-low HCC (Fig. S4E–F), with a similar pattern observed in LICA-FR and LIRI-JP cohorts (Fig. S5A–F). These data suggest that Notch activity marks HCCs with CC features. Notch activity is associated with HCC subtypes of higher cancer stage and poor prognosis. We next compared the transcriptomic signature of Notch-high HCCs to previously defined molecular subtypes. Upregulated genes from Notch-high HCCs are highly enriched in Subclass 1 (S1), a subtype associated with TGF-β signaling, less differentiated tumors and poor survival; conversely, downregulated genes are enriched in S3, a subtype of welldifferentiated tumors defined by activating CTNNB1 mutations38 (Fig. S6A). Notch-high HCCs similarly overlaps with the transcriptional signature of Subgroup 1 (G1), a subtype enriched in genes of developmental signaling and fetal liver39 (Fig. S6B). Notch-high HCCs Zhu et al. Page 6 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript also share molecular homology with a Proliferation subclass characterized by elevated serum a-fetoprotein (AFP) and chromosomal instability32 (Fig. S6C). These data suggested that Notch-high HCCs would be associated with worse outcome. To test this, we applied TCGA iCluster to categorize each tumor. Consistent with our other analyses, Notch-low HCCs primarily fell in iClusters 2 and 3 characterized by CTNNB1 mutations, while > 65% of Notch-high tumors clustered with iCluster 1 (Fig. S6D), which have the worst survival outcomes.9 We performed analogous comparisons of TCGA samples to three different classification signatures, and found that Notch-high HCCs are predicted to show progenitor features40 and poorer survival41 with greater risks of metastasis42 (Fig. S6E–G). Next, we compared the clinical data based on Notch activity. Although subtypematching analyses suggest that patients with Notch-high HCCs would have worse prognosis, and Notch-high group has a higher percentage of advanced stage tumors, with a trend towards higher levels of serum AFP, overall and recurrence-free survival curves were not different (Fig. S6H–K), which may be due to the relatively short follow-up times of TCGA cohort.9 We evaluated another three long-term follow-up HCC cohorts with clinical characterization and microarray-based expression data,40–42 where we observed that higher Notch activity in the tumors (Fig. S7A–C) is associated with worse survival (Fig. 2E–G) and more frequent relapse (Fig. S7D–F), and Notch-high group in general displayed higher tumor stages, larger tumor sizes, and higher AFP levels (Fig. S7G–N). Overall, these data from human HCC datasets demonstrate that Notch-high HCCs adopt a CC-like histology with stromal fibrosis, and are associated with advanced tumor stages and generally poorer outcomes. On the contrary, Notch-inactive HCCs are well-differentiated neoplasms, and typically associated with CTNNB1 mutations that portend better prognosis. 25 Notch signaling in hepatocytes promotes their loss of identity but enhances stromal activity. We next asked whether hepatocyte Notch activation, known to activate hepatic stellate cells (HSC) to trigger liver fibrosis in NASH,13 can directly contribute to NASH-induced HCC development. First, using a Notch-Venus reporter mouse,22 we found that liver Notch signaling is generally confined to nonparenchymal cells (Fig. S8A), but activated in ~15% hepatocytes13 in diet-induced NASH (Fig. 3A). To understand the impact of this, we isolated Notch-active and -inactive hepatocyte populations from NASH diet-fed reporter mice (Fig. 3B) and performed RNA sequencing. Notch-active hepatocytes showed increased expression of pathways that mediate interactions with ECM, consistent with non-cell autonomous activation of stromal cells; conversely, metabolic pathways and β-catenin signaling were repressed, consistent with previous finding43 that Notch reprograms hepatocytes to progenitor/cholangiocyte-like cells (Fig. 3C, D and Fig. S8B). Intriguingly, the transcriptional signatures of Notch-active hepatocytes strongly overlapped with human Notch-high HCCs (Fig. 3E and Fig. S8C) and mouse data from Notch overexpression studies (Fig. S8D), suggesting that Notch-induced change of hepatocyte identity and remodeling of ECM may represent a mechanism that persists during the tumorigenic progression from NASH to HCC. Zhu et al. Page 7 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript Notch signaling is active in a subset of tumors in a mouse model of HCC induced by carcinogen and NASH diet. To test whether Notch is active in mouse HCC, we combined the carcinogen diethylnitrosamine (DEN) with NASH diet-feeding (henceforth, DEN-NASH) in wildtype (WT) mice, similar to prior studies,44 which led to progressive development of multifocal liver tumors (Fig. S9A, B). These tumors were histologically HCC, with varying morphologies including steatotic and steatohepatitic features which are more commonly described in patients with NASH45 (Fig. S9C). We applied this experimental paradigm to Notch-Venus reporter mice (Fig. 4A), which prompted robust tumorigenesis (Fig. 4B) with elevated HCC markers (Fig. 4C). By co-staining with A6 antiserum that marks some mouse HCC,46 we observe both Notch-inactive and -active HCCs (Fig. 4D). Given our data showing Notch and β-catenin mutual exclusivity observed in human HCC, we predicted this would hold true in this mouse model. We found that over half of the tumors are Notchactive, while immunoreactivity of tumors with β-catenin target glutamine synthetase (GS, encoded by Glul )46 was low in DEN-NASH HCC (~10%), and Notch/Venus and GS signals are mutually exclusive (Fig. 4E), even in non-tumor areas (Fig. S10A). We confirmed this result with Osteopontin (Opn, encoded by Spp1 ) staining, a sensitive readout of Notch activity in hepatocytes13 (Fig. S10B–D). These data suggest that similar to human HCC, Notchand β-catenin-active HCCs are fundamentally different in mice. Notch activity is sufficient to drive tumorigenesis in the context of NASH. To determine if hepatocyte Notch activation is sufficient to drive HCC and model Notchhigh HCCs in humans, we transduced young adult male RosaNICD mice23 with AAV8-TbgCre (henceforth, L-NICD mice) to overexpress NICD specifically in hepatocytes,47 followed by 4-month exposure to the NASH diet (Fig. 5A). NASH diet-fed L-NICD mice developed fully penetrant, multifocal large liver tumors in both male (Fig. 5B–D) and female mice (Fig. S11A–C). In contrast, L-NICD mice fed regular chow or a similarly obesogenic high-fat diet did not develop HCC within the same timeframe (Fig. S11D). Histologically, most of the tumors were well-differentiated and steatotic HCCs (Figure 5E), with loss of reticulin fibers (Fig. 5F). Tumors stained positive for the hepatocyte marker HNF4α; consistently, taking advantage of Cre-inducible nuclear eGFP gene in the same locus as NICD for lineage tracing,23 we confirmed that tumor cells derived from hepatocytes (Fig. S11E). L-NICD tumors also stained positive for the liver progenitor antigen A6 (Fig. 5G) and proliferation marker Ki67 (Fig. 5H), but negative for AFP (Fig. S11F), consistent with well-differentiated and proliferative HCCs with some immature features. Finally, consistent with the mutual exclusivity between Notch and β-catenin signaling, Glul expression was downregulated in L-NICD livers and further decreased in the tumors (Fig. S11G–H), leading to the complete lack of GS staining in the tumors (Fig. 5I). Overall, these findings indicate that constitutive Notch activation in hepatocytes leads to Notch-active, β-catenin-deficient HCCs, but only in the fibrotic niche evoked by NASH diet-feeding, consistent with human HCC development in the context of chronic liver diseases. Zhu et al. Page 8 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript NASH diet-fed L-NICD mice recapitulates major aspects of human Notch-active HCC. To determine if HCCs in NASH diet-fed L-NICD mice model Notch-high human HCC, we performed RNA sequencing of the mouse tumors. Compared to surrounding non-tumor (NT) liver, HCCs from L-NICD mice showed upregulation of pathways related to cell cycle and proliferation, with increased expression of the minichromosome maintenance protein complex (MCM) family of genes that are critical for DNA replication and dysregulated in cancers,48 and downregulation of liver metabolic functions (Fig. 6A, B and Fig. S12A). In order to select appropriate human Notch-high HCCs from TCGA for comparison, we first calculated the Notch scores in the surrounding normal liver tissues, and selected Notch-high NT samples (Fig. S12B, C) to parallel increased Notch activity in NT and tumor areas in LNICD mice. Differential gene expression analysis of human Notch-high HCCs normalized to NT (Fig. 6C) showed enrichment of altered pathways similar to L-NICD HCCs (Fig. 12D, E). GSEA revealed that L-NICD tumors transcriptionally resemble the signatures from human Notch-high HCCs (Fig. 6D), as well as hepatocytes with endogenous Notch activity isolated from Notch reporter mice fed on NASH diet (Fig. S12F), ruling out the possible supraphysiologic effects due to NICD overexpression in this model. GSEA of established molecular signatures of HCC subtypes found that mouse L-NICD HCCs and human Notch-high HCCs share the signatures of both KRT19-expressing HCC37 (Fig. S13A) and Proliferation subclass32 (Fig. S13B), especially the subgroups (G1-G3 groups) with chromosomal instability and overexpression of cell cycle-related genes39 (Fig. S13C). However, unlike human Notch-high HCCs, L-NICD HCCs resemble S2, a subclass of HCC related to high rate of proliferation, rather than S138 (Fig. S13D), possibly due to lower expression of stromal ECM-related genes (Fig. S13E). Indeed, in comparison to human Notch-high HCCs, tumors from NASH diet-fed L-NICD mice did not show increased intratumoral fibrosis or numbers of hepatic stellate cell-derived myofibroblasts (Fig. S14A–D) as compared to surrounding non-tumor tissues. Finally, we probed the mutational landscape of 9 HCC samples from NASH diet-fed L-NICD mice using the MSK mouse Integrated Mutation Profiling of Actionable Cancer Targets (M-IMPACT) assay.49 M-IMPACT sequencing found no activating mutations of Ctnnb1 , but unique somatic mutations in 27 genes (Fig. S14E and Table S4), including recurrent mutations in epigenetic regulators Kmt2a and Kmt2d ,50 as well as Amer1 which is a negative regulator of β-catenin stability.51 Chronic inhibition of Notch results in switch of molecular subtypes in mouse HCC. These data establish a potential mouse model for Notch-active HCC, but also that hepatocyte Notch activation is sufficient to promote tumorigenesis in the context of NASH. To investigate if Notch signaling is necessary for HCC development, we transduced DENtreated, NASH diet-fed RosaDNMAM mice carrying a dominant-negative Mastermind (DNMAM) allele24 with AAV8-TbgCre to block all Notch receptors-mediated transcription in hepatocytes (L-DNMAM mice) (Fig. 7A). At the end of the experiment, we found no difference in tumor number or size between the two groups (Fig. S15A–C), despite a reduction in liver fibrosis (Fig. S15D, E), maintained expression of the Dnmam transgene (Fig. S15F) and lowered Notch targets (Fig. 7B and Fig. S15G), including lower tumor Opn Zhu et al. Page 9 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript Highlights •Notch activation induces liver fibrosis, an important niche factor for tumorigenesis •Notch is active in 30% of human hepatocellular carcinomas (HCCs) with poor prognosis •Notch-active HCCs are fibrotic, and have low β-catenin activity in mice and humans. •Forced hepatocyte Notch activity in mice induces diet-dependent HCC. Zhu et al. Page 16 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript Fig. 1. Notch signaling is active in a subset of HCC characterized by a lack of activating CTNNB1 mutations. (A) Heatmap of mRNA expression z-scores of 14 Notch pathway-related genes that collectively define the Notch Score, based on which the 373 tumors in the TCGA HCC dataset were sub-divided into Notch-low, -mid and -high groups. (B to G) Frequencies of commonly found somatic mutations and recurrent activating CTNNB1 mutations in (B and C) TCGA, (D and E) LIRI-JP and (F and G) LICA-FR HCC datasets, based on Notch score. *, P < 0.05, **, P < 0.01 and ***, P < 0.001, calculated by Chi-square tests. Zhu et al. Page 17 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript Fig. 2. Notch-high HCCs show unique molecular and histologic features, and lower survival. (A) Volcano plot showing differentially expressed genes (DEGs, adjusted P -value < 0.05) in Notch-low and -high HCCs (n = 110 per group) from the TCGA dataset. (B) Upregulated and downregulated pathways from KEGG 2019 in Notch-high HCCs, ranked by combined enrichment scores. (C) Representative H&E images of Notch-low, -mid and -high tumors. (D) Percentage of cholangiocarcinoma (CC)-like HCC in Notch-low, -mid and -high tumors in the TCGA dataset. (E) Kaplan-Meier survival curves of patients with Notch-low, -mid and -high tumors in NCI, (D) Korean, and (G) Fudan cohorts. P -values calculated by log-rank tests. Zhu et al. Page 18 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript Fig. 3. Notch activation is associated with reduced hepatocyte metabolic function, but a signature of intercellular crosstalk with stromal cells. (A) Representative images of Venus (green) and HNF4α (red) co-staining in chowand NASH diet-fed Notch-Venus reporter mice. (B) Flow cytometry plot and (C) Volcano plot showing DEGs in Notch-active and -inactive hepatocytes (n = 3 per group) isolated from NASH diet-fed Notch-Venus reporter mice. (D) Upregulated and downregulated pathways from KEGG 2019 in Notch-active hepatocytes, ranked by combined enrichment scores. (E) Gene set enrichment analyses (GSEA) of Notch-active hepatocyte signatures in human Notch-high or low HCCs from TCGA. UP, upregulated. DN, downregulated. NES, normalized enrichment scores. FDR, false discovery rate. Zhu et al. Page 19 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript Fig. 4. Notch and β-catenin activity mark distinct carcinogen-induced HCC subsets. (A) Notch-Venus reporter mice were dosed with diethylnitrosamine (DEN, 100 mg/kg), then fed NASH diet for 10 months prior to sacrifice. (B) Representative images of liver tumors, and (C) Non-tumor and Tumor expression of HCC marker genes in Notch-Venus mice (n = 8 per group). (D) Representative images of A6 (red), Venus (green), and DAPI (blue) costaining in Notch-negative and –positive tumors from Notch-Venus mice. (E) Representative images of β-catenin target glutamine synthetase (GS, red) and Venus (green) co-staining in GS+ and Venus+ tumors from Notch-Venus mice, with quantification of GS+ and/or Notch/ Venus+ tumors. *, P < 0.05 and **, P < 0.01 as compared to the non-tumor group by twotailed t tests. AU, arbitrary unit. Data are shown as means ± SEM. Zhu et al. Page 20 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript Fig. 5. Hepatocyte Notch activation induces HCC in NASH diet-fed mice. (A) Male RosaNICD mice were transduced with AAV8-TbgLacZ (control) or AAV8-TbgCre (L-NICD) and then fed with NASH diet for 16 weeks. (B) Liver-to-body weight ratios and (C) tumor number in control and L-NICD mice (n = 9–11 per group). (D) Representative images of livers from control and L-NICD mice, with (E) H&E and (F) reticulin staining in a representative tumor from L-NICD mice. (G) Representative images of HCC antigen A6 (red) co-staining with HNF4α (green), (H) Ki67 (green), and (I) GS (green) in tumors from L-NICD mice. **, P < 0.01 and ***, P < 0.001 as compared to the control group by two-tailed Student’s t tests. Data are shown as means ± SEM. Zhu et al. Page 21 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript Fig. 6. NASH-diet fed L-NICD mice transcriptionally recapitulates human Notch-active HCC. (A) Volcano plot showing DEGs in tumors as compared to non-tumor surrounding tissues (n = 3 per group) in NASH diet-fed L-NICD mice. (B) Upregulated and downregulated pathways from KEGG 2019 in L-NICD HCCs, ranked by combined enrichment scores. (C) Volcano plot showing DEGs and (D) GSEA in human Notch-high HCCs as compared to adjacent non-tumor liver tissues (n = 10 per group) from the TCGA dataset. Zhu et al. Page 22 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript Fig. 7. Chronic Notch inhibition in carcinogen-treated mice leads to β-catenin-dependent HCC. (A) RosaDNMAM mice were dosed with DEN (100 mg/kg), fed NASH diet, then transduced with AAV8-TbgLacZ (control) or AAV8-TbgCre (L-DNMAM) (n = 9–10 per group). (B) L-DNMAM tumors show reduced expression of the Notch target Spp1 , (C) with correspondingly lower tumor co-staining of Opn (red). (D) L-DNMAM tumors show increased tumoral expression of the β-catenin target Glul , (E) with correspondingly increased tumor co-staining with GS (green). (F) Quantification of Notch/Opn+ and GS+ tumors. *, P < 0.05 and **, P < 0.01 as compared to the indicated control by two-tailed t tests. AU, arbitrary unit. Data are shown as means ± SEM. Zhu et al. Page 23 J Hepatol . Author manuscript; available in PMC 2022 March 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript