Multiple TICs in H9M AML
Normal HSPC express c-kit, lack mature lineage markers, and produce c-kit expressing myeloid or lymphoid progenitors that give rise to mature cells (
Fig S1A). Four unique/expanded populations in H9M mice were assessed for TIA: (i) Lin
−kit
+, lacking mature markers and expressing c-kit, analogous to an HSPC, (ii) Lym
+kit
+, co-expressing lymphoid markers and c-kit, analogous to a lymphocyte precursor, (iii) Gr1
+kit
+, co-expressing c-kit and myeloid markers, representing an immature myeloid cell, and (iv) Gr1
+kit
lo, resembling a mature myeloid population (compare WT & H9M immunophenotypes, and
Fig S1E–J). The difference between the hematopoietic continuum of WT and H9M marrow was also visualized using SPADE, a density normalization, agglomerative clustering, and minimum-spanning tree algorithm that distills multidimensional single-cell data to interconnected “relatedness” clusters (with corresponding populations in conventional flow plots) of metacell populations displayed as a 2D tree plot () (
Qiu et al., 2011). The SPADE tree, drawn using composite data from WT-GFP, primary and secondary H9M-GFP cells, facilitated comparison of common and distinct metacells in each condition (size indicates relative frequency; compare ). As expected, in WT-GFP, the rare, c-kit expressing progenitors lacked lineage marker expression, lymphoid metacells largely cluster together and are distinct from the myeloid region (). In contrast, H9M-GFP cells have increased c-kit expression in the progenitor and myeloid regions, and a diminished lymphoid population ().
The Lin
−kit
+, Lym
+kit
+, Gr1
+kit
+, and Gr1
+kit
lo compartments were assessed for TIA. In retroviral transplant models,
in vitro myeloid colony-forming cells are TICs (
Somervaille and Cleary, 2006). Each population was sorted to high purity and colony-forming activity (CFA) was assessed. The Lin
−kit
+, and unexpectedly, Lym
+kit
+ compartments had robust myeloid CFA (>10%; ). In contrast, the Gr1
+kit
+ cells had a low initial CFA (<1%) which increased upon replating; Gr1
+kit
lo cells lacked any CFA. The most stringent test of TIA is transplanting disease with limiting cell numbers
in vivo (
Clarke et al., 2006). All recipients transplanted with 100 primary Lin
−kit
+ or Lym
+kit
+ cells succumbed rapidly to disease while more than 1000 Gr1
+kit
+ and 10,000 Gr1
+kit
lo cells were required to initiate disease in secondary recipients ( &
Fig S2A–B). Therefore, primary H9M AML is organized as a hierarchy in which multiple, phenotypically distinct TICs (Lin
−kit
+, Lym
+kit
+) are enriched for TIA and give rise to differentiated blasts (Gr1
+kit
lo) lacking TIA. The Gr1
+kit
+ compartment in primary AML has low, but seemingly selectable TIA.
Stage-dependent Differentiation Is Separate from TIA
Constitutive GFP expression allowed
in vivo lineage tracking of progeny. Each c-kit
+ compartment recapitulated the entire immunophenotypic spectrum of primary AML in secondary recipients including cells (for the Lym
+kit
+ and Gr1
+kit
+ compartments) of a less differentiated immunophenotype (Lin
−kit
+) as well as cells with a phenotype corresponding to mature cells of the opposite lineage (i.e., Gr1
+kit
+ cells produced cells expressing lymphoid markers and vice versa; ). To determine if this immunophenotypic recapitulation was clonal, single Lym
+kit
+ and Gr1
+kit
+ cells were expanded into colonies that were injected into secondary recipients (). Two of twenty recipients, transplanted with colonies derived from single Lym
+kit
+ and Gr1
+kit
+ cells respectively, developed AML and these leukemias again recapitulated the complete immunophenotypic spectrum (). Moreover, pre-leukemic H9M-transduced progenitors expressing myeloid markers (CD11b
+Gr1
+) gave rise to the immunophenotypic spectrum when transplanted to primary recipients (
Fig S3), confirming that enforced H9M expression permits cells with defined lineage commitment to assume immunophenotypes of “earlier” or alternate lineages, and that traditional phenotypes associated with developmental stage specificity are separable from TIA.
Dynamic Hierarchy in Secondary AML
To determine whether secondary tumors had the same hierarchical tumorigenic capacity as the primary tumors, 100 Lin−kit+, Lym+kit+, Gr1+kit+, or Gr1+kitlo cells from secondary AML were transplanted into tertiary recipients. As in primary AML, all animals transplanted with 100 Lin−kit+ or Lym+kit+ cells from secondary AML succumbed rapidly to disease (). Strikingly, although no animals transplanted with 100 Gr1+kit+ cells from primary AML developed leukemia, all animals transplanted with 100 Gr1+kit+ cells from secondary AML developed fatal AML indistinguishable (except longer latency, 38 v. 23 days, p=0.002, Log-Rank Test) from tumors arising from Lin−kit+ or Lym+kit+ subsets (). As before, Gr1+ kitlo cells failed to reconstitute tumors () demonstrating a continued lack of TIA in this fraction. Thus, while TIA does not remain fixed within immunophenotypic compartments as tumors progress, H9M AML remains hierarchically organized in secondary disease as TIC give rise to terminal blasts without TIA.
To determine whether H9M AML, like normal hematopoiesis, has a directional developmental progression, cohorts of mice were transplanted with each primary TIC and groups of mice were sacrificed sequentially to assess the immunophenotypes of leukemic progeny and capture individual variability of leukemic expansion (). GFP
+ cells were first detectable on day 7-post transplantation. SPADE analysis showed substantial heterogeneity in immunophenotypic reconstitution prior to day 14, after which nearly all mice (11/12) consolidated to a uniform immunophenotypic equilibrium that was nearly identical to primary AML ( and
Fig S4). Thus, while each TIC ultimately recapitulated the immunophenotypic spectrum of primary AML, unlike in normal development, the
in vivo developmental program that did not appear to follow a strict directional progression. Whether the regeneration occurs by a transdifferentiation-like mechanism, or dedifferentiation followed by forward differentiation remains to be determined, as it rapid and beyond the current level of resolution.
TICs Share Genetic, Signaling and Survival Pathways
We hypothesized these phenotypically distinct TIC, which shared biological and developmental potential, would utilize overlapping genetic, signaling and survival pathways, and targeting such shared pathways could affect TIC survival. Microarray analysis and unsupervised hierarchal clustering was performed on TICs from primary H9M AML, populations from MLL-AF9 AML with high or low TIA (Mac1
+kit
+ or Mac1
+kit
lo, respectively) (
Somervaille and Cleary, 2006), and analogous populations from WT-GFP marrow (Lin
−kit
+, Lym
+, and Gr1
+) (). This analysis showed H9M Lin
−kit
+ and Lym
+kit
+ compartments clustered with normal HSPC (GFP:Lin
−kit
+), H9M Gr1
+kit
+ cells clustered with the Mac1
+kit
+ TIC from MLL-AF9 AML, and terminal blasts clustered together (H9M:Gr1
+kit
lo and MLL-AF9:Mac1
+kit
lo). This suggests that TICs in primary H9M AML (Lin
−kit
+ and Lym
+kit
+), share a genetic signature resembling normal HSPC, as occurs in human AML (
Eppert et al., 2011), and is distinct from the program utilized in MLL-AF9 TICs.
To determine if primary H9M TICs shared enrichment for leukemogenesis-associated genes, the top 2500 differentially expressed probes between Lin
−kit
+ and Lym
+kit
+ cells relative to primary Gr1
+kit
+ cells (which lacked TIA) were assessed and common genes were identified (
Table S1). GO-term enrichment analysis of the resultant 151 genes (94 with annotations) showed the major biological functions were cellular metabolism (26%), signal transduction regulation (19%), and transcription/chromatin organization (18%) (). Specifically, this list included targets whose improper regulation has been previously implicated in AML such as Satb1, DNMT3A, and Msi2 (
Kharas et al., 2010;
Ley et al., 2010;
Steidl et al., 2007). Therefore, primary H9M TICs shared expression of genes implicated in leukemogenesis, suggesting a core cellular program can be induced by H9M expression—although the roles of these genes as drivers or passengers remains to be elucidated. Despite this shared program, these populations are unique as 12 of the top 20 genes differentially expressed between them have roles in B-cell development, Lym
+kit
+ cells express B220, and have rearranged the D-J gene segments on the immunoglobulin heavy chain locus (
Fig S5). Therefore, the Lym
+kit
+ progenitor represents a lymphoid-specified AML TIC, potentially analogous to other lymphoid-primed AML TICs (
Deshpande et al., 2006;
Goardon et al., 2011).
To determine if TICs shared signaling networks, mass cytometry was utilized providing a single-cell view of TIC signaling architecture (
Bendall et al., 2011). After measuring potentiated phosphorylation of 14 intracellular nodes in each TIC after cytokine stimulation, analysis of variance (ANOVA) was utilized to test for statistically significant differences in the fold-change values (FCV) between TICs for each stimulation-phospho node pair (). This analysis showed no statistically significant differences in the FCV for 26 of 28 conditions tested (
Fig S6A), demonstrating these distinct TICs have largely conserved signaling networks (unlike WT samples where FCV between compartments were significantly different for 3 of 5 conditions (
Fig S6B)).
To determine which nodes had conserved roles in TIC proliferation, each compartment was sorted into liquid media supplemented with S36GM—which selectively expands TIC but not terminal blasts in vitro—and proliferation was assessed after 7 days in the presence of pharmacological inhibitors targeting a suite of cellular pathways implicated in TIC-maintenance (). For inhibitors targeting the MAP kinase pathway (Mek inhibitor PD-0325901), DNA methyltransferase (5-aza-2'-deoxycytidine), receptor tyrosine kinases (imatinib), and PI3K/Akt (GDC-0941), termed class I inhibitors, the IC50 values for each TIC were overlapping or within 5-fold, suggesting these nodes play conserved roles in TIC survival and proliferation. In contrast, inhibitors targeting Src-family kinases (dasatinib) and mTOR (KU-0063794), class II inhibitors, had IC50 values greater than 10-fold different in the TICs, suggesting disparate requirements for these pathways across TICs. Collectively, these proliferation and signaling data show that overlapping survival and signaling nodes govern phenotypically distinct TICs.
Targeting Conserved Nodes in vivo Increases Survival
We hypothesized that pathways with conserved roles in TIC survival could serve as effective therapeutic targets in vivo while those with disparate roles would not. To determine if conserved signaling pathways were essential for AML maintenance in vivo, secondary recipients transplanted with total GFP+ H9M marrow were treated with either class I (PD-0325901, 5-aza-2'-deoxycytidine) or class II (Dasatinib) inhibitors starting two weeks after transplant to allow for immunophenotypic reconstitution. All mice treated with vehicle or the class II inhibitor dasatinib (with 10-fold differential selectivity across the TICs) succumbed rapidly to AML (median survival 25 days; ). Strikingly, mice treated with the class I inhibitors PD-0325901 and 5-aza-2'-deoxcytidine had statistically significant increases in overall survival (median survival 37 days, p=0.01 and 42 days, p<0.001, respectively, Log-Rank Test). These results show that in vivo pharmacologic targeting of pathways with conserved roles in TIC cell survival in vitro significantly increased survival, suggesting phenotypically distinct TICs rely on conserved survival nodes in vivo and targeting nodes only important for certain classes of TIC in the tumor allows for other TICs to “escape” and for disease to rapidly progress.