Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

Evolutionary history of transformation from chronic lymphocytic leukemia to Richter syndrome

Abstract

Richter syndrome (RS) arising from chronic lymphocytic leukemia (CLL) exemplifies an aggressive malignancy that develops from an indolent neoplasm. To decipher the genetics underlying this transformation, we computationally deconvoluted admixtures of CLL and RS cells from 52 patients with RS, evaluating paired CLL–RS whole-exome sequencing data. We discovered RS-specific somatic driver mutations (including IRF2BP2, SRSF1, B2M, DNMT3A and CCND3), recurrent copy-number alterations beyond del(9p21)(CDKN2A/B), whole-genome duplication and chromothripsis, which were confirmed in 45 independent RS cases and in an external set of RS whole genomes. Through unsupervised clustering, clonally related RS was largely distinct from diffuse large B cell lymphoma. We distinguished pathways that were dysregulated in RS versus CLL, and detected clonal evolution of transformation at single-cell resolution, identifying intermediate cell states. Our study defines distinct molecular subtypes of RS and highlights cell-free DNA analysis as a potential tool for early diagnosis and monitoring.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Developing an analytic framework for detecting RS-specific clones.
Fig. 2: The landscape of putative driver mutations in RS.
Fig. 3: Tracing evolution of RS on targeted agent therapy.
Fig. 4: Molecular mechanisms underlying transformation to RS.
Fig. 5: Transformation to RS at single-cell resolution.
Fig. 6: cfDNA isolated from plasma of RS patients shows evidence of transformation.

Similar content being viewed by others

Data availability

WES, RNA-seq, WGS, scRNA-seq and cfDNA data are available at dbgap (https://www.ncbi.nlm.nih.gov/gap/) using accession number phs002458.v2.p1 (http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs002458.v2.p1). RNA-seq data from validation cohort can be accessed by registering for an EGA account (https://ega-archive.org/) and contacting the Data Access Committee under study EGAS00001005495 and accession number EGAD00001007922. The following figures have associated raw data: Fig. 2b–h, Fig. 4a,f, Fig. 5a,b,d,e, Fig. 6g, Extended Data Fig. 4a–x, Extended Data Fig. 5a–c, Extended Data Fig. 6a, Extended Data Fig. 8, Extended Data Fig. 9, and Extended Data Fig. 10.

Code availability

Code is available for CNVSingle (https://github.com/broadinstitute/CNVsingle) and PhylogicNDT CopyNumber2Tree (https://github.com/broadinstitute/PhylogicNDT).

References

  1. Offin, M. et al. Concurrent RB1 and TP53 alterations define a subset of EGFR-mutant lung cancers at risk for histologic transformation and inferior clinical outcomes. J. Thorac. Oncol. 14, 1784–1793 (2019).

    Article  CAS  Google Scholar 

  2. Volta, A. D. et al. Transformation of prostate adenocarcinoma into small-cell neuroendocrine cancer under androgen deprivation therapy: much is achieved but more information is needed. J. Clin. Oncol. 37, 350–351 (2019).

    Article  CAS  Google Scholar 

  3. Parikh, S. A., Kay, N. E. & Shanafelt, T. D. How we treat Richter syndrome. Blood 123, 1647–1657 (2014).

    Article  CAS  Google Scholar 

  4. Landau, D. A. et al. Mutations driving CLL and their evolution in progression and relapse. Nature 526, 525–530 (2015).

    Article  CAS  Google Scholar 

  5. Puente, X. S. et al. Non-coding recurrent mutations in chronic lymphocytic leukaemia. Nature 526, 519–524 (2015).

    Article  CAS  Google Scholar 

  6. Knisbacher, B. A. et al. Molecular map of chronic lymphocytic leukemia and its impact on outcome. Nat. Genet. 54, 1664–1674 (2022).

  7. Chigrinova, E. et al. Two main genetic pathways lead to the transformation of chronic lymphocytic leukemia to Richter syndrome. Blood 122, 2673–2682 (2013).

    Article  CAS  Google Scholar 

  8. Fabbri, G. et al. Genetic lesions associated with chronic lymphocytic leukemia transformation to Richter syndrome. J. Exp. Med. 210, 2273–2288 (2013).

    Article  CAS  Google Scholar 

  9. Klintman, J. et al. Genomic and transcriptomic correlates of Richter transformation in chronic lymphocytic leukemia. Blood 137, 2800–2816 (2021).

    Article  CAS  Google Scholar 

  10. Rossi, D. et al. The genetics of Richter syndrome reveals disease heterogeneity and predicts survival after transformation. Blood 117, 3391–3401 (2011).

    Article  CAS  Google Scholar 

  11. Taylor-Weiner, A. et al. DeTiN: overcoming tumor-in-normal contamination. Nat. Methods 15, 531–534 (2018).

    Article  CAS  Google Scholar 

  12. Leshchiner, I. et al. Comprehensive analysis of tumour initiation, spatial and temporal progression under multiple lines of treatment. Preprint at bioRxiv 508127 (2018).

  13. Lawrence, M. S. et al. Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature 499, 214–218 (2013).

    Article  CAS  Google Scholar 

  14. Mermel, C. H. et al. GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers. Genome Biol. 12, R41 (2011).

    Article  Google Scholar 

  15. Schmitz, R. et al. Genetics and pathogenesis of diffuse large B-cell lymphoma. N. Engl. J. Med. 378, 1396–1407 (2018).

    Article  CAS  Google Scholar 

  16. Chapuy, B. et al. Genomic analyses of PMBL reveal new drivers and mechanisms of sensitivity to PD-1 blockade. Blood 134, 2369–2382 (2019).

    Article  Google Scholar 

  17. Biran, A. et al. Activation of Notch and Myc Signaling via B-cell-restricted depletion of Dnmt3a generates a consistent murine model of chronic lymphocytic leukemia. Cancer Res. 81, 6117–6130 (2021).

    Article  CAS  Google Scholar 

  18. Mahajan, V. S. et al. B1a and B2 cells are characterized by distinct CpG modification states at DNMT3A-maintained enhancers. Nat. Commun. 12, 2208 (2021).

    Article  CAS  Google Scholar 

  19. Chapuy, B. et al. Molecular subtypes of diffuse large B cell lymphoma are associated with distinct pathogenic mechanisms and outcomes. Nat. Med. 24, 679–690 (2018).

    Article  CAS  Google Scholar 

  20. Challa-Malladi, M. et al. Combined genetic inactivation of beta2-Microglobulin and CD58 reveals frequent escape from immune recognition in diffuse large B cell lymphoma. Cancer Cell 20, 728–740 (2011).

    Article  CAS  Google Scholar 

  21. Sade-Feldman, M. et al. Resistance to checkpoint blockade therapy through inactivation of antigen presentation. Nat. Commun. 8, 1136 (2017).

    Article  Google Scholar 

  22. Gettinger, S. et al. Impaired HLA Class I antigen processing and presentation as a mechanism of acquired resistance to immune checkpoint inhibitors in lung cancer. Cancer Disco. 7, 1420–1435 (2017).

    Article  CAS  Google Scholar 

  23. Singh, K. et al. c-MYC regulates mRNA translation efficiency and start-site selection in lymphoma. J. Exp. Med. 216, 1509–1524 (2019).

    Article  CAS  Google Scholar 

  24. Lee, S. C. et al. Synthetic lethal and convergent biological effects of cancer-associated spliceosomal gene mutations. Cancer Cell 34, 225–241 e228 (2018).

    Article  CAS  Google Scholar 

  25. Edelmann, J. et al. Genomic alterations in high-risk chronic lymphocytic leukemia frequently affect cell cycle key regulators and NOTCH1-regulated transcription. Haematologica 105, 1379–1390 (2020).

    Article  CAS  Google Scholar 

  26. Anderson, M. A. et al. Clinicopathological features and outcomes of progression of CLL on the BCL2 inhibitor venetoclax. Blood 129, 3362–3370 (2017).

    Article  CAS  Google Scholar 

  27. Jain, P. et al. Long-term outcomes for patients with chronic lymphocytic leukemia who discontinue ibrutinib. Cancer 123, 2268–2273 (2017).

    Article  CAS  Google Scholar 

  28. Maddocks, K. J. et al. Etiology of ibrutinib therapy discontinuation and outcomes in patients with chronic lymphocytic leukemia. JAMA Oncol. 1, 80–87 (2015).

    Article  Google Scholar 

  29. Burger, J. A. et al. Clonal evolution in patients with chronic lymphocytic leukaemia developing resistance to BTK inhibition. Nat. Commun. 7, 11589 (2016).

    Article  CAS  Google Scholar 

  30. Guieze, R. et al. Mitochondrial reprogramming underlies resistance to BCL-2 inhibition in lymphoid malignancies. Cancer Cell 36, 369–384 e313 (2019).

    Article  CAS  Google Scholar 

  31. Kasar, S. et al. Whole-genome sequencing reveals activation-induced cytidine deaminase signatures during indolent chronic lymphocytic leukaemia evolution. Nat. Commun. 6, 8866 (2015).

    Article  CAS  Google Scholar 

  32. Gruber, M. et al. Growth dynamics in naturally progressing chronic lymphocytic leukaemia. Nature 570, 474–479 (2019).

    Article  CAS  Google Scholar 

  33. Zhang, N. et al. Overexpression of Separase induces aneuploidy and mammary tumorigenesis. Proc. Natl Acad. Sci. USA 105, 13033–13038 (2008).

    Article  CAS  Google Scholar 

  34. Patel, A. P. et al. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science 344, 1396–1401 (2014).

    Article  CAS  Google Scholar 

  35. Bergen, V., Lange, M., Peidli, S., Wolf, F. A. & Theis, F. J. Generalizing RNA velocity to transient cell states through dynamical modeling. Nat. Biotechnol. 38, 1408–1414 (2020).

    Article  CAS  Google Scholar 

  36. Adalsteinsson, V. A. et al. Scalable whole-exome sequencing of cell-free DNA reveals high concordance with metastatic tumors. Nat. Commun. 8, 1324 (2017).

    Article  Google Scholar 

  37. Nadeu, F. et al. Detection of early seeding of Richter transformation in chronic lymphocytic leukemia. Nat. Med. 28, 1662–1671 (2022).

    Article  CAS  Google Scholar 

  38. Miller, C. R. et al. Near-tetraploidy is associated with Richter transformation in chronic lymphocytic leukemia patients receiving ibrutinib. Blood Adv. 1, 1584–1588 (2017).

    Article  CAS  Google Scholar 

  39. Quinton, R. J. et al. Whole-genome doubling confers unique genetic vulnerabilities on tumour cells. Nature 590, 492–497 (2021).

    Article  CAS  Google Scholar 

  40. Soilleux, E. J. et al. Diagnostic dilemmas of high-grade transformation (Richter’s syndrome) of chronic lymphocytic leukaemia: results of the phase II National Cancer Research Institute CHOP-OR clinical trial specialist haemato-pathology central review. Histopathology 69, 1066–1076 (2016).

    Article  Google Scholar 

  41. Favini, C. et al. Clonally unrelated Richter syndrome are truly de novo diffuse large B-cell lymphomas with a mutational profile reminiscent of clonally related Richter syndrome. Br. J. Haematol. 198, 1016–1022 (2022).

  42. Parikh, A. R. et al. Liquid versus tissue biopsy for detecting acquired resistance and tumor heterogeneity in gastrointestinal cancers. Nat. Med. 25, 1415–1421 (2019).

    Article  CAS  Google Scholar 

  43. McKenna, A. et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20, 1297–1303 (2010).

    Article  CAS  Google Scholar 

  44. Kim, S. et al. Strelka2: fast and accurate calling of germline and somatic variants. Nat. Methods 15, 591–594 (2018).

    Article  CAS  Google Scholar 

  45. Gao, G. F. et al. Tangent normalization for somatic copy-number inference in cancer genome analysis. Bioinformatics 38, 4677–4686 (2022).

    Article  Google Scholar 

  46. Consortium, I.T.P.-C.A.o.W.G. Pan-cancer analysis of whole genomes. Nature 578, 82–93 (2020).

    Article  Google Scholar 

  47. Chen, X. et al. Manta: rapid detection of structural variants and indels for germline and cancer sequencing applications. Bioinformatics 32, 1220–1222 (2016).

    Article  CAS  Google Scholar 

  48. Wala, J. A. et al. SvABA: genome-wide detection of structural variants and indels by local assembly. Genome Res 28, 581–591 (2018).

    Article  CAS  Google Scholar 

  49. Bass, A. J. et al. Genomic sequencing of colorectal adenocarcinomas identifies a recurrent VTI1A-TCF7L2 fusion. Nat. Genet. 43, 964–968 (2011).

    Article  CAS  Google Scholar 

  50. Drier, Y. et al. Somatic rearrangements across cancer reveal classes of samples with distinct patterns of DNA breakage and rearrangement-induced hypermutability. Genome Res 23, 228–235 (2013).

    Article  CAS  Google Scholar 

  51. Sondka, Z. et al. The COSMIC Cancer Gene Census: describing genetic dysfunction across all human cancers. Nat. Rev. Cancer 18, 696–705 (2018).

    Article  CAS  Google Scholar 

  52. Alexandrov, L. B. et al. The repertoire of mutational signatures in human cancer. Nature 578, 94–101 (2020).

    Article  CAS  Google Scholar 

  53. Olshen, A. B., Venkatraman, E. S., Lucito, R. & Wigler, M. Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics 5, 557–572 (2004).

    Article  Google Scholar 

  54. Brunet, J. P., Tamayo, P., Golub, T. R. & Mesirov, J. P. Metagenes and molecular pattern discovery using matrix factorization. Proc. Natl Acad. Sci. USA 101, 4164–4169 (2004).

    Article  CAS  Google Scholar 

  55. Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21 (2013).

    Article  CAS  Google Scholar 

  56. Law, C. W., Chen, Y., Shi, W. & Smyth, G. K. voom: precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. 15, R29 (2014).

    Article  Google Scholar 

  57. Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl Acad. Sci. USA 102, 15545–15550 (2005).

    Article  CAS  Google Scholar 

  58. Zhang, Y., Parmigiani, G. & Johnson, W. E. ComBat-seq: batch effect adjustment for RNA-seq count data. NAR Genom. Bioinform 2, lqaa078 (2020).

    Article  Google Scholar 

  59. Monti, S., Tamayo, P., Mesirov, J. & Golub, T. Consensus clustering: a resampling-based method for class discovery and visualization of gene expression microarray data. Mach. Learn. 52, 91–118 (2003).

    Article  Google Scholar 

  60. Stuart, T. et al. Comprehensive integration of single-cell data. Cell 177, 1888–1902 e1821 (2019).

    Article  CAS  Google Scholar 

  61. McGinnis, C. S., Murrow, L. M. & Gartner, Z. J. DoubletFinder: doublet detection in single-cell RNA sequencing data using artificial nearest neighbors. Cell Syst. 8, 329–337 e324 (2019).

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We thank C. Hahn, E. Ten Hacken, W. Zhang, S. Gohil and L. Werner for helpful discussions. We thank C. Patterson, S. Pollock, O. Olive, C. J. Shaughnessy, F. Dao and H. Lyon for assistance in data collection and organization and S. Belkin and C. Birger for assistance in data storage. We thank T. Lehmberg, M. McDonough, C. Galler and M. Collins for assistance in sample collection and biobanking. We thank the patients, their families and the investigators of the clinical trials for providing samples and clinical data. This study was supported by NIH/NCI P01 CA206978 (to C.J.W. and G.G.) and NCI (1U10CA180861-01) (to C.J.W.). The work is partially supported by the Broad/IBM Cancer Resistance Research Project (I.L., G.G. and L.P.) and a grant from Force Hemato (R.G.). Individual support was provided by DDCF Physician-Scientist Fellowship (E.M.P.), Dana-Farber Flames FLAIR fellowship (E.M.P.), ASCO Conquer Cancer Young Investigator Award (E.M.P.), Fishman Family Fund (R.G. and C.L.), EMBO fellowship ALTF 14-2018 (B.A.K.), NCI Research Specialist Award R50CA251956 (S.L.) and NIH/NCI R21CA267527-01 (S.Y.). Additional research support was provided by NIH R01 CA213442 (J.R.B.), Melton Family Foundation (J.R.B.), NIH/NCI R01-CA236361 (T.J.K.) and the Deutsche Forschungsgemeinschaft (DFG) SFB1074 subprojects B10 (E.T.) and subprojects B1 and B2 (E.T., C.S. and S.S.)

Author information

Authors and Affiliations

Authors

Contributions

E.M.P., I.L., R.G., C.J., G.G., S.S. and C.J.W. designed and performed the experiments, analyzed the data and wrote the manuscript. E.M.P., N.P.-Z., A.J.A., T.H. and S.L. generated single-cell RNA-seq data. C. Lemvigh, E.M.P. and I.L. analyzed single-cell RNA-seq data, along with C. Levovitz, F.U. and K.R. R.G., E.T., C. Schneider, M.S.D., N.J., W.W., L.R., T.J.K, J.B., S.H., P.F., F.C., N.K., S.P., J.R.B. and S.S. provided patient samples. K.J.L. and S.L. designed targeted NGS assay for detecting the NOTCH1 3′-UTR mutation. N.R. performed mapping and analysis of these NGS data. D.R., F.U., C. Levovitz and S.Y. analyzed the RNA-seq data. S.D. analyzed the mutational data under the supervision of E.M.P. and C.J.W. C.M., J.M., J.H., L.L. and C. Stewart analyzed the WGS data. C.J., I.L., B.P., L.E., D.R., A.T-W., A.M., D.L., E.M.P., R.G. and C.J.W. performed sequencing data analyses, assessment of the clonal architecture and inference of phylogenies under the supervision of I.L. and G.G. E.M.P., R.G., L.R., J.B. and S.F. prepared patient samples. I.L. developed the analytic tool for determining somatic copy-number variations from FFPE samples and CNVSingle for detecting copy-number events in single-cell RNA-seq data. D.N. performed and supervised statistical analyses. R.R. performed statistical analyses. G.F. provided graphical representation of the clinical data. B.A.K. performed immunogenetic data analyses. B.P.D., K.R., C. Levovitz and L.P. helped to design and guide the research. B.P.D., R.A.J. and K.S. were involved in managing the project. I.L., C.J. and Z.L. performed cell-free DNA analyses. All authors discussed, interpreted results and approved the manuscript.

Corresponding authors

Correspondence to Gad Getz or Catherine J. Wu.

Ethics declarations

Competing interests

I.L. serves as a consultant for PACT Pharma, Inc.; has stock in, is on the board of and serves as a consultant for ennov1 LLC.; and is on the board of and holds equity in Nord Bio, Inc. C.J.W., G.G., B.A.K. and Z.L. are inventors on a patent: ‘Compositions, Panels, and Methods for Characterizing Chronic Lymphocytic Leukemia’ (PCT/US21/45144). C.J.W., G.G., E.M.P., I.L. and R.G. are named as inventors on U.S. provisional patent application serial number 63/244,625, filed on 15 September 2021, and U.S. provisional patent application serial number 63/291,213, filed on 17 December 2021, both of which are entitled ‘Diagnosis and Prognosis of Richter’s Syndrome.’ G.G. is a founder of and consultant for and holds privately held equity in Scorpion Therapeutics; receives funding support from IBM and Pharmacyclics; and is an inventor on patent applications related to MSMuTect, MSMutSig, MSIDetect, POLYSOLVER and SignatureAnalyzer-GPU. R.G. receives funding support from Abbvie, AstraZeneca, Gilead, Janssen and Roche. M.S.D. served as a consultant for Abbvie, Adaptive Biotechnologies, Ascentage Pharma, AstraZeneca, BeiGene, Bristol-Myers Squibb, Eli Lilly, Genentech/Roche, Janssen, Merck, Ono Pharmaceuticals, Pharmacyclics, Research to Practice, Takeda, TG Therapeutics, Verastem and Zentalis; receives funding support from Ascentage Pharma, AstraZeneca, Genentech/Roche, MEI Pharma, Novartis, Pharmacyclics, Surface Oncology, TG Therapeutics and Verastem; and receives funding for travel from Abbvie, BeiGene, BioAscend, Clinical Care Options, Curio Science, Imedex, ION Solutions, Janssen, MDOutlook, PeerView, PRIME Oncology, Research to Practice and TG Therapeutics. J.R.B. has served as a consultant for Abbvie, Acerta/AstraZeneca, BeiGene, Bristol-Myers Squibb/Juno/Celgene, Catapult, Eli Lilly, Genentech/Roche, Hutchmed, Janssen, MEI Pharma, MorphoSys AG, Novartis, Pfizer, Pharmacyclics and Rigel; and received research funding from Gilead, Loxo/Lilly, SecuraBio, Sun and TG Therapeutics. C.J.W. receives funding support from Pharmacyclics and holds equity in BioNTech, Inc. N.E.K. serves as an advisor for Abbvie, AstraZeneca, Beigene, Behring, Cytomx Therapy, Dava Oncology, Janssen, Juno Therapeutics, Oncotracker, Pharmacyclics and Targeted Oncology; receives funding support from Abbvie, Acerta Pharma, Bristol Meyer Squib, Celgene, Genentech, MEI Pharma, Pharmacyclics, Sunesis, TG Therapeutics and Tolero Pharmaceuticals; and participates on the Data Safety Monitoring Committee for Agios Pharm, AstraZeneca, BMS-Celgene, Cytomx Therapeutics, Dren Bio, Janssen, MorphoSys and Rigel. T.J.K. is on the advisory board and receives funding support from Abbvie and Roche and serves on the speakers’ bureau for Abbvie, Janssen and Roche. E.T. serves as an advisor and is on the speakers’ bureau for Abbvie, Janssen and Roche; and receives funding support from Abbvie and Roche. S.S. is on the advisory board and receives funding, travel support and speakers’ fees from AbbVie, AstraZeneca, BeiGene, BMS, Celgene, Gilead, GSK, Hoffmann-La Roche, Janssen, Novartis and Sunesis. N.J. receives research funding from AbbVie, Adaptive Biotechnologies, ADC Therapeutics, Aprea Therapeutics, AstraZeneca, BMS, Cellectis, Fate Therapeutics, Genentech, Incyte, Loxo Oncology, Medisix, Mingsight, Novalgen, Pfizer, Pharmacyclics, Precision BioSciences, Servier and Takeda; and serves on the advisory board/receives honoraria from AbbVie, Adaptive Biotechnologies, ADC Therapeutics, AstraZeneca, Beigene, BMS, Cellectis, Genentech, Janssen, MEI Pharma, Pharmacyclics, Precision BioSciences, Servier and TG Therapeutics. W.G.W. reports funding from Abbvie, AstraZeneca/Acerta Pharma, Cyclacel, Genentech, Gilead Sciences, GSK/Novartis, Janssen, Juno Therapeutics, KITE Pharma, Loxo Oncology, Inc., Miragen, Oncternal Therapeutics, Inc., Pharmacyclics LLC, Sunesis and Xencor. S.A.P. has received research funding to the institution from AbbVie, Ascentage Pharma, AstraZeneca, Janssen, Merck, Pharmacyclics and TG Therapeutics for clinical studies in which S.A.P. is a principal investigator. S.A.P has received honoraria for participation in consulting activities/advisory board meetings for AbbVie, Adaptive Biotechnologies, Amgen, AstraZeneca, Genentech, GlaxoSmithKline, Merck and Pharmacyclics (no personal compensation). K.J.L. holds equity in Standard BioTools, Inc. D.N. has stock ownership in Madrigal Pharmaceuticals. C.S. serves on a speakers’ bureau for AbbVie and AstraZeneca. D.L. holds stock in and consults for ennov1. N.P is currently an employee at Bristol Meyers Squibb. All remaining authors declare no competing interests.

Peer review

Peer review information

Nature Medicine thanks Daniel Hodson and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editor: Anna Maria Ranzoni, in collaboration with the Nature Medicine team.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Extended Data Fig. 1 Clonal deconvolution process.

a, distinguishing RS from CLL clones after inferring subclonal composition of paired CLL and RS samples. b, inferring phylogenetic tree from cancer cell fraction using PhylogicNDT. c, sample composition d, mapping copy-number variations to clones using CopyNumber2Tree.

Extended Data Fig. 2 Phylogenetic reconstruction and somatic genomic alterations.

For each of the patient trios with WES data, the left panel shows the phylogenetic tree tracing the transformation history from CLL to RS. The magenta frame denotes the Richter clones. The middle top panel represents the subclonal composition inferred after clustering alterations with similar cancer cell fractions as previously reported4. The middle bottom panel indicates the timeline with RS and CLL sampling time and CLL therapeutic lines. (F, fludarabine; C, cyclophosphamide; R, rituximab; P, pentostatin; O/Ofa, ofatumumab; HDMP, high-dose methylprednisolone; A, alemtuzumab; Auto, autologous stem cell transplantation; CLB, chlorambucil; B, bendamustine; CHOP, cyclophosphamide, doxorubicin, vincristine, prednisone; ESHAP, etoposide, methylprednisolone, high-dose cytarabine, cisplatin; CHP, cyclophosphamide, doxorubicin, prednisone; Len, lenalidomide; Ob, obinutuzumab; idela; idelalisisb; D, dexamethasone; Adria, adriamycin). The right panel is composed of allelic fraction plots and allelic copy ratio plots showing clonal assignment of somatic copy-number events to CLL and RS clones. Cases with whole genome doubling in Extended Data Fig. 2 and clonal unrelated cases in Extended Data Fig. 3.

Extended Data Fig. 3 Phylogenetic reconstruction and somatic genomic alterations.

For each of the patient trios with WES data, the left panel shows the phylogenetic tree tracing the transformation history from CLL to RS. The magenta frame denotes the Richter clones. The middle top panel represents the subclonal composition inferred after clustering alterations with similar cancer cell fractions as previously reported4. The middle bottom panel indicates the timeline with RS and CLL sampling time and CLL therapeutic lines. (F, fludarabine; C, cyclophosphamide; R, rituximab; P, pentostatin; O/Ofa, ofatumumab; HDMP, high-dose methylprednisolone; A, alemtuzumab; Auto, autologous stem cell transplantation; CLB, chlorambucil; B, bendamustine; CHOP, cyclophosphamide, doxorubicin, vincristine, prednisone; ESHAP, etoposide, methylprednisolone, high-dose cytarabine, cisplatin; CHP, cyclophosphamide, doxorubicin, prednisone; Len, lenalidomide; Ob, obinutuzumab; idela; idelalisisb; D, dexamethasone; Adria, adriamycin). The right panel is composed of allelic fraction plots and allelic copy ratio plots showing clonal assignment of somatic copy-number events to CLL and RS clones. Cases with whole genome doubling in Extended Data Fig. 2 and clonal unrelated cases in Extended Data Fig. 3.

Extended Data Fig. 4 Putative RS driver genes.

a-x, individual protein mutation maps for selected putative Richter drivers, showing gene mutation subtype (for example, missense), position and evidence of mutational hotspots. Panels were generated by using the cBioPortal for Cancer Genomics tool.

Extended Data Fig. 5 RS sCNAs and genomic clustering.

GISTIC2-defined recurrent copy-number gains (red, left) and losses (blue, right) are visualized for focal events for RS samples (a) and RS clones (b) (RS samples with CLL events subtracted, bottom). Chromosomes are shown on the vertical axis. Green line denotes a near significant q value of 0.25 and significant events (q<0.1) are annotated in text along with putative driver genes contained within the peak (Supplementary Table 5) c, NMF clustering of RS with DLBCL (304 de novo DLBCL samples19 shows clonal related RS clusters separately from DLBCL and closes to DLBCL from C2 (ref. 19). Clonal unrelated RS clusters across DLBCL subtypes and separate from RS. Samples were annotated for clonal relationship (related RS, gray, unrelated RS, black), cohort (DLBCL, light purple; RS, dark purple) and DLBCL clusters (C1, purple; C2, yellow, C3, pink, C4, blue, C5, green)19. d, NMF clustering of RS shows 5 distinct genomic subtypes of transformation.

Extended Data Fig. 6 Transcriptome supports distinct RS molecular subtypes.

a, Supervised clustering of transcriptome data from 36 RS patients by molecular subtype highlights differentially regulated genes in subtype 1 and 3 (Supplementary Table 8). Samples are annotated for cohort (Discovery, pink; Validation, yellow), clonal relationship (unrelated, black, related, white), and sample purity by WES (green gradient). b, Unsupervised consensus clustering of RS transcriptome data (n=36) shows 5 clusters. (Discovery, pink; Validation, yellow), RS molecular subtype (1, purple; 2, blue; 3, orange; 4, green; and 5, pink), and sample purity by WES (green gradient).c, 5 × 5 table showing association between molecular subtype of RS and unsupervised transcriptome clusters (2 sided Fisher’s exact test, P=0.038) d, Kaplan–Meier curve showing OS of clonal unrelated RS compared to clonal related RS. P value is log rank (2 sided Mantel Cox).

Extended Data Fig. 7 Phylogenetic trees showing CLL and RS clones from WGS of paired samples.

a, Phylogenetic tree and CCF plot for 9 patients based on WGS data showing clonal related RS (magenta box). b, Phylogenetic tree and CCF plot for 2 patients based on WGS demonstrating clonally unrelated RS c, Representative phylogenetic trees and CCF plot for 3 patients from UK cohort9 based on WGS.

Extended Data Fig. 8 WGS Circos plots with or without chromothripsis.

a, chromothripsis and kataegis in RS sample (Pt 42) with whole genome doubling. Circos plots showing structural variants (interchromosomal, blue; deletion, red; inversion, yellow; tandem duplication, green; long range, teal), allelic copy number (middle), rainfall plot with kataegis regions (red) and chromosomes (outside). Adjacent rainfall plots show kategis regions (C to G, red; C to T, yellow; C to A, teal) with corresponding allelic copy-number fragmentation. b, Circos plots from RS WGS samples showing structural variants (interchromosomal, blue; deletion, red; inversion, yellow; tandem duplication, green and long range, teal), allelic copy number (middle), rainfall plot with kataegis regions (red) and chromosomes (outside). SVs impacting known genes and translocation partners are labeled (Supplementary Table 7k).

Extended Data Fig. 9 Single cell processing and transcriptome analysis of RS samples at single cell resolution.

a, flow sorting strategy for RS single-cell samples. Flow sorting to separate RS and CLL cells by size for Patient 19 and Patient 41 (lymph node, LN; peripheral blood, PB; bone marrow, BM). Flow sorting viable cells for Pt 43, Pt 4 and Pt10. Representative flow plots below demonstrate CLL and RS cells were included in sorted population. b, B cell receptor (BCR) clonotypes plotted for RS and CLL clusters on UMAP visualization. c, Representative example from patient 10 showing CNVsingle identifies malignant B cell clusters (5 and 6) separate from immune cell clusters (0,1,2,3,4,7,9). d, UMI/cell and Gene/cell plots for CLL and RS single-cell clusters. RS demonstrates higher UMI/cell (P<2.2 × 10–16 see Methods, Supplementary Table 8). e, RNA inference of directional trajectories is shown on UMAP visualization for Pts 43 and 10. f, copy-number variation heatmap inferred in each cluster from scRNA-seq data using our CNVSingle algorithm for Pts 43 and 10 (Methods).

Extended Data Fig. 10 Single-cell transcriptome and copy-number analysis of RS patients.

a, UMAP visualization of single-cells from patient 4 (left) with associated allelic copy-number ratio plot inferred by CNVsingle (top right) and RS WES (bottom right). b, UMAP visualization of CLL and RS cells from Patient 18 (left top panel) with flow-sorting annotations (right top panel). Inferred CNAs from CNVSingle (bottom panel) are shown as heatmap with CLL (green) and RS (pink) events highlighted. c, UMAP visualization of CLL and RS cells from Patient 41 (left top panel) with flow-sorting annotations (right top panel). Inferred CNAs from CNVSingle (bottom panel) with CLL (green) and RS (pink) events highlighted. d, Plasma of patient 44 shows RS-specific sCNVs on chromosome 9 and 13 leading up to RS diagnosis, which are not reflected in circulating CLL e, Plasma of patient 99 at the start of CLL-directed therapy (top) and just ahead of diagnosis of RS (bottom) during CLL response. f, Chromothripsis in post-transplant RS plasma cfDNA at time of relapse (Pt 112). g, Plot showing allele frequency of RS (purple) and CLL (green) mutations in RS WES (bottom) and plasma cfDNA WES (top) for patient 5 (top) and patient 44 (bottom).

Supplementary information

Supplementary Information

Supplementary Figs. 1 and 2.

Reporting Summary

Supplementary Tables

Supplementary Tables 1–10 combined into single Excel document.

Rights and permissions

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Parry, E.M., Leshchiner, I., Guièze, R. et al. Evolutionary history of transformation from chronic lymphocytic leukemia to Richter syndrome. Nat Med 29, 158–169 (2023). https://doi.org/10.1038/s41591-022-02113-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41591-022-02113-6

This article is cited by

Search

Quick links

Nature Briefing: Cancer

Sign up for the Nature Briefing: Cancer newsletter — what matters in cancer research, free to your inbox weekly.

Get what matters in cancer research, free to your inbox weekly. Sign up for Nature Briefing: Cancer