Whole-genome sequencing identifies recurrent mutations in chronic lymphocytic leukaemia - Nature
- ️Campo, Elías
- ️Sun Jun 05 2011
- Letter
- Open access
- Published: 05 June 2011
- Magda Pinyol2,
- Víctor Quesada1,
- Laura Conde3,
- Gonzalo R. Ordóñez1,
- Neus Villamor3,
- Georgia Escaramis4,
- Pedro Jares3,
- Sílvia Beà3,
- Marcos González-Díaz5,
- Laia Bassaganyas4,
- Tycho Baumann6,
- Manel Juan7,
- Mónica López-Guerra3,
- Dolors Colomer3,
- José M. C. Tubío4,8,
- Cristina López3,
- Alba Navarro3,
- Cristian Tornador4,
- Marta Aymerich3,
- María Rozman3,
- Jesús M. Hernández5,
- Diana A. Puente1,
- José M. P. Freije1,
- Gloria Velasco1,
- Ana Gutiérrez-Fernández1,
- Dolors Costa3,
- Anna Carrió3,
- Sara Guijarro3,
- Anna Enjuanes3,
- Lluís Hernández3,
- Jordi Yagüe7,
- Pilar Nicolás9,
- Carlos M. Romeo-Casabona9,
- Heinz Himmelbauer10,
- Ester Castillo10,
- Juliane C. Dohm10,
- Silvia de Sanjosé11,
- Miguel A. Piris12,
- Enrique de Alava5,
- Jesús San Miguel5,
- Romina Royo13,
- Josep L. Gelpí13,
- David Torrents13,
- Modesto Orozco13,
- David G. Pisano14,
- Alfonso Valencia14,
- Roderic Guigó15,
- Mónica Bayés16,
- Simon Heath16,
- Marta Gut16,
- Peter Klatt17,
- John Marshall18,
- Keiran Raine18,
- Lucy A. Stebbings18,
- P. Andrew Futreal18,
- Michael R. Stratton18,
- Peter J. Campbell18,
- Ivo Gut16,
- Armando López-Guillermo6,
- Xavier Estivill4,
- Emili Montserrat6,
- Carlos López-Otín1 na1 &
- …
- Elías Campo3 na1
Nature volume 475, pages 101–105 (2011)Cite this article
-
30k Accesses
-
64 Altmetric
Subjects
Abstract
Chronic lymphocytic leukaemia (CLL), the most frequent leukaemia in adults in Western countries, is a heterogeneous disease with variable clinical presentation and evolution1,2. Two major molecular subtypes can be distinguished, characterized respectively by a high or low number of somatic hypermutations in the variable region of immunoglobulin genes3,4. The molecular changes leading to the pathogenesis of the disease are still poorly understood. Here we performed whole-genome sequencing of four cases of CLL and identified 46 somatic mutations that potentially affect gene function. Further analysis of these mutations in 363 patients with CLL identified four genes that are recurrently mutated: notch 1 (NOTCH1), exportin 1 (XPO1), myeloid differentiation primary response gene 88 (MYD88) and kelch-like 6 (KLHL6). Mutations in MYD88 and KLHL6 are predominant in cases of CLL with mutated immunoglobulin genes, whereas NOTCH1 and XPO1 mutations are mainly detected in patients with unmutated immunoglobulins. The patterns of somatic mutation, supported by functional and clinical analyses, strongly indicate that the recurrent NOTCH1, MYD88 and XPO1 mutations are oncogenic changes that contribute to the clinical evolution of the disease. To our knowledge, this is the first comprehensive analysis of CLL combining whole-genome sequencing with clinical characteristics and clinical outcomes. It highlights the usefulness of this approach for the identification of clinically relevant mutations in cancer.
Similar content being viewed by others
Main
To gain insights into the molecular alterations that cause CLL, we performed whole-genome sequencing of four cases representative of different forms of the disease: two cases, CLL1 and CLL2, with no mutations in the immunoglobulin genes (IGHV-unmutated) and two cases, CLL3 and CLL4, with mutations in these genes (IGHV-mutated) (Supplementary Table 1 and Supplementary Information). We used a combination of whole-genome sequencing and exome sequencing, as well as long-insert paired-end libraries, to detect variants in chromosomal structure (Supplementary Fig. 1 and Supplementary Tables 2–5). We obtained more than 99.7% concordance between whole-genome sequencing calls and genotyping data, indicating that the coverage and parameters used were sufficient to detect most of the sequence variants in these samples (Supplementary Information). We detected about 1,000 somatic mutations per tumour in non-repetitive regions (Fig. 1a, Supplementary Fig. 2 and Supplementary Table 6). These numbers of somatic mutations were lower than the numbers in melanoma and lung carcinoma5,6, but in agreement with previous estimates of less than one mutation per megabase (Mb) for leukaemias7. The most common substitution was the transition G>A/C>T, usually occurring in a CpG context (Fig. 1b and Supplementary Fig. 2). We also detected marked differences in the mutation pattern between CLL samples and these differences were associated with tumour subtype (Fig. 1b). Thus, IGHV-mutated cases showed a higher proportion of A>C/T>G mutations than cases with unmutated IGHV (16 ± 0.2% versus 6.2 ± 0.1%). The base preceding the adenine in A to C transversions showed an over-representation of thymine, when compared to the prevalence expected from its representation in non-repetitive sequences in the wild-type genome (P < 0.001, Fig. 1c), and there were fewer A to C substitutions at GpA dinucleotides than would be expected by chance (P < 0.001). These differences between CLL subtypes might reflect the molecular mechanisms implicated in their respective development. The pattern and context of mutations are consistent with their being introduced by the error-prone polymerase η during somatic hypermutation in immunoglobulin genes8. This indicates that polymerase η could contribute to the high frequency of A > T to C > G transversions in cases with IGHV-mutated. It also extends the differences observed between these two CLL subtypes to the genomic level.
a, Distribution of somatic alterations. For each tumour genome, copy number (solid lines), density of mutations per 5-Mb window (bars) and protein-coding mutations (dots) are shown. The shaded rectangle indicates the location of the 13q14 deletion that was present in three of the four CLL cases. Chromosome numbers are listed below the four profiles. b, Frequency of substitutions in each CLL tumour for the six possible classes of mutation. c, Distribution of the four possible NpA dinucleotides for the A to C transversion in each tumour genome, compared with the expected distribution across the genome. The total number of A to C substitutions per case is indicated at the top (**, P < 0.001).
We classified the somatic mutations into three different classes according to their potential functional effect (Supplementary Information). We also searched for small insertions and deletions (indels) in coding regions: we found and validated five somatic indels, which caused frameshifts in protein-coding regions (Supplementary Table 7). We identified 46 mutations that changed the protein-coding sequences of 45 genes in the four patients analysed (Supplementary Table 7). None of these nucleotide substitutions had been previously linked to CLL and among the five indel mutations, only one, in NOTCH1 (p.P2515Rfs*4), had been previously found in various lymphoid malignancies, including CLL9,10. To determine whether any of these 45 genes was mutated in more than one CLL case, we analysed an initial validation set of 169 CLL patients. We focused on the 26 genes that are expressed at the RNA level in CLL cells (Supplementary Table 7) because mutations in expressed genes are more likely to have a biological effect than those in non-expressed genes. We used a pooled-sequencing strategy that led us to identify four genes with at least one additional mutation in the validation series: these were NOTCH1, MYD88, XPO1 and KLHL6 (Table 1 and Supplementary Information).
Analysis of additional CLL cases revealed that the deletion of a CT dinucleotide in NOTCH1 (p.P2515Rfs*4) was found in 29 of 255 patients and two additional mutations in the same region were also found (p.Q2503* and p.F2482Ffs*2) (Fig. 2a, b). Accordingly, NOTCH1 is mutated in 12% of CLL patients (Supplementary Table 8). These mutations generate a premature stop codon, resulting in a NOTCH1 protein lacking the C-terminal domain, which contains a PEST sequence (a sequence rich in proline, glutamic acid, serine and threonine) (Fig. 2a). Removal of this region results in the accumulation of an active protein isoform in the mutated CLL cells (Fig. 2c and Supplementary Fig. 3). NOTCH1 is constitutively expressed in CLL11, but the NOTCH1 mutations identified herein generate a more stable and active isoform of the protein. Gene expression analysis of ten NOTCH1-mutated and 49 unmutated CLL cases revealed a high number of differentially expressed genes (n = 542, false discovery rate <0.05; Supplementary Table 9). Likewise, in a gene-set analysis, we found that there was significant differential expression of the NOTCH1 signalling pathway12 and two metabolic pathways (oxidative phosphorylation and glycolysis/gluconeogenesis). This is consistent with the NOTCH1-mediated activation of multiple biosynthetic routes in T acute lymphoblastic leukaemia13. When the differential expression of individual genes from the NOTCH1 pathway was analysed, 23 of the 46 genes assigned to this pathway12 showed a significant differential expression (P < 0.05) in NOTCH1-mutated CLL (Fig. 2d). NOTCH1-mutated patients had a more advanced clinical stage at diagnosis, more adverse biological features and an overall survival that was significantly shorter than those with NOTCH1 unmutated (10-yr overall survival: 21% versus 56%, P = 0.03; Fig. 2e, f). NOTCH1-mutated CLL also underwent transformation into diffuse large B-cell lymphoma more frequently than NOTCH1-unmutated CLL (7 of 31 cases, 23%, versus 3 of 224 cases, 1.3%; P < 0.001). The same IGHV clonal rearrangement and NOTCH1 mutation were found in the CLL and corresponding transformed diffuse large B-cell lymphoma of the four cases studied, indicating a clonal relationship of both components.
a, Schematic representation of human NOTCH1, showing the main domains and locations of the three different somatic mutations identified in CLL. NEC, NOTCH1 extracellular subunit; NTM, NOTCH1 transmembrane subunit; ICN, intracellular domain of NOTCH1; LNR, Lin-12 NOTCH repeats; RAM, RAM domain; ANK, ankyrin repeat domain; PEST, PEST domain. b, Electropherogram showing the heterozygous CT deletion recurrently identified in CLL. c, Western blot showing NOTCH1 protein levels in CLL cases with or without the NOTCH1 p.P2515Rfs*4 mutation, and in Jurkat cells as a control. The arrow indicates the band corresponding to the NTM; the large arrowhead indicates the smaller band corresponding to the mutant form. d, Heat map showing the 23 genes of the NOTCH1 pathway that are differentially expressed in NOTCH1-mutated versus non-mutated CLL. e, Distribution of disease stage (Binet), ZAP-70 expression status, CD38 expression status and IGHV mutational status (UM, unmutated IGHV) in patients with or without mutations in NOTCH1 (*, P < 0.02; **, P < 0.01). f, Actuarial probability of overall survival of CLL patients with mutated or unmutated NOTCH1 (*, P = 0.03).
A recurrent mutation (p.L265P) in the MYD88 gene (Fig. 3a, b) was also identified in 9 of 310 CLL patients (2.9%). During revision of this manuscript, the same mutation has been identified in different lymphomas14, highlighting its relevance in the pathogenesis of lymphoid neoplasias. This protein participates in the signalling pathways of interleukin-1 and Toll-like receptors during the immune response15. MyD88 immunoprecipitation from CLL cells with the p.L265P mutation resulted in the co-immunoprecipitation of large amounts of IRAK1, in contrast to cells lacking this mutation (Fig. 3c). Other effectors of this signalling pathway, including STAT3, IκBα and NF-κB p65 subunit, showed higher phosphorylation in MYD88-mutated than in unmutated CLL cells (Fig. 3d, e) and there was an increased DNA-binding activity of NF-κB in MYD88-mutated cells (Supplementary Fig. 4). These data support the hypothesis that the MYD88 p.L265P mutation constitutes an activating mutation of this novel proto-oncogene14,16. Stimulation of interleukin-1 receptor or Toll-like receptors in MYD88-mutated CLL cells induced the secretion of 5-fold to 150-fold higher levels of interleukin 1 receptor antagonist (IL1RN, also known as IL1RA), interleukin 6 and chemokine (C-C motif) ligands 2, 3 and 4 (CCL2, CCL3 and CCL4), when compared to the secretion of these cytokines by MYD88-unmutated CLLs. Cytokine secretion was elevated in MYD88-mutated cells in response to stimulation of at least four of the eight TLRs tested. No response was observed in lymphocytes carrying the inactivating MYD88 mutation E52DEL (Fig. 3f and Supplementary Fig. 5). The high production of these cytokines has been implicated in the recruitment of macrophages and T lymphocytes by CLL cells, creating a favourable niche for their survival17. Moreover, activation of Toll-like receptors in CLL cells promotes the proliferation of tumour cells and protects them from spontaneous apoptosis18. Patients with MYD88-mutated CLL were diagnosed at a younger age than those with wild-type MYD88 (median 43 yr, range 38–63, versus median 63 yr, range 27–94; P < 0.001) and the disease presented with a more advanced clinical stage (Fig. 3g), although no differences were observed in progression or survival rates. Notably, almost all patients with the MYD88 p.L265P mutation (seven of the eight evaluated) belonged to the IGHV-mutated group.
a, Multiple sequence alignment of MyD88 around the mutated residue (arrow) in different species. Cons., degree of conservation. b, Electropherogram showing the recurrent heterozygous p.L265P MYD88 mutation (arrow) detected in CLL. c, Cell extracts from a MYD88-mutated CLL (L265P) and a MYD88-unmutated CLL (WT) were immunoprecipitated with anti-MyD88 antibody. The immunoprecipitated and unbound fractions were analysed by western blot using anti-IRAK1 and anti-MyD88 antibodies. d, Western blot analysis of phosphorylated STAT3 (p-STAT3 (Tyr 705)) and total STAT3 in cell extracts from MYD88-mutated or unmutated CLL tumour cells. β-Actin was used as a control to show equal loading. e, Western blots showing phosphorylated IκBα (p-IκBα), total IκBα, phosphorylated p65 subunit of NF-κB (p-p65) and total p65 subunit of NF-κB in cell extracts from MYD88-mutated or unmutated CLL tumour cells. f, Heat map representing the cytokine levels secreted by B cells from eight individuals after Toll-like receptor stimulation. Only the five cytokines that showed the most significant differences between MYD88-mutated and MYD88-unmutated CLL are shown. ‘E52DEL’ indicates B cells from two patients with an inactivating MYD88 mutation, ‘WT’ corresponds to tumour cells from CLL patients without MYD88 mutation and ‘L265P’ indicates tumour cells from patients carrying a mutated MYD88. The stimulation experiments for each of the Toll-like receptors (TLRs) are represented in different colours. NS, no stimulus. g, Distribution of disease stage (Binet), ZAP-70 expression status, CD38 expression status and IGHV mutational status (UM, unmutated IGHV) in patients according to the presence or absence of p.L265P MYD88 mutation (*, P < 0.03).
We also identified four cases with mutations in the same codon of the exportin 1 gene (XPO1; p.E571K and p.E571G). Exportin 1 is implicated in the nuclear export of proteins and mRNAs in yeast, including members of the MAP kinase pathway19. The fact that the same residue is mutated in four CLL cases and is part of a highly conserved region (Supplementary Fig. 6) indicates that the mutation affects XPO1 activity. Notably, all four cases with mutations in XPO1 belonged to the IGHV-unmutated subtype and two of them also had the p.P2515Rfs*4 mutation in NOTCH1, indicating that both mutations could have synergic effects in CLL development.
We identified three patients carrying a total of six mutations (F49L/L65P, L90F and L58P/T64A/Q81P) in the gene encoding kelch-like protein 6 (KLHL6), which is implicated in the formation of the germinal centre during B cell maturation20. All six mutations were clustered between residues 49 and 90 (Supplementary Fig. 7). The presence of several point mutations in cis, located near the transcriptional start site of a gene that is highly expressed in the germinal centre, is a characteristic feature of somatic hypermutation. In fact, all three patients had CLL with mutated IGHV. Although somatic hypermutation occurs mainly in IGHV regions, other proto-oncogenes, including BCL6, MYC and PIM1, are mutated by somatic hypermutation in different lymphomas21. However, only BCL6 has been previously shown to be hypermutated by this mechanism in CLL21. Our data show that KLHL6 is probably also a target of somatic hypermutation in IGHV-mutated patients, although its precise contribution to the oncogenic process in CLL remains to be determined.
In addition to these four genes, we identified a series of large genomic alterations that were previously reported2. They included the deletion, in three cases, of the 13q14 region22, and a 40-Mb deletion in chromosome 6q14–q22 (Fig. 1a, Supplementary Fig. 1 and Supplementary Table 5). Finally, in one patient we detected a p.P281R mutation in the cyclin D2 gene (CCND2), which resulted in the accumulation of cyclin D2 in tumour cells (Supplementary Fig. 8). This finding, together with the high conservation of this residue and the identification of mutations in the equivalent residue of cyclin D1 (CCND1) in endometrial cancer23, indicates that this CCND2 mutation could be a driver contributing to the development of CLL in this patient. The finding illustrates the putative relevance of non-recurrent mutations for the pathogenesis of CLL.
The International Cancer Genome Consortium project was founded on the concept that sequencing of cancer genomes could reshape our understanding of cancer biology, with direct implications for clinical translation24. Our study of four CLL genomes underscores this transformative potential, although additional studies will be necessary to translate these findings to the clinic. We have identified four recurrently mutated genes and provided novel insights into the mechanisms by which leukaemic cells recruit, instruct and coordinate a tumour microenvironment. Currently, the biological identification of different subgroups of CLL is based on markers such as IGHV mutational status, cytogenetics, ZAP-70 expression or CD38 expression, which are not fundamental agents in the leukaemic process. The classification of patients based on genomic drivers of the disease is conceptually appealing, as shown by our demonstration that NOTCH1 and MYD88 mutations identify distinct subgroups of patients with particular clinical and biological features. Furthermore, we provide functional evidence that both NOTCH1 and MYD88 mutations are activating events and potential therapeutic targets. The potential to personalize therapeutic choices for patients on the basis of the genomic architecture of their cancers is the long-term aspiration for studies such as this, combining whole-genome sequencing, functional studies and clinical analysis of patients with cancer.
Methods Summary
Four patients with CLL, who had given informed consent for sample collection and analysis, were studied. Tumour samples were obtained before treatment and tumour cells were separated from non-tumour cells by immunomagnetic depletion of T cells, natural killer cells, monocytes and granulocytes (Supplementary Information). Tumour cell purity was ≥98% as assessed by flow cytometry. Normal blood cells from the same patient were obtained after treatment, resulting in no detectable, or less than 0.05%, tumour cell contamination, as assessed by flow cytometry. Additional samples from 363 patients were obtained for clinical validation. Protocols for long-insert and short-insert library construction and for massively parallel paired-end sequencing have been described elsewhere (ref. 25 and Supplementary Information). Genotyping and copy number analysis were performed using the Affymetrix SNP6.0, Agilent 1M and Illumina OmniQuad arrays on the same cases used for whole-genome sequencing. For the validation of candidate genes in a set of 169 additional CLL patients, we used a combination of PCR amplification and Illumina sequencing in pooled samples, resulting in efficient identification of germline and somatic mutations (Supplementary Information). Sequencing data were aligned to the human reference genome (GRCh37) using Burrows–Wheeler alignment (BWA)26 and somatic substitutions were identified using Sidrón, a probabilistic binomial model that uses genotyping data to calibrate sequencing error per sample. Functional analyses of the identified mutations were performed using cryopreserved primary tumour cells. For gene expression analysis, RNA was purified from tumour cells and analysed using the HU133 plus 2.0 GeneChip (Affymetrix). For immunoprecipitation and western blotting, CLL cell extracts were prepared and detected using the indicated antibodies (Supplementary Information). For Toll-like receptor stimulation of CLL cells, the Human TLR1–9 agonist kit (InvivoGen) was used.
Accession codes
Data deposits
Sequencing, expression and genotyping array data have been deposited at the European Genome-Phenome Archive (EGA, http://www.ebi.ac.uk/ega/), which is hosted at the European Bioinformatics Institute (EBI), under accession number EGAS00000000092.
References
Rozman, C. & Montserrat, E. Chronic lymphocytic leukemia. N. Engl. J. Med. 333, 1052–1057 (1995)
Zenz, T., Mertens, D., Kuppers, R., Dohner, H. & Stilgenbauer, S. From pathogenesis to treatment of chronic lymphocytic leukaemia. Nature Rev. Cancer 10, 37–50 (2010)
Damle, R. N. et al. Ig V gene mutation status and CD38 expression as novel prognostic indicators in chronic lymphocytic leukemia. Blood 94, 1840–1847 (1999)
Hamblin, T. J., Davis, Z., Gardiner, A., Oscier, D. G. & Stevenson, F. K. Unmutated Ig V(H) genes are associated with a more aggressive form of chronic lymphocytic leukemia. Blood 94, 1848–1854 (1999)
Pleasance, E. D. et al. A comprehensive catalogue of somatic mutations from a human cancer genome. Nature 463, 191–196 (2010)
Pleasance, E. D. et al. A small-cell lung cancer genome with complex signatures of tobacco exposure. Nature 463, 184–190 (2010)
Greenman, C. et al. Patterns of somatic mutation in human cancer genomes. Nature 446, 153–158 (2007)
Spencer, J. & Dunn-Walters, D. K. Hypermutation at A-T base pairs: the A nucleotide replacement spectrum is affected by adjacent nucleotides and there is no reverse complementarity of sequences flanking mutated A and T nucleotides. J. Immunol. 175, 5170–5177 (2005)
Weng, A. P. et al. Activating mutations of NOTCH1 in human T cell acute lymphoblastic leukemia. Science 306, 269–271 (2004)
Sportoletti, P. et al. NOTCH1 PEST domain mutation is an adverse prognostic factor in B-CLL. Br. J. Haematol. 151, 404–406 (2010)
Rosati, E. et al. Constitutively activated Notch signaling is involved in survival and apoptosis resistance of B-CLL cells. Blood 113, 856–865 (2009)
Kanehisa, M., Goto, S., Furumichi, M., Tanabe, M. & Hirakawa, M. KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 38, D355–D360 (2010)
Palomero, T. et al. NOTCH1 directly regulates c-MYC and activates a feed-forward-loop transcriptional network promoting leukemic cell growth. Proc. Natl Acad. Sci. USA 103, 18261–18266 (2006)
Ngo, V. N. et al. Oncogenically active MYD88 mutations in human lymphoma. Nature 470, 115–119 (2011)
O’Neill, L. A. & Bowie, A. G. The family of five: TIR-domain-containing adaptors in Toll-like receptor signalling. Nature Rev. Immunol. 7, 353–364 (2007)
Coste, I. et al. Dual function of MyD88 in RAS signaling and inflammation, leading to mouse and human cell transformation. J. Clin. Invest. 120, 3663–3667 (2010)
Burger, J. A. et al. High-level expression of the T-cell chemokines CCL3 and CCL4 by chronic lymphocytic leukemia B cells in nurselike cell cocultures and after BCR stimulation. Blood 113, 3050–3058 (2009)
Muzio, M. et al. Expression and function of toll like receptors in chronic lymphocytic leukaemia cells. Br. J. Haematol. 144, 507–516 (2009)
Ferrigno, P., Posas, F., Koepp, D., Saito, H. & Silver, P. A. Regulated nucleo/cytoplasmic exchange of HOG1 MAPK requires the importin β homologs NMD5 and XPO1. EMBO J. 17, 5606–5614 (1998)
Kroll, J. et al. The BTB-kelch protein KLHL6 is involved in B-lymphocyte antigen receptor signaling and germinal center formation. Mol. Cell. Biol. 25, 8531–8540 (2005)
Pasqualucci, L. et al. Hypermutation of multiple proto-oncogenes in B-cell diffuse large-cell lymphomas. Nature 412, 341–346 (2001)
Aqeilan, R. I., Calin, G. A. & Croce, C. M. miR-15a and miR-16–1 in cancer: discovery, function and future perspectives. Cell Death Differ. 17, 215–220 (2010)
Moreno-Bueno, G. et al. Cyclin D1 gene (CCND1) mutations in endometrial cancer. Oncogene 22, 6115–6118 (2003)
Hudson, T. J. et al. International network of cancer genome projects. Nature 464, 993–998 (2010)
Bentley, D. R. et al. Accurate whole human genome sequencing using reversible terminator chemistry. Nature 456, 53–59 (2008)
Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009)
Acknowledgements
This work was funded by the Spanish Ministry of Science and Innovation (MICINN) through the Instituto de Salud Carlos III (ISCIII) and Red Temática de Investigación del Cáncer (RTICC) del ISCIII. C.L.-O. is an Investigator of the Botin Foundation and D.T., of the ICREA program. We thank E. Santos for his support of this project, A. Carracedo and J. Benítez for genotyping studies, C. Fortuny for the supply of samples and N. Villahoz and M. C. Muro for their work in the coordination of the CLL-ICGC Consortium. We are also grateful to all patients with CLL who participated in this study.
Author information
Author notes
Carlos López-Otín and Elías Campo: These authors contributed equally to this work.
Authors and Affiliations
Departamento de Bioquímica y Biología Molecular, Instituto Universitario de Oncología, Universidad de Oviedo, 33006 Oviedo, Spain,
Xose S. Puente, Víctor Quesada, Gonzalo R. Ordóñez, Diana A. Puente, José M. P. Freije, Gloria Velasco, Ana Gutiérrez-Fernández & Carlos López-Otín
Unidad de Genómica, Institut d'Investigacions Biomèdiques August Pi i Sunyer (IDIBAPS), 08036 Barcelona, Spain ,
Magda Pinyol
Unidad de Hematopatología, Servicio de Anatomía Patológica, Hospital Clínic, Universitat de Barcelona, IDIBAPS, 08036 Barcelona, Spain ,
Laura Conde, Neus Villamor, Pedro Jares, Sílvia Beà, Mónica López-Guerra, Dolors Colomer, Cristina López, Alba Navarro, Marta Aymerich, María Rozman, Dolors Costa, Anna Carrió, Sara Guijarro, Anna Enjuanes, Lluís Hernández & Elías Campo
Genes and Disease Programme, Center for Genomic Regulation, Pompeu Fabra University (CRG-UPF), 08003 Barcelona, Spain ,
Georgia Escaramis, Laia Bassaganyas, José M. C. Tubío, Cristian Tornador & Xavier Estivill
Servicio de Hematología, Hospital Universitario, Centro de Investigación del Cáncer-IBMCC (USAL-CSIC), Universidad de Salamanca, 37007 Salamanca, Spain ,
Marcos González-Díaz, Jesús M. Hernández, Enrique de Alava & Jesús San Miguel
Servicio de Hematología, Hospital Clínic, IDIBAPS, Universidad de Barcelona, 08036 Barcelona, Spain ,
Tycho Baumann, Armando López-Guillermo & Emili Montserrat
Servicio de Inmunología, Hospital Clínic, IDIBAPS, 08036 Barcelona, Spain ,
Manel Juan & Jordi Yagüe
Hospital Clínico Universitario de Santiago de Compostela, 15706 Santiago de Compostela, Spain ,
José M. C. Tubío
Cátedra Inter-Universitaria de Derecho y Genoma Humano, Universidad de Deusto, Universidad del País Vasco, 48007 Bilbao, Spain ,
Pilar Nicolás & Carlos M. Romeo-Casabona
Ultrasequencing Unit, Center for Genomic Regulation, Pompeu Fabra University (CRG-UPF), 08003 Barcelona, Spain ,
Heinz Himmelbauer, Ester Castillo & Juliane C. Dohm
Unidad de Infecciones y Cáncer, Institut Català de Oncologia-IDIBELL, 08907 Hospitalet de Llobregat, Spain
Silvia de Sanjosé
Molecular Pathology Programme, Spanish National Cancer Research Centre (CNIO), 28029 Madrid, Spain ,
Miguel A. Piris
Programa Conjunto de Biología Computacional, Barcelona Supercomputing Center (BSC), Institut de Recerca Biomèdica (IRB), Spanish National Bioinformatics Institute, Universitat de Barcelona, 08028 Barcelona, Spain ,
Romina Royo, Josep L. Gelpí, David Torrents & Modesto Orozco
Structural Biology and Biocomputing Programme, Spanish National Cancer Research Centre (CNIO), Spanish National Bioinformatics Institute, 28029 Madrid, Spain ,
David G. Pisano & Alfonso Valencia
Genomic Bioinformatics Programme, Center for Genomic Regulation, Pompeu Fabra University (CRG-UPF), Spanish National Bioinformatics Institute, 08003 Barcelona, Spain ,
Roderic Guigó
Centro Nacional de Análisis Genómico, Parc Científic de Barcelona, 08028 Barcelona, Spain ,
Mónica Bayés, Simon Heath, Marta Gut & Ivo Gut
Centro Nacional de Biotecnología, Consejo Superior de Investigaciones Científicas, 28049 Madrid, Spain ,
Peter Klatt
Wellcome Trust Sanger Institute, Hinxton CB10 1SA, UK ,
John Marshall, Keiran Raine, Lucy A. Stebbings, P. Andrew Futreal, Michael R. Stratton & Peter J. Campbell
Authors
- Xose S. Puente
You can also search for this author in PubMed Google Scholar
- Magda Pinyol
You can also search for this author in PubMed Google Scholar
- Víctor Quesada
You can also search for this author in PubMed Google Scholar
- Laura Conde
You can also search for this author in PubMed Google Scholar
- Gonzalo R. Ordóñez
You can also search for this author in PubMed Google Scholar
- Neus Villamor
You can also search for this author in PubMed Google Scholar
- Georgia Escaramis
You can also search for this author in PubMed Google Scholar
- Pedro Jares
You can also search for this author in PubMed Google Scholar
- Sílvia Beà
You can also search for this author in PubMed Google Scholar
- Marcos González-Díaz
You can also search for this author in PubMed Google Scholar
- Laia Bassaganyas
You can also search for this author in PubMed Google Scholar
- Tycho Baumann
You can also search for this author in PubMed Google Scholar
- Manel Juan
You can also search for this author in PubMed Google Scholar
- Mónica López-Guerra
You can also search for this author in PubMed Google Scholar
- Dolors Colomer
You can also search for this author in PubMed Google Scholar
- José M. C. Tubío
You can also search for this author in PubMed Google Scholar
- Cristina López
You can also search for this author in PubMed Google Scholar
- Alba Navarro
You can also search for this author in PubMed Google Scholar
- Cristian Tornador
You can also search for this author in PubMed Google Scholar
- Marta Aymerich
You can also search for this author in PubMed Google Scholar
- María Rozman
You can also search for this author in PubMed Google Scholar
- Jesús M. Hernández
You can also search for this author in PubMed Google Scholar
- Diana A. Puente
You can also search for this author in PubMed Google Scholar
- José M. P. Freije
You can also search for this author in PubMed Google Scholar
- Gloria Velasco
You can also search for this author in PubMed Google Scholar
- Ana Gutiérrez-Fernández
You can also search for this author in PubMed Google Scholar
- Dolors Costa
You can also search for this author in PubMed Google Scholar
- Anna Carrió
You can also search for this author in PubMed Google Scholar
- Sara Guijarro
You can also search for this author in PubMed Google Scholar
- Anna Enjuanes
You can also search for this author in PubMed Google Scholar
- Lluís Hernández
You can also search for this author in PubMed Google Scholar
- Jordi Yagüe
You can also search for this author in PubMed Google Scholar
- Pilar Nicolás
You can also search for this author in PubMed Google Scholar
- Carlos M. Romeo-Casabona
You can also search for this author in PubMed Google Scholar
- Heinz Himmelbauer
You can also search for this author in PubMed Google Scholar
- Ester Castillo
You can also search for this author in PubMed Google Scholar
- Juliane C. Dohm
You can also search for this author in PubMed Google Scholar
- Silvia de Sanjosé
You can also search for this author in PubMed Google Scholar
- Miguel A. Piris
You can also search for this author in PubMed Google Scholar
- Enrique de Alava
You can also search for this author in PubMed Google Scholar
- Jesús San Miguel
You can also search for this author in PubMed Google Scholar
- Romina Royo
You can also search for this author in PubMed Google Scholar
- Josep L. Gelpí
You can also search for this author in PubMed Google Scholar
- David Torrents
You can also search for this author in PubMed Google Scholar
- Modesto Orozco
You can also search for this author in PubMed Google Scholar
- David G. Pisano
You can also search for this author in PubMed Google Scholar
- Alfonso Valencia
You can also search for this author in PubMed Google Scholar
- Roderic Guigó
You can also search for this author in PubMed Google Scholar
- Mónica Bayés
You can also search for this author in PubMed Google Scholar
- Simon Heath
You can also search for this author in PubMed Google Scholar
- Marta Gut
You can also search for this author in PubMed Google Scholar
- Peter Klatt
You can also search for this author in PubMed Google Scholar
- John Marshall
You can also search for this author in PubMed Google Scholar
- Keiran Raine
You can also search for this author in PubMed Google Scholar
- Lucy A. Stebbings
You can also search for this author in PubMed Google Scholar
- P. Andrew Futreal
You can also search for this author in PubMed Google Scholar
- Michael R. Stratton
You can also search for this author in PubMed Google Scholar
- Peter J. Campbell
You can also search for this author in PubMed Google Scholar
- Ivo Gut
You can also search for this author in PubMed Google Scholar
- Armando López-Guillermo
You can also search for this author in PubMed Google Scholar
- Xavier Estivill
You can also search for this author in PubMed Google Scholar
- Emili Montserrat
You can also search for this author in PubMed Google Scholar
- Carlos López-Otín
You can also search for this author in PubMed Google Scholar
- Elías Campo
You can also search for this author in PubMed Google Scholar
Contributions
X.S.P., V.Q., G.R.O., J.M.P.F., G.V., and A.G.-F. developed the bioinformatic algorithms for analysis of sequence data, implemented SMIPS for clinical validation and performed functional studies. M.P., L.C., P.J., M.J., M.L.-G., D.Colomer, S.G., C.L., A.N. and J.Y. were responsible for downstream validation analysis and functional studies. S.B., A.C., G.E., L.B., J.M.C.T., J.C.D. and C.T. studied structural variants. M.G., M.B., S.H., D.A.P., H.H., E.C., J.M., K.R. and L.A.S. were responsible for generating libraries, performing exome capture and running sequencers. P.N., C.M.R-C. and M.A. prepared and supervised the bioethical requirements. S.S. designed the epidemiologic study. P.K. gave conceptual advice and revised the manuscript. N.V., T.B., D.Costa, A.C., A.E., L.H., M.A.P., E.A., J.S.M., E.M., A.L.-G. and E.C. performed clinical and biological studies. M.R., N.V., J.M.H. and M.G.-D. were the pathologists who reviewed and confirmed the diagnoses. R.R., J.L.G., M.O., D.T., D.G.P., A.V. and R.G. were in charge of bioinformatics data management. I.G., P.A.F., M.R.S. and P.J.C. coordinated the sequencing efforts and performed primary data analysis. X.S.P., X.E., A.L-G., C.L.-O. and E.C. directed the research and wrote the manuscript, which all authors have approved.
Corresponding authors
Correspondence to Carlos López-Otín or Elías Campo.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
The file contains Supplementary Text, additional references, Supplementary Figures 1-10 with legends and Supplementary Tables 1-11 (excluding Supplementary Table 6). (PDF 4388 kb)
Supplementary Table 6
The table shows Somatic substitutions identified in CLL by WGS. (XLS 471 kb)
PowerPoint slides
Rights and permissions
This article is distributed under the terms of the Creative Commons Attribution-Non-Commercial-Share Alike licence (http://creativecommons.org/licenses/by-nc-sa/3.0/), which permits distribution, and reproduction in any medium, provided the original author and source are credited. This licence does not permit commercial exploitation, and derivative works must be licensed under the same or similar licence.
About this article
Cite this article
Puente, X., Pinyol, M., Quesada, V. et al. Whole-genome sequencing identifies recurrent mutations in chronic lymphocytic leukaemia. Nature 475, 101–105 (2011). https://doi.org/10.1038/nature10113
Received: 16 November 2010
Accepted: 06 April 2011
Published: 05 June 2011
Issue Date: 07 July 2011
DOI: https://doi.org/10.1038/nature10113
Editorial Summary
Leukaemia-linked gene mutations
Analysis of the genomes of four patients with chronic lymphocytic leukaemia, and validation in more than 300 patients, has identified four genes — NOTCH1, MYD88, XPO1 and KLHL6 — that are recurrently mutated in the condition. Mutations in NOTCH1, MYD88 and XPO1 are thought to contribute to the clinical evolution of the disease. Evidence that NOTCH1 and MYD88 mutations are activating events highlights them as potential therapeutic targets.