β-Catenin mutations as determinants of hepatoblastoma phenotypes in mice

Hepatoblastoma (HB) is the most common pediatric liver cancer. Although long-term survival of HB is generally favorable, it depends on clinical stage, tumor histology, and a variety of biochemical and molecular features. HB appears almost exclusively before the age of 3 years, is represented by seven histological subtypes, and is usually associated with highly heterogeneous somatic mutations in the catenin β1 (CTNNB1) gene, which encodes β-catenin, a Wnt ligand–responsive transcriptional co-factor. Numerous recurring β-catenin mutations, not previously documented in HB, have also been identified in various other pediatric and adult cancer types. Little is known about the underlying factors that determine the above HB features and behaviors or whether non-HB–associated β-catenin mutations are tumorigenic when expressed in hepatocytes. Here, we investigated the oncogenic properties of 14 different HB– and non-HB–associated β-catenin mutants encoded by Sleeping Beauty vectors following their delivery into the mouse liver by hydrodynamic tail-vein injection. We show that all β-catenin mutations, as well as WT β-catenin, are tumorigenic when co-expressed with a mutant form of yes-associated protein (YAP). However, tumor growth rates, histologies, nuclear-to-cytoplasmic partitioning, and metabolic and transcriptional landscapes were strongly influenced by the identities of the β-catenin mutations. These findings provide a context for understanding at the molecular level the notable biological diversity of HB.

Hepatoblastoma (HB) 6 is the most common pediatric liver cancer, typically arising in children less than 3 years of age (1,2). Various staging systems have been proposed to allow the assessment of therapeutic options and prognosis (3)(4)(5)(6)(7). HB has been subclassified into seven distinct histologic categories (8) with the small-cell undifferentiated variety reportedly having a less favorable outcome (3)(4)(5). Other features associated with inferior long-term survival include older age at the time of diagnosis, prematurity, low levels of circulating ␣-fetoprotein, and high levels of and/or more intense nuclear localization of ␤-catenin in the tumor (2,4,5,9). Predisposing genetic disorders such as Beckwith-Wiedeman syndrome and familial adenomatous polyposis (FAP) also significantly increase the relative risk of HB development (10 -12).
Pediatric cancers in general and HBs in particular have mutational burdens as much as 10 4 -fold lower than those of adult cancers (13)(14)(15). Nonetheless, Ͼ80% of HBs contain somatically acquired CTNNB1 gene mutations, and at least one of these has demonstrable in vivo tumorigenicity in mice (3,(15)(16)(17)(18)(19)(20)(21)(22)(23)(24)(25)(26). CTNNB1 encodes ␤-catenin, a transcriptional co-activator that mediates signaling via the Wnt pathway. In quiescent cells, ␤-catenin resides in the cytoplasm in a transcriptionally inert form as a short-lived member of the APC complex wherein it is sequentially phosphorylated on several neighboring N-terminal serine and threonine residues and thus marked for rapid polyubiquitinylation and proteasome-mediated degradation (27)(28)(29). Extracellular Wnt ligand signaling recruits the APC complex to the inner surface of the plasma membrane, causing the release of unphosphorylated ␤-catenin that cannot be recognized by the degradation complex. The now stabilized ␤-catenin translocates to the nucleus, where it associates with Tcf/Lef transcription factor family members, engages in transcriptional cross-talk with the Hippo tumor suppressor pathway, and drives tumorigenesis by deregulating a large array of targets, including such well-known oncogenes as CCND1 and MYC (3,24,30,31). Mutant forms of ␤-catenin, often involving one or more of its phosphorylation sites, prevent its association with the APC complex, its ubiquitination, and its proteasomal degradation. Instead, the now stabilized mutant ␤-catenin is constitutively redirected to the nucleus in a Wnt-independent manner. A small subset of HBs with WT ␤-catenin harbor acquired mutations in other APC complex components, including Axin1, Axin2, and APC itself, thus leading to the same failure of ␤-catenin phosphorylation, ubiquitinylation, and rapid proteasomal turnover (2,3,22,27,32). Indeed, it is the germ-line mutational inactivation of APC in FAP, leading to the development of polyposis and colo-rectal cancer (CRC) later in life, that is responsible for the greatly increased risk of HB tumorigenesis during early childhood in these same individuals (10,12). Whereas the majority of sporadic CRCs possess acquired APC mutations, those that do not often harbor ␤-catenin mutations identical to or distinct from those described in HBs (33,34). These findings underscore the intimate molecular relationship between these seemingly disparate pediatric and adult cancers.
Unlike most oncogenes, whose mutational repertoires tend to be quite restricted, ␤-catenin mutations are highly heterogeneous, comprising both missense substitutions and variably sized in-frame deletions (21,32,35,36) (http://ongene.bioinfominzhao.org/gene_mutation.cgi and https://cancer.sanger.ac. uk/cosmic/gene/analysis?lnϭCTNNB1). 7 Most occur within the ϳ20-amino acid region encompassing the above-mentioned phosphorylation sites. However, a number of less frequent ones reside elsewhere, particularly within the so-called armadillo repeat region (residues ϳ141-664), which engages with important cooperating proteins, such as E-cadherin, Tcf4/ Lef, and APC (37)(38)(39). These observations suggest that ␤-catenin's oncogenic activation may have different mutantspecific consequences and be associated with distinct tumor features, behaviors, and outcomes. Recurrent ␤-catenin mutations are also found in other cancers, including but not restricted to hepatocellular carcinoma (HCC), medulloblastoma, and ovarian and prostate cancer, with many of the mutations not having been previously described in HBs (40 -43).
The Hippo pathway is also deregulated in most HBs, although the precise nature of the defect(s) remains obscure (3,23). Most HBs, HCCs, cholangiocarcinomas, and CRCs show high-level nuclear localization of yes-associated protein (YAP), a terminal effector of the Hippo pathway (3,23,44,45), and a nucleus-localized mutant form of YAP, YAP S127A , is required to promote HB tumorigenesis by mutant ␤-catenin (23). Additional mutations in NFE2L2/NRF2, TERT, TP53, and other genes have been identified in HB but at frequencies of Ͻ10% (13,14,18). More recently, we showed that HBs induced by a single ␤-catenin mutant known as ⌬(90) (23) together with YAP S127A strongly down-regulated the p19 ARF tumor suppressor and that enforced p19 ARF expression totally abrogated tumorigenesis (46). This confirmed previous findings that human and murine HBs seldom contain TP53 mutations (46,47) but showed that upstream compromise of the pathway may nonetheless occur via p19 ARF silencing.
The foregoing observations raise several important questions, including the extent to which different ␤-catenin mutations dictate tumor growth rates and histology. Do they also differentially affect tumor metabolism, which is dramatically altered from that of the normal liver (20,24,25)? More directly, how do these mutations affect ␤-catenin expression, its nuclear: cytoplasmic partitioning, and tumor-specific gene expression profiles (15,19,24,25,48)? Finally, can non-HB-associated ␤-catenin mutants generate liver tumors, and, if so, do they differ in any discernible ways from those caused by HB-associated mutants? Answers to these questions could also explain the degree to which the mutant identities influence tumor stage, chemotherapeutic response, the likelihood of recurrence, and the prospects for long-term survival and cure.
To address these issues, we relied here on the above-described mouse model of HB in which Sleeping Beauty (SB) vectors encoding mutant forms of ␤-catenin and YAP S127A are stably expressed in murine livers following their delivery by hydrodynamic tail-vein injection (HDTVI) (20,(23)(24)(25)(26). We compared tumors generated by each of 14 patient-derived ␤-catenin variants and WT ␤-catenin protein. The selected mutants were a combination of those associated with HBs only, those associated with tumors other than HBs, and those shared by HBs and other tumors. We found tumor growth rates, histologic subtypes, ␤-catenin levels, and nuclear:cytoplasmic distribution, along with metabolic and transcriptional profiles to be profoundly influenced by the identity of the ␤-catenin mutant. Those previously identified only in non-HB tumors as well as WT ␤-catenin were also shown to be oncogenic, although the latter was associated with low-level expression, slower growth rate, and a marked HCC-like histology. Our studies identify the extent to which numerous tumor properties are influenced by the nature of the driver ␤-catenin mutant and provide a rationale for considering ␤-catenin mutational status in the initial clinical evaluation.

␤-Catenin mutants determine HB growth rates, survival, and histology
The previously described ⌬(90) ␤-catenin mutant, encoding an N-terminal, in-frame 90-amino acid deletion was used as the standard against which all other ␤-catenin mutants were compared (23)(24)(25)49). These consisted of 13 additional human ␤-catenin mutants in the same SB vector backbone (23-25, 49) ( Fig. 1A and Table S1). All mutants were patient-derived and recurrent and were selected based on their location within the ␤-catenin molecule, whether they were missense or in-frame deletions, and whether they had been previously described in HBs only, non-HB tumors only, or both HBs and non-HB tumors. ⌬(90) was Myc epitope-tagged, but all others were C-terminally His 6 -tagged to allow for comparisons of their 7 Please note that the JBC is not responsible for the long-term archiving and maintenance of this site or any other third party hosted site.
␤-Catenin mutations and hepatoblastoma phenotypes expression levels. Purified plasmid DNAs, along with another SB vector encoding YAP S127A and a vector encoding SB translocase, were introduced via HDTVI into 10 -15 FVB mice (23-25, 46, 49). The vast majority of animals developed tumors, and Kaplan-Meier survival curves identified three distinct groups ( Fig. 1B and Fig. S1A). ⌬(90)-associated tumors, along with those generated by five additional mutants, comprised "Group 1" tumors with a combined median survival time of 87.4 Ϯ 6.9 days. Group 2 tumors, generated by six other mutants, showed combined median survival of 124.1 Ϯ 6.1 days, and Group 3 tumors, generated by the R582W and ⌬(36 -53) mutants and WT generated tumors with combined median survivals of 210.2 Ϯ 17.7 days. Notably, S33A tumors were associated with significantly shorter median survival than S33Y tumors, and S45A tumors were associated with shorter median survival than S45F and S45P tumors ( Fig. 1B and Fig. S1B). The tumorigenicity of WT ␤-catenin was unsurprising and consistent with its being the primary oncogenic driver in FAP, in most sporadic CRCs, and in HBs harboring APC, Axin1, or Axin2 mutations (2,3,10,12,18,22,27,32,33,36). These initial results indicated that tumor growth rates and survival were dictated by the identities of the ␤-catenin mutants, that mutants not previously associated with HBs could nonetheless promote liver tumorigenesis, and that, in some cases, alternate missense mutations of the same amino acid could initiate tumors with quite different properties. We next examined ␤-catenin protein levels in representative tumors from each group. As seen in Fig. 1C and Fig. S2A, considerable variation was observed, with two of the slowest growing groups, namely ⌬(36 -53) and WT (Group 3, Fig. S1A), expressing the lowest levels. On the other hand, S33A and K335I tumors also expressed low ␤-catenin levels but were among the fastest growing (Group 1, Fig. 1B and Fig. S1A). Finally, the slowgrowing R582W tumors expressed relatively high levels of mutant protein (Fig. 1C). Thus, tumor growth rates and survival did not necessarily correlate with ␤-catenin protein levels.
Finally, to assess the degree to which each FLAG-tagged mutant ␤-catenin was expressed relative to endogenous murine ␤-catenin and to ascertain the level of expression of Myctagged ⌬(90) (23), we performed immunoblotting on total tumor or normal liver lysates using an anti-␤-catenin antibody. These results showed that all mutant ␤-catenin proteins were expressed at levels that were about equivalent to or modestly higher than that of the endogenous murine protein (Fig. S2B). ⌬ (90), which was the only mutant that could be clearly distinguished from its endogenous counterpart as a result of its large deletion, was the most highly expressed.
RNA-Seq data from a subset of nine tumor groups plus livers (see below) were next used to determine whether protein and transcript levels correlated. We quantified both human and endogenous murine ␤-catenin transcripts from a series of five tumors per cohort (Fig. 1, D and E). Whereas endogenous murine ␤-catenin transcripts showed little variation among the different groups, human transcripts showed a considerable amount but not always in ways that correlated with protein levels. For example, the biggest transcript disparities were seen in S33A and WT tumors, which expressed comparably low protein levels. Conversely, K335I and ⌬(20 -32) tumors showed significant differences in human ␤-catenin protein levels but not in transcript levels. Variations in ␤-catenin protein thus appear to be due to differential expression and/or stability of both transcripts and protein. The fact that every mutant transcript was expressed at levels significantly lower than those of WT transcripts suggests that many if not all mutant transcripts are unstable. Many patient-derived mutations also appear to exert opposite effects on transcript and protein levels.
To determine whether any of the above transcript and protein expression variability could be attributed to SB-␤-catenin vector copy number differences, we quantified these in representative HBs from each tumor group using a TaqMan-based assay. On average, intragroup and intergroup variation seldom exceeded 2-3-fold (Fig. 1, F and G). Exceptions to this were seen with S33A, ⌬(36 -53) and ⌬(90) tumors, which on average contained comparable and significantly fewer copy numbers than other tumor groups, and WT tumors, which contained ϳ50% more. However, these variations bore little relationship to RNA or protein levels. For example, S33A and WT tumors showed the greatest differences in integrated vector copy number (ϳ6-fold) and 8 -10-fold differences in ␤-catenin transcript levels. However, they expressed comparably low levels of ␤-catenin protein and grew at vastly different rates (Group 1 and 2, respectively). Thus, whereas some differences in integrated plasmid copy number were observed within and among each of the tumor groups, most of the variability in expression appeared to be due to differences in transcript and/or protein levels.
Subcellular fractionation studies showed the nuclear proportion of ␤-catenin mutants to range from a high of 95% in ⌬(45-58) tumors to a low of 39% in S37A tumors (Fig. 1, H-J). The nuclear compartmentalization of WT ␤-catenin was also high (82%) despite its overall low-level expression ( Fig. 1C and Fig.  S2). This was surprising, given the expectation that WT ␤-catenin would be largely, if not exclusively, confined to the cytoplasm by virtue of its efficient association with the APC  Table S1). The mutations N335I, N387K, and R582W are within the armadillo-repeat region that interacts with known partner proteins such as Tcf4, E-cadherin, and APC (45). B, Kaplan-Meier survival curves of mice with HBs driven by the indicated ␤-catenin mutants. Three tumor groups were identified based on differences in median survival (Fig. S1A). C, immunoblots of total ␤-catenin levels in whole-tumor lysates. See Fig. S2 for additional sets of tumors. Note that, because ⌬(90) is Myc epitope-tagged rather than His 6 -tagged, it is not detected here. D, total transcript levels for human and murine ␤-catenin. Bars, mean levels from five individual livers or tumors Ϯ S.E. (error bars). E, p values of each pair-wise comparison from the graph depicted in D. F, relative copy number of SB-␤-catenin vectors stably integrated into tumor DNAs from each of the indicated tumor groups. Copy numbers were determined using a TaqMan-based assay and were normalized to the nuclear gene encoding apolipoprotein B. G, p values for all pair-wise comparisons shown in F. H, nuclear:cytoplasmic distribution of ␤-catenin mutants and WT ␤-catenin in tumors. Four tumors (T1-T4) were evaluated for each mutant. Except for ⌬ (90), which was detected with an anti-Myc antibody, all ␤-catenins were detected with an anti-His 6 tag antibody. C, cytoplasmic fraction; N, nuclear fraction. I, quantification of the proportion of ␤-catenin from F that is nuclear. 4 -6 tumors from each group were used to obtain the results. J, p values for all pair-wise comparisons shown in I. K, WT ␤-catenin localizes mostly to the cytoplasm in rapidly growing ⌬(90) tumors. Mice were co-injected with equivalent amounts of SB vectors encoding Myc epitope-tagged ⌬(90) ␤-catenin and His 6 epitope-tagged WT ␤-catenin. The resulting tumors, all of which appeared rapidly (Group 1) (B), were fractionated as described in H.

␤-Catenin mutations and hepatoblastoma phenotypes
complex and proteasome. It suggested that tumorigenesis in response to WT ␤-catenin might require the selective outgrowth of a subpopulation of transfected hepatocytes with the highest nuclear ␤-catenin content to compensate for its low overall expression. This was consistent with the observation that, after appearing following their initial lag phase, these tumors grow at rates comparable with those of Group 1 tumors (not shown). To test this directly, we co-injected several mice with SB vectors encoding Myc-tagged ⌬(90) and His 6 -tagged WT. As expected, these tumors developed rapidly, indicating that they were primarily ⌬(90)-driven. Moreover, the majority of WT ␤-catenin was now observed to be cytoplasmic (Fig. 1K). This supported the idea that, in the absence of long-term selective pressure, the majority of WT ␤-catenin remains confined to the cytoplasm.
Among the different tumor groups, the crowded fetal (CF) pattern was the most commonly observed histopathologic subtype (8), although numerous exceptions were noted ( Fig. 2 and Table 1). These included a combination of CF, macrotrabecu-  Table 1.
a CF, crowded fetal; MT, macrotrabecular; BL, blastemal; PF, pleomorphic fetal; HCC, hepatocellular carcinoma. In cases where multiple histologic subtypes were observed, the most prominent one is listed first.

␤-Catenin mutants generate tumors with distinct metabolic features
⌬(90) HBs as well as HCCs generated by c-Myc overexpression undergo profound metabolic reprogramming involving higher rates of glycolysis and Warburg-type respiration in association with increased pyruvate dehydrogenase complex (PDC) activity (20,24,25,46,50). Concurrently, oxidative phosphorylation (Oxphos), particularly that utilizing fatty acids as an energy source, is significantly reduced. This fatty acid oxidation (FAO)-to-glycolysis switch occurs in the presence of and correlates with a pronounced decrease in mitochondrial DNA (mtDNA) content that presumably reflects the reduction in mitochondrial mass seen in numerous other cancers (24,25,(49)(50)(51).
In control livers, mitochondrial Complex I was responsible for ϳ30 -40% of electron transport chain (ETC) activity, whereas Complex II was responsible for the remaining ϳ60 -70% (Fig. 3, A-D) (24,25,46,50). As previously reported, Complex II activity was diminished in ⌬(90) tumors (20,24). Nearly all tumor groups had signature patterns of Complex I and Complex II activities, and in some cases, these appeared to be independent of one another. For example, the Complex I activities of S45A, S45F, and S45P tumors were higher than those of the ⌬(20 -32) tumors, which in turn were similar to livers (Fig. 3, A and B). In contrast, Complex II activities of ⌬(20 -32) tumors were lower than those of S45A, S45F, and S45P tumors (Fig. 3, C and D). Importantly, the Complex II activities for S33A and S45F tumors were also similar to those of livers. These results indicate not only that ␤-catenin mutants differentially affect Complex II activity but that its down-regulation is context-dependent and does not invariably accompany transformation.
The degree to which tumor groups down-regulated FAO and up-regulated PDC activity was also highly variable. For example, S33Y and S37A tumors down-regulated FAO to a lesser degree relative to ⌬(90) tumors, whereas T41A, S45A, and WT tumors did so to a greater degree (Fig. 3, E and F). PDC activity was similarly highly mutant-dependent, as evidenced by both S45A and R582W tumors having activities similar to those of livers, whereas ⌬(36 -53) and S33A tumors had activities 2.5-2.8 times higher, respectively (Fig. 3, G and H). These results indicate that different metabolic activities are controlled to different degrees, and in some cases independently, by the underlying ␤-catenin mutant. They also demonstrate that, as was true for FAO and Complex II activity, the up-regulation of PDC activity is not a prerequisite for transformation. mtDNA DNA content was also evaluated using a TaqManbased approach with PCR primers that amplified the D-loop region of the mitochondrial genome (24,50). This confirmed the previously observed ϳ80% reduction in total mtDNA content among all ⌬(90) tumors (Fig. 3, I and J) (24). Although all other tumor cohorts were associated with a reduced mtDNA content relative to livers, the magnitude of this loss was significantly less pronounced in S45F tumors.
Finally, we transcriptionally profiled mitochondrial ribosomal protein subunit transcripts. In keeping with the overall loss of mtDNA content in all tumors, these were also reduced, although the patterns of these changes were again highly ␤-catenin mutant-specific (Fig. 3K).

␤-Catenin mutants differ in their transactivation potential
␤-Catenin mutants might vary in their transactivation potential for various nonmutually exclusive reasons, including differences in their expression levels, the extent of their nuclear localization, and their ability to engage Tcf family members and other components of the transcriptional machinery. To examine this, we employed a vector (7TGC) that expresses enhanced GFP (eGFP) driven by a Tcf4/␤-catenin-responsive promoter (52). We co-transfected this with each SB-␤-catenin vector into HEK293 cells and quantified eGFP expression 48 h later. Significant differences were observed in the degree to which each ␤-catenin mutant increased eGFP expression (Fig. 4, A and B). For example, ⌬(36 -53) was the most potent activator of eGFP, whereas R582W and WT were the weakest. In parallel studies, we repeated these same transfections but included the SB-YAP S127A vector. YAP S127A exerted little effect on the reporter except in the case of the S45A mutation, where it modestly but significantly suppressed eGFP expression (Fig. 4, A and C).

␤-Catenin mutants drive distinct HB transcriptional profiles
We hypothesized that the above-reported differences in tumor phenotypes could be explained by variations in their global transcript expression profiles, as already seen for the response of 7TGC (Fig. 4). We further surmised that the overall transcriptional profiles of these tumors would resemble one another more closely despite their variable histologic features ( Fig. 2) (18,24,25,48). We chose tumors generated by nine of the ␤-catenin mutants and subjected five of each along with five normal control livers to RNA-Seq. The former included the Group 1 tumors S33A, S45A, K335I, N387K, ⌬(20 -32), and ⌬(90) and the Group 2 and 3 tumors S33Y, ⌬(36 -53), and WT. K335I and N387K were chosen because their mutations resided within the armadillo repeat region, which is far removed from the main phosphorylation domain located between residues 33 and 45 and whose coordinated multisite phosphorylation is critical for cytoplasmic retention by the APC complex and its subsequent degradation by the proteasome (27)(28)(29) (Fig. 1A). WT tumors were also selected because of their pronounced HCC-like histology ( Fig. 2 and Table 1). This allowed a comparison of their transcriptomes with those from tumors generated by other ␤-catenin mutants and HCCs induced by ectopic Myc overexpression (50).
Relative to livers, tumor cohorts showed a mean of 5646 significant gene expression differences (range ϭ 4222-6820, ␤-Catenin mutations and hepatoblastoma phenotypes q Ͻ 0.05) ( Fig. 5A and Table S2). As hypothesized, the differences were much less among individual tumor cohorts, ranging from as many as 1515 in the case of S33A versus S33Y tumors to as few as 1 in the cases of S33A

␤-Catenin mutations and hepatoblastoma phenotypes
5B) and Table S2). Interestingly, S33Y showed the second fewest differences relative to normal liver (4766) but the largest number of differences relative to the eight other tumor groups (mean ϭ 1071) (Table S2). When S33Y was excluded from the comparison, the tumor group showing the smallest number of differences relative to others was ⌬(90) (average ϭ 48, p Ͻ 10 Ϫ9 ).
Despite the pronounced HCC-like histologic appearance of WT tumors and the variable but significant HCC-like content of several other tumor groups ( Fig. 2 and Table 1), all tumors induced as a result of WT or mutant ␤-catenin overexpression resembled one another at the whole-transcriptome level much more closely than they resembled HCCs induced by Myc overexpression (50) (Fig. 5C and Fig. S3).
Ingenuity pathway analysis (Qiagen) was used to identify the top pathways whose component transcripts were most deregulated in at least one tumor cohort relative to that of liver. In keeping with the biochemical and metabolic abnormalities OCRs were determined on tissue homogenates following the addition of pyruvate, malate, glutamate, and ADP (20,24,25,46,49,50). B, significant pair-wise comparisons between individual tumors groups based on the results shown in A. C, Complex II activity. Following the addition of succinate to achieve the maximal combined activities of Complexes I ϩ II, rotenone was added to inhibit Complex I. The residual activity (i.e. that of Complex II) is shown here. D, significant pair-wise comparisons between individual tumor groups based on the results shown in C. E, FAO activity in the indicated tissues. F, significant pair-wise comparisons between individual tumors groups based on the results shown in E. G, PDH activity in the indicated tissues. H, significant pair-wise comparisons between individual tumors groups based on the results shown in G. I, mtDNA quantification. A TaqMan assay was used to amplify a region of the mitochondrial D-loop region as described previously (24,50). Results were normalized to an amplicon from the nuclear gene encoding apolipoprotein B. J, significant pair-wise comparisons between individual tumors groups based on the results shown in I. K, transcript levels of mitochondrial ribosomal protein subunits. The indicated liver and tumor tissues (five per group) were subjected to RNA-Seq analysis. Levels of 78 transcripts encoding each of the mitochondrial ribosomal protein subunits are expressed as mean values and depicted as Z-scores. Error bars, S.E. ␤-Catenin mutations and hepatoblastoma phenotypes characterized above, seven of the top 15 pathways identified by this approach involved mitochondria-associated processes, including Oxphos, the TCA cycle, branched chain amino acid catabolism, and ␤-FAO, and all were predicted to be downregulated in tumors (Fig. S4). The degree to which each component pathway's transcripts were deregulated (i.e. the Z-score) was also found to be highly characteristic of each ␤-catenin mutant.
Gougelet et al. (53) have identified 613 expressed genes whose proximal promoters directly bound ␤-catenin/Tcf4 complexes in primary murine hepatocytes. The differentially expressed genes in each of our tumor groups overlapped with 44.4 -61.0% of these (mean ϭ 324 (53%); Fig. 5D). Thus, on average, the majority of previously identified ␤-catenin/Tcf4 target genes in liver were dysregulated in HBs by different ␤-catenin mutations, with each cohort possessing a distinct expression profile.
Among the most highly up-regulated transcripts in ⌬(90) tumors is Myc, whose hundreds of target genes include many that encode metabolic enzymes and their regulators (23-25, 54, 55). We recently underscored Myc's functional importance for HB by showing that ⌬(90) tumor growth is markedly impaired in myc Ϫ/Ϫ livers (24,25). We speculated that, by altering pathways involving glycolysis, FAO, Oxphos, and ribosomal biogenesis, Myc balances the growth of tumors with their metabolic needs. Large differences in Myc expression were seen among different tumor groups, although nearly all expressed substantially more Myc than livers (Fig. 6A). All three slowly growing, Group 3 tumors were among those with the lowest Myc levels, and five of the six Group 1 tumors were among the highest Myc expressers. For the remaining tumors, however, this correlation was imperfect. T41A and ⌬(20 -32) tumors were among the fastest growing despite expressing quite low levels of Myc (Group 1, Fig. S1A). Four of the six Group 2 tumors expressed low levels of Myc, but two (S45P and ⌬(45-58)) were among the highest expressers. Myc levels also did not always correlate with the ␤-catenin mutants' activation of the 7TGC promoter (Fig.  4, A and B). For example, ⌬(36 -53) generated tumors with low Myc levels despite being the most potent 7TGC activator, whereas the reverse was true for R582W and ⌬(90). These studies suggest that ␤-catenin mutants differentially up-regulate Myc and that some do so independently of their effects on other genes. In keeping with this idea, L-Myc expression was high in control livers, down-regulated in most but not all tumors, and did not necessarily correlate with those of Myc (Fig. 6A).
We next assessed the expression of 634 direct Myc target genes in the liver using data from the Ingenuity Pathway Analysis (IPA, Qiagen, Inc.) Knowledge Base (25,50). Between 204 and 277 of these (mean ϭ 243 (38.3%)) were deregulated in each of the tumor cohorts ( Fig. 6B and Table S3). In keeping with the fact that liver expresses very low levels of Myc (Fig. 6A), most of these tumor transcripts were expressed in the opposite direction. Interestingly, each HB group expressed its own signature profile of these transcripts that may be in part reflect competition between Myc and L-Myc (Table S3). Each tumor group's expression of these individual target genes as a group was unique and in keeping with our results regarding Myc protein expression and 7TGC responses discussed above. This was again consonant with the idea that each ␤-catenin mutation affected each target gene's response in a specific and largely individualized manner that accounts for this distinct transcript patterning.
IPA was again used to identify families of functionally related transcripts that differed among the nine tumor groups. Over half of the top pathways so identified were involved in aspects of metabolism involving activities such as Oxphos, lipid and xenobiotic metabolism, and FAO. Consistent with the overall reduction of Complex II activity and mtDNA content (Fig. 3), transcripts encoding most TCA cycle enzymes, ETC subunits, and FAO-related enzymes were down-regulated, although the patterns of these changes were again highly mutant-specific (Fig. 6, C-E). The vast majority of transcripts encoding enzymes involved in -FAO, a nonenergy producing and largely cytoplasmic process (56), were also down-regulated (Fig. 6F).
Consistent with the general paucity of HB driver mutations other than those involving ␤-catenin (13-15, 18, 47), we failed to identify de novo mutations in any transcripts expressed by the above 45 tumors (not shown). However, we have previously demonstrated that expression of the p19 ARF tumor suppressor is markedly attenuated in ⌬(90) tumors without being mutated and that tumorigenesis was strongly abrogated by SB-mediated ectopic p19 ARF expression (46). To determine whether p19 ARF loss was a general finding, we evaluated its levels in representative samples from each tumor cohort and verified its downregulation in all ⌬(90) tumors relative to that of liver (Fig. 6G). Similar findings were made in most other tumor groups but not all. For example, S45P, ⌬(20 -32) and ⌬(48 -58) tumors (all Group 1 or Group 2, Fig. 1B) showed, at best, only modest suppression of p19 ARF . Thus, whereas p19 ARF down-regulation occurs frequently and reproducibly in certain HB cohorts, it is Each column represents the mean transcript expression level of five tissue samples from that cohort. Gene expression differences were identified from RNA-Seq data using CLC Genomic Workbench (Qiagen), DeSeq2, and EdgeR, and only those transcripts identified by all three methods as being differentially expressed are included. The number above each tumor cohort indicates the differentially expressed transcript differences between that group and normal liver (also see Table S3). The mean number of differences among all groups was 5646. B, gene expression differences among each of the individual cohorts. Numbers above each tumor group indicate the number of gene expression differences compared with the ⌬(90) cohort. See Table S2 for a comparison of gene expression differences among all cohort comparisons. C, gene expression differences between each cohort of ␤-catenin-induced tumors and HCCs. The latter were induced by induction of a doxycycline-regulated human Myc transgene as described previously (50). Numbers below each cohort indicate the number of gene expression differences relative to the HCC group. This along with principle component analysis (PCA; supporting Fig. S3) indicated that, despite their histologic differences, all ␤-catenin-induced tumors were much more HB-like than HCC-like at the molecular level and more closely resembled liver than they did HCC (Fig. S3). D, HBs deregulate a majority of genes identified as direct and expressed targets of ␤-catenin/Tcf4. Transcripts in the "liver" cohort include the entire 613 repertoire of expressed genes from normal liver whose proximal promoters were bound by ␤-catenin/Tcf4 (30). Numbers below the tumor cohorts indicate the number of transcripts that were deregulated relative to those in the liver. See Table S3 for a pairwise comparison of Myc target gene differences among all cohorts.

␤-Catenin mutations and hepatoblastoma phenotypes
not a universal finding. Together with the results presented in Fig. 6A, this suggests that the pathogenesis of HBs by different ␤-catenin mutants does not invariably involve the inactivation of the same tumor suppressor pathways or the activation of the same oncogene pathways.
Cairo et al. (48) previously described a 16-transcript signature that allowed classification of tumors into two groups, C1 and C2, with favorable and unfavorable outcomes, respectively. The latter group is associated with adverse clinical features, such as higher stage at the time of diagnosis and evidence for vascular invasion. A subsequent study of 24 additional HBs by Hooks et al. (19) subdivided the C2 group into two subgroups based on the expression of four transcripts. C1 HBs displayed a predominantly fetal histology, whereas C2 tumors possessed more immature embryonal features, more mitoses, and more intense nuclear ␤-catenin staining. Murine embryonal liver stem cells demonstrated a C1-like profile prior to their ectopic overexpression of ␤-catenin, Myc, or N-Myc and a C2-like profile afterward. C2 patterns were seen in the human HB cell lines Huh6 and HepG2, which reverted to C1 when either Myc or ␤-catenin was silenced (57). When combined with an additional collection of 85 HBs, 79% of tumors contained an assortment of ␤-catenin mutations.
We combined our RNA-Seq data with those of Hooks et al. (19) and assessed each of our tumor groups for the expression of the murine orthologs of the signature 16 transcript collection identified by Cairo et al. (48). We used the CLC Genomics Workbench 12.0 program to cluster these in an unbiased manner into groups with related expression patterns. As seen in Fig.  6H, this approach indicated that murine HBs appeared most closely related to the C1 human molecular subtype, irrespective of their underlying ␤-catenin mutation.

Discussion
Survival in HB is associated with certain clinical and laboratory features including age, stage, histologic subtype, and serum ␣-fetoprotein levels (2, 7-9, 15, 58 -60). Newer transcriptomebased classifications have identified the C1/C2 groups and the C1/C2A/C2B groups as having distinct long-term survival differences based on 16-gene or four-gene signatures, respectively (19,48). Although 80% or more of HBs bear highly heterogeneous ␤-catenin mutations (3, 16 -18, 21-23, 26, 60), these and the tumor's rarity and extensive histologic variability have made it impossible to correlate specific mutations with survival or other features. Early surgical and/or chemotherapeutic intervention further preclude any assessment of the natural history of tumors associated with specific ␤-catenin mutations.
Our results in mice indicate that several HB biological, biochemical, metabolic, and transcriptional behaviors are predetermined by the identity of the ␤-catenin mutation. It remains unclear whether our failure to detect metastases in any cases is the result of the inadvertent selection of mutants that favored the generation of only localized disease or because most endogenously arising murine tumors tend to have low metastatic propensity (61,62).
CTNNB1 gene mutations occur in numerous other human cancers (Table S1) (16 -18, 21, 22, 40 -43). Whereas some of these are shared with HB, others are not, thus raising the questions of whether they are all equally oncogenic in the liver and, if so, whether they differentially influence tumor behaviors. It seemed reasonable to surmise that in HB, where few if any other recurrent driver mutations exist (3,13,14), different ␤-catenin mutants influence not only growth rates but additional phenotypes. Examples of other oncoproteins exhibiting similar behaviors exist, including non-small-cell lung cancer, where survival and response to erlotinib and gefitinib correlate with missense or deletion mutations, respectively, in the epidermal growth factor receptor (63,64). In CRCs, K-Ras codon 12 missense mutations are associated with inferior survival relative to codon 13 or codon 61 mutations (65). In other cancers, PIK3CA and MAP2K(MEK1) mutations correlate with chemo-sensitivity and survival in vitro, in vivo, and in patients (66 -68).
With regard to ␤-catenin specifically, previous transient transfection studies demonstrated differences in the abilities of nine ␤-catenin mutants, including S45F and S45P, to activate a synthetic reporter similar to 7TGC (Fig. 4) (69). These findings were limited, however, in that some of the ␤-catenin mutants were not patient-derived, nuclear:cytoplasmic distribution of the mutants was not correlated with total expression levels, and other target gene responses were not examined. Consistent with our results, a recent in vivo study showed that ⌬(90) and S45Y, in combination with c-Met, induced HCCs with similar growth rates and histologies but distinct gene expression profiles (70). Thus, our finding that each ␤-catenin mutant was associated with signature transcriptional profiles (Figs. 4 -6) is not without precedent, likely explains the distinct features of each tumor group and provides an overarching mechanistic explanation for the phenotypic differences observed among each tumor group.
Two of the fastest growing tumor cohorts, namely S45A and N387K (Group 1, Fig. S1A), were driven by mutants not previously associated with HBs. By comparison, the common HBassociated mutants S33Y and T41A drove more slowly growing  (53). Transcripts are arranged from highest to lowest expression based on levels in liver. Numbers at the bottom indicate the number of significant expression differences versus liver. Table S3 shows the actual number of differentially regulated transcripts among all pairwise comparisons. C, transcripts encoding TCA cycle enzymes. Table S4 shows the number of similarly regulated transcripts between any two tissues relative to those in liver. D, transcripts encoding ETC and Complex V subunits. Table S5 shows the number of transcript differences among all pairwise comparisons. E, differential regulation of transcripts encoding mitochondrial FAO-related transcripts. Table S6 shows the number of transcript differences among all pairwise comparisons. F, differential regulation of (endosomal) FAO (␤-oxidation)-related transcripts. Table S7 shows the number of transcript differences among all pairwise comparisons. G, expression of p19 ARF in liver and tumor samples. H, murine HBs most closely resemble the C1 molecular subtype of human HBs. The CLC Genomics Workbench hierarchical clustering program was used to compare and organize the murine orthologs of the 16-transcript set of Cairo et al. (48) with the RNA-Seq data from Hooks et al. (19). Human tumor results are from individual tumors, whereas murine results represent the mean expression value of five tumors for each group as done for other panels in the figure. Human C1 tumors are indicated by the green bars at the bottom of the heat map, human C2A and B tumors are indicated by the yellow bars, and murine tumors are indicated by magenta bars. The scale bar here indicates the -fold change rather than the Z-score. Group 2 tumors (Fig. 1B and Fig. 1) (35). The implication of these findings is that mutations previously identified in non-HB tumors may be as tumorigenic as "classical" HB-associated mutants. The recurrent nature of certain mutants therefore implies neither that they are the more efficient oncogenes nor that they drive particularly aggressive neoplasms. Rather, their frequency may reflect greater context-dependent mutagenic susceptibility of certain codons. The identities of these need not be fixed and might well vary as a function of age and/or context. Alternatively, some mutations detected thus far only in adult tumors may yet be identified in HBs as testing for these becomes more commonplace.

␤-Catenin mutations and hepatoblastoma phenotypes
Like growth, HB histology was also ␤-catenin mutant-dependent. For example, both R582W and WT HBs were almost exclusively HCC-like, as were to a lesser extent T41A and K335I tumors ( Table 1 and Fig. 2). Yet no other factors appeared to be predictive of these histologies. The lack of greater histopathologic diversity among our tumor groups more akin to that of the seven HB subtypes (8) suggests either that different ␤-catenin mutants and/or other genes are required for such differences to be fully manifested. These could include TP53, TERT, and NFE2L2, which are recurrently mutated in a small minority of HBs (14,15,18,32). Given the inability to detect such cooperating mutations in the tumors studied here using the CLC Genomics Workbench Basic Variant Detection algorithm, it is tempting to speculate that the histologic variability we did encounter is a result of the inherent transcriptional differences of the ␤-catenin mutants and perhaps epigenetic differences that account for the altered expression of important cooperating genes, such as Myc, L-Myc, and p19 ARF (Fig. 6, A and G). It remains unclear to what extent different tumor histologies are interconvertible and whether they are influenced by factors such as the metabolic state or microenvironment, which can vary over quite small intratumoral distances (Table 1 and Fig. 2) (2,5,8,18).
A notable feature of WT and R582W tumors was their undeniable retention of HB-like transcriptomes despite histologically resembling HCCs (Table 1, Figs. 2 and 5 (A and C), and Fig.  S3). This contrasts with the well-known observation that FAPassociated liver tumors, which express WT ␤-catenin, are always HBs (10). The discrepancy may be at least partly explained by the long latency period of WT murine tumors, which arise well into adulthood. HBs are believed to originate from short-lived primary hepatoblasts or fetal-like multipotent progenitors, thus explaining the strict age dependence even in individuals with FAP (3). This population might be depleted and/or altered in mice by the time WT and R582W tumors appear, at which point they more closely resemble HCCs histologically while retaining HB transcriptomic profiles. Early liver progenitors and fetal-like hepatocytes may also degrade WT ␤-catenin more efficiently within the cytoplasmic APC complex, thus delaying the onset of its stabilization and nuclear translocation (71). This was suggested by our finding that ⌬(90) tumors showed more pronounced cytoplasmic retention of coexpressed WT ␤-catenin than did WT tumors (Fig. 1, compare  H with K). Slowly growing R582W tumors appear even more HCC-like than WT tumors (Table 1 and Fig. 2), and four of the five remaining cohorts with less conspicuous HCC-like histol-ogy were also categorized as belonging to the slower-growing Group 2 or 3 (Table 1). A nonmutually exclusive explanation is that YAP S127A , while incapable of HB induction alone (23), can nevertheless induce HCC in older mice and that its contribution becomes increasingly important with age, whereas WT ␤-catenin continues to dominate the transcriptional output (71). There may exist a specific age-dependence for each ␤-catenin mutant such that the histology of early onset tumors is more HB-like and more ␤-catenindependent, whereas later onset tumor are more HCC-like and YAP S127A -dependent. A similar "window" of susceptibility may also apply to the defects that deregulate other facets of the Hippo pathway (72,73). That ⌬(90) generates HBs or HCCs based on whether it is co-expressed with YAP S127A or c-Met lends further credence to the argument that ␤-catenindriven tumorigenesis can be promiscuous and influenced by other co-expressed oncoproteins as well as age-and tissue-specific factors (70). We did not examine the S45Y mutant used by Qiao et al. (70) but think it likely that it would also drive HB tumorigenesis in association with YAP S127A , given the robust tumorigenicity of S45A, S45F, and S45P (Fig. 1B).
Both absolute ␤-catenin level and nuclear localization are of prognostic significance in some HB studies, although the degrees to which these are mutually and mutationally dependent have not been determined (3,58,74). Missense and in-del mutations can profoundly affect mRNA stability by altering folding, splicing, recognition by micro-RNAs, and translational efficiency (75)(76)(77). Independently, steady-state protein levels can be affected by mutation-dependent effects on folding, protein-protein associations, and post-translational modifications (78 -80). This latter is particularly germane in the case of ␤-catenin mutations, which often involve key sites of APC complex-mediated phosphorylation that determine ␤-catenin's stability and nuclear translocation (27)(28)(29). However, our results suggest that the consequences of such mutations are more complex and nuanced than might have been previously surmised. For example, ⌬(20 -32) and ⌬(45-58) possess comparably sized deletions, with the former affecting none of the known APC phosphorylation sites and the latter affecting only Ser-45. Accordingly, ⌬(20 -32) would have been expected to be expressed at a lower level than ⌬ (45)(46)(47)(48)(49)(50)(51)(52)(53)(54)(55)(56)(57)(58), to be cytoplasmically restricted and to be less pathogenic. Yet, whereas it was indeed more cytoplasmically confined, its total expression was equivalent to ⌬(45-58)'s while generating faster-growing tumors (Fig. 1, A, C, and H-J). Other examples were seen with S33A and S45F. Although each bore only a single phosphorylation site mutation, S33A was expressed at a lower level despite being more nucleus-localized. Yet it also generated faster-growing tumors.
␤-Catenin mRNA and protein steady-state levels also did not often correlate. For example, WT transcripts were ϳ10-fold more abundant than those of S33A, although tumor protein levels were similar (Fig. 1, C and D). Conversely, S45A and K335I transcript levels were equivalent, although S45A protein was more highly expressed. Low ␤-catenin protein levels were associated with three of the four slowest growing tumor groups (R582W, ⌬(36 -53), and WT), although this was not always so, as evidenced by the rapid growth of K335I tumors (Group 1, Fig. 1B and Fig. S1A). All of these proteins, however, efficiently localized to the nucleus (Fig. 1, H-J), suggesting that this could serve as a pro-oncogenic mechanism to compensate for overall low ␤-catenin expression. This was supported by our findings with WT ␤-catenin, which was virtually all nuclear in slowly growing tumors but mostly cytoplasmic when co-expressed with ⌬(90) in rapidly growing tumors (Fig. 1, H and K). The long latency of WT tumors may allow clones with higher expression of nucleus-localized ␤-catenin to be selected and become dominant. The nuclear content of ␤-catenin may thus vary over time. Nonetheless, our findings support the idea that the oncogenic potency of any particular ␤-catenin mutant is more influenced by the amount that is nucleus-localized than by the total amount that is expressed. The mutant-dependent nature of ␤-catenin expression could have significant implications for newer therapeutic approaches that abrogate its interaction with the transcriptional machinery (81,82). It remains to be determined whether these assumptions will hold for non-HB ␤-catenin-driven tumors. More precise determinations of mutant ␤-catenin transcript and protein half-lives will require carefully controlled conditions using metabolic inhibitors, such as actinomycin D and cycloheximide, or transient affinity labeltagging methodologies.

␤-Catenin mutations and hepatoblastoma phenotypes
In tumor groups such as S33A, S45A, ⌬(20 -32), ⌬(45-58), and ⌬(90), the number of stably integrated SB-␤-catenin copy numbers was remarkably consistent among individual tumors, and in other cohorts, it seldom differed by more than 2-3-fold (Fig. 1F). Importantly, mice from each cohort were injected at different times and with different DNA preparations during tumor induction. The highly reproducible intragroup ␤-catenin copy number for the above-mentioned cohorts is testament to the fact that neither of these variables exerted significant influence over the final copy number. Interestingly, those tumor groups with the least copy number variation were all associated with rapidly growing tumors, whereas more slowly growing tumors, such as S33Y, S37A, R582W, and WT, showed greater variation.
Each ␤-catenin mutant also impacted metabolic pathways in distinct ways (Figs. 3 and 6 (C-F)). Mechanistically, this likely reflects differences in the degree to which each mutant regulates these pathways' component transcripts (Figs. 3K and 6 (B-F)). Independent of total transcript levels, the expression patterns of small numbers of transcripts in cancer-related pathways have recently been shown to better predict survival across more cancer types (46,49,(83)(84)(85). Such patterns may serve as surrogate reporters of a tumor's transcriptional environment and ultimately of its clinical behavior (49,83,84). Considering only the differences in ␤-catenin, Myc, and L-Myc protein levels among the different tumor groups (Figs. 1C and 6A and Fig. S2), we can envision a model in which ␤-catenin mutants differentially interact with Tcf/Lef family members, each with its own inherent DNA binding affinity and target gene subsets (86 -88). Myc-and L-Myc network members and their target transcripts would also be expected to possess distinct expression (24,25,46,49,83,84,89). Together, our findings suggest that the mechanisms underlying each ␤-catenin mutant's distinct effects on transformed hepatocyte behavior may involve nonmutually exclusive combinations of total protein levels, subcellular par-titioning, interaction with the transcriptional machinery, crosstalk with the Hippo-YAP pathway, the affinity for target genes, and the magnitude and patterns of their response.
In conclusion, we provide the first direct evidence showing that much of the diversity that distinguishes HBs can be attributed to the underlying identity of their underlying ␤-catenin mutations. Those previously identified only in non-HB tumors are also tumorigenic in the liver, thus suggesting that much of ␤-catenin's oncogenic capacity is neither mutant-nor tissuerestricted. Additional work will be necessary to determine whether these mutations can be used in meaningful ways to predict tumor behavior, response to and durability of chemotherapy, and long-term survival in patients.

Animals and HDTVI
HDTVIs of purified plasmid DNAs were performed as described previously (20,24,25,46). Injections, in a total volume of 2 ml in 0.9% sodium chloride, contained 10 g each of SB vectors encoding the indicated ␤-catenin mutant and the previously described YAP S127A and 2 g of a non-SB vector encoding SB translocase. Each injection was performed over 6 -9 s in 6 -7-week-old FVB mice (Jackson Laboratories, Bar Harbor, ME). To minimize variability in technique or plasmid DNA preparations, 4 -5 mice were injected at a time until a total of 12-15 animals/group had been injected. At least two different plasmid preparations were also used for the generation of each tumor, and all plasmids were prepared at the same time. Neither the timing of injections nor DNA preparations were found to exert significant effects on any tumor features. Mice were monitored three times weekly for evidence of tumor growth and daily after tumors became apparent by physical exam. Mice were euthanized when tumors attained a size of 2 cm in any dimension or when obvious signs of stress were noted. All injections, monitoring procedures, and other routine care and husbandry were previously approved by the University of Pittsburgh institutional animal care and use committee.

Plasmid DNAs
A WT human ␤-catenin cDNA, fused at its 3Ј-end with a His 6 epitope tag, was cloned into the previously described EF1a SB expression vector (20, 23, 24) using a Gateway TM cloning strategy (Thermo Fisher Scientific). Mutagenesis was performed with a QuikChange II site-directed mutagenesis kit (Agilent, Inc., Santa Clara, CA), and each mutant was verified directly by dideoxy-DNA sequencing. Plasmid DNAs were purified using Qiagen columns (Qiagen, Inc., Valencia, CA), quantified on a NanoDrop Microvolume Spectrophotometer (Thermo Fisher Scientific), and checked for overall integrity by gel electrophoresis. The mutations were selected based upon their location and whether they had been previously identified only in HBs, only in non-HB tumors, or in both HB and non-HB tumors (Table S1) (http://ongene.bioinfo-minzhao.org/ gene_mutation.cgi and https://cancer.sanger.ac.uk/cosmic/ gene/analysis?lnϭCTNNB1). 7

␤-Catenin mutations and hepatoblastoma phenotypes TaqMan assays
Relative SB-␤-catenin plasmid copy numbers in tumors were quantified using a TaqMan-based assay. Tumor tissues from several different regions of the liver were first combined so as to minimize the possibility that any differences in plasmid copy number were due to clonal variability. Briefly, HB DNAs were extracted using a DNeasy Blood and Tissue Kit according to the directions of the supplier (Qiagen). Two sets of primers were used, the first of which was designed to amplify human ␤-catenin sequences encoded by the SB vector. The "forward" PCR primer was designed to be homologous to the region encoded by nt 1237-1331 of ␤-catenin (GenBank TM no. NM_001904.4) and consisted of the sequence 5Ј-CAGAGTG-CTGAAGGTGCTATC-3Ј. The "reverse" primer was taken from a region of the SB vector just downstream of the ␤-catenin insertion site: 5Ј-GATCTGTCAGGTGAAGTCCTAAAG-3Ј. Amplification of the 95-bp PCR product was detected in real time with the TaqMan probe: 5Ј-6-FAM/CCGGCTA-TTGTAGAAGCTGGTGGAA/3Ј Iowa Black FQ/-3Ј.
The second PCR primer set was used to amplify a segment of the murine nuclear apolipoprotein B gene. The "forward" PCR primer was designed to be homologous to the region encoded by nt 8010487-8010504 (GenBank TM no. NC_000078.6); sequence, 5Ј-CACGTGGGCTCCAGCATT-3Ј). The "reverse" primer was designed to be homologous to nt 8010539 -8010560 (5Ј-TCACCAGTCATTTCTGCCTTTG-3Ј). Amplification of the 74-bp product was detected with the TaqMan probe: 5Ј-/Cy5-CCAATGGTCGGGCACTGCTCAA-Black Hole Quencher 2-3Ј.
PCR was performed using Kappa Probe Fast qPCR Kit under the following conditions: 95°C for 10 s; 40 cycles at 95°C for 15 s, and 60°C for 1 min. The human ␤-catenin copy number in each sample was normalized to that of the murine nuclear apolipoprotein B gene PCR product. All assays were performed on a CFX96 Touch TM real-time PCR detection system (Bio-Rad). All oligonucleotide primers and TaqMan probes were synthesized by IDT, Inc. (Coralville, IA).
TaqMan assays that amplified a region of the mitochondrial DNA "D-loop" region were performed on the above DNAs as described previously (24,49).

Cell fractionation and immunoblotting
At the time of sacrifice, tumors or control livers were harvested, immediately placed on ice, sectioned into small fragments with a scalpel, and then snap-frozen in liquid nitrogen before being stored at Ϫ80°C. Nuclear and cytoplasmic fractionations were performed on small fragments of fresh tissues (ϳ200 mg) using a subcellular protein fractionation kit according to the directions of the supplier (Thermo Fisher Scientific). For whole-cell analyses, tissue lysates were prepared by placing the frozen samples in 0.5-1.0 ml of SDS-PAGE buffer containing a mixture of protease and phosphatase inhibitors as described previously (20,46,50), followed by rapid homogenization in a BBY24 bullet blender (WisBioMed, San Mateo, CA). SDS-PAGE, transfer to polyvinylidene difluoride membranes, and antibody probing of the membranes were performed as described previously (20,46,50). Antibodies used included those specific for the His 6 tag (BD Biosciences, catalog no. 552565), c-Myc (Cell Signaling (Danvers, MA), catalog no. D3N8F), L-Myc (Santa Cruz Biotechnology, Inc., catalog no. sc-790), ␤-catenin (Abcam, catalog no. 16051), GAPDH (Sigma, catalog no. 71.7), histone H3 (Cell Signaling, catalog no. 9715), and p19 ARF (Santa Cruz Biotechnology, catalog no. sc-22784). All antibody dilutions used were those recommended by the respective vendors. Secondary horseradish peroxidase-conjugated antibodies and chemiluminescent detection (SuperSignal West Pico chemiluminescent substrate kit (Thermo Fisher Scientific) were used and performed as described previously (20,25,50).

Metabolic studies
Oxygen consumption rates (OCRs) were performed on fresh, disrupted tissues immediately after harvesting, using an Oroboros Oxygraph 2k instrument (Oroboros Instruments, Inc., Innsbruck, Austria) as described previously (20,24,25,46,50). Baseline OCRs were determined in 2-ml samples of tissue (ϳ50 mg) in Mir05 buffer containing 10 M cytochrome c and again following the sequential addition of malate (2 mM), ADP (5 mM), pyruvate (5 mM final concentration), and glutamate (10 mM). After achieving plateau OCRs, which assessed the maximal capacity of Complex I, succinate was added (10 mM final concentration) to determine the additional contribution provided by Complex II. The Complex I inhibitor rotenone was then added (0.5 M final concentration), thus allowing the proportional contribution of Complexes I and II to be independently verified. All activities were normalized to total protein.
Pyruvate dehydrogenase (PDH) complex activity was determined on 50 -100 mg of fresh, minced tissue as described previously (20,24,25,46,50). Briefly, liver and tumor tissues were placed into ice-cold glucose-, pyruvate-, and glutamine-free Dulbecco's modified minimal essential medium and homogenized by several passages through an 18-gauge needle. The conversion of [1-14 C]pyruvate (PerkinElmer Life Sciences) to 14 CO 2 was then measured for 40 min at 37°C and quantified by scintillation counting.
␤-FAO was also measured as described previously (20,24,25,46,50). Briefly, ϳ50 mg of fresh tissue was minced in 350 l of ice-cold SET buffer and disrupted by five strokes of a Potter-Elvehjem homogenizer. 5 l of the homogenate was then added to 195 l of FAO reaction buffer containing 0.5 Ci per ml of 3 H-labeled palmitate conjugated to BSA (PerkinElmer Life Sciences). Following incubation for 2 h at 37 ºC, the reaction was terminated by adding 40 l of 1 M KOH and then further incubated at 60°C for 1 h to hydrolyze newly synthesized acyl-carnitine esters. 40 l of 4 M perchloric acid was added, followed by an additional incubation on ice for 60 min and organic extraction. The water-soluble 3 H-labeled products were then quantified by scintillation counting and normalized to total protein content.

Cell culture and transient transfection
HEK293 cells were routinely propagated in Dulbecco's modified Eagle's essential medium containing 10% fetal bovine serum, 2 mM glutamine, and penicillin ϩ streptomycin. The day prior to transfection, the cells were trypsinized and seeded into 6-well tissue culture plates so as to be ϳ30% confluent at the ␤-Catenin mutations and hepatoblastoma phenotypes time of transfection. The cells were then transfected for 4 h using Lipofectamine 2000 according to the directions provided by the vendor (Thermo Fisher Scientific). Each transfection was performed in triplicate on several occasions and contained 2 g of the indicated SB-␤-catenin vectors and 1 g of the reporter vector 7TGC (52) (Addgene, Inc.). The latter is a lentiviral vector that contains a minimal promoter flanked by seven Tcf4/ Lef-binding sites that drive the expression of eGFP in response to Wnt/␤-catenin signaling (52). The vector also contains an mCherry expression cassette under the control of a neutral SV40 promoter to allow for normalization of the eGFP signal. For some transfections, 2 g of the SB-YAP S127A expression vector was included. Following transfection, the cells were propagated for an additional 48 h in standard growth medium, at which time they were harvested for flow cytometry. eGFP and mCherry signals were measured on a BD Biosciences LSR-Fortessa flow cytometer using the following conditions: eGFP-B1, 530/30 nm; Cherry-D2, 610/20 nm. All eGFP signals were normalized to those of mCherry.

Transcriptional profiling of HBs
RNA purification from five representative tumors of each group was performed using Qiagen RNAeasy columns (Qiagen) followed by DNase digestion. Sample integrity was measured using an Agilent 2100 Bioanalyzer (Agilent Technologies, Foster City, CA). All samples had RIN values of 8.5-10 prior to any further processing. Samples for sequencing were prepared using a New England Biolabs NEBNext Ultra Directional RNA Library Prep kit according to the directions provided by the vendor (New England Biolabs, Beverly, MA), and sequencing was performed on a NovaSeq 600 instrument (Illumina, Inc., San Diego, CA) by Novagene, Inc. (Sacramento, CA). Raw and processed original data were deposited in the National Center for Biotechnology Information (NCBI) Gene Expression database and are accessible through the Gene Expression Omnibus (GEO) (90) under accession number GSE130178. To identify differentially expressed transcripts, we utilized three different approaches, namely CLC Genomic Workbench version 12.0 (Qiagen), DeSeq2, and EdgeR. For DeSeq2 and EdgeR, the reads from the FASTQ files were mapped against the GRCm38.p6 mouse reference genome by using STAR (https://github.com/ alexdobin/STAR/releases) 7 version 2.5.2. The output files in BAM format were analyzed by featureCounts (http://bioinf. wehi.edu.au/featureCounts/) 7 to quantify the transcript abundance. Only transcripts showing differential expression using all three methods are reported here. IPA (Qiagen) was used to classify transcripts into pathways whose significance was adjusted for false discovery using the Bonferonni-Hochberg correction and a p value Ͻ0.05. The following parameters were recorded for each pathway: p value, ratios of dysregulated transcripts to all transcripts associated with that pathway (ratio), and predicted pathway activation, inhibition, and indeterminate (Z-score).

Statistical analyses
All data were analyzed using GraphPad Prism version 7 (GraphPad Software, Inc., San Diego, CA). For Kaplan-Meier survival curves, the log-rank (Mantel-Cox) test was used to calculate significance between different cohorts. One-way analysis of variance was applied for multiple comparisons using Fisher's least significant difference test. Student's two-tailed t test was used for comparing differences between two groups. All results are presented as individual values with bars indicating mean Ϯ S.D.