Visualization of recombination–mediated damage-bypass by template switching

Template switching (TS) mediates damage-bypass via a recombination-related mechanism involving PCNA polyubiquitylation and Polymerase δ-dependent DNA synthesis. Using two-dimensional gel electrophoresis and electron microscopy, here we characterize TS intermediates arising in Saccharomyces cerevisiae at a defined chromosome locus, identifying five major families of intermediates. Single-stranded DNA gaps in the range of 150-200 nucleotides, and not DNA ends, initiate TS by strand invasion. This causes re-annealing of the parental strands and exposure of the non-damaged newly synthesized chromatid as template for replication by the other blocked nascent strand. Structures resembling double Holliday Junctions, postulated to be central double-strand break repair intermediates, but so far only visualized in meiosis, mediate late stages of TS, before being processed to hemicatenanes. Our results reveal the DNA transitions accounting for recombination-mediated DNA damage tolerance in mitotic cells and for replication under conditions of genotoxic stress.

The molecular mechanism of TS is only partially understood. Studies in yeast and mammalian cells provided evidence that discontinuities are created in the newly synthesized strands during replication in the presence of damage, opposite the lesions 5,6,14 . Repriming can reactivate transiently stalled forks in vitro 15 , potentially providing a molecular explanation for the formation of gaps. Pulse-labeling experiments indicated that the gaps are then filled-in using the newly replicated strand as template 6 . This notion was substantiated by various genetic and molecular studies that further identified the Rad5-pathway involving PCNA polyubiquitylation, as well as homologous recombination activities to be required in this process 4,5,13,16 . Finally, X-shaped structures involving sister chromatid junctions and having the genetic dependencies of TS intermediates were identified using 2D gel electrophoresis during replication of damaged templates 7 . Using this readout for TS, it was shown that a replication step mediated specifically by Polymerase δ is required 17 , in congruence with other observations 6,18,19 .
The RecQ helicase Sgs1 (BLM in mammalian cells), mutated in the cancer-prone Bloom syndrome, participates in error-free DDT by dissolving, in the context of the Sgs1-Top3-Rmi1 (STR) complex, the TS intermediates after their formation 7,19,20 . Persistent TS structures can also be resolved by the Mus81-Mms4 nuclease or by overexpressing other bacterial or human resolvases [21][22][23] . However, as all these enzymes can process a variety of substrates in vitro 24,25 , the structural nature of TS intermediates has remained undefined. Two main models of TS -strand invasion via a recombination-primed mechanism, or annealing-mediated fork regression -propose as intermediates DNA structures in which the sister chromatids are paired 2,26 . In vitro biochemical characterization of known TS factors provided evidence for either of these mechanisms 27,28,29 , and single, double, nicked Holliday Junctions (HJs), as well as hemicatenane-like structures were proposed as potential intermediates 12,17,20,23,30 .
Here we set out to identify the DNA structures and transitions involved in mediating errorfree DDT by TS using a combination of 2D gel electrophoresis and electron microscopy (EM) studies. This question is biologically important for several reasons. First, TS is important for genome stability by preventing incorporation of mutations 1,3,8,31 , which ultimately lead to genome instability and cancer 32 . Indeed, various TS factors function as tumor suppressors 1,33 . Second, DNA damaging agents, such as ultraviolet light, tobacco smoke and chemotherapeutics, which are potent inducers or occasionally target DDT pathways 8,19,34,35 , cause specific mutational signatures associated with genome instability 32 . Third, DNA gaps have been suggested as the predominant driver of recombination in mitotic cells, as opposed to double-strand breaks (DSBs) in meiosis 31,[36][37][38] , but the molecular mechanism through which ssDNA gaps are engaged in recombination remains elusive.

A linear minichromosome system to study TS-mediated DDT
To obtain a homogenous population of TS intermediates for EM studies, we employed a 2μbased 10.9Kb high-copy number minichromosome, YLpFAT7.1 ( Supplementary Fig. 1a) 39 . In Saccharomyces cerevisiae, YLpFAT7.1 minichromosomes are maintained as linear fragments, and are not entangled in concatemers ( Supplementary Fig. 1b). To improve the specificity towards the intermediates of interest, we designed a restriction enzyme mix that reduces the average length of the endogenous genomic DNA to less than 5Kb ( Supplementary Fig. 1c). We synchronized minichromosome-containing wt and sgs1Δ cells and then allowed them to replicate in the absence or presence of the DNA alkylating agent methyl-methanesulfonate, MMS (Fig. 1a). Similarly to endogenous chromosome replication, besides usual intermediates such as replication forks and bubbles, minichromosome replication in the presence of MMS induced the formation of X-shaped DNA structures, which further accumulated in sgs1Δ cells, proximal to the 2μ replication origin ( Fig. 1a and Supplementary Fig. 1a). The migration pattern of these joint molecules (JMs) appeared as two prominent spots on the X-arc when we performed psoralen-crosslinking in vivo prior to genomic extraction, instead of the smooth X-arc signal observed in the absence of psoralencrosslinking (Fig. 1b). TS intermediates arising on endogenous chromosomes showed the same migration pattern (Fig. 1b). These migration patterns are likely due to the fact that psoralen-crosslinking prevents migration of the junction point along the restriction fragment, stabilizing the initial X-shaped intermediates into specific spots 40 . Importantly, we found that the minichromosome-derived JMs shared the same genetic requirements with classical TS intermediates 7,12 , that is, their formation depended on the Rad51 and Rad5 pathways, in both wt and sgs1Δ backgrounds (Fig. 1c).

Isolation and assignment of TS intermediates
To visualize TS intermediates by EM, we psoralen-crosslinked in vivo genomic DNA of wt and sgs1Δ cells 40 , enriched for minichromosome-derived DNA by cleaving the endogenous chromosomes ( Supplementary Fig. 1c), and further digested with HindIII before running preparative 2D gels (Fig. 2a). We then excised the X-shaped signal from the ethidium bromide-stained second dimension gels (Fig. 2a), purified the DNA by electroelution, and confirmed via a subsequent 2D gel the migration pattern as branched molecules (Fig. 2b).
We proceeded with EM analysis of DNA samples obtained from two independent biological replicates. As we observed similar trends for all the crucial parameters described hereafter, in the following, we present the results obtained by pooling the number of JMs analyzed. We identified over 300 JMs for each wt and sgs1Δ cells ( Supplementary Fig. 1d). To ensure strict analysis of minichromosome-derived structures, we measured the length for all the initially identified JMs (649). The assignment procedure, which allowed a 10% tolerance range in measurements, considered that minichromosome-derived TS intermediates should represent two HindIII 5.8Kb inter-connected fragments, and the branches, taken two to two, should be the size of one HindIII fragment 41 . In this way, we unequivocally assigned 43% (wt) and 51% (sgs1Δ) out of the initially identified X-molecules as minichromosomederived intermediates ( Supplementary Fig. 1d). The JMs excluded by this assignment procedure are likely broken minichromosome-derived DNA or endogenous genomic DNA contaminants, caused by partial digestion. We further classified the assigned minichromosome-derived JMs into symmetric molecules-in which the junction point was stabilized by homologous pairing and two to two or all the four branches had the same length-, or asymmetric, with junction points not stabilized by homologous pairing and all four branches of different length or only two equal arms.

Structural and genetic features of symmetric JMs
Based on the structural features of the identified symmetric JMs, which accounted for 84% of the JMs, we classified these intermediates in 5 prominent families, F1 to F5 (Table 1). Besides the F1-F5 families, a small group of single Holliday Junction (SHJ)-like structures (8/306) was observed ( Table 2). We also classified the asymmetric JMs in two families, F1* and F5*. Table 2 summarizes the number and percentage of JMs belonging to each family.
We noted that two of the identified symmetric families were characterized by the presence of a ssDNA region proximal to the junction point, on one or two of the four branches comprising the JM. Because of the wealth of genetic and molecular information associating TS with gap-filling, we reasoned that the early intermediates are likely to contain such discontinuities close to the junction point. We called these families F1 and F2 ( Table 1). The junction point in F1 resembled a cross between two linear dsDNA filaments, and there was a ssDNA stretch on one branch immediately proximal to the junction point ( Fig. 3a and Supplementary Fig. 2), generally of 150-200 nt in length ( Supplementary Fig. 1e). This family of molecules was equally distributed in both wt and sgs1Δ cells ( Table 2 and Fig. 3c). The junction point in the F2 family involved a longer homology region and contained a thick DNA filament (Table 1 and Fig. 3b). Notably, F2 molecules contained a stretch of ssDNA on one or often two branches of the JM close to the junction point ( Fig. 3b and Supplementary Fig. 3) and were four times more prominent in sgs1Δ than in wt (Fig. 3c). We examined how F1 and F2 families were related by measuring the distribution of ssDNA length and other parameters (Fig. 3d). The median value of the ssDNA regions in F2 molecules (considering all molecules identified in both wt and sgs1Δ) was 104 nt, shorter than the median of the ssDNA region of F1 molecules (196 nt) (Fig. 3d). Notably, the median value of the sum of the ssDNA regions and the thick DNA filament present at the junction point in F2 molecules (195 nt) was very similar, in the 10% tolerance range, to the median value of ssDNA regions found in F1 molecules (196 nt) (Fig. 3d). When instead of the median values, we employed the average values for the parameters above the same pattern emerged (Fig. 3d). Moreover, the same was true when we performed the above calculations only for the F1 and F2 molecules found in sgs1Δ (the number of F2 molecules in wt was too low compared to the F1 molecules in the same background to make such analysis meaningful, Fig. 3c) (data not shown). These results indicate that the F2 molecules derive from F1, by partial annealing of the gap region into the homologous duplex, leading to a double Y-like DNA structure (Fig. 3e).
The other symmetric families we identified did not contain ssDNA discontinuities on the branches ( Table 1). The F3 family structurally resembled double HJs (dHJs). The junction point was similar to a bubble and contained two dsDNA-like filaments of roughly the same length, connecting the two HJ-like crosses, occasionally seen as rhomboid structures with an open center at the crossover position ( Fig. 4a and Supplementary Fig. 4). The distance between the two HJ-like centers generally ranged from 0.15 to 0.6 kb, but occasionally was longer (Fig. 4e). We also found molecules having both F2 and F3 features, called here F3*. The junction point in F3* still resembled a bubble, but we could always identify a ssDNA discontinuity on one of the two dsDNA filaments at the junction point ( Supplementary Fig.  5), likely indicating ongoing DNA synthesis on one of the strands. Based on the length of the duplex filaments at the F3 junction (Fig. 4e) and the ssDNA length in F1 molecules (Fig.  3d), we inferred that the DNA synthesis step implicated in TS 6,17,19 often proceeds beyond the initial size of the gap, concomitantly leading to more extensive annealing of the parental strands (as two DNA duplexes are present between the crossover junctions in F3 and one duplex contains the newly replicated strands, it follows that the other duplex must contain the paired parental strands, a notion for which we provide further evidence below). F4 resembles double Y or dHJs 42 (Table 1), and the junction point contains a small, almost closed bubble ( Fig. 4b and Supplementary Fig. 6). A DNA filament, generally around 80 nt, of a thickness greater than the one of a typical dsDNA filament and likely comprised by two dsDNA filaments, was present at the junction point ( Figs. 4b and 4e).
The F5 molecules resemble a cross between two dsDNA filaments ( Table 1, Fig. 4c and Supplementary Fig. 7). Occasionally, we detected at the junction point (and not on the branches as in F1) ssDNA filaments that appeared to connect the two DNA duplexes implicated in the junction ( Supplementary Fig. 7c).
Besides the F1-F5 families, we observed a few symmetric JMs with structural features of single HJs (SHJs) (8/306, Table 2 and Supplementary Fig. 8a). We suspect that they represent the products of recombination events triggered by occasional nicks arising during replication or induced by MMS lesion processing. From a structural standpoint, it is possible that they represent reversed forks, in which the length of the reversed arm is equal to the non-replicated strand. However, besides the low likelihood for the above, we note that six out of the eight identified SHJ-like molecules were in wt, and previous EM studies of total intermediates did not identify reversed forks under these conditions 14 .
As the symmetric families structurally resembled HJs, we set out to investigate their biochemical properties. Incubation at high temperature causes JMs to branch-migrate and be resolved into linear DNA when the branch points reach the ends [42][43][44] . However, branchmigration of canonical HJs, such as synthetic HJs and meiotic recombination intermediates containing HJs, is prevented by addition of magnesium cations 45 , which stabilize HJs in a stacked-X configuration 45,46 . We found that the minichromosome-derived TS intermediates arising in both wt and sgs1Δ cells showed comparable branch-migration activity regardless of the presence of magnesium cations in the buffer (Fig. 5a). Based on this feature, previously reported also for TS intermediates formed on endogenous chromosomes 20,30 , we conclude that the large fraction of TS intermediates are not canonical HJs.
Notably, we observed the same families of DNA structures, both in what regards the symmetric and the asymmetric junctions (see below) in sgs1Δ and wt cells ( Table 2). As sgs1Δ samples were specifically enriched in F2-F3-F4 families ( Table 2, Figs. 3c, 4d, and Supplementary Fig. 1e), we conclude that these types of DNA structures are substrates for Sgs1 in vivo.

Structural and genetic features of asymmetric JMs
Besides the symmetric molecules, 16% of junctions were asymmetric (49/306). We could divide them in two structural categories. The JMs of one of these categories, named F5* hereafter, had four branches of different length, no ssDNA discontinuities on the branches, and the junction points resembled either crosses between two linear DNA fragments or assumed more complex structures in which three ssDNA filaments forming a triangle-like structure were visible ( Fig. 5b and Supplementary Fig. 8b). As the junction point was not stabilized by homologous pairing, we could rule out F5* being SHJs, congruently with the results of the branch-migration assays (Fig. 5a). Structurally, F5* molecules resembled hemicatenanes and the symmetric F5 molecules. The relative percentage of F5 and F5* versus other JMs was increased in wt versus sgs1Δ (Table 2 and Supplementary Fig. 1e), in line with the notion that Sgs1 is not required for the last step of hemicatenane resolution 47 . Moreover, the abundance of F5 and F5* structures in wt in relation to other JMs (Table 2 and Supplementary Fig. 1e) suggests that hemicatenane-processing may not be very efficient in vivo. The time window from the formation of symmetric F5 hemicatenanes to their final processing may give opportunity for slippage of the junction point, leading to the formation of asymmetric hemicatenane-like F5* molecules, wherein two dsDNA filaments are connected by non-homologous pairing ( Fig. 5b and Supplementary Fig. 8b). We labeled this family F5* to illustrate their possible origin from F5 molecules. Substantiating this view, we observed the same nearly 50% decrease in number and relative percentage of F5* molecules in sgs1Δ compared to wt, as we did for F5 (Table 2 and Supplementary Fig. 1e).
We note that a subset of the F5* molecules (6/23 and 5/11 in wt and sgs1Δ, respectively, 11/306 of the total JMs) had two branches of roughly the same length, and therefore might represent reversed forks. The relative scarcity of JMs with structural properties compatible with those of reversed forks (less than 4% of the JMs) found here and their absence in previous EM studies 14 , largely support the postreplicative gap-filling model of TS 6,7,17,19,35 .
The other category of identified asymmetric JMs structurally resembled both hemicatenanes and F1 molecules in that they contain a ssDNA discontinuity, of about 150-200 nt, on one of the branches, proximal to the junction point ( Fig. 5c and Supplementary Fig. 8c). We labeled these intermediates F1* to illustrate the possibility, inferred from their structural features, of them being F1 precursors. We envisage that replication-associated hemicatenanes 44 may become stabilized into symmetric junctions (F1) when encountering a ssDNA gap 48,49 , thus facilitating strand invasion of the gap-region into the homologous sister duplex 48 .

Pseudo-Double Holliday Junctions mediate TS
dHJ-like structures of the F3-F4 families could arise during gap-filling if besides the plectonemic pairing of the newly synthesized DNA filaments, the parental strands also anneal in vivo (paranemic pairing, Fig. 6b, left panel). To analyze if each of the dsDNA filaments present at the junction point in F3 molecules contain two ssDNA filaments paired in vivo (as opposed to them being together because of occasional nicking induced by the spreading procedure), we analyzed similarly isolated JMs under denaturing conditions. In this assay, large multi-protein complexes -mostly nucleosomal histones -tightly bound to DNA molecules protect the ssDNA filaments of DNA duplexes from being psoralencrosslinked and, under conditions of denaturing spreading, their binding sites appear as strings of ssDNA bubbles, connected by duplex regions representing the unbound linker dsDNA that was psoralen-crosslinked 50 . As a proof of principle, we identified several molecules with F3/dHJ features, whose examination revealed that the ssDNA strands composing each of the two duplexes present at the junction point in F3 molecules (see Fig.  4a) are crosslinked and thus paired in vivo (Fig. 6a).
This type of structure is compatible with F3 molecules being either canonical double HJs, in which the two dsDNA filaments connecting the two HJ crosses are stabilized by plectonemic pairing (Fig. 6b right panel, Double Holliday Junction) or with a structure in which the parental ssDNA strands are connected to each other only by paranemic pairing, for instance by the action of Rad51 (Fig. 6b, left panel, Pseudo-Double Holliday Junction). We reasoned that the two types of structures could be distinguished if the chromosomal DNA was prepared without psoralen-crosslinking in vivo: dHJs would structurally remain the same, whereas the paranemic pairing in pseudo-dHJ structures would be dissolved exposing the two parental ssDNA filaments, which upon Mung Bean nuclease treatment should be digested and yield a double Y-like structure (Fig. 6b, left panel). Consistent with this latter prediction and resembling the behavior of TS intermediates arising on endogenous chromosomes 20,30 , Mung Bean nuclease treatment of minichromosome-enriched DNA prepared in the absence of in vivo psoralen-crosslinking led to a clear transition from the Xarc to a double Y structure (Fig. 6c, left panel). In contrast, the psoralen-crosslinked samples were much more resistant in this regard (Fig. 6c, right panel). Based on these biochemical and structural features, we propose that the F3 molecules accumulating in sgs1Δ and mediating TS are pseudo dHJs.

DISCUSSION
HJs are the most discussed hypothetical intermediates for recombination, modeled for DSB repair 51 , but to date only visualized in meiosis 42 . Here we provide evidence that noncanonical HJs (pseudo dHJs) are implicated in recombination-mediated gap-filling by TS. Analysis of the structural features of early intermediates (F1 and F2) revealed that a substantial part of the gap region becomes paired with the homologous duplex containing the newly synthesized filament (see evidence for paranemic pairing in Fig. 6), leading to a three-stranded DNA structure (F2, Figs. 3e and 7). The features of F2 molecules are consistent with a structure in which the gap region begins to be aligned with the homologous region of the sister duplex, but DNA synthesis did not begin (Figs. 3b, 3e and 7). Reannealing of the parental strands in F2 would facilitate the displacement/exposure of the non-damaged newly synthesized filament, which would then become available as a DNA synthesis template for the stalled 3′end proximal to the gap region (Figs. 3b and 7). Subsequent extension of the 3′end would give rise to F3*-like molecules, in which DNA synthesis is ongoing only on one filament of the F3 junction point ( Supplementary Fig. 5 and Fig. 7). As annotated in the features of the observed symmetric structures, measurements of the ssDNA regions and of the thick DNA filaments present at the junction point in F1-F2 molecules fully support this view (Fig. 3).
Alternative to the model and interpretation above, is that the F2 molecules already illustrate DNA synthesis-mediated gap-filling initiated by the invasion of the transiently stalled 3′ end into the homologous sister duplex and the formation of a D-loop-like structure. A caveat with this interpretation is that parental strand annealing is not necessarily expected to happen in this scenario, but it must be invoked to explain the subsequent paired parental strands observed in F3-F4 molecules (see Fig. 6a). Such pairing could happen if elongation of the D-loop-which this model envisages to happen already in F2 molecules, is coupled with annealing of the parental strands. However, such an interpretation is incongruent with the presence of ssDNA filaments on both sides of the junction point in F2 molecules (Fig. 3b). Thus, the data support the gap-mediated model of strand invasion, which is also consistent with other previous genetic deductions 17,19,37,48,52 .
Regardless, however, of whether TS is initiated by the gap region or by the adjacent 3′end close to the gap, completion of TS DNA synthesis can be followed by second-end capture, with formation of dHJ-like structures, or by strand displacement that would lead to linear molecules and non-crossover products. In agreement with the inference that F3 molecules represent dHJ-like molecules that are resolved by the STR complex 20,47,53,54 , we found an increased percentage of F3 molecules in sgs1Δ (15.5%) as compared to wt (4.6%) cells ( Table 2). Besides F3-F4 types, F2 and F3* families also accumulated in sgs1Δ samples ( Table 2), revealing that Sgs1 also potentiates D-loop displacement (Fig. 7), similar to what was proposed for the mammalian RecQ helicases RecQL5 and BLM 55,56 . We thus propose that both dHJ dissolution and strand displacement account for Sgs1 roles in suppressing sister chromatid exchanges and mitotic crossover 54,57 and in facilitating TS intermediate resolution 7,22 .
Dissolution of dHJ-like intermediates by STR is expected to ultimately lead to symmetric hemicatenanes or F5-like molecules (Fig. 7), later resolved by Top3. In addition to F5 structures, we detected intermediates of the F3 dissolution process (F4, Figs. 4b and 7), which still accumulated in sgs1Δ cells (Figs. 2c and 4d). This indicates that the last step of dissolution is slow, likely because of topological constrains, and that Sgs1 facilitates Top3 activity also in this context 58 (Fig. 7).
Previous work on factors that influence TS intermediate resolution failed to reveal a role for factors such as Srs2 that are important for D-loop unwinding 25,54,57 . Why problems in this early step do not translate into X-molecule accumulation is not well understood. Based on these new results, it seems likely that a combination of events-redundancy of D-loop unwinding activities, very efficient pathways of second end capture and subsequent dissolution of dHJ-like intermediates, as well as sensitivity limits of previously conducted 2D gel assays, is the underlying reason. In this new light, it also appears likely that the TSintermediate accumulation in smc5-smc6-mms21 mutants 59,60 is an indication that Smc5-Smc6-Mms21-mediated SUMOylation influences dHJ dissolution, possibly by facilitating STR activity or by enabling the formation of alternate Top3 complexes with similar in vivo activity with STR 59,60 . The method used here and the advances in the structural studies of replication and recombination intermediates promise to markedly contribute at unraveling how different genome integrity pathways-which show similar phenotypic outcomes when individually ablated-are integrated to support basic chromosome functions such as DDT and replication completion.

Cell culture, psoralen crosslinking and genomic DNA extraction for 2D gels
The yeast strains carrying the minichromosome YLpFAT7.1 39 were grown in synthetic media without leucine, arrested in G1 using alpha factor (SIGMA, 4 μg/ml) and released in YPD media or YPD media containing 0.033% MMS. Cells were fixed in 0.1% sodium azide for 40 minutes on ice as described 7,44,49 . The in vivo psoralen-crosslinking, the genomic DNA preparations (CTAB method), and the neutral-neutral 2D gels were performed as previously described 7,44 . In brief, for the in vivo psoralen-crosslinking cellular pellets corresponding to 1-2 × 10 9 cells were washed and resuspended in 5 ml of sterile cold water. Four cycles of 10 minutes each of in vivo psoralen-crosslinking were performed by placing the cell suspensions (transferred in 6 well-plastic tissue culture plates with flat bottom) in a Stratagene Stratalinker 1800 irradiation chamber, under 365 nm (UV). Before each irradiation cycle a fresh 0.3 ml aliquot of psoralen (4, 5′, 8 trimethyl-psoralen (TMP) 0.2 mg/ml dissolved in ethanol) was added to the cell suspension, mixed with cut blue tips and incubated for 5 minutes under aluminum foil. During the entire procedure, the plates containing the cell suspensions were placed on ice. Following this procedure, the cell suspensions were transferred to 50 ml Falcon tubes, the pellets were washed twice with sterile cold water and genomic DNA was extracted as previously described 43 .

Preparative 2D gels, electroelution and EM
For preparative 2D gels, 50 to 100 μg of genomic DNA was digested with a restriction enzyme mix containing PacI, XhoI, SgrAI, BssHII, SpeI, BclI, PmlI and HindIII, precipitated with potassium acetate /isopropanol and suspended in TE 1X or Tris-HCl 10mM pH 8. Digested genomic DNA was run on Owl Separation Systems, Inc Portsmouth, NH, USA, Model A2 electrophoresis chambers (gel tray 27×20cm-2.5L buffer) using TBE 1X not autoclaved and agarose gels of 500 ml. The first dimensions (0.35% agarose gels) were run at 50V for 24 hours at room temperature, and the second dimensions (0.9% agarose gels) for 12 hours at 150V (current limited at 150mA) at 4°C. For electroelution, agarose fragments were run in TBE 1X for 5 hours at 100V at room temperature in an Elutrap chamber (BioRad) with BT-1 and BT-2 membranes following the manufacturer's instructions. Recovered DNA solutions (around 0.8 ml) were concentrated on conical Amicon Ultra centrifugal filters (0.5ml 100K-MWCO 100K) following manufacturer's instructions. DNA samples from the electroelution were subjected to DNA spreading on carbon coated metal grids (4nm thickness) in the presence of uranyl acetate followed by platinum-based rotatory shadowing (0.4 nm without rotation and up to 8nm with rotation) and subjected to EM analysis as described 41 .

Branch migration and Mung Bean nuclease digestion assays
Branch migration and Mung Bean nuclease digestion assays were conducted as previously described 43 . For branch migration assays, the first dimension slices were incubated in TNE buffer (50mM NaCl, 10mM Tris-HCl pH 8, 0.1mM EDTA pH 8) or TNM buffer (50mM NaCl, 10mM Tris-HCl pH 8, 10mM MgCl 2, 0.1mM EDTA pH 8), for 13 hours at 65°C. For Mung Bean nuclease assays, 12 μg of digested genomic DNA was precipitated with potassium acetate/isopropanol and resuspended in Mung Bean nuclease buffer 1X. 40 units of Mung Bean nuclease (New England Biolabs) were added to each sample followed by incubation for 75 minutes at 37°C. DNA was precipitated again, resuspended in Tris-HCl 10mM pH 8 and first and second dimension gels were run as described above.

Assignment procedure for ssDNA regions on the DNA molecules analyzed
The assignment criteria for a ssDNA region on the DNA molecules analyzed by the EM procedure used here was recently described 41 . We note that in order to assign a ssDNA region on a DNA filament it is necessary to identify two points on the DNA molecule that define the borders of the ssDNA region, in which the thickness of the DNA filament decreases close to one half. The thickness of the DNA filament between these two points can slightly vary due to the positioning of the filament during the spreading procedure, but in general it is clearly lower than the average value of the thickness of a dsDNA filament, which in these experimental conditions is distributed around 20 Angstroms. While the difference between double-stranded and single-stranded DNA is promptly detected, it is generally difficult to visually distinguish thickness differences in multi-strand (>2) DNA. There are two technical explanations behind this: the first one is the reduced relative difference in thickness, a factor 2 from single-to double stranded transition, while only a factor of 1.5 from double to three stranded DNA. The second reason relates to the DNA visualization procedure: it was also previously shown that the combination of uranyl acetate staining and platinum rotatory shadowing used in this work is instrumental to appreciate the difference in thickness between dsDNA filaments and ssDNA filaments 41 . This is because the Platinum rotary shadowing only stains double stranded DNA, while the Uranyl atoms positively stain both single-and double-stranded DNA. Thus, the observed difference in thickness of the stained molecules -largely determined by the amount of deposited heavy atoms -generally exceeds the difference in thickness of the biological material itself. For the same reasons, three-stranded DNA (inferred to be present at the junction point in F2 molecules) will not be covered by much more heavy atoms than double-stranded, as the latter is already abundantly recognized by both staining procedures, reducing the relative contribution to thickness of an additional DNA strand.

Length measurements
The length measurements were performed using a conversion factor expressed in graphic units per base pairs/nucleotides (automatically converted in nanometers/nucleotide according to the magnification value at which the picture was taken), obtained by the measurement in pixels of a plasmid molecule of known length used as internal standard, as detailed in 41 . We note that the length of the ssDNA stretches reported here, using the same conversion factors for both dsDNA and ssDNA, might have been underestimated by a maximum experimentally determined factor of 15%, calculated for extreme conditions of long and entirely single-stranded DNA molecules that show reduced anchorage to the carbon layer and have different stretching properties.

Assignment procedure for the identified joint molecules
The assignment procedure of JMs considers that the total size of the X-molecules (the sum of the branches and, when applicable, of the extended homology region) must fall in the interval 10584bp-12953bp, which corresponds to the size of two HindIII 5.88 Kb linear fragments interconnected (see Supplementary Fig. 1a and 41 ). Moreover, the branches of an X-molecule taken two to two, together with the extended homology region for molecules in the category F2, F3 and F4, must have a length that falls in the interval 5292bp-6468bp, which corresponds to the size of one HindIII fragment. TS intermediates formed on YLpFAT7.1 minichromosomes. (a) Genomic DNA samples isolated from wild-type (wt, CY11340) and sgs1Δ (CY11357) strains carrying YLpFAT7.1 minichromosomes and replicating in the presence of MMS 0.033% were analyzed by 2D gel electrophoresis using a minichromosome-specific probe, leu2d. (b) Genomic DNA samples from wt (CY11340), sgs1Δ (CY11357) cells carrying the minichromosome YLpFAT7.1, and the corresponding strains without minichromosomes, wt (CY12486) and sgs1Δ (CY13249), were analyzed by 2D gel electrophoresis using probes for the minichromosome (leu2d) or the early origin of replication on the endogenous chromosome III (ARS305). (c) The genetic dependency of minichromosome-derived X-molecules was analyzed in wt (CY11340), sgs1Δ (CY11357), rad51Δ (CY12388), rad5Δ (CY12777), sgs1Δ rad51Δ (CY12390), sgs1Δ rad5Δ (CY12775).    Branch migration assays and asymmetric families of JMs. (a) Branch migration assay. Genomic DNA from wt (CY11340) and sgs1Δ (CY11357) strains carrying the minichromosome YLpFAT7.1 was extracted, digested with ApaI and BamHI, and run in triplicate. The first dimension gel slices were incubated in buffer with or without magnesium at 65°C 43 , before second dimension and subsequent analysis with the leu2d probe. The experiments were conducted three times with qualitatively identical results. A schematic representation of the expected position of the spot caused by branch migration of Xmolecules to two linear fragments in 2D gels is shown, together with a schematic representation of the restriction fragment analyzed. (b-c) Representative EM pictures of F5* and F1* families of X-molecules. The analyzed molecules derive from two independent biological experiments. Total views of the indicated X-shaped DNA structures and enlarged views with a schematic representation of the junction point (ssDNA regions in red and dsDNA regions in black). Branch sizes until the junction point are reported in kilobases. Length values indicated with a black or blue code indicate the pairs of branches leading to the size of one minichromosome-derived HindIII fragment. The junction points of the molecules are also represented with a black and light grey code, to show the contributions of the two DNA duplexes to the junction. Scale bars are shown. Giannattasio Table 1 Families of symmetric template switch intermediates.

F1
-simple cross -ssDNA gap proximal to the junction point Nat Struct Mol Biol. Author manuscript; available in PMC 2014 October 08.