Recovery of arrested replication forks by homologous recombination is error-prone.

Homologous recombination is a universal mechanism that allows repair of DNA and provides support for DNA replication. Homologous recombination is therefore a major pathway that suppresses non-homology-mediated genome instability. Here, we report that recovery of impeded replication forks by homologous recombination is error-prone. Using a fork-arrest-based assay in fission yeast, we demonstrate that a single collapsed fork can cause mutations and large-scale genomic changes, including deletions and translocations. Fork-arrest-induced gross chromosomal rearrangements are mediated by inappropriate ectopic recombination events at the site of collapsed forks. Inverted repeats near the site of fork collapse stimulate large-scale genomic changes up to 1,500 times over spontaneous events. We also show that the high accuracy of DNA replication during S-phase is impaired by impediments to fork progression, since fork-arrest-induced mutation is due to erroneous DNA synthesis during recovery of replication forks. The mutations caused are small insertions/duplications between short tandem repeats (micro-homology) indicative of replication slippage. Our data establish that collapsed forks, but not stalled forks, recovered by homologous recombination are prone to replication slippage. The inaccuracy of DNA synthesis does not rely on PCNA ubiquitination or trans-lesion-synthesis DNA polymerases, and it is not counteracted by mismatch repair. We propose that deletions/insertions, mediated by micro-homology, leading to copy number variations during replication stress may arise by progression of error-prone replication forks restarted by homologous recombination.

Introduction

Maintenance of genome stability requires the faithful and accurate replication of the genetic material. Genome instability is a hallmark for most types of cancer and it is strongly associated with predisposition to cancer in many human syndromes (for a review, see [1], [2]). Genome instability is manifest at two levels: at the nucleotide level, resulting in base-substitutions, frame-shifts or in micro-insertions/deletions (referred to herein as mutations); and at the chromosomal level, resulting in duplications, deletions, inversions and translocations (referred to herein as gross chromosomal rearrangements or GCRs).

Genome instability during cancer development and in other human genomic disorders may be consequences of failures in chromosome replication (for a review, see [3], [4]). Origin spacing has recently been shown to cause chromosomal fragility at some human fragile sites [5], [6]. Impediments to replication fork movements per se may also cause genome instability [7]–[9]. Indeed, both slowing down and blockages to fork progression can lead to chromosomal fragilities or GCRs in human cells and yeast models [10]–[14]. However, how a blocked replication fork leads to genetic instability remains poorly understood.

In eukaryotes, DNA replication is initiated at numerous origins along linear chromosomes, and impediments to fork progression appear unavoidable during each S-phase (for a review, see [9], [15]). Impediments to fork progression can be caused by DNA lesions, by non-histone proteins tightly bound to DNA, by sequence-caused secondary structures such as cruciform structures and possibly G-quadruplexes, by nucleotide pool imbalance and by conflicts with transcription machinery (for a review, see [16], [17]). In case of failures in fork progression, DNA replication can be completed either by the recovery of the arrested fork by fork-restart mechanisms, or as a result of the progression of a converging fork which can be ensured by activation of dormant origins [7], [15], [18]. Fork restart is presumably essential in unidirectional replication regions, such as the rDNA locus, in regions of low densities of origins, such as some human fragile sites, and when two converging forks are both impeded [5], [19], [20].

To ensure faithful and complete DNA replication, cells coordinate DNA synthesis restart with specific pathways, including DNA replication checkpoint and homologous recombination mechanisms [17]. The integrity of replication forks is guaranteed by the DNA replication checkpoint that maintains the replisome in a replication-competent state to keep DNA polymerases at the site of nucleotide incorporation [21]. It remains unclear how the DNA replication checkpoint modulates replisome activities to maintain its function [21], [22]. The DNA replication checkpoint also regulates nuclease activities (e.g. Exo1 or Mus81) which contribute to preserving the integrity of stalled forks [23], [24]. If replisome function is lost or the replisome dissociates at broken replication forks, the resumption of DNA synthesis appears to require the replisome to be rebuilt. In E.coli, restart of a collapsed fork involves homologous recombination and the PriA helicase that allows replisome components to be loaded de novo on joint-molecule structures [25], [26]. In eukaryotes, the restart of collapsed or broken replication forks is dependent upon homologous recombination, but the mechanism of origin-independent loading of the replisome remains to be described [20], [27]–[30]. It has been proposed that the repair of a double-strand break (DSB) by recombination (break-induced replication, BIR) in budding yeast similarly involves the assembly of a replication fork (for a review, see [30]–[32]). When BIR occurs outside S-phase, recombination-dependent replication fork assembly can synthesise hundreds of kilobases (Kb). However, this DNA synthesis is highly inaccurate due to frequent template switching of nascent-strands and frame-shift mutations [33], [34].

We previously reported a system that displays replication fork arrest at a specific locus in the fission yeast S. pombe. The system is a polar replication fork barrier (RFB) regulated by the Rtf1 protein binding to its RTS1 binding site [35]. The RTS1-RFB causes fork arrest because of a non-histone protein complex binding to the DNA. As proposed for other polar RFBs, the RTS1-RFB is thought to block fork progression by directly (contact between proteins and the replisome) or indirectly (topological constraint) affecting the replicative helicase activity and thereby preventing DNA unwinding [36], [37]. Recovery of the arrested fork occurs by a DSB-independent mechanism and involves the recruitment of recombination proteins at the RTS1-RFB site. We proposed that recombination proteins associate with unwound nascent strands that then anneal with the initial template to allow DNA synthesis to restart [11], [20]. The causative protein barrier then has to be removed either by DNA helicase or by the recombination machinery itself to allow fork-progression to resume [38]–[40]. Occasionally, the unwound nascent strand can mistakenly anneal with a homologous template in the vicinity of the collapsed fork, resulting in the restart of DNA synthesis on non-contiguous template. This incorrect template switch of nascent strands results in inversions and iso-acentric and dicentric chromosomes in ∼2–3% of cells/generation [11], [20]. Error-free template switching between sister-chromatids provides an efficient mechanism for filling-in single-stranded gaps left behind damage-induced stalled forks [41]. Inverted chromosome fusions in yeast and rare-genome rearrangements in human genomic disorders, may both be consequences of template switching between ectopic repeats associated with impeded replication forks [8], [14].

Here, we used the RTS1-RFB to investigate the consequences of fork collapse on genome instability. We report that recovery from a collapsed fork is associated with a high frequency of instability, with a single fork arrest increasing the rates of mutation, deletion and translocation by 10, 40 and 5 fold, respectively. We show that genetic instability associated with fork arrest is dependent on homologous recombination. Fork-arrest-induced GCRs (deletion and translocation) result from inappropriate ectopic recombination at the site of the collapsed fork. We also demonstrate that restoration of fork progression by homologous recombination results in error-prone DNA synthesis due to frequent replication slippage between short tandem repeats. We investigated the molecular mechanisms of this replication slippage and found that post-replication repair, including ubiquitination of PCNA or trans-lesions-synthesis (TLS) DNA polymerases, is not involved in fork-arrest-induced replication slippage. Micro-deletions/insertions flanked by micro-homology associated with copy number variations (CNVs) in cancer cells or in response to replication stress may therefore be scars left following the restoration of forks progression by homologous recombination.

Results

The conditional replication fork barrier RTS1

We generated fork arrest constructs by manipulating the polar RTS1-RFB (Figure 1A). We introduced the RTS1 sequence on the centromere-proximal (cen-proximal) side of the ura4 locus, 5 kb away from the strong replication origin (ori) 3006/7 on chromosome III. This created the t-ura4”refer to the RTS1-barrier and its polarity that is whether it blocks replication forks travelling from the ori 3006/7 towards the telomere or forks travelling from the telomere towards the ori 3006/7, respectively. We previously confirmed that forks moving from ori 3006/7 towards the telomere (tel) are efficiently blocked by the RTS1-RFB at the t-ura435]. In this model system, fork arrest is activated by inducing the expression of rtf1 gene that is under control of the thiamine repressible promoter nmt41. Thus, the RTS1-RFB is inactivated by adding thiamine to the media and it is activated in thiamine-free media. Efficient induction of Rtf1 expression requires incubation for 12–16 hours in thiamine-free media. Replication intermediates were analysed by native 2-dimensional gel electrophoresis (2DGE). In conditions of Rtf1 expression, more than 95% of replication forks were blocked by the RTS1-RFB at the t-ura4ura4Rtf1 induction (Figure 1B, t-ura420].

Figure 1

Conditional replication fork-arrest assays.

A. Diagrams of fork-arrest constructs. Centromere-proximal and telomere-proximal regions are represented in black and grey, respectively. Strong or putative replication origins (ori) and the centromere are indicated by yellow, green and black circles, respectively. Blues arrows indicate the polarity of the RTS1-RFB. The ura4 gene is indicated in red and the arrow indicates its direction of transcription. Representations of the primary arrested fork structure are given for each construct. The name of each fork-arrest construct is given using the following nomenclature: “t” and “ori” refer to the telomere and the replication origin 3006/7, respectively; “<“ and ”>” indicate the RTS1-barrier and its polarity (< blocks replication forks moving from the ori3006/7 towards the telomere, and > blocks replication forks moving from the telomere towards the origin 3006/7. B. Diagrams of replication intermediates (RIs) within the AseI fragment analysed by 2DGE (top panel). Representative RIs analysed by 2DGE in indicated strains in OFF (Rtf1 being repressed) and ON (Rtf1 being expressed) conditions. Signal corresponding to arrested forks, joints-molecules (JMs) and termination structures are indicated by black, red and green arrows, respectively. Note that the t>ura4-ori construct does not result in a strong fork arrest as the RTS1-RFB is not orientated in the main direction of replication (see text for details).

Conditional fork-arrest constructs to investigate fork-arrest-induced genome instability

The RTS1 sequence was inserted on the tel-proximal side of ura4 creating the t<ura4-ori locus. 2DGE analysis of this construct revealed a strong fork arrest signal on the descending large Y arc (Figure 1A and 1B, t<ura4-ori ON). The ura4 gene, used in this system as a reporter to score genetic instability, is located behind the arrested fork when the RTS1-RFB is active at the t<ura4-ori locus and ahead of the arrested fork at the t-ura4RTS1 sequence at the tel-proximal side of ura4 created the t>ura4-ori locus and no fork arrest signal was detected for this construct by 2DGE when Rtf1 was expressed (Figure 1A and 1B, t>ura4-ori ON). Thus, RTS1 behaves as a polar RFB at the ura4 locus, and replication across this locus is strongly unidirectional due to the relative positions of the origins [42].

Introducing a second RTS1 sequence, such that the two RTS1 sequences are inverted repeats (IRs), created t>ura4ura4>ori loci (Figure 1A and 1B, t>ura4ura4>ori ON). Given the orientation of the polar RTS1-RFB in the t<ura4>ori strain, converging forks cannot be blocked. Whereas block of converging forks can virtually occur in the t>ura4cen-proximal side of ura4 were efficiently recovered by recombination before forks are arrested on the tel-proximal side. Indeed, joint-molecules (JMs) resulting from recombination between RTS1 repeats were detected by 2DGE (see red arrows on Figure 1B, t>ura4ura4>ori ON). Resolution of these JMs gives rise to chromosomal rearrangements [20]. In the absence of homologous recombination (i.e. in a rad22-d mutant), JMs were not detected and termination signals accumulated (see green arrow on Figure 1B, t>ura4ura4ura4RTS1-RFB behaves as a hot spot for replication termination [20].

A single fork arrest induces genomic deletions

We investigated fork-arrest-induced genome instability by selecting for cell resistance to 5-FOAR, the result of loss of ura4 function. Inducing fork-arrest at t-ura4ura4 loss 3 fold (Table 1). Rtf1 expression in the t-ura4-ori and t>ura4-ori strains did not cause site-specific fork-arrest at ura4 as assessed by 2DGE and did not increase the rate of ura4 loss. Thus, ura4 loss results from the RTS1-RFB activity and not simply from the presence of RTS1 and/or Rtf1 expression (Table 1). To investigate the nature of this genetic instability, primers were designed to amplify the ura4 coding sequence and, as a control, the essential rng3 gene, mapping 30 kb tel-proximal to ura4, that should not be rearranged (Figure 2A and 2B) [35]. The absence of ura4 amplification was classified as a deletion event; sequencing of amplified ura4 sequence was used to identify point mutation events (Figure 2B).

Table 1

Rates of ura4 loss (including genomic deletion, translocation, and mutation events), calculated using the method of the median.

strainsa	Rate of ura4 lossb		Fold induction by Rtf1 expression (−/+ thiamine)c	Fold induction by the RFB (relative to the t-ura4-ori construct)
	+ thiamine (Rtf1 repressed)	− thiamine (Rtf1 expressed)
t-ura4-ori	4.6±1	4.3±1.8	0.9
t-ura4<ori	4.1±0.9	13.9±2.8	3.4 (p = 0.002)	3.2
t-ura4<ori RTS1-d	4.9±1.4	11.4±4.1	2.3 (p = 0.03)	2.6
t<ura4-ori	3.9±1.3	4.8±2.5	1.2	1.1
t>ura4<ori	5.8±2.1	97.4±51.5	16.8 (p = 0.004)	22.6
t>ura4<ori RTS1-d	4.9±2.2	19.5±3.3	4 (p = 0.006)	4.5
t>ura4-ori	3.2±1.7	3.6±1.8	1.1	0.8
t<ura4>ori	6.6±1.7	50.3±14.2	7.6 (p = 0.014)	11.7

the following nomenclature is used: “t” and “ori” refer to the telomere and the replication origin 3006/7, respectively; “<“ and ”>” indicate the RTS1-RFB and its polarity (< blocks replication forks moving from the ori3006/7 towards the telomere, and > blocks replication forks moving from the telomere towards the origin 3006/7).

event/cell/division ×10−8 ± standard error. Values are means of at least 3 independent rates.

statistical significance was determined using the nonparametric Mann-Whitney U test.

Figure 2

A single fork-arrest induces GCRs that are stimulated by inverted repeats near the site of fork arrest.

A. Diagrams of chromosome II containing or not the RTS1 sequence (blue arrow or RTS1-d) and of chromosome III containing ura4 alone or associated with RTS1-RFB constructs. The RTS1 sequence maps near the mat1 locus where it helps to ensure unidirectional replication [62]. Primers used for amplifying the 1 Kb ura4 fragment or the 650 bp rng3 fragment are depicted in red and grey, respectively. Primers used to amplify the translocation junction (1.2 kb) are represented in orange on chromosome II (TLII) and in black on chromosome III (TLIII). B. Representative PCR-amplifications from 5-FOAR colonies of the indicated strains; ON and OFF refers to the RTS1-RFB being active or not, respectively. PCR products and their sizes are indicated on the figure. C. Effect of intra- and inter-chromosomal recombination between RTS1 repeats on fork-arrest-induced genomic deletion. RTS1-RFB activity and ura4 location with respect to the RFB are given for each construct. The % of deletion events, as determined by the PCR assay, was used to balance the rate of ura4 loss. Then, the RFB-induced deletion rate was calculated by subtracting the rate obtained in the presence of thiamine (Rtf1 being repressed) from the rate obtained in the absence of thiamine (Rtf1 being expressed). The values reported are means of at least 3 independent median rates. Error bars correspond to the standard error (SE). D. Effect of Rqh1 on RFB-induced deletions (left) and translocations (right), as described for panel C. Error bars indicate SE. Statistically significant fold differences between the rqh1-d and the wild-type strains are indicated with an *. E. Representative PCR amplifications from 5-FOAR colonies of the rqh1-d t-ura4

A. Diagrams of chromosome II containing or not the RTS1 sequence (blue arrow or RTS1-d) and of chromosome III containing ura4 alone or associated with RTS1-RFB constructs. The RTS1 sequence maps near the mat1 locus where it helps to ensure unidirectional replication [62]. Primers used for amplifying the 1 Kb ura4 fragment or the 650 bp rng3 fragment are depicted in red and grey, respectively. Primers used to amplify the translocation junction (1.2 kb) are represented in orange on chromosome II (TLII) and in black on chromosome III (TLIII). B. Representative PCR-amplifications from 5-FOAR colonies of the indicated strains; ON and OFF refers to the RTS1-RFB being active or not, respectively. PCR products and their sizes are indicated on the figure. C. Effect of intra- and inter-chromosomal recombination between RTS1 repeats on fork-arrest-induced genomic deletion. RTS1-RFB activity and ura4 location with respect to the RFB are given for each construct. The % of deletion events, as determined by the PCR assay, was used to balance the rate of ura4 loss. Then, the RFB-induced deletion rate was calculated by subtracting the rate obtained in the presence of thiamine (Rtf1 being repressed) from the rate obtained in the absence of thiamine (Rtf1 being expressed). The values reported are means of at least 3 independent median rates. Error bars correspond to the standard error (SE). D. Effect of Rqh1 on RFB-induced deletions (left) and translocations (right), as described for panel C. Error bars indicate SE. Statistically significant fold differences between the rqh1-d and the wild-type strains are indicated with an *. E. Representative PCR amplifications from 5-FOAR colonies of the rqh1-d t-ura4Rtf1 is expressed or not).

A single arrested fork at the t-ura4ura4-ori strain, p = 0.006) (Figure 2C and Figure S1A). Fork-arrest-induced deletion was recombination-dependent. Spontaneously (i.e. when the RTS1-barrier was inactivated), the rate of genomic deletion in rad22-d or rhp51-d strains was higher than that in the wild-type strain (Figure S1B). Nonetheless, no further increase in the rate of genomic deletion was observed in the surviving rad22-d or rhp51-d cells upon activation of the RTS1-barrier (Figure S1B, t-ura414], [43], [44]. Deleting the natural RTS1 sequence from chromosome II abolished deletion events at collapsed forks, indicating that fork-arrest-induced deletion was also mediated by inter-chromosomal recombination (Figure 2C and t-ura4RTS1-d on Figure S1A). Thus, these data are consistent with the view that homologous recombination makes a major contribution to suppressing genome instability, but can occasionally drive non allelic recombination events leading to GCRs [35], [45].

We detected no fork-arrest-induced deletion in the t<ura4-ori strain, in contrast to the t-ura4ura4 marker is located behind and ahead of collapsed forks in the t<ura4-ori and t-ura4

Inverted repeats stimulate fork-arrest-induced deletion by promoting inter- and intra-chromosomal recombination

Inverted repeats (IRs) are structural elements often associated with genome rearrangements [11], [14], [46], [47]. We investigated the effects of IRs in the vicinity of the RTS1-RFB on fork-arrest-induced genomic deletion. We first compared the t>ura4ura4) to the t-ura4RTS1-RFB). The rate of fork-arrest-induced genomic deletion was 200 times higher in the t>ura4ura4RTS1 sequence on the tel-proximal side of ura4 accounted for 99.5% of the genomic deletions observed in the t>ura4ura4>ori). Preventing inter-chromosomal recombination by deleting RTS1 from the chromosome II (t>ura4RTS1-d) abolished 90% of deletion events (Figure 2C and Figure S1A). Thus, genomic deletions induced by fork-arrest near IRs are due to inter- and intra-chromosomal recombination events. In support of this, stimulation of fork-arrest-induced deletion by IRs is mediated by homologous recombination. Indeed, the rate of genomic deletion was not increased upon induction of the RTS1-RFB in the surviving population of t>ura4

We verified that our data were not influenced by the orientation of IRs or by rare blocking of converging forks in the t>ura4ura4>ori construct in which RTS1 repeats are in the opposite orientations relative to the t>ura4ura4 cannot be blocked (Figure 1). The rate of fork-arrest-induced genomic deletion was 1,000 times higher in the t<ura4>ori than that in the t<ura4-ori strain, that does not contain IRs near the RTS1-RFB (p = 0.008, Figure 2C and Figure S1A). Thus, intra-chromosomal recombination, permitted by the RTS1-RFB sequence on the cen-proximal side of ura4, accounted for nearly 100% of the genomic deletions observed in the t<ura4>ori strain (Figure 2C, compare with t<ura4-ori). Preventing inter-chromosomal recombination by deleting RTS1 from the chromosome II (t<ura4>ori RTS1-d) abolished 90% of deletion events (Figure 2C and Figure S1A). Importantly, the deletion rates for the t<ura4>ori and t>ura4

A single collapsed fork induces translocations that are stimulated by IRs

Fork-arrest at t>ura4RTS1 repeats on chromosomes II and III [35]. We investigated the influence of IRs on fork-arrest-induced translocation. We designed primers to amplify the predicted translocation junction between chromosomes II and III (TLII and TLIII on Figure 2A and 2B). A single arrested fork at the t-ura4ura4ura4RTS1-RFB (p = 0.009, Figure 2D and Figure S1C). Thus, intra-chromosomal recombination accounted for nearly 99% of translocation events observed in the t>ura4ura4RTS1 from the chromosome II (t>ura4RTS1-d on Figure 2B). Therefore, as reported for genomic deletions, fork-arrest-induced translocation associated with IRs is due to inter- and intra-chromosomal ectopic recombination. No translocations were detected in the t<ura4>ori strain (data not shown), so we cannot formally exclude the possibility that fork-arrest-induced translocations in the t>ura4Rad22Rad52 or Rhp51Rad51, it is most likely that translocations occur during fork recovery by recombination (Figure 2B and [35]). Overall, our data indicate that recovery of a single collapsed fork causes translocations and IRs near the site of fork-arrest stimulate translocations by promoting inappropriate inter- and intra-chromosomal recombination.

Fork-arrest-induced GCRs are caused by inter- and intra-chromosomal recombination. We noticed a slightly greater contribution of intra- than inter-chromosomal recombination (Figure 2C). This is consistent with ectopic recombination preferentially occurring at the most proximal homologous sequence, as previously reported [48]. Nonetheless, the rate of fork-arrest-induced deletion in the t>ura4ura4RTS1-d strain) and inter-chromosomal recombination events (4 10−9 in the t-ura4>ori strain). Similar reasoning can be applied for the t<ura4>ori strain. Thus, independent intra- and inter-recombination events cannot themselves explain high rate of GCRs induced by fork arrest near IRs. Therefore, we infer that there is interplay between inter- and intra-chromosomal recombination such that fork-arrest-induced GCRs may involve recombination between three homologous sequences (tri-parental recombination).

The RecQ helicase Rqh1 prevents GCRs at collapsed forks

To confirm that fork-arrest-induced GCRs are the result of inappropriate ectopic recombination during fork recovery, we analysed the involvement of the RecQ helicase Rqh1. We previously reported that Rqh1 limits inappropriate template switching of stalled nascent strands without affecting the efficiency of fork restart [20]. In the t-ura4ura4

A single fork arrest induces mutations

We analysed the effects of collapsed forks on the mutation rate. We sequenced the ura4 coding sequence from 5-FOAR isolated cells and identified base-substitutions, frame-shifts and small insertions and duplications between short tandem repeats (Table 2). A single collapsed fork in the t-ura4RTS1 deleted from chromosome II (Figure 3A and Figure S2A). Thus, fork-arrest-induced mutation is not mediated by inappropriate ectopic recombination. Induction of the RTS1-RFB in the t<ura4-ori strain did not increase the mutation rate of the ura4 gene. Thus, as for GCRs, replicated regions behind arrested forks are not prone to mutation. This observation rules out the hypothesis that fork-arrest-induced mutation is a consequence of the accumulation of damaged single-stranded DNA behind collapsed forks (see discussion). Our data suggest that recovery from collapsed forks results in error-prone DNA-synthesis.

Table 2

Mutations spectra in the indicated strains.

	t-ura4-ori a		t>ura4-ori		t-ura4<ori		t>ura4<ori		t<ura4-ori		t-ura4<ori rad22-d
Rtf1 expression	−	+	−	+	−	+	−	+	−	+	−	+
Fork arrest at ura4	−	−	−	−	−	+	−	+	−	+	−	+
Transitionsb
GC→AT	5			1	7	14	6	2		5	5	5
AT→GC	3		1	2	4	9	1	4	4	7	4	6
Transversionsb
GC→TA	3	5	1		45	20	5	2	4	2	3	4
GC→CG	1	3	4	2	2	3	2			1	8	5
AT→CG	3	3	1	1	6	4	2	2		1	4
AT→TA		1	5	3	4	2	4	3	3	1	2	3
Rate of base substitutionsc	4.6±1.1	3.7±1.5	3.3±1.6	4.2±1.3	3.4±0.9	6.4±1.2d	2.6±1.1	6.5±1.5d	2.5±1.3	3.3±1.7	115±34	152±33
+ Frame-shiftb:
in homo-nt		1							1
in mono-nt		1	2	1	2	1	1		1		1
- Frame-shiftb:
in homo-nt					1	4		1			1	2
in mono-nt					1	3			1			1
Rate of frame-shiftc		0.6±0.25	0.5±0.3	0.5±0.1	0.2±0.06	1±0.2d	0.1±0.06	0.5±0.1	0.7±0.4		8.7±2.6	20.4±4.4
Del/Duplb:
Deletions					9	27		4		1		2
Duplications					5	12
Rate of Del/Duplc		<0.3		<0.5	0.7±0.2	4.9±1	<0.1	2.1±0.5	<0.2	0.2±0.1	<3.7	13±2.8
Other (complexe)					1
Total	15	14	14	10	87	99	21	18	14	18	28	28

the following nomenclature is used: “t” and “ori” refer to the telomere and the replication origin 3006/7, respectively; “<“ and ”>” indicate the RTS1-RFB and its polarity (< blocks replication forks moving from the ori3006/7 towards the telomere, and > blocks replication forks moving from the telomere towards the origin 3006/7).

number of events.

event/cell/division ×10−8 ± standard error. The % of each mutation events was used to balance the rate of mutation.

value not significantly different from those for the t-ura4-ori strain (spontaneous events).

Figure 3

Fork-arrest induces mutations.

A. Effect of intra- and inter-chromosomal recombination between RTS1 repeats on fork-arrest-induced mutation rates (base-substitutions, frame-shifts and small insertions or deletions between short tandem repeats). RTS1-RFB activity and ura4 location with respect to the RFB are given for each construct. The % of mutation events, as determined by the PCR assay and sequencing, was used to balance the rate of ura4 loss. Then, the RFB-induced mutation rate was calculated by subtracting the rate obtained in the presence of thiamine (Rtf1 being repressed) from the rate obtained in the absence of thiamine (Rtf1 being expressed). The values reported are means of at least 3 independent median rates. Error bars correspond to SE. (Refer to Figure S2 for corresponding rates of mutation when Rtf1 is expressed or not). B. Rate of mutation for indicated strains; ON and OFF refers to the RTS1-RFB being active or not, respectively. The % of mutation events, as determined by the PCR assay and sequencing, was used to balance the rate of ura4 loss. The values reported are means of at least 3 independent median rates. Error bars correspond to SE. Statistically significant fold differences in mutation rates between the “OFF” and “ON” conditions are indicated with an *. C. Spectra of mutation events in indicated strains upon RFB induction (refer to Table 2 for exact numbers and to Figure S3 for mapping of deletions/duplications and their features). D. Strains harbouring the ura4 alleles with a single base-substitution or frame-shift or duplication of 20 or 22 nt, together with the RTS1-RFB in the t-ura4

A. Effect of intra- and inter-chromosomal recombination between RTS1 repeats on fork-arrest-induced mutation rates (base-substitutions, frame-shifts and small insertions or deletions between short tandem repeats). RTS1-RFB activity and ura4 location with respect to the RFB are given for each construct. The % of mutation events, as determined by the PCR assay and sequencing, was used to balance the rate of ura4 loss. Then, the RFB-induced mutation rate was calculated by subtracting the rate obtained in the presence of thiamine (Rtf1 being repressed) from the rate obtained in the absence of thiamine (Rtf1 being expressed). The values reported are means of at least 3 independent median rates. Error bars correspond to SE. (Refer to Figure S2 for corresponding rates of mutation when Rtf1 is expressed or not). B. Rate of mutation for indicated strains; ON and OFF refers to the RTS1-RFB being active or not, respectively. The % of mutation events, as determined by the PCR assay and sequencing, was used to balance the rate of ura4 loss. The values reported are means of at least 3 independent median rates. Error bars correspond to SE. Statistically significant fold differences in mutation rates between the “OFF” and “ON” conditions are indicated with an *. C. Spectra of mutation events in indicated strains upon RFB induction (refer to Table 2 for exact numbers and to Figure S3 for mapping of deletions/duplications and their features). D. Strains harbouring the ura4 alleles with a single base-substitution or frame-shift or duplication of 20 or 22 nt, together with the RTS1-RFB in the t-ura4RFB “OFF”) or without (RFB “ON”) thiamine. The bottom diagram indicates strain positions and the mutation events required to obtain Ura+ revertants. E. PCR analysis of Ura+ revertants isolated from the t-ura4-dup20ura4) after RFB induction. With the primers used, a 106 bp fragment is amplified from the ura4 strain and a 126 bp fragment is amplified from the t-ura4-dup20ura4-dup22, ura4-dup20, ura4 alleles and corresponding Ura+ revertants. Micro-homologies are indicated in grey and duplicated sequences are underlined in black. The phenotype of each allele is indicated on the figure. G. Map of deletion and duplication events within the ura4 ORF.

Collapsed forks result specifically in replication slippage

We then analysed the spectra of mutations found in the ura4 ORF by sequencing the PCR products. The rates of base-substitutions and frame-shifts were not significantly increased by the RTS1-RFB activity over spontaneous events (i.e. compare to t-ura4-ori strain, Figure 3C and Table 2). In contrast, the rate of deletions and duplications (Del/Dup) flanked by short homology was increased by 7 times over spontaneous events in the t-ura4ura4-ori strain (Figure 3C and Table 1). These data further confirm that fork-arrest does not promote mutation events behind collapsed forks.

We used reversed mutation assays to test if fork-arrest at the RTS1-RFB specifically induced Del/Dup mutations. We made use of strains harbouring a single mutation within the ura4 ORF: either a single base-substitution or a −1 frame-shift in homo-nucleotide (Figure S2B). We also studied strains harbouring either a duplication of 20 or 22 nt flanked by 5 or 4 bp of micro-homology, respectively (defined as ura4-dup20 and ura4-dup22, Figure S2B). These non-functional ura4 alleles were inserted in front of the RTS1-RFB in the t-ura4ura4 + gene. Activation of the RTS1-RFB at ura4 increased the frequency of Ura+ revertants up to 15 and 7 times in strains harbouring ura4-dup22 and ura4-dup20, respectively (Figure 3D and Figure S2B). Thirty Ura+ colonies were studied by PCR and all gave a product of the same size as the wild-type ura4 gene: they had therefore lost the duplication (Figure 3E and data not shown). Sequencing the full ura4 ORF confirmed that Ura+ revertants contained an intact ura4 sequence, showing that the reversion of these alleles was due solely to the precise deletion of 20 or 22 nt (Figure 3F and data not shown). In contrast, activation of the RTS1-RFB did not increase the frequency of Ura+ revertants of strains harbouring ura4 alleles with a single base-substitution or a −1 frame-shift (Figure 3D and Figure S2B). Thus, collapsed forks tend to induce deletion events between short tandem repeats rather than base-substitution or frame-shift mutations.

Among Del/Dup events, deletions represented the two-third of events in the t-ura4ura4-dup20 and ura4-dup22 alleles used in the reverse mutation assay were representative of the Del/Dup events observed. Del/Dup flanked by micro-homology result from intra-molecular template switching mechanisms in which nascent strands dissociate from the template and misalign with the template when restarting the elongation step. This leads to loop formation, either in the nascent strand or in the template, resulting in duplication or deletion events, respectively [49]. Consequently, we will hereafter refer to Del/Dup as replication slippage. Replication slippage was observed all along the ura4 ORF and up to 1.2 kb ahead of the arrested fork, even if a hot spot of deletion was present 500 bp away from the RTS1-RFB (Figure 3G and Figure S3B). Thus, our data suggest that the DNA synthesis is prone to replication slippage at least for the first 1,200 nt synthetized during the recovery of collapsed forks. Inaccuracy of DNA synthesis on further distances was not directly addressed.

Replication slippage results from error-prone DNA synthesis during fork recovery

To confirm that replication slippage occurs as forks recover, and not behind the fork in the DNA already replicated, we inserted the ura4-dup20 or the ura4-dup22 allele either behind (in the t<ura4-ori configuration) or in front of the RTS1-RFB (in the t-ura4ura4RTS1-RFB resulted in a 8 and 16 fold increases in the replication slippage frequency for the ura4-dup20 and ura4-dup22 alleles, respectively (Figure 4A and 4B). Similar increases in the rate of replication slippages were observed (Figure 4C). In contrast, in the t<ura4-ori background, the frequency of replication slippage was induced by only 2–3 fold by the RTS1-RFB (Figure 4B–4C). These data confirm that DNA located ahead of collapsed forks is more prone to replication slippage than replicated DNA adjacent to arrested forks, further evidence that replication slippage arises during fork recovery.

Figure 4

Fork recovery by homologous recombination results in replication slippage.

A. Serial tenfold-dilutions from indicated strains (t-ura4-dup20-ori associated or not with the RTS1-RFB) were spotted onto the indicated media after cell growth with (Rtf1 -, repressed) or without (Rtf1 +, expressed) thiamine. RTS1-RFB activity is given for each construct and condition. B. Frequency of Ura+ revertants in indicated strains after cell growth with (Rtf1 repressed) or without (Rtf1 expressed) thiamine. The RTS1-RFB activity is given for each construct and condition. Values correspond to the mean of at least three independent experiments and error bars correspond to the standard error of the mean (SEM). C. Rate of replication slippage in the indicated strains and conditions. The rate of Ura+ revertants was calculated using the method of the median from at least 11 independent cultures. Values in brackets indicate the 95% confidence interval. Statistical significance was detected using the nonparametric Mann-Whitney U test. D. Serial tenfold-dilutions from the strains indicated spotted onto the media indicated after cell growth without thiamine. RTS1-RFB activity “–” refers to the strain t-ura4-dup20-ori and “+” refers to the strain t-ura4-dup20

A. Serial tenfold-dilutions from indicated strains (t-ura4-dup20-ori associated or not with the RTS1-RFB) were spotted onto the indicated media after cell growth with (Rtf1 -, repressed) or without (Rtf1 +, expressed) thiamine. RTS1-RFB activity is given for each construct and condition. B. Frequency of Ura+ revertants in indicated strains after cell growth with (Rtf1 repressed) or without (Rtf1 expressed) thiamine. The RTS1-RFB activity is given for each construct and condition. Values correspond to the mean of at least three independent experiments and error bars correspond to the standard error of the mean (SEM). C. Rate of replication slippage in the indicated strains and conditions. The rate of Ura+ revertants was calculated using the method of the median from at least 11 independent cultures. Values in brackets indicate the 95% confidence interval. Statistical significance was detected using the nonparametric Mann-Whitney U test. D. Serial tenfold-dilutions from the strains indicated spotted onto the media indicated after cell growth without thiamine. RTS1-RFB activity “–” refers to the strain t-ura4-dup20-ori and “+” refers to the strain t-ura4-dup20Ura+ revertants frequency for the strains indicated as a function of the number of generations after thiamine removal. RTS1-RFB activity “–” refers to the strain t-ura4-dup20-ori and “+” refers to the strain t-ura4-dup20ura4-dup20ura4-dup20-ori) the active RTS1-RFB. The rate was calculated from the slope of the curves presented in panel F. The values reported are means of three independent experiments and error bars correspond to SE.

Replication slippage results from forks restarted by recombination

Replication slippage occurs in DNA in front of (and not behind) the arrested fork, this DNA being replicated only after restart of the fork. Thus, a defect preventing fork recovery would be expected to abolish the error-prone DNA synthesis during restart. We analyzed fork-arrest-induced mutation in recombination mutants in which collapsed forks at the RTS1-RFB cannot recover, resulting in cell death. Induction of the RTS1-RFB did not increase the overall mutation rate in the surviving populations of t>ura4ura4ura4ura4-dup20 allele, to determine the rate of replication slippage induced by the RTS1-RFB over spontaneous events.

Strains carrying mutations in recombination genes grow slowly, so replication slippage was scored as a function of the number of generations following thiamine removal (i.e. generations subject to fork arrest at ura4) (Figure 4D and 4E). In the wild-type strain, fork arrest at the RTS1-RFB resulted in a 10 fold-increase in the frequency of replication slippage, as expected. In recombination mutants (rad50-d, rhp51-d and rad22-d), fork-arrest at the RTS1-RFB increased the frequency of replication slippage by only 2 times over spontaneous events: therefore, replication slippage occurs less frequently in survivors from recombination mutants than those from the wild-type strain (Figure 4D–4F). Based on 2DGE analysis, fork-restart is severely impaired in the absence of Rad22Rad52 (Figure 1B and [20]), such that even the two-fold induction in replication slippage by fork arrest in the rad22-d strain was surprising. The rad22-d strain accumulates suppressors involving the Fbh1 helicase that limits Rhp51Rad51- dependent recombination at blocked replication forks [50], [51]. Therefore, we analyzed replication slippage in the rad22-d rhp51-d double mutant in which no homologous recombination event occurs. In this background, there was no detectable fork-arrest-induced replication slippage (Figure 4D–4F). Thus, complete defect in fork restart results in a complete abolition of fork-arrest-induced replication slippage in the surviving population. Overall, our data establish that replication slippage results from inaccurate DNA synthesis during the restart of collapsed forks by recombination.

Replication stress leading to fork collapse induces replication slippage

We investigated the effects of replication stress, other than the replication block imposed by the RTS1-RFB, on replication slippage. Strains harbouring ura4 alleles (base-substitutions, −1 frame-shift, and ura4-dup20) were exposed to replication-blocking agents or UV-C-induced DNA damages and the frequency of Ura+ revertants was scored. Three hours of treatment with either the topoisomerase I inhibitor camptothecin (CPT) or mitomycin C (MMC), an inter-strand cross-linking agent (ICls), increased the frequency of Ura+ revertants by 3 to 4 fold in the ura4-dup20 strain (Figure 5A and 5B). At equivalent survival (70–90%), DNA-damages induced by a dose of 100 J/m2 of UV-C did not increase the frequency of Ura+ revertants in the ura4-dup20 strain. Increasing the UV-C dose (150 J/m2) resulted in an increased reversion effect. The other ura4 alleles exhibited an opposite behaviour pattern. As expected, UV-C-induced DNA damages, but not CPT or MMC treatment, increased the frequency of Ura+ revertants of the base-substitution and the −1 frame-shift mutants (Figure 5A). Thus, replication slippage, unlike other point mutations, appears to be a mutation event specifically induced by replication stress.

Figure 5

Collapsed forks, but not stalled forks, induce replication slippage.

A. Left panel: the frequency of Ura+ revertants as a function of time-contact with indicated drugs for the indicated ura4 alleles (single base-substitution, frame-shift, duplication of 20 nt). Right panel: the frequency of Ura+ revertants in response to UV-C irradiation as a function of dose for the indicated ura4 alleles. The values reported are means of two independent experiments. Numbers indicate fold difference in the frequency of Ura+ revertants between the treated and untreated control samples. For ura4 alleles containing base-substitutions or frame-shifts, the mutation event required to obtain Ura+ revertants is indicated on the figure. B. Serial tenfold-dilutions from ura4-dup20 strain spotted onto the media indicated after treatment with MMC or CPT as indicated. C. Frequency of Ura+ revertants after the indicated treatments in the ura4-dup20 strain. DMSO (the vehicle) was used as control for CPT treatment. The values reported are means of at least three independent experiments. Error bars correspond to SEM.

Hydroxyurea (HU) that prevents the bulk of dNTP synthesis during S-phase by inhibiting the ribonucleotide reductase, results in a slow-down of fork progression which did not induce replication slippage (Figure 5A). In contrast, CPT and MMC treatments that lead to replication stress by causing fork collapse induced replication slippage. Homologous recombination is repressed during HU treatment and recombination proteins are recruited to collapsed but not stalled forks [52]–[54]. Consistent with this, we found that the rad22-d mutant is highly sensitive to acute exposure to CPT, but not to HU (Figure S4). Thus, acute exposure to HU results in stalled forks that recover without recombination, whereas recombination may be required for restarting forks that have collapsed due to CPT or MMC treatment. We confirmed that CPT-induced replication slippage results from collapsed forks and was thus S-phase specific: the ura4-dup20 strain was synchronized in early S-phase by HU treatment and released into S-phase with or without CPT. HU-synchronization and release into DMSO (used as vehicle for CPT) did not induce replication slippage. In contrast, the release of cells into S-phase in the presence of CPT stimulated replication slippage up to 12 fold (Figure 5C). These data indicate that CPT-induced fork collapse results in error-prone DNA synthesis characterized by replication slippage. These experiments further support the view that replication slippage results from recovery of collapsed forks by recombination and point out that the RTS1-barrier is representative of collapsed forks restarted by homologous recombination.

TLS–DNA polymerases are not involved in fork-arrest-induced replication slippage

To investigate the inaccuracy of the DNA synthesis occurring immediately following the restart of collapsed forks, we analysed the involvement of TLS-DNA polymerases. In fission yeast, TLS pathways require either mono- or poly-ubiquitination of the clamp loader PCNA on lysine 164 [55]. We found that mutating this lysine to arginine residue did not affect replication slippage induced by the RTS1-RFB. None of Rev1, Rev3 or DinB DNA polymerases were required for fork-arrest-induced replication slippage (Figure S5). Therefore, the error-prone DNA synthesis associated with fork recovery by recombination does not rely on TLS DNA polymerases activity.

The mismatch repair pathway does not counteract fork-arrest-induced replication slippage

The mismatch repair (MMR) pathway is temporally coupled to DNA replication, and MMR components are associated with replication centres [56]. The heterodimer Msh2/Msh6 recognises mispaired DNA and Msh2/Msh3 recognises small DNA loops up to 31 bases long, arising from replication slippage [57]. The failure to repair small DNA loops results in more frequent insertions and deletions [58]. Therefore, MMR activity could potentially lead to an underestimation of the extent of fork-arrest-induced replication slippage. However, replication slippage induced by the RTS1-RFB activity was as frequent in msh2-d, msh6-d and msh3-d strains as in wild-type control. Also, spontaneous replication slippage at ura4 (without RTS1-RFB) was unaffected by the absence of MMR proteins (Figure S5). Therefore, there is no evidence that MMR repairs small DNA loops (20 nt) in fission yeast and fork-arrest-induced replication slippage is not counteracted by MMR in our model system.

Discussion

Using conditional fork arrest constructs, we studied the consequences for genome instability of impediments to replication forks progression. A single fork arrest results in large-scale genomic changes and mutations that occur during recombination-mediated fork recovery (Figure 6). Inappropriate ectopic recombination at arrested forks results in GCRs, whereas appropriate restarting of the fork on the initial template results in error-prone DNA synthesis. GCRs and mutations at collapsed forks are genetically separable: Rqh1 limits fork-arrest-induced GCRs but not mutations (Figure 2D and Figure 3B). We demonstrate here that collapsed forks whose progression resumes by recombination lose accuracy during DNA synthesis, resulting in frequent intra-template switches. Thus, homologous recombination contributes to completion of DNA replication when forks progression is impeded but also fuels genome modifications both at the chromosomal and nucleotide level.

Figure 6

Model of replication-stress-induced genetic instability at collapsed forks.

Collapsed forks might arise from torsional stress, fork breakage (i.e. at nick, ICLs), or proteins tightly-bound to DNA. Replisome disassembly at collapsed forks may favour the unwinding of the nascent strand on which Rad51 nucleates. At this initial stage of fork resumption by recombination, homology-driven template exchange can promote intra- or inter-recombination resulting in GCRs. Fork recovery by recombination overcomes the initial replication block and allows an inaccurate replisome to form (see text for details).

Non allelic homologous recombination (NAHR) between low copy number repeats (LCR) contributes to recombination-mediated GCRs in mitosis and meiosis. NAHR is responsible for translocations, deletions, inversions and loss of heterozygosity [59]. Ou et al. predicted 1143 LCR pairs in the human genome liable to mediate recurrent translocations [60]. In budding yeast, a single DSB is sufficient to mediate recombination-dependent translocation [61]. Here, we report that a single collapsed fork increases the rate of genomic deletion 40 fold, and that of translocation 5 fold. Fork-arrest-induced GCRs are mediated by NAHR between heterologous chromosomes. It is not clear whether fork arrest on both homologous repeats contributes to fork-arrest-induced GCRs. We could not address this question in our model system, because the RTS1 sequence close to the mat1 locus on chromosome II has a low RFB activity [62]. Also, the recruitment of recombination proteins at the RTS1 sequence near the mat1 locus is not regulated by the level of Rtf1 expression, showing that regulating Rtf1 expression was in itself insufficient to regulate the RTS1-RFB activity at the mat1 locus [35].

Inverted repeats (IRs) are structural elements that contribute to genome instability. Impediments to replication forks progressing through IRs favour their fusion to generate acentric and dicentric inverted chromosomes [11], [14]. IRs in humans can also trigger the formation of inverted genomic segments and complex triplication rearrangements by a replication-based mechanism [47]. Here, we report that IRs near collapsed forks can increase the rate of GCRs by up to 1,500 fold. This high level of GCRs cannot be explained merely by the addition of independent inter- and intra-chromosomal recombination events. Rather, our analyses suggest that IRs may stimulate tri-parental recombination events induced by template switching of nascent strands at collapsed forks, such that three homologous sequences are involved. Similarly, recombination-dependent translocations induced by a single DSB in budding yeast is proposed to be the consequence of tri-parental recombination [63], [64]. One possible mechanism is that IRs-induced GCRs result from successive template switches initiated by nascent strands at collapsed forks, reminiscent of the multiple template switches during break-induced-replication (BIR) in budding yeast [33]. Interestingly, Rqh1 prevents fork-arrest-induced GCRs by limiting both inter- and intra-chromosomal recombination without affecting fork restart efficiency. Thus, tri-parental recombination might correspond to multiple and successive template switches between homologous repeats.

The high accuracy of DNA replication is compromised by impediments to fork progression, and recombination-mediated fork recovery results in decreased processivity of DNA synthesis (Figure 6). Recombination-induced mutations associated with DSBs or impeding DNA replication have been described previously. The formation of damaged single-stranded DNA during the resection of DSBs favours base-substitutions [65]. We detected the recruitment of the single-strand binding protein RPA up to 1.4 kb behind, but not ahead of the RTS1-arrested fork (data to be published), showing that fork-arrest-induced mutation is not correlated with damaged single-stranded DNA exposed behind collapsed forks. Nevertheless, there were rare replication slippage occurred behind RTS1-arrested forks (in the t<ura4-ori construct), suggesting that resumption of DNA synthesis can in some cases occur at a position behind the site of the collapsed fork. Recombination-dependent DNA synthesis occurring outside S-phase is also highly inaccurate during gene conversion and BIR, resulting in either template switches, base-substitutions or frame-shifts [33], [34], [66]. Elevated dNTP pools, due to activation of the DNA-damage checkpoint in G2 cells, contributes to the generation of mutations when hundreds of kbs are synthesised during BIR [34]. Here, we demonstrate that during normal S-phase progression, a single collapsed fork, restored by recombination, results in replication slippage up to 1.2 Kb away from the initial restarting point.

Recombination-induced replication slippage has been reported previously.

In fission yeast, a defect in pol alpha (swi7-H4 mutant) is associated with a recombination-mediated mutator phenotype characterized by an increased frequency of base-substitutions and Del/Dup between short direct repeats [67]. DNA-polymerase kappa (DinB) and zeta (Rev3) are responsible for the increased base-substitution rate, but the DNA-polymerases involved in Del/Dup mutations were not identified [68]. In budding yeast, a defect in polymerase delta (pol3-t mutant) results in an increased level of replication slippage, mediated by homologous recombination [49]. In contrast, the increased rate of replication slippage in the absence of the accessory subunit of polymerase delta, Pol32, does not depend on a functional recombination pathway [69]. Here, we report that recovery of collapsed forks by recombination is specifically associated with replication slippage. Nonetheless, spontaneous replication slippage events are also increased in strains mutated for recombination genes (Figure 4D–4F). Pol32 is required for BIR and replication-induced template switches leading to segmental duplication [70], [71]. Recombination is responsible for only half of these segmental duplications. Thus, it is possible that fork-restart mechanisms dependent on Pol32 and homologous recombination are prone to replication slippage and that in the absence of these pathways alternative micro-homology mediated mechanisms are revealed.

We suggest that at least two steps of the recombination-based fork recovery mechanism can compromise genome stability (Figure 6). At an initial stage, recruitment of recombination proteins on stalled nascent strands favours both fork recovery and ectopic template switches leading to GCRs. At a later stage, once the replisome has been reconstituted and the fork has resumed its progression, the nascent strands are prone to intra-template switching leading to replication slippage. The reasons for the inaccuracy of DNA synthesis associated with restarted forks during scheduled DNA replication (i.e. in S-phase) remain to be determined. One possibility is that one or more factors are missing in the rebuilt replisome during recovery by recombination. Oncogene-induced replication stress results from unbalanced DNA replication that contributes to genome instability in precancerous cells [12]. Completion of DNA replication in such stress conditions is likely to rely on recombination-mediated fork recovery that in turn generates genome instability. Insertions/deletions flanked by micro-homology, responsible for copy number variations (CNVs), have been identified both in cancer cells and also in response to replication inhibition [72], [73]; their reported sizes are between 1 Kb and several tens of mega bases, but the analysis of these features has been limited by the resolution of array-based genomics approaches. Sub-microscopic insertions/deletions flanked by micro-homology have been also described at loci in which replication origins are scarce, including the human fragile site FRA3B, the instability of which is a consequence of replication stress [74]. Interestingly, homologous recombination contributes to the stability of fragile sites by facilitating complete replication or by repairing gaps and breaks at these sites. Thus, we propose that micro-homology-mediated CNVs could be viewed as scars left by error-prone replication forks restarted by recombination.

Materials and Methods

Standard genetic and molecular biology

Strains used were constructed by standard genetic techniques and are listed in Table S1. 2DGE was performed as previously described [20]. To create ura4-dup20 and ura4-dup22 alleles associated or not with the RTS1-RFB, genomic DNA was isolated from selected 5-FOAR cells containing duplications identified by sequencing. A PCR fragment containing duplications within the ura4 ORF was amplified using the following primers: TTCTGTTCCAACA-CCAATGTTT and TCACGTTTATTTTCAAACATCCA. The PCR products were purified and used to transform strains SL206 (ura4), SL350 (t-ura4ura4-ori). Transformants were selected on 5-FOA-containing plates. Appropriate replacement of ura4+ by ura4-dup20 or ura4-dup22 was verified by PCR and sequencing.

ura4 loss assay

A minimum of 11 independent single colonies from appropriate strains growing with or without thiamine were inoculated in 10 ml of non-selective media (with or without thiamine) and grown to stationary phase. Appropriate dilutions were plated on supplemented YEA to determine plating efficiency and on 5-FOA-containing plates. Colonies were counted after 5–7 days of incubation at 30°C. The rate of ura4 loss was determined with the method of the median and data are presented on Table 1. Statistical significance was detected using the nonparametric Mann-Whitney U test.

PCR assays and sequencing to determine the rates of genomic deletion, translocation, and mutation

At least 200 5-FOAR colonies per strain and condition were subjected to PCR analysis with the following primers: AAAACAAACGCAAACAAGGC and GTTTAACTATGCTTCGTCGG to amplify ura4 ORF, TGAATCCTCCGTTCAGTAGG and AAGGACTGCGTTCTTCTAGC to amplify rng3 and TTTCCTTTCACGGCTAACCC (TLII) and TGTACCCATGAGCAAACTGC (TLIII) to amplify the translocation junction. The amplified ura4 fragments were then sequenced on both strands, with the primers used to amplify the ura4 ORF. Only mutations present on both strands were considered to determine mutation spectra. Deletions, mutations and translocations were scored as percentages of all events and these values were used to balance the rates of ura4 loss to determine the respective rates of deletion, mutation and translocation (see Figure S1 for deletion and translocation rates and Figure S2 for mutation rates).

The fork-arrest-induced deletion, translocation and mutation rates (Figure 2C–2E and Figure 3A) were calculated by subtracting the rate obtained in presence of thiamine (Rtf1 being repressed, OFF) from the rate obtained in the absence of thiamine (Rtf1 being expressed, ON). This method allows the spontaneous instability of IRs and the leakiness of the RTS1-RFB activity to be disregarded to determine strictly the rate of events induced by fork-arrest. The nonparametric Mann-Whitney U test was used to test for statistically significant differences.

Reverse Mutation Assay

Exponentially growing cells were treated with 20 mM of HU, 40 µM of CPT or 1 mM of MMC. At indicated times, samples were taken and appropriate dilutions were plated on supplemented minimal media to determine plating efficiency and on uracil-free plates. Colonies were counted after incubation at 30°C for 5–7 days and the frequency of Ura+ colonies was determined.

Replication Slippage Assay Using Ura4-Dup20 And Ura4-Dup22 Strains

For strains showing a slow growth phenotype (recombination mutants), the frequency of Ura+ revertants was determined as a function of the number of generations experiencing fork arrest at ura4. Exponentially growing 5-FOAR cells were washed twice in water and used to inoculate uracil-containing media without thiamine. Every 24 hours, cells were counted to determine the number of generations, and appropriate dilutions were plated on supplemented minimal media and on uracil-free plates. Colonies were counted after incubation at 30°C for 5–7 days and the frequency of Ura+ colonies was determined. The slope of the curves presented on Figure 4F corresponds to the rate of replication slippage/generation. For strains showing similar growth to wild-type cells, a single 5-FOAR colony was grown on uracil-containing plates with or without thiamine for 2–3 days, and then grown in uracil-containing media with or without thiamine for 2 days at 30°C. Appropriate dilutions were plated on supplemented minimal media and on uracil-free plates. Colonies were counted after incubation at 30°C for 5–7 days and the frequency of Ura+ colonies was determined.

Fork-arrest results in GCRs in a recombination-dependent manner. A. The rate of deletion for indicated strains, in the presence (Rtf1 repressed) and in the absence (Rtf1 expressed) of thiamine. Numbers of RTS1 repeats present in the S. pombe genome and the presence of a visible fork arrest (based on 2DGE presented on Figure 1B) are given for each strain. The % of deletion events, as determined by the PCR assay, was used to balance the rate of ura4 loss. The values reported are means of at least 3 independent median rates. Error bars correspond to the standard error (SE). Statistically significant fold differences in the rate of deletion events between the Rtf1 “repressed” and “expressed” conditions are indicated with an *. B and C. Rate of deletion (B) and translocation (C) for the strains indicated; ON and OFF refers to the RTS1-RFB being active or not, respectively. The % of deletion and translocation events, as determined by the PCR assay, was used to balance the rate of ura4 loss. The values reported are means at least 3 independent median rates. Error bars correspond to SE. Statistically significant fold differences in the rates of deletion or translocation events between the “OFF” and “ON” conditions are indicated with an *. Translocation events (based on the detection of the TLII/TLIII PCR product) were not detected in rad22-d or rhp51-d strains, whatever the conditional fork arrest construct.

(TIF)

Click here for additional data file.

Fork-arrest induces replication slippage. A. The rate of mutation for indicated strains, in the presence (Rtf1 repressed) and in the absence (Rtf1 expressed) of thiamine. Numbers of RTS1 repeats present in the S. pombe genome and the presence of a visible fork arrest (based on 2DGE presented on Figure 1) are given for each strain. The % of mutation events, as determined by the PCR assay and sequencing, was used to balance the rate of ura4 loss. The reported values are means of at least 3 independent median rates. Error bars correspond to SE. Statistically significant fold differences in the rate of mutation events between the Rtf1 “repressed” and “expressed” conditions are indicated with an *. B. The frequency of Ura+ revertants for the indicated strains and conditions. All strains harbour a non-functional ura4 allele due to a single base-substitution or a frame-shift or a duplication of 20 or 22 nt, together with the RTS1-RFB in the t-ura4

(TIF)

Click here for additional data file.

Features of replication slippage induced by fork arrest. A. Table of deletion/duplication and micro-homology features. B. Map of deletion and duplication events observed within the ura4 ORF in the t-ura4

(TIF)

Click here for additional data file.

Sensitivity of rad22-d strain to acute exposure to 20 mM of HU or 20 µM of CPT. The values reported are means of two to four independent experiments. Error bars indicate the standard error of the mean (SEM).

(TIF)

Click here for additional data file.

Fork-arrest-induced replication slippage is independent of the post-replication repair and mismatch repair. A–C. Left panels: Serial tenfold-dilutions of indicated strains cultured in thiamine-free medium spotted onto the medium indicated. RTS1-RFB activity “–” refers to the strain t-ura4-dup20-ori and “+” refers to the strain t-ura4-dup20Ura+ revertants from the strains indicated (t-ura4-dup20-ori associated or not with the RTS1-RFB) in the conditions indicated. The values reported are means of at least three independent experiments and error bars correspond to SEM.

(TIF)

Click here for additional data file.

Strains used in this study.

(DOCX)

Click here for additional data file.