http://orcid.org/0000-0002-7163-9640
Figures
Abstract
For most eukaryotes, sexual reproduction is a fundamental process that requires meiosis. In turn, meiosis typically depends on a reciprocal exchange of DNA between each pair of homologous chromosomes, known as a crossover (CO), to ensure proper chromosome segregation. The frequency and distribution of COs are regulated by intrinsic and extrinsic environmental factors, but much more is known about the molecular mechanisms governing the former compared to the latter. Here we show that elevated temperature induces meiotic hyper-recombination in Arabidopsis thaliana and we use genetic analysis with mutants in different recombination pathways to demonstrate that the extra COs are derived from the major Type I interference sensitive pathway. We also show that heat-induced COs are not the result of an increase in DNA double-strand breaks and that the hyper-recombinant phenotype is likely specific to thermal stress rather than a more generalized stress response. Taken together, these findings provide initial mechanistic insight into how environmental cues modulate plant meiotic recombination and may also offer practical applications.
Author summary
Meiosis is the cell division used by sexually reproducing species to produce sperm and egg cells. During meiosis, programmed Double Strand Breaks (DSBs) occur on each chromosome, which allows DNA to be exchanged between chromosome pairs, resulting in crossovers (COs). COs are necessary to ensure faithful chromosome segregation during meiosis, and thus fertility, but are also an important source of genetic variation. As such, CO formation is tightly regulated. Despite this, CO frequency can be altered by external factors, such as temperature. In Arabidopsis thaliana, COs are formed through two pathways: interference-sensitive (Type I) and interference-insensitive (Type II). An increase in temperature results in an increase in CO frequency. Using a pollen based assay, we show that COs are formed in the Type I pathway, which accounts for approximately 85% of the COs in Arabidopsis. To investigate whether temperature-dependent COs are the result of additional DSBs, we used immunological staining to examine protein foci, which mark the sites of DSBs. We discovered that temperature likely increases CO frequency by shifting alternative repair outcomes, called non-crossovers, to favor additional COs, rather than by increasing DSBs. Lastly, we found that temperature is not a general stress response, as plants subject to salt stress did not exhibit an increase in CO frequency. Our results may prove valuable in aiding plant breeding by enhancing our ability to rapidly introgress suites of elite traits from wild-plants into their crop relatives, a method that is particularly attractive as it does not require genetic modifications.
Citation: Modliszewski JL, Wang H, Albright AR, Lewis SM, Bennett AR, Huang J, et al. (2018) Elevated temperature increases meiotic crossover frequency via the interfering (Type I) pathway in Arabidopsis thaliana. PLoS Genet 14(5): e1007384. https://doi.org/10.1371/journal.pgen.1007384
Editor: Gregory S. Barsh, Stanford University School of Medicine, UNITED STATES
Received: March 12, 2018; Accepted: April 27, 2018; Published: May 17, 2018
Copyright: © 2018 Modliszewski et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was supported by grants from Rijk Zwaan, the National Natural Science Foundation of China (31370347 and 31570314), and by funds from the University of North Carolina at Chapel Hill and the State Key Laboratory of Genetic Engineering, Fudan University. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Sexually reproducing species use a specialized form of cell division known as meiosis to create haploid gametes from diploid progenitor cells (or a similar genomic reduction in polyploids) [1]. A defining feature of meiosis is the exchange, or crossing-over, of DNA between homologous chromosomes. This exchange results in novel allelic combinations not present in either set of parental chromosomes. In most organisms, crossovers (COs) are also critical for stabilizing homologous chromosome pairing, thus ensuring proper segregation of homologous chromosomes during meiosis. In the absence of COs, chromosomes segregate randomly, resulting in imbalances in chromosome numbers in the gametes, and aneuploidy in progeny. Aneuploidy may affect the fertility of the organism and the viability and fertility of its offspring. Perhaps not surprisingly, the frequency and distribution of COs in the genome are genetically regulated.
CO formation is initiated by the creation of Spo11-induced DNA double-strand breaks (DSBs) [2]. In Arabidopsis meioses, approximately 200 DSBs are formed in each meiocyte, but only about 10 are repaired to form COs [3–6]. The remaining DSBs are repaired as non-crossovers (NCOs), presumably through the Synthesis-Dependent Strand Annealing (SDSA) pathway, or conceivably via sister chromatid repair [7]. Perturbations in DSB frequency do not concomitantly alter CO frequency, indicating that, at least to some extent, CO frequency is under homeostatic control [8]. Despite this, CO number and position can be modulated by external factors such as nutrient availability, exposure to environmental toxins, stress, and temperature [9,10]. How these cues are sensed by meiocytes and how they alter processes such as DSB formation and CO regulation are not known.
The ability of temperature to influence CO numbers was noted only four years after the first genetic map was constructed [11,12]. Many plants, including Hordeum vulgare, Vicia faba, Hyacinthus orientalis, and Arabidopsis thaliana have elevated CO frequencies at moderately higher temperatures [9,10], but the molecular mechanisms that mediate CO frequency changes have not been identified. Possible mechanisms for temperature-induced changes in CO frequency include direct effects of temperature on proteins that execute the steps in meiotic recombination, alteration of chromosome axis or synaptonemal complex structure, modulation of chromatin states, and changes in epigenetic regulations, such as DNA methylation.
In Arabidopsis, COs are formed through the Type I and Type II pathways [13]. The majority (~85%) of COs are derived from the Type I pathway in Arabidopsis and are sensitive to the placement of adjacent crossovers (interference sensitive); Type II COs make up most of the remainder of COs and are not sensitive to the placement of adjacent COs [14,15]. Many of the characterized hyper-recombinant mutants in Arabidopsis, including fancm, figl, top3α, and recq4A/B, operate through a shift in designation of recombination intermediates from NCOs to COs in the Type II pathway, resulting in drastic increases in CO frequency [16–18], much like those seen at 28°C. It is not known if temperature-induced COs in Arabidopsis are formed through the Type I, Type II, or both pathways.
To build an understanding of the molecular mechanisms that govern thermal control of CO frequency in plants, we employed a pollen-based visual assay for recombination in Arabidopsis thaliana to demonstrate that thermal stress-induced COs are generated by the interference sensitive Type I pathway and that they act additively with Type II pathway perturbations resulting in an enhanced hyper-recombination phenotype. We also show that this response is temperature-specific rather than a general stress response and that the extra COs occur without increasing the number of meiotic double-strand breaks. We would also like to note that in the process of submitting our manuscript, Lloyd et al., have published an early online manuscript that also shows in increase in Type I COs at elevated temperature in Arabidopsis [19].
Results
Temperature-dependent COs are derived from the Type I pathway
We had previously used a pollen-based, fluorescent tagged line (FTL) system to demonstrate that Arabidopsis grown at standard 20°C conditions experiences elevated CO frequencies when shifted to 28°C [20]. In brief, the FTL system employs pairs of transgenes encoding fluorescent proteins at defined genetic intervals. These transgenes are expressed under the pollen-specific, post-meiotic promoter, LAT52 [21]. These markers are deployed in a qrt1-2 background [22], which causes pollen from individual meioses to be shed as tetrads, allowing COs to be assayed visually by tracking the pattern of fluorescent protein expression in the tetrads.
To determine if temperature-dependent COs are derived from the Type I or Type II pathway, we analyzed mutant lines of msh4 (At4g17380) and mus81 (At4g30870), which disable the Type I and Type II pathways, respectively [14,15]. MSH4 is a homolog of the bacterial mis-match repair gene MutS that has lost its MMR function in Arabidopsis [14]. MSH4 instead functions in the Type I meiotic recombination pathway along with other ZMM proteins such as ZIP4, MSH5, MER3, HEI10, SHOC1 and PTD [23,24], where it is thought to act early in meiotic prophase I to stabilize double-Holliday Junctions (dHJs) [25]. MUS81 is an endonuclease thought to play a role in the resolution of single- and double-HJs in Arabidopsis [26,27] that also mediates CO formation in the Type II pathway [15]. WT and mus81 plants show an increase in meiotic CO frequency when grown at 28°C, whereas msh4 plants do not (Fig 1A and 1B, S1 Table). As a separate method of confirming these results, we immunostained Arabidopsis male meiocytes at diakinesis using a MLH1 antibody. MLH1 is a MutL homolog that co-localizes with MLH3 in meiosis at sites that will become Type I COs [28]. As expected, in WT plants, a significant increase in MLH1 foci at 28°C (average = 11.8, n = 41) was observed relative to 20°C (average = 9.7, n = 57) (Fig 1C and 1D, S2 Table). In mus81 plants, a similar trend in MLH1 foci at 28°C (average = 11.6, n = 34) was found relative to 20°C (average = 10, n = 60) (Fig 1C and 1D, S2 Table). We used a modified Alexander’s stain assay to assess pollen viability under our control and experimental conditions and observed a significant reduction (p-value = 3.110 × 10−16) in the number of viable pollen per anther at 28°C (average = 294.5, n = 28) compared to 20°C (average = 665.7, n = 22) (S3 Table). Together, these data indicate that the increase in COs seen at 28°C is driven by the Type I meiotic recombination pathway in Arabidopsis.
(A, B) Genetic distances measured using FTLs in WT, mus81, and msh4 plants at 20˚C and 28˚C in the I3a (A) and I1a (B) interval. Significantly different values between 20°C and 28°C at α = 0.05 are marked with an asterisk (*). (C, D) MLHI foci counting. (C) Pollen mother cells containing chromosomes (DAPI, blue) and MLHI foci (red). (D) Boxplot of MLH1 foci counts, mean shown as white circle; **** indicates p ≤ 0.0001. Scale bars represent 5 μm.
Temperature-induced hyper-recombination does not compromise interference
CO interference occurs when one CO influences, typically reducing, the likelihood of a second nearby CO [29–31]. Since temperature-induced hyper-recombination acts through the Type I meiotic recombination pathway, which is sensitive to interference, rather than the interference insensitive Type II pathway, we asked whether thermal stress influences interference in Arabidopsis. To compare the strength of interference at normal and high temperatures, we utilized three linked FTL markers and measured whether COs between the first and second marker influenced the frequency of COs between the second and third marker. We calculated the genetic distance in the first interval in the presence of a CO in the adjacent interval (Xwi) and without (Xwo). In the case of no change in interference, an increase in both Xwi and Xwo will be seen, whereas an increase (or decrease) in interference will result in a change in Xwo without a corresponding change in Xwi. Our results indicate a small but non-significant trend of increased interference in plants grown at 28°C (S1 Fig, S4 Table), suggesting that there is either no change in interference or that the effect is subtle.
Temperature-dependent COs are additive with fancm anti-CO activity
Previous reports of hyper-recombination mutants have identified genes that antagonize CO formation in the Type II pathway, such as FANCM, FIGL1, TOP3α and RECQ4A/B [16–18]. The helicase FANCM is thought to process recombination intermediates into NCO products. In its absence, MUS81, but not the ZMM pathway, acts on these recombination intermediates and resolves them as COs, leading to an increase in CO frequency. TOP3α, RECQ4A/B, and FIGL1 also limit COs formation in the Type II pathway, although they do so independently of FANCM. To test if temperature-driven modulation of meiotic CO frequency can operate independently through the Type I pathway in absence of one of these anti-CO factors, we measured CO frequency in WT and fancm in a genetic interval on chromosome 3 (I3a). We observed an additive effective of temperature and fancm in lines grown at 28°C (Fig 2, S1 Table). Temperature-driven modulation of meiotic CO frequency through the Type I pathway thus appears to act autonomously of CO formation in the Type II pathway.
Genetic distances measured using FTLs in WT and fancm plants at 20°C and 28°C in the I3a interval. Significantly different values between 20°C and 28°C at α = 0.05 are marked with an asterisk (*).
Hyper-recombination is not a universal stress phenomenon
In addition to temperature, other factors, such as nutrient availability, developmental stage, and chemicals have also been shown to modulate meiotic CO frequency in plants [9]. In order to test if other stresses induce an increase in CO frequency, we assayed CO frequency in WT plants grown under control and NaCl treatment conditions. No difference in CO frequency was observed between plants grown at 0, 100, and 200 mM NaCl (Fig 3A). To confirm that the NaCl treatment induced a stress response, we assayed the expression of the transcription factor BHLH122, an osmotic and drought stress biomarker [32], and the aldo/keto reductase family protein AKR4C9, an osmotic and salinity stress biomarker [33]. Although both 100 mM and 200 mM NaCl treated plants showed physiological signs of salinity stress (i.e., loss of turgor pressure), a significant increase in expression of both BHLH122 and AKR4C9 was only observed in the 200 mM treated plants. (Fig 3B, S5 Table). Plants treated with 100 mM NaCl showed no change in either BHLH122 expression (fold change = 1.2, adjusted p-value = 0.837), or AKR4C9 expression (fold change = 1.6, adjusted p-value = 0.315). Plants treated with 200 mM NaCl showed significant and marked stress-induced increases in expression of both BHLH122 (fold change = 8.2 adjusted p-value = 4.4 × 10−5) and AKR4C9 (fold change = 10.5, p-value = 1.81 × 10−5). Taken together, these results show that in Arabidopsis salt stress does not elicit the same meiotic hyper-recombination phenotype as thermal stress, suggesting that the latter is a specific rather than general stress response.
(A) Genetic distances measured in the I1a interval using FTL lines, showing SE; neither the 100 mM NaCl or 200 mM NaCl values are different from the control treatment. (B) ΔCT values of AKR4C9 and BHLH122 for 0, 100, and 200 mM NaCl treatments, using TUB4 as an endogenous control. Adjusted p-values are from Tukey’s honest significant difference test, and are indicated as follows: ns = p > 0.5, **** indicates p ≤ 0.0001.
Temperature-dependent COs are the result of a shift in CO designation through the Type I pathway, not an increase in DSBs
CO frequency can be augmented by at least two mechanisms: an increase in the frequency of DSBs, which could concomitantly increase both NCOs and COs, or by shifting the ratio of NCO: COs to favor increased COs without a parallel increase in DSBs. To determine if temperature-induced COs are the result of an increase in DSBs rather than a shift in the NCO: CO ratio, we used immunostaining to count γH2AX and RAD51 foci in male meiocytes at zygotene of WT, mus81, and msh4 plants. H2AX is a variant histone present in approximately 10% of nucleosomes that becomes phosphorylated (γH2AX) at serine 139 in response to DNA damage [34]. RAD51 is a homolog of the bacterial RecA protein that forms a nucleoprotein filament during recombination [35]. Both γH2AX and RAD51 foci mark the sites of meiotic DSBs in Arabidopsis [6,36]. In WT plants, the number of γH2AX foci (average = 201.8, n = 26) and RAD51 foci (average = 180.5, n = 24) at 20°C did not differ significantly from those at 28°C (average = 194.8, n = 16; average = 181.2, n = 14 respectively) (Fig 4, S6 Table, S7 Table). In mus81 and msh4 mutants, a trend towards fewer DSBs was observed at 28°C, although the results were neither highly significant nor seen in all cases (Fig 4, RAD51 in mus81, p-value = 0.042 and γH2AX in msh4, p-value = 0.035). In spo11-1-1 mutants, no difference in γH2AX foci was seen in plants grown at 20°C (average = 19.8, n = 37) and 28°C (average = 20,8, n = 28) (Fig 4). As expected, the spo11-1-1 mutants did exhibit a dramatic decrease in CO frequency relative to WT plants. These results lead us to reject the hypothesis that the extra COs observed at elevated temperatures in Arabidopsis are the result of an increase in DSBs.
(A, B) γH2Ax foci counting. (A) Pollen mother cells containing chromosomes (DAPI, blue) and γH2Ax foci (red). (B) Boxplots of γH2Ax foci at 20˚C and 28˚C; mean shown as white circle. (C, D) RAD51 foci counting. (C) Pollen mother cells containing chromosomes (DAPI, blue) and RAD51 foci (red). (D) Boxplots of RAD51 foci at 20°C and 28°C; mean shown as white circle. p-values are indicated as follows: ns = p > 0.5, * = p ≤ 0.05, ** = p ≤ 0.01. Scale bars represent 5 μm.
Discussion
Despite the necessity for COs to ensure proper chromosome segregation and evidence indicating that CO frequency is under homeostatic control, CO numbers can be altered by external factors, such as temperature. Here, we show that in Arabidopsis, these additional COs are formed through the interference-sensitive Type I pathway, and that the number of DSBs remains unchanged at 28˚C. This indicates that the increase in CO frequency is derived from a shift in the ratio of NCOs to COs. Although it is known that other factors can alter CO frequencies, we demonstrate here that in Arabidopsis an increase in CO frequency is not a ubiquitous response to stress.
Temperature-dependent modulation of meiotic CO frequency mimics the hyper-recombination phenotypes of Type II anti-CO factors mutants, such as fancm. As such, one may expect that temperature-dependent COs are derived from either the Type II pathway or both pathways. Surprisingly, we discovered that the temperature dependent COs are derived from the Type I pathway. These results contrast those seen in barley, where temperature-dependent COs appear to be derived from the Type II pathway and the distribution but not the frequency of Type I COs is changed [37]. At elevated temperatures, COs in barley redistributed along the length of the chromosomes, shifting from a terminalized to more medialized pattern. Our analysis did not provide the necessary resolution to detect a similar pattern in Arabidopsis, but chromosome field redistribution of COs has been observed in met1 DNA methyltransferase mutants [38]. These observations raise the possibility that epigenetic mechanisms may play a role in regulating the frequency and distribution of COs in response to external cues.
Most sexually reproducing species generate an order of magnitude more meiotic DSBs compared to COs. In addition, most species have relatively few COs per chromosome, though there are interesting exceptions including several fungi, the SAR (Stramenopiles-Alveolates-Rhizaria Eukaryote) supergroup, and some insects [39–41]. At a molecular level, several factors work to limit the number of COs. FANCM and RECQ4A/B are both helicases that promote NCO formation by disassembling recombination intermediates that have formed D-loops and consequently directing them through the SDSA pathway, where they go on to form NCOs [16,18]. TOP3α works in conjunction with RECQ4A/B to maintain the recombination intermediates as NCOs, while FIGL inhibits homologous strand invasion, a necessary step for meiotic recombination [17]. Thus, it appears that despite an abundant pool of initiating events, and a default program that, if left unrestrained, will produce more COs, there is a common trend to limit COs. This problem may be particularly acute in polyploids, which must reduce CO numbers to avoid chromosome entanglements during segregation [39]. It is therefore intriguing that exposure to modest temperature increases, at least in laboratory settings, can cause plants to overcome these limiting programs and evoke a hyper-recombinant response. It would be interesting to explore whether polyploids block temperature-dependent hyper-recombination in addition to strengthening the limits on CO numbers. Our findings in this study and future elaboration of the molecular mechanisms used to increase CO frequency in response to external cues will enable such hypotheses to be tested.
DNA double-strand breaks occur in somatic cells as a result of temperature, osmotic, and oxidative stress, UV-radiation, and many other factors [42–45]. These breaks are repaired via non-homologous end-joining or somatic homologous recombination, effectively resulting in a relationship where external stressors increase somatic homologous recombination. In Arabidopsis, the increase in meiotic CO frequency in response to temperature is not a universal response to stress than can be replicated by other factors, such as salt stress. It should be noted that our assays do not exclude the possibility that salinity stress may result in a redistribution of COs along the chromosomes. It is also possible that even within Arabidopsis, the temperature response may be heterochiasmic. In barley, temperature only increases CO frequency in male meiosis, not female meiosis [46]. Furthermore, the distribution and frequency of COs differs between males and females in Arabidopsis [47,48]. Our pollen FTL system assays male meiosis only, thus it is possible that female meiosis may respond differently to elevated temperatures. Given the fluidity of the relationship between stress and CO frequency across many organisms, it seems that temperature-specific heterochiasmy should not be unexpected.
We selected 28°C as our thermal stress condition based on a range of elevated temperatures assayed in experiments used to originally detect the CO response [20]. Lloyd et al. observed an equally potent increase in COs at moderately cold (8°C- 13°C) temperatures as well [19]. All of our experiments have used the Columbia ecotype, which is adapted to a temperate climate. It would be interesting to test whether the thermal stress-induced CO phenotypes we observed also occur in ecotypes adapted to tropical climates, such as Cvi from the Cape Verde Islands [49], or cold-adapted ecotypes [50]. It is possible that temperate plants are particularly responsive to thermal stress cues and that plants from more extreme climates would have meiotic programs adapted to those conditions and would not be as responsive. For example, isolates of the fungus Sordaria fimicola from harsh micro-environments have higher recombination frequencies compared to those from mild micro-environments at the same collection site leading to the suggestion that adaptation to the harsh conditions had selected for increased COs [51], though in this case their response to changing conditions was not tested. Alternatively, plants adapted to more extreme climactic conditions may respond equally robustly but at different temperatures or “set points”.
While it is tempting to speculate that the temperature response observed here may be of adaptive value, it is equally possible that temperature dependent modulation of meiotic CO frequency is a function of physical factors in the cell, such as synaptonemal complex (SC) length. In barley, the increase in CO frequency observed in male meiosis is coupled with an increase in SC length, although the causal direction of the relationship is unclear [37]. In Arabidopsis, although it is unknown how SC length changes in response to temperature, SC length is longer in males, which also exhibit higher CO frequency [47,48]. It is also possible that the ability to modulate meiotic CO frequency in response to environmental cues may have initially been purely mechanistic in nature, but that the direction and magnitude of response may have subsequently been subject to selective pressure.
In addition to the potential for adaptive value, modulation of meiotic CO frequency via temperature is also of practical importance. Temperature-dependent modulation of meiotic CO frequency provides a means of aiding plant breeding without editing the genome or the necessity of working in mutant backgrounds; our observed reduction in pollen viability at high temperatures may require that more modest increases in temperature are utilized. Easily elevating CO frequency will reduce the number of F2 progeny necessary for classic genetic mapping of traits and will aid in disrupting persistent linkage blocks so that desirable traits can be isolated and bred into elite lines. Valuable genetic diversity, including disease resistance, remains locked in wild relatives and meiotic recombination limits the ability to introgress those traits into commercial relatives. Temperature treatment during reproductive stages may offer a simple and cost effective means to improve our ability to tap into these natural genetic resources.
Materials and methods
Plant lines
Seeds for mutant lines were obtained from the Arabidopsis Biological Resource Center; the following T-DNA lines were used: MSH4 (At4g17380, SALK_136296), MUS81 (At4g30870, SALK_107515), FANCM (At1g35530, SALK_120621). For SPO11-1 (At3g13170) the spo11-1-1 mutant [52] was used. FTL lines were generated as described previously [20]. Mutant lines are in the Columbia-0 ecotype background (CS60000) with the exception of the spo11-1-1 line, which is derived from the Wassilewskija ecotype. FTL lines are in the Columbia-3 ecotype background (CS8846). DNA was extracted as described previously [3], and T-DNA lines were genotyped via PCR using primers and conditions provided in S8 Table.
Plant growth and treatment conditions
Seeds were sown on Metromix-360 (Sun-Gro). Unless otherwise noted, plants were grown under 18 hour days at 20°C in a growth room for control conditions and at 28°C under 18 hours days in a Percival chamber for heat treatment conditions. For the heat treatment, flowering plants were placed in the 28°C chamber for five days and tetrads were counted on the fifth day. For the salt treatment, plants with dry soil were bottom watered until saturation with 0mM NaCl, 100mM NaCl, and 200mM NaCl. Tetrads were counted five days after treatment.
Pollen tetrad assay
Crossover frequency was assayed via a visual fluorescent pollen transgene assay as described previously [20]. In the two color experiments (I1a and I3a intervals, S9 Table), tetrads were counted manually using a Nikon Eclipse E1000 epifluorescence microscope. In the three color experiment (I5cd interval, S9 Table), pollen grains were counted by first capturing an image of the entire slide for each sample on a Zeiss LSM 880 confocal laser scanning microscope, under the conditions provided in S1 File. Tetrads were then counted manually using Fiji [53] and a custom built program TetradAnalysis (https://github.com/jmodlis/TetradAnalysis). For all intervals, genetic distance was calculated using the Perkins equation [54], standard errors were calculated via Stahl Lab Online Tools (http://elizabethhousworth.com/StahlLabOnlineTools/), and p-values were calculated using R [55]. Interference was calculated following the method of Malkova et al [56]. The genetic distance in the I5c interval (FTL1963 and FTL1143, CFP and YFP respectively) measured in tetrads without a CO in the adjacent I5d interval (FTL1143 and FTL2450, YFP and dsRED2, respectively) was divided by the genetic distance in the I5c interval measured in tetrads with a CO in the adjacent I5d interval [57]. The difference in interference ratios at 20°C and 28°C was tested by calculating the p-value from the Z-score. The Z-score was calculated using the formula Z = |R20 –R28|/√(VarR20 + VarR28); the variance for each ratio was calculated as in van Kempen & van Vliet [58] and Stuart & Ord [59] with the exception that covariance was assumed to be zero.
RT-PCR
RNA was collected from floral buds used in the salt treatment on the same day that tetrads were counted. RNA was extracted using the RNeasy Plant Mini Kit (Qiagen), followed by TURBO DNA-free DNase Treatment (Ambion). cDNA was generated using a ProtoScript II First Strand cDNA Synthesis Kit (New England Biolabs). One μg of cDNA was used in RT-PCR reactions, which were performed using the primers and conditions in S8 Table and the PowerUp SYBR Green Master Mix (Applied BioSystems) following the manufacturer’s instructions on a QuantStudio 6 Flex Real-Time PCR System (Applied Biosytems). Fold change (2-ΔΔCT) values for BHLH122 (At1g51140) and AKR4C9 (At2g37770) were calculated using the comparative CT method of Schmittgen and Livak [60] with TUB4 (At5g44340) as an endogenous control. To test for significant differences in gene expression between control and NaCl treatments for BHLH122 and AKR4C9, ΔCT values were used, where ΔCT is the change in the expression of each gene relative to the endogenous control (e.g., ΔCT(BHLH122)—ΔCT(TUB4)). For each gene, an ANOVA was conducted to test for differences among ΔCT means for all treatments, followed by a post-hoc Tukey’s honest significant difference test.
Cytology
Immunofluorescence was performed as previously described by Wang et al., with minor modifications [39]. Unopened flower buds were fixed after 28°C and 20°C treatment for at least 24 hours in a manual climatic box (17600 Lux, 85% humidity, 16h light and 8h dark). The slides were incubated overnight with primary antibodies diluted 1:200 (RAD51, γH2AX, MLH1) in blocking buffer (goat serum, AR0009, Bosterbio) at 4°C and then at 37°C for 60 min with secondary antibody (1:1000, Goat anti-Rabbit IgG (H+L) Cross-Adsorbed Secondary Antibody, Alexa Fluor 555, catalog # A-21428, Thermo Fisher Scientific). Images were taken using a Zeiss Axio Imager A1 microscope and processed using Adobe Photoshop CS6. RAD51, γH2AX, and MLH1 foci were counted using Image Tool version 3.0 software (University of Texas Health Science Center, San Antonio, USA). The viability of pollen from control and temperature treated plants, grown and harvested as describe above, was assayed using the method described by Peterson et al [51]. Welch’s t-test was used to test for differences in means at 20°C and 28°C in R.
Supporting information
S1 Fig. Interference at 28°C.
Genetic distances measured using the I5cd FTL interval in WT and mus81. Genetic distance measured in the I5c interval both with and without COs in the adjacent I5d interval. CFP, YFP and dsRED transgenes shown in blue, yellow, and red, respectively.
https://doi.org/10.1371/journal.pgen.1007384.s001
(TIF)
S1 Table. Two color pollen count data from I3a and I1a intervals for WT, mus81, msh4, and fancm lines conducted in separate experiments.
https://doi.org/10.1371/journal.pgen.1007384.s002
(XLSX)
S2 Table. MLH1 foci data for WT, msh4, and mus81 plants.
https://doi.org/10.1371/journal.pgen.1007384.s003
(XLSX)
S3 Table. Number of viable pollen grains per anther.
https://doi.org/10.1371/journal.pgen.1007384.s004
(XLSX)
S4 Table. Tetrad count data from the I5cd interval for WT and mus81 plants grown at 20°C and 28°C.
https://doi.org/10.1371/journal.pgen.1007384.s005
(XLSX)
S5 Table. Raw qRTPCR data taken from QuantStudio 6 Flex Real-Time PCR system (Applied Biosystems), ΔCT, ΔΔCT, and 2-ΔΔCT calculations, and ANOVA results.
https://doi.org/10.1371/journal.pgen.1007384.s006
(XLSX)
S6 Table. RAD51 foci counts for WT, msh4, and mus81 plants.
https://doi.org/10.1371/journal.pgen.1007384.s007
(XLSX)
S7 Table. γH2AX foci counts from WT, msh4, mus81, and spo11-1-1 lines.
https://doi.org/10.1371/journal.pgen.1007384.s008
(XLSX)
S8 Table. PCR conditions and primers used in this study.
https://doi.org/10.1371/journal.pgen.1007384.s009
(XLSX)
S9 Table. Description of FTL intervals used in this study.
https://doi.org/10.1371/journal.pgen.1007384.s010
(XLSX)
S1 File. Microscope settings for the Zeiss LSM 880 confocal laser scanning microscope used in three color (I5cd) experiment.
https://doi.org/10.1371/journal.pgen.1007384.s011
(PDF)
Acknowledgments
We thank Prof. Raphaël Mercier at Institut National de la Recherche Agronomique (INRA) (France) for providing MLH1 antibody and Tony Perdue at UNC for help with confocal microscopy and Savannah Nunnery at UNC for help with genotyping.
References
- 1. Mercier R, Mézard C, Jenczewski E, Macaisne N, Grelon M. The Molecular Biology of Meiosis in Plants. Annu Rev Plant Biol. 2015; 66(1):297–327. pmid:25494464
- 2. Keeney S, Giroux CN, Kleckner N. Meiosis-specific DNA double-strand breaks are catalyzed by Spo11, a member of a widely conserved protein family. Cell. 1997; 88(3):375–84. pmid:9039264
- 3. Copenhaver GP, Browne WE, Preuss D. Assaying genome-wide recombination and centromere functions with Arabidopsis tetrads. Proc Natl Acad Sci U S A. 1998; 95 (1) 247–52. pmid:9419361
- 4. Mercier R, Jolivet S, Vezon D, Huppe E, Chelysheva L, Giovanni M, et al. Two meiotic crossover classes cohabit in Arabidopsis: one is dependent on MER3, whereas the other one is not. Curr Biol. 2005; 15(8):692–701. pmid:15854901
- 5. Chelysheva L, Gendrot G, Vezon D, Doutriaux M-P, Mercier R, Grelon M. Zip4/Spo22 is required for class I CO formation but not for synapsis completion in Arabidopsis thaliana. PLoS Genet. 2007; 3(5):e83. pmid:17530928
- 6. Sanchez-Moran E, Santos J-L, Jones GH, Franklin FCH. ASY1 mediates AtDMC1-dependent interhomolog recombination during meiosis in Arabidopsis. Genes Dev. 2007; 21(17):2220–33. pmid:17785529
- 7. McMahill MS, Sham CW, Bishop DK. Synthesis-dependent strand annealing in meiosis. PLoS Biol. 2007; 5(11):e299. pmid:17988174
- 8. Martini E, Diaz RL, Hunter N, Keeney S. Crossover Homeostasis in Yeast Meiosis. Cell. 2006; 126(2):285–95. pmid:16873061
- 9. Modliszewski JL, Copenhaver GP. Meiotic recombination gets stressed out: CO frequency is plastic under pressure. Curr Opin Plant Biol. 2017; 36:95–102. pmid:28258986
- 10. Morgan CH, Zhang H, Bomblies K. Are the effects of elevated temperature on meiotic recombination and thermotolerance linked via the axis and synaptonemal complex? Philos Trans R Soc Lond B Biol Sci. 2017; 372(1736):20160470. pmid:29109229
- 11. Sturtevant AH. The linear arrangement of six sex-linked factors in Drosophila, as shown by their mode of association. J Exp Zool. 1913; 14(1):43–59.
- 12. Plough HH. The effect of temperature on crossingover in Drosophila. J Exp Zool. 1917; 24(2):147–209.
- 13. Copenhaver GP, Housworth EA, Stahl FW. Crossover interference in Arabidopsis. Genetics. 2002; 160(4):1631–9. pmid:11973316
- 14. Higgins JD, Armstrong SJ, Franklin FCH, Jones GH. The Arabidopsis MutS homolog AtMSH4 functions at an early step in recombination: evidence for two classes of recombination in Arabidopsis. Genes Dev. 2004; 18(20):2557–70. pmid:15489296
- 15. Berchowitz LE, Francis KE, Bey AL, Copenhaver GP. The role of AtMUS81 in interference-insensitive crossovers in A. thaliana. PLoS Genet. 2007; 3(8):e132. pmid:17696612
- 16. Crismani W, Girard C, Froger N, Pradillo M, Santos JL, Chelysheva L, et al. FANCM limits meiotic crossovers. Science. 2012; 336(6088):1588–90. pmid:22723424
- 17. Girard C, Chelysheva L, Choinard S, Froger N, Macaisne N, Lehmemdi A, et al. AAA-ATPase FIDGETIN-LIKE 1 and Helicase FANCM Antagonize Meiotic Crossovers by Distinct Mechanisms. PLoS Genet. 2015;11(7):1–22. pmid:26161528
- 18. Séguéla-Arnaud M, Crismani W, Larchevêque C, Mazel J, Froger N, Choinard S, et al. Multiple mechanisms limit meiotic crossovers: TOP3α and two BLM homologs antagonize crossovers in parallel to FANCM. Proc Natl Acad Sci. 2015; 112(15):4713–8. pmid:25825745
- 19. Lloyd A, Morgan C, Franklin C, Bomblies K. Plasticity of Meiotic Recombination Rates in Response to Temperature in Arabidopsis. Genetics. 2018; 208(4):1409–20. pmid:29496746
- 20. Francis KE, Lam SY, Harrison BD, Bey AL, Berchowitz LE, Copenhaver GP. Pollen tetrad-based visual assay for meiotic recombination in Arabidopsis. Proc Natl Acad Sci U S A. 2007; 104(10):3913–8. pmid:17360452
- 21. Twell D, Yamaguchi J, McCormick S. Pollen-specific gene expression in transgenic plants: coordinate regulation of two different tomato gene promoters during microsporogenesis. Development. 1990; 109(3):705–13. pmid:2401221
- 22. Francis KE, Lam SY, Copenhaver GP. Separation of Arabidopsis pollen tetrads is regulated by QUARTET1, a pectin methylesterase gene. Plant Physiol. 2006; 142(3):1004–13. pmid:16980565
- 23. Osman K, Higgins JD, Sanchez-Moran E, Armstrong SJ, Franklin FCH. Pathways to meiotic recombination in Arabidopsis thaliana. New Phytol. 2011; 190(3):523–44. pmid:21366595
- 24. Chelysheva L, Vezon D, Chambon A, Gendrot G, Pereira L, Lemhemdi A, et al. The Arabidopsis HEI10 is a new ZMM protein related to Zip3. PLoS Genet. 2012; 8(7):e1002799. pmid:22844245
- 25. Franklin FCH, Higgins JD, Sanchez-Moran E, Armstrong SJ, Osman KE, Jackson N, et al. Control of meiotic recombination in Arabidopsis: role of the MutL and MutS homologues. Biochem Soc Trans. 2006; 34(4):542–4. pmid:16856855
- 26. Boddy MN, Gaillard PHL, McDonald WH, Shanahan P, Yates JR, Russell P. Mus81-Eme1 are essential components of a Holliday junction resolvase. Cell. 2001; 107(4):537–48. pmid:11719193
- 27. Hartung F, Suer S, Bergmann T, Puchta H. The role of AtMUS81 in DNA repair and its genetic interaction with the helicase AtRecQ4A. Nucleic Acids Res. 2006; 34(16):4438–48. pmid:16945961
- 28. Lhuissier FGP, Offenberg HH, Wittich PE, Vischer NOE, Heyting C. The mismatch repair protein MLH1 marks a subset of strongly interfering crossovers in tomato. Plant Cell. 2007; 19(3):862–76. pmid:17337626
- 29. Sturtevant AH. A Third Group of Linked Genes in Drosophila Ampelophila. Science. 1913; 37(965):990–2. pmid:17833164
- 30. Sturtevant AH. The behavior of the chromosomes as studied through linkage. Z Indukt Abstamm Vererbungsl. 1915; 13(1):234–87.
- 31. Berchowitz LE, Copenhaver GP. Genetic interference: don’t stand so close to me. Curr Genomics. 2010; 11(2):91–102. pmid:20885817
- 32. Liu W, Tai H, Li S, Gao W, Zhao M, Xie C, et al. bHLH122 is important for drought and osmotic stress resistance in Arabidopsis and in the repression of ABA catabolism. New Phytol. 2014; 201(4):1192–204. pmid:24261563
- 33. Sham A, Moustafa K, Al-Ameri S, Al-Azzawi A, Iratni R, AbuQamar S. Identification of Arabidopsis candidate genes in response to biotic and abiotic stresses using comparative microarrays. PLoS One. 2015; 10(5): e0125666. pmid:25933420
- 34. Rogakou EP, Boon C, Redon C, Bonner WM. Megabase Chromatin Domains Involved in DNA Double-Strand Breaks in Vivo. J Cell Biol. 1999; 146(5):905–16. pmid:10477747
- 35. Sung P, Krejci L, Van Komen S, Sehorn MG. Rad51 Recombinase and Recombination Mediators. J Biol Chem. 2003; 278(44):42729–32. pmid:12912992
- 36. Kurzbauer M-T, Uanschou C, Chen D, Schlogelhofer P. The Recombinases DMC1 and RAD51 Are Functionally and Spatially Separated during Meiosis in Arabidopsis. Plant Cell. 2012; 24(5):2058–70. pmid:22589466
- 37. Phillips D, Jenkins G, Macaulay M, Nibau C, Wnetrzak J, Fallding D, et al. The effect of temperature on the male and female recombination landscape of barley. New Phytol. 2015; 208(2):421–9. pmid:26255865
- 38. Yelina NE, Choi K, Chelysheva L, Macaulay M, de Snoo B, Wijnker E, et al. Epigenetic Remodeling of Meiotic Crossover Frequency in Arabidopsis thaliana DNA Methyltransferase Mutants. PLoS Genet. 2012; 8(8):e1002844. pmid:22876192
- 39. Wang Y, Cheng Z, Lu P, Timofejeva L, Ma H. Molecular cell biology of male meiotic chromosomes and isolation of male meiocytes in Arabidopsis thaliana. Methods Mol Biol. 2014; 1110:217–30. pmid:24395259
- 40. Mancera E, Bourgon R, Brozzi A, Huber W, Steinmetz LM. High-resolution mapping of meiotic crossovers and non-crossovers in yeast. Nature. 2008 Jul 24; 454(7203):479–85. pmid:18615017
- 41. Stapley J, Feulner PGD, Johnston SE, Santure AW, Smadja CM. Variation in recombination frequency and distribution across eukaryotes: patterns and processes. Philos Trans R Soc B Biol Sci. 2017; 372(1736):20160455. pmid:29109219
- 42. Kovalchuk I, Kovalchuk O, Kalck V, Boyko V, Filkowski J, Heinlein M, et al. Pathogen-induced systemic plant signal triggers DNA rearrangements. Nature. 2003; 423(6941):760–2. pmid:12802336
- 43. Molinier J, Ries G, Zipfel C, Hohn B. Transgeneration memory of stress in plants. Nature. 2006; 442(7106):1046–9. pmid:16892047
- 44. Pecinka A, Rosa M, Schikora A, Berlinger M, Hirt H, Luschnig C, et al. Transgenerational stress memory is not a general response in Arabidopsis. PLoS One. 2009; 4(4): e5202. pmid:19381297
- 45. Kathiria P, Sidler C, Golubov A, Kalischuk M, Kawchuk LM, Kovalchuk I. Tobacco Mosaic Virus Infection Results in an Increase in Recombination Frequency and Resistance to Viral, Bacterial, and Fungal Pathogens in the Progeny of Infected Tobacco Plants. Plant Physiol. 2010; 153(4):1859–70. pmid:20498336
- 46. Higgins JD, Perry RM, Barakate A, Ramsay L, Waugh R, Halpin C, et al. Spatiotemporal Asymmetry of the Meiotic Program Underlies the Predominantly Distal Distribution of Meiotic Crossovers in Barley. Plant Cell. 2012; 24(10):4096–109. pmid:23104831
- 47. Drouaud J, Mercier R, Chelysheva L, Bérard A, Falque M, Martin O, et al. Sex-specific crossover distributions and variations in interference level along Arabidopsis thaliana chromosome 4. PLoS Genet. 2007; 3(6): e106. pmid:17604455
- 48. Giraut L, Falque M, Drouaud J, Pereira L, Martin OC, Mézard C. Genome-wide crossover distribution in Arabidopsis thaliana meiosis reveals sex-specific patterns along chromosomes. PLoS Genet. 2011; 7(11): e1002354. pmid:22072983
- 49. Alonso-Blanco C, El-Assal SE-D, Coupland G, Koornneef M. Analysis of natural variation at flowering time loci in the Landsberg erecta and Cape Verde Island ecotypes of Arabidopsis thaliana. Genetics. 1998; 149(2):749–64. pmid:9611189
- 50. Oakley CG, Savage L, Lotz S, Larson GR, Thomashow MF, Kramer DM, et al. Genetic basis of photosynthetic responses to cold in two locally adapted populations of Arabidopsis thaliana. J Exp Bot. 2018; 69(3):699–709. pmid:29300935
- 51. Saleem M, Lamb BC, Nevo E. Inherited differences in crossing over and gene conversion frequencies between wild strains of Sordaria fimicola from “Evolution Canyon.” Genetics. 2001; 159(4):1573–93. pmid:11779798
- 52. Grelon M, Vezon D, Gendrot G, Pelletier G. AtSPO11-1 is necessary for efficient meiotic recombination in plants. EMBO J. 2001; 20:589–600. pmid:11157765
- 53. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: An open-source platform for biological-image analysis. Nat Methods. 2012; 9(7):676–82. pmid:22743772
- 54. Perkins DD. Biochemical Mutants in the Smut Fungus Ustilago Maydis. Genetics. 1949; 34(5):607–26. pmid:17247336
- 55. Team RC. R: A Language and Environment for Statistical Computing [Internet]. R Foundation for Statistical Computing. 2013.
- 56. Malkova A, Swanson J, German M, Mccusker JH, Housworth EA, Stahl FW, et al. Gene Conversion and Crossing Over Along the 405-kb Left Arm of Saccharomyces cerevisiae Chromosome VII. Gene. 2004; 63:49–63.
- 57. Berchowitz LE, Copenhaver GP. Fluorescent Arabidopsis tetrads: a visual assay for quickly developing large crossover and crossover interference data sets. Nat Protoc. 2008; 3(1):41–50. pmid:18193020
- 58. van Kempen GM, van Vliet LJ. Mean and variance of ratio estimators used in fluorescence ratio imaging. Cytometry. 2000; 39(4):300–5. pmid:10738283
- 59.
Stuart A, Ord JK. Kendall’s Advanced Theory of Statistics, Volume 1 Distribution Theory. 6th ed. New York: Halsted Press; 1994.
- 60. Schmittgen TD, Livak K. Analyzing real-time PCR data by the comparative CT method. Nature Protocols. 2008; 3:1101–8.