HOXA9 promotes MYC-mediated leukemogenesis by maintaining gene expression for multiple anti-apoptotic pathways
Abstract
HOXA9 is often highly expressed in leukemias. However, its precise roles in leukemogenesis remain elusive. Here, we show that HOXA9 maintains gene expression for multiple anti-apoptotic pathways to promote leukemogenesis. In MLL fusion-mediated leukemia, MLL fusion directly activates the expression of MYC and HOXA9. Combined expression of MYC and HOXA9 induced leukemia, whereas single gene transduction of either did not, indicating a synergy between MYC and HOXA9. HOXA9 sustained expression of the genes implicated in the hematopoietic precursor identity when expressed in hematopoietic precursors, but did not reactivate it once silenced. Among the HOXA9 target genes, BCL2 and SOX4 synergistically induced leukemia with MYC. Not only BCL2, but also SOX4 suppressed apoptosis, indicating that multiple anti-apoptotic pathways underlie cooperative leukemogenesis by HOXA9 and MYC. These results demonstrate that HOXA9 is a crucial transcriptional maintenance factor that promotes MYC-mediated leukemogenesis, potentially explaining why HOXA9 is highly expressed in many leukemias.
Introduction
Mutations of transcriptional regulators often cause aberrant gene regulation of hematopoietic cells, which leads to leukemia. Structural alterations of the mixed lineage leukemia gene (KMT2A also known as MLL) by chromosomal translocations cause malignant leukemia that often associates with poor prognosis despite the current intensive treatment regimens (Tamai and Inokuchi, 2010). KMT2A encodes a transcriptional regulator termed MLL that maintains segment-specific expression of homeobox (HOX) genes during embryogenesis (Yu et al., 1998), which determines the positional identity within the body (Luo et al., 2019; Deschamps and van Nes, 2005; Wang et al., 2009). During hematopoiesis, MLL also maintains the expression of posterior HOXA genes and MEIS1 (another homeobox gene), which promote the expansion of hematopoietic stem cells and immature progenitors (Jude et al., 2007; Krivtsov et al., 2006; McMahon et al., 2007; Thorsteinsdottir et al., 2002; Yagi et al., 1998). The oncogenic MLL fusion protein constitutively activates its target genes by constitutively recruiting transcription initiation/elongation factors thereto (Yokoyama et al., 2010; Lin et al., 2010; Okuda et al., 2015). Consequently, HOXA9 and MEIS1 are highly transcribed in MLL fusion-mediated leukemia (Krivtsov et al., 2006). Forced expression of HOXA9 (but not MEIS1) immortalize hematopoietic progenitor cells (HPCs) ex vivo (Schnabel et al., 2000; Kroon et al., 1998). Co-expression of HOXA9 with MEIS1 causes leukemia in mice which recapitulates MLL fusion-mediated leukemia (Kroon et al., 1998). Moreover, overexpression of HOXA9 is observed in many non-MLL fusion-mediated leukemias such as those with NPM1 mutation and NUP98 fusion and is associated with poor prognosis (Collins and Hess, 2016). These findings highlight HOXA9 as a major contributing factor in leukemogenesis. Nevertheless, the mechanism by which HOXA9 promotes oncogenesis remains elusive.
HOXA9 is considered to function as a transcription factor, which retains a sequence-specific DNA binding ability. HOX proteins have an evolutionally conserved homeodomain which possesses strong sequence preferences (Berger et al., 2008). HOXA9 associates with other homeodomain proteins such as PBX and MEIS family proteins (Schnabel et al., 2000; Shen et al., 1999). HOXA9 and those HOXA9 cofactors form a stable complex on a DNA fragment harboring consensus sequences for each homeodomain protein (Shen et al., 1999; Chang et al., 1996), suggesting that they form a complex of different combinations in a locus-specific manner depending on the availability of the binding sites. Recently, it has been reported that HOXA9 specifically associates with enhancer apparatuses (e.g. MLL3/4) to regulate gene expression (Sun et al., 2018; Huang et al., 2012; Zhong et al., 2018). However, the mechanisms by which HOXA9 activates gene expression remain largely unclear.
In this study, we reveal the oncogenic roles for HOXA9 and its target gene products in leukemogenesis and its unique mode of function as a transcriptional maintenance factor that preserves an identity of a hematopoietic precursor.
Results
MLL fusion proteins and HOXA9 sustain MYC expression against differentiation-induced transcriptional suppression
To identify the direct target genes of MLL fusion proteins, we first examined the genome-wide localization pattern of MLL fusion proteins by chromatin immunoprecipitation (ChIP) followed by deep sequencing (ChIP-seq), using HB1119 cells, a cell line expressing the MLL-ENL fusion protein. We observed MLL ChIP signals on the MYC, HOXA9, HOXA10, and MEIS1 loci (Figure 1A; Okuda et al., 2017), which were further confirmed by ChIP-quantitative polymerase chain reaction (qPCR) analysis (Figure 1—figure supplement 1A). These ChIP signals can be mostly attributed to MLL-ENL as the knockdown of wild-type MLL did not affect the MLL ChIP signals (Figure 1A and Figure 1—figure supplement 1B,C; Okuda et al., 2017). The distribution of MLL-ENL was enriched at transcription start sites in a genome-wide manner, which is similar to that of wild-type MLL observed in non-MLL fusion cell lines including HEK293T and REH (Okuda et al., 2017; Miyamoto et al., 2020).
To characterize the dynamic changes in MLL target gene expression during differentiation, we isolated bone marrow cells from mice at various differentiation stages ranging from the most immature population (c-Kit+, Sca1+, Lineage–; KSL) containing hematopoietic stem cells to highly differentiated hematopoietic cells (c-Kitlow/Mac1high) by fluorescence-activated cell sorting (FACS) and performed quantitative reverse transcription-polymerase chain reaction (qRT-PCR) analysis. For comparison, we analyzed leukemia cells (LCs) harvested from mice suffering with MLL-AF10-induced leukemia (MLL-AF10-LCs) and immortalized cells (ICs) transformed by MLL-ENL ex vivo (MLL-ENL-ICs). In accordance with previous reports (Krivtsov et al., 2006; Somervaille and Cleary, 2006; Yokoyama et al., 2013), the expression levels of Hoxa9, Hoxa10, and Meis1 were downregulated during normal hematopoietic differentiation but remained high in MLL fusion-expressing cells (Figure 1B). Myc was highly expressed in KSL, common myeloid progenitors (CMPs), and granulocyte/macrophage progenitors (GMPs), which contain actively dividing populations (Passegué et al., 2005), but was completely suppressed at highly differentiated c-kitlow/Mac1high stages. In the two MLL fusion-expressing cell lines, Myc was expressed at comparable levels to those in the progenitor fractions (CMP and GMP; Figure 1B). Mxd1, a differentiation marker, was highly expressed in differentiated populations alone. These results indicate that Myc, Hoxa9, Hoxa10, and Meis1 are intrinsically programmed to be silenced in normal hematopoietic differentiation but are aberrantly maintained by MLL fusion proteins.
To assess the oncogenic potential of MLL target genes, we performed myeloid progenitor transformation assays (Lavau et al., 1997; Okuda and Yokoyama, 2017a), wherein HPCs were isolated from mice, retrovirally transduced with each MLL target gene, and cultured in a semi-solid medium supplemented with cytokines promoting myeloid lineage differentiation (Figure 1C,D). HPCs exogenously expressing MLL-ENL or MLL-AF10 produced a large number of colonies in the third and fourth passages with high mRNA levels of Myc, Hoxa9, Hoxa10, and Meis1 in colonies of the first passage, confirming their potent transforming capacities. These cells were considered ‘immortalized’ as they proliferate indefinitely in this ex vivo culture (Lavau et al., 1997; Okuda and Yokoyama, 2017a). Ectopic expression of MYC immortalized HPCs; therefore, MYC expression is sufficient to induce proliferation of HPCs. Endogenous Myc expression was severely diminished in MYC-ICs, demonstrating that Myc expression is programed to be silenced and cannot be sustained by MYC itself. Ectopic expression of HOXA9 and HOXA10 (but not MEIS1) immortalized HPCs. Endogenous Myc expression was maintained in HOXA9/A10-ICs, indicating that HOXA9/A10 can sustain Myc expression against differentiation-induced transcriptional suppression. It should be noted that the colony size of HOXA9/A10-ICs was relatively small compared with MYC- or MLL fusion-ICs, suggesting a weaker proliferative potential. Knockdown of Myc by shRNA completely repressed the colony-forming ability of MLL-ENL-ICs and HOXA9-ICs (Figure 1E and Figure 1—figure supplement 2). Taken together, these results demonstrate that MLL fusion proteins and HOXA9 maintain Myc expression despite of the differentiation-induced transcriptional suppression, and the maintenance of MYC expression is indispensable for immortalization of HPCs ex vivo.
HOXA9 confers the identity of a hematopoietic precursor while MYC drives anabolic pathways
To identify the genes specifically regulated by HOXA9 but not by MYC, we performed RNA-seq analysis of HOXA9-ICs and MYC-ICs which do not express HOXA9 (Figure 1C). Genes highly expressed in HOXA9-ICs but lowly expressed in MYC-ICs (defined as ‘HOXA9 high signature’) were associated with hematopoietic identity/functions. In contrast, genes highly expressed in MYC-ICs and lowly expressed in HOXA9-ICs (defined as ‘MYC high signature’) were associated with anabolic pathways involved in nucleotide/protein production (Figure 2A,B). Pathways involved in lipid metabolism were commonly upregulated in both HOXA9- and MYC-ICs compared with non-immortalized c-kit-positive HPCs (Figure 2—figure supplement 1A,B). MLL-AF10-ICs, which express endogenous Hoxa9 and Myc at high levels (Figure 1C), expressed both HOXA9 high and MYC high signature genes (Figure 2A–C), indicating that MLL-AF10-ICs possess MYC-mediated highly proliferative potential and HOXA9-mediated hematopoietic identity. These differences between the HOXA9 high and MYC high signatures indicate that HOXA9 maintains the identity of a hematopoietic precursor, while MYC promotes proliferation by upregulating anabolic pathways.
Apoptosis is induced by MYC, and alleviated by HOXA9 and MLL-AF10
Given that excessive MYC activity promotes apoptosis (McMahon, 2014), we evaluated the apoptotic tendencies in immortalized HPCs. MYC-ICs exhibited increased γH2AX, cleaved poly (ADP) ribose polymerase (PARP), and cleaved caspase 3 levels, indicative of a high degree of replication stress and apoptosis (Figure 3A). FACS analysis with Annexin V also showed that MYC-ICs had more significant apoptotic fraction than that of HOXA9-ICs (Figure 3B). MLL-AF10-ICs exhibited weak apoptotic tendencies similarly to HOXA9-ICs, although their MYC expression tended to be higher than HOXA9-ICs (Figures 1C and 3A). Furthermore, most MYC-ICs underwent massive apoptosis 1 day after cytokine removal, whereas HOXA9-ICs and MLL-AF10-ICs showed relative resistance (Figure 3C) and exhibited successful recovery of the live cell population after cytokine reintroduction (Figure 3—figure supplement 1). These results indicate that HOXA9 confers anti-apoptotic properties.
Next, we evaluated the leukemogenic potential of HOXA9 and MYC in vivo. Transplantation of HPCs exogenously expressing MLL-AF10 into syngeneic mice induced leukemia with full penetrance (Figure 3D). In contrast, neither MYC nor its homologue MYCN could initiate leukemia within 200 days under these experimental conditions, indicating that activation of the MYC high signature alone is insufficient for leukemogenesis in vivo. Although one recipient mouse of MYC-transduced HPCs became sick and was sacrificed approximately 150 days after transplantation, it did not exhibit any leukemia-associated signs. The mouse Myc gene also failed to induce leukemia. Thus, it is unlikely that the inability of MYC alone to induce leukemia under these experimental conditions is due to immune suppression of the cells expressing human transgenes. HOXA9 did not induce leukemia within 200 days either, suggesting that the HOXA9 high signature alone is also insufficient to induce leukemia. Taken together, these results suggest that both high MYC activity and HOXA9-mediated resistance to apoptosis are necessary for driving leukemogenesis in vivo.
HOXA9 promotes MYC-mediated leukemogenesis
We next examined the gene expression of patients with leukemia using publicly available microarray data of the Microarray Innovations in LEukemia (MILE) study (Haferlach et al., 2010). Of 108 cases with MLL translocation, 38 were acute myelogenous leukemia [AML], and 70 were acute lymphoblastic leukemia [ALL]. Most of the cases (78.7%) were categorized in the HOXA9high/MEIS1high group, while some were in the HOXA9low/MEIS1high (13.0%) or HOXA9high/MEIS1low (8.3%) group (Figure 4A). MYC was expressed at high levels irrespective of HOXA9 or MEIS1 expression. These data indicate variability in transcriptional profiles among MLL fusion-mediated leukemia cases, with MYC expression remaining consistently high. HOXA9low/MEIS1high leukemia was predominantly found in ALL, likely due to the MLL-AF4 cases, some of which do not express HOXA genes (Lin et al., 2016). HOXA9high/MEIS1low leukemia was mainly found in AML (Figure 4B).
Next, we assessed the effects of the combined expression of MYC, HOXA9, and MEIS1 in mouse leukemia models. The MYC/HOXA9 combination exhibited high clonogenicity ex vivo, with colony-forming capacity and colony size comparable to those of MLL-ENL- and HOXA9/MEIS1-ICs (Figure 4C,D). A weak synergy between MYC and MEIS1 was also observed. HOXA9/MEIS1-ICs showed high Myc expression (Figure 4—figure supplement 1A) and MYC-dependent proliferation (Figure 4—figure supplement 1B), indicating that MYC is essential in HOXA9/MEIS1-mediated leukemic transformation. In the in vivo leukemogenesis assays, MYC, HOXA9, or MEIS1 expression alone did not initiate leukemia, whereas co-expression of HOXA9 and MEIS1 did in all recipient mice as previously reported (Figure 4E; Kroon et al., 1998). The combined expression of MYC/HOXA9 or MYC/MEIS1 induced leukemia in 38% and 18% of recipient mice, respectively, indicating a synergy of MYC with HOXA9 and MEIS1. qRT-PCR analysis showed that LCs maintained the expression patterns of endogenous Hoxa9, Meis1, and Myc similar to those of the respective ICs (Figure 4—figure supplement 2). We observed rapid onset of leukemia with full penetrance in the secondary transplantation for the three combinations tested, confirming the presence of leukemia-initiating cells (Figure 4E). These results indicate a cooperative role for HOXA9 and MEIS1 in MYC-mediated leukemogenesis with stronger synergy between HOXA9 and MYC.
All the MYC/HOXA9- and HOXA9/MEIS1-induced leukemias were positive to a myeloid marker Mac1 but negative to lymphoid marker B220 or CD3e (Figure 4—figure supplement 3A,B), indicating that both MYC/HOXA9- and HOXA9/MEIS1-combinations induced myeloid leukemia. Gene set enrichment analysis (GSEA) of MYC/HOXA9-LCs and HOXA9/MEIS1-LCs revealed that genes involved in protein synthesis (e.g. ribosome biogenesis, tRNA metabolic process) were enriched in MYC/HOXA9-LCs (Figure 4—figure supplement 4A), while tumour necrosis factor α signaling and KRAS pathways were enriched in HOXA9/MEIS1-LCs (Figure 4—figure supplement 4B). These results suggested that these two combinations promoted leukemogenesis in slightly different ways, wherein MYC-mediated anabolic pathways played more prominent roles in leukemogenesis by MYC/HOXA9 than by HOXA9/MEIS1.
HOXA9 functions as a transcription maintenance factor
Some HOXA9 high signature genes, namely Bcl2, Sox4, and Igf1, have been implicated in leukemogenesis (Zhang et al., 2013; Du et al., 2005; Delbridge et al., 2016; Steger et al., 2015), suggesting that they may be responsible for the synergy between HOXA9 and MYC. Indeed, Bcl2, Sox4, and Igf1 were highly transcribed in HPCs transformed by the combinations containing HOXA9 (i.e. MYC/HOXA9-ICs, HOXA9/vector-ICs, and HOXA9/MEIS1-ICs) but not those lacking it (i.e. MYC/vector-ICs and MYC/MEIS1-ICs) (Figure 5A and Figure 5—figure supplement 1A). Conditional loss of function experiments using HOXA9 conjugated with estrogen receptor (HOXA9-ER) confirmed critical regulation of these genes by HOXA9 (Figure 5B). However, stepwise transduction of MYC followed by HOXA9 failed to upregulate these HOXA9 target genes, indicating that HOXA9 could not reactivate its target genes once silenced (Figure 5C and Figure 5—figure supplement 1B). Conditional inactivation and subsequent reactivation of HOXA9 within the same cell population confirmed that HOXA9 cannot reactivate its target genes once silenced (Figure 5—figure supplement 1C). RNA-seq analysis of various human cell lines showed that HOXA9 was highly expressed in all the MLL fusion cell lines tested (Figure 5—figure supplement 2A, highlighted in blue). BCL2 and SOX4 were also expressed but at different levels depending on the cell lines. In accord, most of the MLL fusion-mediated leukemia patient samples in the MILE study with high HOXA9 expression expressed BCL2 and SOX4 at high levels (Figure 5—figure supplement 2B; Haferlach et al., 2010). IGF1 was only expressed in MV4-11 cells despite comparable HOXA9 expression among the MLL fusion cell lines, supporting the notion that HOXA9 expression itself cannot trigger the expression of silenced HOXA9 target genes. These results indicate that HOXA9 is a transcriptional maintenance factor that may be involved in the maintenance of chromatin structure previously activated by other transcriptional/epigenetic factors.
BCL2 and SOX4 promote MYC-mediated leukemogenesis by alleviating apoptosis
To identify the roles for BCL2 and SOX4 in leukemic transformation, we evaluated the leukemogenic potential of the combined expression of BCL2 or SOX4 with MYC. In myeloid progenitor transformation assays, SOX4 by itself showed weak immortalization capacity as previously reported (Zhang et al., 2013), while BCL2 did not transform HPCs (Figure 6A). Co-expression of BCL2 or SOX4 with MYC led to a substantial increase in colony-forming capacity. The combined expression of BCL2 with MYC induced leukemia in vivo with a penetrance similar to that with the MYC/HOXA9 combination (Figures 4E and 6B, and Figure 6—figure supplement 1A), in accord with previous reports (Beverly and Varmus, 2009; Luo et al., 2005). SOX4 also promoted MYC-mediated leukemogenesis in vivo (Figure 6B and Figure 6—figure supplement 1A), while neither MYC, BCL2, nor SOX4 alone induced leukemia within 200 days. FACS analysis indicated that all of MYC/SOX4-induced leukemias were of myeloid lineage (Figure 6—figure supplement 1B). In contrast, all MYC/BCL2-induced leukemias were of lymphoid lineage, the half of which was B-cell type, and the other half was T-cell type, consistent with a previous report (Luo et al., 2005). It should be noted that single gene transduction of MYC and SOX4 induced leukemia by others in different settings, albeit with low penetrance (Du et al., 2005; Luo et al., 2005). These differences are possibly due to the differences of virus titers and viral genome integration-related gene activation. Furthermore, the combined expression of HOXA9, BCL2, or SOX4 with MYC alleviated apoptotic tendencies (Figure 6C,D). Although SOX4 is reported to modulate transcription of pro/anti-apoptotic genes (Ramezani-Rad et al., 2013), their expression was not drastically altered by SOX4 in this context (Figure 6E,F). Thus, the mechanism underlying the anti-apoptotic properties of SOX4 in this setting is currently unclear. Taken together, the results indicate that multiple anti-apoptotic pathways mediated by BCL2 and SOX4 promote MYC-mediated leukemic transformation.
Endogenous BCL2 and SOX4 support the initiation and maintenance of leukemia
To examine the roles for endogenous BCL2 and SOX4 in the initiation and maintenance of LCs, we conducted myeloid progenitor transformation and in vivo leukemogenesis assays using Bcl2- and Sox4-knockout HPCs. We transduced various oncogenes into HPCs isolated from fetal livers of Bcl2- and Sox4-knockout embryos (Kamada et al., 1995; Schilham et al., 1996) and cultured them ex vivo. Neither Bcl2 nor Sox4 deletion affected the immortalization of HPCs by HOXA9, MYC, or MLL-AF10, indicating that BCL2 and SOX4 are dispensable for proliferation ex vivo (Figure 7—figure supplement 1). On the other hand, Bcl2 deletion delayed the onset of leukemia induced by MLL-AF10 and the HOXA9/MEIS1 combination in vivo (Figure 7A and Figure 7—figure supplement 2A). Sox4 deletion also delayed the onset of HOXA9/MEIS1-induced leukemia, although it did not affect MLL-AF10-induced leukemia (Figure 7B and Figure 7—figure supplement 2A). We also analyzed the steady-state apoptosis of LCs harvested from moribund mice with apoptotic markers (cleaved caspase, γH2AX and Annexin V). HOXA9/MEIS1-LCs were slightly more apoptotic in the absence of Bcl2 and Sox4 (Figure 7—figure supplement 2B,C). However, we did not observe substantial differences in MLL-AF10-LCs between the WT and Bcl2/Sox4 KO genotypes. These results are consistent with the different kinetics of leukemia onset (Figure 7A,B), and the minor differences in apoptotic tendencies in full-blown leukemia cells suggest that these LCs had acquired adequate anti-apoptotic properties at the time of disease presentation. These results indicate that endogenous BCL2 and SOX4 partially contribute to initiating leukemia in vivo.
Next, we examined the roles for endogenous BCL2 and SOX4 in maintaining leukemia-initiating cells using a mouse leukemia cell line that we previously established with MLL-ENL (Okuda et al., 2017). CRISPR/Cas9-mediated sgRNA competition assays indicated the negligible contribution of Bcl2 and Sox4 to MLL-ENL-LC proliferation ex vivo (Figure 7—figure supplement 3). To test their roles in vivo, we transplanted Bcl2- or Sox4-deficient LCs into syngeneic mice, where we observed delayed onset and reduced incidence rate for both Bcl2- and Sox4-knockout LCs (Figure 7C). Forced expression of sgRNA-resistant cDNAs of BCL2 and SOX4 restored leukemogenic potential, validating the on-target effects of sgRNAs (Figure 7D). These results suggest that MLL-ENL-LCs are partially dependent on endogenous BCL2 and SOX4 in vivo. Taken together, our results indicate that MLL fusion-mediated LCs depend on the expression of multiple anti-apoptotic genes via HOXA9 to varying degrees to achieve survival advantages for disease initiation and maintenance (Figure 7E).
Discussion
In this study, we found that HOXA9 regulates a variety of genes to maintain hematopoietic precursor identity and its associated anti-apoptotic properties. In leukemic transformation, MLL fusion proteins exploit both HOXA9 and MYC downstream pathways. Accordingly, HOXA9 and MYC synergistically induce leukemia in mouse models. Thus, we propose that MLL fusion proteins employ two arms to promote oncogenesis: MYC-mediated proliferation and HOXA9-mediated resistance to differentiation/apoptosis.
It is widely accepted that two types of mutations need to occur before leukemia onset; class I mutations that confer proliferative advantages and class II mutations that block differentiation (Gilliland, 2002). However, it was unclear how MLL mutations fit this theory because MLL fusion-mediated leukemia does not require additional mutations in many cases (Andersson et al., 2015). A comparison of gene expression profiles of HOXA9- and MYC-transformed HPCs demonstrated that HOXA9 maintains the expression of a wide range of genes associated with hematopoietic precursor identity. Thus, HOXA9 appears to maintain the intrinsic transcriptional programs of immature HPCs, which are programed to be silenced during differentiation. On the other hand, MYC upregulates genes involved with de novo nucleotide/protein synthesis. This is in line with the known involvement of MYC in proliferation and associated anabolic processes (Ji et al., 2011). Furthermore, MLL-AF10 activated both the HOXA9 high and MYC high signature genes to induce leukemia while ectopic expression of HOXA9 and MYC synergistically induced leukemia. Thus, our findings suggest that MLL fusion proteins activate both HOXA9- and MYC-dependent programs as alternative mechanisms to the combination of class I and II mutations.
Although HOXA9 maintained MYC expression to immortalize myeloid progenitors ex vivo, it did not induce leukemia in vivo. We speculate that HOXA9 alone cannot maintain MYC expression at a level sufficient to confer leukemogenic ability in vivo. MLL fusion- and HOXA9/MEIS1-transduced cells, which are capable of inducing leukemia in vivo, expressed Myc 20–100% more than HOXA9-ICs (Figure 1C and Figure 4—figure supplement 1A). This additional MYC activity may be required to achieve sufficient leukemogenic ability, therefore MYC/HOXA9-transduced cells were able to induce leukemia in vivo (Figure 4E). It should be noted that MYC/HOXA9-ICs and -LCs expressed endogenous Myc at a lower level than other HOXA9-expressing cells (i.e. MLL-AF10- and HOXA9/MEIS1-ICs/LCs) (Figure 4—figure supplement 2), indicating that there is a threshold of MYC expression a cell can endure. We speculate that HOXA9 uplifts the threshold of MYC expression by conferring anti-apoptotic properties but cannot hyper-activate MYC expression by itself. Thus, additional MYC activation mediated by MLL fusions and MEIS1 confers proliferative advantages sufficient to induce leukemia in vivo.
The therapeutic efficacy of BCL2 inhibitor has been reported in AML including MLL fusion-mediated leukemia (Pan et al., 2014). Because there is a correlation between the expression levels of HOX proteins and the sensitivity to BCL2 inhibitor in AML patient samples, the aberrant expression of HOXA9 is the potential mechanism for BCL2-dependence of MLL leukemia (Brumatti et al., 2013; Kontro et al., 2017). Oncogenic MYC expression often leads to apoptosis, which needs to be alleviated by additional genetic events for leukemic cell survival (McMahon, 2014). Indeed, co-expression of HOXA9, BCL2, or SOX4 promoted MYC-mediated leukemogenesis, while MYC alone was insufficient to induce leukemia in vivo. Although the downstream mechanisms could not be addressed, SOX4 also exhibited anti-apoptotic effects on MYC-expressing cells. BCL2 was highly expressed in MLL fusion-mediated leukemia (Figure 5—figure supplement 2A,B), and many human MLL fusion leukemia cell lines were sensitive to a BCL-2 inhibitor (i.e. ABT-199)(Pan et al., 2014). High expression of SOX4 is correlated to poor survival in AML patients (Lu et al., 2017). These results support the key roles of BCL2 and SOX4 in the development of MLL leukemia. Ectopic expression of MYC and BCL2 induced lymphoid leukemia, whereas the MYC/SOX4 combination exclusively induced myeloid leukemia (Figure 6—figure supplement 1B), suggesting that BCL2 confers relatively stronger survival advantages in the lymphoid lineage than in the myeloid lineage in vivo, while SOX4 does in the myeloid lineage. Importantly, there was a partial effect of single-gene knockout of Bcl2 and Sox4 on leukemia initiation and maintenance. This indicates that MLL fusion proteins exploit multiple anti-apoptotic pathways and that blocking a single anti-apoptotic pathway may be insufficient to completely abrogate leukemic potential. Thus, simultaneously blocking multiple anti-apoptotic pathways may be required for efficient molecularly targeted therapy of HOXA9-expressing leukemia.
Our results also provide an insight into the mode of function of HOXA9. HOX genes are known to be expressed in a position-specific manner, conferring a positional identity to a cell (Wang et al., 2009). Our results suggest that HOXA9 unlikely functions as a major upstream factor determining tissue-specific gene expression by turning a silenced chromatin into transcriptionally active chromatin. Rather, HOX proteins likely play a supportive role in maintaining gene expression, which was activated by other transcriptional regulators. Our observation that HOXA9 cannot reactivate gene expression once silenced fits this hypothesis. Accordingly, HOXA9 maintains a subset of genes related to hematopoietic identity when expressed in hematopoietic precursors (Figure 2B). Recently, it has been reported that HOXA9 may recruit enhancer apparatuses, as it colocalizes with active enhancer mark (i.e. acetylated histone H3 lysine 27) (Sun et al., 2018). We speculate that HOXA9 may support the maintenance of an active enhancer, but unlikely establishes it on silenced chromatin. Further functional analysis of HOXA9 is required to understand how HOXA9 regulates gene expression.
In summary, our results describe the oncogenic roles for HOXA9 as a transcriptional maintenance factor for multiple anti-apoptotic genes, which are necessary to promote MYC-mediated leukemogenesis. In the case of MLL fusion-mediated leukemia, MLL fusion proteins directly activate both MYC and HOXA9, while HOXA9 maintains expression of MYC, BCL2, and SOX4, achieving high MYC activity and anti-apoptotic properties simultaneously (Figure 7E). Thus, MLL fusion-mediated LCs possess highly proliferative potentials and survival advantages using HOXA9 as a critical mediator.
Materials and methods
Vector constructs
Request a detailed protocolFor protein expression vectors, cDNAs obtained from Kazusa Genome Technologies Inc (Nagase et al., 2008) were modified by PCR-mediated mutagenesis and cloned into the pMSCV vector (for virus production) or pCMV5 vector (for transient expression) by restriction enzyme digestion and DNA ligation. The MSCV-neo MLL-ENL and MLL-AF10 vectors have been previously described (Okuda et al., 2017). sgRNA-expression vectors were constructed using the pLKO5.sgRNA.EFS.GFP vector (RRID:Addgene_57822) (Heckl et al., 2014). shRNA-expression vectors were purchased from Dharmacon. The target sequences are listed in Supplementary file 1.
Cell lines
Request a detailed protocolHEK293T cells were a gift from Michael Cleary and were authenticated by the JCRB Cell Bank in 2019 (Key Resources Table). HEK293TN cells were purchased from System Biosciences (RRID:CVCL_UL49). The cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and penicillin-streptomycin (PS). The Platinum-E (PLAT-E) ecotropic virus packaging cell line—a gift from Toshio Kitamura (RRID:CVCL_B488)(Morita et al., 2000)—was cultured in DMEM supplemented with 10% FBS, puromycin, blasticidin, and PS. The human leukemia cell lines including HB1119, (RRID:CVCL_8227), K562 (RRID:CVCL_0004), RS4-11 (RRID:CVCL_0093), THP1 (RRID:CVCL_0006), and EOL1 (RRID:CVCL_0258) were gifts from Michael Cleary (Tkachuk et al., 1992; Yokoyama et al., 2004), and were cultured in RPMI 1640 medium supplemented with 10% FBS and PS. The MV4-11 cell line (ATCC, RRID:CVCL_0064) was cultured in Iscove’s modified Dulbecco’s medium (IMDM) supplemented with 10% FBS and PS. The ML-2 cell line (DSMZ, RRID:CVCL_1418) was cultured in RPMI 1640 medium supplemented with 10% FBS and PS. MLL-ENL LCs (MLL-ENLbm0713) were described previously (Okuda et al., 2017). All the cell lines except HEK293TN, PLAT-E, and HB1119 were authenticated using short tandem repeat (STR)-PCR method by the JCRB Cell Bank. Cells were cultured in the incubator at 37°C and 5% CO2, and routinely tested for mycoplasma using the MycoAlert Mycoplasma detection kit (Lonza).
Animal models
Request a detailed protocolFor the in vivo leukemogenesis assay, 8-week-old female C57BL/6JJcl (C57BL/6J) or C. B-17/Icr-scid/scidJcl (SCID) mice were purchased from CLEA Japan (Tokyo, Japan). Bcl2-knockout mice were a gift from Yoshihide Tsujimoto and provided via RIKEN BRC (Kamada et al., 1995). Sox4-knockout mice were a gift from Hans Clevers and provided via RIKEN BRC (Schilham et al., 1996).
Western blotting
View detailed protocolWestern blotting was performed as previously described (Yokoyama et al., 2002). Antibodies used in this study are listed in Key Resources Table.
Virus production
Request a detailed protocolEcotropic retrovirus constructed in pMSCV vectors was produced using PLAT-E packaging cells (Morita et al., 2000). Lentiviruses were produced in HEK293TN cells using the pMDLg/pRRE (RRID:Addgene_12251), pRSV-rev (RRID:Addgene_12253), and pMD2.G (RRID:Addgene_12259) vectors, all of which were gifts from Didier Trono (Dull et al., 1998). The virus-containing medium was harvested 24–48 hr following transfection and used for viral transduction.
Myeloid progenitor transformation assay
Request a detailed protocolThe myeloid progenitor transformation assay was conducted as previously described (Lavau et al., 1997; Okuda and Yokoyama, 2017a). Bone marrow cells were harvested from the femurs and tibiae of 5-week-old female C57BL/6J mice. c-Kit+ cells were enriched using magnetic beads conjugated with an anti-c-Kit antibody (Miltenyi Biotec, RRID:AB_2753213), transduced with a recombinant retrovirus by spinoculation, and then plated (4 × 104 cells/ sample) in a methylcellulose medium (IMDM, 20% FBS, 1.6% methylcellulose, and 100 µM β-mercaptoethanol) containing murine stem cell factor (mSCF), interleukin 3 (mIL-3), and granulocyte-macrophage colony-stimulating factor (mGM-CSF; 10 ng/mL each). During the first culture passage, G418 (1 mg/mL) or puromycin (1 μg/mL) was added to the culture medium to select for transduced cells. Hoxa9 expression was quantified by qRT-PCR after the first passage. Cells were then re-plated once every 4–6 days with fresh medium; the number of plated cells for the second, third, and fourth passages was 4 × 104, 2 × 104, and 1 × 104 cells/well, respectively. CFUs were quantified per 104 plated cells at each passage.
In vivo leukemogenesis assay
Request a detailed protocolIn vivo leukemogenesis assays were conducted as previously described (Lavau et al., 1997; Okuda and Yokoyama, 2017b). c-Kit+ cells (2 × 105) prepared from the femurs and tibiae of 5-week-old female C57BL/6J mouse were transduced with retrovirus by spinoculation and intravenously transplanted into sublethally irradiated (5–6 Gy) C57BL/6J mice. For secondary leukemia, LCs (2 × 105) cultured ex vivo for more than three passages were transplanted. As for knockouts of Bcl2 and Sox4, mice heterozygous for Bcl2 or Sox4 were crossed, and c-Kit+ cells were isolated from fetal livers at E14–15 (for Bcl2) or E13 (for Sox4). The next day, cells were transduced with the retroviruses for MLL-AF10 or HOXA9/MEIS1 and transplanted intravenously into sublethally irradiated (2.5 Gy) SCID mice [2 × 105 (for Bcl2) or 1 × 105 cells/mouse (for Sox4)].
qRT-PCR
Request a detailed protocolTotal RNA was isolated using the RNeasy Mini Kit (Qiagen) and reverse-transcribed using the Superscript III First Strand cDNA Synthesis System (Thermo Fisher Scientific) with oligo (dT) primers. Gene expression was analyzed by qPCR using TaqMan probes (Thermo Fisher Scientific). Relative expression levels were normalized to those of GAPDH/Gapdh, ACTB/Actb, or TBP/Tbp and determined using a standard curve and the relative quantification method, according to manufacturer’s instructions (Thermo Fisher Scientific). Commercially available/custom made PCR probes used are listed in Supplementary file 1.
ChIP-qPCR and ChIP-seq
Request a detailed protocolChIP was performed as previously described (Okuda et al., 2017), using the fanChIP method (Miyamoto and Yokoyama, 2021). DNA was precipitated with glycogen, dissolved in TE buffer, and analyzed by qPCR (ChIP-qPCR) or deep sequencing (ChIP-seq). The qPCR probe/primer sequences are listed in Supplementary file 1. Deep sequencing was performed using the TruSeq ChIP Sample Prep Kit (Illumina) and HiSeq2500 (Illumina) at the core facility of Hiroshima University and described in our previous publication (Okuda et al., 2017).
RNA-seq
Request a detailed protocolTotal RNA was prepared using the RNeasy Kit (Qiagen) and analyzed using a Bioanalyzer (Agilent Technologies). Deep sequencing was performed using a SureSelect Strand Specific RNA Library Prep Kit (Agilent Technologies) and GAIIx (Illumina) with 36 bp single-end reads or HiSeq2500 (Illumina) with 51 bp single-end reads at the core facility of Hiroshima University. Sequenced reads were mapped to the human genome assembly hg19 or the mouse genome assembly mm9 using CASAVA 1.8.2 (Illumina, RRID:SCR_001802), and read counts were normalized as reads per kilo base of exon per million mapped (RPKM). To define HOXA9 high and MYC high signature genes, data were trimmed by removing lowly expressed genes whose RPKM values were less than 2, and the relative expression between HOXA9-ICs and MYC-ICs was estimated. The top 50 HOXA9 high and MYC high signature genes were visualized as a heatmap using the Complex Heatmap package (RRID:SCR_017270)(Gu et al., 2016). To define HOXA9-MYC common target genes, lowly expressed genes were removed, and genes with more than two-fold expression compared to c-Kit+ cells were determined for each of HOXA9-ICs and MYC-ICs. The overlapping genes were defined as ‘HOXA9-MYC common target genes’. GSEA analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed using the DAVID (RRID:SCR_001881)(Jiao et al., 2012), GSEA (Subramanian et al., 2005), and Metascape (Zhou et al., 2019).
sgRNA competition assay
Request a detailed protocolCas9 was introduced via lentiviral transduction using the pKLV2-EF1a-Cas9Bsd-W vector (RRID:Addgene_68343) (Tzelepis et al., 2016). Cas9-expressing stable lines were established with blasticidin (10–30 µg/mL) selection. The targeting sgRNA was co-expressed with GFP via lentiviral transduction using pLKO5.sgRNA.EFS.GFP vector (RRID:Addgene_57822)(Heckl et al., 2014). Percentages of GFP+ cells were initially determined by FACS analysis at 2 or 3 days after sgRNA transduction and then measured once every 3–5 days.
FACS analysis and sorting
Request a detailed protocolTo detect apoptosis, two to five million cells were suspended in 200 µL of reaction buffer (140 mM NaCl, 10 mM HEPES, 2.5 mM CaCl2, and 0.1% BSA) and incubated with APC-Annexin V (RRID:AB_2868885) for 15 min and Propidium iodide for 5 min at room temperature. The cells were then centrifuged, resuspended in fresh reaction buffer, and analyzed with FACS Melody (BD Bioscience). FACS sorting of mouse bone marrow cells was performed with fluorophore-conjugated antibodies listed in Key Resources Table as previously described (Yokoyama et al., 2013).
Accession numbers
Request a detailed protocolDeep sequencing data used in this study have been deposited in the DNA Data Bank of Japan (DDBJ) Sequence Read Archive under the accession numbers listed in Supplementary file 1.
Statistics
Statistical analyses were performed using GraphPad Prism seven software (RRID:SCR_002798). Data are presented as the mean with standard deviation (SD). Comparisons between two groups were analyzed by unpaired two-tailed Student’s t-test, while multiple comparisons were performed by ordinary one-way analysis of variance (ANOVA) followed by Dunnett’s test or two-way ANOVA. Mice transplantation experiments were analyzed by the log-rank test and Bonferroni correction was applied for multiple comparisons. p Values < 0.05 were considered statistically significant. n.s.: p>0.05, *: p ≤ 0.05, **: p ≤ 0.01, ***: p ≤ 0.001, and ****: p ≤ 0.0001.
Study approval
Request a detailed protocolAll animal experimental protocols were approved by the National Cancer Center (Tokyo Japan) Institutional Animal Care and Use Committee.
Appendix 1
Data availability
ChIP-seq data have been deposited to the DDBJ archive and have been published (accession number: DRA004871). RNA-seq data have been deposited to the DDBJ archive and have been published (accession number: DRA010090, DRA012078, DRA004874, DRA012079).
-
DDBJ GEAID E-GEAD-360. Expression profiles of murine myeloid progenitors immortalized by various oncogenes.
-
DDBJ GEAID E-GEAD-435. Expression profiles of murine myeloid progenitors immortalized by various oncogenes.
-
DDBJ GEAID E-GEAD-436. Expression profiles of human cell lines.
-
DDBJ DRAID DRA010090. Sequence reads of murine myeloid progenitors immortalized by various oncogenes.
-
DDBJ DRAID DRA012078. Sequence reads of murine myeloid progenitors immortalized by various oncogenes.
-
DDBJ DRAID DRA004874. Sequence reads of various factors/modificaions in HB1119 cells.
-
DDBJ DRAID DRA012079. Sequence reads of human cell lines.
-
DDBJ GEAID E-GEAD-321. Expression profiles of HB1119 and 293T cell lines.
References
-
Pbx modulation of hox homeodomain amino-terminal arms establishes different DNA-binding specificities across the hox locusMolecular and Cellular Biology 16:1734–1745.https://doi.org/10.1128/MCB.16.4.1734
-
Thirty years of BCL-2: translating cell death discoveries into novel Cancer therapiesNature Reviews Cancer 16:99–109.https://doi.org/10.1038/nrc.2015.17
-
A third-generation Lentivirus vector with a conditional packaging systemJournal of Virology 72:8463–8471.https://doi.org/10.1128/JVI.72.11.8463-8471.1998
-
Molecular genetics of human leukemias: new insights into therapySeminars in Hematology 39:6–11.https://doi.org/10.1053/shem.2002.36921
-
bcl-2 deficiency in mice leads to pleiotropic abnormalities: accelerated lymphoid cell death in Thymus and spleen, polycystic kidney, hair Hypopigmentation, and distorted small intestineCancer Research 55:354–359.
-
MYC and the control of apoptosisCold Spring Harbor Perspectives in Medicine 4:a014407.https://doi.org/10.1101/cshperspect.a014407
-
Cooperative gene activation by AF4 and DOT1L drives MLL-rearranged leukemiaJournal of Clinical Investigation 127:1918–1931.https://doi.org/10.1172/JCI91406
-
Myeloid progenitor transformation assayBio-Protocol 7:475–480.https://doi.org/10.21769/BioProtoc.2626
-
Global analysis of proliferation and cell cycle gene expression in the regulation of hematopoietic stem and progenitor cell fatesJournal of Experimental Medicine 202:1599–1611.https://doi.org/10.1084/jem.20050967
-
HOXA9 forms triple complexes with PBX2 and MEIS1 in myeloid cellsMolecular and Cellular Biology 19:3051–3061.https://doi.org/10.1128/MCB.19.4.3051
-
11q23/MLL acute leukemia : update of clinical aspectsJournal of Clinical and Experimental Hematopathology 50:91–98.https://doi.org/10.3960/jslrt.50.91
-
Integrative genomics viewer (IGV): high-performance genomics data visualization and explorationBriefings in Bioinformatics 14:178–192.https://doi.org/10.1093/bib/bbs017
-
Regeneration, repair and remembering identity: the three rs of hox gene expressionTrends in Cell Biology 19:268–275.https://doi.org/10.1016/j.tcb.2009.03.007
Article and author information
Author details
Funding
Japan Society for the Promotion of Science (16H05337)
- Akihiko Yokoyama
Japan Society for the Promotion of Science (19H03694)
- Akihiko Yokoyama
Japan Society for the Promotion of Science (19K16791)
- Ryo Miyamoto
The Yamagata Prefectural Government (Research grant)
- Akihiko Yokoyama
The City of Tsuruoka (Research grant)
- Akihiko Yokoyama
Dainippon Sumitomo Pharma Co., Ltd. (Research grant)
- Akihiko Yokoyama
The Friends of Leukemia Research Fund (Research grant)
- Akihiko Yokoyama
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Yuzo Sato, Boban Stanojevic, Makiko Okuda, Megumi Nakamura, Etsuko Kanai, Aya Nakayama, Hagumu Sato, Ikuko Yokoyama, Kanae Ito, and Ayako Yokoyama for technical assistance. We thank Drs. Yoshihide Tsujimoto and Hans Clevers for providing us the knockout mouse lines of Bcl2 and Sox4, respectively. These mouse lines were provided by the RIKEN BRC through the National BioResource Project of the MEXT, Japan. We also thank all members of the Shonai Regional Industry Promotion Center for their administrative support. This work was supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI grants (16H05337 and 19H03694 to AY; 19K16791 to RM) and in part by research funds from the Yamagata prefectural government, the City of Tsuruoka, Dainippon Sumitomo Pharma Co. Ltd., and the Friends of Leukemia Research Fund.
Ethics
Animal experimentation: All animal experimental protocols were approved by the National Cancer Center (Tokyo Japan) Institutional Animal Care and Use Committee.
Copyright
© 2021, Miyamoto 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,774
- views
-
- 245
- downloads
-
- 17
- 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
- Computational and Systems Biology
This study investigates the variability among patients with non-small cell lung cancer (NSCLC) in their responses to immune checkpoint inhibitors (ICIs). Recognizing that patients with advanced-stage NSCLC rarely qualify for surgical interventions, it becomes crucial to identify biomarkers that influence responses to ICI therapy. We conducted an analysis of single-cell transcriptomes from 33 lung cancer biopsy samples, with a particular focus on 14 core samples taken before the initiation of palliative ICI treatment. Our objective was to link tumor and immune cell profiles with patient responses to ICI. We discovered that ICI non-responders exhibited a higher presence of CD4+ regulatory T cells, resident memory T cells, and TH17 cells. This contrasts with the diverse activated CD8+ T cells found in responders. Furthermore, tumor cells in non-responders frequently showed heightened transcriptional activity in the NF-kB and STAT3 pathways, suggesting a potential inherent resistance to ICI therapy. Through the integration of immune cell profiles and tumor molecular signatures, we achieved an discriminative power (area under the curve [AUC]) exceeding 95% in identifying patient responses to ICI treatment. These results underscore the crucial importance of the interplay between tumor and immune microenvironment, including within metastatic sites, in affecting the effectiveness of ICIs in NSCLC.
-
- Cancer Biology
- Genetics and Genomics
A new approach helps examine the proportion of cancerous and healthy stem cells in patients with chronic myeloid leukemia and how this influences treatment outcomes.