A Sleeping Beauty mutagenesis screen reveals a tumor suppressor role for Ncoa2/Src-2 in liver cancer

The Sleeping Beauty (SB) transposon mutagenesis system is a powerful tool that facilitates the discovery of mutations that accelerate tumorigenesis. In this study, we sought to identify mutations that cooperate with MYC, one of the most commonly dysregulated genes in human malignancy. We performed a forward genetic screen with a mouse model of MYC-induced liver cancer using SB-mediated mutagenesis. We sequenced insertions in 63 liver tumor nodules and identified at least 16 genes/loci that contribute to accelerated tumor development. RNAi-mediated knockdown in a liver progenitor cell line further validate three of these genes, Ncoa2/Src-2, Zfx, and Dtnb, as tumor suppressors in liver cancer. Moreover, deletion of Ncoa2/Src-2 in mice predisposes to diethylnitrosamine-induced liver tumorigenesis. These findings reveal genes and pathways that functionally restrain MYC-mediated liver tumorigenesis and therefore may provide targets for cancer therapy.

The Sleeping Beauty (SB) transposon mutagenesis system is a powerful tool that facilitates the discovery of mutations that accelerate tumorigenesis. In this study, we sought to identify mutations that cooperate with MYC, one of the most commonly dysregulated genes in human malignancy. We performed a forward genetic screen with a mouse model of MYC-induced liver cancer using SBmediated mutagenesis. We sequenced insertions in 63 liver tumor nodules and identified at least 16 genes/loci that contribute to accelerated tumor development. RNAi-mediated knockdown in a liver progenitor cell line further validate three of these genes, Ncoa2/ Src-2, Zfx, and Dtnb, as tumor suppressors in liver cancer. Moreover, deletion of Ncoa2/Src-2 in mice predisposes to diethylnitrosamineinduced liver tumorigenesis. These findings reveal genes and pathways that functionally restrain MYC-mediated liver tumorigenesis and therefore may provide targets for cancer therapy.
cancer gene identification | mouse models of cancer | Myc oncogene | hepatocellular carcinoma T ransposable elements (TEs) are powerful genetic tools that are widely used in insertional mutagenesis screens because of the facile identification of transposon-induced mutations. The application of transposon-based approaches to cancer gene identification provides opportunities to study the consequences of specific mutations in the context of tumor cell initiation, progression, and maintenance in well-defined, genetically engineered mouse models. Sleeping Beauty (SB), a member of the Tc1/mariner superfamily of DNA transposons, is highly active in mammalian cells (1). A growing body of evidence has demonstrated that SB is an efficient tool for cancer gene discovery when used in forward genetic screens in mice (2,3).
The SB system is based on the use of two transgenic mouse lines, one harboring a transposase and the other with a concatemerized mutagenic transposon containing sequences that can disrupt gene function through either gain-of-function or loss-offunction mechanisms. The transposase binds to the transposon ends and catalyzes its mobilization to new sites. The effect of ubiquitous transposition in mice is the development of T-cell leukemia and brain tumors (3,4). In contrast, restricted expression of SB transposase accelerates fibrosarcoma development in p19 Arf −/− mice (2).
Recently, this approach has been modified through the conditional activation of the SB transposase in specific tissues using Cre/LoxP technology (5,6). In one of these studies, an SB screen identified 19 genes that accelerate liver cancer induced by expression of a dominant-negative Trp53 transgene. However, the SB system has not been used previously to identify genes that cooperate with oncogenes that initiate liver tumorigenesis.
Hepatocellular carcinoma (HCC) is the fifth most common solid tumor worldwide and is the third leading cause of death from cancer (7,8). In the majority of cases, it occurs in a setting of chronic inflammation or cirrhosis. Although numerous genomic alterations have been documented in liver cancer, the key genetic alterations that drive hepatocellular transformation remain poorly understood. Nevertheless, a unifying feature of human HCC is amplification or overexpression of the MYC oncogene (8,9). The need to identify genetic abnormalities and to understand better how they drive cellular transformation is underscored by the fact that HCC often is diagnosed at an advanced stage when no longer amenable to curative therapies. Thus, the identification of critical genes and pathways dysregulated in HCC may accelerate the development of therapeutic strategies.
We sought to identify mutations that cooperate with MYC, one of the most commonly dysregulated oncogenes in human malignancy. We bred mice containing the SB DNA transposon vector transgene and an SB transposase transgene to a well-characterized model of MYC-induced liver cancer (10)(11)(12). We hypothesized that mice harboring the active SB transposon would accumulate mutations that accelerate the rate of tumorigenesis in this sensitized liver tumor model. A cohort of quadruple-transgenic mice was generated, each containing the SB transposon harboring a gene trap element (T2/Onc), a ubiquitously expressed transposase (Rosa26-SB11), a tetracycline (tet)-repressible MYC transgene (tet-o-MYC), and a liver-specific tet-transactivator (LAPtTA) protein. Triple-transgenic mice lacking either the transposon or transposase served as a control cohort. We observed a significant acceleration of tumorigenesis in quadruple-transgenic mice compared with control animals. Sequencing of tumors from quadruple-transgenic animals identified at least 16 genes/loci associated with accelerated liver tumor development. Validation studies in mouse liver progenitor cells (13)(14)(15) confirmed that knockdown of selected genes promotes tumorigenesis in xenograft assays. Finally, we further characterized one of the genes, Nuclear receptor coactivator2 (Ncoa2/Src-2, hereafter referred to as "Ncoa2"), which facilitates transcription of Glucose 6 phosphatase (G6pc) and thus plays a critical role in mammalian glucose availability. We find that low levels of NCOA2 and G6PC expression in HCC patients are associated with poor survival. Moreover, deletion of Ncoa2 in mice promotes diethylnitrosamine (DEN)-induced liver tumorigenesis. These findings reveal genes that function to restrain MYC-induced liver tumorigenesis in multiple model systems, thereby revealing pathways that may be attractive targets for therapy.

SB-Induced Mutations Accelerate Tumor Development in tet-o-MYC;
LAPtTA Mice. We used a previously described model of liver cancer in which mice harboring tet-o-MYC are crossed with mice expressing the tet-transactivator protein (tTA) driven by the liver-activator protein (LAP) promoter (10)(11)(12). Upon removal of doxycycline, double-transgenic animals express MYC specifically in the liver and subsequently develop liver tumors that resemble HCC. These mice served as a sensitized model of HCC.
Before proceeding with the transposon-mediated forward genetic screen, we performed a control experiment to determine tumor latency in the MYC liver model. The purpose of this experiment was twofold. First, the control experiment allowed us to distinguish early-developing tumors that form because of SB insertions from tumors that develop at baseline in tet-o-MYC; LAPtTA mice. Second, this experiment also provided a control for genetic background, because the SB mice are maintained on a mixed FVB/N and C57BL/6J background. We first crossed tet-o-MYC; LAPtTA mice, which are maintained on an FVB/N background, with C57BL/6J WT mice. Because of the positioning of the MYC transgene on the Y chromosome, only males are assessed in this model. MYC was induced when progeny reached 6 wk of age; mice then were dissected and scored for tumors each week from week 11 to week 22 (Fig. S1). Double-transgenic mice developed liver tumors no earlier than 13 wk on this mixed genetic background. Based on these results, we elected to analyze all mice before 16 wk of age, a time point at which liver tumor burden is low in tet-o-MYC; LAPtTA mice ( Fig. S1).
We bred mice containing the SB DNA transposon transgene to the MYC-induced liver tumor model. (The crossing scheme is depicted in Fig. S2A and the mutagenic transposon is depicted in Fig. S2B.)Transposase expression in liver had been demonstrated previously in Rosa26-SB11 mice (4). We hypothesized that mice harboring an active SB transposon would accumulate mutations that accelerate rates of liver tumorigenesis. An experimental cohort of quadruple-transgenic mice was generated, each containing the SB transposon harboring a mutagenic gene trap (T2/Onc), a ubiquitously expressed transposase (Rosa26-SB11), tet-o-MYC, and LAPtTA. We also generated a tripletransgenic control cohort lacking either the transposon or transposase. Doxycycline was withdrawn at 6 wk, and animals were monitored by palpation every 2 or 3 d for early-developing tumors. All animals were euthanized before 16 wk, by which time only 29% of triple-transgenic animals had developed liver tumors ( Fig. 1 A and B and Tables S1 and S2). In contrast, 65% of the quadruple-transgenic mice developed liver tumors by this time point (P < 0.003) (Fig. 1B).
To confirm T2/Onc transgene concatemer mobilization, we performed a PCR excision assay on liver tumor genomic DNA (gDNA) isolated from quadruple-transgenic and control mice. We verified excision of the transgene in all liver tumors examined (Fig. 1C). In contrast, excision of T2/Onc was not observed in a triple-transgenic control animal lacking Rosa26-SB11 transposase (Fig. 1C). We dissected 63 liver tumors from 24 quadrupletransgenic animals (representative pictures of multifocal tumors are shown in Fig. 1D). The gross morphology as well as the age of onset of these tumors is distinct from those observed in Trp53 −/− mice by Keng et al. (5). Histological analysis confirmed that these tumors resembled human HCC (Fig. 1D).
Identification of Common Insertion Sites in Liver Tumors. To identify mutations that cooperate with MYC in the development of accelerated liver tumors, common insertion sites (CISs) were identified by sequencing tumors from quadruple-transgenic mice. We first performed ligation-mediated PCR with 63 liver tumor gDNA samples (Fig. S3) followed by 454 pyrosequencing. More than 90,000 individual sequencing reads were generated, of which 98.6% matched to barcodes that were incorporated into the ligation-mediated PCR. SB has a strong tendency to favor transposition to sites adjacent to the donor concatemer, a phenomenon termed "local hopping" (16). To control stringently for this bias, we initially excluded insertions that mapped to the same chromosome as the T2/ Onc donor concatemer (chromosome 15). We further eliminated insertions that did not map to the canonical TA dinucleotide SB insertion site. We then merged all insertions that mapped to the same nucleotide from the same tumor and analyzed a total of 2,090 nonredundant insertions.
The CISs represent potential tumor-modifier genes and are defined by the presence of SB insertions at these sites significantly more frequently than would be predicted by random chance. Using   Monte Carlo simulation analysis (5, 6), we defined a statistically significant CIS as a region in the genome with three or more SB insertions within a 30-kb window or four or more insertions within 155 kb. Initially, we identified 23 CISs that fit these criteria. Nine of these CISs occurred at the same nucleotide in different tumors from the same mouse and seven were eliminated from consideration as early-occurring passenger insertions less likely to be related to tumor development. Two of these CISs, Zfp189 and Bcl9, remain on the list, because they exhibited concordance with human tumor data. The remaining 14 genes/loci defined by CISs in MYC-induced liver tumors are shown in Table 1, along with two others discovered as outlined below. Consistent with previous studies, we find that genome-level integration of SB with respect to genes and intergenic regions is fairly random overall (Fig. S4A) (17)(18)(19). As expected, we observe the greatest number of insertions clustering around the T2/ Onc transgene concatemer on chromosome 15 (Fig. S4B). Upon closer inspection of transposon integrations on this chromosome, we identified two additional genes as potential CISs that initially were excluded. These CISs correspond to the Growth Hormone Receptor (Ghr) and Nuclear Factor (erythroid-derived 2) genes ( Fig. S4B and Table 1). Because these genes are located many megabases from the donor transgene, these transposon insertions cannot be explained solely by local hopping. Moreover, GHR was implicated previously in human liver cancer (20).
Insertions in candidate liver cancer genes identified in our screen are shown in Fig. 2A and Fig. S5A. Two examples are the Nuclear receptor coactivator2 (Ncoa2) and the Zinc finger transcription factor (Zfx). Insertions in these candidates are distributed randomly throughout the gene and occur in either a sense or antisense orientation. Thus, these insertions are likely to disrupt gene function. To confirm transposon integrations in Ncoa2, we performed PCR assays with a gene-specific primer and a T2/Onc primer. In each case, we observed the SB insertion in liver tumor gDNA but not in tail gDNA from a C57BL/6J wild-type animal (Fig. 2B).
It has been demonstrated previously that Ncoa2, a member of the p160 family of transcriptional coactivators, modulates expression of Glucose 6 phosphatase (encoding G6Pase) and other important genes by acting in concert with nuclear receptors such as RORα (21,22). We measured mRNA expression levels of Ncoa2 and G6pc in four liver tumors with Ncoa2 insertions. Tumors with SB insertions exhibited a reduction of both Ncoa2 and its target G6pc mRNA levels (Fig. 2C). In contrast, tumors lacking SB insertions in Ncoa2 had normal or elevated expression of this gene (Fig. S6), demonstrating that Ncoa2 levels are not naturally lower in the cell of origin of HCC in this model.

SB Insertions in Lymphomas Are Distinct from Liver Tumor Insertions.
Twenty-nine percent of the quadruple-transgenic animals (11 of 37) developed lymphomas before 16 wk of age, whereas none of the triple-transgenic animals exhibited signs of hematopoietic disease. This result was not unexpected, because T2/Onc; Rosa26-SB11 mice are known to develop lymphoma and high-grade gliomas on a wild-type genetic background (4). Seven of these 11 quadrupletransgenic animals developed both liver tumors and lymphomas. To distinguish liver tumor-specific insertions from lymphoma-specific insertions in these seven mice, we performed ligation-mediated PCR with spleen gDNA and then sequenced SB insertions in these samples. Similar to our analysis with liver tumor SB insertions, we searched for CISs that occur at a greater frequency than would be predicted by random chance. We identified 11 genes defined as CISs in these lymphomas (Table S3). Despite the small number of lymphomas analyzed, two of the genes with SB insertions, Myb and Erg, had been found previously as CISs that contribute to leukemias/lymphomas in a larger cohort of 58 T2/Onc; Rosa26-SB11 animals (4). In contrast, none of the lymphoma CIS genes from the Collier et al. study (4) overlapped with liver CIS genes identified in the present study. Interestingly, the only CIS that overlapped between our lymphoma-and liver tumor-sequencing datasets was Ncoa2 (Fig. S5B). This finding underscores the potential significance of this gene in MYC-induced tumorigenesis in multiple settings. Consistent with this, we observed a robust increase in MYC protein levels in SB-induced lymphomas compared with normal spleen samples from control animals (Fig. S7). Although the source of elevated MYC in these tumors remains to be determined, it may derive from either the endogenous or transgenic loci.
Functional Validation of CIS Genes Implicates Ncoa2, Zfx, and Dtnb as Tumor Suppressors. To validate functionally the tumor-suppressor activity of several of the CISs that we identified in our SB screen Nfe2 Nuclear factor erythroid-derived 2 4708 7 7 15 Ghr Growth hormone receptor 259059 7 7 Common insertion sites in liver tumors from quadruple-transgenic mice defined using Monte Carlo simulation analysis. Independent insertions refer to nonredundant insertions. The CISs below the black line represent genes with insertions at the same nucleotide in tumors from the same mouse or were located on chromosome 15. These CISs remain on the list because their expression showed concordance with HCC data. in an independent system, we used hepatoblasts isolated from Trp53 −/− mice and immortalized by MYC. These cells do not form tumors in immunocompromised mice but provide a sensitized background in which one additional genetic alteration suffices to drive tumorigenesis (13)(14)(15). Furthermore, the Myc retrovirus used to immortalize these liver progenitor cells coexpresses GFP, allowing fluorescence imaging of subcutaneous tumors. We modified these cells through the introduction of shRNAs corresponding to candidate genes, followed by transplantation into recipient mice (Fig. 3A). Using this approach, we assess directly whether loss of function of candidate genes influences tumorigenesis and, by extension, whether these genes suppress tumor initiation. As a positive control, we knocked down the Wnt signaling pathway component Apc, a known tumor suppressor in this system (Fig. S8A). We confirmed that shRNAs targeting Apc in Trp53 −/− ; Myc hepatocytes gave rise to tumors within 1 mo after injection into nude mice (Fig. S8 B and C).
We selected five genes from our SB screen to evaluate in this assay. We chose to focus on genes that had not been linked previously to liver cancer and in which the mouse shRNAs were readily available from the RNAi Consortium. We tested four or five different shRNAs per gene and chose the most potent shRNAs for our in vivo validation studies ( Fig. 3B and Fig. S9 A and B). We then infected liver progenitor cells with lentiviruses expressing these shRNAs and assayed for tumor formation. Knockdown of Zfx, one of the top hits from our screen, promotes tumor formation of liver progenitor cells in nude mice (Fig. 3 C and D). We also validated Ncoa2 and β-Dystrobrevin (Dtnb) as putative tumor-suppressor genes using this approach (Fig. 3D). In contrast, the knockdown of two other genes represented by CISs, Zfp189, and Rc3h1, did not result in tumor formation in this assay.    tensive mRNA expression profiling dataset from a large cohort of human liver tumor samples (23). We analyzed expression of 11 orthologs of the CIS genes identified in our screen in a cohort of 139 human HCC samples. Of these, five genes were expressed at statistically significantly lower levels in human tumors than in normal human liver tissue: ZNF608, NCOA2, CTNNA3, SFI1, and DTNB (Fig. 4 A and B). Three genes exhibited a statistically significant increase in expression in HCC samples: SND1, ZNF189, and Wnt pathway component BCL9 (Fig. 4 A and B). Furthermore, we observed an association between the expression of NCOA2 and its target G6PC and overall survival of HCC patients. Kaplan-Meier plot analysis and a log-rank test for overall survival in HCC based on the expression levels of these genes either independently or together revealed a statistically significant reduction in survival of HCC patients with tumors that exhibit low NCOA2 and/or G6PC expression (Fig. 4C and Fig.  S10). These findings demonstrate that our transposon-mediated forward genetic screening approach identified clinically relevant genes that participate in the pathogenesis of human liver cancer.
We next assessed whether the NCOA2 gene accumulated sequence variants in a smaller panel of 14 human liver tumor samples. Paired normal liver biopsy material was available for a subset of these HCC samples. We sequenced all coding exons of NCOA2 and found one potential case of loss of heterozygosity specific to the tumor (Fig. S11A). Consistent with deletion of one allele of NCOA2 in this tumor, expression of NCOA2 was reduced by 50% in this sample relative to its paired normal liver (Fig. S11B). In contrast, NCOA2 expression was not decreased to the same extent in the other 13 liver tumor samples examined relative to matched normal controls (Fig. S11C). These results suggest that hemizygosity or other alterations affecting the NCOA2 gene may contribute to human liver tumorigenesis.
Loss of NCOA2 Enhances DEN-Induced Liver Tumorigenesis. To test directly whether NCOA2 functions as a putative tumor suppressor in vivo, we treated previously described Ncoa2 +/+ and Ncoa2 −/− mice (24) with the chemical carcinogen DEN. Exposure to DEN causes liver damage and hepatocellular carcinogenesis in mice (25,26). We chose to analyze males because DEN-induced tumorigenesis is known to be more effective in males than in females and because males were analyzed in the original transposon screen. pups were injected with 25 mg/kg DEN on postnatal day 15 and then were euthanized at 3 or 6 mo postinjection. For each mouse, the liver was isolated, weighed, and photographed macroscopically ( Fig. S12A and Table S4). DEN-treated livers were evaluated further to determine the total tumor number as well as the average and maximal tumor size. We observed that Ncoa2 −/− male mice exhibited an approximately fourfold increase in tumor number and a statistically significant increase in maximal and average tumor size by 6 mo after DEN treatment (P < 0.05) ( Fig.  5A and Fig. S12B). These findings demonstrate that genetic loss of NCOA2 promotes liver tumorigenesis in a mouse model that closely mimics human HCC pathogenesis.

Discussion
Understanding molecular mechanisms that contribute to tumor development may provide therapeutic approaches for cancer therapy. Screens using transposons facilitate the discovery of genes and pathways contributing to tumor development in mice (5,6,27,28). Thus, forward genetics provides a complementary approach to cancer genome-sequencing studies by illuminating the functional relevance of genes mutated in human tumors. We used transgenic mouse lines that express a ubiquitous transposase and a mutagenic transposon present at low copy number to screen for insertions that accelerate MYC-induced HCC. We found a dramatic increase in the number of quadruple-transgenic animals that developed liver tumors compared with triple-transgenic controls that also expressed the tet-o-MYC and LAPtTA transgenes but lacked an active SB transposon. Sequencing of insertions from liver tumors identified at least 16 CISs/loci that contribute to accelerated tumor development.
Several of the CISs identified in our genetic screen have human orthologs dysregulated in human liver cancer. The genes ZNF608, NCOA2, CTNNA3, SFI1, and DTNB were all expressed at statistically significantly lower levels in human HCC samples. Our validation studies in p53 −/− ; Myc liver progenitor cells further demonstrated that inhibition of Ncoa2 and Dtnb can promote tumor formation in immunocompromised mice. DTNB, a member of the dystrobrevin protein family originally identified as components of the dystrophin-glycoprotein complex, is expressed in nonmuscle tissues including the liver, kidney, and brain (29). Dtnb may have previously unrecognized antitumorigenic activities. Interestingly, DTNB interacts with cAMP-dependent PKA and is proposed to play a scaffolding role in intracellular signaltransduction pathways (30). Atbf1, which encodes AT-motif binding factor-1 protein, also was identified in our screen. Also known as "zinc-finger homeobox 3" (ZFHX3), this protein has been shown previously to negatively regulate the expression of the liver tumor marker α-fetoprotein (AFP) in hepatocytes (31). Moreover, it was demonstrated previously that ATBF1 expression was significantly reduced in 55 of 76 HCC samples (32).
We also validated the tumor-suppressor activity of Zfx, encoding an X-linked pluripotency-associated zinc finger transcription factor that controls the self-renewal of embryonic and hematopoietic stem cells (33), as well as B-cell proliferation (34), in our functional studies. Although ZFX has not been linked previously to human liver tumors, and we did not observe altered expression of this gene in HCC samples, a potential link between the transcriptional networks coregulated by MYC and ZFX was discovered recently. A strong correlation in target gene occupancy was identified between these transcription factors (35). Transcriptional targets of these proteins are enriched for genes involved in the cell cycle, cell death, and cancer. Although future studies are necessary to elucidate the functional relationship between ZFX and MYC, previous findings and those reported here suggest that ZFX may directly oppose the activity of MYC at key transcriptional targets. Loss of ZFX may therefore remove an important negative regulator of MYC-mediated transactivation of oncogenic targets.
Genes with human orthologs that show increases in the level of mRNA expression in human liver tumor samples include SND1, ZNF189, and the Wnt pathway component BCL9. Dysregulation of the Wnt signaling pathway is a key molecular lesion in liver cancer (36)(37)(38). More than 60% of liver tumors exhibit accumulation of β-catenin, a hallmark of activated Wnt signaling. Interestingly, BCL9 was shown previously to be required for β-catenin-mediated transcriptional activity in human cancer cell lines (39). Recent findings suggest that BCL9 is an oncogene, because it has been shown to increase cell proliferation, migration, invasion, and metastatic potential of tumor cells (40). One additional gene identified in our screen, Ctnna3, also interacts with the Wnt pathway (41). Further analysis of all CIS genes using the COSMIC database analysis revealed the presence of additional mutations and/or genomic alterations that may be relevant to the pathogenesis of diverse human tumor types (Table S5). In light of the data presented here, these variants are worthy of further investigation.
An independent SB mutagenesis screen using a dominant-negative Trp53 transgene recently identified 19 genes associated with liver cancer (5). Interestingly, the only CIS gene in common between this study and the Keng study was Sfi1 (5). A meta-analysis of several SB screens indicates that this gene scores very highly in many such screens. Recent work suggests that this finding reflects an artifact caused by previously unrecognized amplification of this gene through segmental duplications and/or other mechanisms in the C57BL/6 mouse genome (42). Nevertheless, we found that SFI1 mRNA expression was down-regulated in a panel of 139 human HCC samples. Consistent with this result, Keng et al. (5) found distinct copy-number losses of SFI1 in 20 of 28 human HCC samples examined. SFI1 is a centrosomal protein with well-studied homologs present in yeast. Although SFI1 has been proposed to Gene-expression studies have demonstrated that NCOA2 has many additional targets, including genes involved in cell-cycle regulation, signal transduction, apoptosis, immunity, and transport, which also may contribute to tumorigenesis. play a role in the formation of centrosome-associated contractile fibers (43), and the data presented here and in Keng et al. (ref. 5) are most consistent with a potential tumor-suppressive role for SFI1, additional studies are needed to determine the role of this gene in chromosome stability and to elucidate how its reduced expression may contribute to tumorigenesis. Transposon insertions in Egfr were very common in liver tumors isolated from mice harboring a dominant-negative Trp53 transgene and a conditional SB transposase (5). In fact, SB insertions were detected most frequently in intron 24 of Egfr, giving rise to a truncated Egfr protein that has a strongly oncogenic effect. Using a PCR genotyping assay, we screened all liver tumors isolated in our genetic screen and did not find a single occurrence of SB insertions in this intron. Presumably, MYC overexpression performs a function that obviates selection for Egfr activation in this model, perhaps because EGFR signaling results in MYC upregulation through the RAS-RAF-MEK-MAPK pathway (44,45). Moreover, this finding underscores the importance of performing screens of this type using multiple sensitized models, because the most relevant pathways in a given tumor may depend on the initiating event that drives tumorigenesis.
Perhaps most intriguing, our findings provide evidence that the transcriptional coactivator NCOA2 acts as a tumor suppressor in liver cancer. Four lines of data support this concept. (i) Recurrent transposon insertions in SB-induced liver tumors result in decreased mRNA expression of both Ncoa2 and its previously characterized target G6pc. (ii) Inhibition of Ncoa2 using multiple independent shRNAs promotes tumor formation of liver progenitor cells in immunocompromised mice. (iii) DEN treatment of Ncoa2 −/− mice enhances liver tumorigenesis compared with Ncoa2 +/+ littermates. (iv) Decreased expression of NCOA2 is observed in human liver tumor samples.
Previous studies of Ncoa2 illuminate our mechanistic understanding of its tumor-suppressor activity. Mice lacking Ncoa2 develop glycogen storage disease type 1 (Von Gierke's disease) and exhibit decreased expression and enzymatic activity of G6pase (21). Similarly, a genetic deficiency of G6PC leads to Von Gierke's disease both in humans and in mice (46,47). Several observations demonstrate that loss of function of G6Pase also promotes liver tumorigenesis. For example, enzymatic activity of G6Pase is reduced in human and rodent liver tumors (48,49), and G6pc −/− mice develop hepatic tumors (47). Furthermore, a significant fraction of patients with glycogen storage disease type 1 develop hepatic adenomas and have an additional risk of undergoing malignant transformation (50). Finally, inhibition of hexokinase, which has the equivalent functional effect as increasing G6Pase activity, is antitumorigenic (51,52).
Based on these findings, loss of Ncoa2 may promote liver tumorigenesis at least in part through a subsequent reduction in G6Pase activity (Fig. 5B). This reduction, in turn, may increase the intracellular pool of glucose-6-phosphate that can enter the glycolytic pathway. Our data indicate that loss of Ncoa2 is particularly protumorigenic in settings of MYC hyperactivity, which is consistent with the known ability of MYC to stimulate glycolytic metabolism potently through the direct induction of expression of numerous genes in the glycolytic pathway (53)(54)(55). Despite the potential relevance of this mechanism to tumorigenesis, we were unable to detect expression of G6Pase in the hepatoblast cells, suggesting that NCOA2 tumor-suppressor activity is independent of G6Pase in the hepatoblast system. Although this finding does not rule out the potential significance of G6Pase as a relevant target of NCOA2 in vivo, we acknowledge that mechanisms of NCOA2-mediated tumor suppression are complex and likely involve multiple targets and pathways. Indeed, a previously reported microarray analysis of gene expression in livers of Ncoa2 −/− mice revealed a signature of increased energy expenditure including the stimulation of genes involved in fatty acid degradation and suppression of genes involved in cholesterol and fatty acid biosynthesis (22). Thus, Ncoa2 deficiency likely contributes to a broader protumorigenic reprogramming of metabolism beyond the glycolytic pathway. Moreover, loss of function of Ncoa2 influences gene-expression programs relevant to other important aspects of tumor biology including cell-cycle regulation, signal transduction, and apoptosis. Therefore future studies are warranted to further dissect the mechanism(s) by which dysregulated expression of NCOA2 contributes to tumorigenesis. PCR Genotyping. Genomic DNA was isolated from tail clippings using the DNeasy Blood and Tissue Kit (Qiagen) according to the manufacturer's instructions. PCR genotyping was performed using 10-100 ng of genomic DNA as template. Primers used were tet-o-MYC, forward 5′-TAGTGAAC-CGTCAGATCGCCTG-3′ and reverse 5′-TTTGATGAAGGTCTCGTCGTCC-3′ (amplicon 500 bp); LAPtTA, forward 5′-GCTGCTTAATGAGGTCGG-3′ and reverse 5′-CTCTGCACCTTGGTGATC-3′ (amplicon 500 bp); T2/Onc, forward 5′-CCTTGCAAAATGGCGTTACT-3′ and reverse 5′-GCCGCACTGGTTGTTAGCAA-3′ (amplicon 200 bp); Rosa26-SB11 WT, forward 5′-CTGTTTTGGAGGCAGGAA-3′, reverse 5′-CCCCAGATGACTACCTATCCTCCC-3′ (amplicon 420 bp); Rosa26-SB11 SB reverse 5′-CTAAAAGGCCTATCACAAAC-3′ (SB knock in amplicon 317 bp ); Hprt, forward 5′-CTTTCTTGTCACTCCCACTTTTCC-3′ and reverse 5′-CACATCCTTCATTCAGGTGTCACT-3′ (amplicon 340 bp). PCR conditions for Ex Taq (Takara) were used according to the manufacturer's instructions. PCR products were separated on a 1.2% agarose gel.

Methods
Liver Tumor Analysis. The whole liver was dissected from the euthanized animal, washed, and placed in ice-cold PBS. The number of liver tumor nodules was counted for each animal. Tumor nodules then were dissected away from surrounding normal liver tissue. When large enough, tumor nodules were split into three pieces for gDNA isolation, RNA isolation, and histological analysis. DNA isolation was performed using the All-in-One Purification Kit (Norgen) or by phenol-chloroform extraction. Total RNA was isolated from tumor and normal tissue using TRIzol (Invitrogen) followed by additional cleanup and DNase digestion using the RNeasy kit (Qiagen) according to the manufacturers' protocols. For histological analysis, tissues were fixed in 10% formalin, embedded in paraffin, and sectioned. Tissue section slides then were stained with H&E using standard protocols.
Ncoa2 Insertion Validation Assay. For each insertion that was validated, 1 kb of sequence surrounding each insertion position was downloaded (500 bp on each side of the insertion), and an ∼800-bp amplicon was designed. To confirm the orientation of the SB insertion, different combinations of Ncoa2 and T2/Onc primers were used in the PCR. Primers used to confirm SB insertion at position 13272606 were forward 5′-CACACACACACACAGCA-GTAGTGAAGCA-3′ and reverse 5′-CGCCAGAAAACTGATTTTGAAGTTGCAG-3′; to confirm the SB insertion at position 13300371, the primers used were forward 5′-GAGCGATGGTTTCTACTCCCACAGACCT-3′ and reverse 5′-CTGTA-CCGTGGTGAGCATTGGGAGTAAT-3′.
Sequence Processing and Annotation. The 47,343 initial 454 sequencing reads were processed and analyzed using a custom semiautomated processing pipeline first described in ref. 6. Briefly, Fasta-formatted sequences from 454 runs were assigned to libraries allowing one mismatch in the 10-bp barcode. All raw sequences were scanned for inverted repeat/direct repeat (IRDR) and linker recognition sequences using EMBOSS Vectorstrip (56) with successively less stringent parameters (10%, 15%, and finally 20% mismatches allowed) until the maximum number of constructs (i.e., IRDR: GTGTATGTAAACTTCC-GACTTCAACTG and linker sequence: GTCCCTTAAGCGGAGCCCT) were recognized and trimmed off. Only 44,294 sequences that had a matched IRDR with an insert sequence of at least 16 bp were carried further for mapping to the mouse genome. These trimmed insert sequences were mapped to the mouse genome (National Center for Biotechnology Information build 37) using BLASTN (DeCypher's TeraBLASTN; Active Motif), requiring query sequences to align within 1 bp of the end of the LTR sequence that was trimmed. Additionally, the query sequence was required to match with at least 95% identity. Ambiguous sequences that mapped to multiple genomic loci were removed, leaving 2,090 uniquely mappable nonredundant insertions. Unambiguously mapped nonredundant insertions were assigned to clusters of CISs if the local density of insertions in a given window size exceeded the density that would be expected by chance, as determined by Monte Carlo simulation. A CIS is defined as statistically significant if a region in the genome has three or more SB insertions with a 30-kb window or four or more insertions within 155 kb. CISs that occurred at the same nucleotide in different tumors from the same mouse were removed with the exceptions of Zfp608 and Bcl9. Nonredundant insertion positions were annotated using the EnsEMBL API (57) with the name of the gene whose start site was closest to the position. CIS positions were annotated similarly using the median insertion position within the CIS as a reference point.
RT-PCR Analysis of CIS Genes. All quantitative real-time PCR analysis was performed using the Step One Plus Real Time PCR system (Applied Biosystems). Analysis of Ncoa2 in SB liver tumors was performed using a predeveloped TaqMan assay with One-Step RT-PCR Master Mix Reagents (Applied Biosystems) or by SYBR green analysis of cDNA samples using primers that cross an exon-exon junction. The Ncoa2 TaqMan assay was performed with Mm00500749_m1 (Applied Biosystems). Primers used for SYBR green were forward 5′-TCACTGCATTGGCTCTTCTG-3′ and reverse 5′-TGT-TTTCTCCCATCCCACTC-3′. Analysis of G6Pase mRNA in SB liver tumors was performed with a TaqMan assay using the Universal Probe #26 (Roche) and forward 5′-GAGGAAAGAAAAAGCCAACGTA-3′ and reverse 5′-GGGACAGA-CAGACGTTCAGC-3′ primers. Expression was normalized to 18SrRNA with the Ribosomal RNA control assay (Applied Biosystems) or Actin. All PCR assays were performed in triplicate.
Xenograft Assay. Trp53 −/− ; Myc cells (2 × 10 6 ) expressing shRNA lentiviruses were injected subcutaneously into the right flank of 6-to 8-wk-old female nude mice (National Cancer Institute-Frederick). Tumor volume was measured with calipers every 3-5 d until the average tumor mass reached 2 cm 3 . Tumor volume was calculated with the formula (length × width 2 )/2. A total of four mice were injected per shRNA tested. At least two independent tumorigenesis experiments were performed for each gene.
GFP Imaging. Whole-body GFP imaging was performed using the IVIS Spectrum Optical Imaging system (Xenogen) according to the manufacturer's instructions.
Microarray Analysis of Human HCC Samples. The mRNA expression data of 12 genes in human HCC were obtained from Gene Expression Omnibus with accession nos. GSE4024 and GSE1898 (23). Hierarchical cluster analysis was conducted using Cluster 3.0, and the image analysis was performed with TreeView 1.60 software (Michael Eisen Laboratory, Lawrence Berkeley National Laboratory and University of California, Berkeley, CA; http://rana.lbl. gov/eisen/). We applied Kaplan-Meier plot analysis and a log rank test for overall survival in HCC by using Prism v. 5.03 (GraphPad Software, Inc.).

DEN Tumorigenesis Experiment.
Fifteen-day-old male Ncoa2 +/+ and Ncoa2 −/− littermate pups were injected intraperitoneally with a 25 mg/kg dose of DEN (Sigma) in 0.9% sterile saline as previously described (58). DEN-treated mice were fed normal chow (Teklad 2920×) and euthanized at 3 or 6 mo postinjection. Mice were weighed, and liver was isolated, weighed, and photographed macroscopically. Livers were evaluated further to determine the total number of tumors as well as average and maximal tumor size using a handheld micrometer. Samples of liver tissue were prepared for histological analysis as previously described (59).