Juxtaposition of heterozygous and homozygous regions causes reciprocal crossover remodelling via interference during Arabidopsis meiosis

Abstract
eLife digest
Introduction
Results
Discussion
Materials and methods
References
Article and author information
Metrics

Abstract

During meiosis homologous chromosomes undergo crossover recombination. Sequence differences between homologs can locally inhibit crossovers. Despite this, nucleotide diversity and population-scaled recombination are positively correlated in eukaryote genomes. To investigate interactions between heterozygosity and recombination we crossed Arabidopsis lines carrying fluorescent crossover reporters to 32 diverse accessions and observed hybrids with significantly higher and lower crossovers than homozygotes. Using recombinant populations derived from these crosses we observed that heterozygous regions increase crossovers when juxtaposed with homozygous regions, which reciprocally decrease. Total crossovers measured by chiasmata were unchanged when heterozygosity was varied, consistent with homeostatic control. We tested the effects of heterozygosity in mutants where the balance of interfering and non-interfering crossover repair is altered. Crossover remodeling at homozygosity-heterozygosity junctions requires interference, and non-interfering repair is inefficient in heterozygous regions. As a consequence, heterozygous regions show stronger crossover interference. Our findings reveal how varying homolog polymorphism patterns can shape meiotic recombination.

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

eLife digest

The genomes of plants and animals consist of several long DNA molecules that are called chromosomes. Most organisms carry two copies of each chromosome: one inherited from each parent. This means that an individual has two copies of each gene. Some of these gene copies may be identical (known as ‘homozygous’), but other gene copies will have sequence differences (or be ‘heterozygous’).

The sex cells (eggs and sperm) that pass half of each parent's genes on to its offspring are made in a process called meiosis. Before the pairs of each chromosome are separated to make two new sex cells, sections of genetic material can be swapped between a chromosome-pair to produce chromosomes with unique combinations of genetic material.

The ‘crossover’ events that cause the genetic material to be swapped are less likely to happen in sections of chromosomes that contain heterozygous genes. However, in a whole population of organisms, the exchange of genetic material between pairs of chromosomes tends to be higher when there are more genetic differences present.

Here, Ziolkowski et al. sought to understand these two seemingly contradictory phenomena by studying crossover events during meiosis in a plant known as Arabidopsis. The plants were genetically modified to carry fluorescent proteins that mark when and where crossovers occur. Ziolkowski et al. cross-bred these plants with 32 other varieties of Arabidopsis. The experiments show that some of these ‘hybrid’ plants had higher numbers of crossover events than plants produced from two genetically identical parents, but other hybrid plants had lower numbers of crossovers.

Ziolkowski et al. found that crossovers are more common between heterozygous regions that are close to homozygous regions on the same chromosome. The boundaries between these identical and non-identical regions are important for determining where crossovers take place. The experiments also show that the heterozygous regions have higher levels of interference—where one crossover event prevents other crossover events from happening nearby on the chromosome. In future, using chromosomes with varying patterns of heterozygosity may shed light on how this interference works.

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

Introduction

Sexual reproduction via meiosis is highly conserved within eukaryotes and allows recombination of genetic variation within populations (Barton and Charlesworth, 1998). During meiosis homologous chromosomes pair and undergo crossover recombination, which together with independent chromosome segregation and gamete fusion increases genetic diversity between progeny (Barton and Charlesworth, 1998; Villeneuve and Hillers, 2001). Meiotic crossovers form via the repair of DNA double-strand breaks (DSBs) generated by the SPO11 endonuclease (Bergerat et al., 1997; Keeney et al., 1997). Nucleolytic resection of DSBs generates 3′ single-stranded DNA (ssDNA), which is bound by the RAD51 and DMC1 recombinases (Bishop et al., 1992; Shinohara et al., 1992). The resulting nucleoprotein filament then invades a homologous chromatid to form a heteroduplex intermediate (Hunter and Kleckner, 2001). The invading ssDNA 3′-ends undergo DNA synthesis using the homologous duplex as a template and after second-end capture forms double Holliday junctions (dHJs) (Szostak et al., 1983; Schwacha and Kleckner, 1995). The dHJs can then be resolved as crossovers, which are cytologically evident as chiasmata (Page and Hawley, 2003; Janssens et al., 2012). Chiasmata hold chromosomes together and ensure that homologous pairs segregate to opposite cell poles, so that gametes inherit a balanced chromosome number (Page and Hawley, 2003).

Crossover numbers are under tight control, with many eukaryote species experiencing 1–2 per chromosome, despite large variation in genome size (Villeneuve and Hillers, 2001; Smukowski and Noor, 2011; Henderson, 2012; Mercier et al., 2014). In Arabidopsis ∼200 DSBs form per meiosis and proceed to form strand invasion intermediates, of which ∼10 are repaired as crossovers, with the excess being repaired as non-crossovers, or via intersister recombination (Giraut et al., 2011; Ferdous et al., 2012; Lu et al., 2012; Sun et al., 2012; Yang et al., 2012; Drouaud et al., 2013; Wijnker et al., 2013; Qi et al., 2014). 80–85% of wild type crossovers are dependent on the ZMM pathway (MSH4, MSH5, MER3, HEI10, ZIP4, SHOC1, PTD) and show interference, that is, they are spaced further apart than expected at random (Copenhaver et al., 2002; Higgins et al., 2004, 2008a; Chen et al., 2005; Mercier et al., 2005; Chelysheva et al., 2007, 2010, 2012; Macaisne et al., 2008). The remaining minority of crossovers are non-interfering and require MUS81 (Berchowitz et al., 2007; Higgins et al., 2008b). However, as chiasmata are still observed in msh4 mus81 double mutants, additional crossover pathways must exist (Higgins et al., 2008b). The majority of interhomolog strand invasion intermediates are dissolved by the FANCM helicase, which acts with the MHF1 and MHF2 co-factors (Crismani et al., 2012; Knoll et al., 2012; Girard et al., 2014). Mutations in FANCM, MHF1 and MHF2 cause dramatic increases in non-interfering crossovers (Crismani et al., 2012; Knoll et al., 2012; Girard et al., 2014). It is presently unclear whether non-interfering crossovers occurring in fancm are generated by the same pathway as in wild type, as a direct test of MUS81 dependence is precluded by fancm mus81 lethality (Crismani et al., 2012; Knoll et al., 2012). Both crossovers and non-crossovers can be accompanied by gene conversion events, which in the case of non-crossovers form via the synthesis-dependent strand annealing pathway (Allers and Lichten, 2001; McMahill et al., 2007; Lu et al., 2012; Sun et al., 2012; Yang et al., 2012; Drouaud et al., 2013; Wijnker et al., 2013; Qi et al., 2014).

Meiotic recombination is sensitive to DNA polymorphism between homologous chromosomes, that is, heterozygosity. For example, insertion-deletion (indel) and single nucleotide polymorphisms (SNPs) suppress crossovers at the scale of hotspots (kb) in fungi, plants and mammals (Dooner, 1986; Borts and Haber, 1987; Jeffreys and Neumann, 2005; Baudat and de Massy, 2007; Cole et al., 2010). This is thought to occur due to heteroduplex base-pair mismatches inhibiting recombination, following interhomolog strand invasion. Large scale chromosome rearrangements, such as inversions or translocations, also suppress crossovers (Schwander et al., 2014; Thompson and Jiggins, 2014). Despite the inhibitory effects of polymorphism on crossovers, nucleotide diversity and population-scaled recombination estimates are positively correlated in many plant and animal genomes (Begun and Aquadro, 1992; Hellmann et al., 2003; Spencer et al., 2006; Gore et al., 2009; Paape et al., 2012; Cutter and Payseur, 2013). For example, linkage disequilibrium-based crossover estimates and sequence diversity (π) are positively correlated in Arabidopsis at varying physical scales (Figure 1A and Table 1) (Cao et al., 2011; Choi et al., 2013). Multiple processes contribute to these relationships. For example, positive or negative directional selection can reduce diversity at linked sites, with a greater effect in regions of low recombination, known as hitchhiking and background selection (Hill and Robertson, 1966; Hudson and Kaplan, 1995; Nordborg et al., 1996; Smith and Haigh, 2007; Cutter and Payseur, 2013; Campos et al., 2014). These phenomena will cause regions of low recombination under selection to have low diversity, consistent with data in Drosophila (Aguade et al., 1989; Begun and Aquadro, 1992; Wiehe and Stephan, 1993; Campos et al., 2014). Recombination may also be mutagenic and increase diversity, for example via mismatch repair enzymes showing a mutational bias for A:T > G:C transversions (Duret and Galtier, 2009; Webster and Hurst, 2012; Glémin et al., 2014).

Figure 1

Download asset Open asset

Testing for crossover modification by Arabidopsis natural variation.

(A) Historical crossover frequency (red, cM/Mb) and sequence diversity (π, blue) along the physical length of the *Arabidopsis thaliana* chromosomes (Mb) (Cao et al., 2011; Choi et al., 2013). Mean values are indicated by horizontal dotted lines and centromeres by vertical dotted lines. The fluorescent crossover intervals analysed are indicated by solid vertical lines and coloured triangles. (B) Map showing the geographical origin of the *Arabidopsis* accessions studied, indicated by red points. (C) Genetic diagram illustrating the experimental approach with a single chromosome shown for simplicity. Fluorescent crossover reporters (triangles) were generated in the Col background (black) and crossed to accessions of interest (red) to generate F₁ heterozygotes. Following meiosis the proportion of parental:crossover gametes from F₁ heterozygotes was analysed to measure genetic distance (cM) between the fluorescent protein encoding transgenes.

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

Table 1

Correlations between historical recombination and sequence diversity at varying physical scales

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

Scale (π)	Chr1	Chr2	Chr3	Chr4	Chr5
5 kb	0.521	0.301	0.545	0.575	0.541
10 kb	0.556	0.305	0.565	0.602	0.562
50 kb	0.657	0.381	0.579	0.692	0.619
100 kb	0.699	0.563	0.601	0.744	0.646
500 kb	0.741	0.528	0.615	0.841	0.653
1 Mb	0.639	0.504	0.683	0.846	0.624

Scale (θ)	Chr1	Chr2	Chr3	Chr4	Chr5
5 kb	0.537	0.298	0.557	0.585	0.553
10 kb	0.569	0.303	0.576	0.610	0.572
50 kb	0.662	0.382	0.592	0.699	0.623
100 kb	0.710	0.573	0.617	0.752	0.650
500 kb	0.754	0.534	0.635	0.844	0.655
1 Mb	0.647	0.504	0.697	0.849	0.635

Spearman's rank correlation between historical crossover frequency estimates from LDhat and sequence diversity (θ and π) at varying physical scales (Cao et al., 2011; Choi et al., 2013). Adjacent windows of the indicated physical size were used for correlations.

Here we use natural variation in Arabidopsis to directly investigate the influence of heterozygosity on meiotic recombination. Extensive evidence exists for cis and trans modification of crossover frequency by plant genetic variation (Barth et al., 2001; Yao and Schnable, 2005; Yandeau-Nelson et al., 2006; Esch et al., 2007; McMullen et al., 2009; López et al., 2012; Salomé et al., 2012; Bauer et al., 2013). We define trans modifiers as loci encoding diffusible molecules that control recombination on other chromosomes, and elsewhere on the same chromosome, as exemplified by mammalian PRDM9 (Baudat et al., 2010; Berg et al., 2010; Myers et al., 2010; Parvanov et al., 2010; Fledel-Alon et al., 2011; Sandor et al., 2012; Kong et al., 2013). We define cis modification as variation that influences recombination only on the same chromosome, for example, the inhibitory effects of high SNP density, inversions and translocations (Dooner, 1986; Borts and Haber, 1987; Jeffreys and Neumann, 2005; Baudat and de Massy, 2007; Cole et al., 2010; Schwander et al., 2014; Thompson and Jiggins, 2014). Regional patterns of chromatin and epigenetic information can also cause significant cis effects, for example loss of either H2A.Z deposition or DNA methylation alters crossover frequency in Arabidopsis (Colomé-Tatché et al., 2012; Melamed-Bessudo and Levy, 2012; Mirouze et al., 2012; Yelina et al., 2012; Choi et al., 2013).

In this study we crossed Arabidopsis lines carrying fluorescent crossover reporters generated in a common background (Col-0) to 32 diverse accessions. We observed extensive variation in F₁ hybrid recombination rates, with both significantly higher and lower crossovers than homozygous backgrounds. We further analysed Col × Ct F₂ recombinant populations using three independent crossover reporter intervals (420, CEN3 and I2f). We did not detect trans modifiers in these crosses, but observed a novel cis modification effect caused by heterozygosity. Specifically, juxtaposition of heterozygous and homozygous regions is associated with increased crossover frequency in the heterozygous region and a reciprocal decrease in the homozygous region. To investigate this phenomenon mechanistically we repeated analysis in mutants where the balance of interfering and non-interfering crossover repair is altered (fancm, zip4 and fancm zip4). This analysis demonstrates that remodelling of crossovers across heterozygosity/homozygosity junctions is dependent on interference. We also show that the non-interfering repair is less efficient in heterozygous regions. As a consequence, interference measurements are stronger in heterozygous regions. Our findings show how varying polymorphism patterns can differentially influence meiotic recombination along chromosomes.

Results

Heterozygosity extensively modifies crossover frequency in Arabidopsis

To test the effect of heterozygosity on meiotic recombination we crossed transgenic Arabidopsis with fluorescent crossover reporters generated in the Col-0 background to 32 diverse accessions that represent global genetic diversity within this species (Figure 1, Tables 2, 3) (Melamed-Bessudo et al., 2005; Berchowitz and Copenhaver, 2008; Yelina et al., 2013). The 5 intervals tested (I1b, I1fg, I2f, 420 and CEN3) range from 0.67–5.40 megabases (Mb), represent 11.5% of the genome (14.34 Mb) in total and are located in sub-telomeric, interstitial and centromeric regions (Figure 1A and Table 2). The intervals vary in experimental recombination rate, with the centromeric interval CEN3 being the lowest (2.11 cM/Mb) and the sub-telomeric interval I2f being the highest (13.02 cM/Mb) (Table 2). As Arabidopsis male meiosis shows elevated sub-telomeric recombination, this likely contributes to the high I2f crossover frequency, which is measured in pollen (Giraut et al., 2011). Low recombination in CEN3 is also expected, as the centromeric regions are heterochromatic and known to show suppressed crossover frequency (Figure 1A) (Copenhaver et al., 1999; Giraut et al., 2011; Salomé et al., 2012; Yelina et al., 2012). To asses relative heterozygosity levels we analysed pairwise sequence differences relative to Col-0 using the 19 genomes dataset, which was generated from a subset of the accessions used in our crosses (Gan et al., 2011). CEN3 shows the highest heterozygosity levels, followed by the interstitial and sub-telomeric intervals (Table 2). Therefore, the regions analysed represent diverse chromosomal environments with varying levels of recombination and inter-accession sequence polymorphism.

Table 2

Fluorescent crossover reporter intervals

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

Interval	Chr	Method	T-DNA 1	T-DNA 2	Mb	Location	cM/Mb (Col-0)	cM/Mb (F₁)	Heterozygosity
I1b	1	Pollen	3,905,441-YFP	5,755,618-dsRed2	1.85	Interstitial	4.25	4.05	1.93 (3.16)
I1c	1	Pollen	5,755,618-dsRed2	9,850,022-CFP	4.09	Interstitial	4.55	N/A	2.80 (3.16)
I1fg	1	Pollen	24,645,163-YFP	25,956,590-dsRed2	1.31	Interstitial	6.20	6.02	2.52 (3.16)
I2a	2	Pollen	12,640,092-CFP	13,226,013-YFP	0.59	Interstitial	5.19	N/A	2.33 (3.30)
I2b	2	Pollen	13,226,013-YFP	14,675,407-dsRed2	1.45	Interstitial	3.09	N/A	1.53 (3.30)
I2f	2	Pollen	18,286,716-dsRed2	18,957,093-YFP	0.67	Sub-telomeric	13.02	17.41	1.43 (3.30)
420	3	Seed	256,516-GFP	5,361,637-dsRed2	5.11	Sub-telomeric	3.70	2.93	1.19 (3.37)
CEN3	3	Pollen	11,115,724-YFP	16,520,560-dsRed2	5.40	Centromeric	2.11	2.38	6.69 (3.37)
I3b	3	Pollen	498,916-CFP	3,126,994-YFP	2.63	Sub-telomeric	5.99	N/A	1.11 (3.37)
I3c	3	Pollen	3,126,994-YFP	4,319,513-dsRed2	1.19	Sub-telomeric	4.01	N/A	1.64 (3.37)
I5c	5	Pollen	2,372,623-CFP	3,760,756-YFP	1.39	Interstitial	4.01	N/A	1.01 (3.27)
I5d	5	Pollen	3,760,756-YFP	5,497,513-dsRed2	1.74	Interstitial	3.20	N/A	1.56 (3.27)

The interval name is listed together with chromosome, method of scoring and location of the flanking T-DNAs together with the fluorescent proteins they encode. Interval cM/Mb values from Col-0 homozygous are listed (Col-0), in addition to the mean cM/Mb observed across all F₁ crosses (F₁). Heterozygosity values were calculated using pairwise comparison of polymorphism data from the 19 genomes project to the Col reference (Gan et al., 2011), and the mean value for the interval shown, in addition to the mean chromosome value in parentheses.

Table 3

Genetic distance in F₁ heterozygotes

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

Accession	Location	I1b	I1fg	I2f	420	CEN3	Total	P
Tsu-0	Tsushima, Japan	6.6	6.3	6.9	14.5	9.4	43.7	<2.00 × 10⁻¹⁶
Hi-0	Hilversum, Netherlands	6.8	6.9	6.9	13.6	9.6	43.8	<2.00 × 10⁻¹⁶
Wil-2	Vilnius, Lithuania	6.1	6.9	6.1	15.9	10.1	45.0	<2.00 × 10⁻¹⁶
Kn-0	Kaunas, Lithuania	7.4	6.6	8.0	15.5	8.7	46.2	<2.00 × 10⁻¹⁶
Ler-0	Gorzow, Poland	6.6	8.2	7.6	12.3	11.9	46.6	<2.00 × 10⁻¹⁶
Ws-0	Vassilyevichy, Belarus	6.7	7.7	10.2	13.0	9.0	46.6	<2.00 × 10⁻¹⁶
No-0	Nossen, Germany	7.4	7.9	6.7	14.1	11.4	47.4	<2.00 × 10⁻¹⁶
Wu-0	Wurzburg, Germany	7.6	6.3	9.5	14.0	11.4	48.8	<2.00 × 10⁻¹⁶
Zu-0	Zurich, Switzerland	7.5	7.1	13.4	12.2	9.9	50.1	0.0438
Po-0	Poppelsdorf, Germany	7.2	7.9	9.1	15.8	10.9	51.0	0.000484
Ct-1	Catania, Italy	7.8	8.7	7.2	15.9	12.1	51.7	9.27 × 10⁻⁰⁸
Oy-0	Oystese, Norway	7.7	8.4	8.5	15.7	12.5	52.8	0.969
Rsch-4	Rschew, Russia	7.9	6.8	10.7	15.2	12.4	53.1	0.505
Col-0	Columbia, USA	8.0	8.2	8.8	18.0	11.5	54.5	–
Sf-2	San Feliu, Spain	8.2	8.8	7.4	18.6	12.3	55.3	0.724
Kas	Kashmir, India	6.9	8.6	13.2	13.8	13.3	55.8	<2.00 × 10⁻¹⁶
Kond	Pugus, Tajikistan	7.1	8.1	15.8	13.7	11.4	56.2	<2.00 × 10⁻¹⁶
Edi-0	Edinburgh, Scotland	8.0	8.0	13.4	13.3	13.6	56.3	<2.00 × 10⁻¹⁶
Bay-0	Bayreuth, Germany	8.6	8.3	11.3	18.6	11.5	58.3	<2.00 × 10⁻¹⁶
Mt-0	Martuba, Libya	9.6	7.8	13.2	20.6	9.6	60.8	<2.00 × 10⁻¹⁶
Sha	Pamiro-Alaya, Tajikistan	7.8	7.5	20.0	7.0	18.6	60.9	0.0012
C24	Columbia, USA	8.8	8.5	18.5	12.1	14.1	61.9	<2.00 × 10⁻¹⁶
Bur-0	Burren, Ireland	6.7	9.1	21.9	14.7	17.8	70.2	<2.00 × 10⁻¹⁶
Cvi-0	Cape Verde Islands	9.1	10.0	11.3	12.6	27.6	70.7	<2.00 × 10⁻¹⁶
Can-0	Las Palmas, Canary Isles	7.8	8.5	22.1	12.4	31.4	82.2	<2.00 × 10⁻¹⁶
Co	Coimbra, Portugal	–	–	–	11.1	13.8	–	–
Nw-0	Neuweilnau, Germany	–	–	–	14.7	14.4	–	–
Mh-0	Szczecin, Poland	–	–	–	14.9	10.1	–	–
Wl-0	Wildbad, Germany	–	–	–	17.0	9.5	–	–
Bu-0	Burghaun, Germany	–	–	–	28.9	8.8	–	–
CIBC5	Ascot, United Kingdom	–	–	–	13.2	11.3	–	–
RRS7	North Liberty, USA	–	–	–	17.2	11.7	–	–
	F₁ cM mean	7.6	7.9	11.5	15.0	12.9	54.8
	cM StDev	0.8	0.9	4.8	3.6	4.9	9.3

The accessions crossed to are listed with their geographic location. Genetic distance (cM) data is shown for the five fluorescent intervals, in addition to a summed total. Also shown are the mean and standard deviation for all F₁s. A generalized linear model (GLM) was used to test for significant differences between total recombinant vs non-recombinant counts between replicate groups of Col-0 homozygotes and F₁ heterozygotes. Tests were performed for genotypes where data from all five tested intervals had been collected.

The crossover reporter systems utilize fluorescent proteins encoded by linked, heterozygous transgenes that are expressed from the pollen-specific LAT52, or seed-specific NapA promoters (Melamed-Bessudo et al., 2005; Francis et al., 2006; Yelina et al., 2013). Fluorescent measurements of gametes or progeny are used to asses segregation of the transgenes through meiosis and thereby measure crossover rates (Melamed-Bessudo et al., 2005; Berchowitz and Copenhaver, 2008; Yelina et al., 2013). Previously, we developed flow cytometry protocols to increase scoring-throughput using fluorescent pollen, allowing up to 80,000 gametes to be scored per individual plant (Yelina et al., 2012, 2013). To increase throughput when measuring fluorescent seed we adapted CellProfiler image analysis software, allowing us to rapidly score ∼2000 seed per individual (Figure 2A–F) (Carpenter et al., 2006). This method gives recombination measurements not significantly different from manually collected data (Figure 2F, Figure 2—source data 1) (generalized linear model (GLM), hereafter GLM, p = 0.373). To test for significant differences between recombinant and non-recombinant counts using replicate groups we used a GLM assuming a binomial count distribution. Replicate heterozygous F₁ individuals were analysed for each cross and 13,264,943 gametes were scored in total, to provide an extensive survey of the influence of polymorphism heterozygosity on crossover frequency (Figure 3 and Table 3).

Figure 2 with 1 supplement see all

Download asset Open asset

High-throughput measurement of crossover frequency using image analysis of fluorescent seed.

(A) Combined red and green, red alone and green alone fluorescent micrographs of seed from a self-fertilized *420/++* plant. (B) CellProfiler output showing identification of seed objects by green lines and scoring of red and green fluorescence shown by shading. Blue shading shows an absence of colour. (C–D) Histograms of seed object fluorescence intensities, with coloured and non-coloured seed divided by vertical dotted lines. (E) Plot of seed object red vs green fluorescence intensities, with each point representing an individual seed. The red and green dashed lines show the colour vs non-colour divisions indicated in (C–D). The formula used for cM calculation is printed below. (F) *420* cM measurements from replicate plants of the indicated genotypes (Col/Col F₁, Col/Ler F₁, Col/Sha F₁) are shown by black dots with mean values indicated by red dots. Data generated by automatic and manual scoring are plotted alongside one another. Measurements made by the different methods were not significantly different as tested using generalized linear model (GLM). See Figure 2—source data 1.

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

Figure 2—source data 1 420 crossover frequency measured via manual or automated scoring of seed fluorescence.: https://doi.org/10.7554/eLife.03708.008
Download elife-03708-fig2-data1-v2.docx

Figure 3

Download asset Open asset

Variation in F₁ hybrid crossover frequency.

(A–E) Genetic distance (cM) measurements for fluorescent crossover intervals *I1b*, *I1fg*, *I2f*, *420* and *CEN3* with individual replicates (black dots) and mean values (red dots) for the crosses labelled on the x-axis. See Figure 3—source data 1–5. (F) Heatmap summarising crossover frequency data for F₁ crosses with data from all five intervals. Accessions are listed as rows and fluorescent intervals listed as columns. The heatmap is ordered according to ascending ‘Total’ cM (red = highest, blue = lowest), which is the sum of the individual interval genetic distances. GLM testing for significant differences between total recombinant vs non-recombinant counts between replicate groups of Col-0 homozygotes and F₁ heterozygotes was performed, for genotypes where data from all five tested intervals were collected (Table 3). Col/Col homozygous data are labelled and highlighted with an arrow in each plot.

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

Figure 3—source data 1 I1b F₁ flow cytometry count data.: https://doi.org/10.7554/eLife.03708.011
Download elife-03708-fig3-data1-v2.docx
Figure 3—source data 2 I1b F₁ flow cytometry count data.: https://doi.org/10.7554/eLife.03708.012
Download elife-03708-fig3-data2-v2.docx
Figure 3—source data 3 I1b F₁ flow cytometry count data.: https://doi.org/10.7554/eLife.03708.013
Download elife-03708-fig3-data3-v2.docx
Figure 3—source data 4 I1b F₁ flow cytometry count data.: https://doi.org/10.7554/eLife.03708.014
Download elife-03708-fig3-data4-v2.docx
Figure 3—source data 5 CEN3 F₁ flow cytometry count data.: https://doi.org/10.7554/eLife.03708.015
Download elife-03708-fig3-data5-v2.docx

We observed substantial variation in crossovers between F₁ crosses, although the interstitial intervals varied less than those in sub-telomeric and centromeric locations (Figure 3A–E, Figure 3—source data 1–5). F₁ heterozygotes showed both significantly higher and lower total recombination compared to Col homozygotes (Figure 3 and Table 3) (GLM with 113° of freedom p < 2.0 × 10⁻¹⁶). F₁ genetic distances and polymorphism levels within the intervals were poorly correlated, consistent with previous observations (Table 4) (Barth et al., 2001; Gan et al., 2011; Salomé et al., 2012). This weak correlation may be partially explained by unknown structural rearrangements. For example, the Shahdara (Sha) accession has a sub-telomeric inversion (3–5.1 Mb) on chromosome 3 relative to Col (Loudet et al., 2002; Simon et al., 2008; Salomé et al., 2012), and Col/Sha F₁s show consistently low crossovers in 420, which overlaps the inversion (Figure 3D and Table 3). Hence the contribution of unknown structural polymorphisms to variation in recombination rates could be significant. Further evidence of the complex effect of polymorphism is evident from the CEN3 interval, which spans the repetitive and structurally diverse centromeric region of chromosome 3 (Figure 1A) (Copenhaver et al., 1999; Clark et al., 2007; Ito et al., 2007; Cao et al., 2011; Gan et al., 2011; Horton et al., 2012; Long et al., 2013), and showed high variability in F₁ crossover frequency (Figure 3E and Table 3). Unexpectedly, some of the most diverged crosses, for example two accessions from Atlantic islands Cvi-0 and Can-0, showed highest CEN3 crossovers (Figure 3E and Table 3) (Ito et al., 2007). 10 of 26 F₁s showed significantly higher summed crossover frequency compared with Col homozygotes, consistent with previous reports that recombination can increase in heterozygous backgrounds in Arabidopsis (Barth et al., 2001) (Figure 3F and Table 3). Both cis and trans modification of crossovers by genetic variation has been observed in plants (Barth et al., 2001; Yao and Schnable, 2005; Yandeau-Nelson et al., 2006; Esch et al., 2007; McMullen et al., 2009; López et al., 2012; Salomé et al., 2012; Bauer et al., 2013). Therefore, the variation in F₁ crossover frequency observed here is likely caused by complex interactions between cis and trans modifying effects.

Table 4

F₁ heterozygosity levels relative to Col-0

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

Accession	Chr 1	I1b	I1fg	Chr 2	I2f	Chr 3	420	CEN3	Chr 4	Chr 5
Bur-0	3.35	1.86	3.62	3.60	1.51	3.58	1.57	6.20	3.89	3.16
Can-0	3.75	2.99	3.51	3.92	0.92	3.98	1.02	8.27	5.34	4.24
Ct-1	2.62	1.67	2.29	2.61	1.85	3.35	0.96	6.91	3.23	3.36
Edi-0	3.30	1.91	3.64	3.26	0.91	3.05	1.34	5.48	3.42	3.81
Hi-0	2.43	1.59	1.87	1.80	1.50	2.58	1.07	4.62	2.69	2.46
Kn-0	3.15	1.78	2.85	3.35	2.18	3.58	1.49	6.69	3.76	3.40
Ler-0	3.10	1.61	2.66	3.62	2.24	3.43	1.13	7.39	3.87	3.53
Mt-0	3.02	1.77	1.16	3.49	1.57	3.17	1.07	5.70	3.95	2.71
No-0	3.25	2.28	2.71	3.36	1.27	3.52	1.21	7.14	3.51	3.56
Oy-0	3.48	1.68	2.10	3.05	0.58	2.94	1.23	6.16	2.95	2.72
Po-0	2.45	1.78	1.19	2.36	0.67	2.87	0.79	5.99	2.53	2.59
Rsch-4	2.94	1.84	1.17	3.36	1.22	3.09	1.05	5.37	3.89	3.22
Sf-2	3.61	1.94	4.24	3.54	2.06	3.74	1.30	8.24	3.81	3.58
Tsu-0	3.37	1.68	2.36	3.69	1.39	3.98	1.14	8.78	3.69	3.05
Wil-2	3.56	1.99	2.45	3.77	2.11	3.81	1.56	7.55	4.44	3.34
Ws-0	3.25	1.93	3.54	3.68	1.58	3.30	1.30	6.65	3.70	3.41
Wu-0	3.13	2.53	1.95	3.14	0.67	3.50	1.22	7.41	3.36	3.15
Zu-0	3.10	1.85	2.02	3.83	1.43	3.19	0.96	5.84	3.38	3.64
Mean	3.16	1.93	2.52	3.30	1.43	3.37	1.19	6.69	3.63	3.27
Correlation (cM)	–	0.13 (p = 0.61)	0.47 (p = 0.05)	–	−0.29 (p = 0.23)	–	0.06 (p = 0.81)	0.28 (p = 0.26)	–	–

Accessions sequenced as part of the 19 genomes project were analysed (Gan et al., 2011) and heterozygosity calculated as the sum of SNPs and indel lengths divided by the length of region (kb). Correlations were between heterozygosity within the interval measured and F₁ cM measurements.

Modification of crossover frequency by juxtaposition of heterozygosity and homozygosity

To investigate the extent of cis and trans modification of crossover frequency by heterozygosity we generated a 420 Col × Ct recombinant F₂ population (n = 139) (Figure 4A). We selected F₂ individuals that were heterozygous for linked T-DNAs expressing red and green fluorescent proteins and Col/Ct heterozygous within 420, but genetically mosaic elsewhere in the genome (Figure 4A,E). The 420/++ Col/Ct F₂ population showed significantly greater variation in recombination rates than Col/Col homozygotes (F-test p = 0.0129) (Figure 4D, Figure 4—source data 1). We genotyped 51 Col/Ct markers throughout the genome and tested for their association with 420 crossover frequency using QTL analysis. We detected no association using markers on chromosomes 1, 2, 4 or 5 (Figure 4B). However, on chromosome 3 itself homozygosity (Col/Col or Ct/Ct) outside of 420 was associated with high recombination (FDR corrected chi-square test p = 3.29 × 10⁻³¹) (Figure 4B,E–F and Table 5). For each marker we used the heterozygous and homozygous counts in the hottest quartile vs the coldest quartile to construct 2 × 2 contingency tables and performed chi-square tests, followed by FDR correction for multiple testing (Table 5).

Figure 4 with 1 supplement see all

Download asset Open asset

Modification of crossover frequency by juxtaposition of heterozygosity and homozygosity.

(A) Diagram illustrating chromosome 3 genotypes (black = Col, red = Ct) in *RG/++* F₁ individuals and their F₂ progeny. A single chromosome is shown for simplicity. Gametes or progeny are analysed for patterns of fluorescence following meiosis to measure genetic distance. (B) The program Rqtl was used to test for association between Col/Ct genotypes and *420* cM in a *420/++* F₂ population. The logarithm of odds (LOD) score is plotted along the 5 chromosomes with the positions of markers shown along the x-axis by ticks. The red horizontal line shows the 5% genome-wide significance threshold calculated with Hayley-Knott regression and by running 10,000 permutations. (C) As for (B) but analyzing Col/Ct markers on chromosomes 2 and 3 for an *I2f/++* F₂ population. (D) *420* cM measurements from Col/Ct *420/++* F₂ (black), Col/Col homozygotes (red) and Col/Ct F₁ (blue) individuals. Mean values are indicated by horizontal dotted lines. See Figure 4—source data 1. (E) Chromosome 3 genotypes shown for *420/++* F₂ individuals ranked by crossover frequency. Each horizontal row represents a single F₂ individual. X-axis ticks show marker positions, and which are coloured red when they showed significantly higher homozygosity in the hottest vs coldest quartiles (FDR-corrected chi square test). Fluorescent T-DNAs are indicated by triangles, in addition to the centromere (Cen). (F) Heterozygosity along chromosome 3 in the hottest (red), coldest (blue) *420* F₂ quartiles and the mean (green). The locations of reporter T-DNAs and the centromeres are indicated by vertical dashed lines. (G–I) As for (D–F) but for interval *I2f*. See Figure 4—source data 2.

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

Figure 4—source data 1 420 Col/Ct F₂ fluorescent seed count data.: https://doi.org/10.7554/eLife.03708.018
Download elife-03708-fig4-data1-v2.docx
Figure 4—source data 2 I2f Col/Ct F₂ fluorescent seed count data.: https://doi.org/10.7554/eLife.03708.019
Download elife-03708-fig4-data2-v2.docx
Figure 4—source data 3 CEN3 Col/Ct F₂ flow cytometry count data.: https://doi.org/10.7554/eLife.03708.020
Download elife-03708-fig4-data3-v2.docx

Table 5

Chromosome 3 genotype counts from hot and cold quartile 420/++ Col/Ct F₂ individuals

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

Marker coordinates (bp)	Hot quartile HET	Hot quartile HOM	Cold quartile HET	Cold quartile HOM	FDR p value
259000	34	0	34	0	1
2718000	34	0	34	0	1
5352000	34	0	34	0	1
6375000	21	13	34	0	4.36 × 10⁻⁰⁴
6948000	17	17	33	1	1.05 × 10⁻⁰⁴
7674000	15	19	33	1	2.12 × 10⁻⁰⁵
8495000	12	22	34	0	3.65 × 10⁻⁰⁷
9404000	8	26	33	1	3.79 × 10⁻⁰⁸
10695000	8	26	30	4	1.36 × 10⁻⁰⁶
11649000	11	23	27	7	4.36 × 10⁻⁰⁴
12356000	11	23	27	7	4.36 × 10⁻⁰⁴
15949000	13	21	23	11	4.48 × 10⁻⁰²
19165000	17	17	21	13	0.591
23040000	13	21	17	17	0.591

The number of 420/++ Col/Ct F₂ individuals showing Col homozygosity (HOM) or Col/Ct heterozygosity (HET) for the indicated marker positions, in either the hottest or coldest F₂ quartile. The p value was obtained by performing a chi square test between homozygous and heterozygous marker genotype counts in the hottest and coldest quartiles (2x2 contingency table), followed by FDR correction for multiple testing.

To test an additional chromosome for the effect of heterozygosity/homozygosity juxtaposition we measured recombination in an I2f Col × Ct F₂ population (n = 78) (Figure 4G–I). The I2f interval is 0.67 Mb and located sub-telomerically on the long arm of chromosome 2 (Figure 1A and Table 2). The I2f/++ Col/Ct F₂ population also showed significantly greater variation in recombination rates than Col/Col homozygotes (F-test, p = 0.04) (Figure 4G, Figure 4—source data 2). We performed QTL analysis for Col/Ct markers on chromosomes 2 and 3 and again observed a significant effect on the same chromosome and no trans effect from chromosome 3. An identical trend to that seen for 420 was observed, with the highest recombination F₂ quartile showing significantly greater marker homozygosity (both Col/Col and Ct/Ct) outside I2f on chromosome 2 (FDR corrected chi-square test p = 1.44 × 10⁻¹⁰) (Figure 4C,G–I and Table 6). The most distal marker showing a significant difference between hot and cold quartiles was of comparable megabase distance for 420 (10.60 Mb) and I2f (10.12 Mb).

Table 6

Chromosome 2 genotype counts from hot and cold quartile I2f/++ Col/Ct F₂ individuals

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

Marker coordinates (bp)	Hot quartile HET	Hot quartile HOM	Cold quartile HET	Cold quartile HOM	FDR p value
132,000	9	11	8	12	1
2,346,000	7	13	8	12	1
4,748,000	8	12	9	11	1
6,789,000	7	13	11	9	0.63
11,443,000	5	15	20	0	6.26 × 10⁻⁰⁵
13,036,000	7	13	20	0	3.32 × 10⁻⁰⁴
14,117,000	9	11	20	0	1.30 × 10⁻⁰³
15,240,000	9	11	20	0	1.30 × 10⁻⁰³
16,909,000	13	7	20	0	0.0262
17,439,000	16	4	20	0	0.238
18,287,000	20	0	20	0	1
18,960,000	20	0	20	0	1
19,311,000	18	2	20	0	0.764

The number of I2f/++ Col/Ct F₂ individuals showing Col homozygosity (HOM) or Col/Ct heterozygosity (HET) for the indicated markers, in either the hottest or coldest F₂ quartile. The p value was obtained by performing a chi square test between homozygous and heterozygous marker genotype counts in the hottest and coldest quartiles (2 × 2 contingency table), followed by FDR correction for multiple testing.

To test whether the effect of heterozygosity/homozygosity juxtaposition is dependent on chromosomal location we measured crossovers in a CEN3 Col × Ct F₂ population (n = 121) (Figures 4A and 5C, Figure 4—figure supplement 1, Figure 4—source data 3). As for 420 and I2f, CEN3 F₂ recombination rates were significantly more variable than Col/Col homozygotes (F-test p = 0.01268) (Figure 4A, Figure 4—figure supplement 1). We genotyped 9 Col/Ct markers on chromosome 3 and observed that 5 markers in proximity to CEN3 were significantly more homozygous in the hottest compared to the coldest F₂ quartile (FDR corrected chi-square test p = 1.14 × 10⁻⁰⁷) (Figure 4D–F, Figure 4—figure supplement 1 and Table 7). The physical extent of the effect was less (2.62 Mb) on the long arm of chromosome 3 for CEN3 than observed for 420 and I2f, potentially due to heterozygosity effects acting independently from both arms across the centromere. Together this shows that juxtaposition of heterozygous and homozygous regions in various chromosomal locations can modify local crossover frequency.

Figure 5 with 1 supplement see all

Download asset Open asset

Juxtaposition of heterozygous and homozygous regions triggers reciprocal crossover remodelling.

(A) Schematic diagram illustrating the physical location of *420* and *I3bc* transgenes expressing fluorescent proteins in seed and pollen. Beneath are diagrams illustrating the locations of Col/Col homozygous (red) and Col/Ct heterozygous (black) regions along chromosome 3. Positions of Col/Ct genotyping markers are indicated by blue ticks along the axis of the chromosome. Printed alongside are formulae for the calculation of genetic distance (cM) and crossover interference using *I3bc*. Counts of pollen with different combinations of fluorescence are indicated. For example, N_BYR indicates the number of pollen with blue, yellow and red fluorescence. (B) *I3b* and *I3c* genetic distance (cM) measured in HOM-HOM and HET-HOM plants as illustrated in (A). See Figure 5—source data 1. (C) As for (B) but showing calculation of crossover interference (1-CoC). See Figure 5—source data 2.

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

Figure 5—source data 1 Three colour I3bc FTL flow cytometry count data.: https://doi.org/10.7554/eLife.03708.025
Download elife-03708-fig5-data1-v2.docx
Figure 5—source data 2 Three colour I3bc FTL flow cytometry count data–measurement of crossover interference.: https://doi.org/10.7554/eLife.03708.026
Download elife-03708-fig5-data2-v2.docx

Table 7

Chromosome 3 genotype counts from hot and cold quartile CEN3/++ Col/Ct F₂ individuals

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

Marker coordinates (bp)	Hot quartile HET	Hot quartile HOM	Cold quartile HET	Cold quartile HOM	FDR P
259000	16	14	17	13	1
2718000	16	14	18	12	1
5352000	19	11	17	13	1
7674000	20	10	12	18	0.129
8495000	23	7	13	17	0.0389
9404000	26	4	16	14	0.0308
11115724	30	0	30	0	1
16520560	30	0	30	0	1
21008000	27	3	14	16	0.00477
22076000	23	7	12	18	0.0308
23040000	24	6	10	20	0.00477

The number of CEN3/++ Col/Ct F₂ individuals showing Col homozygosity (HOM) or Col/Ct heterozygosity (HET) for the indicated markers, in either the hottest or coldest quartile. The p value was obtained by performing a chi square test between homozygous and heterozygous marker genotype counts in the hottest and coldest quartiles (2 × 2 contingency table), followed by FDR correction for multiple testing.

Juxtaposed heterozygous and homozygous regions show reciprocal changes in crossover frequency

We reasoned that if heterozygous regions increase recombination when juxtaposed with homozygous regions, then the linked homozygous regions may show compensatory decreases, due to crossover interference (Copenhaver et al., 2002; Zhang et al., 2014a). To test this idea we constructed a three-colour pollen FTL interval termed I3bc that overlaps the 420 seed interval on chromosome 3 (Figure 5 and Table 2). Three-colour FTL configurations allow simultaneous measurement of crossover frequency in adjacent intervals and measurement of crossover interference (Berchowitz and Copenhaver, 2008; Yelina et al., 2013) (Figure 5—figure supplement 1). To calculate interference, the observed double crossover (DCO) classes (N_-Y- + N_B-R) are compared to the number expected in the absence of interference: (I3b cM/100) × (I3c cM/100) × N_total (Figure 5A). The Coefficient of Coincidence (CoC) is calculated by dividing Observed DCOs by Expected DCOs, and interference strength calculated as 1-CoC (Figure 5A).

I3bc wild type genetic distance was greater than that measured from 420 self-fertilization data, as expected due to increases observed in sub-telomeric regions in male meiosis (Table 2—Figure 5—source data 1) (Giraut et al., 2011). I3b crossover frequency was also higher than I3c, consistent with a telomeric gradient in male crossover frequency (Figure 5B and Table 2) (Giraut et al., 2011). We compared crossovers in plants that were entirely Col homozygous (HOM-HOM) vs plants that were Col/Ct heterozygous within I3b, but Col/Col homozygous in I3c and for the rest of chromosome 3 (HET-HOM) (Figure 5A). Dense genotyping markers were used to confirm the location of homozygous and heterozygous regions (Figure 5A). We observed that I3b crossovers significantly increased in HET-HOM compared to HOM-HOM plants, and there was a reciprocal decrease in I3c crossovers (Figure 5B, Figure 5—source data 2) (both GLM p < 2.0 × 10⁻¹⁶). Together this is consistent with reciprocal crossover changes in juxtaposed heterozygous and homozygous regions being driven by crossover interference.

Reciprocal crossover remodeling across heterozygosity/homozygosity junctions requires interference

The effect of heterozygosity/homozygosity juxtaposition on crossovers extends over megabase distances, which is similar to the scale of crossover interference in Arabidopsis (Copenhaver et al., 2002; Giraut et al., 2011; Salomé et al., 2012). We therefore next used mutations in meiotic recombination pathways to analyse the genetic requirements of these effects. Specifically, we generated plants carrying the linked chromosome 3 fluorescent crossover reporters 420 and CEN3 (420-CEN3), with varying Col/Ct genotype and that were wild type, fancm or fancm zip4 (Figure 6—Figure 6—figure supplement 1). Crossover frequency in 420 and CEN3 can be scored in the same individuals, as these intervals use fluorescent proteins expressed in seed and pollen respectively. In fancm DSBs that would normally be repaired as non-crossovers enter a non-interfering pathway leading to a substantial increase in crossovers, although the interfering pathway remains active (Crismani et al., 2012). In fancm zip4 only non-interfering crossovers occur, due to mutation of the ZMM gene ZIP4 (Chelysheva et al., 2007; Crismani et al., 2012). In wild type, both interfering and non-interfering pathways are active, but interfering crossovers predominate and constitute ∼85% of total crossovers (Copenhaver et al., 2002; Higgins et al., 2004; Mercier et al., 2005). Therefore, by comparing genetic distances in wild type, fancm and fancm zip4, where the relative proportions of interfering and non-interfering repair vary dramatically, we can investigate the sensitivity of different recombination pathways to heterozygosity.

Figure 6 with 1 supplement see all

Download asset Open asset

Genetic requirements of crossover remodelling via juxtaposition of heterozygous and homozygous regions.

(A–D) Replicate measurements of *420* (red) and *CEN3* (blue) genetic distances (cM) are plotted in wild type, *fancm* and *fancm zip4*. See Figure 6—source data 1, 2. Chromosome 3 genotypes of the plants analysed are indicated above the plots (green = Col and red = Ct), for example, HET-HOM indicates heterozygous within *420* and homozygous outside.

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

Figure 6—source data 1 420 fluorescent seed count data from wild type, fancm and fancm zip4 individuals with varying heterozygosity.: https://doi.org/10.7554/eLife.03708.030
Download elife-03708-fig6-data1-v2.docx
Figure 6—source data 2 CEN3 flow cytometry count data from wild type, fancm and fancm zip4 individuals with varying heterozygosity.: https://doi.org/10.7554/eLife.03708.031
Download elife-03708-fig6-data2-v2.docx

When chromosome 3 is Col/Col homozygous (HOM-HOM) genetic distance in the 420 interval significantly increased in fancm and fancm zip4 mutants compared with wild type (both GLM p < 2.0 × 10⁻¹⁶) (Figure 6A, Figure 6—source data 1), consistent with repair of the majority of DSBs via a non-interfering crossover pathway (Crismani et al., 2012). However, the CEN3 interval experienced a smaller yet significant increase in genetic distance in fancm and decreased in fancm zip4 (both GLM p < 2.0 × 10⁻¹⁶), indicating that non-interfering crossover repair is less efficient in this region (Figure 6A, Figure 6—source data 2). We next generated plants that were Col/Ct heterozygous (HET-HET) on chromosome 3 and observed that the previous increase in 420 crossovers was strongly suppressed in fancm and fancm zip4 (GLM p = 1.24 × 10⁻⁰⁶ and p < 2.0 × 10⁻¹⁶), whereas wild type Col/Ct were slightly but significantly higher than wild type Col/Col (GLM p = 0.0126) (Figure 6A–B). CEN3 crossovers were also significantly suppressed by Col/Ct heterozygosity in fancm and nearly eliminated in fancm zip4 compared to Col/Col (both GLM p < 2.0 × 10⁻¹⁶) (Figure 6A–B). Together this indicates that the non-interfering crossover repair pathway that predominates in fancm and fancm zip4 is less efficient in heterozygous regions and particularly within the centromeric region, which shows high polymorphism levels (Table 2).

We next tested the effect of juxtaposing heterozygous and homozygous regions in fancm and fancm zip4 mutants. We first generated lines that were Col/Ct heterozygous within 420 and Col/Col homozygous outside (HET-HOM) (Figure 6—figure supplement 1). As expected, wild type HET-HOM lines show a significant increase in 420 and a reciprocal decrease in CEN3 crossovers compared to wild type HOM-HOM (both GLM p < 2.0 × 10⁻¹⁶) (Figure 6A,C), indicating compensatory changes between the two intervals in the HET-HOM lines. As the HET-HOM lines are heterozygous within 420, this again inhibited crossovers in fancm compared to fancm HOM-HOM (GLM p = 2.38 × 10⁻¹⁵) (Figure 6A,C). HET-HOM lines in fancm zip4 showed lower 420 crossovers than wild type HOM-HOM (GLM p < 2.0 × 10⁻¹⁶), which demonstrates that the interfering pathway is required for the heterozygosity-homozygosity juxtaposition effect (Figure 6A,C). We also generated HOM-HET lines that were homozygous within 420 and heterozygous outside, which significantly reduced 420 crossovers compared to wild type HOM-HOM as expected (GLM p = 7.60 × 10⁻¹¹) (Figure 6A,D). HOM-HET lines in fancm and fancm zip4 showed high 420 crossovers comparable to HOM-HOM, as the non-interfering crossover repair active in these backgrounds is efficient in homozygous regions (Figure 6A,D). CEN3 genetic distance was again strongly suppressed in fancm and fancm zip4 HOM-HET lines compared with HOM-HOM (both GLM p < 2.0 × 10⁻¹⁶), consistent with heterozygosity inhibiting non-interfering crossover repair (Figure 6A,D). Together these data demonstrate that juxtaposition of heterozygous and homozygous regions causes reciprocal changes in crossover frequency via interference.

Total chiasmata are maintained when heterozygosity is varied

As we observed regional changes in crossover frequency with varying patterns of heterozygosity, we next sought to test whether total recombination events were different. When homologous chromosomes align on the metaphase-I plate, crossovers can be cytologically visualized as chiasmata (Sanchez-Moran et al., 2002). To estimate the number of crossovers per meiotic nucleus we performed chromosome spreads of pollen mother cells (PMCs), followed by fluorescence in situ hybridization using a 45S rDNA probe (Figure 7, Figure 7—source data 1). We counted total chiasmata in metaphase-I nuclei in Col/Col homozygotes, Ct/Ct homozygotes and Col/Ct F₁ heterozygotes. In addition, we counted chiasmata in recombinant 420-CEN3 lines showing high (HET-HOM, 27.96 cM) and low (HOM-HET, 13.83 cM) 420 crossover frequency (Figure 7C,D). Adjacent chiasmata count categories were combined to give a minimum expected value of five for the purposes of a chi-square test with 8° of freedom. This test gave no significant differences in chiasmata between the genotypes (p = 0.3365) (Figure 7). Together this is consistent with homeostatic maintenance of crossover numbers, despite local crossover changes caused by juxtaposition of heterozygous and homozygous regions.

Figure 7

Download asset Open asset

Total chiasmata frequencies are stable between Col, Ct and recombinant lines.

(A–E) Metaphase-I chromosome spreads from anthers from (A) Col/Col *420*, (B) Ct/Ct, (C) Col × Ct F₁, (D) a Col × Ct *420* (HOM-HET) cold recombinant line and (E) a Col × Ct *420* hot (HET-HOM) recombinant line. DNA is stained with DAPI (blue) and labelled with a *45S* rDNA probe (green). Scale bars = 10 μM. (F) Boxplot showing total number of chiasmata per nucleus for each genotype. See Figure 7—source data 1.

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

Figure 7—source data 1 Chiasmata count data.: https://doi.org/10.7554/eLife.03708.034
Download elife-03708-fig7-data1-v2.docx

Crossover interference increases in heterozygous regions

Our analysis of 420-CEN3 recombination rates implicated interference as driving crossover changes across homozygosity/heterozygosity junctions. We therefore sought to directly measure interference in lines with varying heterozygosity. We generated I3bc lines that varied in Col/Ct genotype and that were wild type, fancm, zip4 or fancm zip4 (Figure 8—figure supplement 1). We first compared I3bc plants that were Col/Col homozygous (HOM-HOM) with Col/Ct heterozygotes (HET-HET). In wild type, genetic distances did not significantly change between HOM-HOM and HET-HET (GLM p = 0.352 and p = 0.666), but crossover interference significantly increased (GLM p < 2.0 × 10⁻¹⁶) (Figure 8A,B, Figure 8—source data 1). Consistent with previous observations, fancm and fancm zip4 showed a significant reduction and an absence of interference respectively, in a HOM-HOM background (GLM p < 2.0 × 10⁻¹⁶ and p = 4.94 × 10⁻¹⁶) (Figure 8A, Figure 8—source data 2) (Crismani et al., 2012; Yelina et al., 2013). In HET-HET plants the crossover frequency increases seen in fancm and fancm zip4 were again greatly suppressed, or eliminated, relative to HOM-HOM, as observed previously for 420-CEN3 (GLM both p < 2.0 × 10⁻¹⁶) (Figure 8B). Unexpectedly, interference measurements significantly increased in both fancm and fancm zip4 mutants in a HET-HET background compared to HOM-HOM (GLM p < 2.0 × 10⁻¹⁶ and p = 4.94 × 10⁻¹⁶) (Figure 8B). We propose that in the absence of the ZMM pathway alternative repair pathways exist which are differentially sensitive to polymorphism and interference. Multiple, redundant repair pathways are consistent with the residual crossovers observed in msh4 mus81 double mutants (Higgins et al., 2008b). Finally, we measured I3bc cM in zip4 mutants alone (HOM-HOM) and observed significantly decreased crossovers compared with wild type HOM-HOM (GLM p < 2.0 × 10⁻¹⁶) (Figure 8E, Figure 8—source data 1). Importantly, zip4 genetic distances were further significantly reduced when comparing HOM-HOM to HET-HET backgrounds (GLM p = 1.79 × 10⁻¹⁰ and p = 1.53 × 10⁻⁹) (Figure 8E). This provides additional evidence that the non-interfering repair pathway remaining in zip4 is inefficient in heterozygous regions. Interference measurements using I3bc are reliant on the relatively rare double crossover classes (N_-Y- + N_B-R) (Figure 5A). Due to low zip4 fertility it was difficult to obtain sufficient DCO counts to make reliable interference measurements, although the observed counts are consistent with an absence of interference in this mutant (Figure 8—source data 4).

Figure 8 with 1 supplement see all

Download asset Open asset

Crossover interference increases when heterozygous and homozygous regions are juxtaposed.

(A–D) Replicate measurements of *I3b* and *I3c* genetic distances (cM), and *I3bc* crossover interference are plotted in wild type, *fancm, fancm zip4* and *zip4*. Black dots represent replicate measurements with mean values indicated by red dots. Chromosome 3 genotypes of the plants analysed are indicated above the plots (green = Col and red = Ct), for example, HET-HOM indicates heterozygous within *I3bc* and homozygous outside. See Figure 8—source data 1, 2. (E) *I3b* and *I3c* genetic distances (cM) are plotted in wild type and *zip4* mutants with varying patterns of heterozygosity, labelled as for (A–D). Mean values between samples are connected with red lines. See Figure 8—source data 3, 4.

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

Figure 8—source data 1 I3bc fluorescent seed count data from wild type, fancm and fancm zip4 individuals with varying heterozygosity.: https://doi.org/10.7554/eLife.03708.036
Download elife-03708-fig8-data1-v2.docx
Figure 8—source data 2 Calculation of I3bc interference from wild type, fancm and fancm zip4 individuals with varying heterozygosity.: https://doi.org/10.7554/eLife.03708.037
Download elife-03708-fig8-data2-v2.docx
Figure 8—source data 3 I3bc fluorescent seed count data from wild type and zip4 individuals with varying heterozygosity.: https://doi.org/10.7554/eLife.03708.038
Download elife-03708-fig8-data3-v2.docx
Figure 8—source data 4 Calculation of I3bc interference in wild type and zip4.: https://doi.org/10.7554/eLife.03708.039
Download elife-03708-fig8-data4-v2.docx

To test the effects of heterozygosity/homozygosity juxtaposition we next generated lines that were Col/Ct heterozygous within I3bc and Col/Col homozygous outside (HET-HOM). As expected, wild type I3b and I3c genetic distances both significantly increase in HET-HOM lines relative to HOM-HOM (GLM both p < 2.0 × 10⁻¹⁶), consistent with our previous 420 experiments, and this was associated with a significant increase in crossover interference (GLM p < 2.0 × 10⁻¹⁶) (Figure 8A,C). As shown earlier, we observed that Col/Ct (HET-HOM) heterozygosity suppressed the crossover increases seen in fancm and fancm zip4 (GLM p < 2.0 × 10⁻¹⁶), with the same significant increases in crossover interference strength (GLM p < 2.0 × 10⁻¹⁶ and p = 4.94 × 10⁻¹⁶) (Figure 8A,C). The reciprocal situation was observed in HOM-HET plants where I3bc was Col/Col homozygous and the rest of the chromosome Col/Ct heterozygous. I3b and I3c genetic distances were significantly decreased in wild type HOM-HET compared with wild type HOM-HOM plants (GLM both p < 2.0 × 10⁻¹⁶) (Figure 8A,D). HOM-HET fancm and fancm zip4 plants showed high crossovers, as the non-interfering pathway is efficient in the homozygous I3bc interval (Figure 8A,D). We also generated HET-HOM zip4 lines, which unlike wild type showed significantly lower I3b and I3c cM than HOM-HOM zip4 (GLM both P= p < 2.0 × 10⁻¹⁶) (Figure 8E). This again demonstrates that crossover remodelling at heterozygosity/homozygosity junctions requires interference and that non-interfering repair is inefficient in heterozygous regions.

As an independent test of the effect of heterozygosity on crossover interference we analysed four three-colour FTL intervals distributed throughout the genome (Figure 1A and Table 2). We measured crossover frequency and interference in Col/Col homozygotes vs Col/Ler F₁ heterozygotes using meiotic pollen tetrads (Tables 8, 9). This approach is possible as the FTL crossover reporter system was generated in the qrt1-2 mutant background, where sister pollen grains remain physically attached as meiotic tetrads (Berchowitz and Copenhaver, 2008). We scored a total of 49,801 tetrads for Col/Col (an average of 6225 per interval) and 42,422 tetrads for Col/Ler (an average of 5302 per interval) (Tables 8, 9). Compared to Col/Col, genetic distance significantly decreased in Col/Ler for six of the eight intervals measured and the remaining two intervals were not significantly changed (Table 8). To calculate interference strength we compared cM values in each interval from tetrads that had a crossover in the adjacent interval, to the same intervals in tetrads lacking a crossover in the adjacent interval, and detected significant positive interference in all cases (Table 9) (Berchowitz and Copenhaver, 2008). The resulting interference ratios were then compared between Col/Col and Col/Ler using Fisher's combined probability test, which revealed a significant increase in interference strength in Col/Ler (χ²_.001[16] = 39.26) (Table 9). Therefore, the effect of heterozygosity increasing the interference strength is evident in both Col × Ct and Col × Ler crosses.

Table 8

Tetrad FTL cM data in Col/Col and Col/Ler backgrounds

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

	Col/Col				Col/Ler
Interval	PD	NPD	T	cM*	PD	NPD	T	cM*
1b	3976	3	742	8.05 ± 0.29	4395	2	652	6.58 ± 0.25†
1c	3022	11	1695	18.62 ± 0.04	3156	18	1891	19.73 ± 0.04
2a	6787	2	430	3.06 ± 0.15	5920	0	283	2.28 ± 0.13†
2b	6582	2	635	4.48 ± 0.18	5796	0	407	3.28 ± 0.16†
3b	4363	22	2557	19.37 ± 0.35	2758	2	1056	13.99 ± 0.38†
3c	6185	5	736	5.53 ± 0.21	3576	2	238	3.28 ± 0.22†
5c	5356	1	666	5.58 ± 0.21	5458	0	676	5.51 ± 0.20
5d	5358	1	664	5.56 ± 0.21	5540	2	594	4.94 ± 0.20†

*

Map distance in cM (±S.E.).
†

Significant difference in map distance in the heterozygous Col/Ler background compared to the same interval in the Col/Col homozygous background.

Table 9

Tetrad FTL crossover interference data in Col/Col and Col/Ler backgrounds

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

	Col/Col			Col/Ler
Interval	W/o adj. CO*	w/ adj. CO*	R1†	W/o adj. CO*	w/ adj. CO*	R2†
1b	10.69 ± 0.40	3.31 ± 0.30‡	3.23	9.78 ± 0.37	1.22 ± 0.18‡	8.04§
1c	20.61 ± 0.45	7.92 ± 0.76‡	2.6	22.13 ± 0.46	3.52 ± 0.50‡	6.29§
2a	3.20 ± 0.16	1.18 ± 0.30‡	2.75	2.42 ± 0.14	0.37 ± 0.21‡	6.55
2b	4.65 ± 0.19	1.74 ± 0.44‡	2.68	3.41 ± 0.16	0.53 ± 0.30‡	6.44
3b	20.84 ± 0.37	6.95 ± 0.82‡	2.3	14.73 ± 0.40	2.92 ± 0.76‡	5.05
3c	7.65 ± 0.30	1.90 ± 22‡	4.03	4.28 ± 0.30	0.66 ± 0.18‡	6.46
5c	5.87 ± 0.23	3.23 ± 0.47‡	1.82	5.85 ± 0.22	2.35 ± 0.43‡	2.49
5d	5.85 ± 0.23	3.22 ± 0.48‡	1.82	5.29 ± 0.22	2.07 ± 0.38‡	2.56

*

Map distances in cM (±S.E.) for intervals with and without adjacent crossovers (CO).
†

Ratios of map distances for intervals with and without adjacent crossovers in homozygous Col/Col (R1) and heterozygous Col/Ler (R2) backgrounds.
‡

Significant difference in map distances in intervals when adjacent interval does or doesn't have a CO.
§

Significant difference between R2 and R1.

Discussion

We demonstrate reciprocal crossover increases and decreases when heterozygous and homozygous regions are juxtaposed and further demonstrate that this process requires crossover interference. The mechanism of interference is presently unclear, but a Beam-Film model has been developed where crossovers are patterned via forces similar to mechanical stress and which predicts experimental data (Kleckner et al., 2004; Zhang et al., 2014a, 2014b). In this model each chromosome begins with an array of precursor interhomolog strand invasion events, one of which becomes crossover designated via a stress-related force (Designation Driving Force DDF). This causes a local reduction and redistribution of stress in both directions that dissipates with increasing distance (Kleckner et al., 2004; Zhang et al., 2014a, 2014b). At the point where stress increases sufficiently precursor events can again become crossover designated. Any remaining precursors then mature into other fates including non-crossovers and non-interfering crossovers (Kleckner et al., 2004; Zhang et al., 2014a, 2014b).

We considered the effect of juxtaposition of heterozygous/homozygous regions in the context of the Beam-Film model (Kleckner et al., 2004; Zhang et al., 2014a, 2014b). Detection of heterozygosity most likely occurs downstream of interhomolog strand invasion and the formation of base pair mismatches. Therefore, we assume that the initial distribution of meiotic DSBs is unchanged in homozygous or heterozygous states. Mismatches are observed to have a local inhibitory effect on meiotic crossovers (Dooner, 1986; Borts and Haber, 1987; Jeffreys and Neumann, 2005; Baudat and de Massy, 2007; Cole et al., 2010). Therefore, one possibility is that mismatched precursors in heterozygous regions are slowed in maturation and trigger feedback mechanisms that cause further DSBs, for example via ATM/ATR kinase signalling (Carballo et al., 2008; Lange et al., 2011; Zhang et al., 2011; Kurzbauer et al., 2012; Garcia et al., 2015). As a consequence, heterozygous regions would receive more ‘late’ DSBs, leading to more precursors and a higher chance of receiving a crossover designation event. An increased chance of crossover designation would lead to spreading of interference into adjacent homozygous regions causing reciprocal crossover decreases. An alternative model is that mismatched precursors are more sensitive to crossover designation and thus heterozygous regions have a higher chance of an interfering crossover, leading to similar effects. These potential models could be distinguished by measurement of non-crossover (NCO) levels, which should increase in heterozygous regions if more DSBs occur. Our data also demonstrate that non-interfering repair is less efficient in heterozygous regions, which will further contribute to the changes we see across homozygosity/heterozygosity junctions.

Sequence polymorphism has been observed to suppress crossover recombination at the hotspot (kilobase) scale in diverse eukaryotes (Dooner, 1986; Borts and Haber, 1987; Jeffreys and Neumann, 2005; Baudat and de Massy, 2007; Cole et al., 2010). For example, at the mouse A3 hotspot an indel polymorphism within an inverted repeat overlaps a crossover refractory zone (Cole et al., 2010). However, this zone forms significant numbers of non-crossovers, indicating that the repeat/indel does not inhibit DSB formation, but inhibits downstream progression to crossover recombination (Cole et al., 2010). In yeast addition of SNPs to the MAT-URA3 hotspot decreased crossovers and increased the frequency of gene conversions, further indicating that polymorphism can inhibit crossovers at fine-scale (Borts and Haber, 1987). Finally, intragenic mapping of the maize Bronze hotspot demonstrated that transposon insertions suppress crossovers more strongly than single nucleotide changes (Dooner, 1986; Fu et al., 2001; Dooner and He, 2008), again consistent with progression to crossover repair being inhibited by local sequence polymorphisms. Several heteroduplex joint molecules with distinct properties form during meiosis, including displacement-loops and dHJs (Keeney and Neale, 2006). It is possible that these joint molecules and their interactions with recombinases are sensitive to base-pair mismatches. The mismatch repair protein MutS directly recognizes mismatched base-pairs and serves as a paradigm for this type of function (Lamers et al., 2000; Obmolova et al., 2000).

The reciprocal crossover changes we observe when heterozygous regions are juxtaposed with homozygous regions are reminiscent of other homeostatic effects characterized during meiosis (Hillers and Villeneuve, 2003; Martini et al., 2006; Robine et al., 2007; Libuda et al., 2013; Thacker et al., 2014). For example, multiple levels of interference have been detected in mice (de Boer et al., 2006; Cole et al., 2012), Zip3 foci with distinct timing and properties are observed in budding yeast (Serrentino et al., 2013), and ‘upstream’ DSB patterns are altered in ‘downstream’ ZMM mutants (Thacker et al., 2014). As plants, fungi and mammals share the presence of interfering and non-interfering crossover repair pathways similar effects over heterozygosity/homozygosity junctions may be generally important (Stahl et al., 2004). However, when assessing the significance of such effects it is also important to consider how outcrossing vs selfing will influence patterns of homozygosity and heterozygosity within different species. Together our data show how varying patterns of sequence polymorphism along chromosomes can have a significant effect on distributions of meiotic recombination.

Materials and methods

Measuring crossovers using two-colour fluorescence microscopy of seed and flow cytometry of pollen

Request a detailed protocol

Flow cytometry of pollen can be used to rapidly measure meiotic segregation of heterozygous transgenes encoding distinct colours of fluorescent protein (Yelina et al., 2012, 2013). cM were calculated from flow cytometry data using the formula:

cM = 100 \times (R 5 / (R 3 + R 5)),

Where R5 is a number of green-alone fluorescent pollen grains and R3 is a number of green and red fluorescent pollen grains (Yelina et al., 2012, 2013). We previously observed that the number of red-alone pollen exceeded that of green-alone pollen when lines heterozygous for both eYFP and dsRed (eYFPDsRed/++) were analysed (Yelina et al., 2012, 2013). Using pulse-width/SSC (side scatter) analysis and back-gating we demonstrated that the excess counts come primarily from non-hydrated pollen (Yelina et al., 2012, 2013). Therefore to avoid this artifact we multiply the green-alone counts by two to obtain the number of recombinant pollen.

To increase measurement throughput using fluorescent seed we adapted CellProfiler image analysis software (Carpenter et al., 2006) (Figure 2). This program identifies seed boundaries in micrographs and assigns a RFP and GFP fluorescence intensity to each seed object (Figure 2A–B). Three pictures of the seed are acquired at minimum magnification (×0.72) using a charge coupled device (CCD) camera; (i) brightfield, (ii) UV through a dsRed filter and (iii) UV through a GFP filter (Figure 2A). As seed are diploid, there are nine possible fluorescent genotypes when a RFP-GFP/++ heterozygote is self-fertilized, in contrast to four possible states for haploid pollen (Yelina et al., 2013) (Figure 2E). Histograms of seed fluorescence can be used to classify fluorescent and non-fluorescent seed for each colour (Figure 2C–D). Although it is possible to distinguish seed with one vs two T-DNA copies, there is greater overlap between the groups (Figure 2C–E). Therefore, we use fluorescent vs non-fluorescent seed counts for crossover measurement. Using this method it is possible to score 2000–6000 meioses per self-fertilized individual. When plants have been self-fertilized, genetic distance is calculated using the formula:

cM = 100 \times (1 - {[1 - 2 (N_{G} + N_{R}) / N_{T}]}^{1 / 2}),

Where N_G is a number of green-alone fluorescent seeds, N_R is a number of red-alone fluorescent seed and N_T is the total number of seeds counted. During generation of 420/++ F₂ populations we selected for individuals that are heterozygous for transgenes expressing red and green fluorescent proteins (RFP-GFP/++). The majority of these individuals receive a chromosome with linked RFP and GFP transgenes over a non-transgenic chromosome (RFP-GFP/++) (Figure 2—figure supplement 1). In a minority of cases F₂ plants receive recombined RFP-+ and +-GFP chromosomes (Figure 2—figure supplement 1). In the progeny of these individuals the fluorescent seed classes representing parental and crossover genotypes are reversed (Figure 2—figure supplement 1). As R+/+G plants also have variable heterozygosity/homozygosity patterns within 420 depending on crossover positions we excluded these plants from further analysis.

To test whether recombinant and non-recombinant counts were significantly different between replicate groups we used a GLM. We assumed the count data is binomially distributed:

Y_{i} \sim B (n_{i}, p_{i}),

where Y_i represents the recombinant counts, $n_{i}$ are the total counts, and we wish to model the proportions Y_i/n_i. Then:

E (Y_{i} / n_{i}) = p_{i},

and

var (Y_{i} / n_{i}) = \frac{p_{i (1 - p_{i})}}{n_{i}} .

Thus, our variance function is:

V (μ_{i}) = μ_{i} (1 - μ_{i}),

and our link function must map from (0,1) → (−∞, ∞). We used a logistic link function which is:

g (μ_{i}) = logit (μ_{i}) = \log \frac{μ_{i}}{1 - μ_{i}} = β X + ε_{i},

where $ε_{ι} \sim N (0, σ^{2})$ . Both replicates and genotypes are treated as independent variables (X) in our model. We considered p values less than 0.05 as significant.

Measuring crossovers and interference using three-colour flow cytometry of pollen

Request a detailed protocol

Measurements of interference within the I3bc interval were carried out as described previously with minor modifications (Yelina et al., 2013). Inflorescences were collected in polypropylene tubes and pollen was extracted by vigorous shaking in 30 ml of freshly prepared pollen sorting buffer (PSB: 10 mM CaCl₂, 1 mM KCl, 2 mM MES, 5% wt/vol sucrose, 0.01% Triton X-100, pH 6.5). The pollen suspension was filtered through a 70 µM cell strainer to a fresh 50 ml polypropylene tube and centrifuged at 450×g for 3 min. The supernatant was removed and the pollen pellet washed once with 20 ml of PSB without Triton. The pollen suspension was centrifuged at 450×g for 3 min and the supernatant discarded and the pollen pellet resuspended in 500 µl of PSB without Triton. A CyAn ADP Analyser (Beckman Coulter, California, USA) equipped with 405 nm and 488 nm lasers and 530/40 nm, 575/25 nm and 450/50 nm band-pass filters was used to analyse the samples. Polygons were used for gating pollen populations and for each sample eight pollen class counts were obtained (Figure 5—figure supplement 1). I3b and I3c genetic distances were calculated using the following formula:

N_{total} = (N_{- Y -} + N_{B - R} + N_{- YR} + N_{B - -} + N_{BY -} + N_{- - R} + N_{BYR} + N_{- - -})

I 3 b cM = (N_{- Y -} + N_{B - R} + N_{- YR} + N_{B - -}) / N_{total}

I 3 c cM = (N_{- Y -} + N_{B - R} + N_{BY -} + N_{- - R}) / N_{total},

where N_-Y-, N_B-R, N_-YR, N_B--, N_BY-, N_--R, N_BYR, and N_--- are pollen grain counts in each of the eight populations (Figure 5—figure supplement 1). For example, N_BYR is the number of pollen that were blue, yellow and red fluorescent.

Crossover interference was calculated using the following formulas:

Observed DCOs = (N_{- Y -} + N_{B - R}),

Expected DCOs = (I 3 b cM / 100) \times (I 3 c cM / 100) \times N_{total},

Coefficient of Coincidence = Observed DCOs / Expected DCOs,

Interference = 1 - CoC.

At least three biological replicates, constituting 3–5 individual plants were analysed for each sample (Yelina et al., 2013). Statistical tests for genetic distances were performed as described above using a GLM. To test for significant differences in interference we compared observed and expected double crossovers using the same approach.

Generation of fancm and fancm zip4 Col/Ct mapping populations with varying heterozygosity

Request a detailed protocol

Col-0 420 and Ct-1 lines were crossed to fancm-1 zip4-2 double mutant lines in the Col-0 background (Crismani et al., 2012) (Figure 6—figure supplement 1). The resulting F₁ plants were crossed together and progeny identified that were fancm zip4 heterozygous, and 420/++ Col/Ct heterozygous on chromosome 3. Chromosome 3 genotypes were tested in all cases using 13 Col/Ct indel markers (Supplementary file 1). These plants were self-fertilized and 420 homozygous individuals identified (all seed were red and green fluorescent) that were also Ct homozygous outside of 420 and that were fancm zip4 heterozygous (Figure 6—figure supplement 1 and Figure 8—figure supplement 1). These plants were then crossed to CEN3 or I3bc in wild type, fancm and fancm zip4 mutants to obtain scorable progeny with a HOM-HET genotype (Figure 6—figure supplement 1). The selfed progeny of 420/++ Col/Ct fancm zip4 heterozygous plants were also selected for plants with no fluorescent T-DNAs and either chromosome 3 in a Ct homozygous state, or with Ct homozygosity within 420 and Col homozygosity outside (Figure 6—figure supplement 1). These plants were crossed with doubly marked 420-CEN3 or I3bc lines in either wild type, fancm or fancm zip4 mutant backgrounds to obtain HET-HET and HET-HOM scorable plants respectively (Figure 6—figure supplement 1 and Figure 8—figure supplement 1). Equivalent genetic crosses were performed during analysis of I3bc (Figure 8—figure supplement 1). At least three independent lines were generated and analysed for each combination, apart from HOM-HET 420-CEN3 where two were analysed.

To genotype zip4-2 (Salk_068052) the following primers were used:

zip4-2-F 5′-TTGCTACCTTGGGCTCTCTC-3′

zip4-2-R 5′-ATTCTGTTCTCGCTTTCCAG-3′

LBb1.3 5′-ATTTTGCCGATTTCGGAAC-3′

The resulting PCR products were ∼680 bp for wild type (zip4-2-F + zip4-2-R) and ∼340 bp for zip4-2 mutant (zip4-2-F + Lbb1.3) (Crismani et al., 2012).

To genotype the fancm mutation we amplified using the following primers:

fancm1dCAPsF1 5′-ACAATATATGTTTCGTGCAGGTAAGACATTGGAAG-3′

fancm1dCAPsR1 5′-CACCAATAGATGTTGCGACAAT-3′

The resulting PCR product was digested with MboII, which yields a ∼215 bp product for wild type and ∼180 bp for fancm (Crismani et al., 2012).

Chiasmata counting

Request a detailed protocol

Chiasmata counting was performed as previously described (Sanchez-Moran et al., 2002).

References

(1989)
Reduced variation in the yellow-achaete-scute region in natural populations of Drosophila melanogaster

Genetics 122:607–615.
- Google Scholar
1. Allers T
2. Lichten M
(2001) Differential timing and control of noncrossover and crossover recombination during meiosis
Cell 106:47–57.

https://doi.org/10.1016/S0092-8674(01)00416-0
- Google Scholar
(2001) Influence of genetic background and heterozygosity on meiotic recombination in Arabidopsis thaliana
Genome 44:971–978.

https://doi.org/10.1139/g01-094
- Google Scholar
1. Barton NH
2. Charlesworth B
(1998) Why sex and recombination?
Science 281:1986–1990.

https://doi.org/10.1126/science.281.5385.1986
- Google Scholar
1. Baudat F
2. Buard J
3. Grey C
4. Fledel-Alon A
5. Ober C
6. Przeworski M
7. Coop G
8. de Massy B
(2010) PRDM9 is a major determinant of meiotic recombination hotspots in humans and mice
Science 327:836–840.

https://doi.org/10.1126/science.1183439
- Google Scholar
1. Baudat F
2. de Massy B
(2007) Cis- and trans-acting elements regulate the mouse Psmb9 meiotic recombination hotspot
PLOS Genetics 3:e100.

https://doi.org/10.1371/journal.pgen.0030100
- Google Scholar
1. Bauer E
2. Falque M
3. Walter H
4. Bauland C
5. Camisan C
6. Campo L
7. Meyer N
8. Ranc N
9. Rincent R
10. Schipprack W
11. Altmann T
12. Flament P
13. Melchinger AE
14. Menz M
15. Moreno-González J
16. Ouzunova M
17. Revilla P
18. Charcosset A
19. Martin OC
20. Schön CC
(2013) Intraspecific variation of recombination rate in maize
Genome Biology 14:R103.

https://doi.org/10.1186/gb-2013-14-9-r103
- Google Scholar
1. Begun DJ
2. Aquadro CF
(1992) Levels of naturally occurring DNA polymorphism correlate with recombination rates in D. melanogaster
Nature 356:519–520.

https://doi.org/10.1038/356519a0
- Google Scholar
1. Berchowitz LE
2. Copenhaver GP
(2008) Fluorescent Arabidopsis tetrads: a visual assay for quickly developing large crossover and crossover interference data sets
Nature Protocols 3:41–50.

https://doi.org/10.1038/nprot.2007.491
- Google Scholar
(2007) The role of AtMUS81 in interference-insensitive crossovers in A. thaliana
PLOS Genetics 3:10.

https://doi.org/10.1371/journal.pgen.0030132
- Google Scholar
1. Berg IL
2. Neumann R
3. Lam K-WG
4. Sarbajna S
5. Odenthal-Hesse L
6. May CA
7. Jeffreys AJ
(2010) PRDM9 variation strongly influences recombination hot-spot activity and meiotic instability in humans
Nature Genetics 42:859–863.

https://doi.org/10.1038/ng.658
- Google Scholar
(1997) An atypical topoisomerase II from Archaea with implications for meiotic recombination
Nature 386:414–417.

https://doi.org/10.1038/386414a0
- Google Scholar
1. Bishop DK
2. Park D
3. Xu L
4. Kleckner N
(1992) DMC1: a meiosis-specific yeast homolog of E. coli recA required for recombination, synaptonemal complex formation, and cell cycle progression
Cell 69:439–456.

https://doi.org/10.1016/0092-8674(92)90446-J
- Google Scholar
1. de Boer E
2. Stam P
3. Dietrich AJ
4. Pastink A
5. Heyting C
(2006) Two levels of interference in mouse meiotic recombination
Proceedings of the National Academy of Sciences of USA 103:9607–9612.

https://doi.org/10.1073/pnas.0600418103
- Google Scholar
1. Borts RH
2. Haber JE
(1987) Meiotic recombination in yeast: alteration by multiple heterozygosities
Science 237:1459–1465.

https://doi.org/10.1126/science.2820060
- Google Scholar
(2014) The relation between recombination rate and patterns of molecular evolution and variation in Drosophila melanogaster
Molecular Biology and Evolution 31:1010–1028.

https://doi.org/10.1093/molbev/msu056
- Google Scholar
1. Cao J
2. Schneeberger K
3. Ossowski S
4. Günther T
5. Bender S
6. Fitz J
7. Koenig D
8. Lanz C
9. Stegle O
10. Lippert C
11. Wang X
12. Ott F
13. Müller J
14. Alonso-Blanco C
15. Borgwardt K
16. Schmid KJ
17. Weigel D
(2011) Whole-genome sequencing of multiple Arabidopsis thaliana populations
Nature Genetics 43:956–963.

https://doi.org/10.1038/ng.911
- Google Scholar
(2008) Phosphorylation of the axial element protein Hop1 by Mec1/Tel1 ensures meiotic interhomolog recombination
Cell 132:758–770.

https://doi.org/10.1016/j.cell.2008.01.035
- Google Scholar
1. Carpenter AE
2. Jones TR
3. Lamprecht MR
4. Clarke C
5. Kang IH
6. Friman O
7. Guertin DA
8. Chang JH
9. Lindquist RA
10. Moffat J
11. Golland P
12. Sabatini DM
(2006) CellProfiler: image analysis software for identifying and quantifying cell phenotypes
Genome Biology 7:R100.

https://doi.org/10.1186/gb-2006-7-10-r100
- Google Scholar
(2007) Zip4/Spo22 is required for class I CO formation but not for synapsis completion in Arabidopsis thaliana
PLOS Genetics 3:e83.

https://doi.org/10.1371/journal.pgen.0030083
- Google Scholar
(2010) An easy protocol for studying chromatin and recombination protein dynamics during Arabidopsis thaliana meiosis: immunodetection of cohesins, histones and MLH1
Cytogenetic and Genome Research 129:143–153.

https://doi.org/10.1159/000314096
- Google Scholar
1. Chelysheva L
2. Vezon D
3. Chambon A
4. Gendrot G
5. Pereira L
6. Lemhemdi A
7. Vrielynck N
8. Le Guin S
9. Novatchkova M
10. Grelon M
(2012) The Arabidopsis HEI10 is a new ZMM protein related to Zip3
PLOS Genetics 8:e1002799.

https://doi.org/10.1371/journal.pgen.1002799
- Google Scholar
1. Chen C
2. Zhang W
3. Timofejeva L
4. Gerardin Y
5. Ma H
(2005) The Arabidopsis ROCK-N-ROLLERS gene encodes a homolog of the yeast ATP-dependent DNA helicase MER3 and is required for normal meiotic crossover formation
The Plant Journal 43:321–334.

https://doi.org/10.1111/j.1365-313X.2005.02461.x
- Google Scholar
1. Choi K
2. Zhao X
3. Kelly KA
4. Venn O
5. Higgins JD
6. Yelina NE
7. Hardcastle TJ
8. Ziolkowski PA
9. Copenhaver GP
10. Franklin FC
11. McVean G
12. Henderson IR
(2013) Arabidopsis meiotic crossover hot spots overlap with H2A.Z nucleosomes at gene promoters
Nature Genetics 45:1327–1336.

https://doi.org/10.1038/ng.2766
- Google Scholar
1. Clark RM
2. Schweikert G
3. Toomajian C
4. Ossowski S
5. Zeller G
6. Shinn P
7. Warthmann N
8. Hu TT
9. Fu G
10. Hinds DA
11. Chen H
12. Frazer KA
13. Huson DH
14. Schölkopf B
15. Nordborg M
16. Rätsch G
17. Ecker JR
18. Weigel D
(2007) Common sequence polymorphisms shaping genetic diversity in Arabidopsis thaliana
Science 317:338–342.

https://doi.org/10.1126/science.1138632
- Google Scholar
1. Cole F
2. Kauppi L
3. Lange J
4. Roig I
5. Wang R
6. Keeney S
7. Jasin M
(2012) Homeostatic control of recombination is implemented progressively in mouse meiosis
Nature Cell Biology 14:424–430.

https://doi.org/10.1038/ncb2451
- Google Scholar
1. Cole F
2. Keeney S
3. Jasin M
(2010) Comprehensive, fine-scale dissection of homologous recombination outcomes at a hot spot in mouse meiosis
Molecular Cell 39:700–710.

https://doi.org/10.1016/j.molcel.2010.08.017
- Google Scholar
1. Colomé-Tatché M
2. Cortijo S
3. Wardenaar R
4. Morgado L
5. Lahouze B
6. Sarazin A
7. Etcheverry M
8. Martin A
9. Feng S
10. Duvernois-Berthet E
11. Labadie K
12. Wincker P
13. Jacobsen SE
14. Jansen RC
15. Colot V
16. Johannes F
(2012) Features of the Arabidopsis recombination landscape resulting from the combined loss of sequence variation and DNA methylation
Proceedings of the National Academy of Sciences of USA 109:16240–16245.

https://doi.org/10.1073/pnas.1212955109
- Google Scholar
(2002)
Crossover interference in Arabidopsis

Genetics 160:1631–1639.
- Google Scholar
1. Copenhaver GP
2. Nickel K
3. Kuromori T
4. Benito MI
5. Kaul S
6. Lin X
7. Bevan M
8. Murphy G
9. Harris B
10. Parnell LD
11. McCombie WR
12. Martienssen RA
13. Marra M
14. Preuss D
(1999) Genetic definition and sequence analysis of Arabidopsis centromeres
Science 286:2468–2474.

https://doi.org/10.1126/science.286.5449.2468
- Google Scholar
1. Crismani W
2. Girard C
3. Froger N
4. Pradillo M
5. Santos JL
6. Chelysheva L
7. Copenhaver GP
8. Horlow C
9. Mercier R
(2012) FANCM limits meiotic crossovers
Science 336:1588–1590.

https://doi.org/10.1126/science.1220381
- Google Scholar
1. Cutter AD
2. Payseur BA
(2013) Genomic signatures of selection at linked sites: unifying the disparity among species
Nature Reviews Genetics 14:262–274.

https://doi.org/10.1038/nrg3425
- Google Scholar
1. Dooner HK
(1986)
Genetic fine structure of the BRONZE locus in maize

Genetics 113:1021–1036.
- Google Scholar
1. Dooner HK
2. He L
(2008) Maize genome structure variation: interplay between retrotransposon polymorphisms and genic recombination
The Plant Cell 20:249–258.

https://doi.org/10.1105/tpc.107.057596
- Google Scholar
1. Drouaud J
2. Khademian H
3. Giraut L
4. Zanni V
5. Bellalou S
6. Henderson IR
7. Falque M
8. Mézard C
(2013) Contrasted patterns of crossover and non-crossover at Arabidopsis thaliana meiotic recombination hotspots
PLOS Genetics 9:e1003922.

https://doi.org/10.1371/journal.pgen.1003922
- Google Scholar
1. Duret L
2. Galtier N
(2009) Biased gene conversion and the evolution of mammalian genomic landscapes
Annual Review of Genomics and Human Genetics 10:285–311.

https://doi.org/10.1146/annurev-genom-082908-150001
- Google Scholar
(2007) Using crossover breakpoints in recombinant inbred lines to identify quantitative trait loci controlling the global recombination frequency
Genetics 177:1851–1858.

https://doi.org/10.1534/genetics.107.080622
- Google Scholar
1. Ferdous M
2. Higgins JD
3. Osman K
4. Lambing C
5. Roitinger E
6. Mechtler K
7. Armstrong SJ
8. Perry R
9. Pradillo M
10. Cuñado N
11. Franklin FC
(2012) Inter-homolog crossing-over and synapsis in Arabidopsis meiosis are dependent on the chromosome axis protein AtASY3
PLOS Genetics 8:e1002507.

https://doi.org/10.1371/journal.pgen.1002507
- Google Scholar
(2011) Variation in human recombination rates and its genetic determinants
PLOS ONE 6:e20321.

https://doi.org/10.1371/journal.pone.0020321
- Google Scholar
(2006) Separation of Arabidopsis pollen tetrads is regulated by QUARTET1, a pectin methylesterase gene
Plant Physiology 142:1004–1013.

https://doi.org/10.1104/pp.106.085274
- Google Scholar
1. Fu H
2. Park W
3. Yan X
4. Zheng Z
5. Shen B
6. Dooner HK
(2001) The highly recombinogenic bz locus lies in an unusually gene-rich region of the maize genome
Proceedings of the National Academy of Sciences of USA 98:8903–8908.

https://doi.org/10.1073/pnas.141221898
- Google Scholar
1. Gan X
2. Stegle O
3. Behr J
4. Steffen JG
5. Drewe P
6. Hildebrand KL
7. Lyngsoe R
8. Schultheiss SJ
9. Osborne EJ
10. Sreedharan VT
11. Kahles A
12. Bohnert R
13. Jean G
14. Derwent P
15. Kersey P
16. Belfield EJ
17. Harberd NP
18. Kemen E
19. Toomajian C
20. Kover PX
21. Clark RM
22. Rätsch G
23. Mott R
(2011) Multiple reference genomes and transcriptomes for Arabidopsis thaliana
Nature 477:419–423.

https://doi.org/10.1038/nature10414
- Google Scholar
1. Garcia V
2. Gray S
3. Allison RM
4. Cooper TJ
5. Neale MJ
(2015) Tel1(ATM)-mediated interference suppresses clustered meiotic double-strand-break formation
Nature 520:114–118.

https://doi.org/10.1038/nature13993
- Google Scholar
1. Girard C
2. Crismani W
3. Froger N
4. Mazel J
5. Lemhemdi A
6. Horlow C
7. Mercier R
(2014) FANCM-associated proteins MHF1 and MHF2, but not the other Fanconi anemia factors, limit meiotic crossovers
Nucleic Acids Research 42:9087–9095.

https://doi.org/10.1093/nar/gku614
- Google Scholar
1. Giraut L
2. Falque M
3. Drouaud J
4. Pereira L
5. Martin OC
6. Mézard C
(2011) Genome-wide crossover distribution in Arabidopsis thaliana meiosis reveals sex-specific patterns along chromosomes
PLOS Genetics 7:e1002354.

https://doi.org/10.1371/journal.pgen.1002354
- Google Scholar
(2014) GC content evolution in coding regions of angiosperm genomes: a unifying hypothesis
Trends in Genetics 30:263–270.

https://doi.org/10.1016/j.tig.2014.05.002
- Google Scholar
1. Gore MA
2. Chia J-M
3. Elshire RJ
4. Sun Q
5. Ersoz ES
6. Hurwitz BL
7. Peiffer JA
8. McMullen MD
9. Grills GS
10. Ross-Ibarra J
11. Ware DH
12. Buckler ES
(2009) A first-generation haplotype map of maize
Science 326:1115–1117.

https://doi.org/10.1126/science.1177837
- Google Scholar
(2003) A neutral explanation for the correlation of diversity with recombination rates in humans
American Journal of Human Genetics 72:1527–1535.

https://doi.org/10.1086/375657
- Google Scholar
1. Henderson IR
(2012) Control of meiotic recombination frequency in plant genomes
Current Opinion in Plant Biology 15:556–561.

https://doi.org/10.1016/j.pbi.2012.09.002
- Google Scholar
(2004) The Arabidopsis MutS homolog AtMSH4 functions at an early step in recombination: evidence for two classes of recombination in Arabidopsis
Genes & Development 18:2557–2570.

https://doi.org/10.1101/gad.317504
- Google Scholar
(2008b) Expression and functional analysis of AtMUS81 in Arabidopsis meiosis reveals a role in the second pathway of crossing-over
The Plant Journal 54:152–162.

https://doi.org/10.1111/j.1365-313X.2008.03403.x
- Google Scholar
1. Higgins JD
2. Vignard J
3. Mercier R
4. Pugh AG
5. Franklin FCH
6. Jones GH
(2008a) AtMSH5 partners AtMSH4 in the class I meiotic crossover pathway in Arabidopsis thaliana, but is not required for synapsis
The Plant Journal: for Cell and Molecular Biology 55:28–39.

https://doi.org/10.1111/j.1365-313X.2008.03470.x
- Google Scholar
1. Hill WG
2. Robertson A
(1966) The effect of linkage on limits to artificial selection
Genetics Research 8:269–294.

https://doi.org/10.1017/S0016672300010156
- Google Scholar
1. Hillers KJ
2. Villeneuve AM
(2003) Chromosome-wide control of meiotic crossing over in C. elegans
Current Biology 13:1641–1647.

https://doi.org/10.1016/j.cub.2003.08.026
- Google Scholar
1. Horton MW
2. Hancock AM
3. Huang YS
4. Toomajian C
5. Atwell S
6. Auton A
7. Muliyati NW
8. Platt A
9. Sperone FG
10. Vilhjálmsson BJ
11. Nordborg M
12. Borevitz JO
13. Bergelson J
(2012) Genome-wide patterns of genetic variation in worldwide Arabidopsis thaliana accessions from the RegMap panel
Nature Genetics 44:212–216.

https://doi.org/10.1038/ng.1042
- Google Scholar
1. Hudson RR
2. Kaplan NL
(1995)
Deleterious background selection with recombination

Genetics 141:1605–1617.
- Google Scholar
1. Hunter N
2. Kleckner N
(2001) The single-end invasion: an asymmetric intermediate at the double-strand break to double-holliday junction transition of meiotic recombination
Cell 106:59–70.

https://doi.org/10.1016/S0092-8674(01)00430-5
- Google Scholar
1. Ito H
2. Miura A
3. Takashima K
4. Kakutani T
(2007) Ecotype-specific and chromosome-specific expansion of variant centromeric satellites in Arabidopsis thaliana
Molecular Genetics and Genomics 277:23–30.

https://doi.org/10.1007/s00438-006-0172-2
- Google Scholar
(2012) The chiasmatype theory. A new interpretation of the maturation divisions. 1909
Genetics 191:319–346.

https://doi.org/10.1534/genetics.112.139725
- Google Scholar
1. Jeffreys AJ
2. Neumann R
(2005) Factors influencing recombination frequency and distribution in a human meiotic crossover hotspot
Human Molecular Genetics 14:2277–2287.

https://doi.org/10.1093/hmg/ddi232
- Google Scholar
(1997) Meiosis-specific DNA double-strand breaks are catalyzed by Spo11, a member of a widely conserved protein family
Cell 88:375–384.

https://doi.org/10.1016/S0092-8674(00)81876-0
- Google Scholar
1. Keeney S
2. Neale MJ
(2006) Initiation of meiotic recombination by formation of DNA double-strand breaks: mechanism and regulation
Biochemical Society Transactions 34:523–525.

https://doi.org/10.1042/BST0340523
- Google Scholar
1. Kleckner N
2. Zickler D
3. Jones GH
4. Dekker J
5. Padmore R
6. Henle J
7. Hutchinson J
(2004) A mechanical basis for chromosome function
Proceedings of the National Academy of Sciences of USA 101:12592–12597.

https://doi.org/10.1073/pnas.0402724101
- Google Scholar
1. Knoll A
2. Higgins JD
3. Seeliger K
4. Reha SJ
5. Dangel NJ
6. Bauknecht M
7. Schröpfer S
8. Franklin FCH
9. Puchta H
(2012) The Fanconi anemia ortholog FANCM ensures ordered homologous recombination in both somatic and meiotic cells in Arabidopsis
The Plant Cell 24:1448–1464.

https://doi.org/10.1105/tpc.112.096644
- Google Scholar
(2013) Common and low-frequency variants associated with genome-wide recombination rate
Nature Genetics 46:11–16.

https://doi.org/10.1038/ng.2833
- Google Scholar
(2012) The recombinases DMC1 and RAD51 are functionally and spatially separated during meiosis in Arabidopsis
The Plant Cell 24:2058–2070.

https://doi.org/10.1105/tpc.112.098459
- Google Scholar
(2000) The crystal structure of DNA mismatch repair protein MutS binding to a G x T mismatch
Nature 407:711–717.

https://doi.org/10.1038/35037523
- Google Scholar
1. Lange J
2. Pan J
3. Cole F
4. Thelen MP
5. Jasin M
6. Keeney S
(2011) ATM controls meiotic double-strand-break formation
Nature 479:237–240.

https://doi.org/10.1038/nature10508
- Google Scholar
(2013) Meiotic chromosome structures constrain and respond to designation of crossover sites
Nature 502:703–706.

https://doi.org/10.1038/nature12577
- Google Scholar
1. Long Q
2. Rabanal FA
3. Meng D
4. Huber CD
5. Farlow A
6. Platzer A
7. Zhang Q
8. Vilhjálmsson BJ
9. Korte A
10. Nizhynska V
11. Voronin V
12. Korte P
13. Sedman L
14. Mandáková T
15. Lysak MA
16. Seren Ü
17. Hellmann I
18. Nordborg M
(2013) Massive genomic variation and strong selection in Arabidopsis thaliana lines from Sweden
Nature Genetics 45:884–890.

https://doi.org/10.1038/ng.2678
- Google Scholar
1. López E
2. Pradillo M
3. Oliver C
4. Romero C
5. Cuñado N
6. Santos JL
(2012) Looking for natural variation in chiasma frequency in Arabidopsis thaliana
Journal of Experimental Botany 63:887–894.

https://doi.org/10.1093/jxb/err319
- Google Scholar
(2002) Bay-0 x Shahdara recombinant inbred line population: a powerful tool for the genetic dissection of complex traits in Arabidopsis
TAG 104:1173–1184.

https://doi.org/10.1007/s00122-001-0825-9
- Google Scholar
1. Lu P
2. Han X
3. Qi J
4. Yang J
5. Wijeratne AJ
6. Li T
7. Ma H
(2012) Analysis of Arabidopsis genome-wide variations before and after meiosis and meiotic recombination by resequencing Landsberg erecta and all four products of a single meiosis
Genome Research 22:508–518.

https://doi.org/10.1101/gr.127522.111
- Google Scholar
1. Macaisne N
2. Novatchkova M
3. Peirera L
4. Vezon D
5. Jolivet S
6. Froger N
7. Chelysheva L
8. Grelon M
9. Mercier R
(2008) SHOC1, an XPF endonuclease-related protein, is essential for the formation of class I meiotic crossovers
Current Biology 18:1432–1437.

https://doi.org/10.1016/j.cub.2008.08.041
- Google Scholar
1. Martini E
2. Diaz RL
3. Hunter N
4. Keeney S
(2006) Crossover homeostasis in yeast meiosis
Cell 126:285–295.

https://doi.org/10.1016/j.cell.2006.05.044
- Google Scholar
1. McMullen MD
2. Kresovich S
3. Villeda HS
4. Bradbury P
5. Li H
6. Sun Q
7. Flint-Garcia S
8. Thornsberry J
9. Acharya C
10. Bottoms C
11. Brown P
12. Browne C
13. Eller M
14. Guill K
15. Harjes C
16. Kroon D
17. Lepak N
18. Mitchell SE
19. Peterson B
20. Pressoir G
21. Romero S
22. Oropeza Rosas M
23. Salvo S
24. Yates H
25. Hanson M
26. Jones E
27. Smith S
28. Glaubitz JC
29. Goodman M
30. Ware D
31. Holland JB
32. Buckler ES
(2009) Genetic properties of the maize nested association mapping population
Science 325:737–740.

https://doi.org/10.1126/science.1174320
- Google Scholar
(2007) Synthesis-dependent strand annealing in meiosis
PLOS Biology 5:e299.

https://doi.org/10.1371/journal.pbio.0050299
- Google Scholar
1. Melamed-Bessudo C
2. Levy AA
(2012) Deficiency in DNA methylation increases meiotic crossover rates in euchromatic but not in heterochromatic regions in Arabidopsis
Proceedings of the National Academy of Sciences of USA 109:E981–E988.

https://doi.org/10.1073/pnas.1120742109
- Google Scholar
(2005) A new seed-based assay for meiotic recombination in Arabidopsis thaliana
The Plant Journal 43:458–466.

https://doi.org/10.1111/j.1365-313X.2005.02466.x
- Google Scholar
1. Mercier R
2. Jolivet S
3. Vezon D
4. Huppe E
5. Chelysheva L
6. Giovanni M
7. Nogué F
8. Doutriaux MP
9. Horlow C
10. Grelon M
11. Mézard C
(2005) Two meiotic crossover classes cohabit in Arabidopsis: one is dependent on MER3, whereas the other one is not
Current Biology 15:692–701.

https://doi.org/10.1016/j.cub.2005.02.056
- Google Scholar
(2014) The molecular biology of meiosis in plants
Annual Review of Plant Biology.

https://doi.org/10.1146/annurev-arplant-050213-035923
- Google Scholar
(2012) Loss of DNA methylation affects the recombination landscape in Arabidopsis
Proceedings of the National Academy of Sciences of USA 109:5880–5885.

https://doi.org/10.1073/pnas.1120841109
- Google Scholar
1. Myers S
2. Bowden R
3. Tumian A
4. Bontrop RE
5. Freeman C
6. MacFie TS
7. McVean G
8. Donnelly P
(2010) Drive against hotspot motifs in primates implicates the PRDM9 gene in meiotic recombination
Science 327:876–879.

https://doi.org/10.1126/science.1182363
- Google Scholar
(1996) The effect of recombination on background selection
Genetics Research 67:159–174.

https://doi.org/10.1017/S0016672300033619
- Google Scholar
1. Obmolova G
2. Ban C
3. Hsieh P
4. Yang W
(2000) Crystal structures of mismatch repair protein MutS and its complex with a substrate DNA
Nature 407:703–710.

https://doi.org/10.1038/35037509
- Google Scholar
1. Paape T
2. Zhou P
3. Branca A
4. Briskine R
5. Young N
6. Tiffin P
(2012) Fine-scale population recombination rates, hotspots, and correlates of recombination in the Medicago truncatula genome
Genome Biology and Evolution 4:726–737.

https://doi.org/10.1093/gbe/evs046
- Google Scholar
1. Page SL
2. Hawley RS
(2003) Chromosome choreography: the meiotic ballet
Science 301:785–789.

https://doi.org/10.1126/science.1086605
- Google Scholar
(2010) Prdm9 controls activation of mammalian recombination hotspots
Science 327:835.

https://doi.org/10.1126/science.1181495
- Google Scholar
1. Qi J
2. Chen Y
3. Copenhaver GP
4. Ma H
(2014) Detection of genomic variations and DNA polymorphisms and impact on analysis of meiotic recombination and genetic mapping
Proceedings of the National Academy of Sciences of USA 111:10007–10012.

https://doi.org/10.1073/pnas.1321897111
- Google Scholar
1. Robine N
2. Uematsu N
3. Amiot F
4. Gidrol X
5. Barillot E
6. Nicolas A
7. Borde V
(2007) Genome-wide redistribution of meiotic double-strand breaks in Saccharomyces cerevisiae
Molecular and Cellular Biology 27:1868–1880.

https://doi.org/10.1128/MCB.02063-06
- Google Scholar
1. Salomé PA
2. Bomblies K
3. Fitz J
4. Laitinen RAE
5. Warthmann N
6. Yant L
7. Weigel D
(2012) The recombination landscape in Arabidopsis thaliana F2 populations
Heredity 108:447–455.

https://doi.org/10.1038/hdy.2011.95
- Google Scholar
(2002)
Variation in chiasma frequency among eight accessions of Arabidopsis thaliana

Genetics 162:1415–1422.
- Google Scholar
1. Sandor C
2. Li W
3. Coppieters W
4. Druet T
5. Charlier C
6. Georges M
(2012) Genetic variants in REC8, RNF212, and PRDM9 influence male recombination in cattle
PLOS Genetics 8:e1002854.

https://doi.org/10.1371/journal.pgen.1002854
- Google Scholar
1. Schwacha A
2. Kleckner N
(1995) Identification of double Holliday junctions as intermediates in meiotic recombination
Cell 83:783–791.

https://doi.org/10.1016/0092-8674(95)90191-4
- Google Scholar
(2014) Supergenes and complex phenotypes
Current Biology 24:R288–R294.

https://doi.org/10.1016/j.cub.2014.01.056
- Google Scholar
(2013) Differential association of the conserved SUMO ligase Zip3 with meiotic double-strand break sites reveals regional variations in the outcome of meiotic recombination
PLOS Genetics 9:e1003416.

https://doi.org/10.1371/journal.pgen.1003416
- Google Scholar
(1992) Rad51 protein involved in repair and recombination in S. cerevisiae is a RecA-like protein
Cell 69:457–470.

https://doi.org/10.1016/0092-8674(92)90447-K
- Google Scholar
1. Simon M
2. Loudet O
3. Durand S
4. Bérard A
5. Brunel D
6. Sennesal FX
7. Durand-Tardif M
8. Pelletier G
9. Camilleri C
(2008) Quantitative trait loci mapping in five new large recombinant inbred line populations of Arabidopsis thaliana genotyped with consensus single-nucleotide polymorphism markers
Genetics 178:2253–2264.

https://doi.org/10.1534/genetics.107.083899
- Google Scholar
1. Smith JM
2. Haigh J
(2007) The hitch-hiking effect of a favourable gene
Genetical Research 89:391–403.

https://doi.org/10.1017/S0016672308009579
- Google Scholar
1. Smukowski CS
2. Noor MA
(2011) Recombination rate variation in closely related species
Heredity 107:496–508.

https://doi.org/10.1038/hdy.2011.44
- Google Scholar
1. Spencer CC
2. Deloukas P
3. Hunt S
4. Mullikin J
5. Myers S
6. Silverman B
7. Donnelly P
8. Bentley D
9. McVean G
(2006) The influence of recombination on human genetic diversity
PLOS Genetics 2:e148.

https://doi.org/10.1371/journal.pgen.0020148
- Google Scholar
1. Stahl FW
2. Foss HM
3. Young LS
4. Borts RH
5. Abdullah MF
6. Copenhaver GP
(2004) Does crossover interference count in Saccharomyces cerevisiae?
Genetics 168:35–48.

https://doi.org/10.1534/genetics.104.027789
- Google Scholar
1. Sun Y
2. Ambrose JH
3. Haughey BS
4. Webster TD
5. Pierrie SN
6. Muñoz DF
7. Wellman EC
8. Cherian S
9. Lewis SM
10. Berchowitz LE
11. Copenhaver GP
(2012) Deep genome-wide measurement of meiotic gene conversion using tetrad analysis in Arabidopsis thaliana
PLOS Genetics 8:e1002968.

https://doi.org/10.1371/journal.pgen.1002968
- Google Scholar
(1983) The double-strand-break repair model for recombination
Cell 33:25–35.

https://doi.org/10.1016/0092-8674(83)90331-8
- Google Scholar
1. Thacker D
2. Mohibullah N
3. Zhu X
4. Keeney S
(2014) Homologue engagement controls meiotic DNA break number and distribution
Nature 510:241–246.

https://doi.org/10.1038/nature13120
- Google Scholar
1. Thompson MJ
2. Jiggins CD
(2014) Supergenes and their role in evolution
Heredity 113:1–8.

https://doi.org/10.1038/hdy.2014.20
- Google Scholar
1. Villeneuve AM
2. Hillers KJ
(2001) Whence meiosis?
Cell 106:647–650.

https://doi.org/10.1016/S0092-8674(01)00500-1
- Google Scholar
1. Webster MT
2. Hurst LD
(2012) Direct and indirect consequences of meiotic recombination: implications for genome evolution
Trends in Genetics 28:101–109.

https://doi.org/10.1016/j.tig.2011.11.002
- Google Scholar
1. Wiehe TH
2. Stephan W
(1993)
Analysis of a genetic hitchhiking model, and its application to DNA polymorphism data from Drosophila melanogaster

Molecular Biology and Evolution 10:842–854.
- Google Scholar
1. Wijnker E
2. Velikkakam James G
3. Ding J
4. Becker F
5. Klasen JR
6. Rawat V
7. Rowan BA
8. de Jong DF
9. de Snoo CB
10. Zapata L
11. Huettel B
12. de Jong H
13. Ossowski S
14. Weigel D
15. Koornneef M
16. Keurentjes JJ
17. Schneeberger K
(2013) The genomic landscape of meiotic crossovers and gene conversions in Arabidopsis thaliana
eLife 2:e01426.

https://doi.org/10.7554/eLife.01426
- Google Scholar
(2006) Effects of trans-acting genetic modifiers on meiotic recombination across the a1-sh2 interval of maize
Genetics 174:101–112.

https://doi.org/10.1534/genetics.105.049270
- Google Scholar
1. Yang S
2. Yuan Y
3. Wang L
4. Li J
5. Wang W
6. Liu H
7. Chen JQ
8. Hurst LD
9. Tian D
(2012) Great majority of recombination events in Arabidopsis are gene conversion events
Proceedings of the National Academy of Sciences of USA 109:20992–20997.

https://doi.org/10.1073/pnas.1211827110
- Google Scholar
1. Yao H
2. Schnable PS
(2005) Cis-effects on meiotic recombination across distinct a1-sh2 intervals in a common Zea genetic background
Genetics 170:1929–1944.

https://doi.org/10.1534/genetics.104.034454
- Google Scholar
1. Yelina NE
2. Choi K
3. Chelysheva L
4. Macaulay M
5. de Snoo B
6. Wijnker E
7. Miller N
8. Drouaud J
9. Grelon M
10. Copenhaver GP
11. Mezard C
12. Kelly KA
13. Henderson IR
(2012) Epigenetic remodeling of meiotic crossover frequency in Arabidopsis thaliana DNA methyltransferase mutants
PLOS Genetics 8:e1002844.

https://doi.org/10.1371/journal.pgen.1002844
- Google Scholar
1. Yelina NE
2. Ziolkowski PA
3. Miller N
4. Zhao X
5. Kelly KA
6. Muñoz DF
7. Mann DJ
8. Copenhaver GP
9. Henderson IR
(2013) High-throughput analysis of meiotic crossover frequency and interference via flow cytometry of fluorescent pollen in Arabidopsis thaliana
Nature Protocols 8:2119–2134.

https://doi.org/10.1038/nprot.2013.131
- Google Scholar
(2011) Meiotic double-strand breaks occur once per pair of (sister) chromatids and, via Mec1/ATR and Tel1/ATM, once per quartet of chromatids
Proceedings of the National Academy of Sciences of USA 108:20036–20041.

https://doi.org/10.1073/pnas.1117937108
- Google Scholar
(2014a) Crossover patterning by the beam-film model: analysis and implications
PLOS Genetics 10:e1004042.

https://doi.org/10.1371/journal.pgen.1004042
- Google Scholar
1. Zhang L
2. Wang S
3. Yin S
4. Hong S
5. Kim KP
6. Kleckner N
(2014b) Topoisomerase II mediates meiotic crossover interference
Nature 511:551–556.

https://doi.org/10.1038/nature13442
- Google Scholar

Article and author information

Author details

Piotr A Ziolkowski
1. Department of Plant Sciences, University of Cambridge, Cambridge, United Kingdom
2. Adam Mickiewicz University, Poznań, Poland
Contribution
PAZ, Designed and performed all experiments and analysed data, except Col x Ler tetrad and chiasmata counting, and participated in preparation of the manuscript

Competing interests
The authors declare that no competing interests exist.
Luke E Berchowitz
1. Department of Biology and the Carolina Center for Genome Sciences, University of North Carolina at Chapel Hill, Chapel Hill, United States
2. Lineberger Comprehensive Cancer Center, University of North Carolina School of Medicine, Chapel Hill, United States
Contribution
LEB, Designed, performed and analysed Col x Ler tetrad experiments, read and approved the final manuscript

Competing interests
The authors declare that no competing interests exist.
Christophe Lambing

Department of Plant Sciences, University of Cambridge, Cambridge, United Kingdom

Contribution
CL, Performed the chiasmata experiments, read and approved the final manuscript

Competing interests
The authors declare that no competing interests exist.
Nataliya E Yelina

Department of Plant Sciences, University of Cambridge, Cambridge, United Kingdom

Contribution
NEY, Participated in data acquisition and analysis of CEN3 F₁ data, read and approved the final manuscript

Competing interests
The authors declare that no competing interests exist.
Xiaohui Zhao

Department of Plant Sciences, University of Cambridge, Cambridge, United Kingdom

Contribution
XZ, Performed statistical analyses, read and approved the final manuscript

Competing interests
The authors declare that no competing interests exist.
Krystyna A Kelly

Department of Plant Sciences, University of Cambridge, Cambridge, United Kingdom

Contribution
KAK, Performed statistical analyses, read and approved the final manuscript

Competing interests
The authors declare that no competing interests exist.
Kyuha Choi

Department of Plant Sciences, University of Cambridge, Cambridge, United Kingdom

Contribution
KC, Contributed to data analysis, read and approved the final manuscript

Competing interests
The authors declare that no competing interests exist.
Liliana Ziolkowska

Department of Plant Sciences, University of Cambridge, Cambridge, United Kingdom

Contribution
LZ, Aided PAZ in data acquisition, read and approved the final manuscript

Competing interests
The authors declare that no competing interests exist.
Viviana June

Department of Plant Sciences, University of Cambridge, Cambridge, United Kingdom

Contribution
VJ, Aided PAZ in data acquisition, read and approved the final manuscript

Competing interests
The authors declare that no competing interests exist.
Eugenio Sanchez-Moran

School of Biosciences, University of Birmingham, Birmingham, United Kingdom

Contribution
ES-M, Performed the chiasmata experiments, read and approved the final manuscript

Competing interests
The authors declare that no competing interests exist.
Chris Franklin

School of Biosciences, University of Birmingham, Birmingham, United Kingdom

Contribution
CF, Performed the chiasmata experiments, read and approved the final manuscript

Competing interests
The authors declare that no competing interests exist.
Gregory P Copenhaver
1. Department of Biology and the Carolina Center for Genome Sciences, University of North Carolina at Chapel Hill, Chapel Hill, United States
2. Lineberger Comprehensive Cancer Center, University of North Carolina School of Medicine, Chapel Hill, United States
Contribution
GPC, Designed, performed and analysed Col x Ler tetrad experiments, read and approved the final manuscript

Competing interests
The authors declare that no competing interests exist.

"This ORCID iD identifies the author of this article:" 0000-0002-7962-3862
Ian R Henderson

Department of Plant Sciences, University of Cambridge, Cambridge, United Kingdom

Contribution
IRH, Designed experiments, collected F₁ data, analysed data and wrote the manuscript

For correspondence
irh25@cam.ac.uk

Competing interests
The authors declare that no competing interests exist.

Funding

Royal Society (University Research Fellowship)

Ian R Henderson

Gatsby Charitable Foundation (2962)

Ian R Henderson

National Science Foundation (NSF) (MCB-1121563)

Gregory P Copenhaver

Ministry of Science and Higher Education, Republic of Poland (605/MOB/2011/0)

Piotr A Ziolkowski

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

Research was supported by a Royal Society University Research Fellowship and Gatsby Charitable Foundation Grant 2962 to IRH, and United States National Science Foundation grant MCB-1121563 to GPC. We thank Raphaël Mercier for providing fancm and zip4 mutations and genotyping information and Avi Levy for the 420 line. PAZ was supported by Polish Mobility Plus Fellowship 605/MOB/2011/0. We thank the editor and reviewers for insightful comments.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.