Clonal heterogeneity of FLT3-ITD detected by high-throughput amplicon sequencing correlates with adverse prognosis in acute myeloid leukemia

In acute myeloid leukemia (AML), internal tandem duplications (ITDs) of FLT3 are frequent mutations associated with unfavorable prognosis. At diagnosis, the FLT3-ITD status is routinely assessed by fragment analysis, providing information about the length but not the position and sequence of the ITD. To overcome this limitation, we performed cDNA-based high-throughput amplicon sequencing (HTAS) in 250 FLT3-ITD positive AML patients, treated on German AML Cooperative Group (AMLCG) trials. FLT3-ITD status determined by routine diagnostics was confirmed by HTAS in 242 out of 250 patients (97%). The total number of ITDs detected by HTAS was higher than in routine diagnostics (n = 312 vs. n = 274). In particular, HTAS detected a higher number of ITDs per patient compared to fragment analysis, indicating higher sensitivity for subclonal ITDs. Patients with more than one ITD according to HTAS had a significantly shorter overall and relapse free survival. There was a close correlation between FLT3-ITD mRNA levels in fragment analysis and variant allele frequency in HTAS. However, the abundance of long ITDs (≥75nt) was underestimated by HTAS, as the size of the ITD affected the mappability of the corresponding sequence reads. In summary, this study demonstrates that HTAS is a feasible approach for FLT3-ITD detection in AML patients, delivering length, position, sequence and mutational burden of this alteration in a single assay with high sensitivity. Our findings provide insights into the clonal architecture of FLT3-ITD positive AML and have clinical implications

Mediation analysis reveals common mechanisms of RUNX1 point mutations and RUNX1/RUNX1T1 fusions influencing survival of patients with acute myeloid leukemia

Alterations of RUNX1 in acute myeloid leukemia (AML) are associated with either a more favorable outcome in the case of the RUNX1/RUNX1T1 fusion or unfavorable prognosis in the case of point mutations. In this project we aimed to identify genes responsible for the observed differences in outcome that are common to both RUNX1 alterations. Analyzing four AML gene expression data sets (n = 1514), a total of 80 patients with RUNX1/RUNX1T1 and 156 patients with point mutations in RUNX1 were compared. Using the statistical tool of mediation analysis we identified the genes CD109, HOPX, and KIAA0125 as candidates for mediator genes. In an analysis of an independent validation cohort, KIAA0125 again showed a significant influence with respect to the impact of the RUNX1/RUNX1T1 fusion. While there were no significant results for the other two genes in this smaller validation cohort, the observed relations linked with mediation effects (i.e., those between alterations, gene expression and survival) were almost without exception as strong as in the main analysis. Our analysis demonstrates that mediation analysis is a powerful tool in the identification of regulative networks in AML subgroups and could be further used to characterize the influence of genetic alterations

A 29-gene and cytogenetic score for the prediction of resistance to induction treatment in acute myeloid leukemia

Primary therapy resistance is a major problem in acute myeloid leukemia treatment. We set out to develop a powerful and robust predictor for therapy resistance for intensively treated adult patients. We used two large gene expression data sets (n=856) to develop a predictor of therapy resistance, which was validated in an independent cohort analyzed by RNA sequencing (n=250). In addition to gene expression markers, standard clinical and laboratory variables as well as the mutation status of 68 genes were considered during construction of the model. The final predictor (PS29MRC) consisted of 29 gene expression markers and a cytogenetic risk classification. A continuous predictor is calculated as a weighted linear sum of the individual variables. In addition, a cut off was defined to divide patients into a high-risk and a low-risk group for resistant disease. PS29MRC was highly significant in the validation set, both as a continuous score (OR=2.39, P=8.63·10−9, AUC=0.76) and as a dichotomous classifier (OR=8.03, P=4.29·10−9); accuracy was 77%. In multivariable models, only TP53 mutation, age and PS29MRC (continuous: OR=1.75, P=0.0011; dichotomous: OR=4.44, P=0.00021) were left as significant variables. PS29MRC dominated all models when compared with currently used predictors, and also predicted overall survival independently of established markers. When integrated into the European LeukemiaNet (ELN) 2017 genetic risk stratification, four groups (median survival of 8, 18, 41 months, and not reached) could be defined (P=4.01·10−10). PS29MRC will make it possible to design trials which stratify induction treatment according to the probability of response, and refines the ELN 2017 classification

ZBTB7A mutations in acute myeloid leukaemia with t(8;21) translocation

The t(8;21) translocation is one of the most frequent cytogenetic abnormalities in acute myeloid leukaemia (AML) and results in the RUNX1/RUNX1T1 rearrangement. Despite the causative role of the RUNX1/RUNX1T1 fusion gene in leukaemia initiation, additional genetic lesions are required for disease development. Here we identify recurring ZBTB7A mutations in 23% (13/56) of AML t(8;21) patients, including missense and truncating mutations resulting in alteration or loss of the C-terminal zinc-finger domain of ZBTB7A. The transcription factor ZBTB7A is important for haematopoietic lineage fate decisions and for regulation of glycolysis. On a functional level, we show that ZBTB7A mutations disrupt the transcriptional repressor potential and the anti-proliferative effect of ZBTB7A. The specific association of ZBTB7A mutations with t(8;21) rearranged AML points towards leukaemogenic cooperativity between mutant ZBTB7A and the RUNX1/RUNX1T1 fusion