The correlation pattern of acquired copy number changes in 164 ETV6/RUNX1-positive childhood acute lymphoblastic leukemias

Optical Genome Mapping Identifies Novel Recurrent Structural Alterations in Childhood ETV6::RUNX1+ and High Hyperdiploid Acute Lymphoblastic Leukemia

The mutational landscape of B-cell precursor acute lymphoblastic leukemia (BCP-ALL), the most common pediatric cancer, is not fully described partially because commonly applied short-read next generation sequencing has a limited ability to identify structural variations. By combining comprehensive analysis of structural variants (SVs), single-nucleotide variants (SNVs), and small insertions-deletions, new subtype-defining and therapeutic targets may be detected. We analyzed the landscape of somatic alterations in 60 pediatric patients diagnosed with the most common BCP-ALL subtypes, ETV6::RUNX1+ and classical hyperdiploid (HD), using conventional cytogenetics, single nucleotide polymorphism (SNP) array, whole exome sequencing (WES), and the novel optical genome mapping (OGM) technique. Ninety-five percent of SVs detected by cytogenetics and SNP-array were verified by OGM. OGM detected an additional 677 SVs not identified using the conventional methods, including (subclonal) IKZF1 deletions. Based on OGM, ETV6::RUNX1+ BCP-ALL harbored 2.7 times more SVs than HD BCP-ALL, mainly focal deletions. Besides SVs in known leukemia development genes (ETV6, PAX5, BTG1, CDKN2A), we identified 19 novel recurrently altered regions (in n ≥ 3) including 9p21.3 (FOCAD/HACD4), 8p11.21 (IKBKB), 1p34.3 (ZMYM1), 4q24 (MANBA), 8p23.1 (MSRA), and 10p14 (SFMBT2), as well as ETV6::RUNX1+ subtype-specific SVs (12p13.1 (GPRC5A), 12q24.21 (MED13L), 18q11.2 (MIB1), 20q11.22 (NCOA6)). We detected 3 novel fusion genes (SFMBT2::DGKD, PDS5B::STAG2, and TDRD5::LPCAT2), for which the sequence and expression were validated by long-read and whole transcriptome sequencing, respectively. OGM and WES identified double hits of SVs and SNVs (ETV6, BTG1, STAG2, MANBA, TBL1XR1, NSD2) in the same patient demonstrating the power of the combined approach to define the landscape of genomic alterations in BCP-ALL.

Disruption to the FOXO-PRDM1 axis resulting from deletions of chromosome 6 in acute lymphoblastic leukaemia

A common problem in the study of human malignancy is the elucidation of cancer driver mechanisms associated with recurrent deletion of regions containing multiple genes. Taking B-cell acute lymphoblastic leukaemia (B-ALL) and large deletions of 6q [del(6q)] as a model, we integrated analysis of functional cDNA clone tracking assays with patient genomic and transcriptomic data, to identify the transcription factors FOXO3 and PRDM1 as candidate tumour suppressor genes (TSG). Analysis of cell cycle and transcriptomic changes following overexpression of FOXO3 or PRDM1 indicated that they co-operate to promote cell cycle exit at the pre-B cell stage. FOXO1 abnormalities are absent in B-ALL, but like FOXO3, FOXO1 expression suppressed growth of TCF3::PBX1 and ETV6::RUNX1 B-ALL in-vitro. While both FOXOs induced PRDM1 and other genes contributing to late pre-B cell development, FOXO1 alone induced the key transcription factor, IRF4, and chemokine, CXCR4. CRISPR-Cas9 screening identified FOXO3 as a TSG, while FOXO1 emerged as essential for B-ALL growth. We relate this FOXO3-specific leukaemia-protective role to suppression of glycolysis based on integrated analysis of CRISPR-data and gene sets induced or suppressed by FOXO1 and FOXO3. Pan-FOXO agonist Selinexor induced the glycolysis inhibitor TXNIP and suppressed B-ALL growth at low dose (ID50 < 50 nM).

The Landscape of Secondary Genetic Rearrangements in Pediatric Patients with B-Cell Acute Lymphoblastic Leukemia with t(12;21)

The most frequent chromosomal rearrangement in childhood B-cell acute lymphoblastic leukemia (B-ALL) is translocation t(12;21)(p13;q22). It results in the fusion of the ETV6::RUNX1 gene, which is active in the regulation of multiple crucial cellular pathways. Recent studies hypothesize that many translocations are influenced by RAG-initiated deletions, as well as defects in the RAS and NRAS pathways. According to a "two-hit" model for the molecular pathogenesis of pediatric ETV6::RUNX1-positive B-ALL, the t(12;21) translocation requires leukemia-causing secondary mutations. Patients with ETV6::RUNX1 express up to 60 different aberrations, which highlights the heterogeneity of this B-ALL subtype and is reflected in differences in patient response to treatment and chances of relapse. Most studies of secondary genetic changes have concentrated on deletions of the normal, non-rearranged ETV6 allele. Other predominant structural changes included deletions of chromosomes 6q and 9p, loss of entire chromosomes X, 8, and 13, duplications of chromosome 4q, or trisomy of chromosomes 21 and 16, but the impact of these changes on overall survival remains unclarified. An equally genetically diverse group is the recently identified new B-ALL subtype ETV6::RUNX1-like ALL. In our review, we provide a comprehensive description of recurrent secondary mutations in pediatric B-ALL with t(12;21) to emphasize the value of investigating detailed molecular mechanisms in ETV6::RUNX1-positive B-ALL, both for our understanding of the etiology of the disease and for future clinical advances in patient treatment and management.

A somatic UBA2 variant preceded ETV6-RUNX1 in the concordant BCP-ALL of monozygotic twins

Genetic analysis of leukemic clones in monozygotic twins with concordant acute lymphoblastic leukemia (ALL) has proved a unique opportunity to gain insight into the molecular phylogenetics of leukemogenesis. Using whole-genome sequencing, we characterized constitutional and somatic single nucleotide variants/insertion-deletions (indels) and structural variants in a monozygotic twin pair with concordant ETV6-RUNX1+ B-cell precursor ALL (BCP-ALL). In addition, digital PCR (dPCR) was applied to evaluate the presence of and quantify selected somatic variants at birth, diagnosis, and remission. A shared somatic complex rearrangement involving chromosomes 11, 12, and 21 with identical fusion sequences in leukemias of both twins offered direct proof of a common clonal origin. The ETV6-RUNX1 fusion detected at diagnosis was found to originate from this complex rearrangement. A shared somatic frameshift deletion in UBA2 was also identified in diagnostic samples. In addition, each leukemia independently acquired analogous deletions of 3 genes recurrently targeted in BCP-ALLs (ETV6, ATF7IP, and RAG1/RAG2), providing evidence of a convergent clonal evolution only explained by a strong concurrent selective pressure. Quantification of the UBA2 deletion by dPCR surprisingly indicated it persisted in remission. This, for the first time to our knowledge, provided evidence of a UBA2 variant preceding the well-established initiating event ETV6-RUNX1. Further, we suggest the UBA2 deletion exerted a leukemia predisposing effect and that its essential role in Small Ubiquitin-like Modifier (SUMO) attachment (SUMOylation), regulating nearly all physiological and pathological cellular processes such as DNA-repair by nonhomologous end joining, may hold a mechanistic explanation for the predisposition.

Identification of new ETV6 modulators through a high-throughput functional screening

ETV6 transcriptional activity is critical for proper blood cell development in the bone marrow. Despite the accumulating body of evidence linking ETV6 malfunction to hematological malignancies, its regulatory network remains unclear. To uncover genes that modulate ETV6 repressive transcriptional activity, we performed a specifically designed, unbiased genome-wide shRNA screen in pre-B acute lymphoblastic leukemia cells. Following an extensive validation process, we identified 13 shRNAs inducing overexpression of ETV6 transcriptional target genes. We showed that the silencing of AKIRIN1, COMMD9, DYRK4, JUNB, and SRP72 led to an abrogation of ETV6 repressive activity. We identified critical modulators of the ETV6 function which could participate in cellular transformation through the ETV6 transcriptional network.

•

We develop a genome-wide shRNAs screen for ETV6 modulators

•

The screen uncovered 13 novel putative ETV6 modulator genes

•

The modulators demonstrated a broad impact on the ETV6 transcriptional network

•

T-ALL cells results suggest modulators are conserved in other cellular contexts

Copy Number Changes and Allele Distribution Patterns of Chromosome 21 in B Cell Precursor Acute Lymphoblastic Leukemia

Chromosome 21 is the most affected chromosome in childhood acute lymphoblastic leukemia. Many of its numerical and structural abnormalities define diagnostically and clinically important subgroups. To obtain an overview about their types and their approximate genetic subgroup-specific incidence and distribution, we performed cytogenetic, FISH and array analyses in a total of 578 ALL patients (including 26 with a constitutional trisomy 21). The latter is the preferred method to assess genome-wide large and fine-scale copy number abnormalities (CNA) together with their corresponding allele distribution patterns. We identified a total of 258 cases (49%) with chromosome 21-associated CNA, a number that is perhaps lower-than-expected because ETV6-RUNX1-positive cases (11%) were significantly underrepresented in this array-analyzed cohort. Our most interesting observations relate to hyperdiploid leukemias with tetra- and pentasomies of chromosome 21 that develop in constitutionally trisomic patients. Utilizing comparative short tandem repeat analyses, we were able to prove that switches in the array-derived allele patterns are in fact meiotic recombination sites, which only become evident in patients with inborn trisomies that result from a meiosis 1 error. The detailed analysis of such cases may eventually provide important clues about the respective maldistribution mechanisms and the operative relevance of chromosome 21-specific regions in hyperdiploid leukemias.

Evidence-based review of genomic aberrations in B-lymphoblastic leukemia/lymphoma: Report from the cancer genomics consortium working group for lymphoblastic leukemia

Clinical management and risk stratification of B-lymphoblastic leukemia/ lymphoma (B-ALL/LBL) depend largely on identification of chromosomal abnormalities obtained using conventional cytogenetics and Fluorescence In Situ Hybridization (FISH) testing. In the last few decades, testing algorithms have been implemented to support an optimal risk-oriented therapy, leading to a large improvement in overall survival. In addition, large scale genomic studies have identified multiple aberrations of prognostic significance that are not routinely tested by existing modalities. However, as chromosomal microarray analysis (CMA) and next-generation sequencing (NGS) technologies are increasingly used in clinical management of hematologic malignancies, these abnormalities may be more readily detected. In this article, we have compiled a comprehensive, evidence-based review of the current B-ALL literature, focusing on known and published subtypes described to date. More specifically, we describe the role of various testing modalities in the diagnosis, prognosis, and therapeutic relevance. In addition, we propose a testing algorithm aimed at assisting laboratories in the most effective detection of the underlying genomic abnormalities.

SLX4IP Antagonizes Promiscuous BLM Activity during ALT Maintenance

Cancer cells acquire unlimited proliferative capacity by either re-expressing telomerase or inducing alternative lengthening of telomeres (ALT), which relies on telomere recombination. Here, we show that ALT recombination requires coordinate regulation of the SMX and BTR complexes to ensure the appropriate balance of resolution and dissolution activities at recombining telomeres. Critical to this control is SLX4IP, which accumulates at ALT telomeres and interacts with SLX4, XPF, and BLM. Loss of SLX4IP increases ALT-related phenotypes, which is incompatible with cell growth following concomitant loss of SLX4. Inactivation of BLM is sufficient to rescue telomere aggregation and the synthetic growth defect in this context, suggesting that SLX4IP favors SMX-dependent resolution by antagonizing promiscuous BLM activity during ALT recombination. Finally, we show that SLX4IP is inactivated in a subset of ALT-positive osteosarcomas. Collectively, our findings uncover an SLX4IP-dependent regulatory mechanism critical for telomere maintenance in ALT cancer cells.

Insights of Acute Lymphoblastic Leukemia with Development of Genomic Investigation

Treatment outcomes for acute lymphoblastic leukemia (ALL), especially pediatric ALL, have greatly improved due to the risk-adapted therapy. Combination of drug development, clinical practice, as well as basic genetic researches has brought the survival rate of ALL from less than 10% to more than 90% today, not only increasing the treatment efficacy but also limiting adverse drug reactions (ADRs). In this review, we summarized the landscape identification of ALL genetic alterations, which provided the opportunity to increase the survival rate and especially minimize the relapse risk of ALL, and highlighted the importance of the development of new technologies of genomic investigation for translational medicine.

Genome wide mapping of ETV6 binding sites in pre-B leukemic cells

Genetic alterations in the transcriptional repressor ETV6 are associated with hematological malignancies. Notably, the t(12;21) translocation leading to an ETV6-AML1 fusion gene is the most common genetic alteration found in childhood acute lymphoblastic leukemia. Moreover, most of these patients also lack ETV6 expression, suggesting a tumor suppressor function. To gain insights on ETV6 DNA-binding specificity and genome wide transcriptional regulation capacities, we performed chromatin immunoprecipitation experiments coupled to deep sequencing in a t(12;21)-positive pre-B leukemic cell line. This strategy led to the identification of ETV6-bound regions that were further associated to gene expression. ETV6 binding is mostly cell type-specific as only few regions are shared with other blood cell subtypes. Peaks localization and motif enrichment analyses revealed that this unique binding profile could be associated with the ETV6-AML1 fusion protein specific to the t(12;21) background. This study underscores the complexity of ETV6 binding and uncovers ETV6 transcriptional network in pre-B leukemia cells bearing the recurrent t(12;21) translocation.

Translocations involving ETS family proteins in human cancer

The ETS transcription factors regulate expression of genes involved in normal cell development, proliferation, differentiation, angiogenesis, and apoptosis, consisting of 28 family members in humans. Dysregulation of these transcription factors facilitates cell proliferation in cancers, and several members participate in invasion and metastasis by activating certain gene transcriptions. ETS1 and ETS2 are the founding members of the ETS family and regulate transcription by binding to ETS sequences. Three chimeric genes involving ETS genes have been identified in human cancers, which are EWS-FLI1 in Ewing's sarcoma, TMPRSS2-ERG in prostate cancer, and ETV6-RUNX1 in acute lymphocytic leukemia. Although these fusion transcripts definitely contribute to the pathogenesis of the disease, the impact of these fusion transcripts on patients' prognosis is highly controversial. In the present review, the roles of ETS protein translocations in human carcinogenesis are discussed.

ETV6/RUNX1-positive childhood acute lymphoblastic leukemia (ALL): The spectrum of clonal heterogeneity and its impact on prognosis

The prognostic significance of the ETV6/RUNX1-fusion and of the accompanying aberrations is disputable; whether co-existing sub-clones are responsible for delayed MRD-clearance and thus, moderate outcome, remains to be clarified. We studied, in a paediatric cohort of 119 B-ALLs, the relation between the ETV6/RUNX1 aberration and the co-existing subclones with (a) presenting clinical/biological features, (b) early response to treatment(MRD) and (c) long-term outcome over a 12-year period. Patients were homogeneously treated according to BFM-based-protocols. 27/119 patients (22.7%) were ETV6/RUNX1-positive; 19/27 (70.4%) harbored additional genetic abnormalities while 9/19 (33.3%) presented with clonal heterogeneity. The most common abnormalities were del12p13 (37%), 3-6×21q22 (22.2%), del9p21 (18.5%) and 2-3xETV6/RUNX1 (18.5%). MRDd15-positivity (≥10-3) was detected in 44% of the cohort; the corresponding MRD among patients carrying subclones rises to 88.9%. Common features of all relapses were sub-clonal diversity, FCM-MRDd15-positivity and additional del(9p21) while there were no censored relapses among ETV6/RUNX1-positive patients with sole translocation and absence of additional aberrations, within a median follow-up time of 90 months. In our study, the presence of clonal heterogeneity and impaired FCM-MRD clearance among ETV6/RUNX1-positive patients, ultimately influenced prognosis. Longer follow-up is needed in order to further validate these initial results.

Pathogenesis of ETV6/RUNX1-positive childhood acute lymphoblastic leukemia and mechanisms underlying its relapse

ETV6/RUNX1 (E/R) is the most common fusion gene in childhood acute lymphoblastic leukemia (ALL). Multiple lines of evidence imply a “two-hit” model for the molecular pathogenesis of E/R-positive ALL, whereby E/R rearrangement is followed by a series of secondary mutations that trigger overt leukemia. The cellular framework in which E/R arises and the maintenance of a pre-leukemic condition by E/R are fundamental to the mechanism that underlies leukemogenesis. Accordingly, a variety of studies have focused on the relationship between the clones giving rise to the primary and recurrent E/R-positive ALL. We review here the most recent insights into the pathogenic mechanisms underlying E/R-positive ALL, as well as the molecular abnormalities prevailing at relapse.

RUNX-mediated growth arrest and senescence are attenuated by diverse mechanisms in cells expressing RUNX1 fusion oncoproteins

RUNX gene over‐expression inhibits growth of primary cells but transforms cells with tumor suppressor defects, consistent with reported associations with tumor progression. In contrast, chromosomal translocations involving RUNX1 are detectable in utero, suggesting an initiating role in leukemias. How do cells expressing RUNX1 fusion oncoproteins evade RUNX‐mediated growth suppression? Previous studies showed that the TEL‐RUNX1 fusion from t(12;21) B‐ALLs is unable to induce senescence‐like growth arrest (SLGA) in primary fibroblasts while potent activity is displayed by the RUNX1‐ETO fusion found in t(8;21) AMLs. We now show that SLGA potential is suppressed in TEL‐RUNX1 but reactivated by deletion of the TEL HLH domain or mutation of a key residue (K99R). Attenuation of SLGA activity is also a feature of RUNX1‐ETO9a, a minor product of t(8;21) translocations with increased leukemogenicity. Finally, while RUNX1‐ETO induces SLGA it also drives a potent senescence‐associated secretory phenotype (SASP), and promotes the immortalization of rare cells that escape SLGA. Moreover, the RUNX1‐ETO SASP is not strictly linked to growth arrest as it is largely suppressed by RUNX1 and partially activated by RUNX1‐ETO9a. These findings underline the heterogeneous nature of premature senescence and the multiple mechanisms by which this failsafe process is subverted in cells expressing RUNX1 oncoproteins.

Whole-Exome Sequencing of ETV6/RUNX1 in Four Childhood Acute Lymphoblastic Leukaemia Cases

Background:

ETV6/RUNX1 gene fusion is the most frequently seen chromosomal abnormality in childhood acute lymphobastic leukamia (ALL). However, additional genetic changes are known to be required for the development of this type of leukaemia. Therefore, we here aimed to assess the somatic mutational profile of four ALL cases carrying the ETV6/RUNX1 fusion gene using whole-exome sequencing.

Methods:

DNA was isolated from bone marrow samples using a QIAmp DNA Blood Mini kit and subsequently sequenced using the Illumina MiSeq system.

Results:

We identified 12,960 to17,601 mutations in each sample, with a total of 16,466 somatic mutations in total. Some 15,533 variants were single nucleotide polymorphisms (SNPs), 129 were substitutions, 415 were insertions and 389 were deletions. When taking into account the coding region and protein impact, 1,875 variants were synonymous and 1,956 were non-synonymous SNPs. Among non-synonymous SNPs, 1,862 were missense, 13 nonsense, 35 frameshifts, 11 nonstop, 3 misstart, 15 splices disrupt and 17 in-frame indels. A total of 86 variants were located in leukaemia-related genes of which 32 variants were located in the coding regions of GLI2, SP140, GATA2, SMAD5, KMT2C, CDH17, CDX2, FLT3, PML and MOV10L1.

Conclusions:

Detection and identification of secondary genetic alterations are important in identifying new therapeutic targets and developing rationally designed treatment regimens with less toxicity in ALL patients.

Role of RUNX1 in hematological malignancies

RUNX1 is a member of the core-binding factor family of transcription factors and is indispensable for the establishment of definitive hematopoiesis in vertebrates. RUNX1 is one of the most frequently mutated genes in a variety of hematological malignancies. Germ line mutations in RUNX1 cause familial platelet disorder with associated myeloid malignancies. Somatic mutations and chromosomal rearrangements involving RUNX1 are frequently observed in myelodysplastic syndrome and leukemias of myeloid and lymphoid lineages, that is, acute myeloid leukemia, acute lymphoblastic leukemia, and chronic myelomonocytic leukemia. More recent studies suggest that the wild-type RUNX1 is required for growth and survival of certain types of leukemia cells. The purpose of this review is to discuss the current status of our understanding about the role of RUNX1 in hematological malignancies.

Mechanism of ETV6-RUNX1 Leukemia

The t(12;21)(p13;q22) translocation is the most frequently occurring single genetic abnormality in pediatric leukemia. This translocation results in the fusion of the ETV6 and RUNX1 genes. Since its discovery in the 1990s, the function of the ETV6-RUNX1 fusion gene has attracted intense interest. In this chapter, we will summarize current knowledge on the clinical significance of ETV6-RUNX1, the experimental models used to unravel its function in leukemogenesis, the identification of co-operating mutations and the mechanisms responsible for their acquisition, the function of the encoded transcription factor and finally, the future therapeutic approaches available to mitigate the associated disease.

The RUNX Genes as Conditional Oncogenes: Insights from Retroviral Targeting and Mouse Models

The observation that the Runx genes act as targets for transcriptional activation by retroviral insertion identified a new family of dominant oncogenes. However, it is now clear that Runx genes are 'conditional' oncogenes whose over-expression is growth inhibitory unless accompanied by another event such as concomitant over-expression of MYC or loss of p53 function. Remarkably, while the oncogenic activities of either MYC or RUNX over-expression are suppressed while p53 is intact, the combination of both neutralises p53 tumour suppression in vivo by as yet unknown mechanisms. Moreover, there is emerging evidence that endogenous, basal RUNX activity is important to maintain the viability and proliferation of MYC-driven lymphoma cells. There is also growing evidence that the human RUNX genes play a similar conditional oncogenic role and are selected for over-expression in end-stage cancers of multiple types. Paradoxically, reduced RUNX activity can also predispose to cell immortalisation and transformation, particularly by mutant Ras. These apparently conflicting observations may be reconciled in a stage-specific model of RUNX involvement in cancer. A question that has yet to be fully addressed is the extent to which the three Runx genes are functionally redundant in cancer promotion and suppression.

Characterization of a novel acquired der(1)del(1)(p13p31)t(1;15)(q42;q15) in a high risk t(12;21)-positive acute lymphoblastic leukemia

The t(12;21)(p13;q22) with ETV6-RUNX1 fusion occurs in 25% of cases of B-cell precursor acute lymphoblastic leukemia (BCP-ALL); and is generally associated with favorable prognosis. However, 15-20% of the t(12;21)-positive cases are associated with high-risk disease due to for example slow early responses to therapy. It is well-known that development of overt leukemia in t(12;21)-positive ALL requires secondary chromosomal aberrations although the full spectrum of these cytogenetic alterations is yet unsettled, and also, how they may be associated with disease outcome. This report describes the case of an adolescent male with t(12;21)-positive ALL who displayed a G-banded karyotype initially interpreted as del(1)(p22p13) and del(15)(q15). The patient was treated according to NOPHO standard risk protocol at diagnosis, but had minimal residual disease (MRD) at 6,4% on day 29 as determined by flow cytometric immunophenotyping. Because of MRD level>0.1% he was then assigned as a high risk patient and received intensified chemotherapy accordingly. Further molecular cytogenetic studies and oligo-based aCGH (oaCGH) analysis characterized the acquired complex structural rearrangements on chromosomes 1 and 15, which can be described as der(1)del(1)(p13.1p31.1)t(1;15)(q42;q15) with concurrent deletions at 1q31.2-q31.3, 1q42.12-q43, and 15q15.1-q15.3. The unbalanced complex rearrangements have not been described previously. Extended locus-specific FISH analyses showed that the three deletions were on the same chromosome 1 homologue that was involved in the t(1;15), and that the deletion on chromosome 15 also was on the same chromosome 15 homologue as involved in the t(1;15). Together these findings show the great importance of the combined usage of molecular cytogenetic analyses and oaCGH analysis to enhance characterization of apparently simple G-banded karyotypes, and to provide a more complete spectrum of secondary chromosomal aberrations in high risk t(12;21)-positive BCP-ALLs.

CLIC5: a novel ETV6 target gene in childhood acute lymphoblastic leukemia

The most common rearrangement in childhood precursor B-cell acute lymphoblastic leukemia is the t(12;21)(p13;q22) translocation resulting in the ETV6-AML1 fusion gene. A frequent concomitant event is the loss of the residual ETV6 allele suggesting a critical role for the ETV6 transcriptional repressor in the etiology of this cancer. However, the precise mechanism through which loss of functional ETV6 contributes to disease pathogenesis is still unclear. To investigate the impact of ETV6 loss on the transcriptional network and to identify new transcriptional targets of ETV6, we used whole transcriptome analysis of both pre-B leukemic cell lines and patients combined with chromatin immunoprecipitation. Using this integrative approach, we identified 4 novel direct ETV6 target genes: CLIC5, BIRC7, ANGPTL2 and WBP1L To further evaluate the role of chloride intracellular channel protein CLIC5 in leukemogenesis, we generated cell lines overexpressing CLIC5 and demonstrated an increased resistance to hydrogen peroxide-induced apoptosis. We further described the implications of CLIC5's ion channel activity in lysosomal-mediated cell death, possibly by modulating the function of the transferrin receptor with which it colocalizes intracellularly. For the first time, we showed that loss of ETV6 leads to significant overexpression of CLIC5, which in turn leads to decreased lysosome-mediated apoptosis. Our data suggest that heightened CLIC5 activity could promote a permissive environment for oxidative stress-induced DNA damage accumulation, and thereby contribute to leukemogenesis.

Genetic alterations in glucocorticoid signaling pathway components are associated with adverse prognosis in children with relapsed ETV6/RUNX1-positive acute lymphoblastic leukemia

The ETV6/RUNX1 gene fusion defines the largest genetic subgroup of childhood ALL with overall rapid treatment response. However, up to 15% of cases relapse. Because an impaired glucocorticoid pathway is implicated in disease recurrence we studied the impact of genetic alterations by SNP array analysis in 31 relapsed cases. In 58% of samples, we found deletions in various glucocorticoid signaling pathway-associated genes, but only NR3C1 and ETV6 deletions prevailed in minimal residual disease poor responding and subsequently relapsing cases (p<0.05). To prove the necessity of a functional glucocorticoid receptor, we reconstituted wild-type NR3C1 expression in mutant, glucocorticoid-resistant REH cells and studied the glucocorticoid response in vitro and in a xenograft mouse model. While these results prove that glucocorticoid receptor defects are crucial for glucocorticoid resistance in an experimental setting, they do not address the essential clinical situation where glucocorticoid resistance at relapse is rather part of a global drug resistance.

The mutational landscape in pediatric acute lymphoblastic leukemia deciphered by whole genome sequencing

Genomic characterization of pediatric acute lymphoblastic leukemia (ALL) has identified distinct patterns of genes and pathways altered in patients with well-defined genetic aberrations. To extend the spectrum of known somatic variants in ALL, we performed whole genome and transcriptome sequencing of three B-cell precursor patients, of which one carried the t(12;21)ETV6-RUNX1 translocation and two lacked a known primary genetic aberration, and one T-ALL patient. We found that each patient had a unique genome, with a combination of well-known and previously undetected genomic aberrations. By targeted sequencing in 168 patients, we identified KMT2D and KIF1B as novel putative driver genes. We also identified a putative regulatory non-coding variant that coincided with overexpression of the growth factor MDK. Our results contribute to an increased understanding of the biological mechanisms that lead to ALL and suggest that regulatory variants may be more important for cancer development than recognized to date. The heterogeneity of the genetic aberrations in ALL renders whole genome sequencing particularly well suited for analysis of somatic variants in both research and diagnostic applications.

Somatic drivers of B-ALL in a model of ETV6-RUNX1; Pax5(+/-) leukemia

BACKGROUND

B-cell precursor acute lymphoblastic leukemia (B-ALL) is amongst the leading causes of childhood cancer-related mortality. Its most common chromosomal aberration is the ETV6-RUNX1 fusion gene, with ~25% of ETV6-RUNX1 patients also carrying PAX5 alterations.

METHODS

We have recreated this mutation background by inter-crossing Etv6-RUNX1 (Etv6 (RUNX1-SB)) and Pax5(+/-) mice and performed an in vivo analysis to find driver genes using Sleeping Beauty transposon-mediated mutagenesis and also exome sequencing.

RESULTS

Combination of Etv6-RUNX1 and Pax5(+/-) alleles generated a transplantable B220 + CD19+ B-ALL with a significant disease incidence. RNA-seq analysis showed a gene expression pattern consistent with arrest at the pre-B stage. Analysis of the transposon common insertion sites identified genes involved in B-cell development (Zfp423) and the JAK/STAT signaling pathway (Jak1, Stat5 and Il2rb), while exome sequencing revealed somatic hotspot mutations in Jak1 and Jak3 at residues analogous to those mutated in human leukemias, and also mutation of Trp53.

CONCLUSIONS

Powerful synergies exists in our model suggesting STAT pathway activation and mutation of Trp53 are potent drivers of B-ALL in the context of Etv6-RUNX1;Pax5(+/-).

CD200/BTLA deletions in pediatric precursor B-cell acute lymphoblastic leukemia treated according to the EORTC-CLG 58951 protocol

DNA copy number analysis has been instrumental for the identification of genetic alterations in B-cell precursor acute lymphoblastic leukemia. Notably, some of these genetic defects have been associated with poor treatment outcome and might be relevant for future risk stratification. In this study, we characterized recurrent deletions of CD200 and BTLA genes, mediated by recombination-activating genes, and used breakpoint-specific polymerase chain reaction assay to screen a cohort of 1154 cases of B-cell precursor acute lymphoblastic leukemia uniformly treated according to the EORTC-CLG 58951 protocol. CD200/BTLA deletions were identified in 56 of the patients (4.8%) and were associated with an inferior 8-year event free survival in this treatment protocol [70.2% ± 1.2% for patients with deletions versus 83.5% ± 6.4% for non-deleted cases (hazard ratio 2.02; 95% confidence interval 1.23-3.32; P=0.005)]. Genetically, CD200/BTLA deletions were strongly associated with ETV6-RUNX1-positive leukemias (P<0.0001), but were also identified in patients who did not have any genetic abnormality that is currently used for risk stratification. Within the latter population of patients, the presence of CD200/BTLA deletions was associated with inferior event-free survival and overall survival. Moreover, the multivariate Cox model indicated that these deletions had independent prognostic impact on event-free survival when adjusting for conventional risk criteria. All together, these findings further underscore the rationale for copy number profiling as an important tool for risk stratification in human B-cell precursor acute lymphoblastic leukemia. This trial was registered at www.ClinicalTrials.gov as #NCT00003728.

ETV6/RUNX1 induces reactive oxygen species and drives the accumulation of DNA damage in B cells

The t(12;21)(p13;q22) chromosomal translocation is the most frequent translocation in childhood B cell precursor-acute lymphoblastic leukemia and results in the expression of an ETV6/RUNX1 fusion protein. The frequency of ETV6/RUNX1 fusions in newborns clearly exceeds the leukemia rate revealing that additional events occur in ETV6/RUNX1-positive cells for leukemic transformation. Hitherto, the mechanisms triggering these second hits remain largely elusive. Thus, we generated a novel ETV6/RUNX1 transgenic mouse model where the expression of the fusion protein is restricted to CD19(+) B cells. These animals harbor regular B cell development and lack gross abnormalities. We established stable pro-B cell lines carrying the ETV6/RUNX1 transgene that allowed us to investigate whether ETV6/RUNX1 itself favors the acquisition of second hits. Remarkably, these pro-B cell lines as well as primary bone marrow cells derived from ETV6/RUNX1 transgenic animals display elevated levels of reactive oxygen species (ROS) as tested with ETV6/RUNX1 transgenic dihydroethidium staining. In line, intracellular phospho-histone H2AX flow cytometry and comet assay revealed increased DNA damage indicating that ETV6/RUNX1 expression enhances ROS. On the basis of our data, we propose the following model: the expression of ETV6/RUNX1 creates a preleukemic clone and leads to increased ROS levels. These elevated ROS favor the accumulation of secondary hits by increasing genetic instability and double-strand breaks, thus allowing preleukemic clones to develop into fully transformed leukemic cells.

A new workflow for whole-genome sequencing of single human cells

Unbiased amplification of the whole-genome amplification (WGA) of single cells is crucial to study cancer evolution and genetic heterogeneity, but is challenging due to the high complexity of the human genome. Here, we present a new workflow combining an efficient adapter-linker PCR-based WGA method with second-generation sequencing. This approach allows comparison of single cells at base pair resolution. Amplification recovered up to 74% of the human genome. Copy-number variants and loss of heterozygosity detected in single cell genomes showed concordance of up to 99% to pooled genomic DNA. Allele frequencies of mutations could be determined accurately due to an allele dropout rate of only 2%, clearly demonstrating the low bias of our PCR-based WGA approach. Sequencing with paired-end reads allowed genome-wide analysis of structural variants. By direct comparison to other WGA methods, we further endorse its suitability to analyze genetic heterogeneity.

Blocking ETV6/RUNX1-induced MDM2 overexpression by Nutlin-3 reactivates p53 signaling in childhood leukemia

ETV6/RUNX1 (E/R) is the most common fusion gene in childhood acute lymphoblastic leukemia. It is responsible for the initiation of leukemia but also indispensable for disease maintenance and propagation, although its function in these latter processes is less clear. We therefore investigated the effects of the perceived p53 pathway alterations in model cell lines and primary leukemias and, in particular, how E/R upregulates MDM2, the predominant negative regulator of p53. We found that E/R transactivates MDM2 in both p53(+/+) and p53(-/-) HCT116 cells by binding to promoter-inherent RUNX1 motifs, which indicates that this activation occurs in a direct and p53-independent manner. Treatment of E/R-positive leukemic cell lines with Nutlin-3, a small molecule that inhibits the MDM2/p53 interaction, arrests their cell cycle and induces apoptosis. These phenomena concur with a p53-induced expression of p21, pro-apoptotic BAX and PUMA, as well as caspase 3 activation and poly ADP-ribose polymerase cleavage. The addition of DNA-damaging and p53-activating chemotherapeutic drugs intensifies apoptosis. Moreover, Nutlin-3 exposure leads to an analogous p53 accumulation and apoptotic surge in E/R-positive primary leukemic cells. Our findings clarify the role of p53 signaling in E/R-positive leukemias and outline the potential basis for its therapeutic exploitation in this setting.

Frequent and sex-biased deletion of SLX4IP by illegitimate V(D)J-mediated recombination in childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) accounts for ∼25% of pediatric malignancies. Of interest, the incidence of ALL is observed ∼20% higher in males relative to females. The mechanism behind the phenomenon of sex-specific differences is presently not understood. Employing genome-wide genetic aberration screening in 19 ALL samples, one of the most recurrent lesions identified was monoallelic deletion of the 5' region of SLX4IP. We characterized this deletion by conventional molecular genetic techniques and analyzed its interrelationships with biological and clinical characteristics using specimens and data from 993 pediatric patients enrolled into trial AIEOP-BFM ALL 2000. Deletion of SLX4IP was detected in ∼30% of patients. Breakpoints within SLX4IP were defined to recurrent positions and revealed junctions with typical characteristics of illegitimate V(D)J-mediated recombination. In initial and validation analyses, SLX4IP deletions were significantly associated with male gender and ETV6/RUNX1-rearranged ALL (both overall P < 0.0001). For mechanistic validation, a second recurrent deletion affecting TAL1 and caused by the same molecular mechanism was analyzed in 1149 T-cell ALL patients. Validating a differential role by sex of illegitimate V(D)J-mediated recombination at the TAL1 locus, 128 out of 1149 T-cell ALL samples bore a deletion and males were significantly more often affected (P = 0.002). The repeatedly detected association of SLX4IP deletion with male sex and the extension of the sex bias to deletion of the TAL1 locus suggest that differential illegitimate V(D)J-mediated recombination events at specific loci may contribute to the consistent observation of higher incidence rates of childhood ALL in boys compared with girls.

Copy number genome alterations are associated with treatment response and outcome in relapsed childhood ETV6/RUNX1-positive acute lymphoblastic leukemia

The clinical heterogeneity among first relapses of childhood ETV6/RUNX1-positive acute lymphoblastic leukemia indicates that further genetic alterations in leukemic cells might affect the course of salvage therapy and be of prognostic relevance. To assess the incidence and prognostic relevance of additional copy number alterations at relapse of the disease, we performed whole genome array comparative genomic hybridization of leukemic cell DNA from 51 patients with first ETV6/RUNX1-positive relapse enrolled in and treated according to the relapse trials ALL-REZ of the Berlin-Frankfurt-Münster Study Group. Within this cohort of patients with relapsed ETV6/RUNX1-positive acute lymphoblastic leukemia, the largest analyzed for genome wide DNA copy number alterations to date, alterations were present in every ETV6/RUNX1-positive relapse and a high proportion of them occurred in recurrent overlapping chromosomal regions. Recurrent losses affected chromosomal regions 12p13, 6q21, 15q15.1, 9p21, 3p21, 5q and 3p14.2, whereas gains occurred in regions 21q22 and 12p. Loss of 12p13 including CDKN1B was associated with a shorter remission duration (P=0.009) and a lower probability of event-free survival (P=0.001). Distribution of X-chromosomal copy number alterations was gender-specific: whole X-chromosome loss occurred exclusively in females, gain of Xq only in males. Loss of the glucocorticoid receptor gene NR3C1 (5q31.3) was associated with a poor response to induction treatment (P=0.003), possibly accounting for the adverse prognosis of some of the ETV6/RUNX1-positive relapses.

Genome-wide analysis of cytogenetic aberrations in ETV6/RUNX1-positive childhood acute lymphoblastic leukaemia

The chromosomal translocation t(12;21) resulting in the ETV6/RUNX1 fusion gene is the most frequent structural cytogenetic abnormality among patients with childhood acute lymphoblastic leukaemia (ALL). We investigated 62 ETV6/RUNX1-positive childhood ALL patients by single nucleotide polymorphism array to explore acquired copy number alterations (CNAs) at diagnosis. The mean number of CNAs was 2·82 (range 0-14). Concordance with available G-band karyotyping and comparative genomic hybridization was 93%. Based on three major protein-protein complexes disrupted by these CNAs, patients could be categorized into four distinct subgroups, defined by different underlying biological mechanisms relevant to the aetiology of childhood ALL. When recurrent CNAs were evaluated by an oncogenetic tree analysis classifying their sequential order, the most common genetic aberrations (deletions of 6q, 9p, 13q and X, and gains of 10 and 21) seemed independent of each other. Finally, we identified the most common regions with recurrent gains and losses, which comprise microRNA clusters with known oncogenic or tumour-suppressive roles. The present study sheds further light on the genetic diversity of ETV6/RUNX1-positive childhood ALL, which may be important for understanding poor responses among this otherwise highly curable subset of ALL and lead to novel targeted treatment strategies.

Whole-exome sequencing of pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL), the most common malignant disorder in childhood, is typically associated with numerical chromosomal aberrations, fusion genes or small focal deletions, thought to represent important pathogenetic events in the development of the leukemia. Mutations, such as single nucleotide changes, have also been reported in childhood ALL, but these have only been studied by sequencing a small number of candidate genes. Herein, we report the first unbiased sequencing of the whole exome of two cases of pediatric ALL carrying the ETV6/RUNX1 (TEL/AML1) fusion gene (the most common genetic subtype) and corresponding normal samples. A total of 14 somatic mutations were identified, including four and seven protein-altering nucleotide substitutions in each ALL. Twelve mutations (86%) occurred in genes previously described to be mutated in other types of cancer, but none was found to be recurrent in an extended series of 29 ETV6/RUNX1-positive ALLs. The number of single nucleotide mutations was similar to the number of copy number alterations as detected by single nucleotide polymorphism arrays. Although the true pathogenetic significance of the mutations must await future functional evaluations, this study provides a first estimate of the mutational burden at the genetic level of t(12;21)-positive childhood ALL.

High frequency of BTG1 deletions in acute lymphoblastic leukemia in children with down syndrome

Previous cytogenetic studies of myeloid and acute lymphoblastic leukemias in children with Down syndrome (ML-DS and DS-ALL) have revealed significant differences in abnormality patterns between such cases and acute leukemias in general. Also, certain molecular genetic aberrations characterize DS-related leukemias, such as GATA1 mutations in ML-DS and deregulation of the CRLF2 gene in DS-ALL. Whether microdeletions/microduplications also vary between DS and non-DS cases is presently unclear. To address this issue, we performed single nucleotide polymorphism array analyses of eight pediatric ML-DS and 17 B-cell precursor DS-ALL. In the ML-DS cases, a total of 29 imbalances (20 gains and nine losses) and two partial uniparental isodisomies (pUPDs) were detected. None of the 11 small (defined as <10 Mb) imbalances were recurrent, nor were the pUPDs, whereas of the 18 large aberrations, three were recurrent-dup(1q), +8 and +21. In contrast, several frequent changes were identified in the DS-ALL cases, which harbored 82 imbalances (30 gains and 52 losses) and four pUPDs. Of the 40 large changes, 28 were gains and 12 losses, with +X, dup(Xq), dup(1q), del(7p), dup(8q), del(9p), dup(9p), del(12p), dup(17q), and +21 being recurrent. Of the 40 microdeletions identified, several targeted specific genes, with the following being repeatedly deleted: BTG1 and CDKN2A/B (29% of cases), ETV6, IKZF1, PAX5 and SERP2 (18%), and BTLA, INPP4B, P2RY8, and RB1 (12%). Loss of the SERP2 and INPP4B genes, encoding the stress-associated endoplasmic reticulum protein family member 2 and the inositol polyphosphate 4-phosphatase-II, respectively, has previously never been implicated in leukemia. Although deletions of the other genes have been associated with ALL, the high frequency of BTG1 loss is a novel finding. Such deletions may characterize a clinical subgroup of DS-ALL, comprising mainly boys with a high median age. In conclusion, ML-DS and DS-ALL are genetically distinct, with mainly gains in ML-DS and deletions in DS-ALL. Furthermore, DS-ALL is characterized by several recurrent gene deletions, with BTG1 loss being particularly frequent.

Modeling the evolution of ETV6-RUNX1-induced B-cell precursor acute lymphoblastic leukemia in mice

The t(12;21) translocation that generates the ETV6-RUNX1 (TEL-AML1) fusion gene, is the most common chromosomal rearrangement in childhood cancer and is exclusively associated with B-cell precursor acute lymphoblastic leukemia (BCP-ALL). The translocation arises in utero and is necessary but insufficient for the development of leukemia. Single-nucleotide polymorphism array analysis of ETV6-RUNX1 patient samples has identified multiple additional genetic alterations; however, the role of these lesions in leukemogenesis remains undetermined. Moreover, murine models of ETV6-RUNX1 ALL that faithfully recapitulate the human disease are lacking. To identify novel genes that cooperate with ETV6-RUNX1 in leukemogenesis, we generated a mouse model that uses the endogenous Etv6 locus to coexpress the Etv6-RUNX1 fusion and Sleeping Beauty transposase. An insertional mutagenesis screen was performed by intercrossing these mice with those carrying a Sleeping Beauty transposon array. In contrast to previous models, a substantial proportion (20%) of the offspring developed BCP-ALL. Isolation of the transposon insertion sites identified genes known to be associated with BCP-ALL, including Ebf1 and Epor, in addition to other novel candidates. This is the first mouse model of ETV6-RUNX1 to develop BCP-ALL and provides important insight into the cooperating genetic alterations in ETV6-RUNX1 leukemia.

Clonal origins of relapse in ETV6-RUNX1 acute lymphoblastic leukemia

B-cell precursor childhood acute lymphoblastic leukemia with ETV6-RUNX1 (TEL-AML1) fusion has an overall good prognosis, but relapses occur, usually after cessation of treatment and occasionally many years later. We have investigated the clonal origins of relapse by comparing the profiles of genomewide copy number alterations at presentation in 21 patients with those in matched relapse (12-119 months). We identified, in total, 159 copy number alterations at presentation and 231 at relapse (excluding Ig/TCR). Deletions of CDKN2A/B or CCNC (6q16.2-3) or both increased from 38% at presentation to 76% in relapse, suggesting that cell-cycle deregulation contributed to emergence of relapse. A novel observation was recurrent gain of chromosome 16 (2 patients at presentation, 4 at relapse) and deletion of plasmocytoma variant translocation 1 in 3 patients. The data indicate that, irrespective of time to relapse, the relapse clone was derived from either a major or minor clone at presentation. Backtracking analysis by FISH identified a minor subclone at diagnosis whose genotype matched that observed in relapse ∼ 10 years later. These data indicate subclonal diversity at diagnosis, providing a variable basis for intraclonal origins of relapse and extended periods (years) of dormancy, possibly by quiescence, for stem cells in ETV6-RUNX1(+) acute lymphoblastic leukemia.

New strategies in acute lymphoblastic leukemia: translating advances in genomics into clinical practice

B-precursor acute lymphoblastic leukemia (B-ALL) is the most common childhood malignancy and remains a leading cause of death in children and young adults. Current therapeutic approaches involve intensive combination chemotherapy, which fails in up to one quarter of patients. New treatment approaches directed against rational therapeutic targets are required. Recent genomic profiling of ALL has identified several genetic alterations associated with a high risk of treatment failure. Deletion or sequence mutation of the lymphoid transcription factor gene IKZF1 (IKAROS) is associated with a high rate of leukemic relapse, and testing for IKZF1 alterations at diagnosis may aid risk stratification. A subset of B-ALL patients with IKZF1 alterations have a transcriptional profile similar to BCR-ABL1-positive ALL, and these patients commonly have novel rearrangements and mutations resulting in aberrant cytokine receptor signaling and activation of kinase signaling cascades, including rearrangement of CRLF2 and activating mutations of Janus kinases (JAK1 and JAK2). JAK inhibitor therapy is under investigation in children with relapsed and refractory malignancies, including leukemia.