CoChat
Sign Up
Cited by 178Chronic Lymphocytic Leukemia ResearchLymphoma Diagnosis and TreatmentCancer Genomics and Diagnostics

The DOI matches its Crossref record; 53 of the 53 references it cites were found in a registry or an index.

The evolutionary landscape of chronic lymphocytic leukemia treated with ibrutinib targeted therapy

Open access
Landau, Dan A., Sun, Clare, Rosebrock, Daniel, Herman, Sarah E. M., Fein, Joshua, Sivina, Mariela, Underbayev, Chingiz, Liu, Delong, Hoellenriegel, Julia, Ravichandran, Sarangan, Farooqui, Mohammed Z. H., Zhang, Wandi, Cibulskis, Carrie, Zviran, Asaf, Neuberg, Donna S., Livitz, Dimitri, Bozic, Ivana, Leshchiner, Ignaty, Getz, Gad, Burger, Jan A., Wiestner, Adrian, Wu, Catherine J.·2017·Nature Communications
DOI 10.1038/s41467-017-02329-yPMC PMC5736707
Open this paper in CoChat — free accountOpen PDF

Abstract

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.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41467-017-02329-y

Cites

53 of 53 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1 Putative driver gene mutations and copy number alterations at treatment initiation. a Treatment schema, absolute lymphocyte count (ALC) and number of samples per cohort that underwent whole-exome sequencing (WES) and RNA-sequencing at the indicated time points. Box plot shows the median, the interquartile range (IQR), and Tukey whiskers (±1.5 times IQR). b Distribution of clonal (black) and subclonal (tan) putative CLL driver mutations and copy number alterations (blue) across the 61 patients
Fig. 1 Putative driver gene mutations and copy number alterations at treatment initiation. a Treatment schema, absolute lymphocyte count (ALC) and number of samples per cohort that underwent whole-exome sequencing (WES) and RNA-sequencing at the indicated time points. Box plot shows the median, the interquartile range (IQR), and Tukey whiskers (±1.5 times IQR). b Distribution of clonal (black) and subclonal (tan) putative CLL driver mutations and copy number alterations (blue) across the 61 patientsCC BY · PubMed Central
Fig. 2 Early clonal shifts during treatment with ibrutinib. a 31% of CLLs showed significant clonal shifts during the first year of therapy, excluding relapse timepoints to quantify early clonal shifts in the absence of disease progression. Examples show an evolved (left) and unevolved (right) CLL. For each of the two CLLs, n-dimensional clustering across multiple mutations and timepoints (left panel) allows the derivation of the clonal cancer cell fraction (CCF) for each subclone at each timepoint, integrating the CCF information of each of the mutations harbored by the subclone (top panel). Comparison of clonal CCF between pre-treatment and the latest pre-relapse samples enables to determine the presence of significant clonal shifts over time (middle panel) and phylogenetic tree inference (bottom). Candidate CLL driver mutations6 are highlighted in yellow. b The magnitude of change in CCF of the largest rising or falling clone within each patient over the duration studied is displayed as a function of the CCF change. Red data points indicate significant (q < 0.1, see Methods section for P value estimation) changes in subclone size and blue indicating non-significant changes. Median values of CCF change are shown as dotted lines. c Association between clinical characteristics and the risk of early clonal shifts (n = 61, except for U-IGHV (n = 57) and persistent lymphocytosis at 6 months (n = 60) due to data availability). FISH fluorescent in situ hybridization, U-IGHV unmutated IGHV gene, CI confidence interval
Fig. 2 Early clonal shifts during treatment with ibrutinib. a 31% of CLLs showed significant clonal shifts during the first year of therapy, excluding relapse timepoints to quantify early clonal shifts in the absence of disease progression. Examples show an evolved (left) and unevolved (right) CLL. For each of the two CLLs, n-dimensional clustering across multiple mutations and timepoints (left panel) allows the derivation of the clonal cancer cell fraction (CCF) for each subclone at each timepoint, integrating the CCF information of each of the mutations harbored by the subclone (top panel). Comparison of clonal CCF between pre-treatment and the latest pre-relapse samples enables to determine the presence of significant clonal shifts over time (middle panel) and phylogenetic tree inference (bottom). Candidate CLL driver mutations6 are highlighted in yellow. b The magnitude of change in CCF of the largest rising or falling clone within each patient over the duration studied is displayed as a function of the CCF change. Red data points indicate significant (q < 0.1, see Methods section for P value estimation) changes in subclone size and blue indicating non-significant changes. Median values of CCF change are shown as dotted lines. c Association between clinical characteristics and the risk of early clonal shifts (n = 61, except for U-IGHV (n = 57) and persistent lymphocytosis at 6 months (n = 60) due to data availability). FISH fluorescent in situ hybridization, U-IGHV unmutated IGHV gene, CI confidence intervalCC BY · PubMed Central
Fig. 3 Genotype specific kinetics and disease progression on treatment. a Absolute CCF change per day (d|CCF|/dt) in the first period (first 30 days) and second period (beyond 30 days and up to 365 days) of significantly changing clones in patients with early clonal shifts. Box plot shows the median, the interquartile range (IQR), and Tukey whiskers (±1.5 times IQR, paired Welch’s t test). b The deceleration of change in CCF over the duration studied (dv|CCF|/dt) of these clones. c Measured circulating clonal sizes are shown for cases A35 and B04 (black dots represent the mode of CCF distribution multiplied by the number of circulating CLL cells, bars show 95% CI, CLL cell number is given in thousands of cells per microliter). The color curve represents fit with exponential growth with associated R 2 value. Total tumor burden is represented by the black line (ALC × sample purity) and insets show the corresponding changes in CCF for each CLL. d CCF change over time is shown for each putative CLL driver (somatic copy number alterations (SCNAs) and somatic single-nucleotide variations/indels (SSNVs)) represented by at least six instances across the patient cohorts. The bar on the right side of each plot reflects total number of clones harboring the indicated putative driver, with red, blue and gray representing the rate of significant CCF increase, decrease and no change, respectively. e The rate of change in CCF over time (dCCF/dt) does not demonstrate a clear trend by subclonal driver, as shown here in both the first time period (left) and second time period (right). f Kaplan–Meier plot of time-to-progression separated by the presence (red) or absence (black) of early clonal shifts (CCF change >0.1 with Q < 0.1)
Fig. 3 Genotype specific kinetics and disease progression on treatment. a Absolute CCF change per day (d|CCF|/dt) in the first period (first 30 days) and second period (beyond 30 days and up to 365 days) of significantly changing clones in patients with early clonal shifts. Box plot shows the median, the interquartile range (IQR), and Tukey whiskers (±1.5 times IQR, paired Welch’s t test). b The deceleration of change in CCF over the duration studied (dv|CCF|/dt) of these clones. c Measured circulating clonal sizes are shown for cases A35 and B04 (black dots represent the mode of CCF distribution multiplied by the number of circulating CLL cells, bars show 95% CI, CLL cell number is given in thousands of cells per microliter). The color curve represents fit with exponential growth with associated R 2 value. Total tumor burden is represented by the black line (ALC × sample purity) and insets show the corresponding changes in CCF for each CLL. d CCF change over time is shown for each putative CLL driver (somatic copy number alterations (SCNAs) and somatic single-nucleotide variations/indels (SSNVs)) represented by at least six instances across the patient cohorts. The bar on the right side of each plot reflects total number of clones harboring the indicated putative driver, with red, blue and gray representing the rate of significant CCF increase, decrease and no change, respectively. e The rate of change in CCF over time (dCCF/dt) does not demonstrate a clear trend by subclonal driver, as shown here in both the first time period (left) and second time period (right). f Kaplan–Meier plot of time-to-progression separated by the presence (red) or absence (black) of early clonal shifts (CCF change >0.1 with Q < 0.1)CC BY · PubMed Central
Fig. 4 Transcriptomic changes on treatment. a RNA-seq principal component analysis showing PC1 and PC2 across three time points. b Heat map of differentially expressed genes (fold change >2; FDR < 0.1; paired Student t test). 155 genes increased on drug; 498 genes deceased on drug. c Mean (±SEM) change in gene expression compared to pre-treatment across genes upregulated (top) and downregulated (bottom) on ibrutinib therapy. P values determined by a paired Student t test d Mean (±SEM) change in gene expression in experimentally derived gene sets that were enriched in genes downregulated on ibrutinib across patients. Asterisks represent significant difference in gene expression between 1 month and 6 months as determined by a paired Student t test. Genes sets in blue represent immune-receptor signaling, red represent cytokine signaling and brown represent general cellular responses. In all panels, green color is assigned to pre-treatment samples, purple to 1 month samples and orange 6 months samples in 14 patients totaling 42 samples
Fig. 4 Transcriptomic changes on treatment. a RNA-seq principal component analysis showing PC1 and PC2 across three time points. b Heat map of differentially expressed genes (fold change >2; FDR < 0.1; paired Student t test). 155 genes increased on drug; 498 genes deceased on drug. c Mean (±SEM) change in gene expression compared to pre-treatment across genes upregulated (top) and downregulated (bottom) on ibrutinib therapy. P values determined by a paired Student t test d Mean (±SEM) change in gene expression in experimentally derived gene sets that were enriched in genes downregulated on ibrutinib across patients. Asterisks represent significant difference in gene expression between 1 month and 6 months as determined by a paired Student t test. Genes sets in blue represent immune-receptor signaling, red represent cytokine signaling and brown represent general cellular responses. In all panels, green color is assigned to pre-treatment samples, purple to 1 month samples and orange 6 months samples in 14 patients totaling 42 samplesCC BY · PubMed Central
Fig. 5 PLCG2 mutations in ibrutinib resistant CLL. a The relapse characteristics are provided for the entire cohort. Patients without progressive disease (PD), are shown in gray with the time from treatment initiation to the last follow-up. For patients with PD, in addition to the time-to-progression, we provide the resistant genotype information. b Map of the PLCG2 gene with mutations identified in cases of ibrutinib resistance14,16,17,29. Red circles denote the number of patients with indicated mutations identified in the current study. Gray bars denote the regions covered by targeted sequencing. Domains PH Pleckstrin homology, X-box phosphatidylinositol-specific phospholipase C X domain, SH2 1 C-terminal Src homology 2, SH2 2 N-terminal Src homology 2, SH3 Src homology 3, Y-box phosphatidylinositol-specific phospholipase C Y domain, C2 calcium-binding motif. c Detailed information is presented for the two cases in which WES revealed additional PLCG2 mutations. Top panel shows the absolute lymphocyte count (ALC) over the patient’s clinical course, as well as changes in CCF of subclones as depicted in the inferred phylogenetic tree. Bottom panel shows the inferred growth kinetics of the different subclones, including measurements with corresponding 95% CI, as well as the exponential growth curves with 95% CI as shaded area. The calculated growth or decline rates from the exponential growth curves as well as the corresponding R 2 fit with exponential growth dynamics is listed in the table in the bottom panel. *We note that clones with R 2 = 1.0 merely reflects that only two data points were available for fitting
Fig. 5 PLCG2 mutations in ibrutinib resistant CLL. a The relapse characteristics are provided for the entire cohort. Patients without progressive disease (PD), are shown in gray with the time from treatment initiation to the last follow-up. For patients with PD, in addition to the time-to-progression, we provide the resistant genotype information. b Map of the PLCG2 gene with mutations identified in cases of ibrutinib resistance14,16,17,29. Red circles denote the number of patients with indicated mutations identified in the current study. Gray bars denote the regions covered by targeted sequencing. Domains PH Pleckstrin homology, X-box phosphatidylinositol-specific phospholipase C X domain, SH2 1 C-terminal Src homology 2, SH2 2 N-terminal Src homology 2, SH3 Src homology 3, Y-box phosphatidylinositol-specific phospholipase C Y domain, C2 calcium-binding motif. c Detailed information is presented for the two cases in which WES revealed additional PLCG2 mutations. Top panel shows the absolute lymphocyte count (ALC) over the patient’s clinical course, as well as changes in CCF of subclones as depicted in the inferred phylogenetic tree. Bottom panel shows the inferred growth kinetics of the different subclones, including measurements with corresponding 95% CI, as well as the exponential growth curves with 95% CI as shaded area. The calculated growth or decline rates from the exponential growth curves as well as the corresponding R 2 fit with exponential growth dynamics is listed in the table in the bottom panel. *We note that clones with R 2 = 1.0 merely reflects that only two data points were available for fittingCC BY · PubMed Central
Fig. 6 Disease progression without BTK or PLCG2 mutations. As in Fig. 5, detailed information is presented for two cases in which WES revealed no PLCG2 or BTK mutations. Top panel shows the absolute lymphocyte count (ALC) over the patient’s clinical course, as well as changes in CCF of subclones as depicted in the inferred phylogenetic tree. Bottom panel shows the inferred growth kinetics of the different subclones, including measurements with corresponding CI, as well as the exponential growth curves with CI as shaded area. The calculated growth or decline rates from the exponential growth curves as well as the corresponding R 2 fit with exponential growth dynamics is listed in the table in the bottom panel. *We note that R 2 = 1.0 merely reflects that only two data points were available for fitting
Fig. 6 Disease progression without BTK or PLCG2 mutations. As in Fig. 5, detailed information is presented for two cases in which WES revealed no PLCG2 or BTK mutations. Top panel shows the absolute lymphocyte count (ALC) over the patient’s clinical course, as well as changes in CCF of subclones as depicted in the inferred phylogenetic tree. Bottom panel shows the inferred growth kinetics of the different subclones, including measurements with corresponding CI, as well as the exponential growth curves with CI as shaded area. The calculated growth or decline rates from the exponential growth curves as well as the corresponding R 2 fit with exponential growth dynamics is listed in the table in the bottom panel. *We note that R 2 = 1.0 merely reflects that only two data points were available for fittingCC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 53 references

Cited by

Sorted by citations

Items 101–178 of 178

Related

See all

Metadata

Authors
Landau, D. A., Sun, C., Rosebrock, D., Herman, S. E. M., Fein, J., Sivina, M., Underbayev, C., Liu, D., Hoellenriegel, J., Ravichandran, S., Farooqui, M. Z. H., Zhang, W., Cibulskis, C., Zviran, A., Neuberg, D. S., Livitz, D., Bozic, I., Leshchiner, I., Getz, G., Burger, J. A., Wiestner, A., & Wu, C. J.
Affiliations
  • Broad Institute (US)
  • National Institutes of Health (US)
  • Tel Aviv University (IL)
  • The University of Texas MD Anderson Cancer Center (US)
  • Frederick National Laboratory for Cancer Research (US)
  • Dana-Farber Cancer Institute (US)
  • New York Genome Center (US)
Corresponding author
  • Wiestner, Adrian
  • Wu, Catherine J.
ORCID
Venue
Nature Communications
Publisher
Springer Science and Business Media LLC
Published
2017-12-19
Volume
8
Issue
1
Pages
2185-2185
ISSN
2041-1723 (ISSN-L)
DOI
10.1038/s41467-017-02329-y
PMC
PMC5736707
OpenAlex
W2773946416
Cited by
176
References
53
Citations by year
2026: 3 · 2025: 10 · 2024: 19 · 2023: 19 · 2022: 16 · 2021: 30 · 2020: 28 · 2019: 34 · 2018: 16
FWCI
12.41
Type
journal article
Language
en
Open access
Licence
Locations
  • Nature Communications (publishedVersion, cc-by)
  • Europe PMC (PubMed Central) (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • Digital Access to Scholarship at Harvard (DASH) (Harvard University) (submittedVersion, cc-by)
  • RePEc: Research Papers in Economics (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
Funders
  • American Association for Cancer Research
  • Burroughs Wellcome Fund
  • Leukemia and Lymphoma Society
  • Entertainment Industry Foundation
  • Broad Institute (U54HG003067)
  • National Institutes of Health
  • Stand Up To Cancer
  • National Heart, Lung, and Blood Institute (R01HL116452, R01HL131768)
  • National Human Genome Research Institute (U54HG003067)
  • National Cancer Institute (U10CA180861, R01CA184922, R01CA216273, R01CA182461, P01CA206978)
  • National Institute of Environmental Health Sciences (K01ES025431)
Topics
  • Health Sciences › Medicine › Hematology › Chronic Lymphocytic Leukemia Research
  • Health Sciences › Medicine › Pathology and Forensic Medicine › Lymphoma Diagnosis and Treatment
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Cancer Research › Cancer Genomics and Diagnostics
Keywords
chronic lymphocytic leukemia, ibrutinib, clonal evolution, drug resistance, BTK mutations, PLCG2 mutations, B cell receptor signaling, energy metabolism, cell cycle, transcriptome, exome sequencing, treatment outcomes
MeSH terms
Rituximab, Exome Sequencing, Agammaglobulinaemia Tyrosine Kinase, Adenine, Adult, Aged, Aged, 80 and over, Antineoplastic Combined Chemotherapy Protocols, Female, Humans, Longitudinal Studies, Male, Middle Aged, Mutation, Piperidines, Prognosis, Protein-Tyrosine Kinases, Pyrazoles, Pyrimidines, Signal Transduction, Leukemia, Lymphocytic, Chronic, B-Cell, Down-Regulation, Gene Expression Regulation, Neoplastic, Treatment Outcome, Disease Progression, Drug Resistance, Neoplasm, Phospholipase C gamma, Molecular Targeted Therapy, Clonal Evolution

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review