While splicing changes caused by somatic mutations in SF3B1 are known, identifying full-length isoform changes may better elucidate the functional consequences of these mutations. We report nanopore sequencing of full-length cDNA from CLL samples with and without SF3B1 mutation, as well as normal B cell samples, giving a total of 149 million pass reads. We present FLAIR (Full-Length Alternative Isoform analysis of RNA), a computational workflow to identify high-confidence transcripts, perform differential splicing event analysis, and differential isoform analysis. Using nanopore reads, we demonstrate differential 3' splice site changes associated with SF3B1 mutation, agreeing with previous studies. We also observe a strong downregulation of intron retention events associated with SF3B1 mutation. Full-length transcript analysis links multiple alternative splicing events together and allows for better estimates of the abundance of productive versus unproductive isoforms. Our work demonstrates the potential utility of nanopore sequencing for cancer and splicing research.

Chronic lymphocytic leukaemia (CLL) is a malignancy of CD5⁺ B cells that is characterized by the accumulation of small, mature-appearing lymphocytes in the blood, marrow and lymphoid tissues. Signalling via surface immunoglobulin, which constitutes the major part of the B cell receptor, and several genetic alterations play a part in CLL pathogenesis, in addition to interactions between CLL cells and other cell types, such as stromal cells, T cells and nurse-like cells in the lymph nodes. The clinical progression of CLL is heterogeneous and ranges from patients who require treatment soon after diagnosis to others who do not require therapy for many years, if at all. Several factors, including the immunoglobulin heavy-chain variable region gene (IGHV) mutational status, genomic changes, patient age and the presence of comorbidities, should be considered when defining the optimal management strategies, which include chemotherapy, chemoimmunotherapy and/or drugs targeting B cell receptor signalling or inhibitors of apoptosis, such as BCL-2. Research on the biology of CLL has profoundly enhanced our ability to identify patients who are at higher risk for disease progression and our capacity to treat patients with drugs that selectively target distinctive phenotypic or physiological features of CLL. How these and other advances have shaped our current understanding and treatment of patients with CLL is the subject of this Primer.

Combining single-cell cytometry datasets increases the analytical flexibility and the statistical power of data analyses. However, in many cases the full potential of co-analyses is not reached due to technical variance between data from different experimental batches. Here, we present cyCombine, a method to robustly integrate cytometry data from different batches, experiments, or even different experimental techniques, such as CITE-seq, flow cytometry, and mass cytometry. We demonstrate that cyCombine maintains the biological variance and the structure of the data, while minimizing the technical variance between datasets. cyCombine does not require technical replicates across datasets, and computation time scales linearly with the number of cells, allowing for integration of massive datasets. Robust, accurate, and scalable integration of cytometry data enables integration of multiple datasets for primary data analyses and the validation of results using public datasets.

Lineage-tracing methods have enabled characterization of clonal dynamics in complex populations, but generally lack the ability to integrate genomic, epigenomic and transcriptomic measurements with live-cell manipulation of specific clones of interest. We developed a functionalized lineage-tracing system, ClonMapper, which integrates DNA barcoding with single-cell RNA sequencing and clonal isolation to comprehensively characterize thousands of clones within heterogeneous populations. Using ClonMapper, we identified subpopulations of a chronic lymphocytic leukemia cell line with distinct clonal compositions, transcriptional signatures and chemotherapy survivorship trajectories; patterns that were also observed in primary human chronic lymphocytic leukemia. The ability to retrieve specific clones before, during and after treatment enabled direct measurements of clonal diversification and durable subpopulation transcriptional signatures. ClonMapper is a powerful multifunctional approach to dissect the complex clonal dynamics of tumor progression and therapeutic response.

How the genomic features of a patient's cancer relate to individual disease kinetics remains poorly understood. Here we used the indolent growth dynamics of chronic lymphocytic leukaemia (CLL) to analyse the growth rates and corresponding genomic patterns of leukaemia cells from 107 patients with CLL, spanning decades-long disease courses. We found that CLL commonly demonstrates not only exponential expansion but also logistic growth, which is sigmoidal and reaches a certain steady-state level. Each growth pattern was associated with marked differences in genetic composition, the pace of disease progression and the extent of clonal evolution. In a subset of patients, whose serial samples underwent next-generation sequencing, we found that dynamic changes in the disease course of CLL were shaped by the genetic events that were already present in the early slow-growing stages. Finally, by analysing the growth rates of subclones compared with their parental clones, we quantified the growth advantage conferred by putative CLL drivers in vivo.

Abstract Previous studies have revealed a critical role of methylation deregulation in the onset and progression of chronic lymphocytic leukemia (CLL). In mammalian cells, DNA methylation is dynamically established by the DNA methyltransferase 3 (DNMT3) family of de novo methyltransferases DNMT3A. Although mutations of DNMT3A are rarely observed in CLL, our RNA-sequencing (RNA-seq) analysis of 107 human CLLs show that low DNMT3A expression is associated with more aggressive disease, and supports a driving role of DNMT3A loss in CLL. To test this hypothesis, we generated a conditional knock-out mouse model with B cell-restricted deletion of Dnmt3a. Homozygous Dnmt3a depletion in B cells resulted in the development of CD5+ B cell leukemia mimicking human CLL with 100% penetrance at a median age of onset of 5.3 months, and heterozygous Dnmt3a depletion yielded a disease penetrance of 89% with a median onset at 18.5 months, confirming its role as a haplo-insufficient tumor suppressor. Given the known role of Dnmt3a as a de novo methyltransferase, we first evaluated the impact of Dnmt3a depletion on global DNA methylation in non-leukemic CD5+ B cells isolated from the peritoneal cavity by cell sorting (i.e. B1a cells) using reduced representation bisulfite sequencing (RRBS). We identified a set of differentially methylated regions (DMRs) (difference>0.2), mostly hypomethylated, in Dnmt3afl/fl versus WT B1a cells (473 hypomethylated, 19 hypermethylated). Genes with dysregulated methylation were highly enriched in pathways involved in immune response (e.g., Interferon-α signaling, JAK/STAT3 signaling) and proliferation (Wnt Signaling and Notch signaling). Given the prominent hypomethylation changes observed in Dnmt3a depleted B1a cells, we investigated whether these would lead to altered gene transcript expression. Using RNA-seq, we detected 460 downregulated and 168 upregulated genes in the Dnmt3afl/fl B1a cells compared to WT B1a cells (FDR<0.05, fold change >2). Consistent with the methylation data, differentially expressed genes were likewise enriched for JAK/STAT3 signaling, Wnt Signaling and Notch signaling, supporting a direct influence of dysregulated methylation on downstream signaling cascades. We investigated the changes in methylomes of the CLL cells arising from the Dnmt3afl/fl animals. Compared to WT B1a cells, Dnmt3afl/fl CLL cells generated 1335 hypomethylated and 2369 hypermethylated DMRs in. Focusing on genes that were hypomethylated in CLL compared to WT B1a cells, we found that these were highly enriched for several oncogenic signaling pathways including Notch signaling and Wnt Signaling, consistent with the pre-leukemia findings. RNA-seq analysis identified more upregulated (n=2801) than downregulated (n=1244) genes in CLL cells compared to WT B1a cells (FDR<0.05, FC>2), supporting a role of Dnmt3a depletion in transcriptional activation. We observed a general upregulation of Notch signaling genes and the downstream Notch targets, implicating Notch activation in this CLL mouse model. Of note, we showed Dnmt3a-depleted CLL cells to be highly sensitive to Notch inhibitor DAPT both in vitro and in a transplantable mouse model. Consistently, primary human CLL cells with low DNMT3A expression were more sensitive to DAPT than those with higher DNMT3A expression (P=0.005, Spearman correlation), despite similar sensitivity to ibrutinib and venetoclax. Together, our results have confirmed the causal role of Dnmt3a downregulation in CLL generation. We provide evidence in support of the interaction between Dnmt3a-dependent methylation changes and hyperactivation of Notch signaling in transcriptional reprogramming and transformation of B1a cells into CLL. Furthermore, we demonstrate differential sensitivity of DNMT3A high and low expressing primary CLLs to Notch inhibition, indicative of consistent dependencies of human and murine CLLs. Thus, the Dnmt3a models provides a unique opportunity for the study of non-mutational Notch activation, and a useful platform for the study of Notch-signaling targeted therapeutics. Disclosures Kipps: Abbott Laboratories: Consultancy, Research Funding; Celgene Corporation: Consultancy, Honoraria, Research Funding; Pharmacyclics LLC, an Abbvie Company: Consultancy, Honoraria, Other: Travel, Accommodations, Expenses, Research Funding, Speakers Bureau; Genentech, Inc.: Honoraria, Research Funding, Speakers Bureau; Gilead Sciences, Inc.: Honoraria, Research Funding; GlaxoSmithKline: Research Funding; MedImmune Inc: Research Funding; Moores Cancer Center: Current Employment; Oncternal Therapeutics, Inc.: Current holder of stock options in a privately-held company, Other: Stock or other ownership, Patents & Royalties: Cirmtuzumab was developed by Thomas J. Kipps in the Thomas J. Kipps laboratory and licensed by the University of California to Oncternal Therapeutics, Inc., which provided stock options and research funding to the Thomas J. Kipps laboratory., Research Funding; AbbVie: Consultancy, Honoraria, Other, Speakers Bureau; DAVAOncology: Consultancy, Honoraria, Other; DAVA Pharmaceuticals: Speakers Bureau; Bionest Partner: Other; Celgene: Consultancy, Honoraria, Other, Research Funding; Genetech: Honoraria, Other; Genentech-Roche: Consultancy; Gilead Sciences: Consultancy, Honoraria, Other, Speakers Bureau; Janssen: Consultancy, Honoraria, Other, Research Funding, Speakers Bureau; Roche: Honoraria, Other; MD Anderson Cancer Center: Research Funding; Velos: Research Funding; CRIM: Research Funding; Indy Hematology Review: Other; TG Therapeutics: Other; Verstem: Other, Speakers Bureau; University of California, San Diego: Current Employment; Pharmacyclics/AbbVie: Honoraria, Research Funding; Breast Cancer Research Foundation: Research Funding; SCOR - The Leukemia and Lymphoma Society: Research Funding; National Cancer Institute/NIH: Honoraria, Research Funding; Genentech/Roche: Honoraria; European Research Initiative on CLL (ERIC): Honoraria. Neuberg: Madrigal Pharmaceuticals: Other: Stock ownership; Pharmacyclics: Research Funding. Letai: Flash Therapeutics: Other: equity holding member of the scientific advisory board; Dialectic Therapeutics: Other: equity holding member of the scientific advisory board; Zentalis Pharmaceuticals: Other: equity holding member of the scientific advisory board. Wu: BioNTech: Current equity holder in publicly-traded company; Pharmacyclics: Research Funding.

While cancers evolve during disease progression and in response to therapy, temporal dynamics remain difficult to study in humans due to the lack of consistent barcodes marking individual clones in vivo. We employ mitochondrial single-cell assay for transposase-accessible chromatin with sequencing to profile 163,279 cells from 9 patients with chronic lymphocytic leukemia (CLL) collected across disease course and utilize mitochondrial DNA (mtDNA) mutations as natural genetic markers of cancer clones. We observe stable propagation of mtDNA mutations over years in the absence of strong selective pressure, indicating clonal persistence, but dramatic changes following tight bottlenecks, including disease transformation and relapse posttherapy, paralleled by acquisition of copy-number variants and changes in chromatin accessibility and gene expression. Furthermore, we link CLL subclones to distinct chromatin states, providing insight into nongenetic sources of relapse. mtDNA mutations thus mirror disease history and provide naturally occurring genetic barcodes to enable patient-specific study of cancer subclonal dynamics.

Significance

Single-cell multi-omic profiling of CLL reveals the utility of somatic mtDNA mutations as in vivo barcodes, which mark subclones that can evolve over time along with changes in accessible chromatin and gene expression profiles to capture dynamics of disease evolution. See related commentary by Hilton and Scott, p. 2965. This article is highlighted in the In This Issue feature, p. 2945.

We conducted a multicenter pilot investigation of the safety and feasibility of bone marrow transplantation (BMT) in adults with severe sickle cell disease (SCD) (NCT 01565616) using a reduced toxicity preparative regimen of busulfan (13.2 mg/kg), fludarabine (175 mg/m² ) and thymoglobulin (6 mg/kg) and cyclosporine or tacrolimus and methotrexate for graft-vs-host disease (GVHD) prophylaxis. Twenty-two patients (median age 22 years; range 17-36) were enrolled at eight centers. Seventeen patients received marrow from an HLA-identical sibling donor and five patients received marrow from an 8/8 HLA-allele matched unrelated donor. Before BMT, patients had stroke, acute chest syndrome, recurrent pain events, were receiving regular red blood cell transfusions, or had an elevated tricuspid regurgitant jet (TRJ) velocity, which fulfilled eligibility criteria. Four patients developed grades II-III acute GVHD (18%) and six developed chronic GVHD (27%) that was moderate in two and severe in one patient. One patient died of intracranial hemorrhage and one of GVHD. Nineteen patients had stable donor chimerism, 1-year post-transplant. One patient who developed secondary graft failure survives disease-free after a second BMT. The one-year overall survival and event-free survival (EFS) are 91% (95% CI 68%-98%) and 86% (95% CI, 63%-95%), respectively, and 3-year EFS is 82%. Statistically significant improvements in the pain interference and physical function domains of health-related quality of life were observed. The study satisfied the primary endpoint of 1-year EFS ≥70%. This regimen is being studied in a prospective clinical trial comparing HLA-matched donor BMT with standard of care in adults with severe SCD (NCT02766465).

SF3B1 is recurrently mutated in chronic lymphocytic leukemia (CLL), but its role in the pathogenesis of CLL remains elusive. Here, we show that conditional expression of Sf3b1-K700E mutation in mouse B cells disrupts pre-mRNA splicing, alters cell development, and induces a state of cellular senescence. Combination with Atm deletion leads to the overcoming of cellular senescence and the development of CLL-like disease in elderly mice. These CLL-like cells show genome instability and dysregulation of multiple CLL-associated cellular processes, including deregulated B cell receptor signaling, which we also identified in human CLL cases. Notably, human CLLs harboring SF3B1 mutations exhibit altered response to BTK inhibition. Our murine model of CLL thus provides insights into human CLL disease mechanisms and treatment.

Most human cancers converge to a deregulated methylome with reduced global levels and elevated methylation at select CpG islands. To investigate the emergence and dynamics of the cancer methylome, we characterized genome-wide DNA methylation in pre-neoplastic monoclonal B cell lymphocytosis (MBL) and chronic lymphocytic leukemia (CLL), including serial samples collected across disease course. We detected the aberrant tumor-associated methylation landscape at CLL diagnosis and found no significantly differentially methylated regions in the high-count MBL-to-CLL transition. Patient methylomes showed remarkable stability with natural disease and post-therapy progression. Single CLL cells were consistently aberrantly methylated, indicating a homogeneous transition to the altered epigenetic state, and a distinct expression profile together with MBL cells compared to normal B cells. Our longitudinal analysis reveals the cancer methylome to emerge early, which may provide a platform for subsequent genetically-driven growth dynamics and together with its persistent presence suggests a central role in the normal-to-cancer transition.