Letters to Blood To the editor: RNA-sequencing analysis of core binding factor AML identifies recurrent ZBTB7A mutations and defines RUNX1-CBFA2T3 fusion signature Vincent-Philippe Lavall ´ ee, 1,2 S´ ebastien Lemieux, 1,3 Genevi ` eve Boucher, 1 Patrick Gendron, 1 Isabel Boivin, 1 Richard N. Armstrong, 1 Guy Sauvageau, 1,2,4,5 and Jos´ ee H´ ebert 1,2,4,5 1 The Leucegene project at Institute for Research in Immunology and Cancer, Universit ´ e de Montr´ eal, Montr´ eal, QC, Canada; 2 Division of Hematology, Maisonneuve-Rosemont Hospital, Montr ´ eal, QC, Canada; 3 Department of Computer Science and Operations Research, Universit ´ e de Montr´ eal, Montr ´ eal, QC, Canada; 4 Quebec Leukemia Cell Bank, Maisonneuve-Rosemont Hospital, Montr ´ eal, QC, Canada; and 5 Department of Medicine, Faculty of Medicine, Universit ´ e de Montr´ eal, Montr ´ eal, QC, Canada RUNX1 (also known as AML1 or CBFA2) and CBFB encode the a and b subunits of a heterodimeric core binding transcription factor complex involved in the development of normal hematopoiesis (reviewed by de Bruijn and Speck 1 ). Both genes are rearranged in acute myeloid leukemia (AML) with t(8;21)(q22;q22);RUNX1-RUNX1T1 and inv(16) (p13.1q22)/t(16;16)(p13.1;q22);CBFB-MYH11, which are collectively called core binding factor (CBF) AML. Targeted mutational studies have shown that both CBF AML subgroups are characterized by recurrent mutations in KIT, FLT3, NRAS, and KRAS. 2-5 Additional mutations in ASXL2 and ASXL1 were also recently described in t(8;21) AML. 6,7 Mutation analyses employing untargeted approaches in large CBF AML cohorts are still lacking. We previously described comparative transcriptomic approaches leading to the comprehensive description of the mutational and transcriptomic landscape of MLL, 8 EVI1, 9 and NUP98-NSD1 10 AML subgroups. Using this same methodology and cohort, we now report the results of the 48 CBF AML specimens in- cluded in our collection of 415 specimens. We identified novel mutations and differentially expressed genes and demonstrated that the RUNX1- CBFA2T3 AML sample is characterized by a gene expression profile highly similar to that of RUNX1-RUNX1T1 AML. This study is part of the Leucegene project, an initiative approved by the research ethics boards of Universit´ e de Montr´ eal and Maisonneuve- Rosemont Hospital. All AML samples were collected with an informed consent between 2001 and 2015 according to Quebec Leukemia Cell Bank (BCLQ) procedures. Workflow for sequencing, mutation analysis, and transcripts quantification have been described previously 8-10 and are complemented in the supplemental Methods (available on the Blood Web site). Forty-eight CBF AML specimens, including 28 samples with inv(16) and 20 with t(8;21), and 367 control AML specimens were part of this analysis. Patient characteristics are described in Figure 1A. The cytogenetic distribution of the entire cohort is shown in supplemental Figure 1. Using the most minimally and differentially expressed genes, we identified signatures of 145 and 127 genes that best characterize t(8;21) and inv(16) subgroups, respectively (Figure 1B-C; supplemental Tables 1 and 2). Fusion partner genes, RUNX1T1 and MYH11, are among the single most differentially expressed genes in each corresponding group. Previously reported candidates such as POU4F1 [t(8;21)] and ST18 [inv(16)] were among the most discriminatory genes identified by our analysis. Other CBF microarray data sets were readily enriched in Gene Set Enrichment Analysis studies 11-13 (supplemental Tables 1 and 2). Importantly, ;80% of genes identified in our CBF AML signatures have not been previously described in those data sets. For example, ADARB2-AS1 and LINC00958 are typical for t(8;21) AML and MEGF10 and APLN for inv(16) specimens in our collection. Our signatures shared 50% and 25% of the most significantly overex- pressed or underexpressed genes with pediatric t(8;21) and inv(16) AML cohorts, respectively (supplemental Tables 1 and 2). 14 This suggests that similar networks are at play in pediatric and adult CBF AML. Using the subgroup-specific gene signatures and performing principal component analyses (PCAs), each CBF subgroup homoge- neously clustered together (Figure 1D-E). Most interestingly, 1 sample harboring a t(16;21)(q24;q22);RUNX1-CBFA2T3 unambiguously clustered with t(8;21) specimens, suggesting that the transcriptional network is shared between these 2 entities (Figure 1D; supplemental Figure 2). In agreement with this observation, a KIT D817V mutation was detected in this RUNX1-CBFA2T3 sample. The RUNX1-CBFA2T3 fusion is a rare but recurrent gene rearrangement in AML (recently reviewed by Athanasiadou et al 15 ). In contrast to other RUNX1 fusion partners, CBFA2T3 (also known as MTG16) shows high sequence identity with RUNX1T1 leading to a RUNX1-CBFA2T3 chimeric protein that shares similar structural characteristics to RUNX1-RUNX1T1. 16 Further analyses are needed to determine if these observations translate into a similar susceptibility to chemotherapy and clinical outcome. Genes mutated in t(8;21) and inv(16) cohorts are shown in Figure 2A and detailed in supplemental Table 3. As previously reported, the most frequent mutations in both CBF subgroups were found in activated signaling genes (Figure 2A). Notably, 15 (31%) CBF samples contained 2 to 5 mutations in activated signaling genes. In the majority of cases, the sum of VAF did not exceed ;50% (Figure 2B), hence suggesting that they occurred in different subclones. The loss of 2 signaling mutations and expansion of a third one in a relapse specimen further support this concept (Figure 2C). Sixteen different genes were mutated in t(8;21) AML cohort: KIT (8/20, 40%); ASXL2 (5/20, 25%); FLT3, ASXL1 (4/20, 20% each); ZBTB7A, NRAS, TET2, and SMC1A (3/20, 15%); DNMT3A (2/20, 10%); and KDM6A, KMT2C, SMC3, STAG2, WT1, JAK2, and CSF3R (1/20, 5% each) (Figure 2A). No association was found between muta- tions and additional cytogenetic aberrations and clinical or laboratory characteristics (Figure 2A). As previously described, 17 an association was observed between t(8;21) and del(9q) or -Y. The activated signaling genes (14/20, 70%) were the most frequently mutated, followed by chromatin modifier (10/20, 50%), cohesin (5/20, 25% each), and DNA methylation (4/20, 20%) genes (Figure 2A). Mutations in the chromatin modifier ASXL2 were largely restricted to t(8;21) AML 2498 BLOOD, 19 MAY 2016 x VOLUME 127, NUMBER 20 Downloaded from http://ashpublications.org/blood/article-pdf/127/20/2498/1393512/2498.pdf by guest on 08 June 2022