Comprehensive CRISPR-Cas9 screens identify genetic determinants of drug responsiveness in multiple myeloma.

The introduction of new drugs in the past years has substantially improved outcome in multiple myeloma (MM). However, the majority of patients eventually relapse and become resistant to one or multiple drugs. While the genetic landscape of relapsed/ resistant multiple myeloma has been elucidated, the causal relationship between relapse-specific gene mutations and the sensitivity to a given drug in MM has not systematically been evaluated. To determine the functional impact of gene mutations, we performed combined whole-exome sequencing (WES) of longitudinal patient samples with CRISPR-Cas9 drug resistance screens for lenalidomide, bortezomib, dexamethasone, and melphalan. WES of longitudinal samples from 16 MM patients identified a large number of mutations in each patient that were newly acquired or evolved from a small subclone (median 9, range 1-55), including recurrent mutations in TP53, DNAH5, and WSCD2. Focused CRISPR-Cas9 resistance screens against 170 relapse-specific mutations functionally linked 15 of them to drug resistance. These included cereblon E3 ligase complex members for lenalidomide, structural genes PCDHA5 and ANKMY2 for dexamethasone, RB1 and CDK2NC for bortezomib, and TP53 for melphalan. In contrast, inactivation of genes involved in the DNA damage repair pathway, including ATM, FANCA, RAD54B, and BRCC3, enhanced susceptibility to cytotoxic chemotherapy. Resistance patterns were highly drug specific with low overlap and highly correlated with the treatment-dependent clonal evolution in patients. The functional association of specific genetic alterations with drug sensitivity will help to personalize treatment of MM in the future.

Deregulated expression of circular RNAs in acute myeloid leukemia.

Circular RNAs (circRNAs) are dynamically regulated during differentiation and show cell type-specific expression, which is altered in cancer and can have a direct impact on its various hallmarks. We hypothesized that circRNA expression is deregulated in acute myeloid leukemia (AML) and that circRNA candidates might contribute to the pathogenesis of the disease. To identify leukemia-associated and differentiation-independent changes in circRNA expression, we determined the circular RNAome of 61 AML patients and 16 healthy hematopoietic stem and progenitor cell (HSPC) samples using ribosomal RNA-depleted RNA sequencing. We found hundreds of circRNAs that were differentially expressed between AML and healthy HSPCs. Gene set analysis found that many of these circRNAs were transcribed from genes implicated in leukemia biology. We discovered a circRNA derived from the T-cell transcription factor gene B cell CLL/lymphoma 11B, circBCL11B, which was exclusively expressed in AML patients, but not detected in healthy HSPCs, and associated with a T-cell-like gene expression signature. We were able to validate this finding in an independent cohort of 332 AML patients. Knockdown of circBCL11B had a negative effect on leukemic cell proliferation and resulted in increased cell death of leukemic cells, thereby suggesting circBCL11B as a novel functionally relevant candidate in AML pathogenesis. In summary, our study enables comprehensive insights into circRNA expression changes upon leukemic transformation and provides valuable information on the biology of leukemic cells and potential novel pathway dependencies that are relevant for AML therapy.

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 alterations in high-risk chronic lymphocytic leukemia frequently affect cell cycle key regulators and NOTCH1-regulated transcription.

To identify genomic alterations contributing to the pathogenesis of high-risk chronic lymphocytic leukemia (CLL) beyond the well-established role of TP53 aberrations, we comprehensively analyzed 75 relapsed/refractory and 71 treatment-naïve high-risk cases from prospective clinical trials by single nucleotide polymorphism arrays and targeted next-generation sequencing. Increased genomic complexity was a hallmark of relapsed/refractory and treatment-naïve high-risk CLL. In relapsed/refractory cases previously exposed to the selective pressure of chemo(immuno)therapy, gain(8)(q24.21) and del(9)(p21.3) were particularly enriched. Both alterations affect key regulators of cell-cycle progression, namely MYC and CDKN2A/B While homozygous CDKN2A/B loss has been directly associated with Richter transformation, we did not find this association for heterozygous loss of CDKN2A/B Gains in 8q24.21 were either focal gains in a MYC enhancer region or large gains affecting the MYC locus, but only the latter type was highly enriched in relapsed/refractory CLL (17%). In addition to a high frequency of NOTCH1 mutations (23%), we found recurrent genetic alterations in SPEN (4% mutated), RBPJ (8% deleted) and SNW1 (8% deleted), all affecting a protein complex that represses transcription of NOTCH1 target genes. We investigated the functional impact of these alterations on HES1, DTX1 and MYC gene transcription and found derepression of these NOTCH1 target genes particularly with SPEN mutations. In summary, we provide new insights into the genomic architecture of high-risk CLL, define novel recurrent DNA copy number alterations and refine knowledge on del(9p), gain(8q) and alterations affecting NOTCH1 signaling. This study was registered at ClinicalTrials.gov with number NCT01392079.

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.

Resting cells rely on the DNA helicase component MCM2 to build cilia.

Minichromosome maintenance (MCM) proteins facilitate replication by licensing origins and unwinding the DNA double strand. Interestingly, the number of MCM hexamers greatly exceeds the number of firing origins suggesting additional roles of MCMs. Here we show a hitherto unanticipated function of MCM2 in cilia formation in human cells and zebrafish that is uncoupled from replication. Zebrafish depleted of MCM2 develop ciliopathy-phenotypes including microcephaly and aberrant heart looping due to malformed cilia. In non-cycling human fibroblasts, loss of MCM2 promotes transcription of a subset of genes, which cause cilia shortening and centriole overduplication. Chromatin immunoprecipitation experiments show that MCM2 binds to transcription start sites of cilia inhibiting genes. We propose that such binding may block RNA polymerase II-mediated transcription. Depletion of a second MCM (MCM7), which functions in complex with MCM2 during its canonical functions, reveals an overlapping cilia-deficiency phenotype likely unconnected to replication, although MCM7 appears to regulate a distinct subset of genes and pathways. Our data suggests that MCM2 and 7 exert a role in ciliogenesis in post-mitotic tissues.

Circular RNAs of the nucleophosmin (NPM1) gene in acute myeloid leukemia.

In acute myeloid leukemia, there is growing evidence for splicing pattern deregulation, including differential expression of linear splice isoforms of the commonly mutated gene nucleophosmin (NPM1). In this study, we detect circular RNAs of NPM1 and quantify circRNA hsa_circ_0075001 in a cohort of NPM1 wild-type and mutated acute myeloid leukemia (n=46). Hsa_circ_0075001 expression correlates positively with total NPM1 expression, but is independent of the NPM1 mutational status. High versus low hsa_circ_0075001 expression defines patient subgroups characterized by distinct gene expression patterns, such as lower expression of components of the Toll-like receptor signaling pathway in high hsa_circ_0075001 expression cases. Global evaluation of circRNA expression in sorted healthy hematopoietic controls (n=10) and acute myeloid leukemia (n=10) reveals circRNA transcripts for 47.9% of all highly expressed genes. While circRNA expression correlates globally with parental gene expression, we identify hematopoietic differentiation-associated as well as acute myeloid leukemia subgroup-specific circRNA signatures.

Partitioned learning of deep Boltzmann machines for SNP data.

MOTIVATION: Learning the joint distributions of measurements, and in particular identification of an appropriate low-dimensional manifold, has been found to be a powerful ingredient of deep leaning approaches. Yet, such approaches have hardly been applied to single nucleotide polymorphism (SNP) data, probably due to the high number of features typically exceeding the number of studied individuals. RESULTS: After a brief overview of how deep Boltzmann machines (DBMs), a deep learning approach, can be adapted to SNP data in principle, we specifically present a way to alleviate the dimensionality problem by partitioned learning. We propose a sparse regression approach to coarsely screen the joint distribution of SNPs, followed by training several DBMs on SNP partitions that were identified by the screening. Aggregate features representing SNP patterns and the corresponding SNPs are extracted from the DBMs by a combination of statistical tests and sparse regression. In simulated case-control data, we show how this can uncover complex SNP patterns and augment results from univariate approaches, while maintaining type 1 error control. Time-to-event endpoints are considered in an application with acute myeloid leukemia patients, where SNP patterns are modeled after a pre-screening based on gene expression data. The proposed approach identified three SNPs that seem to jointly influence survival in a validation dataset. This indicates the added value of jointly investigating SNPs compared to standard univariate analyses and makes partitioned learning of DBMs an interesting complementary approach when analyzing SNP data. AVAILABILITY AND IMPLEMENTATION: A Julia package is provided at 'http://github.com/binderh/BoltzmannMachines.jl'. CONTACT: binderh@imbi.uni-freiburg.de. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

BiTrinA--multiscale binarization and trinarization with quality analysis.

MOTIVATION: When processing gene expression profiles or other biological data, it is often required to assign measurements to distinct categories (e.g. 'high' and 'low' and possibly 'intermediate'). Subsequent analyses strongly depend on the results of this quantization. Poor quantization will have potentially misleading effects on further investigations. We propose the BiTrinA package that integrates different multiscale algorithms for binarization and for trinarization of one-dimensional data with methods for quality assessment and visualization of the results. By identifying measurements that show large variations over different time points or conditions, this quality assessment can determine candidates that are related to the specific experimental setting. AVAILABILITY AND IMPLEMENTATION: BiTrinA is freely available on CRAN. CONTACT: hans.kestler@leibniz-fli.de or hans.kestler@uni-ulm.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.