Research ArticleCancer

Overcoming mutational complexity in acute myeloid leukemia by inhibition of critical pathways

See allHide authors and affiliations

Science Translational Medicine  25 Oct 2017:
Vol. 9, Issue 413, eaao1214
DOI: 10.1126/scitranslmed.aao1214

The right treatments for the right mutations

A variety of mutations have been observed in cancer cells from patients with acute myeloid leukemia, but it can be difficult to know which of these mutations contribute to tumorigenesis and should therefore be targeted. To address this issue, Saito et al. isolated subpopulations of leukemic cells with specific mutations and monitored their leukemogenic capacity in immunosuppressed mice. By combining this approach with genomic analysis, the authors were able to identify mutations that drive the evolution of leukemia and figure out effective approaches to target them.

Abstract

Numerous variant alleles are associated with human acute myeloid leukemia (AML). However, the same variants are also found in individuals with no hematological disease, making their functional relevance obscure. Through NOD.Cg-PrkdcscidIl2rgtmlWjl/Sz (NSG) xenotransplantation, we functionally identified preleukemic and leukemic stem cell populations present in FMS-like tyrosine kinase 3 internal tandem duplication–positive (FLT3-ITD)+ AML patient samples. By single-cell DNA sequencing, we identified clonal structures and linked mutations with in vivo fates, distinguishing mutations permissive of nonmalignant multilineage hematopoiesis from leukemogenic mutations. Although multiple somatic mutations coexisted at the single-cell level, inhibition of the mutation strongly associated with preleukemic to leukemic stem cell transition eliminated AML in vivo. Moreover, concurrent inhibition of BCL-2 (B cell lymphoma 2) uncovered a critical dependence of resistant AML cells on antiapoptotic pathways. Co-inhibition of pathways critical for oncogenesis and survival may be an effective strategy that overcomes genetic diversity in human malignancies. This approach incorporating single-cell genomics with the NSG patient-derived xenograft model may serve as a broadly applicable resource for precision target identification and drug discovery.

INTRODUCTION

Acute myeloid leukemia (AML) is a biologically and clinically heterogeneous entity. Recent studies using deep DNA and RNA sequencing have demonstrated intrapatient heterogeneity (14). Mutations in genes, such as NPM1, TET2, WT1, IDH1/2, and DNMT3A, are commonly found in AML (58), and next-generation DNA sequencing (NGS) suggests that certain mutations occur earlier than others based on variant allele frequencies (911). Preleukemic stem cells carrying these early somatic mutations may contribute to leukemogenesis and disease relapse (9, 10). On the other hand, large-scale population-based sequencing studies have revealed that hematopoietic cells in 5 to 18.4% of elderly subjects with nonmalignant conditions, such as diabetes mellitus and cardiovascular diseases, harbored somatic mutations in genes including ASXL1, DNMT3A, and TET2 (1214). In these subjects, the mutations were associated with a 0.5 to 1% annual rate of progression to hematological malignancies. This raised some fundamental questions regarding leukemogenesis and treatment strategies. First, which among AML-associated recurrent mutations contribute to leukemogenesis? Second, which of the many mutations and pathways must be targeted for greatest clinical efficacy? To address these questions, we examined mutational profiles of phenotypically and functionally defined human AML cell populations to link mutations with in vivo fates. We then used a functional single-cell genomic approach to identify critical targets, allowing in vivo elimination of human AML cells with multiple coexisting mutations.

RESULTS

In vivo fates of human AML cells are linked with distinct mutational profiles through NSG xenotransplantation

We obtained bone marrow (BM) or peripheral blood (PB) samples from 27 patients with FMS-like tyrosine kinase 3 internal tandem duplication–positive (FLT3-ITD+) AML (table S1). Most of the patients had poor prognostic factors, such as complex chromosomal abnormalities in addition to FLT3-ITD mutation, and/or had a known aggressive disease (for example, primary resistance or relapse after multiple stem cell transplantations). Because AML-associated hematopoiesis consists of both normal and malignant cells, we profiled patterns of recurrent mutations in patient-derived cell populations purified according to cell surface phenotype that defines hematopoietic stem cells (HSCs), multipotent progenitor cells, multilymphoid progenitors, and mature lymphoid and myeloid cells (15). To link these mutations with in vivo fates, we transplanted the cell populations in newborn NSG mouse recipients (Fig. 1). If human lymphoid and myeloid subsets were engrafted in NSG recipients (multilineage human hematopoietic repopulation), then the transplanted subpopulation contained normal HSCs and/or preleukemic stem cells. If NSG recipients developed leukemia with uncontrolled proliferation of myeloid blasts and without lymphoid differentiation, then the transplanted subpopulation contained leukemia-initiating cells (LICs). As expected, frequencies of hematopoietic subpopulations and their population-level mutational profiles varied among patients, and frequencies of mutated alleles varied among subpopulations within individual patients (representative data from six patients in Fig. 2 and fig. S1; sequence information in table S2). Upon transplantation, subpopulations with similar surface phenotypes isolated from different patients showed distinct behaviors in vivo. For instance, the patient 21–derived CD34+CD38CD90CD45RA cell population initiated AML in NSG mice and therefore contained LICs (Fig. 2A). In contrast, in patients 20, 23, and 24, the CD34+CD38CD90CD45RA cell population reconstituted multilineage human hematopoiesis in NSG mice and therefore contained multilineage-engrafting HSCs or preleukemic stem cells, whereas LICs were present in the CD34+CD38CD90CD45RA+ cell population (Fig. 2B and fig. S1, A and B). In patient 13, the CD34+CD38 population reconstituted multilineage human hematopoiesis, and the CD34CD33+ population contained LICs, whereas patient 1–derived CD34+CD38 cells initiated AML in vivo (Fig. 2C and fig. S1C). These observations are consistent with recent reports showing variable cell surface phenotype of LICs (16). We next examined mutational profiles in these subpopulations with defined in vivo fates. In patient 20, the same set of mutations (FLT3-ITD, DNMT3A, and WT1) was present in CD34+CD38CD90CD45RA and CD34+CD38CD90CD45RA+ subpopulations, but these subpopulations showed distinct in vivo fates (Fig. 2B). This functional difference may be due to uneven distribution of mutations identified in bulk AML cell populations in single-cell clones, resulting in disparate combinations of mutations and divergent in vivo fates. Therefore, we performed single-cell DNA sequencing of functionally defined preleukemic stem cell– and LIC-containing subpopulations along with human multilineage hematopoietic cells and leukemia cells they generate in vivo to define clonal structures and identify mutation(s) associated with leukemia-initiating versus multilineage-engrafting function.

Fig. 1. In vivo fates of patient-derived AML cells defined by mutational profile.

(A) Somatic mutation profiles were identified in patient cell subpopulations defined by surface phenotype based on developmental hierarchy of human hematopoiesis. (B) The in vivo fate of each subpopulation was determined through transplantation into newborn NSG mice. If repopulation by multilineage hematopoiesis occurred, then the transplanted subpopulation contained hematopoietic or preleukemic stem cells; if AML engraftment occurred, then the transplanted subpopulation contained LICs.

Fig. 2. Mutational profiles and in vivo engraftment of patient-derived subpopulations.

Three additional patients are shown in fig. S1. In each patient sample, human CD45+CD3+CD19 T cells and human CD45+CD3CD19+ B cells were identified. Within CD3CD19 non-T, non-B cells, subpopulations were identified on the basis of CD34, CD38, CD90, and CD45RA surface expression. These populations underwent polymerase chain reaction (PCR) for FLT3-ITD mutation and DNA sequencing (DNA-seq) for the other genes indicated. Variant allele frequencies are shown as heat maps. In patients 21 (A) and 20 (B), the CD34+CD38CD90CD45RA and CD34+CD38CD90CD45RA+ subpopulations were identified. The in vivo fates of CD34+CD38CD90CD45RA subpopulations differed between patients 20 and 21, showing engraftment with multilineage human hematopoiesis in patient 20 (indicated by green rectangles) but initiation of AML in patient 21 (indicated by red rectangles). In patient 13 (C), the CD34+CD38 subpopulation showed multilineage repopulation, whereas the CD34CD33+ subpopulation with additional FLT3-ITD and NPM1 mutations initiated AML. AML-engrafted recipients showed no B cell engraftment (indicated by gray dashed outlines on flow cytometry plots). Detailed information on variants found in each patient is shown in table S2.

Functional genomic approach combining patient-derived xenograft model and single-cell DNA sequencing distinguishes leukemogenic from permissive mutations

We examined DNMT3A, TET2, NPM1, and WT1 mutations and FLT3-ITD among single cells isolated from multilineage-engrafting cell– and LIC-containing patient-derived populations and their in vivo progeny (Fig. 3A). Patient-derived multilineage-engrafting preleukemic stem cells showed mutational heterogeneity at the single-cell level, carrying combinations of multiple mutations (Fig. 3, B and C). DNMT3A mutation was identified both in multilineage-engrafting CD34+CD38CD90CD45RA preleukemic cells (patients 20 and 24) and AML-initiating CD34+CD38CD90CD45RA cells (patient 21) at the single-cell level (Fig. 3B). In vivo–generated single B cells from patient 20– and patient 24–derived preleukemic stem cells harbored DNMT3A mutation (and WT1 mutation in the case of patient 20), whereas there were no FLT3-mutated B cell clones. There were no FLT3-ITD mutations in patient 24–derived preleukemic stem cells and in vivo–generated single B cells. Although there were FLT3-ITD–mutated single cells in the patient 20–derived preleukemic stem cell population, those FLT3-ITD–mutated subclones did not contribute to normal lymphopoiesis. These findings indicate that DNMT3A and WT1 mutations are permissive and can coexist in a single cell without hindering human multilineage differentiation. In contrast, FLT3-ITD was identified in substantial proportions of patient 21–, patient 20–, and patient 24–derived LICs and engrafted AML cells at the single-cell level. Likewise, CD34+CD38 cells derived from patients 1 and 13 exhibited distinct in vivo fates: At the single-cell level, the latter carried wild-type (WT) FLT3, and the former harbored FLT3-ITD (Fig. 3C). Note that, although FLT3-ITD+ single cells were a minority among the LIC population in patient 13, every engrafted AML cell harbored FLT3-ITD mutation. In addition, FLT3-ITD+ single cells were enriched among LIC-containing CD34CD33+ population at the time of relapse in patient 13. These findings suggest that acquisition of the FLT3-ITD mutation acts as a critical trigger for leukemia initiation, working in cooperation with accumulated mutations in DNMT3A, TET2, NPM1, and/or WT1.

Fig. 3. Mutations contributing to distinct in vivo cell fates identified by single-cell functional genomics.

(A) Patient-derived subpopulations with defined in vivo fates and their in vivo progeny underwent single-cell mutation profiling. Through this strategy, mutations present in patient-derived preleukemic stem cells (pre-LSCs) and LIC clones were tracked and linked to in vivo fates. (B and C) Using samples from five patients, patient-derived multilineage-engrafting and LIC-containing population and engrafted B cells and AML cells were subjected to single-cell DNA sequencing for variants detected in each indicated gene in each patient. FLT3-ITD sequences with highly variable repeated sequence patterns were detected by single-cell PCR. In (B) and (C), each column of rectangles represents an individual cell. The presence or absence of mutations in each gene is shown by colors of rectangles, as indicated.

Kinase inhibition effectively targets human AML with mutational diversity

Therapeutic efficacy of targeting such a mutation among multiple coexisting mutations is an important question for clinical translation. We addressed this by in vivo inhibition of the FLT3 pathway in an NSG patient-derived xenograft (PDX) model. To serve as a realistic platform for in vivo therapeutic testing, a PDX model must reflect the mutational diversity of human AML cells. At the single-cell level, we detected AML cells with various patterns of mutations in engrafted recipients. In addition, engrafted AML cells retained the mutations present in patient-derived LIC-containing cells from 12 patients examined, indicating that engrafted AML cells reflect mutational diversity of patient-derived AML cells (Fig. 4). Therefore, we went on to examine the effect of FLT3 pathway inhibition in human AML cells with multiple coexisting mutations by using RK-20449, a pyrrolo-pyrimidine derivative inhibitor of Src family kinase HCK and FLT3 that we had previously identified (17). By treating 56 NSG mice that were engrafted with FLT3-ITD+ AML from 19 patients, we found significant responses to the single agent RK-20449 in vivo in all 19 cases (data and associated P values are shown in table S3). For five patients, RK-20449 completely eliminated AML cells in the BM, spleen, and PB of all recipient mice treated, despite the presence of mutations not directly targeted by RK-20449 (patient 1: DNMT3A, NPM1, and TET2; patient 2: IDH1; patient 16: FLT3 D835H point mutation and NRAS) (Fig. 5A and table S2). For 11 additional patients, RK-20449 treatment alone resulted in complete responses in the spleen in most of the recipient mice tested (Fig. 5B). Although statistically significant (P < 0.05) treatment effects were observed in groupwise comparisons, residual RK-20449–resistant AML cells were present in the BM of at least one treated mouse for 14 of 19 patients (Fig. 5, B and C).

Fig. 4. Profiles of AML-associated mutations in nine genes in patient- and recipient-derived AML cells from 12 AML cases.

Variant allele frequencies of indicated genes are represented as heat maps. Patient, LIC-containing population from the patient; 1′, human AML cells from primary recipients; 2′, human AML cells from secondary recipients. All patient-derived and recipient-derived leukemia populations were positive for FLT3-ITD by PCR. Non-ITD FLT3 mutations were identified by sequencing. Information on variants is shown in table S2.

Fig. 5. Induction of apoptosis via enhanced BCL-2 dependence in FLT3-ITD+ AML cells with diverse coexisting somatic mutations by in vivo kinase inhibition.

Overall, treatment with RK-20449 resulted in significant reduction of AML cells in the BM, spleen, and PB of recipients (P < 1 × 10−19 for each; data are tabulated in table S3). To document patient-to-patient variability, RK-20449 responses were classified as follows: complete response (A), if all recipients treated showed residual BM human CD45+ chimerism of <5%; good response (B), if the case did not meet the criterion for complete responder but all recipients showed <50% residual BM human CD45+; partial response (C), if at least one of the recipients showed >50% residual BM human CD45+. For each response group, PB time course of human AML chimerism (leftmost panels) for RK-20449–treated recipients and final BM (middle panels) and spleen (rightmost panels) human AML chimerism of RK-20449–treated and untreated recipients are shown. Pretreatment PB human AML cell chimerism is shown at week 0. The numbers of recipients for each patient/each treatment group and pre-/posttreatment AML chimerism are shown in table S3. In all response groups, BM, spleen, and PB chimerism was significantly reduced with RK-20449 treatment. For the BM and spleen, final chimerism for recipients in each treatment group was compared. For the PB, pre- and posttreatment chimerism for RK-20449–treated recipients was compared. P < 5.8 × 10−5 by two-tailed t test for all comparisons. In each scatter graph, dotted lines are drawn at 0, 5, and 50%. (D) Dynamics of apoptotic response to BIM peptide in the presence of RK-20449 was measured for six AML cases. Bars represent the BIM IC50 of cytochrome c loss in RK-20449–treated human CD45+ cells as percentages of IC50 in cells treated with dimethyl sulfoxide (DMSO) alone. Enhancement of apoptotic response to BIM by RK-20449 showed substantial patient-to-patient variability. (E) Dynamics of apoptotic response to BH3-only peptides BAD, HRK, and NOXA in the presence of RK-20449 or DMSO alone was measured for seven AML cases. Bars represent the percentage reduction of IC50 in the presence of RK-20449 compared with DMSO alone. RK-20449 enhanced apoptotic response to BAD and HRK peptides with substantial patient-to-patient variability, whereas apoptotic response to NOXA peptide was not substantially enhanced by RK-20449. (F) Apoptotic response to a BCL-2 inhibitor ABT-199 was enhanced by RK-20449 in seven AML cases. Bars show increased cytochrome c loss in human AML cells treated with ABT-199 and RK-20449 compared with those treated with DMSO alone.

BCL-2 inhibition enhances mitochondrial priming to apoptosis induced by kinase inhibition

Resistance to RK-20449 may be mediated through pre-existing or newly acquired somatic mutations (18). However, enrichment or new acquisition of AML-associated somatic mutations was not identified in five AML cases after in vivo RK-20449 treatment (fig. S2). In addition, the WT FLT3 gene was not identified in resistant cells, indicating that emergence of FLT3 WT cells is not a substantial mechanism for resistance. This is consistent with mutational profiles of patient samples serially obtained at primary presentation and at relapse (fig. S3). In six cases examined, neither emergence of FLT3 WT cells nor substantial increase in frequencies of somatic mutations was detected, with the exception of FLT3 D835H mutation–positive cells in patient 16 emerging at the time of relapse. In vivo RK-20449 treatment resulted in transcriptional up-regulation of S100A8 (associated with drug resistance in leukemia), HSPA5 (promotes cell survival under endoplasmic reticulum stress and suppresses ferroptosis), and IFI6 (negatively regulates apoptosis) (fig. S4) (1923). Therefore, we functionally assessed dependence on antiapoptotic mechanisms in RK-20449–resistant human AML cells by dynamic BCL-2 (B cell lymphoma 2) homology domain 3 (BH3) profiling (24). Some human malignancies are dependent on specific antiapoptotic proteins for survival and are therefore sensitive to the small-molecule antagonists of those proteins (2529). Dynamic BH3 profiling determines how “primed” cells are to apoptotic cell death and how changing conditions (such as exposure to drugs) affect baseline priming by quantifying mitochondrial cytochrome c release in response to BH3-only peptides that activate proapoptotic effectors BAX and BAK. Despite patient-to-patient variability, RK-20449 treatment lowered the half maximal inhibitory concentration (IC50) of BIM peptide for mitochondrial cytochrome c release, indicating enhanced proapoptotic signaling in FLT3-ITD+ human AML cells (Fig. 5D). This was consistent with our finding that RK-20449 alone completely eliminated AML cells with DNMT3A, TET2, IDH1, and/or WT1 mutations. In addition, exposure to RK-20449 facilitated mitochondrial cytochrome c release in response to BAD and HRK peptides, indicating enhanced mitochondrial sensitivity to BCL-2 and BCL-XL inhibition (Fig. 5E). In addition, dynamic BH3 profiling showed that a selective BCL-2 inhibitor, ABT-199, enhanced RK-20449–induced apoptosis in FLT3-ITD+ human AML cells (Fig. 5F). This suggested that co-inhibition of antiapoptotic signal, specifically BCL-2, might result in a more robust eradication of human AML with diverse coexisting mutations.

Combined inhibition of kinase and antiapoptosis pathways eliminates AML with genetic complexity and clonal diversity

To examine whether combining RK-20449 and ABT-199 eliminates resistant AML cells, we chose cases that were not completely eradicated in vivo by RK-20449 alone. In most of the cases, the BCL-2 inhibitor ABT-199 alone resulted in transient and limited responses. In contrast, in all 12 cases, combination treatment significantly reduced human AML chimerism in the PB, BM, and spleen (Fig. 6, A and B, and tables S4 to S6; P < 0.05 in all comparisons). Combination treatment completely eliminated human AML cells in vivo without targeting coexisting mutations in 9 of 12 cases (Fig. 6B; cases 3, 5 to 10, 13, and 14). We compared residual leukemia-initiating capacity by transplanting a predetermined fraction of viable human CD45+ AML cells from treated mice into secondary NSG recipients (Fig. 6C). We found that BM treated with RK-20449 alone or ABT-199 alone contained enough viable LICs to initiate AML in every secondary recipient transplanted, whereas AML cells remaining after combination treatment engrafted in only 1 of 23 mice, indicating that combination treatment more effectively reduced the frequency of LICs in vivo (Fig. 6C). On the other hand, combination treatment did not have significant effects on human leucocytes, T cells, B cells, and myeloid cells in human cord blood (CB) HSC–engrafted NSG recipients (table S7).

Fig. 6. Eradication of FLT3-ITD+ AML cells in vivo through combined inhibition of kinase and antiapoptotic pathways.

(A and B) Mice engrafted with AML derived from 12 patients associated with good or partial responses to RK-20449 alone received four different treatments (no treatment, ABT-199 alone, RK-20449 alone, and combination). The numbers of recipients for each patient/each treatment group and pre-/posttreatment AML chimerism are shown in table S4. (A) PB human CD45+ chimerism is shown over time. Recipients were phlebotomized weekly, and pretreatment PB human CD45+ AML chimerism is shown at time 0. Mean PB human CD45+ chimerism for each patient/each treatment group and statistics comparing the treatment groups are shown in table S5. (B) Final BM and spleen human AML chimerism is shown for mice engrafted with AML derived from 12 patients. Nine cases showed complete responses, and three cases showed good responses to combination treatment. Each circle represents an AML-engrafted recipient. Mean BM/spleen human CD45+ chimerism for each patient/each treatment group and statistics comparing the treatment groups are shown in table S6. In each scatter graph, dotted lines are drawn at 0, 5, and 50%. (C) Residual human AML initiation capacity of human CD45+ cells after in vivo treatment was assessed by serial transplantation for four treatment groups. To compare the amount of residual LICs in recipients after in vivo treatment, each secondary recipient was transplanted with human CD45+ cells sorted from 2.5% (by cell number) of viable BM cells remaining in treated recipients. The mean and SEM for human CD45+ AML cell chimerism in the BM of secondary recipients are shown. Each circle represents a secondary recipient.

DISCUSSION

Multilineage-engrafting preleukemic HSCs carrying somatic mutations are thought to be poised for leukemogenesis and may act as a reservoir for leukemic progression and relapse (9, 10). To prevent such events, early-occurring somatic mutations or founder mutations with high variant allele frequencies have been considered as critical therapeutic targets (30). However, some of these somatic mutations were found in individuals with advanced age with no apparent hematological disease, associated with a relatively low rate of progression to hematological malignancies (fig. S5, top left) (31). Here, we demonstrated that the contribution of these somatic mutations to normal hematopoietic differentiation or to leukemogenesis could inform therapeutic target selection in AML. By integrating population- and single-cell–level genomics and in vivo functional assessment in PDXs, we found that relatively late acquisition of FLT3-ITD on the background of permissive earlier-occurring mutations in DNMT3A, TET2, NPM1, and WT1, alone or in combination, triggered in vivo leukemogenesis (fig. S5, top right). Moreover, acquisition of FLT3-ITD triggered leukemogenesis along a broad spectrum of hematopoietic differentiation. Through single-cell sequencing, we found substantial mutational heterogeneity in patient-derived preleukemic cells and NSG-engrafted human lymphoid cells, whereas patient-derived LICs were less mutationally heterogeneous at the single-cell level, with substantial enrichment of FLT3-ITD+ single cells. This is consistent with the hypothesis that late acquisition of FLT3-ITD mutation results in selective clonal expansion. Single-cell RNA sequencing may help fully characterize clonal structures of preleukemic and leukemic stem cells.

In a recent clinical study, midostaurin, a small-molecule inhibitor of FLT3, improved both 5-year event-free and overall survivals in FLT3-ITD+ AML patients when combined with standard chemotherapy, making it the first U.S. Food and Drug Administration–approved targeted agent in AML (3234). Our findings are consistent with this clinical finding: Targeting FLT3-ITD alone resulted in effective elimination of AML in vivo despite coexisting earlier-occurring mutations (fig. S5, bottom left and right). Furthermore, dynamic BH3 profiling identified an additional target, BCL-2, in kinase inhibitor–resistant AML cells. Co-inhibition of BCL-2 resulted in complete eradication of human AML cells resistant to kinase inhibition. Note that this approach was effective in cases with poor prognoses and highly aggressive clinical courses, multiple coexisting mutations, and/or complex chromosomal abnormalities. With the advent of high-throughput sequencing technologies and the genomic information gained from each patient, precision medicine is becoming more and more feasible. By functionally connecting genomic information with in vivo fate and behavior of patient-derived cells at the level of single cells through a PDX model, we could identify therapeutic targets with improved precision to promote more effective drug discovery in genetically complex human malignancies.

MATERIALS AND METHODS

Study design

The overall objective of this study was to explore strategies for eliminating FLT3-ITD+ human AML by (i) tracing profiles of known somatic mutations associated with AML and other malignancies in single-cell clones of patient-derived hematopoietic cells to human hematopoietic cells engrafting in NSG xenotransplantation recipients to identify mutations associated with leukemia initiation, (ii) examining whether the antiapoptosis pathway can be a therapeutic target independent of leukemogenic somatic mutations, and (iii) testing whether targeting these pathways by small-molecule inhibitors would result in efficient elimination of human AML in vivo. To do so, we isolated various subpopulations of patient-derived hematopoietic cells by fluorescence-activated cell sorting (FACS) and performed NGS of bulk and single-cell genomic DNA for known malignancy-associated somatic mutations in parallel with NSG xenotransplantation. By performing NGS of bulk and single-cell DNA in human hematopoietic cells that were engrafted in NSG recipients, we traced and identified somatic mutations affiliated with in vivo leukemia-initiating AML cell clones and preleukemic cell clones that reconstituted nonmalignant human hematopoiesis in vivo. To determine the degree to which human AML cells resistant to kinase inhibition were dependent on the antiapoptotic pathway, we performed dynamic BH3 profiling. To test the efficacy of kinase inhibition and inhibition of the antiapoptotic pathway, we treated human FLT3-ITD+ AML-engrafted NSG recipients with a dual kinase inhibitor, RK-20449, and a BCL-2 inhibitor, ABT-199, and assessed human AML cell chimerism by flow cytometry in the PB weekly during treatment and in the BM, spleen, and PB at the time the animal was sacrificed.

We did not predetermine the numbers of mice in each treatment group. To ensure that each comparison was sufficiently powered, we performed power calculation for each comparison (either two-tailed t test or paired two-tailed t test) using StatMate (GraphPad). We found that the comparisons deemed statistically significant (P < 0.05) were powered at 85 to 99% in detecting the differences observed using the statistical test used. PB sampling was carried out before treatment and every week after start of treatment. Recipients were checked daily and sacrificed when they showed signs of progressive disease such as ruffled fur and weakness. No randomization and blinding were performed, and there were no exclusions.

For in vivo treatment experiments, recipients with similar extent of PB human AML engraftment chimerism were chosen as a set for indicated treatments. Pretreatment human CD45 chimerism and statistical data are shown in tables S3 to S5. There were no significant differences found in pretreatment human AML engraftment between treatment groups.

Because it is logistically difficult and undesirable from the standpoint of animal health and comfort to obtain replicate samples from highly immunosuppressed NSG recipients, PB was sampled only once every week for all recipients in all experiments. Human AML chimerism obtained at the time of sacrifice was also evaluated once in the BM, spleen, and PB of the recipients. Overall experimental replication for in vivo treatment studies was ensured by the numbers of patient samples tested and the numbers of recipients treated for each group. For each patient sample, treatment experiments were performed as a set of untreated control, treatment A, treatment B, and/or treatment C, with A, B, and C being ABT-199 alone, RK-20449 alone, or combination, respectively. Therefore, there were at least three experimental replicates for each patient sample.

Ethical statements

Written informed consent was obtained from all patients. The study was performed with authorization from the Institutional Review Board for Human Research at RIKEN and Toranomon Hospital, in accordance with the ethical standards of responsible committees on human experimentation at each institution. CB samples were purchased from Lonza. All experiments using NSG mice (35, 36) were performed with authorization from and according to guidelines established by the Institutional Animal Committees at RIKEN and the Jackson Laboratory.

Mice

Mice were bred and maintained under defined flora at the animal facility of RIKEN and at the Jackson Laboratory. Both female and male NSG mice received transplants at 2 days of age. Treatment studies were conducted when sufficient engraftment was observed at about 6 weeks of age.

Flow cytometry

The following monoclonal antibodies (mAbs) were used for flow cytometry: mouse anti-human CD19 (catalog nos. 562653, 555412, and 341093), CD3 (catalog nos. 563800, 562426, and 555341), CD33 (catalog nos. 562854 and 555450), CD34 (catalog no. 348791), CD38 (catalog no. 340439 and 555459), CD4 (catalog no. 555348), and CD45 (catalog nos. 563204, 641399, and 563204); and rat anti-mouse Ter119 (catalog no. 557915) and CD45 (catalog nos. 563891 and 563410) (BD Biosciences). Analyses were performed with FACSAria III and FACSCanto II (BD Biosciences). To obtain cells for xenogeneic transplantation, BV786 (Brilliant Violet 786)–conjugated anti-CD3, BV605-conjugated anti-CD19, BV421-conjugated anti-CD33, PE-Cy7 (phycoerythrin–cyanin 7)–conjugated anti-CD34, allophycocyanin-conjugated anti-CD38, fluorescein isothiocyanate–conjugated anti-CD90, and PE-conjugated anti-CD45RA mAbs were used. For single-cell sorting, a 100-μm nozzle was used.

Transplantation

NSG newborn mice received 1.5-gray total body irradiation, followed by intravenous injection of purified human cells. For primary transplantation shown in Figs. 2 to 4, cells of the indicated phenotype were sorted from BM or PB cells obtained from the patients, and 5 × 102 to 2 × 105 cells were transplanted per recipient, depending on the frequency of the cell population. For in vivo treatment experiments, the known LIC population from each patient was transplanted to create human AML-engrafted recipients. Donor cells were purified according to their cell surface phenotype using mAbs against human CD34, CD38, CD90, CD45RA, CD3, CD19, and CD33. The extent of engraftment of human cells in the NSG recipients was assessed by retro-orbital phlebotomy and flow cytometry.

Genome analysis

DNA was extracted from human cells purified from patient samples or recipient organs using the DNeasy Blood and Tissue kit (Qiagen). PCR detection of FLT3-ITD was performed using the TaKaRa PCR FLT3/ITD Mutation Detection set (Takara Bio).

The bulk DNA sequences were determined by NGS. After shearing with a Covaris S220 (Covaris), the fragmented genomic DNA (10 ng) was converted to an NGS sequencing library with a KAPA Hyper Prep kit (Kapa Biosystems) according to the protocol provided by the supplier. Targeted sequencing of AML-related genes was carried out by a hybridization capture method with xGen AML Cancer Panel v1.0 (Integrated DNA Technologies) according to the protocol provided by the supplier. The hybridization-captured DNA library was subjected to NGS in a paired-end read mode (200 cycles) with an Illumina HiSeq 1500 (Illumina). The obtained DNA sequences were mapped onto human genome sequence (hg19) using Burrows-Wheeler Aligner (BWA) (v0.7.12; http://bio-bwa.sourceforge.net/) and then realigned with a Realigner Target Creator in the Genome Analysis Toolkit (v1.6-13; https://software.broadinstitute.org/gatk/index.php). After treatment with Fix Mate Information in Picard (v1.119; http://broadinstitute.github.io/picard/) and Quality Score Covariate and Table Recalibration in the Genome Analysis Toolkit (v1.6-13), variants were detected with VarScan (v2.3.6; read depth, >10; http://varscan.sourceforge.net/somatic-calling.html).

Single-cell variation analysis was carried out for single cells sorted on a BD FACSAria into 96-well plates. Single-cell whole-genome amplification by Multiple Annealing and Looping-Based Amplification Cycles (MALBAC) (37), associated with a substantially low allelic dropout rate at ~1% and a false-positive rate of ~4 × 10−5, was used (3739). After single-cell genome amplification with a MALBAC Single Cell WGA kit (Yikon Genomics), target regions of genes of interest were PCR-amplified with primers including well indexes by PCR. The first-round gene-specific PCR was conducted using gene-specific primers by 25-cycle PCR with Gflex DNA Polymerase (Takara Bio) in a single-plex mode. The second-round PCR was performed to link Illumina anchor sequences at both sides of the first-round PCR products. PCR conditions were as in the first-round PCR, except that PCR primers were replaced with those for attaching Illumina anchoring sequences. Because of low PCR efficiency, the first-round PCRs for WT1 and CEBPA were carried out as follows: For WT1, PCR was performed over 40 cycles using Taq DNA polymerase (Qiagen) under a three-step thermal cycling of 94°C for 30 s, 55°C for 30 s, and 72°C for 30 s; for CEPBA, PCR was performed over 36 cycles using Taq DNA polymerase and Q-solution (Qiagen) under a two-step thermal cycling of 95°C for 1 min and 68°C for 6 min. Primer sequences are described in table S8. The PCR products were sequenced in a paired-end read mode (300 cycles) on an Illumina MiSeq. The obtained DNA sequences were mapped onto the human genome sequence (hg19) with BWA (v0.7.12), and paired-end reads were merged with SAMtools (v1.0; read depth, >100; http://samtools.sourceforge.net/). Variant detection and frequency calculation were carried out with mpileup in SAMtools. Here, variation genotype was assigned on an assumption that the sequencing error rate is lower than 2% and an apparent variation with allele frequency lower than 2% was regarded as a WT. Three independent amplification and sequencing rounds were performed before calling a locus WT. No germline tissue was available for evaluation of somatic status of mutations. CEBPA, DNMT3A, FLT3, IDH1, TET2, and WT1 variants were included or excluded according to gene-specific characteristics, as described by Lindsley et al. (40).

For RNA sequencing (RNA-seq) analysis, total RNA was extracted from FACS-purified cells treated with TRIzol (Thermo Fisher Scientific). RNA-seq libraries were prepared using the SureSelect Strand-Specific RNA Library Preparation kit (Agilent Technologies) according to the manufacturer’s protocol and were sequenced by a HiSeq 1500 (Illumina). The sequence reads were mapped to the human genome (hg19) using TopHat2 software (v2.0.8). Cufflinks (v2.1.1) was run with the same reference annotation with TopHat2 to generate FPKM (fragments per kilobase of transcript per million mapped reads) values. Statistical evaluation of gene expression change was performed using the edgeR algorithm with read counts on exons determined using the R program.

In vivo treatment

In vivo treatment experiments were performed with AML-engrafted NSG recipients using RK-20449 (17) and ABT-199 (41, 42). The recipients were treated with RK-20449 (30 mg/kg) intraperitoneally twice a day, ABT-199 (70 mg/kg) orally once a day, or both RK-20449 and ABT-199. The mice were euthanized when they became moribund or after 4 to 6 weeks of treatment, and human AML chimerism in the BM, spleen, and PB was determined using flow cytometry. In secondary transplantation experiments, each mouse received 7-aminoactinomycin D(−) viable human CD45+ cells from 2.5% of total BM that remained in AML-engrafted recipients at the time of sacrifice to simulate relapse occurring from residual viable AML cells. All treated recipients and their pre- and posttreatment engraftment data are tabulated in tables S3 and S4. No sample size pre-estimation was performed. To ensure that each comparison was sufficiently powered, we performed power calculation for each comparison (either two-tailed t test or paired two-tailed t test) using StatMate (GraphPad). Each comparison deemed significant (P < 0.05) was powered at 85 to 99% for detecting the differences observed.

Dynamic BH3 profiling

Cells were harvested from the BM of recipients engrafted with AML derived from patients with FLT3-ITD+ AML, and BH3 profiling was performed using the plate-based assay previously described (43, 44). For dynamic BH3 profiling, 2 × 106 harvested cells were exposed to 500 nM RK-20449 in 2 ml of hematopoietic growth medium supplemented with stem cell factor (50 ng/ml), FLT3 ligand (50 ng/ml), and thrombopoietin (50 ng/ml) for 4 hours at 37°C in humidified atmosphere containing 5% CO2. After surface labeling with anti-human CD45 and dead cell exclusion with Zombie NIR (BioLegend), the cells were permeabilized and exposed to BH3 peptides (0.39 μM BIM, 80 μM BAD, 80 μM HRK, or 80 μM NOXA) or 1 μM ABT-199, and retained intracellular cytochrome c was measured by flow cytometry using anti–cytochrome c antibody.

Statistical analysis

For in vivo treatment experiments, numerical data are presented as means ± SEM. The differences were examined with two-tailed t test (GraphPad Prism, GraphPad). Statistical methods for genomic analyses are included in the “Genome analysis” section.

For in vivo treatment experiments, pre- and posttreatment data were obtained from the PB of each mouse. To detect differences between pre- and posttreatment data obtained from the same mouse, paired two-tailed t test (pairing pre- and posttreatment values of each mouse) was used. Because pretreatment data from the BM/spleen of the mice are not available, differences in the BM and spleen of the mice were evaluated by unpaired two-tailed t test. We designed the experiment such that there were three to six engrafted recipients per treatment group from a given patient. Therefore, in comparisons restricted to mice engrafted with AML from a particular patient, n was insufficient (three to six in each group) for meaningful tests of normality or variation. In all such comparisons, performing nonparametric tests did not yield contradicting results on the significance of the differences detected, and the comparisons were sufficiently powered, although n was small, because the sizes of the differences detected were sufficiently large and the values observed within each group were sufficiently tightly distributed (SD was sufficiently low). For comparisons of treatment groups including all mice (across patient samples), we chose to use the t test because parametric tests on non-Gaussian data are robust as long as sample sizes are sufficient. In these cases, using nonparametric tests also did not yield contradicting results.

Because n was relatively small (<20) for groups restricted to a particular patient sample, independent data points were plotted in all relevant figures, and all data are tabulated in tables S1 to S8. In comparisons of whole treatment groups (containing experiments using samples from multiple patients), normality testing using two different tests (D’Agostino-Pearson omnibus test and Shapiro-Wilk test) yielded variable results. However, because n is relatively large, detected differences are relatively large, and deviations among data points within each group are relatively small; parametric test (two-tailed t test) should be robust. Performing nonparametric tests on these data sets did not contradict the results of the parametric two-tailed t test.

In comparisons between groups restricted to individual patients, n was insufficient to obtain a meaningful estimate on the variance. For aggregate data of mice engrafted with AML from multiple patients, n was sufficiently large; however, in such comparisons, we did not expect variances to be similar because patients have highly heterogeneous disease characteristics and biology. As expected, the F test used to compare variances showed similar variances between some groups but not in others. In comparisons between groups restricted to individual patients, n was insufficient to obtain a meaningful estimate on the variance. We did not expect variances to be similar for aggregate data across multiple patient samples and comparisons made on mice engrafted with AML from multiple patients with highly heterogeneous disease characteristics and biology.

SUPPLEMENTARY MATERIALS

www.sciencetranslationalmedicine.org/cgi/content/full/9/413/eaao1214/DC1

Fig. S1. Identification of hematopoietic subpopulations in patient samples by surface phenotype and in vivo function in three additional patients.

Fig. S2. Mutational profiles of AML cells with and without in vivo RK-20449 treatment.

Fig. S3. Mutational profiles of AML patient samples obtained serially during clinical course.

Fig. S4. Altered transcription profiles of human AML cells with in vivo RK-20449 treatment.

Fig. S5. Overcoming genetically complex AML by targeting leukemia-initiating mutation and antiapoptotic BCL-2 pathway.

Table S1. Patient characteristics (provided as an Excel file).

Table S2. Information of identified variants (provided as an Excel file).

Table S3. Human AML chimerism in the PB, BM, and spleen of AML-engrafted recipients treated with RK-20449 alone (provided as an Excel file).

Table S4. Human AML chimerism in the PB, BM, and spleen of AML-engrafted recipients treated with ABT-199 alone, RK-20449 alone, or combination (provided as an Excel file).

Table S5. Statistics comparing pre- and posttreatment PB human AML chimerism in recipients treated with ABT-199 alone, RK-20449 alone, or combination (provided as an Excel file).

Table S6. Statistics comparing BM and spleen human AML chimerism in recipients treated with ABT-199 alone, RK-20449 alone, or combination (provided as an Excel file).

Table S7. Effect of in vivo exposure to combined RK-20449 and ABT-199 on human multilineage hematopoiesis (provided as an Excel file).

Table S8. PCR primers for targeted sequencing of MALBAC products (provided as an Excel file).

Reference (45)

REFERENCES AND NOTES

Acknowledgments: We would like to acknowledge and thank the patients and the clinical care team at the Toranomon Hospital, whose cooperation and effort made this study possible. We thank R. C. Lindsley and J. S. Garcia for the critical review of the manuscript. Funding: This study was supported by the Basic Science and Platform Technology Program for Innovative Biological Medicine (to F.I.) and the Project for Development of Innovative Research on Cancer Therapeutics (to F.I.) from the Japanese Ministry of Education, Culture, Sports, Science, and Technology, the Japan Agency for Medical Research and Development (to F.I.), and the Maine Cancer Foundation and the NIH (grants CA034196 and CA17983 to L.D.S.). Author contributions: F.I., Y.S., O.O., and A.G.L. conceptualized and designed the study. F.I., O.O., Y.M., and T.W. developed the methodology. I.O., L.H., K.S., A.K., and H.K. acquired the data. F.I., Y.S., O.O., A.G.L., Y.M., and T.W. analyzed and interpreted the data. F.I., Y.S., O.O., A.G.L., and L.D.S. wrote, reviewed, and/or revised the manuscript. T.F. provided administrative, technical, and material support. F.I. and Y.S. supervised the study. S. Takagi, N.U., and S. Taniguchi provided the patient samples and clinical information. Competing interests: Y.S. and F.I. are inventors on patent application (62/394871) submitted by RIKEN that covers the combinatory use of RK-20449 and BCL-2 inhibitors for AML. All other authors declare that they have no competing interests. Data and materials availability: Raw sequence data reported in this paper are available at the National Bioscience Database Center Human database (Japan) (accession no. hum0116). NSG mice are available from the Jackson Laboratory under a material transfer agreement.
View Abstract

Navigate This Article