Whole-genome sequencing reveals activation-induced cytidine deaminase signatures during indolent chronic lymphocytic leukaemia evolution

Abstract

Patients with chromosome 13q deletion or normal cytogenetics represent the majority of chronic lymphocytic leukaemia (CLL) cases, yet have relatively few driver mutations. To better understand their genomic landscape, here we perform whole-genome sequencing on a cohort of patients enriched with these cytogenetic characteristics. Mutations in known CLL drivers are seen in only 33% of this cohort, and associated with normal cytogenetics and unmutated IGHV. The most commonly mutated gene in our cohort, IGLL5, shows a mutational pattern suggestive of activation-induced cytidine deaminase (AID) activity. Unsupervised analysis of mutational signatures demonstrates the activities of canonical AID (c-AID), leading to clustered mutations near active transcriptional start sites; non-canonical AID (nc-AID), leading to genome-wide non-clustered mutations, and an ageing signature responsible for most mutations. Using mutation clonality to infer time of onset, we find that while ageing and c-AID activities are ongoing, nc-AID-associated mutations likely occur earlier in tumour evolution.

ARTICLE

Received 2 Jul 2015 | Accepted 8 Oct 2015 | Published 7 Dec 2015

S. Kasar1,2,*, J. Kim3,*, R. Improgo1,2, G. Tiao3, P. Polak3, N. Haradhvala3, M.S. Lawrence3, A. Kiezun3,S.M. Fernandes1, S. Bahl3, C. Sougnez3, S. Gabriel3, E.S. Lander3, H.T. Kim4, G. Getz3,5,6,** & J.R. Brown1,2,**

Patients with chromosome 13q deletion or normal cytogenetics represent the majority of chronic lymphocytic leukaemia (CLL) cases, yet have relatively few driver mutations. To better understand their genomic landscape, here we perform whole-genome sequencing on a cohort of patients enriched with these cytogenetic characteristics. Mutations in known CLL drivers are seen in only 33% of this cohort, and associated with normal cytogenetics and unmutated IGHV. The most commonly mutated gene in our cohort, IGLL5, shows a mutational pattern suggestive of activation-induced cytidine deaminase (AID) activity. Unsupervised analysis of mutational signatures demonstrates the activities of canonical AID (c-AID), leading to clustered mutations near active transcriptional start sites; non-canonical AID (nc-AID), leading to genome-wide non-clustered mutations, and an ageing signature responsible for most mutations. Using mutation clonality to infer time of onset, we nd that while ageing and c-AID activities are ongoing, nc-AID-associated mutations likely occur earlier in tumour evolution.

1 Department of Medical Oncology, Dana-Farber Cancer Institute, Brigham and Womens Hospital, Harvard Medical School, Boston, Massachusetts 02215, USA. 2 Department of Medicine, Brigham and Womens Hospital, Harvard Medical School, Boston, Massachusetts 02215, USA. 3 Cancer Program, Broad Institute of Harvard and MIT, Cambridge, Massachusetts 02142, USA. 4 Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts 02215, USA. 5 Department of Pathology and Cancer Center, Massachusetts General Hospital, Boston, Massachusetts 02114, USA. 6 Department of Pathology, Harvard Medical School, Boston, Massachusetts 02215, USA. * These authors contributed equally to this work. ** These authors jointly supervised this work. Correspondence and requests for materials should be addressed to G.G. (email: mailto:[email protected]

Web End [email protected] ) or to J.R.B. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 6:8866 | DOI: 10.1038/ncomms9866 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

DOI: 10.1038/ncomms9866 OPEN

Whole-genome sequencing reveals activation-induced cytidine deaminase signatures during indolent chronic lymphocytic leukaemia evolution

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9866

Chronic lymphocytic leukaemia (CLL) is a clinically heterogeneous incurable malignancy of CD5 CD19 B lymphocytes1. Among the strongest predictors of

outcome are the disease-associated chromosome abnormalities, with 17p and 11q deletion and trisomy 12 associated with more aggressive disease, while 13q deletion (incidence 5060%) and normal cytogenetics (incidence 1520%) are lower risk according to Dohners cytogenetic classication2. Interestingly, the recurrent coding mutations identied to date in CLL have been associated with the higher-risk cytogenetic abnormalities, and are less commonly seen in CLLs with a lower-risk cytogenetic prole. We therefore set out to explore the genetic basis of the lower-risk cytogenetic group by whole-genome sequencing, as clues to the genetic basis of disease in this more indolent group may lie elsewhere in the genome.

Whole-genome sequencing provides unique information not available from prior studies with whole-exome sequencing, including data on translocations, complex rearrangements and genome-wide mutational patterns. However, relatively higher sequencing costs have limited the number of whole-genome studies (n 4, Puente et al.3; n 28, Alexandrov et al.4, with only

signature analysis reported without detailed cohort description) and to date, most studies involved larger exome data sets, which were likely the major driver of the primary ndings. Here we present a comprehensive analysis of structural rearrangements and somatic mutations in 30 CLL whole genomes having low-risk cytogenetic aberrations. We deliberately balanced our cohort to evenly represent higher- and lower-risk IGHV cases, since different driver events might be relevant to these subgroups, as in fact turned out to be the case.

Recently developed techniques using Non-negative Matrix Factorization (NMF)5 to perform unsupervised analysis of somatic mutation data has enabled the unbiased discovery of genome-wide mutational patterns in multiple tumour types4,6,7. One such study, by Alexandrov et al., analysed 28 CLL WGS and 103 whole-exome sequencing samples and found that CLL mutations comprise three mutational signatures: (i) ageing-related mutations (C4T at CpG mutations due to spontaneous deamination4; signature 1B); (ii) APOBEC signature (signature 2); and (iii) an activation-induced cytidine deaminase (AID)-related signature (signature 9). During B-cell development, AID induces deamination of cytosine to uracil. Resolution of these lesions by the error-prone DNA polymerase Z (eta) results in A to

C mutations at WA (W A or T) motifs; described as signature 9

by Alexandrov et al.4 However, as noted by Alexandrov et al.4, signature 9 does not exhibit the known mutation features of canonical AID (c-AID) (C to T/G at WRCY motifs, W A or T,

R purine, Y pyrimidine)8, and is therefore referred to as a

non-canonical AID (nc-AID) signature throughout this paper. Note that previously the c-AID signature was not separated as a distinct one in CLL4. However, prior experimental evidence has suggested that somatic hypermutation could be ongoing in a limited number of CLLs9. In addition, supervised mutation analysis8,10,11 did identify c-AID mutations in the immunoglobulin heavy chain locus in CLL8,9, as well as in multiple myeloma10 and diffuse large B-cell lymphoma (DLBCL)11. To the best of our knowledge, at the time of manuscript submission, genome-wide unsupervised discovery of c-AID signatures had not been performed in CLL.

In this study, our data set of 30 whole genomes provides the opportunity to perform an unsupervised analysis of the mutational patterns giving rise to indolent CLL. Here we present a modied Bayesian NMF algorithm that we have developed to analyse the mutation spectrum of CLL and show that it can successfully delineate both canonical and nc-AID signatures in an unsupervised genome-wide manner. In the

context of known CLL and AID biology, our results support a model of differential activities of the two AID signatures and the ageing signature throughout tumour evolution.

ResultsStructural rearrangements in CLL reveal chromoplexy events. To assess the degree of genomic structural instability in the 30 CLLs genomes (see Supplementary Tables 1 and 2 for patient characteristics), we rst analysed rearrangements, and identied a total of 92 events using the dRanger12 algorithm followed by BreakPointer12. This result corresponds to a median of 2.5 (range 015) rearrangements per genome, signicantly fewer than most other cancers (Fig. 1a), underscoring the relative stability of these CLL genomes. Interestingly, deletion of 13q often occurred by an inter-chromosomal unbalanced translocation (6/16, 37.5% of 13q deletion cases, Fig. 1a) rather than a simple deletion (see Supplementary Data 1 for rearrangement partners). Hruba et al.13 reported a similar frequency of inter-chromosomal 13q rearrangements by uorescence in situ hybridization. Apart from chr13, three other chromosomes were rearranged in X20% cases;

chr2 (12 cases, 40%), chr14 (9 cases, 30%) and chr1 (6 cases, 20%). Chr14 accounted for 14% of rearrangements (13/92), mostly representing deletions. Eleven of these 13 rearrangements (9 cases) had break points at the 50IgH region (chr14q32.33). Patients with chr14q32.33 rearrangements had a shorter time to next treatment post sampling (TTNT) (P 0.019, log-rank test, Fig. 1b). Although

deletions in this locus have been frequently reported in a variety of B-cell neoplasms including CLL, they have not been detected in normal B cells, indicating that they are not a by-product of normal immunoglobulin rearrangements14,15. Single-nucleotide polymorphism (SNP) array analysis of the 30 CLLs revealed a median of 1 somatic copy-number alteration (sCNA) per case (range 06), similar to our previous report16. Other than 13q loss, we detected previously described sCNAs such as focal amplications in 3q25.33 (ref. 15) in two cases (which include PIK3CA) and a focal deletion at 1q42.2 that was reported by Pfeifer et al.17

Next, we looked for complex structural rearrangements such as chromothripsis18 and chromoplexy19 since these may have disrupted multiple genes in a single event. While chromothripsis typically involves multiple focal deletions in a single chromosome (thought to occur during metaphase)20,21, chromoplexy is dened as a series of inter-dependent rearrangements among multiple chromosomes (most likely during interphase)19. Although we did not nd evidence of chromothripsis, three of the 30 cases had evidence for chromoplexy (detected by ChainFinder19an algorithm that links close rearrangements to balanced chains of events). Two of the cases had a single chain (both with three rearrangements) and one had two chains (with three and eight rearrangements) (Fig. 1c; Supplementary Data 2). In one case each, the chain included known common CLL copy-number changes, namely 13q deletion and 14q32 deletion. Interestingly, all three of these patients were untreated before sampling but underwent therapy shortly thereafter, suggesting that these events may indicate poor outcome (P 0.02, log-rank test),

although this nding needs to be conrmed in a larger cohort. Our results indicate, for the rst time, that chromoplexy events occur in CLL (involving 17 of the total 92 rearrangements). Earlier cytogenetic reports of chained translocations may also have reected this phenomenon, albeit at much lower resolution. Statistical analysis of the copy-number and structural rearrangement data by the ChainFinder algorithm suggests that these events likely occurred at the same time and hence adds additional information beyond the previous cytogenetic studies. Taken together, our data suggest that in this indolent cohort a subset of 13q deletions may occur by inter-chromosomal

2 NATURE COMMUNICATIONS | 6:8866 | DOI: 10.1038/ncomms9866 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9866 ARTICLE

a b

Chr14q32.33 rearrangements

Present

Absent

P=0.019

x

y

1

22

21

20

19 18

17

16

15

14

13

12

11

10

2

4

1.0

0.8

0.6

0.4

0.2

0.0

0 2 4 6 8 Years post sampling

3

Probability

5

6

8

7

6

c

1

1

22

21

20

19 18

17

16

15

14

13

12

11

10

x

y

22

21

20

19 18

17

16

15

14

13

12

11

10

x

y

x

y

1

2

2

22

21

20

19 18

17

16

15

14

13

12

11

10

2

3

3

3

4

4

4

6

8

7

6

5

6

8

7

6

5

6

8

7

6

5

Figure 1 | Summary of structural rearrangements. (a) Circos plot representing the structural rearrangements observed across 30 CLL genomes. Purple lines indicate inter-chromosomal rearrangements, grey lines indicate intra-chromosomal rearrangements; red arrows point to inter-chromosomal rearrangements giving rise to 13q deletion. (b) KaplanMeier curve showing the relationship between time to next treatment post sampling and rearrangements in Chr14q32.33 (50IGH) in the vicinity of KIAA0125. (c) Circos plots depicting the presence of chained rearrangements detected by the

ChainFinder algorithm. Red arrows indicate deletion bridges and inter-dependent chains. Left1 chain near LPIN1, TRIB2 and TMEM194B geneson chr2; middle1 chain near KIAA0125 on chr14 and ANP32A on chr15; rightchain 1 (blue) near SNAR-H, REG3G and CTNNA2 on chr2and chain 2 (green) near ARMCX6, SAGE1, ZCCHC5 and ITM2A on chr23 and CYSLTR2, EBPL, RNASEH2B and KPNA3 on chr13. The genes listed here either fall within a deletion or are within 25kb of a chained rearrangement breakpoint.

rearrangement or even more involved chromoplexy events. Future larger studies are needed to correlate such events more denitively with clinical outcome.

Increased subclonal mutation rate with age. Turning our attention to mutational patterns, we identied an average of 3,055 mutations per genome (Fig. 2a; Methods). The average genome-wide mutation frequency of 1.10.4 per Mb (range 0.42.1, n 30, data are shown as means.d.) is lower than that of many

other haematological malignancies and solid tumours7,21,22. The pattern of mutation densities (intergenic4intronic4untranslated region (UTR)4exonic) in different genomic regions was similar to other WGS studies23 (Fig. 2b).

Although CLL is largely considered to be a geriatric malignancy, about one-third of patients develop the disease much earlier. We have previously reported that disease diagnosed at an older age is associated with a higher number of clonal mutations (in coding regions24), but not subclonal mutations. Clonal mutations are likely to have occurred during life before onset of malignancy, while subclonal mutations arise subsequent to transformation, after the last selective sweep, and are therefore in only a subset of cells24. In this cohort, we deliberately selected older and younger patients matched for other disease characteristics, so as to better associate mutational patterns with age of diagnosis. We conrmed the expected increase in clonal mutations with age, but we were also able to observe a clearly higher number of subclonal mutations with older age (Fig. 2cg), even though the disease duration was comparable in the young and old cohorts (Fig. 2d). Thus, in addition to the well-described higher clonal mutation burden acquired before disease development in older patients, we also nd

a higher ongoing rate of subclonal mutations, which may reect more clonal evolution and heterogeneity.

Somatic mutational landscape of indolent CLL by WGS. Focusing on the specic somatic mutations, we observe that only 10 (33%) out of 30 patients displayed at least one mutation in a previously reported CLL driver. In comparison, 57% (91/160) of the cohort in our previous whole-exome study24 had at least one mutation in a CLL driver, indicating that this cohort does indeed capture a different biology (Fishers exact test, P 0.027). Those

patients with at least one driver mutation in a previously reported CLL cancer gene were more likely to have unmutated IGHV (P 0.014) and normal cytogenetics as compared with 13q

deletion (P 0.033, Fishers exact test). These patients also had a

higher risk of progressing to next treatment (Hazard Ratio, HR 5.71, P 0.0076), as expected by their IGHV status.

Interestingly, 7/30 cases (23%) did not harbour mutations in any gene previously associated with cancer (CLL drivers24, COSMIC or PanCancer25 mutations) (Fig. 3a); these patients all carried 13q deletion (Fishers exact test, P 0.009) and ve of them had

mutated IGHV. The number of nonsilent mutations per tumour in these seven cases was also signicantly lower than the rest of the cohort (135.5 versus 208.6, P 0.0027, MannWhitney U-

test, data are shown as means.d.). No difference in the number of rearrangements and sCNAs was seen in these two groups.

Interestingly, the 13q-deleted subgroup was enriched in 50UTR and coding-region mutations in IGLL5 (P 0.04, Fishers exact

test), the gene carrying the most frequent coding-region mutations in our cohort (Fig. 3b,c). These mutations were also more common in IGHV-mutated cases (P 0.013).

NATURE COMMUNICATIONS | 6:8866 | DOI: 10.1038/ncomms9866 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9866

a

b

c

d

Mutation density per MB

1.5

1.0

0.5

0.0

Age at diagnosis (years)

100

80 60 40 20

0

Time from Dx to sampling (years)

20

15

10

5

0

P=NS

Total=91,672

0.60% 3UTR0.20% 5UTR0.00% De_novo_Start_InFrame0.00% De_novo_Start_OutOfFrame57.27% IGR31.25% Intron6.22% lincRNA0.54% Missense_Mutation0.03% Nonsense_Mutation0.00% Nonstop_Mutation3.64% RNA0.21% Silent0.05% Splice_Site0.00% Start_Codon_SNP

IGR

Intron

Noncoding

UTR

Exon

Old Young

Old Young

e

f

g

Total

Coding

All

Old Young

P=0.0032

P=0.0028

No. of mutations per sample

40

30

20

10

P=0.025

P=0.0001

Mutation rate per MB

1.5

1.0

0.5

0.0

No. of mutations per sample

4,000

3,000

2,000

1,000

0

P=0.0006

P=0.019

Total Coding

Old Young Old Young

0 Old Young Old Young

Clonal Subclonal

Clonal Subclonal

Figure 2 | Overview of somatic mutational landscape. (a) Pie chart depicting the percentage of different types of sSNVs detected in our cohort genome-wide. (b) Bar chart of average mutation densities across different regions of the genome. n 30, error bars indicate s.e.m. (c) Dot plot of

age at diagnosis in the older versus younger cohort. The horizontal line indicates median age. (d) Dot plot of time from diagnosis to sampling in the older versus younger cohort. (e) Bar chart comparing the mutation rate per MB genome-wide (total) and in coding regions in the entire cohort and in younger (n 13) versus older (n 17) patients. (f,g) Dot plot of average number of clonal and subclonal mutations total (f) and in coding regions

(g) in younger versus older subgroups is shown. Error bars indicate s.e.m., P values were calculated using the MannWhitney U-test. NS, not signicant(i.e. P40.05).

Little is known about the function of the IGLL5 gene, but it is homologous to IGLL1 (lambda5), which is critical for B-cell development. Furthermore, IGLL5 has been reported to be recurrently mutated in diffuse large B-cell lymphoma26. The mutation pattern in IGLL5 was suggestive of off-target AID activity, with clustering of mutations near the transcription start site (TSS) through the rst intron, as well as biallelic mutations (Supplementary Fig. 1a). The mutations included non-synonymous coding-region mutations (n 4), 50UTR mutations

(n 4) and one patient with both (total 9/30 30%), as well as 15

samples with mutations in the rst intron (total 15/30 50%).

The 50UTR and coding mutations in IGLL5 were enriched in subclonal mutations, whereas the intronic mutations were mostly clonal (P 0.006, Fishers exact test), suggesting that the 50UTR/

coding mutations were acquired later than the intronic mutations, after the last selective sweep. The presence of 50UTR and rst exon mutations was conrmed by Sanger sequencing (n 7/8 cases, we did not have additional DNA from the 8th

patient after sequencing). In addition, expression of the mutant alleles in the IGLL5 coding region was also conrmed using matched RNA sequencing (RNA-seq) data (Supplementary Fig. 1b shows a representative Integrative Genomics Viewer [IGV] screenshot). IGLL5 mutants showed a trend towards reduced transcript levels as compared with wild type (Supplementary Fig. 1c). Comparing the fraction of reads supporting the mutated allele in the WGS and RNA-seq data showed higher mutation allele fraction in the RNA-seq data in the coding mutations (P 0.0078), whereas the 50UTR mutations had

similar allele fractions (P 0.16) (Supplementary Fig. 1d). These

data suggest a potentially different functional role for coding and 50UTR mutations in IGLL5, but future experiments will be required to determine their true role if any in CLL pathogenesis. Given that IGLL5 was the most commonly mutated gene in our cohort and the mutational pattern suggested the potential involvement of AID activity, we were interested in exploring

more broadly the mutagenic processes, including AID, that give rise to the somatic mutations that lead to CLL27,28.

Unsupervised discovery of mutational signatures. Normal B cells undergo somatic hypermutation in the germinal centre29a process that is mediated by AID and induces clustered mutations in immunoglobulin loci and some off-target regions10,30. Following AID-induced deamination of cytosine to uracil, different repair processes lead to different mutational signatures, called either c-AID or nc-AID. Specically, direct replication over the AID-induced U:G lesions or removal of the uracil by UNG (uracil DNA glycosylase) followed by replication accounts for the mutations of the c-AID signature (C to T/G mutation at WRCY motifs, W A or T,

R purine, Y pyrimidine; reviewed in ref. 31). Alternatively,

processing of the AID-induced lesions by the mismatch repair pathway that recruits the error-prone DNA polymerase Z gives rise to the nc-AID-related mutations (A to C/G at WA motifs28,32; reviewed in refs 31,33).

Given our ndings with IGLL5, and due to the known clustered nature of c-AID mutations34, we considered the nearest mutation distance (NMDsee Methods for details) as a parameter to stratify somatic mutations. We observed a bimodal distribution of mutation distance that enabled a partitioning of the mutations into two groups: (i) a clustered group (NMDo1,000 nt)

consisting of 7% of mutations and (ii) a non-clustered group (NMD41,000 nt) with the remaining 93% of mutations (Fig. 4a).

Comparing the mutational spectra of these two groups revealed a marked increase of C4T/G at GCT motifs in the group of clustered mutations. This pattern of mutations matches the known c-AID signature7, suggesting that this process contributes to the mutational load in CLL (Supplementary Fig. 2). Although AID expression has been reported in only 0.012% of quiescent circulating CLL cells35, our nding is consistent with previous

4 NATURE COMMUNICATIONS | 6:8866 | DOI: 10.1038/ncomms9866 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9866 ARTICLE

a

Cytogenetics

JB-0025

JB-0006

JB-0073

JB-0007

JB-0010

JB-0011

JB-0005

JB-0012

JB-0031

Zap70

IGHV

Age_at_Dx

b

13q del Normal Other

Neg Pos unk

unk

Mut

unmut

Old Young

CW176

CW114

JB-0014

JB-0004

CW127

CW126

JB-0017

JB-0030

JB-0032

JB-0018

JB-0026

JB-0009

JB-0008

JB-0003

JB-0016

JB-0024

JB-0029

JB-0023

JB-0013

JB-0072

CW119

cdsStart +274

+170

TSS

cdsEnd +7,915

% Incidence 20 10

Cytogenetics Zap70IGHV Age_at_Dx

+429

+6,039, splice

***

ATM BCOR DDX3X EGR2 FBXW7 KRAS MED12

NOTCH1 SF3B1

TP53 XPO1 ZMYM3

BCL11A BCL11B BIRC3 BRAF CBL CLTCL1 CNOT3 DICER1 EBF1

FUBP1 HERPUD1 KDM6A

MET

NIN

NONO PAX5 PTPRC TCL1A TERT FAT1

MGA MXRA5 OR4A16 POU2F2

RPS15 SACS SMC1A XIRP2

IGLL5 nonsilent IGLL5 5UTR

Rearrangement >0.2 MB CNA

Driver Driver Driver Driver Driver Driver Driver Driver Driver Driver Driver Driver

COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC COSMIC

PanCan PanCan PanCan PanCan PanCan PanCan PanCan PanCan PanCan

23,230,439

23,229,959 23,235,879 23,235,998 23,238,005

c

L28L A32T A44S G110Splice

209

214

35

Signal

G19S L38QA44V

P52R

J

109

115

Ig-like-C1-type

1

98

Figure 3 | Distribution of mutations in selected genes. (a) Heatmap showing the presence of non-synonymous mutations in genes specied onthe right. In the heatmap, white box = one event, black box two events and gray no events. The genes are classied based on potential functional

signicance as shown in the rightmost column. The top panel shows the clinical characteristics of each sample. The bar chart on the left indicatesthe percentage of cases with at least one mutation in the gene on the right. The bottom four rows in the heatmap represent the presence of mutations in IGLL5, rearrangement events and copy number alterations. The black box highlights samples with mutations in known CLL driver genes; the bluebox highlights cases with no mutations in known cancer-associated genes. (b) Graphical representation of 50UTR and coding mutations in the

IGLL5 transcript. ***indicates mutations concentrated in the rst intron. (c) Graphical representation of IGLL5-coding mutation alterations at the protein level.

studies that have reported AID activity in CLL cells by analysing intra-clonal IGHV diversity36,37 and induction of de novo somatic hypermutation in vitro8,38.

Recently, Alexandrov et al.4 identied, in an unsupervised manner, 21 mutational signatures across 30 different tumour types by applying NMF to the mutation counts across the 96 available trinucleotide mutation contexts. Here we characterized the mutational signatures operating in our 30 CLL cases using a related Bayesian NMF method5 considering NMD as an additional feature. Thus, instead of analysing a 96-by-30 matrix of mutation counts, we partitioned the mutations in each tumour into two groups of clustered and non-clustered mutations, giving rise to a 96-by-60 matrix (Methods). This partitioning enabled the discovery of mutational signatures unique to the clustered and non-clustered mutations. Our analysis identied three mutational processes, only two of which were reported by Alexandrov et al.4: an ageing signature characterized by increased C4T transitions at CpG sites (analogous to their signature 1B); a nc-AID signature, dominated by A4C at WA motifs (analogous to their signature 9) (Figs 4b,c and 5a); and a third signature that matches the c-AID signature (C to T/G mutation at WRCY motifs, W A or T, R purine (A or G), Y pyrimidine (C or

T)) that was not reported by Alexandrov et al.4 To further validate the nding of the c-AID signature, we reanalysed the 28 WGS CLL samples from Alexandrov et al. using our method and were able to validate both the c-AID and nc-AID signatures in their data (Supplementary Fig. 3), although the c-AID signal was not as strong as in our cohort. Thus, our analysis provides denitive evidence that c-AID activity in CLL is strong enough to be discovered in an unsupervised analysis of genome-wide mutational patterns.

Next, we calculated for each mutation, m, the probability (pms)

that it was generated by each of the three mutational signatures, s, and assigned it to a signature if that probability (pms) was greater

than 0.75 (Methods). As expected, plotting the NMD along the genome for each signature revealed that the c-AID-associated mutations form distinct clusters, whereas the nc-AID- and ageing-associated mutations are scattered more evenly (Supplementary Fig. 4). From this analysis, we were able to determine that the three signatures exhibit differential contribution to the overall mutational landscape of each patient. The ageing signature was predominant across all cases and the number of ageing-related mutations was signicantly higher in patients with older age at diagnosis, as might be expected (Fig. 5b, P 0.004, Wilcoxons rank-sum test). However, 70% of cases had

at least 10% of mutations due to AID activities (Fig. 5a). The number of mutations due to c-AID and nc-AID was signicantly higher in IGHV-mutated CLLs (Fig. 4c and Fig. 5c,d; c-AID P 0.0004, nc-AID Po0.0001, Wilcoxons rank-sum test).

Consistent with this, 7/12 IGHV-unmutated cases showed 495% ageing signature. Among coding mutations, 95% were associated with the ageing signature (pm,ageing40.75), whereas

only 1.9 and 2.4% were associated with c-AID and nc-AID, respectively. Therefore, the ageing signature is likely to be the primary contributor to driver mutations in coding regions in CLL. Interestingly, the seven samples with no mutations in known CLL drivers or other cancer genes showed a lower number of ageing-associated mutations (P 0.021, MannWhitney test).

c-AID signature exhibits classical features of SHM. Apart from the WRCY recognition motif, other previously described

NATURE COMMUNICATIONS | 6:8866 | DOI: 10.1038/ncomms9866 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9866

a

Distribution of NMDs in CLL 30 WGS

Frequency

2,500

1,500

500

0

1 2 3 4 5 6 7 Log10 (nearest mutation distance)

b

Mutation contributions

3,000 2,000 1,000

0 200 150 100

500

1,000

500

0

Tri-nucleotide motifs

c

sSNVs associated with signatures discovered in 30 CLL WGS

Samples

NMD 1,000 nt NMD >1,000 nt

Mutation contributions

Mutation signatures discovered in 30 CLL WGSA->C A->G A->T C->A C->G C->T

Ageingc-AIDnc-AID

4,000

3,000

2,000

1,000

0

300 200 100

0 2,000

1,500

1,000

500

0

0003

0004

Ageingc-AIDnc-AID

0005

0006

0007*

0008*

0009*

0010*

0013*

0014*

0016*

0017*

0018*

0023*

0024*

0011

0012

0025

0026

0029*

0030*

0031

0032

0072

0073

CW127*

CW126*

CW119*

CW114*

CW176

0003

0004

0005

0006

0010*

0009*

0008*

0007*

0011

0012

0024*

0023*

0018*

0017*

0016*

0014*

0013*

0025

0026

0031

0032

0072

0073

0029*

0030*

CW176

CW127*

CW126*

CW119*

CW114*

Figure 4 | Analysis of mutational signature in CLL. (a) Frequency histogram of nearest mutation distance (NMD) shows bimodal distribution. (b) Estimated mutation contributions of the indicated mutational signatures detected upon inclusion of NMD as a factor in Bayesian NMF. (c) Number of clustered mutations (left) and non-clustered mutations (right) associated with canonical AID (green), ageing (blue) and non-canonical AID (purple) signature across samples. * Indicates cases with mutated IGHV.

characteristics of off-target c-AID activity include (1) multiple and biallelic mutations (evident from the observed clustering and the mutation pattern in c-AID targets such as IGLL5) and (2) preferential targeting 2 kb from the TSS of highly transcribed genes (reviewed in ref. 31). We observe that the c-AID-associated mutation rate was increased 2.5-fold within 2 kb of TSS as compared with the genome-wide rate (Fig. 6a). To conrm the preference for highly transcribed genes, we divided the genes into four quartiles based on their expression levels determined by RNA-seq, and compared the contribution of the three signatures with the mutation rate in genes in each quartile. As expected, the overall rate of mutations decreases with higher expression levels due to transcription-coupled repair. However, the c-AID mutation rate was found to be the highest compared with the other two signatures in the genes in the quartile with the highest expression (Q4) (Fig. 6b).

Identication of genome-wide targets of the AID signatures. Next, we focused on the contribution of c-AID and nc-AID to the mutational density in individual genes (including UTRs and introns), to nd specic target genes unique to each of the AID processes. Specically, we identied non-overlapping sets of trinucleotide sequence contexts that distinguish the c-AID and nc-AID signatures (Supplementary Fig. 5a) (using only mutations with pms40.75). For each AID signature, we then compared the observed mutation density in every gene with at least one signature-associated mutation (pms40.75, 281 c-AID and 809 nc-AID genes) to the context-specic background mutation density in these genes (Methods, Supplementary Fig. 5). We then corrected for multiple hypotheses and identied genes associated with each signature using a q-value cutoff of 0.1 (Supplementary Fig. 5b).

6 NATURE COMMUNICATIONS | 6:8866 | DOI: 10.1038/ncomms9866 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9866 ARTICLE

a b

100%

90% 80% 70% 60% 50% 40% 30% 20% 10%

0%

Ageing

P=0.004

nc-AID

Ageing

c-AID

Mutation count

Mutation count

4,000

3,000

2,000

1,000

0

5,000

Normalized signature activity

IGHV Mut

Old

IGHV Unmut

Young

JB-0025

JB-0006

JB-0073

JB-0007

JB-0010

CW176

JB-0011

JB-0005

JB-0012

JB-0031

JB-0014

CW114

CW127

JB-0004

CW126

JB-0017

JB-0030

JB-0032

JB-0018

JB-0026

JB-0009

JB-0008

JB-0003

JB-0016

JB-0024

JB-0029

JB-0023

JB-0013

CW119

JB-0072

c d

Non-canonical AID 3,000

2,000

1,000

0

Canonical AID

Mutation count

400

300

200

100

0

P<0.0001

P=0.0004

IGHV Mut

IGHV Unmut

Old

IGHV Mut

IGHV Unmut

Young

Old

Young

Figure 5 | Association of signatures with clinical characteristics. (a) Percentage contribution of each of the mutational signatures to the overall mutation spectrum across samples. (bd) Dot plots showing total mutation counts associated with the indicated signatures in relation to age at diagnosis (younger versus older) and IGHV mutation status (mut versus unmut). Error bars indicates.e.m. P values were calculated using the Wilcoxons Rank Sum Test.

1.5

1

0.5

0

0.5

1

1.5

a b

c-AID Ageing nc-AID

Log 2ratio

Ratio of mutation density within 2kb of TSS to whole-genome density 1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0

Normalized mutation frequency

c-AID Ageing nc-AID Increasing expression

Q4

Q1

Q2

Q3

Figure 6 | c-AID mutations exhibit classical features of SHM. (a) Ratio of mutation frequency within 2 kb of transcription start site (TSS) to the genome-wide mutation rate for each signature. (b) Genes were divided into four quartiles, Q1 through Q4, in order of increasing expression. Bar graph showing normalized mutation density of each signature per quartile.

For c-AID, we detected 34 associated genes with qo0.1 (Supplementary Data 3). Consistent with known AID biology, 24 (70%) of these genes were located in the cytobands with the three immunoglobulin loci (14q32.33, 22q11.22 and 2p11.2). Unsurprisingly, IGLL5 was one of the most signicant genes associated with the c-AID signature in this analysis (qo10 90,

Supplementary Data 3). The list also included BCL6 and LTB, other known off-targets of AID in post-germinal centre B-cell malignancies39,40 and another haematologic malignancy-related gene, TRIP11 (ref. 41) (Supplementary Data 3).

For nc-AID, we discovered 14 genes that were specically targeted (qo0.1, Supplementary Data 4). This list includes two genes in the immunoglobulin cytobands as well as cancer-related

genes such as CADM2 (renal cell carcinoma)42, CHRM3 (colon cancer)43, LPHN3 (panCancer analysis)44 and ROBO1 (breast cancer)45 (Supplementary Data 4). The biologic basis of this signature selectivity and its relevance to cancer development will need to be claried in future studies.

Ageing and c-AID activities are ongoing in CLL. We analysed mutation clonality to assess the activity of these three mutational processes over time during the life history of the CLLs. The clonality of a given mutation can be used to deduce the time of its onset in relation to the most recent selective sweep, with clonal mutations being earlier events and subclonal mutations occurring

NATURE COMMUNICATIONS | 6:8866 | DOI: 10.1038/ncomms9866 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9866

a

Signature association > 0.75 Signature association > 0.75

*

30,000

Clonal

Subclonal

Mutation counts

Fractions

1.0

0.8

0.6

0.4

0.2

0.0

20,000

10,000

0 c-AID Ageing nc-AID Other

c-AID Ageing Total

c-AID Ageing nc-AID Other

Signatures Signatures

b

c-AID > 0.75

c-AID > 0.75

Ageing > 0.75

nc-AID > 0.75

nc-AID > 0.75

Clonal Subclonal

Clonal Subclonal

Clonal Subclonal

600

400

200

Counts

Counts

15,000

10,000

5,000

Fraction

Counts

3,000

2,000

1,000

0 0.025

0.22

0.38

0.52

0.68

0.82

0.98

0 0.025

0.22

0.38

0.52

0.68

0.82

0.98

0 0.025

0.22

0.38

0.52

0.68

0.82

0.98

CCF

CCF

CCF

Ageing > 0.75

Fraction

0.20

0.15

0.10

0.05

0.00

0.8

0.6

0.4

0.2

0.0

Fraction

0.15

0.10

0.05

0.00

0.025

0.22

0.38

0.52

0.68

0.82

0.98

0.025

0.22

0.38

0.52

0.68

0.82

0.98

0.025

0.22

0.38

0.52

0.68

0.82

0.98

CCF

CCF

CCF

c

IGHV mutated IGHV unmutated

P=0.001 NS

P=0.04

5

4

3

2

1

0

Subclonal:clonal mutations

Figure 7 | Chronological order of mutational processes. (a) Bar graph showing absolute number (left) and ratio (right) of clonal and subclonal mutations in the indicated categories (pms40.75). Other includes mutations that were not assigned to any of the three signatures. *Po0.000001, P value was calculated using the Chi-square Test. (b) Distribution of CCF of mutations assigned to each signature; total number (top) and the fraction (bottom) of mutations for given CCF. (c) Ratio of subclonal:clonal mutations among mutations associated with either c-AID or ageing, compared with total mutations, shown divided by IGHV status. Note that pms40.5 is shown here (P 0.001) since pms40.75 had a low overall n, but a similar trend was observed with

pms40.75, P 0.055. We only considered cases with at least ve c-AID-associated mutations, resulting in N(IGHV

mut)

17 and N(IGHV unmut) 8 for

pms40.5. NS, not signicant. (i.e. Po0.05). Error bars indicates.e.m. P values were calculated using the Mann Whitney U Test.

later. We used ABSOLUTE46 to assess clonal versus subclonal status (Methods) and examined the proportion of clonal and subclonal mutations associated with each mutational process (pms40.75) (Fig. 7a). Overall, we found a signicant association between clonality and the mutational processes (Po0.00001, w2-test). Next, we tested each mutational process independently and found that each signature had a different proportion of clonal mutations. The nc-AID signature had the highest proportion of clonal mutations (Po0.00001, Fishers exact test).

On the other hand, c-AID-associated mutations were equally distributed between clonal and subclonal populations (P 0.26,

Fishers exact test), and ageing was enriched in subclonal mutations (Po0.00001, Fishers exact test).

Given this difference in the proportion of clonal mutations in each signature, we were interested in using these data to infer the time of onset of each of the mutational processes in the life history of the CLL. To do this, we looked at the distribution of the proportion of tumour cells bearing a signature-associated mutation, namely, the cancer cell fraction (CCF) for each mutation. We plotted the fraction of mutations associated with each signature as a function of their CCF. This analysis showed that 54% of the nc-AID mutations were clonal (high CCF). As

8 NATURE COMMUNICATIONS | 6:8866 | DOI: 10.1038/ncomms9866 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9866 ARTICLE

expected, therefore, the fraction of nc-AID-associated mutations declined sharply at low CCF values, indicating relatively few subclonal mutations, and suggesting that nc-AID was more active at earlier stages of tumour evolution (Fig. 7b; Supplementary Fig. 6). In contrast, roughly 40% of the c-AID mutations were clonal, and c-AID mutations showed a constant proportion across CCF values, suggesting both early and ongoing c-AID activity. Ageing-associated mutations were least likely to be clonal, with only 36% clonal, and, consistent with that, ageing-associated mutations showed a slight increase towards lower CCF, that is, more representation in subclonal mutations. These data therefore suggest that nc-AID occurred mostly at earlier times, whereas c-AID and ageing are continuing to operate even after the last selective sweep.

Ongoing c-AID activity is enriched in unmutated IGHV cases. Being a strong mutator, AID activity is tightly regulated in cells47. It is expressed in a very small fraction of circulating CLL cells, likely those in the proliferative fraction. Interestingly, and unexpectedly given that more somatic hypermutation is present in mutated IGHV CLLs, this expression of AID (encoded by the AICDA gene) among circulating CLL cells has been shown to be enriched among unmutated IGHV cases48,49. In our present cohort, although AID mRNA expression is very low overall, its expression is signicantly higher in unmutated IGHV (Supplementary Fig. 7, P 0.001, MannWhitney U-test).

Hence, we hypothesized that the ongoing c-AID activity evident in the low-CCF subclonal c-AID-associated mutations would be enriched in IGHV-unmutated patients. In fact, the ratio of c-AID-associated subclonal:clonal mutations was higher in the unmutated compared with mutated IGHV CLLs (Fig. 7c, pms40.5, P 0.001;

pms40.75, P 0.055, MannWhitney U-test, Supplementary Data

5 and 6), even though the overall frequency of c-AID-associated mutations was higher in mutated IGHV CLLs (Fig. 5d). These data suggest that the ongoing c-AID activity is enriched in unmutated IGHV CLLs, even though the sum total of c-AID activity across all of tumour evolution is enriched in mutated IGHV CLLs. As a control, we assessed whether the ratio of ageing-associated subclonal:clonal mutations was associated with IGHV status, and found no association (Fig. 7c).

DiscussionIn summary, we describe here the results of whole-genome sequencing of a CLL cohort comprised of low-risk cytogenetic subgroups in which we nd that a small subset have complex rearrangements that may be associated with more aggressive disease, while a signicant number have only 13q deletion as an obvious CLL driver.

Unlike previous studies24, we nd signicant enrichment in not just clonal but also subclonal mutations with age. Acquisition of subclonal mutations or clonal evolution has been associated with worsening disease24,50. While this supports the paradigm of acquisition of passenger mutations with age, it also points towards a more heterogeneous tumour in older patients and/or a faster ongoing acquisition of new mutations within the tumour. Thus, the age-associated increase in clonal diversication may be a key factor promoting worse disease outcomes in older patients.

We discovered recurrent mutations in IGLL5 that were previously undescribed in CLL. Interestingly, these mutations segregate independently of the known CLL driver genes and thus seem to be a unique feature of low-risk CLL. The pattern of IGLL5 mutation is suggestive of off-target AID activity, which is more prominent in lower-risk IGHV-mutated CLL. The mutations are expressed and were associated with a trend towards lower overall gene expression. Although the complete

functional characterization of this protein is beyond the scope of this manuscript, we have presented several indicators that suggest IGLL5 mutations may be of biological importance. Taken together, these ndings point towards a potential functional role of IGLL5 perturbation in low-risk CLL. However, further experimental work is required to conrm any such role.

Systematic analysis of mutational signatures gives insights into key mutagenic processes governing the developmental history of a cancer cell. Using a novel signature discovery method that uses information on both sequence context and mutation distance, we were able to identify three mutational signatures operative in CLL, including two distinct AID processes (c-AID and nc-AID) that represent a greater fraction of mutational activity in mutated IGHV cases, and the ageing-related signature. Somatic hypermutation is a critical physiological process in B-cell development responsible for afnity maturation of antibodies51. This process is initiated by AID, a 24-kDa protein that catalyses cytosine deamination to produce uracil, thereby creating U:G mismatches52. Repair of AID-induced lesions give rise to C to T/G mutations at WRCY motifstermed as c-AID, and A to C/G mutations at WA motifstermed as nc-AID51. Although AID activity is tightly regulated to primarily target the immunoglobulin-variable region genes, off-target AID activity can cause oncogenic mutations and chromosomal instability51. Despite scattered data suggesting that c-AID activity is present in CLL, the c-AID mutational signature was not identied in a previous unsupervised analysis of mutational signatures for CLL4. Here we demonstrate that the c-AID activity is separable as a distinct mutational signature using genome-wide unsupervised analysis considering the mutation distance as an additional feature.

Recently, Pettersen et al. reported c-AID-induced mutations in kataegis regions in CLL using a supervised motif discovery method7. A similar supervised motif discovery in four CLL whole genomes by Rebhandl et al. had previously implicated APOBEC activity in CLL53. Although APOBEC is widely expressed in CLL, our unsupervised signature discovery did not yield an APOBEC mutational footprint in our cohort. We applied the supervised motif analysis, as per Rebhandl et al.53, to our cohort, and were unable to detect evidence for APOBEC activity in clustered and non-clustered mutations at either immunoglobulin or non-immunoglobulin loci.

Analysis of whether a mutation is clonal or subclonal can be used to infer the time of occurrence of that mutation in relation to initial malignant transformation, with clonal mutations occurring earlier. We therefore examined the clonal fraction of mutations associated with each of our signatures. The ageing signature activity was enriched in patients with late-onset disease and enriched in subclonal mutations, which occur later, therefore suggesting that the ageing signature is a source of ongoing mutagenesis in CLL. This nding is consistent with our observation of not just increased clonal, but also increased subclonal mutations with age. Interestingly, the nc-AID-associated mutations were more clonal, suggesting that this process primarily occurred before the last selective sweep and perhaps even before cancer initiation. The mutation clonality analysis suggested that c-AID activity represents both an early and an ongoing process in the CLL life cycle. The higher proportion of newer subclonal c-AID-related mutations in IGHV-unmutated CLL, suggests that ongoing c-AID activity is higher in this subgroup. These ndings are consistent with prior work showing that in mature circulating CLL cells, AID activity is more easily induced in unmutated IGHV patients, and hence more likely to create newer mutations8. It should be noted that the sum total of all AID-related mutations is signicantly higher in the mutated IGHV cases, as also reported by Alexandrov et al.4

NATURE COMMUNICATIONS | 6:8866 | DOI: 10.1038/ncomms9866 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9866

Circulating CLL cells are mostly in the G0/G1 phase54,55, although recent work has demonstrated that a small pool of proliferating cells is always present56, likely arising from tissue niches. Our data suggest that DNA-polymerase-Z-mediated repair of AID-induced genetic lesions, which results in nc-AID mutations, occurs predominantly earlier in the CLL life cycle, perhaps even before transformation, while UNG-mediated repair, which results in c-AID mutations, is ongoing early and later. Little is known about the factors governing the relative activity and timing of pol-Z and UNG in repairing AID-mediated strand breaks. DNA pol-Z which contributes to nc-AID mutations, is more active in the S phase57,58. UNG, which contributes to c-AID mutations, is most abundant in the G1/S transition and S phase59, but it is possible that in somatic hypermutation and class-switch recombination, UNG exerts its function in G1, similar to AID47,60. Specically, using a mouse model, Sharbeen et al.61 have shown that the mutagenic activity of UNG on AID-induced lesions was exclusively restricted to the G1 phase. These data may suggest that pol-Z is more active in the S phase, while

UNG is more active in G1; how this paradigm applies during the early development and transformation of a B cell to a CLL cell is yet unclear, but our data suggest that pol-Z activity is earlier in this process.

In summary, by characterizing a lower-risk CLL cohort with WGS, we were able to show for the rst time the operation of a distinct c-AID signature in CLL using unsupervised genome-wide analysis, and to demonstrate that this signature is in fact more abundant in cases with lower-risk mutated IGHV. These cases have fewer driver mutations and their key causative events beyond deletion 13q are still unclear. We are continuing to analyse noncoding and promoter regions from these whole genomes, as well as to correlate these results with epigenetic analyses, in an effort to identify the key driving events in these indolent CLL patients. Meanwhile, our new mutational signature detection method can be extended to other cancers to better elucidate signatures associated with clustered mutations.

Methods

Sample preparation. Matched peripheral blood (tumour) and saliva (normal) samples were collected after obtaining informed consent to a tissue banking protocol approved by the Institutional Review Board at Dana-Farber Cancer Institute (protocol no. 99224). For samples with white blood cell count o25 K or absolute lymphocyte count o20 K, B cells were puried using the Easy Sep Human

B cell Enrichment Kit (StemCell Technologies Inc., Vancouver, Canada) according to the manufacturers instructions before viably freezing. Tumour and saliva DNA were extracted using QIAamp Blood DNA (Qiagen Inc., Valencia, CA) and Oragene DNA (Oragene, Ontario, Canada) kits, respectively, according to the manufacturers directions.

Whole-genome sequencing. Puried DNA was submitted to the Genomics Platform at the Broad Institute (Cambridge, MA) for high-throughput whole-genome sequencing. All samples were subjected to in-house quality control (QC) procedures such as Picogreen-based double-stranded DNA quantication (Life Technologies, Carlsbad, CA) and ngerprinting to conrm the match between a tumour and its intended normal, before library preparation.

For a subset of samples, starting with 3 mg of genomic DNA, library construction was performed as described by Fisher et al.62 Another subset of samples, however, was prepared using the protocol by Fisher et al., with some slight modications. Initial genomic DNA input into shearing was reduced from 3 mg to 100 ng in 50 ml of solution. In addition, for adapter ligation, Illumina paired-end adapters were replaced with palindromic forked adapters with unique eight-base index sequences embedded within the adapter. For a subset of samples, size selection was performed using gel electrophoresis, with a target insert size of either 340 or 370 bp10%. Multiple gel cuts were taken for libraries that required high sequencing coverage. For another subset of samples, size selection was performed using Sages Pippin Prep.

Following sample preparation, libraries were quantied using quantitative PCR (kit purchased from KAPA biosystems) with probes specic to the ends of the adapters. Cluster amplication was performed according to the manufacturers protocol (Illumina) using either HiSeq 2000 v2, or HiSeq v3 cluster chemistry and owcells. For a subset of samples, after cluster amplication, SYBR Green dye was added to all owcell lanes, and a portion of each lane was visualized using a light

microscope, to conrm target cluster density. Flowcells were sequenced on HiSeq 2000 using HiSeq 2000 v2 or v3 Sequencing-by-Synthesis kits, then analysed using RTA v1.10.15. or RTA v.1.12.4.2.

Mean target coverage of 30X and 60X was achieved for the tumour and normal samples, respectively. Pilot analysis of two normal saliva samples to determine the percentage of bacterial DNA contamination suggested that 60X coverage would be adequate to achieve 30X human DNA coverage in our samples. Average length of the paired-end reads was 101 bp with an 8-bp index. The raw sequencereads were processed and aligned to the hg19 human reference genome using the Picard pipeline, followed by QC using Firehose tools developed at the Broad22 (https://www.broadinstitute.org/cancer/cga/Firehose

Web End =https://www.broadinstitute.org/cancer/cga/Firehose). The QC parameters tested include lane cross-check ngerprinting for sample identity, tumour normal cross-contamination measured using ContEst63 and coverage statistics. All samples passed the QC check.

Identication of somatic mutations. High-condence somatic mutation calls were made by applying MuTect64 to whole-genome sequencing data from tumours and patient-matched normal samples. Refer to Cibulskis et al.64 for more details. In addition, commonly occurring germline variants were ltered out using a panel of normals. The somatic mutation calls were further subjected to a realignment lter to remove remaining false-positive calls (Supplementary Data 7 for genome-wide somatic single nucleotide variants (sSNV) calls).

Estimation of clonality using ABSOLUTE. Tumour samples are frequently contaminated with normal cells. ABSOLUTE46 infers the purity and ploidy of this heterogeneous population using copy-number and mutation data. ABSOLUTE also estimates local copy number in the cancer cells and the CCF of each mutation (that is, the fraction of cancer cells harbouring the mutation). We followed the same procedure as described in Landau et al.24 (Supplementary Data 8). Specically, mutations with probabilityX0.5 of having CCFX0.95 were classied as clonal, and the rest were classied as subclonal. Mutations with CCF o0.1 were ltered out due to low power.

Discovery of structural rearrangements. Clusters of discordant read pairs were used to infer the presence of structural rearrangements using the dRanger11 and BreakPointer11 algorithms. Mapped distance between pairs that is greater than that expected, based on library insert-size distribution, indicated the presence of a deletion. Inter-chromosomal rearrangements were identied as mate-pairs with each end mapping to different chromosomes. Tandem duplications were identied as pairs with same orientation, as well as an unexpected insert size. dRanger uses a panel of 177 whole-genome-sequenced normals to lter known germline rearrangements and artefacts. The algorithm assigns a nal score based on number of supporting read pairs and a series of ltering matrices described in greater detail previously11. A score cutoff X4 was selected, as previous work has shown that it yields at least 85%

true positives in a large-scale PCR-based validation study11. Breakpointer can be downloaded at https://www.broadinstitute.org/cancer/cga/breakpointer

Web End =https://www.broadinstitute.org/cancer/cga/breakpointer . See Supplementary Data 1 for a list of structural rearrangements.

Analysis of SNP array data. A minimum of 250 ng of tumour and matched normal DNA was used to run Affymetrix Genome-Wide Human SNP Array 6.0 containing 906,600 SNPs and more than 946,000 probes for the detection of copy-number variation on a single genotyping array. The Genome-Wide Human SNP Array 6.0 uses a Birdsuite calling pipeline that delivers SNP as well as CNA calls. Germline CNAs and artefacts were removed by normalizing against a panel of normals. The resultant copy number segments le (seg le) with log2 copy-number ratios was used for further analysis. The number of sCNAs per sample was calculated manually using the following parameters, followed by visual inspection in IGV and comparison with germline: segment length X0.2 MB, amplication threshold X0.1, deletion threshold p 0.1.

The structural rearrangement data and SNP array data were modelled together using the ChainFinder18 algorithm to detect inter-dependent events, as described in detail by Baca et al.18 The algorithm is available at https://www.broadinstitute.org/cancer/cga/chainfinder

Web End =https://www.broadinstitute. https://www.broadinstitute.org/cancer/cga/chainfinder

Web End =org/cancer/cga/chainnder . See Supplementary Data 2 for ChainFinder output.

Discovering mutational signatures. The mutation signatures discovery is a de-convolution process of the somatic mutation counts in each tumour, stratied by mutation contexts and potentially other biologically meaningful parameters, into a set of characteristic mutational signatures. Here we applied the Bayesian NMF algorithm (BayesNMF)5 to infer the number of mutational signatures and their sample-specic contributions. In addition to raw mutation counts stratied by 96 base substitutions in trinucleotide sequence contexts, we also considered the clustering information of mutations as an additional feature in the signature discovery. We considered the NMDs, a minimum genomic distance to all other mutations on the same chromosome in the same patient, as a parameter to stratify mutations, and partitioned them into clustered (NMDp1,000 nt) and non-clustered groups (NMD41,000 nt) (Fig. 4a). The comparison of overall mutation spectrum between clustered and non-clustered mutation groups (Supplementary Fig. 4) revealed a signicant elevation of C4T/G at GCT context and A4G at WA

10 NATURE COMMUNICATIONS | 6:8866 | DOI: 10.1038/ncomms9866 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9866 ARTICLE

(W A/T) in the clustered mutation group, corresponding to the known AID

mutation motifs. On the basis of this observation, we separately counted clustered and non-clustered mutations across 96 trinucleotide mutation contexts in each sample. We split mutations in each tumour into two columns representing clustered and non-clustered mutational groups, giving rise to the mutation count matrix X (96 by 2M, M the number of samples). This mutation count matrix was ingested as an

input for the BayesNMF and factored into two matrices, W0 (96 by K) and H0 (K by 2M), approximating X by W0H0. It should be noted that clustered and non-clustered mutations from the same patient were separately handled to capture a characteristic signal from clustered mutations. While the conventional NMF requires the number of signatures K a priori, BayesNMF automatically prunes away irrelevant components that do not contribute to explaining X and effectively determines the appropriate number of K. We ran BayesNMF 50 times with exponential priors for W0 and H0 and 41 out of the 50 runs converged to the solution of K 3, while 9 runs

converged to the solution of K 4. We used the three-signature solution (K 3)

with the maximum posterior for the downstream analysis. To enumerate the number of mutations associated with each mutation signature, we performed a scaling transformation, XBW0H0 WH, W W0U 1 and H UH0, where U is a K-by-K

diagonal matrix with the element corresponding to the 1-norm of column vectors of W0, resulting in the nal signature matrix W and the activity matrix H. Note that the kth column vector of W (wk) represents a normalized mutability of 96 trinucleotide mutation contexts in the kth signature and the kth row vector of H (hk) dictates the estimation of clustered and non-clustered mutations associated with the kth signature across samples (Fig. 4c).

Signature-enrichment analysis. Using the determined W and H from the BayesNMF, we annotated each mutation with the probability (likelihood of association) that it was generated by each of the discovered mutational signatures, pms,

where m denoted a mutation and s refers to the signature. More specically, the likelihood of association to the kth signature for a set of mutations corresponding to ith mutation context and jth clustered or non-clustered mutation group was dened as [wkhk/Skwkhk]ij, where wk and hk correspond to the kth column vector and kth row vector of W and H, respectively.

For the gene-level signature-enrichment analysis, we rst attempted to identify a hotspot mutation motif out of 96 contexts in each signature by considering mutations only with pms 40.75. Note that keeping mutations with a higher pms, ltered out mutations shared by multiple signatures and enabled the discovery of more distinct mutation motifs unique to each signature. By considering contributions more than the third quintile, we were able to extract characteristic mutation motifs to each signature19 hotspot mutation motifs in c-AID and ve hotspot mutation motifs in nc-AID (Supplementary Fig. 5a). To take into account sequence composition variation across the genome, we enumerated all available trinucleotide contexts across genes having non-zero mutations with pms40.75 in

each signature. This information was used to estimate the background mutation rates at the hotspot motifs in each signature, resulting in rcAID 0.27 per Mb and

rncAID 0.58 per Mb for c-AID and nc-AID signatures, respectively. Then, for

given mutation counts, x, at hotspot motifs and available sequence context, n, in each gene, we performed a binomial test with the estimated background mutation rate to assess the signicance of the enrichment of each signature across 281 genes for c-AID and 809 genes for nc-AID having non-zero mutations with pms 40.75

(Supplementary Fig. 5b; Supplementary Data 3 and 4). We corrected for multiple hypotheses and identied genes that are associated with each signature using a q-value cutoff of 0.1 (see QQ plots in Supplementary Fig. 5c).

Two threshold values (0.5 and 0.75, Supplementary Data 5 and 6) for pms were

utilized for the clonality analysis to dichotomize the signature association of mutations.

RNA sequencing and analysis. RNA was extracted using the Qiagen RNeasy kit and the RNA integrity number was measured using Agilent Bioanalyzer at the Harvard BioPolymers Facility to assess the quality of the extracted RNA. Only samples with a RNA Integrity Number 48 were submitted for sequencing. Poly-A-selected RNA was used for library construction using the Illumina TruSeq Paired End Strand-specic kit according to the manufacturers protocol and sequenced using Illumina HiSeq. The RNA-seq BAMs were aligned to the hg19 genome using TopHat65 (Gencode gtf used for annotation). QC analysis was performed using the metrics described by DeLuca et al.66 Gene-level expression data represented as fragments per kilobase of exons mapped was obtained using Cufinks67.

Statistical analysis. Statistical analysis was performed using with SAS version 9.2 (SAS Institute, Cary, NC) and R version 2.15.2 (the CRAN project). Categorical variables were compared using the Fishers exact test or a w2-test as appropriate, and continuous variables were compared using the Wilcoxons rank-sum test. TTNT was dened as the time of sampling to the rst treatment after sampling or death, whichever occurs rst. Patients who did not receive a treatment after sampling were censored at the date last known alive and without any treatment. TTNT was estimated using the Kaplan and Meier method, and the difference was tested using the log-rank test. In addition, univariable Cox modelling was performed for known CLL risk factors as well as exploratory factors presented in this paper. Due to the limited number of events, multivariable Cox modelling was not

explored. The linearity assumption for continuous variables was examined using restricted cubic spline estimates of the relationship between the continuous variable and log-relative hazard, and the cutoff points of these variables were based on the change of the log-relative hazards. All P values are two sided and considered signicant at the 0.05 level. Due to the exploratory nature, multiple comparisons were not adjusted in the signicance level.

References

1. Hallek, M. et al. Guidelines for the diagnosis and treatment of chronic lymphocytic leukemia: a report from the International Workshop on Chronic Lymphocytic Leukemia updating the National Cancer InstituteWorking Group 1996 guidelines. Blood 111, 54465456 (2008).

2. Dohner, H. et al. Genomic aberrations and survival in chronic lymphocytic leukemia. N. Engl. J. Med. 343, 19101916 (2000).

3. Puente, X. S. et al. Whole-genome sequencing identies recurrent mutations in chronic lymphocytic leukaemia. Nature 475, 101105 (2011).

4. Alexandrov, L. B. et al. Signatures of mutational processes in human cancer. Nature 500, 415421 (2013).

5. Tan, V. Y. & Fevotte, C. Automatic relevance determination in nonnegative matrix factorization with the beta-divergence. IEEE Trans. Pattern Anal. Mach. Intell. 35, 15921605 (2013).

6. Nik-Zainal, S. et al. Mutational processes molding the genomes of 21 breast cancers. Cell 149, 979993 (2012).

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

8. Pettersen, H. S. et al. AID expression in B-cell lymphomas causes accumulation of genomic uracil and a distinct AID mutational signature. DNA Repair (Amst) 25, 6071 (2015).

9. Patten, P. E. et al. IGHV-unmutated and IGHV-mutated chronic lymphocytic leukemia cells produce activation-induced deaminase protein with a full range of biologic functions. Blood 120, 48024811 (2012).

10. Roberts, S. A. et al. Clustered mutations in yeast and in human cancers can arise from damaged long single-strand DNA regions. Mol. Cell 46, 424435 (2012).

11. Pasqualucci, L. et al. Hypermutation of multiple proto-oncogenes in B-cell diffuse large-cell lymphomas. Nature 412, 341346 (2001).

12. 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, 228235 (2013).

13. Hruba, M., Dvorak, P., Weberova, L. & Subrt, I. Independent coexistence of clones with 13q14 deletion at reciprocal translocation breakpoint and 13q14 interstitial deletion in chronic lymphocytic leukemia. Leuk. Lymphoma 53, 20542062 (2012).

14. Reindl, L. et al. Biological and clinical characterization of recurrent 14q deletions in CLL and other mature B-cell neoplasms. Br. J. Haematol. 151, 2536 (2010).

15. Quintero-Rivera, F., Nooraie, F. & Rao, P. N. Frequency of 5IGH deletions in B-cell chronic lymphocytic leukemia. Cancer Genet. Cytogenet. 190, 3339 (2009).

16. Brown, J. R. et al. Integrative genomic analysis implicates gain of PIK3CA at 3q26 and MYC at 8q24 in chronic lymphocytic leukemia. Clin. Cancer Res. 18, 37913802 (2012).

17. Pfeifer, D. et al. Genome-wide analysis of DNA copy number changes and LOH in CLL using high-density SNP arrays. Blood 109, 12021210 (2007).18. Stephens, P. J. et al. Massive genomic rearrangement acquired in a single catastrophic event during cancer development. Cell 144, 2740 (2011).

19. Baca, S. C. et al. Punctuated evolution of prostate cancer genomes. Cell 153, 666677 (2013).

20. Korbel, J. O. & Campbell, P. J. Criteria for inference of chromothripsis in cancer genomes. Cell 152, 12261236 (2013).

21. Zhang, C.-Z. et al. Chromothripsis from DNA damage in micronuclei. Nature 522, 179184 (2015).

22. Morin, R. D. et al. Mutational and structural analysis of diffuse large B-cell lymphoma using whole-genome sequencing. Blood 122, 12561265 (2013).

23. Chapman, M. A. et al. Initial genome sequencing and analysis of multiple myeloma. Nature 471, 467472 (2011).

24. Landau, D. A. et al. Evolution and impact of subclonal mutations in chronic lymphocytic leukemia. Cell 152, 714726 (2013).

25. Lawrence, M. S. et al. Discovery and saturation analysis of cancer genes across 21 tumour types. Nature 505, 495501 (2014).

26. de Miranda, N. F. et al. Exome sequencing reveals novel mutation targets in diffuse large B-cell lymphomas derived from Chinese patients. Blood 124, 25442553 (2014).

27. Alexandrov, L. B., Nik-Zainal, S., Wedge, D. C., Campbell, P. J. & Stratton, M.R. Deciphering signatures of mutational processes operative in human cancer. Cell Rep. 3, 246259 (2013).28. Roberts, S. A. & Gordenin, D. A. Hypermutation in human cancer genomes: footprints and mechanisms. Nat. Rev. Cancer 14, 786800 (2014).

NATURE COMMUNICATIONS | 6:8866 | DOI: 10.1038/ncomms9866 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 11

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9866

29. Victora, G. D. & Nussenzweig, M. C. Germinal centers. Annu. Rev. Immunol.

30, 429457 (2012).

30. Yamane, A. et al. Deep-sequencing identication of the genomic targets of the cytidine deaminase AID and its cofactor RPA in B lymphocytes. Nat. Immunol. 12, 6269 (2011).

31. Di Noia, J. M. & Neuberger, M. S. Molecular mechanisms of antibody somatic hypermutation. Annu. Rev. Biochem. 76, 122 (2007).

32. Betz, A. G., Rada, C., Pannell, R., Milstein, C. & Neuberger, M. S. Passenger transgenes reveal intrinsic specicity of the antibody hypermutation mechanism: clustering, polarity, and specic hot spots. Proc. Natl Acad. Sci. USA 90, 23852388 (1993).

33. Stavnezer, J. Complex regulation and function of activation-induced cytidine deaminase. Trends Immunol. 32, 194201 (2011).

34. Canugovi, C., Samaranayake, M. & Bhagwat, A. S. Transcriptional pausing and stalling causes multiple clustered mutations by human activation-induced deaminase. FASEB J. 23, 3444 (2009).

35. Albesiano, E. et al. Activation-induced cytidine deaminase in chronic lymphocytic leukemia B cells: expression as multiple forms in a dynamic, variably sized fraction of the clone. Blood 102, 33333339 (2003).

36. Brown, J. R. et al. Next-generation sequencing reveals clonal evolution at the immunoglobulin loci in chronic lymphocytic leukemia. Blood 124, 3302 (2014).

37. Messmer, B. T., Albesiano, E., Messmer, D. & Chiorazzi, N. The pattern and distribution of immunoglobulin VH gene mutations in chronic lymphocytic leukemia B cells are consistent with the canonical somatic hypermutation process. Blood 103, 34903495 (2004).

38. Huemer, M. et al. AID induces intraclonal diversity and genomic damage in CD86( ) chronic lymphocytic leukemia cells. Eur. J. Immunol. 44, 37473757

(2014).39. Capello, D. et al. Distribution and pattern of BCL-6 mutations throughout the spectrum of B-cell neoplasia. Blood 95, 651659 (2000).

40. Khodabakhshi, A. H. et al. Recurrent targets of aberrant somatic hypermutation in lymphoma. Oncotarget 3, 13081319 (2012).

41. Ranzani, M., Annunziato, S., Adams, D. J. & Montini, E. Cancer gene discovery: exploiting insertional mutagenesis. Mol. Cancer Res. 11, 11411158 (2013).

42. He, W. et al. Aberrant methylation and loss of CADM2 tumor suppressor expression is associated with human renal cell carcinoma tumor progression. Biochem. Biophys. Res. Commun. 435, 526532 (2013).

43. Raufman, J.-P. et al. Muscarinic receptor subtype-3 gene ablation and scopolamine butylbromidetreatment attenuate small intestinal neoplasia in apcmin/ mice. Carcinogenesis 32, 13961402 (2011).

44. Kan, Z. et al. Diverse somatic mutation patterns and pathway alterations in human cancers. Nature 466, 869873 (2010).

45. Chang, P. H. et al. Activation of Robo1 signaling of breast cancer cells by Slit2 from stromal broblast restrains tumorigenesis via blocking PI3K/Akt/beta-catenin pathway. Cancer Res. 72, 46524661 (2012).

46. Carter, S. L. et al. Absolute quantication of somatic DNA alterations in human cancer. Nat. Biotechnol. 30, 413421 (2012).

47. Zan, H. & Casali, P. Regulation of Aicda expression and AID activity. Autoimmunity 46, 83101 (2013).

48. Oppezzo, P. et al. Chronic lymphocytic leukemia B cells expressing AID display dissociation between class switch recombination and somatic hypermutation. Blood 101, 40294032 (2003).

49. Heintel, D. et al. High expression of activation-induced cytidine deaminase (AID) mRNA is associated with unmutated IGVH gene status and unfavourable cytogenetic aberrations in patients with chronic lymphocytic leukaemia. Leukemia 18, 756762 (2004).

50. Greaves, M. & Maley, C. C. Clonal evolution in cancer. Nature 481, 306313 (2012).

51. Keim, C., Kazadi, D., Rothschild, G. & Basu, U. Regulation of AID, the B-cell genome mutator. Genes Dev. 27, 117 (2013).

52. Muramatsu, M. et al. Specic expression of activation-induced cytidine deaminase (AID), a novel member of the RNA-editing deaminase family in germinal center B cells. J. Biol. Chem. 274, 1847018476 (1999).

53. Rebhandl, S. et al. APOBEC3 signature mutations in chronic lymphocytic leukemia. Leukemia 28, 19291932 (2014).

54. Gaidano, G., Foa, R. & Dalla-Favera, R. Molecular pathogenesis of chronic lymphocytic leukemia. J. Clin. Invest. 122, 34323438 (2012).

55. Lanasa, M. C. Novel insights into the biology of CLL. Hematology Am. Soc. Hematol. Educ. Program 2010, 7076 (2010).

56. Messmer, B. T. et al. In vivo measurements document the dynamic cellular kinetics of chronic lymphocytic leukemia B cells. J. Clin. Invest. 115, 755764 (2005).57. Zeng, X. et al. DNA polymerase eta is an A-T mutator in somatic hypermutation of immunoglobulin variable genes. Nat. Immunol. 2, 537541 (2001).

58. Yamada, A., Masutani, C., Iwai, S. & Hanaoka, F. Complementation of defective translesion synthesis and UV light sensitivity in xeroderma pigmentosum variant cells by human and mouse DNA polymerase eta. Nucleic Acids Res. 28, 24732480 (2000).

59. Visnes, T. et al. Uracil in DNA and its processing by different DNA glycosylases. Philos. Trans. R. Soc. Lond. B Biol. Sci. 364, 563568 (2009).

60. Li, S., Zhao, Y. & Wang, J. Y. Analysis of Ig gene hypermutation in Ung(-/-) Polh(-/-) mice suggests that UNG and A:T mutagenesis pathway target different U:G lesions. Mol. Immunol. 53, 214217 (2013).

61. Sharbeen, G., Yee, C. W., Smith, A. L. & Jolly, C. J. Ectopic restriction of DNA repair reveals that UNG2 excises AID-induced uracils predominantly or exclusively during G1 phase. J. Exp. Med. 209, 965974 (2012).

62. Fisher, S. et al. A scalable, fully automated process for construction of sequence-ready human exome targeted capture libraries. Genome Biol. 12, R1 (2011).

63. Cibulskis, K. et al. ContEst: estimating cross-contamination of human samples in next-generation sequencing data. Bioinformatics 27, 26012602 (2011).64. Cibulskis, K. et al. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat. Biotechnol. 31, 213219 (2013).

65. Trapnell, C., Pachter, L. & Salzberg, S. L. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25, 11051111 (2009).

66. DeLuca, D. S. et al. RNA-SeQC: RNA-seq metrics for quality control and process optimization. Bioinformatics 28, 15301532 (2012).

67. Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufinks. Nat. Protoc. 7, 562578 (2012).

Acknowledgements

We thank the Melton Family Fund for CLL Research and the National Human Genome Research Institute Grant U54HG003067 for funding the whole-genome sequencing and SNP arrays, respectively. B.J.R. is a clinical scholar of the Leukemia Lymphoma Society and is supported by the American Cancer Society (RSG-13-002-01-CCE) and the Leukemia Lymphoma Society (TRP#6289-13). G.G. is the Paul C. Zamecnik, MD, Chair in Oncology at Massachusetts General Hospital. We thank the Genome Sequencing Platform and Firehose Team at the Broad Institute for their help. We also thank Drs. Nicholas Chiorazzi, Catherine Wu and Dan-Avi Landau for critical review of this manuscript and Dr. Sylvan Baca for reviewing the ChainFinder results.

Author contributions

K.S. performed the experiments, analysed the data and wrote the manuscript. K.J. developed the modied Bayesian NMF algorithm, analysed the data and wrote the manuscript. I.R. helped with sample selection and edited the manuscript. T.G., P.P., H.N. and K.A. analysed the data. L.M.S. provided critical insights for bioinformatics analyses. F.S.M. collected and processed patient samples. B.S. and S.C. were the project managers for the sequencing samples submitted to the Genomics Platform at Broad Institute. G.S. and L.E.S. are the co-PIs of the NHGRI grant U54HG003067 used for obtaining the SNP array data. K.H.T. performed the statistical analyses. G.G. and B.J.R. conceived and directed the study and wrote the manuscript. B.J.R. funded the sequencing project.

Additional information

Accession codes: The sequence data have been deposited in dbGAP under the accession code phs000879.v1.p1 http://npg.nature.com/reprintsandpermissions/

Web End =http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi? http://npg.nature.com/reprintsandpermissions/

Web End =study_id=phs000879.v1.p1 .

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Kasar, S. et al. Whole-genome sequencing reveals activation-induced cytidine deaminase signatures during indolent chronic lymphocytic leukaemia evolution. Nat. Commun. 6:8866 doi: 10.1038/ncomms9866 (2015).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

12 NATURE COMMUNICATIONS | 6:8866 | DOI: 10.1038/ncomms9866 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications