Covering Hierarchical Dirichlet Mixture Models on binary data to enhance genomic stratifications in onco-hematology.

Onco-hematological studies are increasingly adopting statistical mixture models to support the advancement of the genomically-driven classification systems for blood cancer. Targeting enhanced patients stratification based on the sole role of molecular biology attracted much interest and contributes to bring personalized medicine closer to reality. In onco-hematology, Hierarchical Dirichlet Mixture Models (HDMM) have become one of the preferred method to cluster the genomics data, that include the presence or absence of gene mutations and cytogenetics anomalies, into components. This work unfolds the standard workflow used in onco-hematology to improve patient stratification and proposes alternative approaches to characterize the components and to assign patient to them, as they are crucial tasks usually supported by a priori clinical knowledge. We propose (a) to compute the parameters of the multinomial components of the HDMM or (b) to estimate the parameters of the HDMM components as if they were Multivariate Fisher's Non-Central Hypergeometric (MFNCH) distributions. Then, our approach to perform patients assignments to the HDMM components is designed to essentially determine for each patient its most likely component. We show on simulated data that the patients assignment using the MFNCH-based approach can be superior, if not comparable, to using the multinomial-based approach. Lastly, we illustrate on real Acute Myeloid Leukemia data how the utilization of MFNCH-based approach emerges as a good trade-off between the rigorous multinomial-based characterization of the HDMM components and the common refinement of them based on a priori clinical knowledge.

AML with complex karyotype: extreme genomic complexity revealed by combined long-read sequencing and Hi-C technology.

Acute myeloid leukemia with complex karyotype (CK-AML) is associated with poor prognosis, which is only in part explained by underlying TP53 mutations. Especially in the presence of complex chromosomal rearrangements, such as chromothripsis, the outcome of CK-AML is dismal. However, this degree of complexity of genomic rearrangements contributes to the leukemogenic phenotype and treatment resistance of CK-AML remains largely unknown. Applying an integrative workflow for the detection of structural variants (SVs) based on Oxford Nanopore (ONT) genomic DNA long-read sequencing (gDNA-LRS) and high-throughput chromosome confirmation capture (Hi-C) in a well-defined cohort of CK-AML identified regions with an extreme density of SVs. These rearrangements consisted to a large degree of focal amplifications enriched in the proximity of mammalian-wide interspersed repeat elements, which often result in oncogenic fusion transcripts, such as USP7::MVD, or the deregulation of oncogenic driver genes as confirmed by RNA-seq and ONT direct complementary DNA sequencing. We termed this novel phenomenon chromocataclysm. Thus, our integrative SV detection workflow combing gDNA-LRS and Hi-C enables to unravel complex genomic rearrangements at a very high resolution in regions hard to analyze by conventional sequencing technology, thereby providing an important tool to identify novel important drivers underlying cancer with complex karyotypic changes.

Clonal evolution of acute myeloid leukemia with FLT3-ITD mutation under treatment with midostaurin.

In the international randomized phase 3 RATIFY (Randomized AML Trial In FLT3 in patients less than 60 Years old) trial, the multikinase inhibitor midostaurin significantly improved overall and event-free survival in patients 18 to 59 years of age with FLT3-mutated acute myeloid leukemia (AML). However, only 59% of patients in the midostaurin arm achieved protocol-specified complete remission (CR), and almost half of patients achieving CR relapsed. To explore underlying mechanisms of resistance, we studied patterns of clonal evolution in patients with FLT3-internal tandem duplications (ITD)-positive AML who were entered in the RATIFY or German-Austrian Acute Myeloid Leukemia Study Group 16-10 trial and received treatment with midostaurin. To this end, paired samples from 54 patients obtained at time of diagnosis and at time of either relapsed or refractory disease were analyzed using conventional Genescan-based testing for FLT3-ITD and whole exome sequencing. At the time of disease resistance or progression, almost half of the patients (46%) became FLT3-ITD negative but acquired mutations in signaling pathways (eg, MAPK), thereby providing a new proliferative advantage. In cases with FLT3-ITD persistence, the selection of resistant ITD clones was found in 11% as potential drivers of disease. In 32% of cases, no FLT3-ITD mutational change was observed, suggesting either resistance mechanisms bypassing FLT3 inhibition or loss of midostaurin inhibitory activity because of inadequate drug levels. In summary, our study provides novel insights into the clonal evolution and resistance mechanisms of FLT3-ITD-mutated AML under treatment with midostaurin in combination with intensive chemotherapy.

Genomic heterogeneity in core-binding factor acute myeloid leukemia and its clinical implication.

Core-binding factor (CBF) acute myeloid leukemia (AML) encompasses AML with inv(16)(p13.1q22) and AML with t(8;21)(q22;q22.1). Despite sharing a common pathogenic mechanism involving rearrangements of the CBF transcriptional complex, there is growing evidence for considerable genotypic heterogeneity. We comprehensively characterized the mutational landscape of 350 adult CBF-AML [inv(16): n = 160, t(8;21): n = 190] performing targeted sequencing of 230 myeloid cancer-associated genes. Apart from common mutations in signaling genes, mainly NRAS, KIT, and FLT3, both CBF-AML entities demonstrated a remarkably diverse pattern with respect to the underlying cooperating molecular events, in particular in genes encoding for epigenetic modifiers and the cohesin complex. In addition, recurrent mutations in novel collaborating candidate genes such as SRCAP (5% overall) and DNM2 (6% of t(8;21) AML) were identified. Moreover, aberrations altering transcription and differentiation occurred at earlier leukemic stages and preceded mutations impairing proliferation. Lasso-penalized models revealed an inferior prognosis for t(8;21) AML, trisomy 8, as well as FLT3 and KIT exon 17 mutations, whereas NRAS and WT1 mutations conferred superior prognosis. Interestingly, clonal heterogeneity was associated with a favorable prognosis. When entering mutations by functional groups in the model, mutations in genes of the methylation group (ie, DNMT3A, TET2) had a strong negative prognostic impact.

Clonal evolution patterns in acute myeloid leukemia with NPM1 mutation.

Mutations in the nucleophosmin 1 (NPM1) gene are considered founder mutations in the pathogenesis of acute myeloid leukemia (AML). To characterize the genetic composition of NPM1 mutated (NPM1mut) AML, we assess mutation status of five recurrently mutated oncogenes in 129 paired NPM1mut samples obtained at diagnosis and relapse. We find a substantial shift in the genetic pattern from diagnosis to relapse including NPM1mut loss (n = 11). To better understand these NPM1mut loss cases, we perform whole exome sequencing (WES) and RNA-Seq. At the time of relapse, NPM1mut loss patients (pts) feature distinct mutational patterns that share almost no somatic mutation with the corresponding diagnosis sample and impact different signaling pathways. In contrast, profiles of pts with persistent NPM1mut are reflected by a high overlap of mutations between diagnosis and relapse. Our findings confirm that relapse often originates from persistent leukemic clones, though NPM1mut loss cases suggest a second "de novo" or treatment-associated AML (tAML) as alternative cause of relapse.