Learn more: PMC Disclaimer | PMC Copyright Notice
Topoisomerase I suppresses genomic instability by preventing interference between replication and transcription
Associated Data
Abstract
Topoisomerase I (Top1) is a key enzyme acting at the interface between DNA replication, transcription and mRNA maturation. Here, we show that Top1 suppresses genomic instability in mammalian cells by preventing conflicts between transcription and DNA replication. Using DNA combing and ChIP-on-chip, we found that Top1-deficient cells accumulate stalled replication forks and chromosome breaks in S phase and that breaks occur preferentially at gene-rich regions of the genome. Strikingly, these phenotypes were suppressed by preventing the formation of RNA-DNA hybrids (R-loops) during transcription. Moreover, these defects could be mimicked by depletion of the splicing factor ASF/SF2, which interacts functionally with Top1. Taken together, these data indicate that Top1 prevents replication fork collapse by suppressing the formation of R-loops in an ASF/SF2-dependent manner. We propose that interference between replication and transcription represents a major source of spontaneous replication stress, which could drive genomic instability during early stages of tumorigenesis.
Introduction
DNA replication is a complex process that involves the coordinated activation of thousands of replication origins distributed along the chromosomes 1. Replication forks progressing from these origins frequently stall when they encounter obstacles such as DNA lesions or tightly-bound protein complexes 2, 3. Arrested forks are unstable structures, which must be readily processed to avoid inappropriate recombination and genomic instability 3, 4. It has been recently proposed that activated oncogenes increase the rate of replication fork stalling and chromosome rearrangement at common fragile sites (CFS) in precancerous lesions 5, 6. This chronic replication stress promotes the bypass of anticancer barriers such as checkpoints and senescence 7. However, the origin of this replication stress and the nature of sites that impede fork progression in the human genome are currently unknown.
Approximately 1400 natural replication pause sites have been identified in the yeast genome, which include centromeres, telomeres, inactive replication origins and highly-expressed genes 2, 8. Replication through active genes has been implicated in genomic instability through a process called transcription-associated recombination (TAR) 4. TAR does not result from frontal collisions between RNA and DNA polymerases but is rather due to RNA-DNA hybrids (R-loops) that form when the assembly of mRNA-particle complexes (mRNPs) is perturbed 4. In agreement with this model, several mutations affecting the maturation of mRNPs increase TAR in budding yeast 4. TAR-related events have also been reported in human cells 9, 10. It is therefore likely that transcription represents a potential source of replication stress and genomic instability in human cells, as it is the case in budding yeast.
DNA topoisomerase I (Top1) is a ubiquitous enzyme that plays multiple biological functions at the crossroads between replication, transcription and mRNA maturation 11. Top1 relaxes DNA supercoiling generated by transcription, replication and chromatin remodelling 12. It is essential for viability in metazoan 13, 14 and Top1 depletion by RNA interference induces DNA damage in human cells 15. Besides its DNA relaxation activity, Top1 is also implicated in the regulation of mRNA splicing in higher eukaryotes, presumably through the direct phosphorylation of splicing factors of the Serine/Arginine (SR)-rich family 16–18. Since one of these factors, ASF/SF2, prevents DNA breaks by avoiding the formation of R-loops 19, we asked whether Top1 suppresses genomic instability in higher eukaryotes by coordinating replication and transcription.
Results
Top1-deficient cells accumulate chromosome breaks in S phase
We first monitored genomic instability in murine B lymphoma-derived cells (P388), Top1-deficient subclone (45/R) and 45/R cells complemented with human Top1-GFP (21/P; Fig. 1a). Using comet assay, a 4.4-fold increase of DNA breaks was detected in 45/R cells relative to control cells (Fig. 1b, c, Table S1). We also observed a sharp increase of histone H2AX phosphorylation (γ-H2AX) 20 in BrdU-positive Top1-deficient cells, indicating that chromosome breaks occur in S phase (Fig. 1d, e, Table S1). Importantly, both DNA breaks and γ-H2AX foci were suppressed upon complementation of 45/R cells with Top1-GFP (Fig. 1c–e). Moreover, M-FISH analysis revealed a 6-fold increase of chromosome rearrangements in Top1- cells (Fig. 1f, g). This increased genomic instability is reminiscent of the phenotype of human HCT116 cells expressing Top1 shRNA (shTop1), which accumulate chromosome breaks (Table S1) and display activation of S-phase checkpoints 15. Interestingly, Top1 depletion also induced a 2- to 3-fold increase of the frequency of chromosome breaks at the common fragile sites (CFSs) FRA3B, FRA16D and FRAXD (Fig. 1h, i), which are frequently rearranged in cancer cells 21, 22. These data indicate that Top1 plays an important role in S phase by preventing DNA breaks and chromosomal aberrations.
Replication forks are slower in the absence of Top1
Since CFSs break more frequently upon replication stress, we asked whether the chronic genomic instability in Top1- cells results from replication defects. Murine cells were pulse-labelled with BrdU to identify newly-replicated regions and the length of BrdU tracks was monitored along individual DNA fibres by DNA combing 23, 24. This analysis revealed that BrdU tracks are shorter in Top1-deficient cells (17.9 kb) than in control and complemented cells (28.5 and 28.2 kb; Fig. 2a, b, S1). Forks are therefore ~50% slower in the absence of Top1. Interestingly, we also observed a concomitant increase of the initiation rate in Top1- deficient cells, as centre-to-centre distances between BrdU tracks were significantly shorter in Top1-deficient cells (54.7 kb) than in complemented and control cells (94.5 kb and 111.6 kb; Fig. 2c, S1).
In human shTop1 cells, BrdU tracks length and centre-to-centre distances were also found to be shorter than in control cells (Fig. S2). This effect was slightly lower in than in Top1-deficient murine cells, probably because of differences in residual levels of Top1. In order to determine elongation rates, human cells were pulse-labelled with IdU and CldU and the distance covered by individual forks during the pulse was determined. This analysis revealed that forks move at 0.7 kb/min in shTop1 cells versus 1.1 kb/min in control HCT116 cells (Fig. S2a). Similar results were obtained upon transient depletion of Top1 with siRNA (Fig. 2d, 2e, S3). We therefore conclude that Top1 is required for normal fork progression in mammalian cells and that backup origins fire in Top1-deficient cells to compensate for slower forks.
Shorter BrdU tracks observed in Top1-deficient cells could be due to slower forks or to increased fork stalling. To discriminate between these two possibilities, progression of sister replication forks was analysed by DNA combing (Fig. 3). In control cells, sister forks progress at a similar rate from a given origin 25 and generate symmetrical patterns of IdU/CldU incorporation (Fig. 3a). In contrast, more than 50% of asymmetrical patterns were detected in murine Top1- cells (Fig. 3b). Analysis of the ratio of the longest to the shortest IdU signals for each pair of sister replication forks also revealed a 3-fold increase in fork asymmetry in Top1-deficient cells (42%) relative to complemented and control cells (17.2 and 16.4 %; Fig. 3c, S4a). Similarly, a significant increase of sister fork asymmetry was detected in Top1-depleted human cells (Fig. 3d–f, S4b–c). Together, these results indicate that replication forks are not only slower, but also pause or stall more frequently in the absence of Top1.
Top1 modulates the function of ASF/SF2 to prevent genomic instability in S phase
Next, we asked whether the replication stress observed in Top1-deficient cells results from the accumulation of supercoiled DNA or whether it reflects defects in the regulation of RNA splicing. Indeed, Top1 has been implicated in the regulation of mRNP assembly, presumably through the binding and the phosphorylation of splicing factors of the SR family 16–18. Since the THO/TREX complex prevents fork collapse in yeast by promoting mRNP assembly and preventing the formation of R-loops 4, we reasoned that Top1 could play a similar role in metazoan by regulating RNA splicing. In order to address this possibility, we have complemented Top1-deficient human cells with the S. cerevisiae Top1 enzyme (ScTop1), which is required for normal fork progression but does not have a kinase domain and does not regulate fork pausing (Fig. S6). The relaxation activity of Top1 is highly conserved between yeast and human cells 26. Moreover, it has been reported that ScTop1 is targeted to the nucleus in mammalian cells and is fully able to relax supercoiled DNA 27. Here, we found that although transfection with ScTop1 significantly reduced DNA breaks in human Top1- cells (Fig. 4a, S5), this suppression was much less efficient than complementation with human Top1-GFP (Fig. 1c). These data indicate that accumulation of supercoiled DNA is not the only cause of DNA breaks in Top1- cells.
To test whether RNA splicing is altered in Top1- cells, we next monitored the subnuclear localization of SR proteins by indirect immunofluorescence. Indeed, it has been reported that SR proteins form speckles in the nucleus and are recruited to transcription sites by phosphorylation 28. Interestingly, we found that the distribution of phospho-SR proteins was severely perturbed in Top1-deficient cells, the absence of speckles correlating with the formation of γ-H2AX foci (Fig. 4b). We therefore depleted the SR protein ASF/SF2 using RNA interference (Fig. S7a) and monitored the effect of this depletion on genomic instability. A sharp increase of both DNA damage (Fig. 4c) and chromosome breaks (Fig. 4d, e) was detected in ASF/SF2-depleted cells, which is reminiscent of Top1-deficient cells. Importantly, depletion of both Top1 and ASF/SF2 did not further increase DNA breaks, which indicates that both proteins work in the same pathway to prevent genomic instability in S phase. This view is also supported by the fact that Diospyrin D1, an inhibitor of Top1 kinase 29, induced a 10-fold increase of DNA breaks in Top1-proficient cells (Fig. 4f, S7b, Table S2), in a replication-dependent manner (Fig. 4g, Table S2). Finally, we found that both Diospyrin and ASF/SF2 depletion induced a sharp increase of sister fork asymmetry in human and murine cells (Fig. 4h, i, S7c, d), to a level comparable to Top1 depletion (Fig. 3). We therefore propose that mammalian Top1 suppresses replication stress in S phase both by relaxing DNA supercoiling and by regulating mRNA splicing.
RNaseH1 suppresses fork asymmetry and chromosome breaks in Top1- cells
The above results suggest that transcription interferes with replication fork progression in Top1-deficient cells through the accumulation of R-loops. To test this possibility, Top1- cells were treated with Cordycepin, a potent inhibitor of RNA chain elongation, and sister forks progression was monitored by DNA combing. Cordycepin suppressed sister-fork asymmetry in murine and human Top1- cells (Fig. 5a, b, Table S3), supporting the view that transcription interferes with fork progression in Top1- cells. To confirm the implication of R-loops in this process, human cells were transfected with a vector expressing RNaseH1 to degrade DNA-RNA hybrids. Remarkably, RNaseH1 induced a reduction of both sister fork asymmetry and DNA breaks (Fig. 5c–e, Table S3) in shTop1 cells. In contrast, control cells were not affected by this treatment.
Since Cordycepin and RNaseH1 suppress sister fork asymmetry in Top1-deficient cells, we next asked whether these treatments restore normal fork progression in the absence of Top1. Conversely, we checked whether increased fork stalling induced by Diospyrin or ASF/SF2 depletion also reduced fork rates in control cells. To this end, fork speed was measured in cells exposed to Cordycepin, RNaseH1, Diospyrin and siASF (Fig. 5f, S8, Table 1). We found that Diospyrin and siASF induced a 0.33 kb/min reduction of fork rate in control cells (1.2 kb/min), which represents ~70% of the delay measured in Top1-deficient cells (−0.45 kb/min relative to control cells). Conversely, 46% of the delay measured in Top1-deficient cells was suppressed upon treatment with Cordycepin or RNaseH1 (+0.21 kb/min; Fig. 5f, S8, Table 1). Together, these data indicate that the slow fork phenotype of Top1-deficient cells is largely caused by transcription-dependent replication fork pausing.
Table 1
Fork rate (kb/min) | P388 | 45/R | HCT116 | shTop1 | |||||||
---|---|---|---|---|---|---|---|---|---|---|---|
Ctrl | D1 | siASF | Ctrl | Cord. | Ctrl | D1 | siASF | Ctrl | Cord. | RNaseH | |
Nb. values | 137 | 146 | 103 | 141 | 112 | 122 | 133 | 108 | 149 | 118 | 103 |
Median | 1.22 | 0.96 | 0.78 | 0.73 | 0.90 | 1.18 | 0.85 | 0.91 | 0.77 | 1.01 | 0.95 |
Mean | 1.25 | 0.94 | 0.88 | 0.79 | 0.90 | 1.16 | 0.89 | 0.92 | 0.83 | 1.02 | 1.01 |
St. Dev. | 0.41 | 0.34 | 0.36 | 0.29 | 0.29 | 0.31 | 0.28 | 0.27 | 0.34 | 0.32 | 0.39 |
P-value | - | *** | *** | - | ** | - | *** | *** | - | *** | *** |
- | 0.06 | - | - | - | 0.44 | - | 0.34 |
γ-H2AX is preferentially detected at transcribed regions in Top1- cells
To test whether DNA breaks observed in Top1-deficient cells occur preferentially at transcribed loci, we monitored the distribution of γ-H2AX by ChIP-on-chip on high-resolution tiling arrays encompassing human chromosomes 1 and 6 (Fig. 6a). This analysis identified 435 loci in shTop1 and 149 loci in shCtrl cells (Fig. 6b) averaging ~1 kb in length (Fig. 6c). Interestingly, 50% of γ-H2AX loci were found within 2 kb of the 5′- or 3′-boundary of annotated genes in shTop1 cells, a frequency that is significantly higher than expected for a random distribution (11.2%; Fig. 6d). In contrast, no significant enrichment was detected in shCtrl cells (8.3%; Fig. 6d). This is best illustrated for SFRS3 (Fig. 6e), one of the most highly expressed genes on chromosome 6. Significant γ-H2AX enrichment near promoters, termination sites and at converging genes was also detected by calculating average γ-H2AX profiles for the complete set of protein-coding genes on human chromosomes 1 and 6 (Fig. 6f).
We next asked whether γ-H2AX enrichment at annotated genes correlates with transcriptional activity in S phase. To this end, we focused on replication-dependent histone genes, whose expression is linked to DNA replication by transcriptional and posttranscriptional mechanisms 30. We found that about 70% of the replication-dependent histone genes were enriched in γ-H2AX specifically in shTop1 cells (Fig. 6g, S9). Moreover, the five major histone H4 genes (representing more than 80% of histone H4 mRNAs 30) were highly enriched in γ-H2AX in shTop1 cells (Fig. S8). We found a positive correlation between the level of expression of H4 genes and their probability to overlap with γ-H2AX domains (r2=0.71; Fig. 6h). Collectively, these data indicate that, in the absence of Top1, DNA breaks form at highly-expressed genes during S phase and that the rate of DNA damage correlates with the level of gene expression.
Discussion
It is generally believed that Top1 relaxes DNA supercoiling that accumulates ahead of replication forks. However, this view has been recently challenged by reports showing that Top2 relaxes chromatin templates more efficiently than Top1 31 and that Top1 is dispensable for normal DNA replication in budding yeast 32. Here, we have used DNA combing to monitor fork progression in top1Δ yeast mutants and in mammalian cells presenting low to undetectable levels of Top1. We found that forks are ~50% slower in these cells, confirming thereby that Top1 plays a major role in replication elongation. We also found that Top1- cells compensate for slow forks by increasing the rate of initiation. These data are consistent with other studies showing that dormant replication origins fire in response to various types of fork impediments 33–35. They also explain why yeast top1Δ mutants have a normal S phase 32 although their forks are 50% slower.
We also show that, besides slow replication forks, Top1-deficient mammalian cells display an increased rate of fork stalling. Our data indicate that fork arrest in Top1-deficient cells is largely due to the accumulation of R-loops during transcription. Since DNA supercoiling accumulates at transcription sites in the absence of Top1, an attractive possibility could be that nascent RNA chains reanneal with the DNA template in Top1-deficient cells and interfere with fork progression. This view is supported by studies in E. coli showing that Top1 prevents the formation of abnormal RNA-DNA hybrids by relieving negative supercoiling behind RNA polymerases 36. This model is also consistent with the fact that DNA breaks in Top1-human cells are at least partially suppressed by S. cerevisiae Top1. Alternatively, Top1 could prevent R-loop formation by promoting the ASF/SF2-dependent assembly of mRNPs 16, 17, 37. In agreement with this model, we found that the subnuclear organization of ASF/SF2 speckles is profoundly altered in Top1-deficient cells. Moreover, inhibition of Top1 kinase with Diospyrin or depletion of its target ASF/SF2 induces fork arrest and chromosome breaks to a similar extent as in Top1-cells. Importantly, no additive effect of ASF/SF2- and Top1- depletion was detected, which indicates that both proteins act in the same pathway to prevent replication stress. We therefore conclude that Top1 avoids conflicts between replication and transcription both by relaxing DNA supercoiling and by promoting mRNP assembly (Fig. 7).
A simple prediction of the model described above is that chromosome breaks occur preferentially at gene-rich regions in the absence of Top1. To test this possibility, we mapped γ-H2AX enriched loci by ChIP-on-chip. This analysis identified a significant overlap between γ-H2AX-enriched regions and annotated genes, with a preferential localization at 5′- and 3′-boundaries. We also observed a strong accumulation of γ-H2AX at replication-dependent histone genes, which is proportional to their level of expression. γ-H2AX-enriched regions encompass ~1 kb on average which is strikingly smaller than the megabase-length domains that form at DNA double-strand breaks induced by ionizing radiation 20. They are reminiscent of the γ-H2AX tracks that form at natural replication pause sites in budding yeast (our unpublished results and D. Durocher, personal communication) and probably correspond to checkpoint signalling of stalled forks. These observations support the view that transient breaks accumulate in Top1-deficient cells as the consequence of conflicts between replication and transcription.
It is now well established that TAR represents a prominent source of genomic instability in budding yeast 4. Together with earlier studies 9, 10, 19 our results indicate that TAR-related events also occur in mammals and are suppressed in a Top1- and ASF/SF2-dependent manner. Intron-containing genes represent only a small fraction of the human genome, which may question the significance of TAR as a general source of genomic instability. However, other types of transcriptional activity could promote TAR in human cells. Indeed, stalled RNA polymerases have been detected at the promoter of most inactive genes 38 and small nuclear RNAs complementary to the 5′- and 3′-ends of genes have also been identified 39. Moreover, unbiased transcriptome analyses have recently shown that most of the human genome is transcribed 40. Since ASF/SF2 is permanently associated with RNA pol II and is required for optimal elongation 37, 41, it is tempting to speculate that pervasive transcription interferes with fork progression and is regulated in a similar Top1-dependent manner. Since it is believed that R-loops affect DNA replication by generating marks that persist after transcription 4, gene expression could therefore interfere with replication even if these processes do not occur simultaneously 42, 43.
Recent evidence indicates that the DNA damage checkpoint is constitutively induced in precancerous lesions due to persistent replication defects 5, 6 and acts as a selective pressure to inactivate p53 7. However, the origin of replication stress in precancerous cells is currently unknown. Here, we show that transcription represents a source of replication fork stalling in mammalian cells. Since transcription units act as polar replication barriers, alterations of origin usage in oncogene-activated cells 44, 45 could modify the direction of fork progression and increase conflicts with transcription. It was recently proposed that transcription at highly-expressed genes is co-oriented with the replication fork in a large fraction of the human genome, presumably to minimise these conflicts 46. This view is supported by our ChIP-on-chip data, which show that γ-H2AX is preferentially detected outside of these organised domains in Top1- cells (data not shown).
The human genome contains a large excess of replication origins that can be used as backup initiation sites when fork progression is impeded 47, 48. Likewise, we found here that Top1-depleted cells activate new origins to rescue stalled forks (Fig. 7). Common fragile sites (CFSs) are specific regions of the genome that are prone to breakage upon replication stress and induce chromosomal rearrangements from the early stages of tumorigenesis 7, 21, 22. These loci end their replication very late in S phase and are frequently associated with very large genes 49, 50. We observed a significant increase of chromosome breaks at CFSs in Top1- cells. It is therefore tempting to speculate that regions of the genome that are either devoid of dormant origins or that replicate very late in S phase are hypersensitive to fork stalling induced by transcription, this sensitivity further increasing in cancer cells due to perturbations of the replication program.
Methods
Cell lines and culture conditions
Cells lines and culture conditions for P388 murine leukaemia cell line (control), 45/R cells (CPT-resistant subclone which do not display Top1 activity) and 21/P cells (45/R cells complemented with human Top1-GFP) were described previously 17. HCT116 and derived HCT116-shTop1 cells were grown in DMEM, with 500 μg/ml Hygromycin B as described 15.
RNAi, transient transfections and drugs
The following siRNA duplexes designed to repress ASF/SF2: siASF1 (5′-GUAUUGACCUUAUACUAAAdTdT-3′) and siASF2 (5′-GGGUAGCAAUGCCAGUAAAdTdT-3′) were transfected using INTERFERin (Polyplus-transfection) at 1 nM and 10 nM, respectively. Depletion was verified by western blot with antibodies against SF2/ASF (Zymed Laboratories) 24 to 48 hrs after transfection. For Top1 transient downregulation, cells were transfected with the following siRNA (5′-GGACTCCATCAGATACTAT-3′) and were analysed 72 hrs after transfection as described 15. For RNaseH1 expression studies, 5×104 HCT116 cells were seeded 24 hrs before transfection and were transfected with 1 μg of pCMV6-XL5-RNaseH1 vector (Origene, Rockville, MD) and 1 μg of pHygEGFP reporter vector (Clontech) using jetPEI reagents (Polyplus-transfection). Cells were analysed 24 hrs after transfection. Cells were analysed 48 hrs after transfection. Transcription was inhibited by exposing murine cells for 6 hrs to 100 μM Cordycepin and human cells for 3 hrs to 50 μM Cordycepin. Inhibition of Top1 kinase activity was performed by treating cells with the following concentrations Diospyrin (IC50): shCtrl HCT116 (35 μM), shTop1 HCT116 (10 μM), P388 (10 μM), 45R (5 μM), 21P (7.5 μM).
Comet Assay
DNA breaks were monitored using the CometAssay Reagent Kit for Single Cell Gel Electrophoresis Assay (Trevigen, Inc., Gaithersburg, MD) according to the manufacturer’s instructions with minor modifications. Briefly, 75 μl of cells resuspended in LMP agarose (3×105 cells/ml) were pipetted onto Trevigen CometSlides. After electrophoresis with a horizontal apparatus (30 min at 1 V/cm in 1×TBE), slides were stained with 50 μl SYBR green dye (Trevigen; 1/10,000 in TE) and viewed using a Leica epifluorescence microscope. Imaging was performed using the Comet Imager Software V 2.0.105 (Metasystem, Germany). Tail moment (TM) considers both the tail length (TL) and the fraction of DNA in the comet tail (TM = %DNA in tail x TL/100). A total of 100 cells were analysed per slide.
Analysis of metaphase spreads, FISH and M-FISH
Cells were incubated for 2 hrs with 10 μM nocodazole, harvested then spread on slides and analysed as described 21. FISH with BACs RP11-147N7 (FRA3B), RP11-22N7 (FRA16D) RP11-323J20 (FRAXB) was performed as described 21 For each cell line, the frequency of chromosome breaks were scored on more than 100 metaphases. Multiplex-Fluorescence In Situ Hybridization (M-FISH) analysis of metaphase spreads was performed as described 51. For each experiment, 20 metaphase spreads were acquired by using a Leica DM RXA epifluorescence microscope equipped with a Sensys CCD camera (Photometrics) and the Q-FISH software (Leica). Images were processed using MCK software (Leica Microsystems Imaging Solutions).
Indirect immunofluorescence
Cells were labelled for 10 min with 10 μM BrdU (Sigma). Murine leukaemia cells were resuspended in PBS (106 cells/ml) and cytospined (Shandon Cytospin 3, Thermo Scientific). Cells were fixed for 15 min with 4% formaldehyde at room temperature, washed twice with PBS and incubated with 4N HCl for 30 min. Cells were permeabilised with 0.5% Triton X100 in PBS for 10 min and then incubated in PBS, 3% BSA (Sigma-Aldrich) for 30 min at 37°C. BrdU was detected with a rat monoclonal antibody (BU1/75, AbCys) and a secondary antibody-coupled Alexa 488 (Molecular Probes). γ-H2AX (Ser139) was detected with a rabbit antibody (2577, Cell Signalling) and an anti-rabbit IgG conjugated with TRITC (Jackson ImmunoResearch). Phospho-SR proteins were detected with the monoclonal antibody Mab104 17 and an anti-mouse IgG conjugated with TRITC (Jackson ImmunoResearch). Slides were mounted with Vectashield containing 1.5 μg/ml DAPI.
DNA combing
DNA combing was performed as described 23, 52. Briefly, cells DNA fibres were extracted in agarose plugs immediately after BrdU labelling and were stretched on silanized coverslips. BrdU was detected with a rat monoclonal antibody (BU1/75, AbCys; 1/20) and a secondary antibody coupled to Alexa 488 (A11006, Molecular Probes; 1/50). DNA molecules were counterstained with an anti-ssDNA antibody (MAB3034, Chemicon; 1/500) and an anti-mouse IgG coupled to Alexa 546 (A11030, Molecular Probes, 1/50). CldU and IdU were detected with BU1/75 (AbCys, 1/20) and BD44 (Becton Dickinson, 1/20) anti-BrdU antibodies, respectively. DNA fibres were analysed on a Leica DM6000B microscope equipped with a CoolSNAP HQ CCD camera (Roper Scientifics). Data acquisition was performed with MetaMorph (Roper Scientifics). Representative images of DNA fibres were assembled from different fields of view and were processed as described 53.
Graphs and statistical analysis
Box-and-whiskers graphs were plotted with Prism v5.0 (GraphPad Software). For all graphs, top and bottom of the box correspond to the 25th and 75th percentile (the lower and upper quartiles, respectively) and the line near the middle of the box marks the median (50th percentile). Whiskers correspond to the 10–90 percentiles. Data not included between the whiskers are plotted as outliers (dots). Statistical analysis was performed in Prism v5.0 (GraphPad Software) using the non-parametric Mann-Whitney rank sum test. Normality of distributions was assessed with the Kolmogorov-Smirnov test. Linear regressions and r2 were computed with Prism v5.0 (GraphPad Software).
Cloning and expression of S. cerevisiae Topoisomerase I
cDNA of scTop1 cDNA was amplified by PCR with primers 5′-GGGAGCTCATGACTATTGCTGATGCTTCC-3′ (introducing a SacI site) and 5′-GGCCCGGGTTAAAACCTCCAATTTTCATCTACC-3′ (introducing a SmaI site). PCR products were cloned into pIRES2-AcGFP1 (Clontech Laboratories) using SacI and SmaI cloning sites. Sequence of ScTop1 was verified with primers 5′-TTTCTCGTTGCGATTTCAC-3′, 5′-AAACTACATGCCGGGATT-3′ and (‘-GCAAATGGGCGGTAGGCGTG-3′ (CMV promoter). 3×105 HCT116 cells were seeded 24 hrs before transfection with 6 μg of pIRES2-AcScTop1-GFP1 vector and 6 μg of pIRES2-AcGFP1 reporter vector using jetPEI reagents (Polyplus-transfection). Expression of ScTop1 was checked by RT-PCR with a Light Cycler LC480 (Roche) and cells were analysed 48 hrs after transfection.
ChIP-on-chip
Asynchronous cultures of HCT116 cells were harvested with Trypsin, washed with PBS, and fixed with 1% formaldehyde in PBS during 10 minutes at room temperature. The ChIP experiments were performed as previously described 54 with the following modifications. Briefly, 5 μg of anti γ-H2AX antibody (Abcam ab2893) was first bound to 20 μL of Dynabeads Protein A (Dynal, Invitrogen) then added to 200 μL of lysate, normalised to 2.106 cells. After purification, the immunoprecipitated DNA was amplified using the WGA2 GenomePlex Complete Whole Genome Amplification kit (SIGMA) with addition of 0.1 mM dUTP. Seven μg of amplified DNA was fragmented and labelled using the GeneChip WT Double-Stranded DNA Terminal Labelling kit (Affymetrix, PN 900812) following manufacturer’s recommendations. DNA was hybridised to Affymetrix GeneChip Human Tiling 2.0R A Array Set (Chromosome 1, 6), using the GeneChip Hybridization, Wash, and Stain Kit (Affymetrix, PN 900720). Tiling arrays were scanned with the GeneChip scanner 3000 7G. Enrichment was calculated with the Tiling Analysis Software v1.1 (Affymetrix) and profiles were generated with the Affymetrix Integrated Genome Browser (v5.01). Data from two biological replicates were combined and loci enriched in γ-H2AX with a P-value lower than 0.01 were scored. Enrichment analyses were performed with the Affymetrix Tiling Analysis Software (v. 1.1.02).
Supplementary Material
Supplementary Figure S1
Single-molecule analysis of DNA replication in Top1-deficient murine leukemia cells. P388 (control), 45/R (Top1-) and 21/P (Top1-GFP) cells were pulse-labelled for 15 min with BrdU and DNA fibres were stretched on silanized coverslips. (a) Distribution of BrdU tracks length. (b) Distribution of centre-to-centre distances between BrdU tracks. (c) Statistics for BrdU tracks length and centre-to-centre distances in P388, 45/R and 21/P cells. See Fig. 2a-c for experimental details.
Supplementary Figure S2
Single-molecule analysis of DNA replication in shCtrl and shTop1 HCT116 cells. (a) Box plots of fork velocity, BrdU tracks length and centre-to-centre distances between BrdU tracks in shCtrl and shTop1 HCT116 cells. (b) Frequency distribution of data displayed in above panels. (c) Statistics of fork velocity, BrdU tracks length and centre-to-centre distances. Fork velocity was determined as described in Fig. 2d. Analysis of BrdU tracks length and centreto- centre distances between BrdU tracks was performed as described for P388 cells (Fig. 2a-c).
Supplementary Figure S3
Single-molecule analysis of DNA replication in human HCT116 cells transfected with siRNA against Top1. (a) Frequency distribution of BrdU tracks length and centre-to-centre distances in siCtrl and siTop1 HCT116 cells. (b) Statistical analysis of DNA combing data. See Fig. 2d,e for experimental details. (c) Western blot analysis of Top1 levels in HCT116 cells transfected with control or Top1 siRNAs.
Supplementary Figure S4
DNA combing analysis of replication fork stalling in Top1-depleted cells. Cells were pulselabelled for 15 min with IdU and 15 min with CldU and processed as described (Fig. 3) (a) Statistical analysis of fork asymmetry in P388, 45/R and 21/P cells. See Fig. 3b,c for experimental details. (b) Statistical analysis of fork asymmetry in shCtrl and shTop1 HCT116 cells. (c) Frequency of fork arrest in cells transiently transfected with siTop1 and siCtrl, calculated as described previously (Conti et al., 2007).
Supplementary Figure S5
Complementation of Top1-deficient human cells with S. cerevisiae TOP1. (a) Frequency distribution of comet tail moments in shCtrl and shTop1 HCT116 cells, transfected (ScTop1) or not (empty vector) with a vector expressing the S. cerevisiae TOP1 gene. Note that ScTop1 overexpression is not toxic in HCT116 cells as it did not significantly increase the incidence of DNA breaks. (b) Statistical analysis of comet tail moment in shCtrl and shTop1 HCT116 cells complemented or not with yeast TOP1. (c) RT-PCR analysis of ScTop1 mRNA levels in shCtrl and shTop1 cells transfected with an empty vector (-) or with a vector expressing the S. cerevisiae TOP1 gene (+ ScTop1). ScTop1 mRNA levels were normalized to human HPRT and 18S RNAs.
Supplementary Figure S6
Single-molecule analysis of DNA replication in S. cerevisiae wild type and top1Δ cells. (a) Frequency distribution of BrdU tracks length and inter-origin distance in isogenic wild type (wt) and top1Δ S. cerevisiae cells. Asynchronous cultures were pulse-labelled for 20 minutes with BrdU and newly-replicated tracks were analysed by DNA combing. BrdU tracks are shorter in top1Δ mutants (19.7 kb) than in wild type cells (37.2 kb) but the density of BrdU tracks is higher in top1Δ cells (b) Asynchronous cultures of wild type and top1Δ cells were pulse-labelled for 10 minutes with IdU and for 15 minutes with CldU. Fork asymmetry was determined in wild type and top1Δ cells as described in fig 3. Medians are indicated. No statistical difference is observed between wt and top1Δ cells. (c) Representative DNA fibres from wild type and top1Δ cells. Red: IdU, Green: CldU, Blue in merged image: DNA. ORI: position of replication origin. Bar: 20 kb. (d) Statistical analysis of BrdU tracks length, interorigin distance and fork asymmetry in wild type and top1Δ cells.
Supplementary Figure S7
ASF/SF2 depletion increases replication fork pausing and DNA breaks in Top1-proficient cells, but not in Top1-deficient cells. P388, HCT116 and shTop1 HCT116 cells were treated with Diospyrin or with siRNAs against ASF/SF2 (siASF1 or siASF2). (a) Western blot analysis of ASF/SF2 depletion in HCT116 shCtrl cells. (b) Analysis of DNA damage (comet assay) in human HCT116 cells treated (+) or not (-) with diospyrin (D1). (c, d) Fork asymmetry was monitored by DNA combing. (e, f) DNA breaks were quantitated by comet assay, as described in Fig. 4.
Supplementary Figure S8
Frequency distribution of replication fork rates in P388, 45/R, HCT116 shCtrl and shTop1 cells exposed or not to Dyospyrin (D1), siASF1, Cordycepin (Cord.) and RNaseH. Median fork rate is indicated in red. See Fig. 5 for experimental details.
Supplementary Figure S9
ChIP-on-chip analysis of ɣ-H2AX enrichment at the three major histone gene clusters on human chromosome 1 and 6. Chromatin immunoprecipitation and microarray analysis were performed as described in Fig. 6 for two biological replicates. The position of histone genes on (+) and (-) strands is shown in red. Non-histone genes are indicated in black. ɣ-H2AX-enriched loci (p<0.01) are indicated. Red: shTop1. Gray: shCtrl. The five nonallelic histone H4 genes contributing to more than 80% of histone H4 mRNAs are indicated with an asterisk.
Supplementary Table S1
Accumulation of chromosome breaks in Top1-deficient cells. Numerical values for the experiments shown in Fig. 1c and 1e are indicated, together with statistical differences between Top1-deficient and control cells. See Fig. 1 for experimental details. Comet tail moment in shTop1 and shCtrl HCT116 cells was measured as described for P388 cells (see Methods section).
Supplementary Table S2
Analysis of chromosome breaks in HCT116 and P388 cells exposed to Diospyrin or to an siRNA against ASF/SF2. See Fig. 4 for experimental details.
Supplementary Table S3
Analysis of fork asymmetry and DNA breaks in Top1-deficient cells exposed to Cordycepin and to RNaseH. See Fig. 5 for experimental details.
Acknowledgments
We thank Michelle Debatisse and Marcel Méchali for discussions and for critical comments on the manuscript. We are also grateful to D. Durocher for communicating data prior to publication. We thank Jacques Piette for his initial contribution to this project. We also thank Etienne Schwob and the DNA combing facility of Montpellier for providing silanized coverslips. Work in the PP lab is supported by FRM (Equipe FRM), ANR, INCa and the EMBO Young Investigator Programme. AC is supported by the Ligue contre le cancer (Comité de l’Hérault), ARC (Association pour la Recherche contre le Cancer) and INCa. ST, LC and HT were recipients of fellowships from ARC, EMBO and CNRS, respectively.
Footnotes
Accession Numbers
Microarray data presented in this article can be obtained from Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) with accession number . GSE17552
Author Contributions
AC and PP designed the project and wrote the paper; ST, LC, CC and HT performed most of the experiments; MFISH studies were done by HHG and AJ; VP and JDV performed microarray experiments. Bioinformatic analyses were done by AT; CT, YP and JT provided important experimental materials and advices. All authors contributed to the interpretation of the results and edited the manuscript.