Inferring clonal heterogeneity in cancer using SNP arrays and whole genome sequencing.

MOTIVATION: Clonal heterogeneity is common in many types of cancer, including chronic lymphocytic leukemia (CLL). Previous research suggests that the presence of multiple distinct cancer clones is associated with clinical outcome. Detection of clonal heterogeneity from high throughput data, such as sequencing or single nucleotide polymorphism (SNP) array data, is important for gaining a better understanding of cancer and may improve prediction of clinical outcome or response to treatment. Here, we present a new method, CloneSeeker, for inferring clinical heterogeneity from sequencing data, SNP array data, or both. RESULTS: We generated simulated SNP array and sequencing data and applied CloneSeeker along with two other methods. We demonstrate that CloneSeeker is more accurate than existing algorithms at determining the number of clones, distribution of cancer cells among clones, and mutation and/or copy numbers belonging to each clone. Next, we applied CloneSeeker to SNP array data from samples of 258 previously untreated CLL patients to gain a better understanding of the characteristics of CLL tumors and to elucidate the relationship between clonal heterogeneity and clinical outcome. We found that a significant majority of CLL patients appear to have multiple clones distinguished by copy number alterations alone. We also found that the presence of multiple clones corresponded with significantly worse survival among CLL patients. These findings may prove useful for improving the accuracy of prognosis and design of treatment strategies. AVAILABILITY AND IMPLEMENTATION: Code available on R-Forge: https://r-forge.r-project.org/projects/CloneSeeker/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Time-to-progression after front-line fludarabine, cyclophosphamide, and rituximab chemoimmunotherapy for chronic lymphocytic leukaemia: a retrospective, multicohort study.

BACKGROUND: Fludarabine, cyclophosphamide, and rituximab (FCR) has become a gold-standard chemoimmunotherapy regimen for patients with chronic lymphocytic leukaemia. However, the question remains of how to treat treatment-naive patients with IGHV-unmutated chronic lymphocytic leukaemia. We therefore aimed to develop and validate a gene expression signature to identify which of these patients are likely to achieve durable remissions with FCR chemoimmunotherapy. METHODS: We did a retrospective cohort study in two cohorts of treatment-naive patients (aged ≥18 years) with chronic lymphocytic leukaemia. The discovery and training cohort consisted of peripheral blood samples collected from patients treated at the University of Texas MD Anderson Cancer Center (Houston, TX, USA), who fulfilled the diagnostic criteria of the International Workshop on Chronic Lymphocytic Leukemia, had received at least three cycles of FCR chemoimmunotherapy, and had been treated between Oct 10, 2000, and Oct 26, 2006 (ie, the MDACC cohort). We did transcriptional profiling on samples obtained from the MDACC cohort to identify genes associated with time to progression. We did univariate Cox proportional hazards analyses and used significant genes to cluster IGHV-unmutated samples into two groups (intermediate prognosis and unfavourable prognosis). After using cross-validation to assess robustness, we applied the Lasso method to standardise the gene expression values to find a minimum gene signature. We validated this signature in an external cohort of treatment-naive patients with IGHV-unmutated chronic lymphocytic leukaemia enrolled on the CLL8 trial of the German Chronic Lymphocytic Leukaemia Study Group who were treated between July 21, 2003, and April 4, 2006 (ie, the CLL8 cohort). FINDINGS: The MDACC cohort consisted of 101 patients and the CLL8 cohort consisted of 109 patients. Using the MDACC cohort, we identified and developed a 17-gene expression signature that distinguished IGHV-unmutated patients who were likely to achieve a long-term remission following front-line FCR chemoimmunotherapy from those who might benefit from alternative front-line regimens (hazard ratio 3·83, 95% CI 1·94-7·59; p<0·0001). We validated this gene signature in the CLL8 cohort; patients with an unfavourable prognosis versus those with an intermediate prognosis had a cause-specific hazard ratio of 1·90 (95% CI 1·18-3·06; p=0·008). Median time to progression was 39 months (IQR 22-69) for those with an unfavourable prognosis compared with 59 months (28-84) for those with an intermediate prognosis. INTERPRETATION: We have developed a robust, reproducible 17-gene signature that identifies a subset of treatment-naive patients with IGHV-unmutated chronic lymphocytic leukaemia who might substantially benefit from treatment with FCR chemoimmunotherapy. We recommend testing the value of this gene signature in a prospective study that compares FCR treatment with newer alternative therapies as part of a randomised clinical trial. FUNDING: Chronic Lymphocytic Leukaemia Global Research Foundation and the National Institutes of Health/National Cancer Institute.

Trisomy 12 chronic lymphocytic leukemia expresses a unique set of activated and targetable pathways.

Although trisomy 12 (+12) chronic lymphocytic leukemia (CLL) comprises about 20% of cases, relatively little is known about its pathophysiology. These cases often demonstrate atypical morphological and immunophenotypic features, high proliferative rates, unmutated immunoglobulin heavy chain variable region genes, and a high frequency of NOTCH1 mutation. Patients with +12 CLL have an intermediate prognosis, and show higher incidences of thrombocytopenia, Richter transformation, and other secondary cancers. Despite these important differences, relatively few transcriptional profiling studies have focused on identifying dysregulated pathways that characterize +12 CLL, and most have used a hierarchical cytogenetic classification in which cases with more than one recurrent abnormality are categorized according to the abnormality with the poorest prognosis. In this study, we sought to identify protein-coding genes whose expression contributes to the unique pathophysiology of +12 CLL. To exclude the likely confounding effects of multiple cytogenetic abnormalities on gene expression, our +12 patient cohort had +12 as the sole abnormality. We profiled samples obtained from 147 treatment-naïve patients. We compared cases with +12 as the only cytogenetic abnormality to cases with only del(13q), del(11q), or diploid cytogenetics using independent discovery (n=97) and validation (n=50) sets. We demonstrate that CLL cases with +12 as the sole abnormality express a unique set of activated pathways compared to other cytogenetic subtypes. Among these pathways, we identify the NFAT signaling pathway and the immune checkpoint molecule, NT5E (CD73), which may represent new therapeutic targets.

LDOC1 mRNA is differentially expressed in chronic lymphocytic leukemia and predicts overall survival in untreated patients.

We previously identified LDOC1 as one of the most significantly differentially expressed genes in untreated chronic lymphocytic leukemia (CLL) patients with respect to the somatic mutation status of the immunoglobulin heavy-chain variable region genes. However, little is known about the normal function of LDOC1, its contribution to the pathophysiology of CLL, or its prognostic significance. In this study, we have investigated LDOC1 mRNA expression in a large cohort of untreated CLL patients, as well as in normal peripheral blood B-cell (NBC) subsets and primary B-cell lymphoma samples. We have confirmed that LDOC1 is dramatically down-regulated in mutated CLL cases compared with unmutated cases, and have identified a new splice variant, LDOC1S. We show that LDOC1 is expressed in NBC subsets (naive > memory), suggesting that it may play a role in normal B-cell development. It is also expressed in primary B-cell lymphoma samples, in which its expression is associated with somatic mutation status. In CLL, we show that high levels of LDOC1 correlate with biomarkers of poor prognosis, including cytogenetic markers, unmutated somatic mutation status, and ZAP70 expression. Finally, we demonstrate that LDOC1 mRNA expression is an excellent predictor of overall survival in untreated CLL patients.

microRNA fingerprinting of CLL patients with chromosome 17p deletion identify a miR-21 score that stratifies early survival.

Aberrant expression of microRNAs (miRNAs) has been associated with clinical outcome in patients with chronic lymphocytic leukemia (CLL). To identify a powerful and easily assessable miRNA bio-marker of prognosis and survival, we performed quantitative reverse-transcription polymerase chain reaction (qRT-PCR) profiling in 104 CLL patients with a well-defined chromosome 17p status, and we validated our findings with miRNA microarray data from an independent cohort of 80 patients. We found that miR-15a, miR-21, miR-34a, miR-155, and miR-181b were differentially expressed between CLLs with chromosome 17p deletion and CLLs with normal 17p and normal karyotype, and that miR-181b was down-regulated in therapy-refractory cases. miR-21 expression levels were significantly higher in patients with poor prognosis and predicted overall survival (OS), and miR-181b expression levels significantly predicted treatment-free survival. We developed a 21FK score (miR-21 qRT-PCR, fluorescence in situ hybridization, Karyotype) to stratify patients according to OS and found that patients with a low score had a significantly longer OS time. When we evaluated the relative power of the 21FK score with the most used prognostic factors, the score was the most significant in both CLL cohorts. We conclude that the 21FK score represents a useful tool for distinguishing between good-prognosis and poor-prognosis CLL patients.

The B cell antigen receptor in atypical chronic lymphocytic leukemia with t(14;19)(q32;q13) demonstrates remarkable stereotypy.

The t(14;19)(q32;q13) is a recurrent chromosomal translocation reported in a variety of B-cell leukemias and lymphomas, including chronic lymphocytic leukemia (CLL). CLL cases associated with t(14;19) often have atypical morphologic and immunophenotypic features and unmutated immunoglobulin heavy chain (IGH) variable region (V) genes, associated with an aggressive clinical course. We analyzed IGHV somatic mutation status and gene use in 11 patients with t(14;19)-positive CLL. All cases were unmutated, and the IGHV genes in 10 cases showed minimal deviation from germline sequences. In 7 of 11 patients, we found homologous heavy chain rearrangements using IGHV4-39; light chain analysis revealed identical IGKV1-39 use. Corresponding V-(D)-J sequences demonstrated remarkable stereotypy of the immunoglobulin heavy and kappa light chain complementarity determining region 3 (H/K CDR3) genes. These findings raise the possibility that specific antigen drive is involved in the clonal development and/or selection of t(14;19)(q32;q13)-positive CLL cells. Our findings support the hypothesis that stimulatory signals through specific antigen receptors may promote the expansion of either CLL precursor cells or CLL clones that harbor distinct chromosomal abnormalities.

Relevance of the immunoglobulin VH somatic mutation status in patients with chronic lymphocytic leukemia treated with fludarabine, cyclophosphamide, and rituximab (FCR) or related chemoimmunotherapy regimens.

Although immunoglobulin V(H) mutation status (IgV(H) MS) is prognostic in patients with chronic lymphocytic leukemia (CLL) who are treated with alkylating agents or single-agent fludarabine, its significance in the era of chemoimmunotherapy is not known. We determined the IgV(H) somatic mutation status (MS) in 177 patients enrolled in a phase 2 study of fludarabine, cyclophosphamide, and rituximab (FCR) and in 127 patients treated with subsequent chemoimmunotherapy protocols. IgV(H) MS did not impact significantly on the complete remission (CR) rate of patients receiving FCR or related regimens. However, CR duration was significantly shorter in patients with CLL that used unmutated IgV(H) than those whose CLL used mutated IgV(H) (TTP 47% vs 82% at 6 years, P < .001). In a multivariate model considering all baseline characteristics, IgV(H) MS emerged as the only determinant of remission duration (hazard ratio 3.8, P < .001). Our results suggest that postremission interventions should be targeted toward patients with unmutated IgV(H) status.

Identification and validation of biomarkers of IgV(H) mutation status in chronic lymphocytic leukemia using microfluidics quantitative real-time polymerase chain reaction technology.

To develop a model incorporating relevant prognostic biomarkers for untreated chronic lymphocytic leukemia patients, we re-analyzed the raw data from four published gene expression profiling studies. We selected 88 candidate biomarkers linked to immunoglobulin heavy-chain variable region gene (IgV(H)) mutation status and produced a reliable and reproducible microfluidics quantitative real-time polymerase chain reaction array. We applied this array to a training set of 29 purified samples from previously untreated patients. In an unsupervised analysis, the samples clustered into two groups. Using a cutoff point of 2% homology to the germline IgV(H) sequence, one group contained all 14 IgV(H)-unmutated samples; the other contained all 15 mutated samples. We confirmed the differential expression of 37 of the candidate biomarkers using two-sample t-tests. Next, we constructed 16 different models to predict IgV(H) mutation status and evaluated their performance on an independent test set of 20 new samples. Nine models correctly classified 11 of 11 IgV(H)-mutated cases and eight of nine IgV(H)-unmutated cases, with some models using three to seven genes. Thus, we can classify cases with 95% accuracy based on the expression of as few as three genes.

Biological validation of differentially expressed genes in chronic lymphocytic leukemia identified by applying multiple statistical methods to oligonucleotide microarrays.

Oligonucleotide microarrays are a powerful tool for profiling the expression levels of thousands of genes. Different statistical methods for identifying differentially expressed genes can yield different results. To our knowledge, no experimental test has been performed to decide which method best identifies genes that are truly differentially expressed. We applied three statistical methods (dChip, t-test on log-transformed data, and Wilcoxon test) to identify differentially expressed genes in previously untreated patients with chronic lymphocytic leukemia (CLL). We used a training set of Affymetrix Hu133A microarray data from 11 patients with unmutated immunoglobulin (Ig) heavy chain variable region (VH) genes and 8 patients with mutated Ig VH genes. Differential expression was validated using semiquantitative real-time polymerase chain reaction assays and by validating models to predict the somatic mutation status of an independent test set of nine CLL samples. The methods identified 144 genes that were differentially expressed between cases of CLL with unmutated compared with mutated Ig VH genes. Eighty genes were identified by Wilcoxon test, 60 by t-test, and 65 by dChip, but only 11 were identified by all three methods. Greater agreement was found between the t-test and the Wilcoxon test. Differential expression was validated by semiquantitative real-time polymerase chain reaction assays for 83% of individual genes, regardless of the statistical method. However, the Wilcoxon test gave the most accurate predictions on new samples, and dChip, the least accurate. We found that all three methods were equally good for finding differentially expressed genes, but they found different genes. The genes selected by the nonparametric Wilcoxon test are the most robust for predicting the status of new cases. A comprehensive list of all differentially expressed genes can only be obtained by combining the results of multiple statistical tests.

B-Lymphoblastic Leukemia in Patients With Chronic Lymphocytic Leukemia: A Report of Four Cases.

OBJECTIVES: B-lymphoblastic leukemia (B-LBL) arising in patients with chronic lymphocytic leukemia (CLL) is exceedingly rare and poorly characterized. METHODS: We describe four patients with CLL and concurrent or subsequent B-LBL diagnosed by morphologic, immunophenotypic, cytogenetic, and molecular analysis and reviewed the literature. RESULTS: In three patients, B-LBL followed CLL by 5 to 15 years, and in one patient, B-LBL was diagnosed simultaneously with CLL. In all cases, the CLL had a typical immunophenotype, and the B-LBL blasts showed an immature B-cell immunophenotype with expression of CD10, CD19, and TdT and absence of surface immunoglobulin. In two patients, B-LBL blasts harbored t(9;22)(q34;q11.2)/BCR-ABL1. We sequenced the IGHV genes in both CLL and B-LBL in two patients and showed that IGHV usage differed. CONCLUSIONS: Our data suggest that at least some cases of B-LBL arising in patients with CLL are independent, secondary neoplasms rather than a manifestation of histologic transformation.

Genomic variation by whole-genome SNP mapping arrays predicts time-to-event outcome in patients with chronic lymphocytic leukemia: a comparison of CLL and HapMap genotypes.

Genomic abnormalities, such as deletions in 11q22 or 17p13, are associated with poorer prognosis in patients with chronic lymphocytic leukemia (CLL). We hypothesized that unknown regions of copy number variation (CNV) affect clinical outcome and can be detected by array-based single-nucleotide polymorphism (SNP) genotyping. We compared SNP genotypes from 168 untreated patients with CLL with genotypes from 73 white HapMap controls. We identified 322 regions of recurrent CNV, 82 of which occurred significantly more often in CLL than in HapMap (CLL-specific CNV), including regions typically aberrant in CLL: deletions in 6q21, 11q22, 13q14, and 17p13 and trisomy 12. In univariate analyses, 35 of total and 11 of CLL-specific CNVs were associated with unfavorable time-to-event outcomes, including gains or losses in chromosomes 2p, 4p, 4q, 6p, 6q, 7q, 11p, 11q, and 17p. In multivariate analyses, six CNVs (ie, CLL-specific variations in 11p15.1-15.4 or 6q27) predicted time-to-treatment or overall survival independently of established markers of prognosis. Moreover, genotypic complexity (ie, the number of independent CNVs per patient) significantly predicted prognosis, with a median time-to-treatment of 64 months versus 23 months in patients with zero to one versus two or more CNVs, respectively (P = 3.3 × 10(-8)). In summary, a comparison of SNP genotypes from patients with CLL with HapMap controls allowed us to identify known and unknown recurrent CNVs and to determine regions and rates of CNV that predict poorer prognosis in patients with CLL.

A two-gene signature, SKI and SLAMF1, predicts time-to-treatment in previously untreated patients with chronic lymphocytic leukemia.

We developed and validated a two-gene signature that predicts prognosis in previously-untreated chronic lymphocytic leukemia (CLL) patients. Using a 65 sample training set, from a cohort of 131 patients, we identified the best clinical models to predict time-to-treatment (TTT) and overall survival (OS). To identify individual genes or combinations in the training set with expression related to prognosis, we cross-validated univariate and multivariate models to predict TTT. We identified four gene sets (5, 6, 12, or 13 genes) to construct multivariate prognostic models. By optimizing each gene set on the training set, we constructed 11 models to predict the time from diagnosis to treatment. Each model also predicted OS and added value to the best clinical models. To determine which contributed the most value when added to clinical variables, we applied the Akaike Information Criterion. Two genes were consistently retained in the models with clinical variables: SKI (v-SKI avian sarcoma viral oncogene homolog) and SLAMF1 (signaling lymphocytic activation molecule family member 1; CD150). We optimized a two-gene model and validated it on an independent test set of 66 samples. This two-gene model predicted prognosis better on the test set than any of the known predictors, including ZAP70 and serum ß2-microglobulin.

Chronic lymphocytic leukemia With t(2;14)(p16;q32) involves the BCL11A and IgH genes and is associated with atypical morphologic features and unmutated IgVH genes.

The t(2;14)(p16;q32) has been reported previously in only 12 cases of chronic lymphocytic leukemia/small lymphocytic lymphoma (CLL/SLL). The clinicopathologic features have been incompletely described. We describe 6 new cases of CLL/SLL with t(2;14)(p16;q32). All had marrow involvement, 4 had absolute lymphocytosis, 4 had lymphadenopathy, and 3 had hepatosplenomegaly. All showed atypical lymphocyte morphologic features with plasmacytoid differentiation and irregular nuclei; 3 had increased prolymphocytes. Flow cytometry demonstrated typical immunophenotypes in 5 and an atypical immunophenotype in 1. All expressed ZAP70; 5 assessed showed unmutated IgV(H) genes. Karyotyping identified t(2;14)(p16;q32) as the sole abnormality in 1, primary abnormality in 2, and part of a complex karyotype in 3. Fluorescence in situ hybridization analysis revealed BCL11A/IgH rearrangement in all. After chemotherapy, 3 patients died of disease and 3 were alive with disease (median follow-up, 80 months). We conclude that CLL/SLL with t(2;14) (p16;q32) and BCL11A/IgH rearrangement is characterized by atypical morphologic features and unmutated IgV(H) genes.

Array CGH analysis of chronic lymphocytic leukemia reveals frequent cryptic monoallelic and biallelic deletions of chromosome 22q11 that include the PRAME gene.

We used BAC array-based CGH to detect genomic imbalances in 187 CLL cases. Submicroscopic deletions of chromosome 22q11 were observed in 28 cases (15%), and the frequency of these deletions was second only to loss of the 13q14 region, the most common genomic aberration in CLL. Oligonucleotide-based array CGH analysis showed that the 22q11 deletions ranged in size from 0.34 Mb up to approximately 1 Mb. The minimally deleted region included the ZNF280A, ZNF280B, GGTLC2, and PRAME genes. Quantitative real-time PCR revealed that ZNF280A, ZNF280B, and PRAME mRNA expression was significantly lower in the 22q11 deletion cases compared to non-deleted cases.

High expression of activation-induced cytidine deaminase (AID) and splice variants is a distinctive feature of poor-prognosis chronic lymphocytic leukemia.

In chronic lymphocytic leukemia (CLL), analysis of immunoglobulin heavy chain variable regions for somatic hypermutation identifies 2 prognostic subsets, mutated and unmutated. Investigators have postulated that unmutated and mutated CLL arises from malignant transformation of pre- and post-germinal center (GC) B cells, respectively. Alternatively, unmutated cases may arise from B cells stimulated by T-cell-independent antigens or from GC B cells with inactive somatic hypermutation. Activation-induced cytidine deaminase (AID), a protein essential for somatic hypermutation, is expressed by GC B cells in which this process occurs. We investigated AID mRNA expression in 20 CLL cases. In 8 cases we detected high expression of wild-type AID mRNA and 2 splice variants; in 12 cases and 5 normal peripheral blood B-cell samples we detected no expression using standard conditions. Of 8 CLL cases that highly expressed AID, 7 were unmutated, suggesting that this subset may arise from GC-experienced B cells with inactive somatic hypermutation, and may predict prognosis.