Search

Targeted-sequencing analysis of PDX and PairTree predicted mutational population clusters

dPDX and dRI-PDX leukemia-initiating cell frequencies

Disease recurrence causes significant mortality in B-progenitor acute lymphoblastic leukemia (B-ALL). Genomic analysis of matched diagnosis and relapse samples shows relapse often arising from minor diagnosis subclones. However, why therapy eradicates some subclones while others survive and progress to relapse remains obscure. Elucidation of mechanisms underlying these differing fates requires functional analysis of isolated subclones. Here, large-scale limiting dilution xenografting of diagnosis and relapse samples, combined with targeted sequencing, identified and isolated minor diagnosis subclones that initiate an evolutionary trajectory toward relapse [termed diagnosis Relapse Initiating clones (dRI)]. Compared with other diagnosis subclones, dRIs were drug-tolerant with distinct engraftment and metabolic properties. Transcriptionally, dRIs displayed enrichment for chromatin remodeling, mitochondrial metabolism, proteostasis programs, and an increase in stemness pathways. The isolation and characterization of dRI subclones reveals new avenues for eradicating dRI cells by targeting their distinct metabolic and transcriptional pathways before further evolution renders them fully therapy-resistant.Significance:Isolation and characterization of subclones from diagnosis samples of patients with B-ALL who relapsed showed that relapse-fated subclones had increased drug tolerance and distinct metabolic and survival transcriptional programs compared with other diagnosis subclones. This study provides strategies to identify and target clinically relevant subclones before further evolution toward relapse.See related video: https://vimeo.com/442838617See related article by E. Waanders et al.

Patient characteristics and patient genomic analysis (WES and SNP 6.0)

Leukemia-initiating cell frequencies of paired diagnosis and relapse patient samples

RNA-sequencing data from Patient 9 PDX

RNA-sequencing, pathway enrichment (GSEA) reports and GSVA results (including gene list HSC vs B) of PDX and paired patient samples

PDX targeted-sequencing tissue concordance

Purpose: In many children with cancer and characteristics suggestive of a genetic predisposition syndrome, the genetic cause is still unknown. We studied the yield of pathogenic mutations by applying whole-exome sequencing on a selected cohort of children with cancer.Experimental Design: To identify mutations in known and novel cancer-predisposing genes, we performed trio-based whole-exome sequencing on germline DNA of 40 selected children and their parents. These children were diagnosed with cancer and had at least one of the following features: (1) intellectual disability and/or congenital anomalies, (2) multiple malignancies, (3) family history of cancer, or (4) an adult type of cancer. We first analyzed the sequence data for germline mutations in 146 known cancer-predisposing genes. If no causative mutation was found, the analysis was extended to the whole exome.Results: Four patients carried causative mutations in a known cancer-predisposing gene: TP53 and DICER1 (n = 3). In another 4 patients, exome sequencing revealed mutations causing syndromes that might have contributed to the malignancy (EP300-based Rubinstein–Taybi syndrome, ARID1A-based Coffin–Siris syndrome, ACTB-based Baraitser–Winter syndrome, and EZH2-based Weaver syndrome). In addition, we identified two genes, KDM3B and TYK2, which are possibly involved in genetic cancer predisposition.Conclusions: In our selected cohort of patients, pathogenic germline mutations causative or likely causative of the cancer phenotype were found in 8 patients, and two possible novel cancer-predisposing genes were identified. Therewith, our study shows the added value of sequencing beyond a cancer gene panel in selected patients, to recognize childhood cancer predisposition. Clin Cancer Res; 24(7); 1594–603. ©2018 AACR.

Supplementary Table S1 - Cancer gene panel Supplementary Table S2 - Patients with multiple malignancies Supplementary Table S3 - Tumor types of index patients Supplementary Table S4 - Candidate genes Supplementary Table S5 - Overview of clinical features and mutations found in category 1 (ID and/or congenital anomalies) Supplementary Table S6 - Overview of clinical features and mutations found in category 2 (cancer twice) Supplementary Table S7 - Overview of clinical features and mutations found in category 3 (family history) Supplementary Table S8 - Overview of clinical features and mutations found in category 4 (adult type of cancer) Supplementary Table S9 - Overview of clinical features and mutations found in category 5 (multiple reasons for inclusion) Supplementary Figure S1 - Detection of mosaic mutation in ARID1A Supplementary Figure S2 - Family tree of case #04 Supplementary Figure S3 - Family tree of case #05 Supplementary Figure S4 - Family tree of case #07 Supplementary Figure S5 - Family tree of case #10 Supplementary Figure S6 - Family tree of case #12 Supplementary Figure S7 - Family tree of case #17 Supplementary Figure S8 - Family tree of case #21 Supplementary Figure S9 - Family tree of case #22 Supplementary Figure S10 - Family tree of case #24 Supplementary Figure S11 - Family tree of case #27 Supplementary Figure S12 - Family tree of case #28 Supplementary Figure S13 - Family tree of case #29 Supplementary Figure S14 - Family tree of case #34 Supplementary Figure S15 - Family tree of case #36 Supplementary Figure S16 - Family tree of case #39

This file contains all mutations identified in all patients.

B-cell acute lymphoblastic leukemia (B-ALL) readily disseminates to the leptomeninges of the central nervous system (CNS). CNS involvement is more frequent in certain poor prognosis subgroups including patients with MLL-AF4 translocations, and late CNS relapse is often lethal. The biology and clonal history of CNS leukemia are poorly defined and consequently therapies exploiting drivers of metastasis are lacking. To characterize leptomeningeal leukemia we performed targeted DNA sequencing, SNP copy number analysis, RNA sequencing, and functional analysis on cells isolated from the bone marrow (BM) and CNS of xenografts generated from a cohort of paired diagnosis and relapse samples from 14 B-ALL patients. The majority of patient samples disseminated to the CNS following intrafemoral injection into irradiated NSG mice, with greater CNS involvement in xenografts derived from relapse patient samples. Secondary transplantation of both BM- and CNS-purified cells demonstrated their capacity to re-engraft BM, CNS, and spleen. Targeted-sequencing results were analyzed using a Bayesian clustering method to determine the clonal composition of matched BM and CNS, demonstrating discordance in subclonal prevalence in nearly half the xenografts tested. Xenografts derived from two patient samples demonstrated recurrent enrichment of a particular subclone in the CNS versus BM. Similarly, copy number analysis identified frequent discordance between BM and CNS tissues within individual mice. All xenografts from one patient exhibited chromosome 6p and 17p hemi-deletions that were exclusive to CNS cells. While these data suggest that individual B-ALLs harbor subclones with CNS tropism, there were no recurrently enriched single nucleotide mutations or copy number alterations across all patients. RNA-sequencing of 45 BM and CNS pairs from primary xenografts demonstrated that CNS-isolated cells were consistently distinct from their matched BM. GSEA analysis of xenografts generated from patients with MLL-AF4 translocations (MLL) (n=2 patients, 26 mice), identified CNS cell enrichment of gene sets related to mRNA translation and nascent peptide elongation compared to BM. MLL-CNS cells exhibited altered rates of protein synthesis compared to BM cells from the same mouse. The clinically-approved translation inhibitor omacetaxine mepesuccinate effectively diminished protein translation rates of CNS isolated cells and reduced CNS engraftment by four fold in xenografts derived from two MLL-AF4 patients. These data demonstrate that the CNS microenvironment selects for the outgrowth of B-ALL cells with genetically and/or biologically distinct properties. Moreover, we demonstrate that in MLL-AF4 patients, altered protein synthesis occurs in CNS dissemination and that targeting this process may clinically benefit patients with CNS disease. Citation Format: Robert J. Vanner, Stephanie M. Dobson, Ildiko Grandal, Olga Gan, Jessica McLeod, James Kennedy, Veroniqu Voisin, Abdellatif Daghrach, Erwin M. Schoof, Cynthia Guidos, Jayne Danska, Esme Waanders, Mark Minden, Charles G. Mullighan, John E. Dick. Genetic profiling of central nervous system dissemination of B-acute lymphoblastic leukemia reveals clonal selection and therapeutic vulnerability [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5173.

Current cancer therapies are directed at molecular markers or dominant pathways present in the bulk of neoplastic cells. However, recent studies have identified many genetically distinct subclones co-existing within a single neoplasm. In over 50% of patients with relapsed acute lymphoblastic leukemia (ALL), the genetic clones present at relapse are not the dominant clone present at diagnosis, but have evolved from a minor or ancestral clone (Mullighan et al., Science, 2008). Previous work has shown that this subclonal diversity in B-ALL exists at the level of the leukemia-initiating cells (L-IC) capable of generating patient derived xenografts (Notta et al., Nature, 2011). In order to investigate the functional consequences of genetic clonal evolution during disease progression, we performed in-depth genomic and functional analysis of 14 paired diagnosis/relapse samples from adult and pediatric B-ALL patients of varying cytogenetics. Patient samples were subjected to whole exome sequencing (WES), SNP analysis and RNA sequencing. Diagnosis-specific, relapse-specific, and shared variants at both clonal and subclonal frequencies were identified. Limiting dilution analysis by transplantation of CD19+ leukemic blasts into 870 immune-deficient mice (xenografts) identified no significant trend in enrichment in L-IC frequency between paired patient samples with a median frequency of 1 in 2691. Despite similar frequencies of L-IC, functional differences within identically sourced patient xenografts were observed, including increased leukemic dissemination of relapse cells to distal sites such as the central nervous system (CNS), differences in engraftment levels and differences in immunophenotypes. Targeted-sequencing and copy number analysis of the xenografts, in comparison to the patient sample from which they were derived, has uncovered clonal variation and the unequivocal identification of minor subclones ancestral to the relapse in xenografts transplanted with the diagnostic sample from 8 patients. Some of these subclones are rare and were not captured through standard WES analysis of the patient samples, highlighting the value of xenografting to functionally identify and viably isolate subclones for further study. Interrogation of the therapeutic responses of the ‘relapse-like’ diagnosis subclones in secondary xenografts displayed differential resistance to standard chemotherapeutic agents (vincristine and L-asparaginase) pre-existing in the patient diagnosis samples prior to treatment. Furthermore, investigation of different sites of leukemic infiltration in the xenografts provided evidence of distinct clonal selection in the CNS, a known site of disease relapse, in comparison to the bone marrow. Using this data we can begin to draw the evolutionary paths to relapse. We have shown evidence that minor subclones at diagnosis, ancestral to the relapsing clone, possess functional advantages over other diagnostic clones. Overall, this work provides a substantial advance in connecting genetic diversity to functional consequences, thereby furthering our understanding of the heterogeneity identified in B-ALL and its contributions to therapy failure and disease recurrence. Citation Format: Stephanie M. Dobson, Robert Vanner, Esmé Waanders, Olga I. Gan, Jessica McLeod, Ildiko Grandal, Debbie Payne-Turner, Michael Edmonson, Zhaohui Gu, Xioatu Ma, Yiping Fan, Pankaj Gupta, Sagi Abelson, Michael Rusch, Ying Shao, Scott Olsen, Geoffrey Neale, John Easton, Cynthia J. Guidos, Jayne S. Danska, Jinghui Zhang, Mark D. Minden, Charles G. Mullighan, John E. Dick. Evolving functional heterogeneity in B-acute lymphoblastic leukemia. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr LB-341.

Current cancer therapies are directed at molecular markers or dominant pathways present in the bulk of neoplastic cells. However, recent studies have identified many genetically distinct subclones co-existing within a single neoplasm. In over 50% of patients with relapsed acute lymphoblastic leukemia (ALL), the genetic clones present at relapse are not the dominant clone present at diagnosis, but have evolved from a minor or ancestral clone (Mullighan et al., Science, 2008). Previous work has shown that this subclonal diversity in B-ALL exists at the level of the leukemia-initiating cells (L-IC) capable of generating patient derived xenografts (Notta et al., Nature, 2011). In order to investigate the functional consequences of genetic clonal evolution during disease progression, we performed in-depth genomic and functional analysis of 14 paired diagnosis/relapse samples from adult and pediatric B-ALL patients of varying cytogenetics. Patient samples were subjected to whole exome sequencing (WES), SNP analysis and RNA sequencing. Diagnosis-specific, relapse-specific, and shared variants at both clonal and subclonal frequencies were identified. Limiting dilution analysis by transplantation of CD19+ leukemic blasts into 870 immune-deficient mice (xenografts) identified no significant trend in enrichment in L-IC frequency between paired patient samples with a median frequency of 1 in 2691. Despite similar frequencies of L-IC, functional differences within identically sourced patient xenografts were observed, including increased leukemic dissemination of relapse cells to distal sites such as the central nervous system (CNS), differences in engraftment levels and differences in immunophenotypes. Targeted-sequencing and copy number analysis of the xenografts, in comparison to the patient sample from which they were derived, has uncovered clonal variation and the unequivocal identification of minor subclones ancestral to the relapse in xenografts transplanted with the diagnostic sample from 8 patients. Some of these subclones are rare and were not captured through standard WES analysis of the patient samples, highlighting the value of xenografting to functionally identify and viably isolate subclones for further study. Interrogation of the therapeutic responses of the ‘relapse-like’ diagnosis subclones in secondary xenografts displayed differential resistance to standard chemotherapeutic agents (vincristine and L-asparaginase) pre-existing in the patient diagnosis samples prior to treatment. Furthermore, investigation of different sites of leukemic infiltration in the xenografts provided evidence of distinct clonal selection in the CNS, a known site of disease relapse, in comparison to the bone marrow. Using this data we can begin to draw the evolutionary paths to relapse. We have shown evidence that minor subclones at diagnosis, ancestral to the relapsing clone, possess functional advantages over other diagnostic clones. Overall, this work provides a substantial advance in connecting genetic diversity to functional consequences, thereby furthering our understanding of the heterogeneity identified in B-ALL and its contributions to therapy failure and disease recurrence. Citation Format: Stephanie M. Dobson, Robert Vanner, Esme Waanders, Olga I. Gan, Jessica McLeod, Ildiko Grandal, Debbie Payne-Turner, Michael Edmonson, Zhaohui Gu, Xioatu Ma, Yiping Fan, Pankaj Gupta, Sagi Abelson, Michael Rusch, Ying Shao, Scott Olsen, Geoffrey Neale, John Easton, Cynthia J. Guidos, Jayne S. Danska, Jinghui Zhang, Mark D. Minden, Charles G. Mullighan, John E. Dick. Evolving functional heterogeneity in B-acute lymphoblastic leukemia. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr LB-341.