DNA deaminases induce break-associated mutation showers with implication of APOBEC3B and 3A in breast cancer kataegis

  1. Benjamin JM Taylor
  2. Serena Nik-Zainal
  3. Yee Ling Wu
  4. Lucy A Stebbings
  5. Keiran Raine
  6. Peter J Campbell
  7. Cristina Rada
  8. Michael R Stratton
  9. Michael S Neuberger  Is a corresponding author
  1. Medical Research Council Laboratory of Molecular Biology, United Kingdom
  2. Cancer Genome Project, Wellcome Trust Sanger Institute, United Kingdom
4 figures, 1 table and 1 additional file

Figures

Figure 1 with 1 supplement
AID/APOBEC-induced kataegis in yeast.

(A) The total number of mutations detected in canavanine-resistant (CanR) AID/APOBEC yeast transformants, median frequency indicated. (B) Rainfall plot of genome-wide intermutational distances (IMD) in individual CanR AID* transformants. Clone identifier is indicated along the top, mutations shown as dots with the y-axis giving the distance to the next downstream mutation on the same chromosome. For each clone, dots are ordered sequentially along the genome. Dot colours represent: cluster mutations (at C, red; G, black), unclustered mutations (at C, pink; G, grey). Mutations at A:T (14 out of 1078 total), single mutations on individual chromosomes (15% of the database), the most downstream mutation in each chromosome and transformants without any multiply mutated chromosomes (5/40) are not depicted. Supplementary file 1B contains the location of all identified mutations. (C) Observed distribution of IMDs in the AID* dataset compared to a simulation assuming mutations are randomly scattered throughout the genome. (D) IMD plots of APOBEC3A/B/G*-expressing yeast transformants (28/78 transformants harboured no multiply mutated chromosome and are not depicted). (E) Detailed view of AID* mutation clusters. Each line represents an individual cluster with the clone identifier (grey box), number of mutations in the cluster and total mutations in the clone indicated. Mutations are coloured as in (B), a horizontal line indicates mutations that have coalesced, * indicates clusters localising within 10 kb of CAN1. All clusters containing ≥5 mutations are depicted. (F) Mutation clusters identified in yeast APOBEC3 transformants.

https://doi.org/10.7554/eLife.00534.003
Figure 1—figure supplement 1
Characterisation of AID/APOBEC yeast transformants.

(A) CanR frequencies of yeast transformants expressing different AID/APOBEC proteins. (B) AID/APOBEC protein levels in yeast transformants. Whole cell extracts from 24 hr galactose induced yeast transformants were immunoblotted with anti-FLAG (M2; Sigma) and anti-β-Actin (ab8224; Abcam).

https://doi.org/10.7554/eLife.00534.004
Figure 2 with 2 supplements
Yeast kataegic clusters are associated with transversions, are reduced by UNG-deficiency and can be triggered by a double strand DNA break.

(A) IMD plots of AID*/APOBEC transformants reveal preferential association of nucleotide transversions with kataegic clusters. Mutation datasets and presentation are as in Figure 1B but with transition mutations represented by yellow dots and transversions by blue dots. Density plots depict the overall distribution of transition (Tn) and transversion (Tv) mutations at C:G pairs. (B) IMD plots of AID*/APOBEC-expressing ung1Δ yeast transformants, depicted as in Figure 2A. Density plots compare the distributions of IMDs in AID* transformants of ung1Δ and wild type yeast. (C) All mutation clusters identified in AID*/APOBEC3A-transformants of ung1Δ yeast depicted as in Figure 1E. (D) Kataegis localised to a double strand break. Mutations in the vicinity of the CAN1 locus of (I-SceI+APOBEC3G*) transformants of either control cells or of a KanMX-ISceIRS derivative carrying a CAN1-proximal I-SceI recognition sequence. The I-SceI cut site is marked with an arrow. All CAN1 mutations in control cells and 33/36 CAN1 region mutations in KanMX-ISceIRS cells occur at the canonical APOBEC3G CC context. Two-thirds of the CAN1 region mutations in the KanMX-ISceIRS cells were transversions.

https://doi.org/10.7554/eLife.00534.005
Figure 2—figure supplement 1
Canavanine resistance frequencies of yeast transformants carrying an I-SceI recognition sequence (I-SceIRS).

CanR mutation frequency of yeast cells carrying an I-SceIRS 420 bp downstream of the CAN1 polyadenylation site (KanMX-ISceIRS cells) which have been induced for APOBEC3G* and/or I-SceI expression. * indicated p<0.05.

https://doi.org/10.7554/eLife.00534.006
Figure 2—figure supplement 2
IMD plots of mutations from CanR (I-SceI+APOBEC3G*) transformants of control or KanMX-ISceIRS cells IMD plots coloured as in Figure 1B.
https://doi.org/10.7554/eLife.00534.007
Figure 3 with 4 supplements
Comparison of kataegic mutations in yeast AID/APOBEC transformants with those in breast cancers.

(A) Comparison of the length, number of mutations and polarity of yeast kataegic clusters with those in breast cancers. The degree of strand polarity is indicated by colour intensity. The breast cancer data (Nik-Zainal et al., 2012) are a compilation from three tumours (PD4103a, PD4107a, PD4199a) chosen for their large number of clusters. (B) Context of the genome wide mutated C bases in yeast AID/APOBEC transformants with total numbers of mutations in each dataset indicated. (C) Context of the kataegic and singlet mutated C bases in selected breast cancers. Analyses of all sequenced breast cancers are presented in Figure 3—figure supplements 2–4. (D) Similarity of sequence contexts of C mutations in breast cancer kataegic stretches compared to those of deaminase-induced C mutations in yeast. (D1) Identity of the base at the −2 position of TC mutations in cancer kataegic regions and in APOBEC3A/B yeast transformants. The base compositions were normalised to the genomic base composition of the −2 base at TC dinucleotides. (D2) Sequence contexts similarity p-value at positions (−1 plus −2) to the mutated Cs. The contexts of all Cs throughout the yeast and human genomes are included for comparison. Mutation context of wild type versions of AID and APOBEC3G are shown in Figure 3—figure supplement 1. Analysis of additional yeast transformants and breast cancers is shown in Figure 3—figure supplements 2–4.

https://doi.org/10.7554/eLife.00534.009
Figure 3—figure supplement 1
Mutation context of hyperactive APOBEC3G and AID are identical to the wild type proteins.

Mutation context of CanR transformants expressing APOBEC3G and AID.

https://doi.org/10.7554/eLife.00534.010
Figure 3—figure supplement 2
Analysis of kataegic stretches and mutation distributions of 21 breast cancers.

Analysis of tumours PD3851a to PD4086a. For each tumour, somatic mutations are displayed as intermutational distance (IMD) plots (top left), as described in Figure 1B. Dots are coloured according to kataegic mutations (black) or singlet mutations (red). The sequence context of C mutations in unclustered and kataegic locations are shown (lower plots), with the number of mutations indicated, coloured as in Figure 3B. The length and number of mutations within each cluster is displayed (top right), with the dot colour intensity indicating polarity. Data are taken from Nik-Zainal et al. (2012).

https://doi.org/10.7554/eLife.00534.011
Figure 3—figure supplement 3
Analysis of kataegic stretches and mutation distributions of 21 breast cancers.

Analysis of tumours PD4088a to PD4194a. For each tumour, somatic mutations are displayed as intermutational distance (IMD) plots (top left), as described in Figure 1B. Dots are coloured according to kataegic mutations (black) or singlet mutations (red). The sequence context of C mutations in unclustered and kataegic locations are shown (lower plots), with the number of mutations indicated, coloured as in Figure 3B. The length and number of mutations within each cluster is displayed (top right), with the dot colour intensity indicating polarity. Data are taken from Nik-Zainal et al. (2012).

https://doi.org/10.7554/eLife.00534.012
Figure 3—figure supplement 4
Analysis of kataegic stretches and mutation distributions of 21 breast cancers.

Analysis of tumours PD4198a to PD4248a. For each tumour, somatic mutations are displayed as intermutational distance (IMD) plots (top left), as described in Figure 1B. Dots are coloured according to kataegic mutations (black) or singlet mutations (red). The sequence context of C mutations in unclustered and kataegic locations are shown (lower plots), with the number of mutations indicated, coloured as in Figure 3B. The length and number of mutations within each cluster is displayed (top right), with the dot colour intensity indicating polarity. Data are taken from Nik-Zainal et al. (2012).

https://doi.org/10.7554/eLife.00534.013
DNA damage by APOBEC3 family members and expression in breast cancer cell-lines.

(A) Effect of enforced APOBEC expression on cell viability. (A1) Stable transfectants of KBM7 cells that inducibly express APOBEC proteins were incubated with inducer (doxycyclin) and viability was monitored after 72 hr. (A2) Expression of FLAG-tagged APOBECs after 24 hr doxycyclin treatment. (B) Enforced expression of APOBEC3A and APOBEC3B leads to induction of histone γH2AX and of 53BP1 foci. The percentage of cells (B1) positive for histone γH2AX expression quantified by flow cytometry and (B2) exhibiting punctate rather than diffuse 53BP1 staining quantified by immunofluorescence microscopy, with the foci number per cell indicated. (C) Expression of APOBEC family members in six human breast cancer cell-lines as well as in HEK293 cells was analysed by qRT-PCR of total cellular RNA. Expression is shown relative to the average of housekeeping genes HPRT and HMBS. The effect of phorbol ester (PMA) and interferon alpha (INFα) treatment on APOBEC3A and APOBEC3B levels is also shown. In all cases, * indicates p<0.1, ** indicates p<0.0001 compared to control (unpaired t-test).

https://doi.org/10.7554/eLife.00534.014

Tables

Table 1

Pattern of nucleotide substitutions at C:G pairs in AID*/APOBEC yeast transformants. All mutations at C:G pairs were computed as substitutions at C

https://doi.org/10.7554/eLife.00534.008
AID*ung1Δ AID*rev1Δ AID*A3AA3BA3G*
KataegicUnclusteredTotal
C→T (%)46877699.3100798178
C→G (%)4711210.20171620
C→A (%)7230.50432
Total2747901064208877130121560

Additional files

Supplementary file 1

(A) Catalogue of yeast mutations. (B) Singlet and kataegic C:G mutations in yeast transformants. Kataegic mutations in all yeast transformants, are defined as those mutations that form part of a group of at least 5 mutations in which each mutation is separated from its downstream neighbour by ≤8.5 kb. All other mutations (including those that constitute the sole mutation on an individual chromosome in a particular cell) are classed as singlet mutations. (C) Primers used for APOBEC RNA RT-PCR.

https://doi.org/10.7554/eLife.00534.015

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Benjamin JM Taylor
  2. Serena Nik-Zainal
  3. Yee Ling Wu
  4. Lucy A Stebbings
  5. Keiran Raine
  6. Peter J Campbell
  7. Cristina Rada
  8. Michael R Stratton
  9. Michael S Neuberger
(2013)
DNA deaminases induce break-associated mutation showers with implication of APOBEC3B and 3A in breast cancer kataegis
eLife 2:e00534.
https://doi.org/10.7554/eLife.00534