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.

Monoclonal B-cell lymphocytosis (MBL) is a precursor to CLL. Other than age, sex, and CLL family-history, little is known about factors associated with MBL risk. A polygenic-risk-score (PRS) of 41 CLL-susceptibility variants has been found to be associated with CLL risk among individuals of European-ancestry(EA). Here, we evaluate these variants, the PRS, and environmental factors for MBL risk. We also evaluate these variants and the CLL-PRS among African-American (AA) and EA-CLL cases and controls. Our study included 560 EA MBLs, 869 CLLs (696 EA/173 AA), and 2866 controls (2631 EA/235 AA). We used logistic regression, adjusting for age and sex, to estimate odds ratios (OR) and 95% confidence intervals within each race. We found significant associations with MBL risk among 21 of 41 variants and with the CLL-PRS (OR = 1.86, P = 1.9 × 10-29, c-statistic = 0.72). Little evidence of any association between MBL risk and environmental factors was observed. We observed significant associations of the CLL-PRS with EA-CLL risk (OR = 2.53, P = 4.0 × 10-63, c-statistic = 0.77) and AA-CLL risk (OR = 1.76, P = 5.1 × 10-5, c-statistic = 0.62). Inherited genetic factors and not environmental are associated with MBL risk. In particular, the CLL-PRS is a strong predictor for both risk of MBL and EA-CLL, but less so for AA-CLL supporting the need for further work in this population.

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.

Immunoglobulins (Ig) are produced by B lymphocytes as secreted antibodies or as part of the B-cell receptor. There is tremendous diversity of potential Ig transcripts (>1 × 10(12)) as a result of hundreds of germ-line gene segments, random nucleotide incorporation during joining of gene segments into a complete transcript, and the process of somatic hypermutation at individual nucleotides. This recombination and mutation process takes place in the maturing B cell and is responsible for the diversity of potential epitope recognition. Cancers arising from mature B cells are characterized by clonal production of Ig heavy (IGH@) and light chain transcripts, although whether the sequence has undergone somatic hypermutation is dependent on the maturation stage at which the neoplastic clone arose. Chronic lymphocytic leukemia (CLL) is the most common leukemia in adults and arises from a mature B cell with either mutated or unmutated IGH@ transcripts, the latter having worse prognosis and the assessment of which is routinely performed in the clinic. Currently, IGHV mutation status is assessed by Sanger sequencing and comparing the transcript to known germ-line genes. In this paper, we demonstrate that complete IGH@ V-D-J sequences can be computed from unselected RNA-seq reads with results equal or superior to the clinical procedure: in the only discordant case, the clinical transcript was out-of-frame. Therefore, a single RNA-seq assay can simultaneously yield gene expression profile, SNP and mutation information, as well as IGHV mutation status, and may one day be performed as a general test to capture multidimensional clinically relevant data in CLL.

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.

Patients with chronic lymphocytic leukemia (CLL) have high plasma-levels of Wnt5a, which can induce phosphorylation of ERK1/2 and enhance CLL-cell proliferation. Such effects could be inhibited by treatment with an ERK1/2 inhibitor, ERK1/2-specific siRNA, or cirmtuzumab, an anti-ROR1 mAb. The CLL-derived line, MEC1, expresses Wnt5a, but not ROR1. MEC1 cells transfected to express ROR1 (MEC1-ROR1) had higher levels of phosphorylated ERK1/2 than parental MEC1, or MEC1 transfected with ROR1ΔPRD, a truncated ROR1 lacking the cytoplasmic proline-rich domain (PRD), or ROR1P808A a mutant ROR1 with a P→A substitution at 808, which is required for complexing with the Rac-specific-guanine-nucleotide-exchange factor DOCK2 upon stimulation with Wnt5a. We silenced DOCK2 with siRNA and found this repressed the capacity of Wnt5a to induce ERK1/2 phosphorylation in MEC1-ROR1 or CLL cells. CLL cells that expressed ROR1 had higher levels of phosphorylated ERK1/2 or DOCK2 than CLL cells lacking ROR1. Although we found ibrutinib could inhibit the phosphorylation of ERK1/2 and DOCK2 induced by B-cell-receptor ligation, we found that this drug was unable to inhibit Wnt5a-induced, ROR1-dependent phosphorylation of ERK1/2 or DOCK2. This study demonstrates that Wnt5a can induce activation of ERK1/2 and enhance CLL-cell proliferation via a ROR1/DOCK2-dependent pathway independent of BTK.

Although numerous studies highlighted the role of Epstein-Barr Virus (EBV) in B-cell transformation, the involvement of EBV proteins or genome in the development of the most frequent adult leukemia, chronic lymphocytic leukemia (CLL), has not yet been defined. We hypothesized that EBV microRNAs contribute to progression of CLL and demonstrated the presence of EBV miRNAs in B-cells, in paraffin-embedded bone marrow biopsies and in the plasma of patients with CLL by using three different methods (small RNA-sequencing, quantitative reverse transcription PCR [q-RT-PCR] and miRNAs in situ hybridization [miRNA-ISH]). We found that EBV miRNA BHRF1-1 expression levels were significantly higher in the plasma of patients with CLL compared with healthy individuals (p < 0 · 0001). Notably, BHRF1-1 as well as BART4 expression were detected in the plasma of either seronegative or seropositive (anti-EBNA-1 IgG and EBV DNA tested) patients; similarly, miRNA-ISH stained positive in bone marrow specimens while LMP1 and EBER immunohistochemistry failed to detect viral proteins and RNA. We also found that BHRF1-1 plasma expression levels were positively associated with elevated beta-2-microglobulin levels and advanced Rai stages and observed a correlation between higher BHRF1-1 expression levels and shorter survival in two independent patients' cohorts. Furthermore, in the majority of CLL cases where BHRF1-1 was exogenously induced in primary malignant B cells the levels of TP53 were reduced. Our findings suggest that EBV may have a role in the process of disease progression in CLL and that miRNA RT-PCR and miRNAs ISH could represent additional methods to detect EBV miRNAs in patients with CLL.

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.

Cirmtuzumab may enhance the therapeutic activity of ibrutinib by inhibiting ROR1-dependent signaling pathway in patients with chronic lymphocytic leukemia (CLL). Mantle cell lymphoma (MCL) is B-cell malignancy that also expresses ROR1. In this study, we found that the plasma of patients with MCL had high levels of Wnt5a, a ROR1 ligand, that were comparable to those found in patients with CLL; in contrast Wnt5a was virtually undetectable in the plasma of age-matched healthy adults. We also found that Wnt5a induced Rac1 activation in the primary MCL cells. Cirmtuzumab, but not ibrutinib, could inhibit the capacity of Wnt5a to induce primary MCL cells to activate Rac1. Addition of exogenous Wnt5a in vitro significantly enhanced the numbers of MCL cell divisions and the proportion of dividing MCL cells entering S/G2 in MCL cells over time in the presence of CD154 and IL-4/10. Treatment of the MCL cells with cirmtuzumab, but not ibrutinib, blocked Wnt5a-enhanced proliferation of MCL cells. This study indicates that cirmtuzumab and ibrutinib may have complementary activity in the treatment of patients with MCL.

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.