Polyploid cells rewire DNA damage response networks to overcome replication stress-induced barriers for tumour progression

Abstract

Mutations in genes involved in DNA replication, such as flap endonuclease 1 (FEN1), can cause single-stranded DNA breaks (SSBs) and subsequent collapse of DNA replication forks leading to DNA replication stresses. Persistent replication stresses normally induce p53-mediated senescence or apoptosis to prevent tumour progression. It is unclear how some mutant cells can overcome persistent replication stresses and bypass the p53-mediated pathways to develop malignancy. Here we show that polyploidy, which is often observed in human cancers, leads to overexpression of BRCA1, p19arf and other DNA repair genes in FEN1 mutant cells. This overexpression triggers SSB repair and non-homologous end-joining pathways to increase DNA repair activity, but at the cost of frequent chromosomal translocations. Meanwhile, DNA methylation silences p53 target genes to bypass the p53-mediated senescence and apoptosis. These molecular changes rewire DNA damage response and repair gene networks in polyploid tumour cells, enabling them to escape replication stress-induced senescence barriers.

ARTICLE

Received 10 Feb 2012 | Accepted 5 Apr 2012 | Published 8 May 2012 DOI: 10.1038/ncomms1825

Li Zheng1, Huifang Dai1, Mian Zhou1, Xiaojin Li2, Changwei Liu1,3, Zhigang Guo1,, Xiwei Wu2, Jun Wu4, Charles Wang2, John Zhong5, Qin Huang6, Julio Garcia-Aguilar7, Gerd P. Pfeifer1 & Binghui Shen1,3

Mutations in genes involved in DNA replication, such as ap endonuclease 1 (FEN1), can cause single-stranded DNA breaks (SSBs) and subsequent collapse of DNA replication forks leading to DNA replication stresses. Persistent replication stresses normally induce p53-mediated senescence or apoptosis to prevent tumour progression. It is unclear how some mutant cells can overcome persistent replication stresses and bypass the p53-mediated pathways to develop malignancy. Here we show that polyploidy, which is often observed in human cancers, leads to overexpression of BRCA1, p19arf and other DNA repair genes in FEN1 mutant cells. This overexpression triggers SSB repair and non-homologous end-joining pathways to increase DNA repair activity, but at the cost of frequent chromosomal translocations. Meanwhile, DNA methylation silences p53 target genes to bypass the p53-mediated senescence and apoptosis. These molecular changes rewire DNA damage response and repair gene networks in polyploid tumour cells, enabling them to escape replication stress-induced senescence barriers.

1 Department of Cancer Biology, City of Hope National Medical Center and Beckman Research Institute, 1500 East Duarte Road, Duarte, California 91010, USA. 2 Department of Molecular Medicine, City of Hope National Medical Center and Beckman Research Institute, 1500 East Duarte Road, Duarte, California 91010, USA. 3 Institute of Genetics, College of Life Sciences, Zhejiang University, Hangzhou, China. 4 Division of Comparative Medicine, Cityof Hope National Medical Center and Beckman Research Institute, 1500 East Duarte Road, Duarte, California 91010, USA. 5 Department of Pathology, University of Southern California, 1501 San Pablo St, ZNI 529, Los Angeles, 90033, USA. 6 Department of Pathology, City of Hope National Medical Center and Beckman Research Institute, 1500 East Duarte Road, Duarte, California 91010, USA. 7 Department of Surgery, City of Hope National Medical Center and Beckman Research Institute, 1500 East Duarte Road, Duarte, California 91010, USA. Present address: College of Life Sciences, Nanjing Normal University, China. Correspondence and requests for materials should be addressed to L.Z. (email: [email protected]) or B.S. (email: [email protected]).

Polyploid cells rewire DNA damage response networks to overcome replication stress-induced barriers for tumour progression

NATURE COMMUNICATIONS | 3:815 | DOI: 10.1038/ncomms1825 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1825

DNA single-strand breaks (SSBs) frequently occur during DNA replication and repair1. The most well-known example occurs during replication of the lagging DNA strand, a proc

ess in which millions of short DNA fragments (Okazaki fragments) are synthesised. Okazaki fragments are processed by the combined action of DNA polymerase , ap endonuclease 1 (FEN1) and DNA ligase 1, which seal the nicks between each fragment to produce the intact lagging DNA strand24. Mutations in any one of these genes, including FEN1, can impair Okazaki fragment maturation, leading to persistent DNA SSBs in the genome1,5,6. In addition, SSBs or gaps are transiently produced by nuclease cleavage during DNA repair, including mismatch repair, nucleotide excision repair and base excision repair1,5,6. Furthermore, SSBs may also arise directly aer exposure to DNA damaging agents such as ionising radiation1. Mammalian cells have evolved DNA SSB repair pathways and are relatively tolerant to SSBs. Critically, however, unrepaired DNA SSBs can collapse DNA replication forks, leading to replication stresses and one-ended DNA double-strand breaks (DSBs), which are the most challenging lesions for cells to repair1,7,8.

Cells have evolved multiple DNA damage response and repair pathways to deal with collapsed replication forks and DSBs to prevent transmission of DNA damage and genome instabilities to daughter cells. In response to collapsed replication forks, cells activate DNA checkpoints and arrest cell cycle progression to allow the repair of the stalled forks. Accurate repair of the one-ended DSBs and restart of the collapsed replication fork require the homologous recombination (HR) repair machinery and a homologous DNA sequence as the repair template7,8. In rare cases, a one-ended DSB may be misaligned with another DSB and misjoined by the non-homologous end-joining (NHEJ) machinery7. This misjoining results in chromosomal translocations, which are among the most common cytogenetic aberrations in cancer cells9. Because faithful duplication and transmission of the chromosomes into daughter cells is critical for the survival of organisms, it is thought that NHEJ complexes are inhibited at one-ended DNA DSBs7. In the case of continued DNA-replication stress and persistent formation of DSBs, p53-mediated cellular senescence or apoptosis is induced to prevent proliferation of damaged cells1012.

It is unclear how tumour cells breach this replication-stress-mediated barrier and progress into malignant cancers. To address this question, we used a previously established mutant mouse model carrying a point mutation (F343A/F344A, FFAA) in the FEN1 gene that encodes an endonuclease critical for the Okazaki fragment maturation6,13. Mutant FFAA mouse cells display defects in the Okazaki fragment maturation and consequently high levels of unli-gated DNA SSBs and DSBs6. These cells also have high potential to become tetraploid and aneuploid6. Here we show that the near-polyploid aneuploid cancer cells can induce epigenetic changes that lead to the overexpression of BRCA1, p19arf and other DNA repair genes, but downregulate p53 target genes in the senescence and apoptosis pathways. As a consequence, the DNA damage response and repair networks are rewired in the near-polyploid aneuploid cancer cells, leading to reduction of DNA replication stresses and escaping of senescence and apoptosis.

ResultsPolyploid tumour cells overcome ATR-mediated senescence.

We previously showed that heterozygous mutant mice harbouring the FFAA mutation in FEN1, which causes defects in the Okazaki fragment maturation were prone to the development of malignant cancers6,13. To determine how WT/FFAA cells develop malignancy, we examined their response to DNA replication stresses and sought to identify the underlying molecular events. Previously, we showed that primary WT/FFAA mouse embryonic broblasts (MEFs) grown in normal culture conditions have the capacity to form large numbers of foci6. Cells from these foci were thought to be tumour-

initiating cells14,15. To characterise these tumour-initiating cells, we conducted serial expansion of these cells (Fig. 1a). We found that 23% of the foci could expand unlimitedly with the remaining 77%, displaying only limited clonal expansion capability. The unlimitedly expanding cells displayed a fast and unrestricted proliferation, compared with the normal diploid MEF cells, which showed a gradually deceased proliferation (Fig. 1b). Consistently, they had signicantly higher proportions of S- and G2-phase cells than those of primary diploid WT/FFAA MEF cells (Supplementary Fig S1), and formed subcutaneous cancers in the immunodecient NOD/ SCID mice. In contrast, cells capable of only limited expansion showed little growth at 1015 days (Fig. 1b). Previously, the FFAA FEN1 mutation was shown to induce tetraploidy6, which presumably causes more DNA replication stresses due to the high DNA content of the tetraploid cells. To determine if tetraploidy was associated with the expansion phenotypes of WT/FFAA tumour cells, we counted the chromosome number of these two types of tumour cells. Surprisingly, all the unlimitedly expanding clones had near-polyploid aneuploidy (hereaer referred to as aneuploidy; Fig. 1c and Supplementary Fig. S2). In contrast, cells of three limitedly expanding clones were diploid or near-diploid (Fig. 1c).

To determine the activation of DNA damage response pathways and cellular senescence in primary WT/FFAA MEFs, and the limitedly or unlimitedly expanding tumour cells, we immunouorescently stained the cells for the phosphorylated (activated) forms of ATM, ATR, H2AX, Chk1 and Chk2. Previously, we showed that the FFAA mutation caused DNA replication defects that activated ATR and Chk16 (Fig. 1d). A number of the primary WT/FFAA cells were positive for X-gal staining, which was the indicator for cellular senescence (Fig. 1d). Near-diploid tumour cells of limited expansion were positive for phosphorylated ATR-, H2AX- and Chk1 (Fig. 1d), but unexpectedly, were negative for phospho-ATM and -Chk2. X-gal staining also showed that tumour cells in the limitedly expanding clones were predominantly senescent (Fig. 1d). In contrast, very few of the aneuploid tumour cells of unlimited expansion showed activation of these DNA damage response proteins, and they did not undergo cellular senescence (Fig. 1d).

Polyploid tumour cells alter DNA damage response pathways. To search for the molecular mechanisms driving how the aneuploid cancer cells overcome DNA replication stresses and abrogate DNA damage responses and cellular senescence, we used mRNA micro-array analysis to compare genome-wide gene expression proling of primary WT/FFAA MEF cells (control) with that of the tumour cells having limited or unlimited proliferation potential. A total of 79 genes were upregulated and 121 downregulated in the senescent tumour cells, compared with the normal control. Ingenuity pathway analysis indicated cancer as the most relevant disease, conrming that these cells were indeed tumour-related. Furthermore, CXCL2-centered cytokine signalling was the only signicantly upregulated gene network in the limitedly expanding tumour cells (Fig. 2). The most upregulated gene CXCL2 (15-fold) is a ligand of CXCR2, which has a key role in reinforcement of ATM/ATR-mediated onco-genesis-induced senescence16. In contrast to the senescent tumour cells, expression of CXCL2 and other related chemokine-signalling genes was signicantly downregulated in the expanding aneuploid cancer cells (Fig. 2). This suggested that the aneuploid cancer cells may downregulate CXCL2-mediated signalling to attenuate ATR-mediated oncogenesis-induced senescence.

Approximately 1,400 genes were signicantly upregulated in the aneuploid cancer cells compared with the control cells. The highly upregulated gene networks were DNA replication, chromosome organisation and DNA repair genes for the HR, NHEJ, BER and MMR pathways (Fig. 2 and Supplementary Table S1). In addition to DNA repair genes, the cell-cycle checkpoint genes BRCA1 and p19arf were signicantly upregulated. This increased expression was

NATURE COMMUNICATIONS | 3:815 | DOI: 10.1038/ncomms1825 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1825

ARTICLE

a

Grafting-tumour

Senescence

Primary MEF

Colony formation

Colony expansion

Unlimited proliferation (23%)

Limited proliferation (77%)

b Unlimited proliferation

Limited proliferation

Primary MEF

c

80 Primary MEF

Metaphase cell (%)

60

60

Relative cell number

50

40

40

20

30

20

0

40

4149

5069

7079

10

3539

80 8189

9099

100109

110119

120129

130139

140149

150159

>159

0

0 5 10 15 20 Time (days)

50 Limited proliferating clones

Metaphase cell (%)

40

d Unlimited

proliferation

Limited

proliferation

Primary MEF

30

20

10

0

pATR

X-gal-H2AX

pChk1

3539

40

4149

5069

7079

80 8189

9099

100109

110119

120129

130139

140149

150159

>159

50

Unlimited proliferating clones

Metaphase cell (%)

40

30

20

10

0

7079

40

4149

5069

6069

9099

100109

110119

120129

130139

140149

150159

>159

80

8189

Chromosome number

Figure 1 | Selection of near-polyploid aneuploid cells clonal expansion under oncogenesis-induced DNA replication stresses. (a) Cartoon of expansion method for spontaneously formed colonies of WT/FFAA MEFs. Cells from a single colony were cultured in a 48-well plate. Once the cell culture reached 90% conuence, the cells were further expanded in 24-, 12- and 6-well plates, or 10-cm dishes. Colonies that expanded for at least 10 passages without a decrease in proliferation rate were considered unlimited expansion clones, and were conrmed to form cancers in NOD/SCID mice. (b) Growth rates of diploid MEF cells, limited and unlimited expansion colonies. Values are means s.d. of ve independent clones. (c) Chromosome numbers of primary WT/

FFAA MEF cells (P0) and representative clones with limited or unlimited expansion at the 6-well plate stage were counted using metaphase spread. More than 100 mitotic cells were analysed for each spread. (d) Upper panels: primary WT/FFAA MEF cells and tumour cells from limitedly and unlimitedly proliferating clones stained for the phosphorylated forms of ATR, H2AX and Chk1 (red). Nuclei were stained with DAPI (blue). Scale bar equals 25 m.

Lower panels: detection of senescent cells by X-gal staining (blue). Scale bar equals 100 m.

conrmed by quantitative real-time PCR, because the number of mRNAs of BRCA1 and p19arf per gene copy in the aneuploid cancer cells was signicantly higher than in the control diploid cells (Supplementary Fig. S3). BRCA1 is central to the upregulated DNA-repair networks interacting with other DNA-repair proteins such as Rad51 and the Fanconi Anemia proteins, and it is thought to be important for various DSB-repair pathways (Fig. 2)1720. BRCA1 also interacts with the tumour suppressor gene p19arf to stabilise p53 and induce the expression of p21, resulting in cellular senescence2123. Interestingly, we noticed that p21 and other p53 target genes were signicantly downregulated in the aneuploid cancer cells (Fig. 2 and Supplementary Table S2), suggesting that the activation

of BRCA1 and p19arf, and expression of p21 were uncoupled in the aneuploid cancer cells. Therefore, the transcriptome proling analysis suggested the BRCA1- and p19arf-mediated DNA damage response and repair networks, which could facilitate DNA repair but suppress senescence or apoptosis.

Overexpression of BRCA1 and p19arf enhances DNA SSB repair. We next investigated if overexpression of BRCA1 and p19arf was induced by polyploidisation. We used uorescence-activated cell sorting (FACS) to separate the tetraploid cell population from the diploid one, which had similar cell-cycle proling to each other (Supplementary Fig. S4a and 4b). Because the primary tetraploid

NATURE COMMUNICATIONS | 3:815 | DOI: 10.1038/ncomms1825 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1825

Unlimited proliferating clones (near-polyploid cancer cells)

A

Cytokine signalling

p53-pathways

BRCA1-mediated DSBR

Limited proliferating clones (near-diploid senescent cells)

DNA replication

Signaling chromosome organisation

B

O

Primary MEFs (diploid cells)

Cytokine signalling

Figure 2 | Gene expression analysis of near-polyploid aneuploid cancer cells. Microarray gene expression proling analysis of primary WT/FFAA diploid control cells (O), or tumour cells with limited expansion (B) or unlimited expansion (A). Microarray hybridisation data analysed with Ingenuity Pathway analysis identied signicantly relevant disorders and upregulated or downregulated gene networks or pathways, compared with the diploid control. Red circles, upregulated genes; green circles, downregulated genes; solid lines, physical interaction; arrows with dashed lines, induction of gene expression; arrows with solid lines, activation. The degree of the change in expression is qualitatively reected by the darkness of the color. DSBR: DSB repair.

cells are genetically and epigenetically unstable24, it seemed possible that not all tetraploid cells overexpressed BRCA1 and p19arf. A microuidic device was used to capture single cells from the FACS-sorted diploid or tetraploid cell populations, most of which were at the G1/G0 phase (Supplementary Fig. S4). Expression proling of BRCA1 and p19arf was determined at the single-cell level. The diploid WT/FFAA cells, which were previously shown to have increased numbers of DNA SSBs and DSBs6, expressed low levels of BRCA1 or p19arf in a relatively uniform manner (Supplementary Fig. S5), suggesting that DNA damage due to the FFAA FEN1 mutation

is not responsible for the overexpression of BRCA1 and p19arf. In contrast, the tetraploid cells overexpressed BRCA1 or p19arf in a heterogeneous fashion (Supplementary Fig. S5). The heterogeneous overexpression of BRCA1 and p19arf in primary tetraploid cells was conrmed by in situ ViewRNA analysis (Supplementary Fig. S6). Interestingly, all the aneuploid cancer cells uniformly overexpressed both BRCA1 and p19arf (Supplementary Fig. S6). It seems possible, therefore, that tetraploidy could result in the heterogeneous induction of BRCA1 and/or p19arf, and that the cells that overexpress both BRCA1 and p19arf are selected for during clonal expansion.

NATURE COMMUNICATIONS | 3:815 | DOI: 10.1038/ncomms1825 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1825

ARTICLE

a

b

5 20 nt

32 nt 40 nt

80 nt

3

53 5

32 nt 40 nt

80 nt

5

3

3 5

Diploid M

Polyploid

Diploid

M

Polyploid

Time

Time

80

60

40

Ligated

products

Non-ligated

products

80

60

40

Ligated

products

Non-ligated

products

(nt)

(nt)

c

d

5 20 nt

53 5

32 nt 40 nt

80 nt

3

32 nt 40 nt

80 nt

3

53 5

Ig G -Brca1 -p19arf -p19arf

Time (min)

10 20 60 10 20 60 10 20 60

80

60

40

(nt)

Ligated

products

Non-ligated

products

Ig G -Brca1

M

M

10 20 60 10 20 60 10 20 60

Time (min)

80

60

40

Ligated

products

Non-ligated

products

(nt)

e siRNA control siRNA Brca1 siRNA p19arf

4.5

-H2AX foci per nucleus

siRNA Control Brca1 p19arf

4

3.5

3

2.5

-H2Ax

2

1.5

1

0.5

0

Figure 3 | BRCA1 and p19arf are crucial for DNA SSB repair and reduction of DNA DSBs. A gapped DNA substrate (a) or a gapped DNA ap substrate (b) (upper panels) was used to assay the NEs of primary diploid MEFs (P0) or the aneuploid cancer cells. A quantity 2 g of each NE was incubated with the indicated DNA substrates (1 pmol). All reactions were carried under conditions described in methods for 10, 20 and 60 min, then analysed by 15% denaturing PAGE. (c,d) Reactions to assess SSB repair using NEs from the aneuploid cancer cells were carried out as described for (a) and (b), except that mouse IgG (400 ng) against Brca1 or p19arf was added to block their function. Non-specic mouse IgG (400 ng) was added for the controls (e) Aneuploid cancer cells with siRNA-mediated knockdown of BRCA1 or p19arf were stained for -H2AX (red) to evaluate DNA damage. Nuclei were stained with DAPI (blue). Left panel: -H2AX staining. Right panel: quantication of -H2AX foci per nuclei. The total number of -H2AX foci was normalised with the total number of nuclei. Values are mean s.d., n = 3 independent assays. Scale bar equals 25 m.

Next, we investigated the role of overexpression of BRCA1 and p19arf in coping with DNA replication stresses. One possible mechanism is that it promotes the repair of DNA SSBs that arise due to FEN1 FFAA mutation as well as oncogenesis-induced hyper-DNA replication. To evaluate if the aneuploid cancer cells that overexpressed both BRCA1 and p19arf had a greater capacity for repairing DNA SSBs than did the diploid MEFs, nuclear extracts (NEs) were prepared from both cell types and assayed the DNA SSB repair efficiencies, using two gapped DNA substrates representing DNA SSB intermediate structures that occur during Okazaki fragment maturation or long-patch BER (Fig. 3a,b). NEs from the aneuploid cancer cells generated considerably more fully repaired products than did NEs from the primary diploid MEFs (Fig. 3a,b). However, adding BRCA1 or p19arf antibodies to NEs from the aneuploid cancer cells reduced the in vitro SSB repair efficiency by more than 90% (Fig. 3c,d). It indicated that that BRCA1 and p19arf have important roles in stimulating DNA SSB repair in these cells. To further elucidate how BRCA1 and p19arf contribute to SSB repair, the eect of BRCA1 and p19arf on gap lling mediated by Pol and Pol, which are essential steps during DNA SSB repair1,5, was analysed. We found that recombinant human BRCA1 could slightly (approximately two-fold) stimulate human Pol and

Pol to incorporate 32P-dCTP into a gapped DNA duplex, whereas recombinant human p14arf protein, the mouse p19arf homologue, greatly enhanced the gap-lling activity (Supplementary Fig S7a, b). In addition, both BRCA1 and p14arf enhanced FEN1-mediated ap cleavage (Supplementary Fig. S8), which occurs during Okazaki fragment maturation, and can also occur during LB-BER, DNA SSB repair and NHEJ4,5,2527. Small interfering RNA (siRNA) to knockdown BRCA1 or p19arf expression in the aneuploid cancer cells (Supplementary Fig. S9a,b) showed that the number of H2AX-foci per nuclei was greatly increased in the aneuploid cancer cells treated with BRCA1 or p19arf siRNA (Fig. 3e). Knockdown of BRCA1 or p19arf also arrested the growth of the aneuploid cancer cells (Supplementary Fig. S9c).

The NHEJ pathway is stimulated in polyploid cancer cells. We noticed that the levels of some DNA repair genes that are involved in the NHEJ pathway, including DNA-PK, WRN and XRCC47, were increased in the near-polyploid aneuploid cancer cells (Fig. 2). To determine if the aneuploid cancer cells had an enhanced NHEJ activity, we assayed the NHEJ activity of NEs from normal MEFs or aneuploid cancer cells using a synthetic oligo-based substrate with non-compatible 3 ends28. The NHEJ activity of the NEs from the

NATURE COMMUNICATIONS | 3:815 | DOI: 10.1038/ncomms1825 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1825

a

b

32P 5 3 5

82 nt 80 nt

32P

80 nt 82 nt

*

WT/FFAA Diploid (2n = 38XY)

*

5

3 5

2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

1

17 18 19 X Y

Diploid

Polyploid

M

200

(nt)

Ligated

products

Substrates

120

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 X Y

80

60

WT/FFAA Polyploid (63)

1 2 3 4 5 6 7 8 9 10 1116

M M

c

32P

80 nt 82 nt

* 5

3 5

82 nt 80 nt

32P

5

3 5

*

12 13 14 15

17 18 19 X Y

-p19arf 10 20 60 10 20 60 10 20 60

Ig G -Brca1

Time (min)

M

1

4

1

X

X 186 7 8 9 10

11 12 13 14 15

16

MM

17 18 19 X Y

1 2 3 4 5

200

(nt)

Ligated

products

120

80

Substrates

60

Figure 4 | Near-polyploid aneuploid cancer cells promote the NHEJ activity and cause chromosome translocations. (a) NHEJ activity of NEs from the diploid normal and aneuploid cancer cells. A 32P-labeled DNA duplex with a two nucleotide 3-end overhang was prepared (upper panel) and incubated with the indicated NEs (4 g each) under conditions described in methods for 10, 20 and 60 min, and analysed by 6% denaturing PAGE. (b) GTG (upper panels) and SKY (lower panels) karyotypes of chromosomes in primary WT/FFAA MEFs or aneuploid cancer cells (clone C21). Structural chromosome aberrations, including two dicentric chromosomes composed of material from chromosome 1 and chromosome X (Dic(1H1;XA1)), and material from chromosome 4 and chromosome 18 (Dic(4C7, 18E4)), two inversions of chromosome 5 (In(5B-F)), and two additional micro-chromosomes, are indicated by arrows. (c) Reactions to assess NHEJ activities using NEs from the aneuploid cancer cells were carried out as described in (a), except mouse IgG against BRCA1 or p19arf was added to block their function. Non-specic mouse IgG was used for control reactions.

aneuploid cancer cell was increased by more than 20-fold compared to the NEs from primary MEFs, which had little NHEJ activity on these non-compatible DNA ends (Fig. 4a). Although the increase in NHEJ activity should allow the tumour cells to repair DSBs and suppress DSB-induced cellular senescence or apoptosis, at the same time, it could also lead to misjoining the one-ended DNA DSBs with other DSBs. In support of this hypothesis, all of the WT/FFAA aneuploid cancer cells had chromosome translocations, but none in the diploid WT/FFAA cells (Fig. 4b). Adding BRCA1 or p19arf antibodies to the NEs also inhibited the NHEJ activity on non-compatible DNA ends, albeit at a moderate level (Fig. 4c).

Epigenetic silencing of p21 abrogates cellular senescence. BRCA1 and p19arf, on the other hand, may also activate p53, which induces the expression of downstream eectors such as p21, to limit the proliferation of aneuploid cancer cells2123,29. Activation of p53 in WT/ FFAA aneuploid cancer cells was conrmed by western blot analysis showing that the levels of nuclear p53 and phosphorylated-p53 were considerably elevated (Fig. 5a). This result led us to question how

the cancer cells employed increased levels of BRCA1 and p19arf to reduce DNA replication stresses, yet avoided BRCA1-, p19arf- and p53-mediated cellular senescence and apoptosis. Our microarray data indicated that p21 was downregulated (Fig. 2). Because p21 is a key eector gene downstream of the BRCA1p19arfp53 pathway3032, this event might be crucial for enabling the tumour cells to escape BRCA1p19arfp53-mediated cellular senescence. Western blot analysis conrmed that the activation of p53 and expression of p21 were uncoupled in the aneuploid cancer cells (Fig. 5a). We searched for DNA mutations in the coding region and boundary regions between the exons and introns of p53. However, we did not detect any mutations or deletions of the p53 gene in the aneuploid cancer cell lines or in lung tumours, suggesting that other mechanisms, such as epigenetic changes, were allowing the cells to bypass the p53-mediated barriers33. Supporting this hypothesis, we found that DNA methylation occurred in CpG islands in the promoter regions of the p53 target genes p21, PERP and RPRM (Fig. 5b and Supplementary Fig. S10). To further test our hypothesis that DNA methylation contributed to the downregulation of p21,

NATURE COMMUNICATIONS | 3:815 | DOI: 10.1038/ncomms1825 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1825

ARTICLE

a

b

p21 (Cdkn1a)

CpG islands

c

100

DMSO

Polyploid

AzadC

Normalized p21mRNA

levels

28,820,200

28,821,000

3.02.00.0

3.02.00.0

3.02.00.0

28,822,000

P-p53

p53

p21

Diploid

50

50

kd

10

25 15

40

Log 2ratio (MIRA/INPUT)

Polyploid cells

1 Diploid Polyploid

Diploid cells

-actin

d

e

AzadC

f

siRNA control

g

120 0 100 500 ng ml1

120 siRNA p21

300 24 48 72h

Time:

GFP GFP-p21

AzadC

+ +

p21

-actin

Diploid Polyploid

100

100

25 15

40

kd

Relative growth rate (%)

Relative cell growth

60

Relative growth rate (%)

250

80

80

200

60

150

40

40

100

20

20

50

0 Diploid Polyploid

0 Control AzadC 0

Figure 5 | DNA methylation silencing of p21 uncouples p53 activation and p21 expression. (a) Western blot analyses of levels of the nuclear formof p53, phosphorylated p53 (P-p53), and p21 in the normal diploid and aneuploid cancer cells. (b) DNA methylation status in the promoter regions of p21 (Cdkn1a) in aneuploid cancer cells and normal diploid cells. DNA methylation status is indicated by the Log2 value of the hybridisation signal from the enriched methylated genomic DNA (MIRA) normalised with that from the input genomic DNA (INPUT). Black bars, hybridisation signals; grey bar, location of CpG islands. The DNA methylation proling data accession number: GSE18815. (c,d) Near-polyploid aneuploid cancer or diploid normal cells (n = 3 for each type) were treated with 5 AzadC (500 ng ml 1) for 48 h. The mRNA (c) and protein levels (d) of p21 in the 5 AzadC-treated or untreated cells were analysed by quantitative PCR and western blot analysis, respectively. (e) Aneuploid cancer or diploid normal cells were treated with 5 AzadC (0, 100, or 500 ng ml 1; 48 h). Live cells were counted, and the relative growth rate calculated by dividing the number of cells by the number of the untreated cells (mock). Values are means s.d. of three independent assays. (f) Aneuploid cancer cells were pretreated with control or p21 siRNA then cultured in the absence or presence of 500 ng ml 1 5 AzadC for 48 h. The cell number was counted and the relative growth rate was calculated by dividing the number of cells under varying treatments by the number of cells treated with control siRNA oligos only. (g) Near-polyploid aneuploid cancer cells expressing GFP or GFP-p21 were seeded in 96-well plates (4,000 each well). Cell growth was monitored by CellTiter 96 Aqueous One Solution Cell Proliferation Assay kit every 24 h for 72 h. The relative cell number was calculated by comparison of the absorbance value at 570 nm (A570) of each well. The relative cell number of the wells with the cells expressing GFP only at 24 h was set as 100. Values are mean s.d. of n = 5 assays on two independent tetraploid cell lines.

we treated the aneuploid cancer cells and diploid normal cells with the DNA methylation inhibitor 5-aza-2-deoxycytidine (5-AzadC; 500 ng ml 1), and this treatment restored p21 expression in the aneuploid cancer cells (Fig. 5c,d). Consistently, 5 AzadC suppressed the proliferation of aneuploid cancer cells, but not of the primary MEFs (Fig. 5e). To further test if DNA methylation silencing of p21 was indeed a key factor in the ability of the aneuploid cancer cells to evade growth arrest, we pre-treated the cells with siRNA specic to p21 and then cultured them in the absence or presence of 5 AzadC. We found that 5 AzadC signicantly inhibited the growth of the aneuploid cells treated with a control siRNA, but this eect was largely reversed in cells treated with p21 siRNA and 5 AzadC (Fig. 5f). Furthermore, we showed that the expression of exogenous p21 in the aneuploid cancer cells resulted in an obvious growth inhibition (Fig. 5g).

Aneuploid cancer cells were also found to be resistant to exogenous cytotoxic insults, including tumour necrosis factor (TNF) (Fig. 6a), which induces p53-depedent senescence and apoptosis34,35. However, addition of 5 AzadC greatly sensitised the aneuploid cancer cells to TNF-mediated cytotoxic eects (Fig. 6b). This nding provided additional evidence that the epigenetic silencing of p21 and other p53-mediated senescence or apoptosis genes is a

critical mechanism that allows the aneuploid cancer cells to escape p53-mediated proliferation barriers. To investigate if this result was also true in vivo, 5 AzadC was administered to NOD/SCID mice (500 ng g 1 body weight for 5 days) that were injected with WT/FFAA aneuploid cancer cells. Treatments with 5 AzadC signicantly decreased the size of aneuploid-associated tumours compared with untreated NOD/SCID controls (Fig. 6c).

Polyploidisation promotes cancer progression in vivo. We previously showed that the WT/FFAA mice developed high frequency of lung adenoma and adenocarcinoma6. It has been suggested that mutations in oncogenes including EGFR and Kras can initiate lung tumours3638. Supporting this hypothesis, we found that of the ten WT/FFAA lung tumours we examined, two tumours carried the L705S or K738N mutations in the EGFR gene, and another two tumours had the G12D mutation in the Kras gene, which frequently occurs in human and mouse lung cancers3638. However, hyper-DNA replication in the newly initiated tumour cells due to oncogene mutations can cause DNA replication stresses, especially if the cells already carry mutations in DNA replication genes, such as FEN1. To investigate if the formation of polyploidy and subsequent aneuploidy contributed to overcoming such replication

NATURE COMMUNICATIONS | 3:815 | DOI: 10.1038/ncomms1825 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1825

a

Recombinant TNF:

b c

2.5

120 0 10 40

ng ml1

120 AzadC

DMSO

Mean

100

100

2

Relative growth rate (%)

Relative growth rate (%)

80

80

1.5

60

60

Tumour weight (g)

1

40

40

0.5

20

20

0

0 Diploid Polyploid

0 Control TNF

AzadC 0 100 200 500 (g kg1)

Figure 6 | The DNA methylation inhibitor 5 AzadC restores p53-mediated growth arrest of near-polyploid aneuploid cancer cells. (a) Aneuploid cancer cells were resistant to TNF-mediated growth arrest. Proliferation of aneuploid cancer or normal diploid cells in the presence of indicated concentrations of recombinant mouse TNF was assessed. (b) Growth of aneuploid cancer cells treated with 40 ng ml 1 TNF in the absence or presence of 100 ng ml 1 5 AzadC for 72 h. In panels (a) and (b), live cells were counted and the relative growth rate calculated by dividing the number of cells by the number of untreated cells (mock). Values are mean s.d. of three independent assays. (c) 5 AzadC inhibits aneuploid cancer growth in vivo. Tumour weights of NOD/SCID mice that rst received 5106 aneuploid tumour cells (from clone C21), and were then treated with PBS buffer (control) or 5 AzadC (100, 200 or 500 ng g 1 body weight; n = 15 for the control group, n = 8 for each treatment group) for 5 consecutive days.

stresses in vivo, we rst employed a dual-colour uorescence in situ hybridisation (FISH) to detect aneuploid cells carrying three or more copies of both chromosome 2 and chromosome 8 in tissue sections of lung adenoma (n = 8) and adenocarcinoma (n = 10) from WT/FFAA mice. The adenocarcinoma but not adenoma had a signicant amount of such aneuploid cells (Fig. 7a and Supplementary Fig. S11a). To further verify that polyploidisation was associated with tumour expansion and progression in vivo, we isolated single cells from a pool of small (12 mm) or large lung tumours (45 mm), which are predominantly lung adenoma or adenocarcinoma, respectively, according to previous histological analysis6,38. FACS-analysis revealed that the adenoma cells were diploid or near diploidy, whereas the adenocarcinoma cells were near-polyploid aneuploidy (Fig. 7b and Supplementary Table S3). We noticed the G1/G0 peak of adenocarcinoma were relatively broader than that of the normal lung or adenoma, indicating that the population might contain heterogeneous aneuploid cells. In addition, aneuploid adenocarcinoma cells had a higher S/G2 cells than the diploid normal or adenoma cells (Fig. 7b and Supplementary Table S3). Moreover, lung adenocarcinoma had signicantly less TUNEL (terminal deoxynucleotidyl transferase dUTP nick end labelling)-positive cells than lung adenoma (Fig. 7c,d and Supplementary Fig. S11b). The lung adenocarcinoma also had less X-gal-positive cells, which is indicative of cellular senescence39, than the lung adenoma (Fig. 7e,f and Supplementary Fig. S11c).

Disscussion

DNA SSBs arise and persist due to genetic defects in DNA replication and repair genes, including FEN1. The FFAA FEN1 mutation disrupting the FEN1/PCNA interaction causes defects in DNA replication and repair, and results in DNA SSBs in the genome6. Consistently, mutations at the R192 residue, including R192Q, a FEN1 SNP present in the human genome (http://www.ncbi.nlm.nih. gov/projects/SNP/), that also disrupt the FEN1/PCNA interaction, result in DNA strand breaks (Guo et al.2). In addition, we have identied a group of FEN1 mutations, which eliminate the exonuclease activity of FEN1 and are present in human cancers38. More recently, we have shown that mouse cells carrying the E160D mutation that models this class of nuclease-decient mutants, display BER defects and accumulate DNA SSBs in response to the DNA mutagen methylnitrosourea40. It is known that unre-paired SSBs may collapse DNA replication forks, resulting in DNA replication stresses and DSBs1,5,40. Meanwhile, oncogenesis also induces DNA replication stresses, which activate ATM/

ATR and p53, and lead to cellular senescence. These molecular mechanisms have been thought to prevent early tumour progression11,12,41. Therefore, the initiated tumour cells must adapt to or reduce the DNA replication stresses to progress to malignancy. In support of this hypothesis, recent studies indicate that carcinoma has a considerably lower DNA damage response index than the corresponding precancerous lesion11,42,43. Our current studies consistently revealed that in unrestrictedly proliferating tumour cells, ATR-mediated DNA damage responses and cellular senescence were greatly reduced. In addition, WT/FFAA lung adenoma cells, but not adenocarcinoa, were associated with high levels of senescence and apoptosis.

We further suggest the formation of tetraploidy and aneuploidy, which is a hallmark of human cancer24,44,45, is a key molecular mechanism that enables tumour cells to overcome DNA replication stresses and bypass the p53-mediated cell death pathways to acquire unlimited proliferation potential (Fig. 8). In our current study, a very striking cellular phenotype observed in the in vitro expansion assays was that all unlimitedly proliferating WT/FFAA cells were near-polyploid aneuploid. In addition, lung adenoma from WT/FFAA mice was shown to be largely diploid, but lung adenocarcinoma was largely polyploid. These observations suggest that polyploidisation may positively aect cancer progression. Our single-cell gene expression proling studies suggest that polyploidisation produces populations of cells with a wide spectrum of DNA damage response and repair gene networks. This would increase the chance that a subpopulation of the aneuploid tumour cells adapt or overcome replication stresses and bypass the p53-mediated cell death or senescence pathways, thereby evolving into malignant cancer cells. This hypothesis is consistent with the current view that cancer is a micro-evolutionary disease46.

We have shown that WT/FFAA aneuploid cancer cells with upregulated BRCA1- and p19arf-mediated DNA damage response and repair gene networks are selected for clonal expansion, suggesting an important role for these two proteins in coping with DNA replication stresses. BRCA1 interacts with Rad51 to promote DSB repair via the HR pathway8,1720. In addition, BRCA1 has been shown to promote NHEJ activity19. The p19arf has been shown to stabilise p53, which mediates DNA repair, cellular senescence and apoptosis21,47. Here, our data suggest new, critical roles for BRCA1 and p19arf in SSB repair and the reduction of DNA replication stresses. Our in vitro assays using puried recombinant proteins demonstrated that p19arf strongly stimulated both Pol- or Pol-mediated gap lling and FEN1-mediated ap cleavage, and BRCA1

NATURE COMMUNICATIONS | 3:815 | DOI: 10.1038/ncomms1825 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1825

ARTICLE

a

b

c

d

100

100 Adenoma

100 Adenocarcinoma

DNA index: 1.1

DNA index: 1.7

80

60

Adenoma

Adenocarcinoma

counts

80

60

40

20

0

Adenoma Adenocarcinom a

14.22%

Adenoma

Adenocarcinoma

Counts

Counts

40

20

0

80

counts

100

80

60

40

20

0

60

1.21%

40

20

20 40

DNA content

0 0 10 1,000 10,000

FITC

0

60 (K) 100

e

f

Adenoma

Adenocarcinoma

140

-galactosidase activity

(unit per mg cell extracts)

120

100

80

60

40

20

0 Adenoma

Adenocarcinoma

Figure 7 | Malignant adenocarcinoma in WT/FFAA mice was associated with near-polyploid aneuploidy and reduced apoptosis. (a) In vivo levels of aneuploidy of lung adenoma (n = 8) and adenocarcinoma (n = 10) were analysed with dual-colour interphase FISH (probes 2qA1 (red, chromosome 2) and 8qA1 (green, chromosome 8)). Aneuploid cells with three or more copies of chromosomes 2 and 8 are indicated by white arrows. Scale bar equals 20 m.

(b) FACS analysis of cells isolated from normal lung, adenoma (a pool of eight tumours of 12 mm diameters), and adenomcarcinoma (tumours of 35 mm diameters), following a previously published protocol51. The data was processed using the FlowJo software. DNA index was determined by comparingthe position of the G1/G0 peak of adenoma or adenocarcinoma with that of the normal lung cells (diploid). (c) TUNEL staining (red) of apoptotic cellsin tissue sections of the lung adenoma (n = 8) or adenocarcinoma (n = 10). Nuclei were stained with DAPI (blue). Pink spots were counted as TUNEL-positive nuclei. Scale bar equals 25 m. (d) Tumour cells isolated from normal lung (negative control), adenoma (a poll of eight small tumours with diameters of 1-2 mm) or adenocarcinoma (a pool of three large tumours with diameters of 35 mm) were stained with APO-DIRECT kit (TUNEL assay). Cells with FITC intensities that were greater than 500 were considered as apoptotic cells (Supplementary Fig. S12). Dotted lines: Cutoff for quantication of apoptotic cells. (e) X-gal staining (blue) of senescent cells in the tissue sections of lung adenoma (n = 8) or adenocarcinoma (n = 6). The cytoplasm was counterstained with eosin solution (red). Tumours with at least ten blue spots, as shown for adenoma, were considered X-gal-positive. Scale bar equals 50 m. (f) Mammalian -galactosidase activity in the cell extracts from the lung adenoma or adenocarcinoma were assayed by the mammalian -galactosidase activity assay kit (Thermo Scientic). One unit -galactosidase activity was arbitrarily dened as an increase of 0.001 of absorbance405 nm (A405) per hour. Values are mean s.d of assays using three independent cell extracts.

also enhanced ap cleavage by FEN1. Because gap lling and ap cleavage are important steps in both the DNA SSB repair and NHEJ repair pathways1,5,26, BRCA1 and p19arf may promote DNA SSB and NHEJ activity. A decrease in the DNA SSBs would reduce the frequency of collapsed DNA replication forks, thereby attenuating DNA replication stress response. The enhanced NHEJ activity should facilitate the repair of DSBs to suppress DNA damage response pathways that lead to cellular senescence. However, this may also increase improper repair of many one-ended DNA DSBs leading to chromosomal translocations and genome instabilities.

Both BRCA1 and p19arf can activate p53, which induces the expression of p21 and other genes to arrest the cell cycle and trigger senescence or apoptosis2123,29. Therefore, inactivation of these p53-mediated pathways may be critical for cells to progress towards malignancy. Mutations in p53 occur in approximately 50% of human cancers, and it is proposed to be a key mechanism enabling tumour

cells to escape p53-mediated cellular senescence or apoptosis29,48,49. Here we suggest that DNA methylation at the promoter regions of p21 and other p53 target genes can also act to silence the p53 pathways, and results in a unique response to DNA replication stress, allowing the aneuploid cancer cells to promote DNA repair to prevent breakage-mediated cell death, but escaped cellular senescence and apoptosis (Fig. 8). In addition, we suggest that DNA methylation-induced silencing of the p53 pathways allows the aneuploid cancer cells to become resistant to exogenous cytotoxic insults, such as TNF, which is secreted by immune cells to induce senescence and apoptosis via the p53 pathways50. This would in turn, further promote the in vivo progression of cells towards malignancy. Taken together, the results from our current study indicate that aneuploid cancer cells overcome DNA replication stresses by increasing DNA repair activity and escape senescence and apoptosis via epigenetic reprogramming of the p53-mediated senescence and apoptosis pathways rather than through p53 mutation.

NATURE COMMUNICATIONS | 3:815 | DOI: 10.1038/ncomms1825 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1825

Tetraploidisation

Hyper-

proliferation

DNA replication

stresses

Initiated tumour cells

(Fen1 mutations)

Enhanced DNA repair

Epigenetically-silencing of p21

Senescent benign tumour

Polyploid cancer

Figure 8 | A model of tetraploidy formation as a key molecular eventfor overcoming DNA replication stresses and malignant cancer progression. DNA replication stresses, which can occur due to mutations in DNA replication genes such as Fen1, environmental insults or oncogene activation, normally arrest the cell cycle and the senescence and apoptosis pathways serve as barriers to prevent tumour cell progression. At the same time, DNA replication and repair defects frequently result in the formation of tetraploidy and polyploidy. This leads to populations of genetically and epigenetically heterogeneous cells represented by different colours (white, red, pink, green, or blue). We observed subpopulations of polyploid cells with overexpression of BRCA1, p19arf and other DNA repair genes that reduce DNA replication stresses, and DNA methylation silencing of p53 downstream targets to escape senescence and apoptosis, and proposed these cells undergo clonal expansion (red cells).

Methods

Generation of limited and unlimited expansion transformed cells. Spontaneously transformed cells were produced as previously described6. Briey, 5105 mouse primary MEFs (P1, E13.5) were seeded onto a 10-cm dish and cultured in Dulbeccos modied Eagles medium supplemented with 10% fetal bovine serum for 1525 days. Cells from each colony were picked, transferred and cultured in a 48-well plate. Once the culture reached 90% conuence, cell numbers were counted and all cells were transferred to a 24-well plate. Likewise, the cells were passedto a 12-well or 6-well plate, or a 10-cm dish. If cells from a clone continuously expanded with similar or higher proliferation rates aer six additional passages in the 10-cm dish, we considered them to be an unlimited expansion clone. Those that did not continue to proliferate under these conditions were considered to be limited expansion clones.

Tumour grafting in NOD/SCID mice. Aneuploid cancer cells (5107) were subcutaneously injected into NOD/SCID mice (2 months old, n = 3 for each cell line). For the 5 AzadC experiments, NOD/SCID mice that were injected with aneuploid cancer cells were also treated with 5 AzadC daily for 5 days (0, 100, 200, 500 ng g 1 body weight, 100200 l injection volume). All mice were observed for 2 months. The mice were then killed, and subcutaneous tumours were dissected and weighed. All protocols involving animals were approved by the Research Animal Care Committee of City of Hope, in compliance with the Public Health Service Policy of the United States.

Recombinant protein. Recombinant human FEN1, Pol and Pol were expressed and puried as previously described6,52. Puried recombinant human BRCA1 was purchased from Active Motif. To express and purify human p14arf (the human homologue of mouse p19arf), a pcDNA-myc-ARF plasmid, which encodes a c-myc-tagged human p14arf53, was transfected into 293T cells. Aer 48 h additional culturing, the cells were harvested and lysed. The c-myc-tagged p14arf was puried by the affinity purication kit for the c-myc tagged protein, according to the manufactures instruction (MBL International). The eluted c-myc-tagged p14arf was examined by SDSPAGE and veried by western blotting analysis (Supplementary Fig. S13), using an antibody against human p14arf (Santa Cruz Biotchnologies).

In vitro DNA SSB repair and NHEJ assays. SSB repair on the gapped DNA substrate with or without a DNARNA ap was assayed as previously described6,54. Briey, NEs were prepared and mixed with DNA substrates (1 pmol) in reaction buer A (50 mM Hepes-KOH, pH 7.5, 45 mM KCl, 5 mM MgCl2, 1 mM dithiothreitol, 0.1 mM EDTA, 2 mM ATP, 200 units creatine-phosphokinase, 0.5 mM NAD and 5 mM phosphocreatine). Each reaction (15 l) also contained 5 Ci [-32P] dCTP and 50 M each of dATP, dGTP and dTTP. NHEJ was assayed as previously described19,28. A 3 end-32P-labeled oligo-based DNA duplex was prepared. There was a two-nucleotide (GG) overhang at the non-labeled 3-end of the DNA substrate to resemble non-compatible DNA end joining28. NEs were incubated with DNA substrates (1 pmol) in the reaction buer (50 mM triethanolamineHCl, pH

7.5, 5 mM Mg(OAc)2, 80 mM potassium acetate, 2 mM ATP, 1 mM dithiothreitol and 100 g ml 1 BSA) containing 50 M each of the four deoxyribonucleotides. SSB repair or NHEJ reactions were carried out for the indicated times at 37 C, and the product was analysed with 15 or 6% denaturing PAGE and autography.

Metaphase spread preparation and analysis. Cells treated with colcemid (0.1 g ml 1, 4 h) were harvested, washed with PBS buer and incubated with hypotonic solution (75 mM KCl). Cells were then xed with Carnoys solution and dropped onto a glass slide, which was then baked (65 C, overnight), stained with a Giemsa solution, scanned under a microscope (Olympus AX70) for mitotic cells and imaged using Image Pro 6.3 (Media Cybernetics, Inc.). The chromosome number of each mitotic cell was analysed and scored using Image Pro 6.3. Typically, 100200 mitotic cells were analysed.

Immunouorescence staining. To stain for phospho-Chk1, phospho-ATR or H2AX, MEFs grown on glass cover slides were xed with 4% paraformaldehyde, permeabilised with 0.1%. Triton X100, blocked with Image iT FX signal enhancer (Invitrogen) and incubated (1.5 h, room temperature) with the indicated primary antibodies. Antibodies against phospho-Chk1 (S345; 1:400) and H2AX (1:400) were from Abcam, and the antibodies against phospho-ATR (1:200) and phospho-ATM (1:200) were from Cell Signaling Technology. Cells were then washed with PBS buer and incubated (1 h, room temperature) with corresponding secondary antibodies (1:200, Invitrogen). Slides were washed with PBS buer, counter stained with DAPI and analysed with a uorescence microscope (Olympus AX70).

Detection of senescent cells. Senescent cells were detected by staining for X-gal-positive cells as previously described39. The X-gal staining kit was from Cell Signaling. Briey, cells cultured in 12- or 6-well plates were xed (10% formalin in PBS) for 10 min. Aer being washed with PBS buer, xed cells were incubated (overnight, 37 C) with the X-gal staining solution. Cells were then analysed with an inverted light microscope (IX81, Olympus).

mRNA microarray analysis. For microarray analysis, total RNA was extracted from parental MEFs (WT/FFAA, P1, E13.5, diploid normal control) and cells from clones with limited or unlimited expansion (n = 3 independent lines for each group) using the Qiagen RNaeasy kit (Qiagen). The Aymetrix GeneChip Mouse Gene 1.0-ST array (Aymetrix, Santa Clara, CA) was used to dene gene expression proles from the samples. Biotinylated single-stranded cDNA was generated using 100 ng total RNA according to the manufacturers protocol. Hybridisation cocktails containing 2.5 g fragmented, end-labeled cDNA were prepared and applied to the GeneChip and incubated for 16 h. Arrays were washed and stained with the GeneChip Fluidics Station 450 using FS450_0007 script. They were then scanned at 5-m resolution using the Aymetrix GCS 3000 7G and GeneChip Operating Soware v. 1.4 to produce CEL intensity les.

Raw intensity measurements of all probe sets were background-corrected, normalised and converted into expression measurements using Aymetrixs Expression Console v1.1.1. The Bioconductor LIMMA package was then used to identify genes dierentially expressed between the tetraploid and diploid samples. Signicant genes were selected using a cut-o of P < 0.05 and fold change of 2. Data were analysed with DAVID (Database for Annotation, Visualization and Integrated Discovery) functional annotation tools and Ingenuity Pathway analysis programto identify statistically signicant diseases, gene networks and pathways that were upregulated or downregulated in transformed cells with limited or unlimited expansion potential and compared witho the diploid control.

DNA methylation analysis. DNA methylation in aneuploid cancer cells and normal diploid cells was analysed with MIRA55. Briey, sonicated genomic DNA (approximately 500 bp) was incubated with the methylated DNA-binding proteins MBD2B and MBD3L1 to pull down methylated DNA. Enriched methylated and input genomic DNA samples were quantitatively amplied by real-time PCR and labeled with Cy3 and Cy5, respectively. The labeled DNA fragments were hybridised to the Nimblegen HG18 CpG Promoter array (Roche Nimblegen). Normalised hybridisation signals were analysed with SignalMap (Nimblegen). Methylated CpG islands shown in SignalMap were then veried by Bisulte DNA sequencing55.

DNA methylation inhibitor treatment. WT/FFAA aneuploid cancer cells or normal diploid cells were treated as indicated with the DNA methylation inhibitor, 5-aza-2-deoxycytidine (5 AzadC)56. Expression of p21 and CXCL2 was determined by quantitative real-time PCR (Taqman gene expression assays, Applied Biosystems) and western blot. The corresponding cell growth rate was measured by cell counting.

Histopathology and staining. To conduct histopathology analysis, xed (10% formalin) or frozen tissue sections were stained with haematoxylin and eosin.

To determine the ploidy status of cells in tumour tissues, dual-colour FISH for chromosomes 2 and 8 was conducted with FISH probes that corresponded to the mouse chromosomes 2qA1 and 8qA1. Slides were analysed with a BioView Imaging system (BioView, UK). To stain DNA strand breaks, paraffin-embedded tissue sections (5 m) were de-waxed with xylene (5 min) and rehydrated in an ice-cold ethanol series (100, 95, 75 and 50%; 2 min each). Sections were permeabilised with

0

NATURE COMMUNICATIONS | 3:815 | DOI: 10.1038/ncomms1825 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1825

ARTICLE

0.1 M sodium citrate and 0.1% Triton X-100 (8 min). Slides were then incubated (1 h, 37 C) with TUNEL staining solution (Roche), washed with PBS and examined with a uorescence microscope (Olympus AX70). To stain senescent cells in lung tumour tissue, frozen tissue sections (5 m) were xed with 1% formalin and stained with X-gal staining solution (Cell Signaling) at 37 C overnight. Cytoplasm was stained with eosin Y solution (Sigma). All slides for each staining were analysed in a double-blind fashion.

References

1. Caldecott, K. W. Mammalian DNA single-strand break repair: an X-ra(y)ted aair. Bioessays 23, 447455 (2001).

2. Guo, Z. et al. Functional regulation of ap endonuclease 1 by methylation-mediated phosphorylation suppression and PCNA interaction. Nat. Chem. Biol. 6, 766773 (2010).

3. Liu, Y., Kao, H. I. & Bambara, R. A. Flap endonuclease 1: a central component of DNA metabolism. Annu. Rev. Biochem. 73, 589615 (2004).

4. Zheng, L. et al. Functional regulation of FEN1 nuclease and its link to cancer. Nucleic Acids Res. 39, 781794 (2011).

5. Schar, P. Spontaneous DNA damage, genome instability, and cancerwhen DNA replication escapes control. Cell 104, 329332 (2001).

6. Zheng, L. et al. Fen1 mutations that specically disrupt its interaction with PCNA cause aneuploidy-associated cancer. Cell Res. 21, 10521067 (2011).

7. Helleday, T., Lo, J., van Gent, D. C. & Engelward, B. P. DNA double-strand break repair: from mechanistic understanding to cancer treatment. DNA Repair (Amst) 6, 923935 (2007).

8. Khanna, K. K. & Jackson, S. P. DNA double-strand breaks: signaling, repair and the cancer connection. Nat. Genet. 27, 247254 (2001).

9. Nussenzweig, A. & Nussenzweig, M. C. Origin of chromosomal translocations in lymphoid cancer. Cell 141, 2738 (2010).

10. Bartkova, J. et al. DNA damage response as a candidate anti-cancer barrier in early human tumorigenesis. Nature 434, 864870 (2005).

11. Bartkova, J. et al. Oncogene-induced senescence is part of the tumorigenesis barrier imposed by DNA damage checkpoints. Nature 444, 633637 (2006).

12. Di Micco, R. et al. Oncogene-induced senescence is a DNA damage response triggered by DNA hyper-replication. Nature 444, 638642 (2006).

13. Zheng, L., Dai, H., Qiu, J., Huang, Q. & Shen, B. Disruption of the FEN-1/ PCNA interaction results in DNA replication defects, pulmonary hypoplasia, pancytopenia, and newborn lethality in mice. Mol. Cell. Biol. 27, 31763186 (2007).

14. Jacobsen, K., Groth, A. & Willumsen, B. M. Ras-inducible immortalized broblasts: focus formation without cell cycle deregulation. Oncogene 21, 30583067 (2002).

15. Zavada, J. & Macpherson, I. Transformation of hamster cell lines in vitro by a hamster sarcoma virus. Nature 225, 2426 (1970).

16. Acosta, J. C. et al. Chemokine signaling via the CXCR2 receptor reinforces senescence. Cell 133, 10061018 (2008).

17. Zhong, Q. et al. Association of BRCA1 with the hRad50-hMre11-p95 complex and the DNA damage response. Science 285, 747750 (1999).

18. Scully, R. et al. Association of BRCA1 with Rad51 in mitotic and meiotic cells. Cell 88, 265275 (1997).

19. Zhong, Q., Boyer, T. G., Chen, P. L. & Lee, W. H. Decient nonhomologous end-joining activity in cell-free extracts from Brca1-null broblasts. Cancer Res. 62, 39663970 (2002).

20. Garcia-Higuera, I. et al. Interaction of the Fanconi anemia proteins and BRCA1 in a common pathway. Mol. Cell. 7, 249262 (2001).

21. Palmero, I., Pantoja, C. & Serrano, M. p19ARF links the tumour suppressor p53 to Ras. Nature 395, 125126 (1998).

22. Pomerantz, J. et al. The Ink4a tumor suppressor gene product, p19Arf, interacts with MDM2 and neutralizes MDM2s inhibition of p53. Cell 92, 713723 (1998).

23. Somasundaram, K. et al. BRCA1 signals ARF-dependent stabilization and coactivation of p53. Oncogene 18, 66056614 (1999).

24. Ganem, N. J., Storchova, Z. & Pellman, D. Tetraploidy, aneuploidy and cancer. Curr. Opin. Genet. Dev. 17, 157162 (2007).

25. Klungland, A. & Lindahl, T. Second pathway for completion of human DNA base excision-repair: reconstitution with puried proteins and requirement for DNase IV (FEN1). EMBO J 16, 33413348 (1997).

26. Wu, X., Wilson, T. E. & Lieber, M. R. A role for FEN-1 in nonhomologous DNA end joining: the order of strand annealing and nucleolytic processing events. Proc. Natl Acad. Sci. USA 96, 13031308 (1999).

27. Zheng, L. & Shen, B. Okazaki fragment maturation: nucleases take centre stage. J. Mol. Cell. Biol. 3, 2330 (2011).28. Ma, Y. et al. A biochemically dened system for mammalian nonhomologous DNA end joining. Mol. Cell. 16, 701713 (2004).

29. Zambetti, G. P. The p53 tumor suppressor pathway and cancer. Protein Rev. 2, 5380 (2005).

30. Deng, C., Zhang, P., Harper, J. W., Elledge, S. J. & Leder, P. Mice lacking p21CIP1/WAF1 undergo normal development, but are defective in G1 checkpoint control. Cell 82, 675684 (1995).

31. Dulic, V. et al. p53-dependent inhibition of cyclin-dependent kinase activities in human broblasts during radiation-induced G1 arrest. Cell 76, 10131023 (1994).

32. Bunz, F. et al. Requirement for p53 and p21 to sustain G2 arrest aer DNA damage. Science 282, 14971501 (1998).

33. Allan, L. A., Duhig, T., Read, M. & Fried, M. The p21(WAF1/CIP1) promoter is methylated in Rat-1 cells: stable restoration of p53-dependent p21(WAF1/ CIP1) expression aer transfection of a genomic clone containing the p21(WAF1/CIP1) gene. Mol. Cell. Biol. 20, 12911298 (2000).

34. Donato, N. J. & Perez, M. Tumor necrosis factor-induced apoptosis stimulates p53 accumulation and p21WAF1 proteolysis in ME-180 cells. J. Biol. Chem. 273, 50675072 (1998).

35. Rokhlin, O. W. et al. p53 is involved in tumor necrosis factor-alpha-induced apoptosis in the human prostatic carcinoma cell line LNCaP. Oncogene 19, 19591968 (2000).

36. Paez, J. G. et al. EGFR mutations in lung cancer: correlation with clinical response to getinib therapy. Science 304, 14971500 (2004).

37. Slebos, R. J. et al. K-ras oncogene activation as a prognostic marker in adenocarcinoma of the lung. N. Engl. J. Med. 323, 561565 (1990).

38. Zheng, L. et al. Fen1 mutations result in autoimmunity, chronic inammation and cancers. Nat. Med. 13, 812819 (2007).

39. Dimri, G. P. et al. A biomarker that identies senescent human cells in culture and in aging skin in vivo. Proc. Natl Acad. Sci. USA 92, 93639367 (1995).40. Xu, H. et al. Chemical-induced cancer incidence and underlying mechanisms in Fen1 mutant mice. Oncogene 30, 10721081 (2011).

41. Mallette, F. A., Gaumont-Leclerc, M. F. & Ferbeyre, G. The DNA damage signaling pathway is a critical mediator of oncogene-induced senescence. Genes Dev. 21, 4348 (2007).

42. Campisi, J. Suppressing cancer: the importance of being senescent. Science 309, 886887 (2005).

43. Gorgoulis, V. G. et al. Activation of the DNA damage checkpoint and genomic instability in human precancerous lesions. Nature 434, 907913 (2005).

44. Boveri, T. The Origin of Malignant Tumors, Williams and Wilkins, Baltimore, 1929.

45. Pellman, D. Cell biology: aneuploidy and cancer. Nature 446, 3839 (2007).46. Merlo, L. M., Pepper, J. W., Reid, B. J. & Maley, C. C. Cancer as an evolutionary and ecological process. Nat. Rev. Cancer 6, 924935 (2006).

47. Kamijo, T. et al. Tumor suppression at the mouse INK4a locus mediated by the alternative reading frame product p19ARF. Cell 91, 649659 (1997).

48. Soussi, T. & Beroud, C. Assessing TP53 status in human tumours to evaluate clinical outcome. Nat. Rev. Cancer 1, 233240 (2001).

49. Meek, D. W. Tumour suppression by p53: a role for the DNA damage response? Nat. Rev. Cancer 9, 714723 (2009).

50. Wallach, D. et al. Tumor necrosis factor receptor and Fas signaling mechanisms. Annu. Rev. Immunol. 17, 331367 (1999).

51. Heinlein, C., Deppert, W., Braithwaite, A. W. & Speidel, D. A rapid and optimization-free procedure allows the in vivo detection of subtle cell cycle and ploidy alterations in tissues by ow cytometry. Cell Cycle 9, 35843590 (2010).

52. Guo, Z. et al. Human DNA polymerase beta polymorphism, Arg137Gln, impairs its polymerase activity and interaction with PCNA and the cellular base excision repair capacity. Nucleic Acids Res. 37, 34313441 (2009).

53. Zhang, Y., Xiong, Y. & Yarbrough, W. G. ARF promotes MDM2 degradation and stabilizes p53: ARF-INK4a locus deletion impairs both the Rb and p53 tumor suppression pathways. Cell 92, 725734 (1998).

54. Zheng, L. et al. Human DNA2 is a mitochondrial nuclease/helicase for efficient processing of DNA replication and repair intermediates. Mol. Cell. 32, 325336 (2008).

55. Rauch, T., Li, H., Wu, X. & Pfeifer, G. P. MIRA-assisted microarray analysis, a new technology for the determination of DNA methylation patterns, identies frequent methylation of homeodomain-containing genes in lung cancer cells. Cancer Res. 66, 79397947 (2006).

56. Juttermann, R., Li, E. & Jaenisch, R. Toxicity of 5-aza-2-deoxycytidineto mammalian cells is mediated primarily by covalent trapping of DNA methyltransferase rather than DNA demethylation. Proc. Natl Acad. Sci. USA 91, 1179711801 (1994).

Acknowledgements

We thank V. Bedell, D. Bobadilla, J. Collins and M. Slovak in the City of Hope Cytogenetics Core Facility for technical assistance with the karyotype analysis and FISH detection of polyploid cells. We thank M. Lee and B. Armstrong in the City of Hope Microscopy Core Facility for assistance with immunouorescence, TUNEL and X-gal staining. We thank L. Brown in the City of Hope Flow Cytometry Core for assistance with FACS sorting and analysis. We thank Drs. K. Walker, M. Morgan and S. Tseo for editorial assistance. This work was supported by NIH Grants R01 CA085344 and R01 CA073764 to B. H. S, by P30 CA33572 from the National Cancer Institute to City of Hope, and by the Lung Cancer Program of the City of Hope Comprehensive Cancer

NATURE COMMUNICATIONS | 3:815 | DOI: 10.1038/ncomms1825 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1825

Center. This work is also based upon work partially supported by the U. S. Army RDECOM Acquisition Center, Natick Contracting Division, Natick, MA, under Contract No: W911QY-10-C-0176.

Author contributions

L.Z. designed and conducted biochemical and cellular assays, and contributed to manuscript preparation. H.F. characterised diploid normal, polyploid cancer cells and performed mouse phenotype analyses. M.Z. conducted DNA methylation analyses. X.L., X.W., J.Z. and C.W. conducted gene expression proling analysis. C.L. and Z.G. conducted biochemical assays with puried recombinant proteins. J. W. performed tumour graing in NOD/SCID mice. Q.H. and G.J. coordinated and performed histological and cytogenetic analyses. P. G. designed and supervised DNA methylation proling analyses. B.S. supervised the entire project, designed and coordinated most of experiments in this study, and contributed to manuscript preparation.

Additional information

Accession codes: The microarray data have been deposited in the GEO database under accession code GSE18815.

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

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

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

How to cite this article: Zheng L. et al. Polyploid cells rewire DNA damage response networks to overcome replication stress-induced barriers for tumour progression. Nat. Commun. 3:815 doi: 10.1038/ncomms1825 (2012).

NATURE COMMUNICATIONS | 3:815 | DOI: 10.1038/ncomms1825 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.