MLL3 regulates the CDKN2A tumor suppressor locus in liver cancer
Abstract
Mutations in genes encoding components of chromatin modifying and remodeling complexes are among the most frequently observed somatic events in human cancers. For example, missense and nonsense mutations targeting the mixed lineage leukemia family member 3 (MLL3, encoded by KMT2C) histone methyltransferase occur in a range of solid tumors, and heterozygous deletions encompassing KMT2C occur in a subset of aggressive leukemias. Although MLL3 loss can promote tumorigenesis in mice, the molecular targets and biological processes by which MLL3 suppresses tumorigenesis remain poorly characterized. Here, we combined genetic, epigenomic, and animal modeling approaches to demonstrate that one of the mechanisms by which MLL3 links chromatin remodeling to tumor suppression is by co-activating the Cdkn2a tumor suppressor locus. Disruption of Kmt2c cooperates with Myc overexpression in the development of murine hepatocellular carcinoma (HCC), in which MLL3 binding to the Cdkn2a locus is blunted, resulting in reduced H3K4 methylation and low expression levels of the locus-encoded tumor suppressors p16/Ink4a and p19/Arf. Conversely, elevated KMT2C expression increases its binding to the CDKN2A locus and co-activates gene transcription. Endogenous Kmt2c restoration reverses these chromatin and transcriptional effects and triggers Ink4a/Arf-dependent apoptosis. Underscoring the human relevance of this epistasis, we found that genomic alterations in KMT2C and CDKN2A were associated with similar transcriptional profiles in human HCC samples. These results collectively point to a new mechanism for disrupting CDKN2A activity during cancer development and, in doing so, link MLL3 to an established tumor suppressor network.
Editor's evaluation
This paper convincingly shows that MLL3 regulates the CDKN2A tumor suppressor in MYC-driven liver cancers. The use of in vivo models and epigenomic analysis made the findings particularly robust. This work significantly advances our understanding of the function of MLL3 in cancer.
https://doi.org/10.7554/eLife.80854.sa0Introduction
Hepatocellular carcinoma (HCC) is a deadly primary liver cancer with a 5 year survival rate of only 18% (Jemal et al., 2017). HCC is currently the fourth most frequent cause of cancer-related mortality worldwide, and its incidence continues to grow (Llovet et al., 2021). Genomic alterations found in HCC are highly diverse and are characterized by promoter mutations in TERT (telomerase reverse transcriptase), amplifications, or chromosomal gains encompassing the MYC oncogene, activating hotspot mutations in CTNNB1 (β-catenin), and inactivating mutations and deletions in the TP53 and CDKN2A tumor suppressor genes (2017; Schulze et al., 2015).
Among these alterations, genetic gain of MYC and inactivation of tumor suppressor p53 are known to cooperate to drive tumorigenesis in HCC (Molina-Sánchez et al., 2020). Mechanistically, oncogenic MYC activation triggers increased expression of the tumor suppressor ARF, one of two proteins encoded in CDKN2A in alternative reading frames. ARF binds to the E3 ubiquitin ligase MDM2 to prevent p53 degradation, leading to apoptosis to restrain MYC-driven tumorigenesis (Lowe and Sherr, 2003). However, it is unclear how the CDKN2A locus is regulated in response to MYC overexpression.
Beyond these well-studied drivers, HCC frequently harbors mutations in one or more chromatin modifying enzymes, including MLL3 (encoded by KMT2C; Fujimoto et al., 2012; Kan et al., 2013). MLL3 is a component of the COMPASS-like complex that has structural and functional similarities to the developmentally essential Drosophila Trithorax-related complex (Schuettengruber et al., 2017). This multiprotein complex controls gene expression through its histone H3 lysine 4 (H3K4) methyltransferase activity, which establishes chromatin modifications most often associated with transcriptional activation (Shilatifard, 2012). Most studies have shown that MLL3 and its paralog MLL4 (encoded by KMT2D) typically catalyze H3K4 monomethylation (H3K4me1) at enhancers (Herz et al., 2012; Hu et al., 2013), while the MLL1/2 complex is responsible for H3K4 trimethylation (H3K4me3) at promoters and enhancers in a locus-specific manner (Denissov et al., 2014; Rickels et al., 2016; Wang et al., 2009).
While less characterized, MLL3/4 regulation of promoter activity is emerging as an additional mechanism connecting the COMPASS-like complex to gene expression. Some publications report that H3K4me1 enrichment at promoters has been associated with gene repression (Cheng et al., 2014), and MLL3 inactivation decreases H3K4me3 levels at the promoters of metabolism-related genes in normal murine livers (Valekunja et al., 2013) and human liver cancer cells (Ananthanarayanan et al., 2011). Furthermore, a recent study in leukemia cells demonstrated that MLL3 and MLL4, in the absence of MLL 1/2 complex, are capable of binding to promoters to activate tumor suppressor genes (Soto-Feliciano et al., 2023). These divergent results suggest that the genomic binding pattern and functions of MLL3 are highly context dependent.
Notably, HCC also harbors mutations in KMT2D (Cleary et al., 2013), while KDM6A/UTX, an H3K27 demethylase within the COMPASS-like complex, has been functionally established as a potent tumor suppressor in pancreatic and liver cancers (Revia et al., 2022). These observations suggest that epigenetic-based mechanisms of gene regulation controlled by the MLL3 complex may constrain HCC development. However, because chromatin regulators such as ARID1A often exhibit context-specific tumor suppressive and oncogenic roles in liver cancer development (Sun et al., 2017), it is unclear whether MLL3 is a bona fide tumor suppressor in HCC. We therefore employed mouse models of HCC to investigate the molecular targets of MLL3 and the biological processes it affects.
Results
MLL3 is a tumor suppressor in Myc-driven liver cancer
To better understand the functional significance of genes commonly inactivated in HCC, including a number of chromatin regulators, we selected 12 genes with recurrent inactivating mutations in human HCC (Cancer Genome Atlas Research Network, 2017; Ahn et al., 2014; Fujimoto et al., 2012; Figure 1—figure supplement 1A) and performed a CRISPR-based in vivo screen to determine whether they behave as tumor suppressors in HCC. Specifically, the screen tested whether loss of each of these 12 genes would drive hepatic tumorigenesis in cooperation with Myc—one of the most frequently gained and/or amplified oncogenes in HCC (Huang et al., 2014). We applied hydrodynamic tail vein injection (HTVI) in wild-type mice to directly introduce genetic manipulations into adult hepatocytes in vivo (Bell et al., 2007). We introduced both a transposon vector for stable genomic integration of oncogenic Myc cDNA and plasmids designed for transient expression of Cas9 and single guide RNAs (sgRNAs; a mix of two for each gene; Figure 1B; Largaespada, 2009; Moon et al., 2019; Tschaharganeh et al., 2014; Xue et al., 2014). 3 months after HTVI, only sgKmt2c resulted in liver tumor formation with high penetrance (Figure 1—figure supplement 1B), suggesting that MLL3 likely acts as a tumor suppressor to constrain Myc-driven liver cancer. Supporting this idea, KMT2C mutations co-occur with MYC genomic gains and amplifications in human HCC tumors (Figure 1A).
To validate and extend the results from the screen, we applied the same approach to test whether the screen phenotype could be recapitulated with oncogenic Myc and single Kmt2c-targeted sgRNAs (Figure 1B). Mice injected with an Myc cDNA transposon combined with either of two independent Cas9/Kmt2c sgRNAs (Myc; sgKmt2c.1 or Myc; sgKmt2c.2) developed liver tumors, with a slightly later onset and slightly longer survival than mice receiving the Myc transposon combined with an sgRNA targeting Trp53 (Myc; sgTrp53; Figure 1C and D). In contrast, mice injected with Myc and a control sgRNA (sgChrom8) did not succumb to disease over the observation period (Figure 1C). These findings were confirmed in an independent cohort of mice (Figure 1—figure supplement 2b). Analyses of tumor-derived genomic DNA revealed insertions and deletions (indels) in either Kmt2c or Trp53 depending on the genotype of tumor-derived cells (Figure 1—figure supplement 2B). DNA sequencing of the CRISPR-targeted region from two independent Myc; sgKmt2c tumors revealed indels predicted to generate premature stop codons (Figure 1—figure supplement 2C). In one case, the indel was heterozygous, implying that even partial suppression of Kmt2c can promote tumorigenesis. In support of this, GFP-linked Kmt2c shRNAs efficiently cooperated with Myc overexpression to drive liver cancer, producing tumors with 50–80% reduction in Kmt2c mRNA expression (Figure 1—figure supplement 2D–G). shKmt2c.2 resulted in less potent knockdown than shKmt2c.1 yet produced faster tumor formation, suggesting that, as in acute myeloid leukemia (Chen et al., 2014), MLL3 can likely act as a haploinsufficient tumor suppressor in liver cancer (Figure 1—figure supplement 2E, G).
Apart from MYC, CTNNB1 (β-catenin) is a frequently mutated oncogene in human HCC (Rebouissou et al., 2016), although the co-occurrence between CTNNB1 and KMT2C mutations was not statistically significant (Figure 1A). To test whether MLL3 loss can cooperate with oncogenic CTNNB1 to promote liver tumorigenesis, we performed analogous HTVI of a transposon vector expressing constitutively active β-catenin (Ctnnb1-N90; Tward et al., 2007) in combination with Kmt2c- or Trp53-targeted sgRNAs (Figure 1—figure supplement 3A). However, no tumor formation was observed in mice that received the Kmt2c-targeted sgRNAs (Figure 1—figure supplement 3B), indicating that unlike p53, the tumor-suppressive role of MLL3 is specific to the oncogene and contexts.
MLL3 loss alters the chromatin landscape of liver cancer cells
MLL3 and MLL4 are histone methyltransferases that can deposit the H3K4 monomethylation mark at genomic enhancers and intergenic regions during organ development (Hu et al., 2013). However, more studies indicate that MLL3 and MLL4 are also capable of binding to promoter regions (Chen et al., 2014; Dhar et al., 2016; Wang et al., 2010), especially in the context of cancer (Soto-Feliciano et al., 2023). To determine the genomic binding patterns of MLL3 in HCC, we performed MLL3 chromatin immunoprecipitation (ChIP)-sequencing (ChIP-Seq) analysis in Myc; sgKmt2c (sgKmt2c.1 which generates heterozygous or homozygous indels) and Myc; sgTrp53 liver cancer cell lines. Compared to the sgTrp53 cells, sgKmt2c cells had a marked reduction in MLL3 chromatin binding at a subset of genomic loci (Figure 2A). Approximately 40% of the peaks that were selectively lost in Kmt2c-deficient cells occurred at promoter regions, whereas unchanged MLL3 peaks between the two genotypes were more likely to be within intergenic regions (Figure 2B, Figure 2—figure supplement 1). Therefore, our data suggest that, beyond the canonical action of MLL3 at gene enhancers, MLL3 can also occupy promoter regions in Myc-induced liver cancer. Of note, the residual ChIP-seq signal observed in the sgKmt2c cells most likely reflects the binding of MLL4 and/or remnant MLL3 since the antibody used in these experiments can recognize both MLL3 and MLL4 proteins (Dorighi et al., 2017). Nonetheless, the downregulated peak signals in Myc; sgKmt2c cells were specifically due to MLL3 disruption.
Similar to the Drosophila Trithorax-related complex (Schuettengruber et al., 2017), the mammalian MLL3 and MLL4 complexes facilitate gene transcription by establishing permissive modifications on histone H3K4 via the MLL3 and MLL4 methyltransferase (Shilatifard, 2012). To determine whether MLL3 disruption impacts the local or global chromatin landscape of HCC cells, we performed ChIP-Seq analyses for H3K4 methylation and H3K27 acetylation in six independently derived tumor cell lines: three each for Myc; sgKmt2c and Myc; sgTrp53 (Figure 2C). Cluster analysis on genomic areas revealed three clusters of genomic loci that showed enrichment or depletion between Myc; sgKmt2c and Myc; sgTrp53 tumor cells for each tested histone modification (Figure 2C, Figure 2—figure supplement 2A–C). Loci in cluster 1 (reduced H3K4me3, H3K4me1, and H3K27ac in sgKmt2c cells) showed the most pronounced differences in chromatin modifications between the two liver tumor genotypes. In contrast, the loci in cluster 2 showed increased H3K4me1 and H3K27ac marks, most of which mapped to intergenic regions. The loci in cluster 3 showed increased H3K4me3 and included some p53 target genes such as Cdkn1a and Eda2r.
To determine whether these drastic changes in the chromatin landscape were associated with changes in MLL3 binding, we integrated the chromatin modifications results with our MLL3 ChIP-Seq results (Figure 2C). Interestingly, the loci in cluster 1, which displayed the most substantial changes in histone modifications, involved genes that showed enriched MLL3 binding in Myc; sgTrp53 cells compared to the Myc; sgKmt2c genotype. These data support a model whereby MLL3 binding to these loci facilitates the acquisition of a chromatin environment conducive for active gene transcription.
MLL3 regulates specific tumor suppression programs in liver cancer cells
Transcriptional profiling helped hone in on potentially critical targets of MLL3. Specifically, we determined the output of these chromatin landscape changes by transcriptional profiling of the same set of Myc; sgTrp53 and Myc; sgKmt2c liver cancer cell lines described above. Despite the broad binding of MLL3 across the genome, we found only 248 differentially expressed genes (DEGs): 132 significantly upregulated (p<0.05, log2 fold-change >2) and 116 significantly downregulated (p<0.05, log2 fold-change <−2) in Myc; sgKmt2c liver tumor cells compared to Myc; sgTrp53 controls.
As predicted, transcripts encoding p53 and p53 target genes such as Ccng1, Cdkn1a, and Zmat3 (Bieging-Rolett et al., 2020) were upregulated in Myc; sgKmt2c cells, consistent with nonsense-mediated decay of truncated p53 transcripts and a concomitant reduction in p53 effector genes. Strikingly, some of the downregulated genes in Myc; sgKmt2c lines mapped to loci enriched in cluster 1, including Cdkn2a, Bmp6, and Lrp2 (Figure 2C–D, Figure 2—figure supplement 2D, E). Of note, sgKmt2c did not lead to compensatory changes in the transcript levels of other major components of the COMPASS-like complexes, including Mll4 (Kmt2d), Utx (Kdm6a), Mll1 (Kmt2a), and Mll2 (Kmt2b; Figure 3—figure supplement 1A), suggesting that the alterations in MLL3 binding, histone modification, and transcription were specifically attributed to MLL3 disruption.
We reason that the mediators of MLL3 actions in tumor suppression should be within cluster 1 with reduced transcription and MLL3 binding in Kmt2c-deficient cells. To further characterize the gene repertoire directly regulated by MLL3 genomic binding, we integrated the results of MLL3 ChIP-seq and RNA-seq. Specifically, we selected downregulated DEGs that show concordant decreased binding of MLL3 in Myc; sgKmt2c lines and subjected them to gene ontology analysis. Apart from the cluster 1 genes noted above, the integrative analysis revealed multiple MLL3-regulated tumor suppressive programs (Figure 3A, Figure 3—figure supplement 1B), including both cell-autonomous mechanisms (cellular metabolism) and non-autonomous mechanisms (interaction with extracellular matrix and immune system).
KMT2C and CDKN2A mutations result in similar transcriptomes in human HCC
One genomic locus that stood out in our integrative analysis was Cdkn2a, which encompasses both the p16/Ink4a and p19/Arf (p14 in human) tumor suppressors (Gil and Peters, 2006). CDKN2A is located on the human chromosome 9p and is deleted or epigenetically silenced in many cancer types (Sherr, 2012), including HCC (2017). While MLL3 likely regulates a plethora of genes that contribute to its tumor-suppressive potential, the well-defined and potent antitumor functions of Cdkn2a-encoded proteins make them attractive candidates as functionally relevant MLL3 effectors.
Furthermore, in our analysis publicly available genomic data on 1280 HCC samples, we found that CDKN2A alterations, like KMT2C alterations, showed significant co-occurrence with MYC gains and amplifications. However, we were unable to conduct a meaningful test of mutual exclusivity between CDKN2A and KMT2C alterations (Figure 3—figure supplement 2A), given the constraints of sample size and the modest frequencies of alteration in each gene. Further dissection of transcriptional profiling datasets from human and mouse HCCs harboring known gene alterations using gene set enrichment analysis (GSEA) revealed that human tumors with CDKN2A deletions transcriptionally resembled both mouse and human HCC harboring KMT2C alterations (Figure 3B and C) but not those harboring RB1 loss (Figure 3—figure supplement 2B), even though the tumor suppressor RB1 is regulated by CDKN2A/p16INK4A and their genomic alterations exhibit mutual exclusivity in multiple cancer types (Knudsen et al., 2020). While we cannot rule out the possibility that other factors drive these associations, our results support a biologically meaningful relationship between MLL3 and CDKN2A.
Cdkn2a locus is a genomic and transcriptional target of MLL3 in liver cancer
To explore the relationship between MLL3 and Cdkn2a locus in more detail, we tested whether genes encoded by Cdkn2a were direct targets of MLL3-regulated transcription. Indeed, Cdkn2a is a cluster 1 locus that, in Myc; sgKmt2c cancer cells, displays significant reduction in (1) expression, (2) H3K4me1/3 and H3K27ac levels, and (3) MLL3 binding at the Cdkn2a promoter compared with Myc; sgTrp53 cells (Figure 2C–D, Figure 4A). Of note, MLL3 binding peaks were also observed within the gene body of Cdkn2a. The differential expression of Ink4a and Arf was confirmed by qPCR, immunoblotting, and ChIP-qPCR analyses on multiple Myc; sgKmt2c and Myc; sgTrp53 liver cancer lines (Figure 4B, Figure 4—figure supplement 1A, B). These results imply that Cdkn2a locus is a genomic and transcriptional target of MLL3 in liver cancer cells.
Since the Myc; sgTrp53 and Myc; sgKmt2c cells we studied above are not isogenic, we performed a series of additional experiments to demonstrate a direct transcriptional effect of MLL3 on the Cdkn2a locus. Because p53 inactivation can lead to compensatory increases in Ink4a and Arf expression (Stott et al., 1998), representing an alternative possibility accounting for the observed difference of Cdkn2a expression in sgTrp53 vs sgKmt2c cells. However, p53 suppression in sgKmt2c cells produced only a subtle and inconsistent effect on the expression of Ink4a and Arf, whereas Kmt2c suppression in sgTrp53 cells consistently attenuated p16Ink4a and p19Arf protein levels (Figure 4—figure supplement 1C, D). As another means of ruling out the p53 pathway as an explanation for altered Cdkn2a expression in Myc; sgKmt2c cells, we tested the ability of MLL3 to regulate Cdkn2a transcripts in an orthogonal liver cancer model driven by Myc and inactivation of Axin1, which is a well-defined tumor suppressor that negatively regulates β-catenin activity in HCC (Satoh et al., 2000). Liver cancer cells produced by hydrodynamic delivery of the Myc transposon vector and Axin1 sgRNAs displayed reduced Ink4a and Arf expression upon Kmt2c knockdown without targeting p53 (Figure 4—figure supplement 1E,F). Importantly, MLL3 binding peaks at the Cdkn2a locus were also detected in Myc; sgAxin1 liver cancer cells (Figure 4—figure supplement 1G), suggesting that Cdkn2a transcription is directly regulated by MLL3 rather than an indirect outcome of p53 loss. These data imply that MLL3 supports a chromatin environment at the Cdkn2a locus that facilitates the transcription of both Ink4a and Arf and raises the possibility that these factors contribute to the tumor suppressor activity of MLL3 in liver cancer.
We next set out to determine whether MLL3 binding is sufficient to induce transcriptional activation of the CDKN2A locus and, in doing so, extend our analysis to human liver cancer cells. As the KMT2C transcript is too large (14,733 bp) for cDNA transduction, we turned to the CRISPR activation (CRISPRa) system (Chavez et al., 2015) in a human hepatocellular carcinoma cell line (HLE). Following stable integration of the nuclease dead Cas9 fused to the VP64-p65-Rta (VPR) transcriptional activator, cells were transduced with two orthogonal sgRNAs targeting the human KMT2C promoter (or, as control, transduced with sgRNA against GFP; Figure 4C). Cells expressing the KMT2C sgRNAs showed a marked and specific increase in the expression of endogenous KMT2C, but not of KMT2D or TP53 (Figure 4D, Figure 5—figure supplement 1A, B), which was accompanied by an increase in MLL3 binding to the CDKN2A locus (Figure 4E) and transcriptional upregulation of both CDKN2A transcripts (Figure 4F and G). Therefore, MLL3 directly binds and co-activates transcription of the CDKN2A locus in human liver cancer cells.
MLL3 mediates oncogene-induced apoptosis in a Cdkn2a-dependent manner
The above results raise the possibility that the Cdkn2a products, INK4A and ARF, may contribute to the tumor suppressive activity of MLL3. In this regard, Myc overexpression in primary cells (mouse embryonic fibroblasts; MEFs) often triggers apoptosis (Evan et al., 1992), and this in turn limits tumorigenesis in a manner that is dependent on Cdkn2a (Zindy et al., 1998). This pathway also suppresses liver tumorigenesis since concomitant disruption of Ink4a and Arf using CRISPR, or germline deletion of Arf alone, cooperated with Myc overexpression to rapidly promote tumor development (Figure 5—figure supplement 1C). Similarly, Kmt2c suppression also attenuated MYC-induced apoptosis, as shown by tumor histology and apoptosis by TUNEL assay (Negoescu et al., 1997), 5 days after hydrodynamic delivery of transposon vectors encoding Myc together with GFP-linked shRNAs targeting Kmt2c (or Renilla luciferase as a control; Figure 5A and B). This difference in apoptosis correlated with an increase in retention of GFP-shKmt2c expressing cells 10 days after injection (Figure 5—figure supplement 1D, E). Altogether, these results show that Kmt2c suppression impairs Myc-induced apoptosis in vivo in a manner that is reminiscent of the anti-apoptotic effects of Cdkn2a loss in the context of aberrant Myc activation (Eischen et al., 1999; Jacobs et al., 1999; Schmitt et al., 1999).
To model the interaction between Myc overexpression, MLL3 function, and Cdkn2a regulation, we transduced liver progenitor cells (LPCs) with retroviral vectors encoding Myc linked to a reverse tetracycline transactivator (rtTA3), together with doxycycline (dox)-inducible Kmt2c shRNAs to enable reversible Kmt2c silencing (Figure 5—figure supplement 2A). Infection of LPCs with Myc in the presence of MLL3 (i.e. cells infected with Myc-rtTA3 and a dox-inducible shRNA targeting Renilla luciferase) acutely activated INK4A and ARF expression (Figure 5—figure supplement 2B), and these cells could not be maintained in culture. Phenocopying the ability of Myc and Kmt2c suppression to transform liver cells in vivo, combined Myc and shKmt2c expression facilitated the persistent growth of cells maintained on Dox (Figure 5—figure supplement 2C,D). By contrast, dox withdrawal induced Kmt2c mRNA expression and H3K4me3 deposition at the Cdkn2a promoters, ultimately leading to elevations in Arf and Ink4a mRNA and protein (Figure 5C and D), reduced colony formation, and increased apoptosis (Figure 5—figure supplement 2D–F). Furthermore, constitutive shRNA-mediated knockdown of Arf and Ink4a through targeting of the shared exon 2 (shCdkn2a) significantly rescued colony-forming capacity and prevented cell death following Kmt2c restoration as determined by time-lapse microscopy of cells cultured with a fluorescent dye that stains dead cells (Figure 5E–G, Figure 5—figure supplement 2G). These data support a model whereby a prominent tumor suppressive output of MLL3 in liver cancer involves direct upregulation of Cdkn2a that, when impaired, attenuates the MYC-induced apoptotic program and permits tumor progression.
Discussion
Our study combined genetic, epigenomic, and animal modeling approaches to identify Cdkn2a as an important regulatory target of MLL3 in both mouse and human liver cancers. Our results support a model whereby oncogenic stress, herein produced by MYC, leads to an increase in the binding of MLL3 to the CDKN2A locus, an event that is associated with the accumulation of histone marks linked to the biochemical activity of MLL3-containing complexes and conducive to gene activation (Figure 6). Accordingly, these events are accompanied by transcriptional upregulation of two key Cdkn2a gene products, Ink4a and Arf. Moreover, suppression of Kmt2c phenocopies the effects of Cdkn2a inactivation in abrogating MYC-induced apoptosis. Conversely, suppression of Cdkn2a diminishes the anti-proliferative effects of Kmt2c restoration. As such, our results establish a conserved epistatic relationship between the chromatin modifier MLL3 and a well-characterized tumor suppressor network.
The epistatic relationship described above might be expected to lead to mutual exclusivity of KMT2C and CDKN2A alterations; however, we did not observe significant mutual exclusivity in liver cancer, which is likely due to insufficient samples sizes needed to obtain statistical power. Alternatively, other functionally important components linked to the CDKN2A locus could produce CDKN2A-independent forces that drive selection for chromosome 9p deletions, including type I interferon genes, CDKN2B, and MTAP (Barriga et al., 2022). Alternatively, mutual exclusivity between KMT2C and CDKN2A alterations would be expected only under circumstances where CDKN2A action is the most dominant MLL3 effector. Indeed, it seems likely that multiple downstream genes, including factors involved in interactions with stromal and immune populations, contribute to MLL3-driven tumor suppression, and their relative importance may vary between cell and tissue types. Such a variable output in cancer-relevant gene regulation has been noted for other chromatin regulators that, at the extreme, serve as pro-oncogenic factors in some contexts and tumor suppressors in others (Fountain et al., 1992; Schmid et al., 2000; Sun et al., 2017; Xia et al., 2021). Furthermore, our observation that Kmt2c deficiency cooperated with MYC but not CTNNB1 to drive HCC highlights such context specificity and is in line with recent findings that chromatin context could favor particular oncogenic alterations over others (Weiss et al., 2022).
UTX (KDM6A), MLL3 (KMT2C), and MLL4 (KMT2D), the core catalytic components of the COMPASS-like complex, are all considered tumor suppressors, with frequent loss-of-function genomic alterations found in a broad spectrum of human cancers (Revia et al., 2022; Sze and Shilatifard, 2016). While each of these components regulates redundant sets of genes (Hu et al., 2013; Lee et al., 2009), they may exert their tumor suppressive functions through different mechanisms. In liver and pancreas cancer models, UTX can control the expression of negative regulators of mTOR such as DEPTOR, and its disruption prevents their transcription and facilitates tumorigenesis through increased mTORC1 activity (Revia et al., 2022). Additionally, while the mechanisms of MLL4 activity have not been examined in liver cancer, studies suggest that MLL4 suppresses skin carcinogenesis by promoting lineage stability and ferroptosis independently of MLL3 (Egolf et al., 2021). Our study demonstrates that MLL3 is both necessary and sufficient for efficient transcriptional activation of the CDKN2A locus that drives oncogene-induced apoptosis. The molecular basis for this heterogeneity in effector output remains to be determined, but it seems likely that different subsets of target genes are preferentially disabled by haploinsufficiency of individual components and/or subject to compensation by the remaining COMPASS complex activities. Systematic studies comparing the binding, histone modifications, and transcriptional output of cells across a spectrum of allelic configurations of COMPASS complex factors will be needed to achieve a more holistic understanding of their functions and interactions in different contexts.
The most well-established role for MLL3/4-UTX-containing complexes is the control of H3K4 monomethylation at enhancers during development (Herz et al., 2010; Hu et al., 2013). While our ChIP-Seq studies also revealed binding of MLL3/4 to enhancers in liver tumor cells, an even larger fraction of genes—including Cdkn2a—showed MLL3/4 chromatin enrichment at gene promoters, and indeed, transcription of this class of genes was most affected by Kmt2c disruption. Interestingly, Kmt2c suppression preferentially limited the MLL3/4 enrichment at promoters and shifted residual complex binding toward intergenic regions. Such dynamic regulation of distinct cis-acting elements by the MLL3/4 complex has also been observed in other contexts (Cheng et al., 2014; Soto-Feliciano et al., 2023), where the non-canonical binding of MLL3/4 at promoters is a recurrent tumor suppressive mechanism in cancer cells. MLL3/4 has also been observed to bind within the exons and introns, which may enable chromatin looping of enhancers to activate gene expression (Panigrahi and O’Malley, 2021). Further studies into the action and regulation of MLL3/4 complexes at promoters and gene bodies will be informative and may yield new insights into the actions of the COMPASS-like complex in cancer.
While CDKN2A showed a surprisingly dominant role in mediating the tumor-suppressive effects of the broadly acting MLL3 enzyme, there are precedents for a predominant contribution of a single gene to the functional output of chromatin-complex disruption. Indeed, polycomb repressive complexes (PRCs) broadly repress gene expression in different cell types through the coordinated action of PRC1 and PRC2 complexes that deposit and maintain repressive H3K27me3 marks on the enhancers of target genes, including CDKN2A (Bracken et al., 2007; Kotake et al., 2007). Despite these similarly broad effects, CDKN2A is often the most functionally relevant target of PRC-mediated repression, as genetic deletion of either the PRC1 component Bmi1 or the PRC2 component Ezh2, or treatment with small molecule inhibitors of EZH2, can facilitate Cdkn2a induction in normal and tumor cells. This, in turn, triggers anti-proliferative responses that can be rescued by Cdkn2a deletion (Jacobs et al., 1999; Richly et al., 2011). Notably, the COMPASS-like complexes are biochemically and functionally similar to Trithorax complexes in Drosophila, which have an evolutionarily conserved antagonistic relationship with PRC1 and PRC2 that controls epigenetic memory and cell fate during development (Mills, 2010; Piunti and Shilatifard, 2016). Our findings suggest such antagonism extends to tumor suppression in mammalian cells, likely via regulation of Cdkn2a and other tumor suppressor genes (Soto-Feliciano et al., 2023).
Materials availability statement
Source files of all original gels and western blots were provided for the following figures:
RNA sequencing and ChIP-Seq data files that support the findings of this study have been deposited in the Gene Expression Omnibus under the accession code GSE85055, as well as in the Dryad digital repository (doi:10.5061/dryad.7pvmcvdwm; doi:10.5061/dryad.f1vhhmh0h). Sequences of sgRNAs, shRNAs, and primers used in this manuscript are included in the Supplementary file 1.
Materials and methods
Animal experiments
Request a detailed protocol8- to 10-week-old female C57BL/6 animals were purchased from Envigo (formerly Harlan). Each experiment was performed in mice from the same order. Arf-null animals (C57BL/6 background), originally provided by Dr. Charles Sherr, St. Jude Children’s Research Hospital, were maintained in our breeding colony. For HTVI, a sterile 0.9% NaCl solution/plasmid mix was prepared containing oncogene transposons (5 µg DNA of pT3-Myc or 10 µg pT3-Ctnnb1 N90) with either 20 µg of pX330 expressing the indicated sgRNAs or 20 µg of pT3-EF1a-GFP-miRE plasmid together with CMV-SB13 Transposase (1:5 ratio). Mice were randomly assigned to experimental groups and injected with the 0.9% NaCl solution/plasmid mix into the lateral tail vein with a total volume corresponding to 10% of body weight in 5–7 s as described before (Largaespada, 2009; Moon et al., 2019; Tschaharganeh et al., 2014; Xue et al., 2014). Injected mice were monitored for tumor formation by abdominal palpation. All animal experiments were approved by the Memorial Sloan Kettering Cancer Center (MSK) Institutional Animal Care and Use Committee (protocol 11-06-011). Animals were monitored for signs of ill health by veterinary staff at the Research Animal Resource Center at MSK, and efforts were made to minimize suffering.
Vector constructs
Request a detailed protocolThe pT3-Myc vector Addgene (#92046) and pT3-EF1a-GFP-miRE plasmid were described before (Huang et al., 2014). The pT3-Ctnnb1 N90 vector (Tward et al., 2007) was obtained from Addgene (#31785). For CRISPR/Cas9-mediated genome editing, sgRNAs were subcloned into pX330 (Addgene, #42230; Hsu et al., 2013). All shRNA and sgRNA sequences are listed in Supplementary file 1.
Derivation of primary liver tumor cell lines
Request a detailed protocolLiver tumors were resected with sterile instruments, and 10–50 mg of tumor tissue was minced and washed in sterile PBS, incubated in a mix of 1 mg/mL collagenase IV and 3 mg/mL dispase (dissolved in sterile, serum-free DMEM(Dulbecco's Modified Eagle Medium)) with gentle shaking, washed with PBS, incubated for 5 min in 0.05% (w/v) trypsin, and washed and plated in complete DMEM (10% FBS(fetal bovine serum), 1× penicillin/streptomycin) on collagen-coated plates (PurCol, Advanced Biomatrix). Primary cultures were passaged until visibly free from fibroblasts. Cell lines were authenticated on a routine basis using short tandem repeat profiling, as well as tested for mycoplasma contamination and immediately discarded upon a positive test.
Analysis of CRISPR-directed mutations
Request a detailed protocolCRISPR-mediated insertions and deletions were detected by surveyor assay as directed by the manufacturer (Transgenomic/IDT). Briefly, after overnight lysis of primary tumors and cell lines at 37°C in buffer containing 0.4 mg/mL proteinase K, 10 mM Tris, 100 mM NaCl, 10 mM EDTA, and 0.5% SDS, pH 8.0, genomic DNA was extracted by isopropanol precipitation. ~250–500 bp regions flanking predicted CRISPR cleavage sites were PCR amplified with Herculase II taq polymerase, column purified (Qiagen), heated to 95°C, and slowly cooled to promote annealing of heteroduplexes. Following treatment with Surveyor nuclease, products were analyzed by electrophoresis on a 2% polyacrylamide gel. Primers used for surveyor assay are listed in Supplementary file 1. Amplified PCR products were separately gel purified and ligated into blunt-end digested pBlueScript (Stratagene). DNA from 48 transformed colonies was analyzed by Sanger sequencing using a T7 primer.
CRISPR activation
Request a detailed protocolHuman HCC cell line HLE, purchased from JCRB Cell Bank (JCRB0404), was transduced by the lentivirus expressing nuclease-dead Cas9 (dCas9) fused with VPR (Chavez et al., 2015) and sgRNAs against KMT2C (sequence in Supplementary file 1) to generate stable MLL3 CRISPRa HLE line by puromycin selection.
Generation and modification of primary cells
Request a detailed protocolLPCs from E13.5–15.5 C57BL/6 embryos were isolated and grown in hepatocyte growth media (HGM) as previously described (Zender et al., 2005). To simultaneously overexpress Myc and conditionally suppress Kmt2c, LPCs were co-infected with a retroviral construct constitutively expressing both Myc and a reverse tet-transactivator (rtTA) (MSCV-Myc-IRES-rtTA) along with retroviral TRMPV vectors (MSCV-TRE-dsRed-miR30/shRNA-PGK-Venus-IRES-NeoR) (Zuber et al., 2011) expressing ds-Red linked, teint-responsive shRNAs targeting Kmt2c cloned into an optimized mir-30 context (‘mir-E,’ TRPMVe; Zuber et al., 2011). For selection of infected cells and sustained shKmt2c expression, cells were maintained in HGM with neomycin (1 mg/mL) and dox (1 µg/mL) starting 2 days after infection. To introduce constitutively expressed shRNAs in the setting of inducible shKmt2c, retroviral MLPe vectors (MSCV-LTR-mir-E-PGK-Puro-IRES-GFP; Dickins et al., 2005). GFP-linked shRNAs targeting either Cdkn2a or Renilla luciferase (as control) were co-infected with MSCV-Myc-IRES-rtTa and TRMPVe-shKmt2c. Triple-infected cells were maintained in media with neomycin, puromycin (2 µg/mL), and dox 2 days post infection. Infected continuously proliferating cells were transitioned to growth in complete DMEM and maintained on collagen-coated plates.
Colony assays
Request a detailed protocolFor measurement of cell proliferation, 5000 transduced and selected LPCs or MEFs were plated in triplicate in 6-well plates. Tetracycline-inducible shKmt2c–expressing LPCs were grown in the presence or absence of dox, and after 5 days, cells were fixed with formalin and methanol and stained with 0.05% crystal violet. MEFs were fixed after 6 days with formalin and methanol and stained with 0.05% crystal violet.
Apoptosis assays
Request a detailed protocolApoptosis was measured in LPCs via Annexin V staining according to the manufacturer’s instructions (eBiosciences, Annexin-V APC). 25,000 cells were grown with and without dox for 3 days, trypsinized, washed with Annexin-V binding buffer, and ~100,000 cells were incubated with Annexin-V APC and analyzed on an LSRII flow cytometer (BD).
Live imaging
Request a detailed protocolImaging was performed on LPCs immortalized by linked overexpression of Myc and two independent, inducible Kmt2c shRNAs constitutively expressing shRNAs targeting Renilla luciferase or Cdkn2a (generated as detailed above). 1000 cells were plated on collagen-coated, 96 well, clear bottom imaging plates in media supplemented with 300 nM Draq7 (Invitrogen) with and without dox, in triplicate by genotype. 18 hr after plating cells, Venus (marking all plated cells) and Draq7 fluorescence was collected in two, 10× fields of each well every 15 min for 41 hr using an automated, high content microscope (InCell 6000, General Electric).
Chromatin immunoprecipitation
Request a detailed protocolHistone ChIP was performed as previously described (Lee et al., 2006). Briefly, cell samples were cross-linked in 1% formaldehyde for 10 min, and the reaction was stopped by addition of glycine to 125 mM final concentration. Fixed cells were lysed in SDS lysis buffer, and the chromatin was fragmented by sonication (Covaris). Sheared chromatin was incubated with antibodies (final concentration 10 µg/mL) against H3K4me3 (Abcam, ab8580, Lot:GR164706-1), H3K27ac (Abcam, ab4729, Lot:GR200563-1), or H3K4me1 (Abcam; ab8895, Lot:GR114265-2) or with normal rabbit IgG (Abcam, ab46540) at 4°C for overnight. Antibodies were recovered by binding to protein A/G agarose (Millipore), and the eluted DNA fragments were used directly for qPCR or subjected to high-throughput sequencing (ChIP-Seq) using a HiSeq 2000 platform (Illumina). High-throughput reads were aligned to mouse genome assembly NCBI37/mm9 as previously described (Barradas et al., 2009). Reads that aligned to multiple loci in the mouse genome were discarded. The ChIP-Seq signal for each gene was quantified as total number of reads per million in the region 2 kb upstream to 2 kb downstream of the transcription start site (TSS). Primers used for ChIP-qPCR of mouse Cdkn2a promoter (Barradas et al., 2009) are listed in Table S1.
The complete dataset is available at NCBI Gene Expression Omnibus (GSE85055), as well as the Dryad digital repository (doi:10.5061/dryad.7pvmcvdwm).
For the MLL3 ChIP-Seq, the following protocol was used. Cross-linking ChIP in mouse and human HCC cells was performed with 10–20×107 cells per immunoprecipitation. Cells were collected, washed once with ice-cold PBS, and flash-frozen. Cells were resuspended in ice-cold PBS and cross-linked using 1% paraformaldehyde (PFA; Electron Microscopy Sciences) for 5 min at room temperature with gentle rotation. Unreacted PFA was quenched with glycine (final concentration 125 mM) for 5 min at room temperature with gentle rotation. Cells were washed once with ice-cold PBS and pelleted by centrifugation (800 g for 5 min). To obtain a soluble chromatin extract, cells were resuspended in 1 mL of LB1 (50 mM HEPES pH 7.5, 140 mM NaCl, 1 mM EDTA, 10% glycerol, 0.5% NP-40, 0.25% Triton X-100, and 1× complete protease inhibitor cocktail) and incubated at 4°C for 10 min while rotating. Samples were centrifuged (1400 g for 5 min), resuspended in 1 mL of LB2 (10 mM Tris-HCl pH 8.0, 200 mM NaCl, 1 mM EDTA, 0.5 mM EGTA, and 1× complete protease inhibitor cocktail), and incubated at 4°C for 10 min while rotating. Finally, samples were centrifuged (1400 g for 5 min) and resuspended in 1 mL of LB3 (10 mM Tris-HCl pH 8.0, 100 mM NaCl, 1 mM EDTA, 0.5 mM EGTA, 0.1% sodium deoxycholate, 0.5% N-lauroylsarcosine, and 1× complete protease inhibitor cocktail). Samples were homogenized by passing seven to eight times through a 28-gauge needle, then Triton X-100 was added to a final concentration of 1%. Chromatin extracts were sonicated for 14 min using a Covaris E220 focused ultrasonicator. Lysates were centrifuged at 20,000 g for 10 min at 4°C, and 5% of the supernatant was saved as input DNA. Beads were prepared by incubating them in 0.5% BSA in PBS and antibodies overnight (100 μL of Dynabeads Protein A or Protein G [Invitrogen] plus 20 μL of antibody). The antibody was anti-MLL3/4, which was kindly provided by the Wysocka laboratory (Dorighi et al., 2017). Antibody-beads mixes were washed with 0.5% BSA in PBS and then added to the lysates overnight while rotating at 4°C. Beads were then washed six times with RIPA buffer (50 mM HEPES pH 7.5, 500 mM LiCl, 1 mM EDTA, 0.7% sodium-deoxycholate, and 1% NP-40) and once with TE-NaCl Buffer (10 mM Tris-HCl pH 8.0, 50 mM NaCl, and 1 mM EDTA). Chromatin was eluted from beads in Elution buffer (50 mM Tris-HCl pH 8.0, 10 mM EDTA, and 1% SDS) by incubating at 65°C for 30 min while shaking, supernatant was removed by centrifugation, and crosslinking was reversed by further incubating chromatin overnight at 65°C. The eluted chromatin was then treated with RNase A (10 mg/mL) for 1 hr at 37°C and with proteinase K (Roche) for 2 hr at 55°C. DNA was purified by using phenol-chloroform extraction followed with ethanol precipitation. The NEBNext Ultra II DNA Library Prep kit was used to prepare samples for sequencing on an Illumina NextSeq500 (75 bp read length, single-end, or 37 bp read length, and paired-end).
The complete dataset for MLL3 ChIP-Seq is available at the Dryad digital repository (doi:10.5061/dryad.f1vhhmh0h).
Immunoblotting
Request a detailed protocolCell pellets were lysed in Laemmli buffer (100 mM Tris-HCl pH 6.8, 5% glycerol, 2% SDS, and 5% 2-mercaptoethanol). Equal amounts of protein were separated on 12% SDS–polyacrylamide gels and transferred to PVDF(polyvinylidene difluoride) membranes (90 V, 75 min). β-actin was used as a control to ensure equal loading, and images were analyzed using the AlphaView software (ProteinSimple). Immunoblotting was performed using antibodies for MYC (1:1000, Abcam, ab32072), p53 (1:500, Leica Biosystems, NCL-p53-505), p19 (1:250, Santa Cruz Biotechnology, sc-32748), p16 (1:250, Santa Cruz Biotechnology, sc-1207), Axin1 (1:1000, Cell Signaling, #2074), and β-actin (1:10000, Sigma-Aldrich, clone AC-15). Source files of all western blots were provided for Figure 1—figure supplement 2, Figure 4—figure supplement 1, Figure 5—figure supplement 2.
Quantitative RT-PCR
Request a detailed protocolTotal RNA was isolated using RNeasy Mini Kit, QIAshredder Columns, and RNase-Free DNase Set (Qiagen). cDNA synthesis was performed using TaqMan Reverse Transcription Reagents (Thermo Fisher Scientific). Real-time PCR was carried out using Power SYBR Green Master Mix (Thermo Fisher Scientific) and the Life Technologies ViiA 7 machine. Transcript levels were normalized to the levels of mouse or human Actb mRNA expression and calculated using the ΔΔCt method. Each qRT-PCR was performed in triplicate using gene-specific primers (sequences listed in Table S1).
RNA sequencing and differential expression analysis
Request a detailed protocolFor RNA sequencing, total RNA from three independent tumor-derived cell lines (Myc; sgTrp53 and Myc; sgKmt2c) was isolated using RNeasy Mini Kit, QIAshredder Columns and RNase-Free DNase Set (Qiagen). RNA-Seq library construction and sequencing were performed according to protocols used by the integrated genomics operation Core at MSK. 5–10 million reads were acquired per replicate sample. After removing adaptor sequences with Trimmomatic, RNA-seq reads were aligned to GRCm38.91(mm10) with STAR (Dobin et al., 2013). Genome-wide transcript counting was performed by HTSeq to generate an FPKM(Fragments Per Kilobase per Million mapped fragments) matrix (Anders et al., 2015). DEGs were identified by DESeq2 (v.1.8.2, package in R) and plotted in the volcano plot. The complete dataset is available at NCBI Gene Expression Omnibus (GSE85055).
Integrative analyses of RNA-seq and MLL3 ChIP-seq
Request a detailed protocolDifferential peaks from ChIP-Seq data were annotated by assigning all intragenic peaks to that gene while intergenic peaks were assigned using linear genomic distance to the TSS. Genes that were coordinately regulated (fold change >1.5 and adjusted p-value<0.1) in MLL3 ChIP-seq and RNA-seq data were selected for the integrated analysis. Enriched pathways were scored using the enrichGO function with ‘biological process’ in the clusterProfiler R package. Redundant pathways were collapsed using the ‘simplify’ function with a cutoff of 0.7 with the p.adjust metric. Network analysis was performed using differential peaks and genes by running enrichplot::cnetplot in R with default parameters.
Human cancer analyses
Request a detailed protocolRNA sequencing data of selected samples with somatic mutations or homozygous deletions of KMT2C, CDKN2A, TP53, or RB1 in the TCGA HCC dataset were downloaded from Broad Institute TCGA Genome Data Analysis Center. To obtain transcriptional signatures of HCC with genomic mutations and deletions of either KMT2C, CDKN2A, and RB1, differential gene expression analyses were performed by DESeq2 (with TP53-mutated HCCs as controls). The oncoprints of homozygous deletions and somatic mutations of KMT2C, CDKN2A, and TP53, as well as MYC gains and amplifications from human HCC datasets (Cancer Genome Atlas Research Network, 2017, MSK [Harding et al., 2019; Zheng et al., 2018], INSERM [Schulze et al., 2015], RIKEN [Fujimoto et al., 2012], AMC [Ahn et al., 2014], and MERCi [Ng et al., 2022]) were generated by cBioPortal (https://www.cbioportal.org ; Cerami et al., 2012; Gao et al., 2013).
Gene set enrichment analysis
Request a detailed protocolGSEA was performed using the GSEAPreranked tool for conducting GSEA of data derived from RNA-seq experiments (version 2.07) against other signatures. The metric scores (normalized enrichment scores and false discovery rate q-values) were calculated using the sign of the fold change multiplied by the inverse of the p-value (Subramanian et al., 2005). Specifically, transcriptional signatures were derived based on significantly changed genes (p-adjusted<0.05, absolute fold change >2) from RNA-seq of mouse HCC cell lines (Myc; sgKmt2c vs Myc; sgTrp53, n=3 each genotype), and in human HCCs with mutations in KMT2C vs TP53 (p-adjusted<0.05, absolute fold change >2). These signatures were compared to the transcriptional comparison of TCGA human HCCs with genomic inactivation of CDKN2A vs TP53.
Statistical analyses
Request a detailed protocolData are presented as mean ± D or SEM as specified. The statistical comparison between two groups was accomplished with the two-tailed student’s t-test or one-way ANOVA followed by post hoc t-tests among three or more groups. The analyses for co-occurrence or mutual exclusivity of mutations were performed using Fisher Exact test. Comparisons of survival curves were performed by log-rank tests. All statistical tests were performed using the Prism 8 software. All data presented in the manuscript have been replicated in independent cohorts of mice or in at least three biological replicates for in vitro experiments. On the basis of predicted effects of oncogene-tumor suppressor interaction introduced by HTVI in mice, with a power of 0.8 and p<0.05, we calculated a minimum sample size of 5 mice per group. Animals within the same cage were randomly allocated into control and experimental groups, with the group assignment recorded in a master spreadsheet and unmasked only when all samples of the respective experiments were analyzed. Data collection of each experiment was detailed in the respective figures, figure legends, and methods. No data were excluded from studies in this manuscript.
Appendix 1
Data availability
Source files of all original gels and Western Blots were provided for the following figures: Figure 1—figure supplement 2B; Figure 4—figure supplement 1A, C, D, E; Figure 5—figure supplement 2B, F, G. RNA sequencing and ChIP sequencing data files that support the findings of this study have been deposited in the Gene Expression Omnibus under the accession code GSE85055, as well as in the Dryad digital repository (doi:10.5061/dryad.7pvmcvdwm; doi:10.5061/dryad.f1vhhmh0h). Sequences of sgRNAs, shRNAs, and primers used in this manuscript are included in the Supplementary File 1.
-
Dryad Digital RepositoryMll3 suppresses tumorigenesis by activating the Ink4a/Arf locus.https://doi.org/10.5061/dryad.7pvmcvdwm
-
Dryad Digital RepositoryMLL3 ChIP sequencing in murine and human HCC cells.https://doi.org/10.5061/dryad.f1vhhmh0h
-
NCBI Gene Expression OmnibusID GSE85055. Mll3 suppresses tumorigenesis by activating the Ink4a/Arf locus.
References
-
Histone H3K4 Trimethylation by Mll3 as part of ASCOM complex is critical for NR activation of bile acid transporter genes and is downregulated in cholestasisAmerican Journal of Physiology. Gastrointestinal and Liver Physiology 300:G771–G781.https://doi.org/10.1152/ajpgi.00499.2010
-
Histone demethylase Jmjd3 contributes to epigenetic control of Ink4A/ARF by Oncogenic RASGenes & Development 23:1177–1182.https://doi.org/10.1101/gad.511109
-
Highly efficient Cas9-mediated transcriptional programmingNature Methods 12:326–328.https://doi.org/10.1038/nmeth.3312
-
An essential role for UTX in resolution and activation of Bivalent promotersNucleic Acids Research 44:3659–3674.https://doi.org/10.1093/nar/gkv1516
-
Probing tumor phenotypes using stable and regulated synthetic microRNA precursorsNature Genetics 37:1289–1295.https://doi.org/10.1038/ng1651
-
STAR: Ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Disruption of the ARF-Mdm2-P53 tumor Suppressor pathway in Myc-induced LymphomagenesisGenes & Development 13:2658–2669.https://doi.org/10.1101/gad.13.20.2658
-
Regulation of the Ink4B-ARF-Ink4A tumour Suppressor locus: all for one or one for allNature Reviews. Molecular Cell Biology 7:667–677.https://doi.org/10.1038/nrm1987
-
The H3K27Me3 demethylase dUTX is a Suppressor of Notch- and RB-dependent tumors in DrosophilaMolecular and Cellular Biology 30:2485–2497.https://doi.org/10.1128/MCB.01633-09
-
DNA targeting specificity of RNA-guided Cas9 nucleasesNature Biotechnology 31:827–832.https://doi.org/10.1038/nbt.2647
-
The Mll3/Mll4 branches of the COMPASS family function as major Histone H3K4 Monomethylases at enhancersMolecular and Cellular Biology 33:4745–4754.https://doi.org/10.1128/MCB.01181-13
-
Annual report to the nation on the status of cancer, 1975-2014Journal of the National Cancer Institute 109:djx030.https://doi.org/10.1093/jnci/djx030
-
Pan-cancer molecular analysis of the RB tumor Suppressor pathwayCommunications Biology 3:158.https://doi.org/10.1038/s42003-020-0873-9
-
Hepatocellular carcinomaNature Reviews. Disease Primers 7:6.https://doi.org/10.1038/s41572-020-00240-3
-
Tumor suppression by Ink4A-ARF: progress and puzzlesCurrent Opinion in Genetics & Development 13:77–83.https://doi.org/10.1016/s0959-437x(02)00013-8
-
Throwing the cancer switch: reciprocal roles of Polycomb and Trithorax proteinsNature Reviews. Cancer 10:669–682.https://doi.org/10.1038/nrc2931
-
TUNEL: Improvement and evaluation of the method for in situ apoptotic cell identificationBiochemica 2:12–17.
-
Roles of the Polycomb group proteins in stem cells and cancerCell Death & Disease 2:e204.https://doi.org/10.1038/cddis.2011.84
-
Ink4-ARF locus in cancer and agingWiley Interdisciplinary Reviews. Developmental Biology 1:731–741.https://doi.org/10.1002/wdev.40
-
Mll3/Mll4/COMPASS family on epigenetic regulation of enhancer function and cancerCold Spring Harbor Perspectives in Medicine 6:a026427.https://doi.org/10.1101/cshperspect.a026427
-
The Histone demethylase UTX enables RB-dependent cell fate controlGenes & Development 24:327–332.https://doi.org/10.1101/gad.1882610
-
Generation and analysis of genetically defined liver Carcinomas derived from Bipotential liver ProgenitorsCold Spring Harbor Symposia on Quantitative Biology 70:251–261.https://doi.org/10.1101/sqb.2005.70.059
Article and author information
Author details
Funding
National Cancer Institute (P01 CA013106)
- Scott W Lowe
National Cancer Institute (R01 CA233944)
- Scott W Lowe
National Institute of General Medical Sciences (1K99GM140265-01)
- Yadira M Soto-Feliciano
National Cancer Institute (1F32CA257103)
- Changyu Zhu
American Cancer Society (PF-14-066-01-TBE)
- John P Morris
Helmholtz foundation (VH-NG-1114)
- Darjus F Tschaharganeh
National Cancer Institute (P30 CA008748)
- Scott W Lowe
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Charles Sherr and Janet Novak for constructive guidance and advice on all aspects of this study. We thank Ali Shilatifard, Lu Wang, and all members of the Lowe lab for helpful and stimulating discussions. We gratefully thank A Chramiec for excellent technical assistance. We thank Joanna Wysocka (Stanford University) for kindly sharing the anti-MLL3/4 antibody used in our ChIP-Seq experiments. This work was supported by grants to SWL (P01 CA013106 and R01 CA233944) from the NIH/NCI, as well as by the National Center for Tumor Disease, Heidelberg, and grants of the German Research Foundation to DFG (SFB/TRR77). This work was also supported by the NIH/NCI Cancer Center Support Grant to Memorial Sloan Kettering Cancer Center (P30 CA008748). YMSF is supported by a MOSAIC K99/R00 Award from the NIH/NIGMS (1K99GM140265-01). CZ is supported by an F32 Postdoctoral Fellowship (1F32CA257103) from the NIH/NCI. JPM was a recipient of a Postdoctoral Fellowship (PF-14-066-01-TBE) from the American Cancer Society. DFT is supported by a Young Investigator Group (VH-NG-1114) by the Helmholtz foundation. SWL is the Geoffrey Beene Chair for Cancer Biology and an investigator of the Howard Hughes Medical Institute.
Ethics
All animal experiments were approved by the MSKCC Institutional Animal Care and Use Committee (protocol 11-06-011). Animals were monitored for signs of ill-health by veterinary staff at the Research Animal Resource Center (RARC) at MSKCC and efforts were made to minimize suffering.
Copyright
© 2023, Zhu, Soto-Feliciano, Morris et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,618
- views
-
- 223
- downloads
-
- 6
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
- Immunology and Inflammation
Despite major successes with inhibitory receptor blockade in cancer, the identification of novel inhibitory receptors as putative drug targets is needed due to lack of durable responses, therapy resistance, and side effects. Most inhibitory receptors signal via immunoreceptor tyrosine-based inhibitory motifs (ITIMs) and previous studies estimated that our genome contains over 1600 ITIM-bearing transmembrane proteins. However, testing and development of these candidates requires increased understanding of their expression patterns and likelihood to function as inhibitory receptor. Therefore, we designed a novel bioinformatics pipeline integrating machine learning-guided structural predictions and sequence-based likelihood models to identify putative inhibitory receptors. Using transcriptomics data of immune cells, we determined the expression of these novel inhibitory receptors, and classified them into previously proposed functional categories. Known and putative inhibitory receptors were expressed across different immune cell subsets with cell type-specific expression patterns. Furthermore, putative immune inhibitory receptors were differentially expressed in subsets of tumour infiltrating T cells. In conclusion, we present an inhibitory receptor pipeline that identifies 51 known and 390 novel human inhibitory receptors. This pipeline will support future drug target selection across diseases where therapeutic targeting of immune inhibitory receptors is warranted.
-
- Cancer Biology
Gastric cancer (GC) is a major cause of cancer-related mortality worldwide. Despite the widespread recognition of tumor immunotherapy in treating unresectable GC, challenges, including ineffective immunotherapy and drug resistance, persist. Therefore, understanding the regulatory mechanisms of PD-L1, particularly in the context of super-enhancers (SEs) and zinc finger protein 36 ring finger protein-like 1 (ZFP36L1) RNA-binding protein, is crucial. In this study, we performed H3K27ac Cleavage Under Targets and Tagmentation (CUT&Tag) sequencing, investigated the heterogeneity of SEs between two GC subtypes with differential growth patterns, and revealed the immune escape signatures driven by ZFP36L1-SE in infiltrative GC through SEs inhibitors treatment. The regulation of ZFP36L1 to PD-L1 was evaluated by quantitative PCR, western blot, flow cytometry, and immunohistochemistry. Furthermore, we explored its regulatory mechanisms using a combination of molecular biology techniques, including luciferase reporter assay, GST/RNA pull-down, chromatin immunoprecipitation (ChIP)/RIP experiments, and in vivo functional assays. We demonstrated that ZFP36L1, driven by an SE, enhances IFN-γ-induced PD-L1 expression, with SPI1 identified as the specific transcription factor binding to ZFP36L1-SE. Mechanistically, ZFP36L1 binds to the adenylate uridylate-rich element in the 3ʹ untranslated region (3ʹUTR) of HDAC3 mRNA, exacerbating its mRNA decay, and thereby facilitating PD-L1 abnormal transcriptional activation. Collectively, our findings provide mechanistic insights into the role of the SPI1-ZFP36L1-HDAC3-PD-L1 signaling axis in orchestrating immune escape mechanisms in GC, thereby offering valuable insights into the potential targets for immune checkpoint therapy in GC management.