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.

Chronic lymphocytic leukemia (CLL) is characterized by disordered DNA methylation, suggesting these epigenetic changes might play a critical role in disease onset and progression. The methyltransferase DNMT3A is a key regulator of DNA methylation. Although DNMT3A somatic mutations in CLL are rare, we found that low DNMT3A expression is associated with more aggressive disease. A conditional knockout mouse model showed that homozygous depletion of Dnmt3a from B cells results in the development of CLL with 100% penetrance at a median age of onset of 5.3 months, and heterozygous Dnmt3a depletion yields a disease penetrance of 89% with a median onset at 18.5 months, confirming its role as a haploinsufficient tumor suppressor. B1a cells were confirmed as the cell of origin of disease in this model, and Dnmt3a depletion resulted in focal hypomethylation and activation of Notch and Myc signaling. Amplification of chromosome 15 containing the Myc gene was detected in all CLL mice tested, and infiltration of high-Myc-expressing CLL cells in the spleen was observed. Notably, hyperactivation of Notch and Myc signaling was exclusively observed in the Dnmt3a CLL mice, but not in three other CLL mouse models tested (Sf3b1-Atm, Ikzf3, and MDR), and Dnmt3a-depleted CLL were sensitive to pharmacologic inhibition of Notch signaling in vitro and in vivo. Consistent with these findings, human CLL samples with lower DNMT3A expression were more sensitive to Notch inhibition than those with higher DNMT3A expression. Altogether, these results suggest that Dnmt3a depletion induces CLL that is highly dependent on activation of Notch and Myc signaling. SIGNIFICANCE: Loss of DNMT3A expression is a driving event in CLL and is associated with aggressive disease, activation of Notch and Myc signaling, and enhanced sensitivity to Notch inhibition.

Treatment of chronic lymphocytic leukemia (CLL) has shifted from chemo-immunotherapy to targeted agents. To define the evolutionary dynamics induced by targeted therapy in CLL, we perform serial exome and transcriptome sequencing for 61 ibrutinib-treated CLLs. Here, we report clonal shifts (change >0.1 in clonal cancer cell fraction, Q < 0.1) in 31% of patients during the first year of therapy, associated with adverse outcome. We also observe transcriptional downregulation of pathways mediating energy metabolism, cell cycle, and B cell receptor signaling. Known and previously undescribed mutations in BTK and PLCG2, or uncommonly, other candidate alterations are present in seventeen subjects at the time of progression. Thus, the frequently observed clonal shifts during the early treatment period and its potential association with adverse outcome may reflect greater evolutionary capacity, heralding the emergence of drug-resistant clones.

In chronic lymphocytic leukemia (CLL), B-cell receptor signaling, tumor-microenvironment interactions, and somatic mutations drive disease progression. To better understand the intersection between the microenvironment and molecular events in CLL pathogenesis, we integrated bulk transcriptome profiling of paired peripheral blood (PB) and lymph node (LN) samples from 34 patients. Oncogenic processes were upregulated in LN compared with PB and in immunoglobulin heavy-chain variable (IGHV) region unmutated compared with mutated cases. Single-cell RNA sequencing (scRNA-seq) distinguished 3 major cell states: quiescent, activated, and proliferating. The activated subpopulation comprised only 2.2% to 4.3% of the total tumor bulk in LN samples. RNA velocity analysis found that CLL cell fate in LN is unidirectional, starts in the proliferating state, transitions to the activated state, and ends in the quiescent state. A 10-gene signature derived from activated tumor cells was associated with inferior treatment-free survival (TFS) and positively correlated with the proportion of activated CD4+ memory T cells and M2 macrophages in LN. Whole exome sequencing (WES) of paired PB and LN samples showed subclonal expansion in LN in approximately half of the patients. Since mouse models have implicated activation-induced cytidine deaminase in mutagenesis, we compared AICDA expression between cases with and without clonal evolution but did not find a difference. In contrast, the presence of a T-cell inflamed microenvironment in LN was associated with clonal stability. In summary, a distinct minor tumor subpopulation underlies CLL pathogenesis and drives the clinical outcome. Clonal trajectories are shaped by the LN milieu, where T-cell immunity may contribute to suppressing clonal outgrowth. The clinical study is registered at clinicaltrials.gov as NCT00923507.

Intratumoral heterogeneity plays a critical role in tumor evolution. To define the contribution of DNA methylation to heterogeneity within tumors, we performed genome-scale bisulfite sequencing of 104 primary chronic lymphocytic leukemias (CLLs). Compared with 26 normal B cell samples, CLLs consistently displayed higher intrasample variability of DNA methylation patterns across the genome, which appears to arise from stochastically disordered methylation in malignant cells. Transcriptome analysis of bulk and single CLL cells revealed that methylation disorder was linked to low-level expression. Disordered methylation was further associated with adverse clinical outcome. We therefore propose that disordered methylation plays a similar role to that of genetic instability, enhancing the ability of cancer cells to search for superior evolutionary trajectories.

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.

Hotspot mutation of IKZF3 (IKZF3-L162R) has been identified as a putative driver of chronic lymphocytic leukemia (CLL), but its function remains unknown. Here, we demonstrate its driving role in CLL through a B cell-restricted conditional knockin mouse model. Mutant Ikzf3 alters DNA binding specificity and target selection, leading to hyperactivation of B cell receptor (BCR) signaling, overexpression of nuclear factor κB (NF-κB) target genes, and development of CLL-like disease in elderly mice with a penetrance of ~40%. Human CLL carrying either IKZF3 mutation or high IKZF3 expression was associated with overexpression of BCR/NF-κB pathway members and reduced sensitivity to BCR signaling inhibition by ibrutinib. Our results thus highlight IKZF3 oncogenic function in CLL via transcriptional dysregulation and demonstrate that this pro-survival function can be achieved by either somatic mutation or overexpression of this CLL driver. This emphasizes the need for combinatorial approaches to overcome IKZF3-mediated BCR inhibitor resistance.

Unlike many other hematologic malignancies, Richter syndrome (RS), an aggressive B cell lymphoma originating from indolent chronic lymphocytic leukemia, is responsive to PD-1 blockade. To discover the determinants of response, we analyze single-cell transcriptome data generated from 17 bone marrow samples longitudinally collected from 6 patients with RS. Response is associated with intermediate exhausted CD8 effector/effector memory T cells marked by high expression of the transcription factor ZNF683, determined to be evolving from stem-like memory cells and divergent from terminally exhausted cells. This signature overlaps with that of tumor-infiltrating populations from anti-PD-1 responsive solid tumors. ZNF683 is found to directly target key T cell genes (TCF7, LMO2, CD69) and impact pathways of T cell cytotoxicity and activation. Analysis of pre-treatment peripheral blood from 10 independent patients with RS treated with anti-PD-1, as well as patients with solid tumors treated with anti-PD-1, supports an association of ZNF683high T cells with response.

The B cell leukemia/lymphoma 2 (BCL-2) inhibitor venetoclax is effective in chronic lymphocytic leukemia (CLL); however, resistance may develop over time. Other lymphoid malignancies such as diffuse large B cell lymphoma (DLBCL) are frequently intrinsically resistant to venetoclax. Although genomic resistance mechanisms such as BCL2 mutations have been described, this probably only explains a subset of resistant cases. Using 2 complementary functional precision medicine techniques - BH3 profiling and high-throughput kinase activity mapping - we found that hyperphosphorylation of BCL-2 family proteins, including antiapoptotic myeloid leukemia 1 (MCL-1) and BCL-2 and proapoptotic BCL-2 agonist of cell death (BAD) and BCL-2 associated X, apoptosis regulator (BAX), underlies functional mechanisms of both intrinsic and acquired resistance to venetoclax in CLL and DLBCL. Additionally, we provide evidence that antiapoptotic BCL-2 family protein phosphorylation altered the apoptotic protein interactome, thereby changing the profile of functional dependence on these prosurvival proteins. Targeting BCL-2 family protein phosphorylation with phosphatase-activating drugs rewired these dependencies, thus restoring sensitivity to venetoclax in a panel of venetoclax-resistant lymphoid cell lines, a resistant mouse model, and in paired patient samples before venetoclax treatment and at the time of progression.

B cell malignancies frequently colonize the bone marrow. The mechanisms responsible for this preferential homing are incompletely understood. Here we studied multiple myeloma (MM) as a model of a terminally differentiated B cell malignancy that selectively colonizes the bone marrow. We found that extracellular CyPA (eCyPA), secreted by bone marrow endothelial cells (BMECs), promoted the colonization and proliferation of MM cells in an in vivo scaffold system via binding to its receptor, CD147, on MM cells. The expression and secretion of eCyPA by BMECs was enhanced by BCL9, a Wnt-β-catenin transcriptional coactivator that is selectively expressed by these cells. eCyPA levels were higher in bone marrow serum than in peripheral blood in individuals with MM, and eCyPA-CD147 blockade suppressed MM colonization and tumor growth in the in vivo scaffold system. eCyPA also promoted the migration of chronic lymphocytic leukemia and lymphoplasmacytic lymphoma cells, two other B cell malignancies that colonize the bone marrow and express CD147. These findings suggest that eCyPA-CD147 signaling promotes the bone marrow homing of B cell malignancies and offer a compelling rationale for exploring this axis as a therapeutic target for these malignancies.

Mitochondrial apoptosis can be effectively targeted in lymphoid malignancies with the FDA-approved B cell lymphoma 2 (BCL-2) inhibitor venetoclax, but resistance to this agent is emerging. We show that venetoclax resistance in chronic lymphocytic leukemia is associated with complex clonal shifts. To identify determinants of resistance, we conducted parallel genome-scale screens of the BCL-2-driven OCI-Ly1 lymphoma cell line after venetoclax exposure along with integrated expression profiling and functional characterization of drug-resistant and engineered cell lines. We identified regulators of lymphoid transcription and cellular energy metabolism as drivers of venetoclax resistance in addition to the known involvement by BCL-2 family members, which were confirmed in patient samples. Our data support the implementation of combinatorial therapy with metabolic modulators to address venetoclax resistance.

Clinical vignette: A 49-year-old man with stage IV BRAFV600E-driven melanoma was initiated on twice-daily 960 mg of vemurafenib for treatment of progressive and recurrent subcutaneous metastatic disease of the left lower extremity. The patient's melanoma responded well to targeted BRAF inhibition. At treatment onset, hematologic parameters were all within normal limits; however, within three months of initiating therapy, wbc were found to be elevated (to 20 K) with sustained lymphocytosis of mature phenotype. Immunophenotypic analysis was consistent with chronic lymphocytic leukemia (CLL), and FISH results revealed presence of the CLL-associated deletion in chromosome 13q14 as well as in 2p33. Vemurafenib was withdrawn after approximately one year of therapy, and subsequently, his peripheral lymphocytosis resolved and CLL regressed. Nevertheless, a monoclonal B cell population persisted even 732 days after discontinuation of vemurafenib.

Transformation to aggressive disease histologies generates formidable clinical challenges across cancers, but biological insights remain few. We modeled the genetic heterogeneity of chronic lymphocytic leukemia (CLL) through multiplexed in vivo CRISPR-Cas9 B-cell editing of recurrent CLL loss-of-function drivers in mice and recapitulated the process of transformation from indolent CLL into large cell lymphoma [i.e., Richter syndrome (RS)]. Evolutionary trajectories of 64 mice carrying diverse combinatorial gene assortments revealed coselection of mutations in Trp53, Mga, and Chd2 and the dual impact of clonal Mga/Chd2 mutations on E2F/MYC and interferon signaling dysregulation. Comparative human and murine RS analyses demonstrated tonic PI3K signaling as a key feature of transformed disease, with constitutive activation of the AKT and S6 kinases, downmodulation of the PTEN phosphatase, and convergent activation of MYC/PI3K transcriptional programs underlying enhanced sensitivity to MYC/mTOR/PI3K inhibition. This robust experimental system presents a unique framework to study lymphoid biology and therapy.

RNA splicing factor SF3B1 is recurrently mutated in various cancers, particularly in hematologic malignancies. We previously reported that coexpression of Sf3b1 mutation and Atm deletion in B cells, but not either lesion alone, leads to the onset of chronic lymphocytic leukemia (CLL) with CLL cells harboring chromosome amplification. However, the exact role of Sf3b1 mutation and Atm deletion in chromosomal instability (CIN) remains unclear. Here, we demonstrated that SF3B1 mutation promotes centromeric R-loop (cen-R-loop) accumulation, leading to increased chromosome oscillation, impaired chromosome segregation, altered spindle architecture, and aneuploidy, which could be alleviated by removal of cen-R-loop and exaggerated by deletion of ATM. Aberrant splicing of key genes involved in R-loop processing underlay augmentation of cen-R-loop, as overexpression of the normal isoform, but not the altered form, mitigated mitotic stress in SF3B1-mutant cells. Our study identifies a critical role of splice variants in linking RNA splicing dysregulation and CIN and highlights cen-R-loop augmentation as a key mechanism for leukemogenesis.

Major international projects are underway that are aimed at creating a comprehensive catalogue of all the genes responsible for the initiation and progression of cancer. These studies involve the sequencing of matched tumour-normal samples followed by mathematical analysis to identify those genes in which mutations occur more frequently than expected by random chance. Here we describe a fundamental problem with cancer genome studies: as the sample size increases, the list of putatively significant genes produced by current analytical methods burgeons into the hundreds. The list includes many implausible genes (such as those encoding olfactory receptors and the muscle protein titin), suggesting extensive false-positive findings that overshadow true driver events. We show that this problem stems largely from mutational heterogeneity and provide a novel analytical methodology, MutSigCV, for resolving the problem. We apply MutSigCV to exome sequences from 3,083 tumour-normal pairs and discover extraordinary variation in mutation frequency and spectrum within cancer types, which sheds light on mutational processes and disease aetiology, and in mutation frequency across the genome, which is strongly correlated with DNA replication timing and also with transcriptional activity. By incorporating mutational heterogeneity into the analyses, MutSigCV is able to eliminate most of the apparent artefactual findings and enable the identification of genes truly associated with cancer.