Acute promyelocytic leukemias share cooperative mutations with other myeloid-leukemia subgroups

Acute promyelocytic leukemias (APLs) are characterized by the expression of the PML-RARA oncogene, which is a product of the 15;17 chromosomal translocation. Two lines of evidence suggest that this genetic alteration is involved in the initiation of promyelocytic leukemogenesis in vivo : i) the t(15;17) often represents the only cytogenetic abnormality in an otherwise-normal karyotype; ii) PML-RARA is able to initiate leukemogenesis in mouse models on expression in hematopoietic stem cells (HSCs) or progenitors, giving rise to a disease that recapitulates clinical and morphological features of human APLs. 1 APL onset in mice, however, occurs after a long latency and with a low penetrance, suggesting that, like other cancers, APL is a multistep disease, whereby PML-RARA cooperates with secondary mutations for the full development of the leukemia phenotype. Notably, in the pre-leukemic phase of the disease, PML-RARA expression induces DNA damage in HSCs/progenitors and activates DNA repair, thus imparting a mutator phenotype that might contribute to disease progression. 2,3 Next-generation sequencing represents a powerful tool for the discovery of genetic alterations at high resolution. To identify gene mutations that might cooperate with PML-RARA in the leukemogenic process, we performed whole-exome sequencing of 5 leukemias that developed in PML-RARA transgenic mice (mouse APLs; mAPLs) and 11 patients’ leukemias expressing PML-RARA (human APLs; hAPLs) (Supplementary Table 1). For the scope of this investigation, we only considered non-synonymous single-nucleotide variants (SNVs) and small insertions/deletions (indels) occurring in the tumor DNA with an allelic

Acute promyelocytic leukemias (APLs) are characterized by the expression of the PML-RARA oncogene, which is a product of the 15;17 chromosomal translocation. Two lines of evidence suggest that this genetic alteration is involved in the initiation of promyelocytic leukemogenesis in vivo: i) the t(15;17) often represents the only cytogenetic abnormality in an otherwisenormal karyotype; ii) PML-RARA is able to initiate leukemogenesis in mouse models on expression in hematopoietic stem cells (HSCs) or progenitors, giving rise to a disease that recapitulates clinical and morphological features of human APLs. 1 APL onset in mice, however, occurs after a long latency and with a low penetrance, suggesting that, like other cancers, APL is a multistep disease, whereby PML-RARA cooperates with secondary mutations for the full development of the leukemia phenotype. Notably, in the pre-leukemic phase of the disease, PML-RARA expression induces DNA damage in HSCs/progenitors and activates DNA repair, thus imparting a mutator phenotype that might contribute to disease progression. 2,3 Next-generation sequencing represents a powerful tool for the discovery of genetic alterations at high resolution. To identify gene mutations that might cooperate with PML-RARA in the leukemogenic process, we performed whole-exome sequencing of 5 leukemias that developed in PML-RARA transgenic mice (mouse APLs; mAPLs) and 11 patients' leukemias expressing PML-RARA (human APLs; hAPLs) (Supplementary Table 1). For the scope of this investigation, we only considered non-synonymous single-nucleotide variants (SNVs) and small insertions/deletions (indels) occurring in the tumor DNA with an allelic frequency X25% (corresponding to a 450% frequency of cells carrying the mutation in the tumor sample, assuming the majority of mutations to be heterozygous). These experimental conditions allow identification of SNVs with a validation rate of 100%, as shown by an independent sequencing approach (Sanger sequencing; unpublished). We identified a total of 18 mutations in the 5 mAPLs (16 SNVs and 2 indels) and 73 mutations in the 11 hAPLs (59 SNVs and 14 indels) ( Table 1 and Supplementary Table 2). We validated all the 16 mAPL SNVs by Sanger sequencing (the 2 indels were not amplifiable for technical reasons), 18/18 randomly selected human SNVs and 1/1 indel (Supplementary Table 2).
We next combined the results of our sequencing analyses with those previously published for APLs (1 mAPL and 24 hAPLs), 4-7 obtaining a data set of 41 APL samples (6 mAPLs and 35 hAPLs). Analysis of this data set showed a total of 270 mutations affecting 248 human genes (231 SNVs, 34 indels and 5 splice junction mutations-SJMs), with a low frequency of mutations per case (B7.1 in hAPLs and B3.5 in mAPLs). In three cases (two mAPLs and one hAPL), we found no SNVs or indels with a frequency higher than 25%. Surprisingly, two hAPL samples showed several hundreds of low-frequency mutations, whose significance remains unclear (not shown).
Of the 248 mutated genes, 9 were found in X2 patients (recurrent mutations), and 8 showed a significantly higher mutation rate (qo0.005; FLT3; WT1; KRAS; CALR; CSMD1; DDR2; REV3L and TCERG1L; Table 2 and Supplementary Table 3). FLT3 and KRAS have been already described as cooperators of PML-RARA in mouse models of APLs, 8 whereas WT1 is infrequently mutated in hAPLs. 9 Of the remaining five, the DDR2 (discoidin domain receptor 2) tyrosine kinase is mutated in a small subset of squamous cell lung cancer. Notably, DDR2 mutations are critical oncogenic events for these tumors and confer high sensitivity to the multi-targeted kinase inhibitor dasatinib. 10 We found a total of 25 mutations affecting these eight genes (APL driver mutations), with FLT3 and WT1 being the most frequently involved ( Table 2).
AML mutations were separately analyzed in each of the AML subgroups. We found a total of 1360 mutated genes (B9 mutations per case), of which 153 were recurrently mutated (X2 patients in the same subgroup) and 40 showed a significantly higher mutation rate (qo0.005; Table 2 and  Supplementary Table 3). Notably, included in the 40 genes were 21 of the 22 significantly mutated genes identified in a recent analysis of 200 AMLs (also part of our data set 7 ), and additional mutated genes critical for AMLs (BCOR, ASXL1, GATA2, SUZ12 and DDX41) or for selected epithelial cancers (CTCF, PLCE1 and CHD4). The most frequently mutated genes were also significantly associated with specific AML subgroups: FLT3, IDH1, DNMT3A and PTPN11 with NPM1 þ AMLs; RUNX1, CEBPA and GATA2 with NK-NPM1 À AMLs; TP53 with AMLs with complex karyotypes; KIT with CBFB/MYH11 AMLs; IDH1 with AMLs with trisomy 8 (Table 2). Moreover, mutations in SF3B1, PTPN11, DNAH9, are present in both human and mouse leukemias.
We then analyzed the distribution of the significantly mutated APL and AML genes (n ¼ 44) across all samples (n ¼ 239). Twentyeight genes (B64%) were mutated in more than one cytogenetic subgroup, covering 383 of the 416 mutations identified in all samples (B92%). In the remaining 16 genes (B36%), mutations were instead associated with a specific subgroup, corresponding to just 33 of the identified mutations (B8%) and suggesting that subgroup specificity might be due to their low frequency. Indeed, we found a significant correlation between frequency of mutations per gene and numbers of subgroups where it is mutated (Spearman's coefficient value of 0.93). As regards APLs, 5 of the 8 significantly mutated genes were also found in other AML cytogenetic subgroups (FLT3; WT1; KRAS; CSMD1 and DDR2), covering a total of 100 mutations, while 3 (CALR; REV3L and TCERG1L) were only found in APLs, covering 6 mutations. Three genes that were found significantly mutated in other AML subgroups were also mutated in APLs (PHF6, FAM5C and PTPN11).