Recombination Rate Heterogeneity within Arabidopsis Disease Resistance Genes

Structural variation and DNA methylation shape the centromere-proximal meiotic crossover landscape in Arabidopsis

BACKGROUND

Centromeres load kinetochore complexes onto chromosomes, which mediate spindle attachment and allow segregation during cell division. Although centromeres perform a conserved cellular function, their underlying DNA sequences are highly divergent within and between species. Despite variability in DNA sequence, centromeres are also universally suppressed for meiotic crossover recombination, across eukaryotes. However, the genetic and epigenetic factors responsible for suppression of centromeric crossovers remain to be completely defined.

RESULTS

To explore the centromere-proximal meiotic recombination landscape, we map 14,397 crossovers against fully assembled Arabidopsis thaliana (A. thaliana) genomes. A. thaliana centromeres comprise megabase satellite repeat arrays that load nucleosomes containing the CENH3 histone variant. Each chromosome contains a structurally polymorphic region of ~3-4 megabases, which lack crossovers and include the satellite arrays. This polymorphic region is flanked by ~1-2 megabase low-recombination zones. These recombination-suppressed regions are enriched for Gypsy/Ty3 retrotransposons, and additionally contain expressed genes with high genetic diversity that initiate meiotic recombination, yet do not crossover. We map crossovers at high-resolution in proximity to CEN3, which resolves punctate centromere-proximal hotspots that overlap gene islands embedded in heterochromatin. Centromeres are densely DNA methylated and the recombination landscape is remodelled in DNA methylation mutants. We observe that the centromeric low-recombining zones decrease and increase crossovers in CG (met1) and non-CG (cmt3) mutants, respectively, whereas the core non-recombining zones remain suppressed.

CONCLUSION

Our work relates the genetic and epigenetic organization of A. thaliana centromeres and flanking pericentromeric heterochromatin to the zones of crossover suppression that surround the CENH3-occupied satellite repeat arrays.

Fine-mapping and evolutionary history of R-BPMV, a dominant resistance gene to Bean pod mottle virus in Phaseolus vulgaris L

KEY MESSAGE

R-BPMV is located within a recently expanded TNL cluster in the Phaseolus genus with suppressed recombination and known for resistance to multiple pathogens including potyviruses controlled by the I gene. Bean pod mottle virus (BPMV) is a comovirus that infects common bean and legumes in general. BPMV is distributed throughout the world and is a major threat on soybean, a closely related species of common bean. In common bean, BAT93 was reported to carry the R-BPMV resistance gene conferring resistance to BPMV and linked with the I resistance gene. To fine map R-BPMV, 182 recombinant inbred lines (RILs) derived from the cross BAT93 × JaloEEP558 were genotyped with polymerase chain reaction (PCR)-based markers developed using genome assemblies from G19833 and BAT93, as well as BAT93 BAC clone sequences. Analysis of RILs carrying key recombination events positioned R-BPMV to a target region containing at least 16 TIR-NB-LRR (TNL) sequences in BAT93. Because the I cluster presents a suppression of recombination and a large number of repeated sequences, none of the 16 TNLs could be excluded as R-BPMV candidate gene. The evolutionary history of the TNLs for the I cluster were reconstructed using microsynteny and phylogenetic analyses within the legume family. A single I TNL was present in Medicago truncatula and lost in soybean, mirroring the absence of complete BPMV resistance in soybean. Amplification of TNLs in the I cluster predates the divergence of the Phaseolus species, in agreement with the emergence of R-BPMV before the separation of the common bean wild centers of diversity. This analysis provides PCR-based markers useful in marker-assisted selection (MAS) and laid the foundation for cloning of R-BPMV resistance gene in order to transfer the resistance into soybean.

Natural variation in Arabidopsis responses to Plasmodiophora brassicae reveals an essential role for Resistance to Plasmodiophora brasssicae 1 (RPB1)

Despite the identification of clubroot resistance genes in various Brassica crops our understanding of the genetic basis of immunity to Plasmodiophora brassicae infection in the model plant Arabidopsis thaliana remains limited. To address this issue, we performed a screen of 142 natural accessions and identified 11 clubroot-resistant Arabidopsis lines. Genome-wide association analysis identified several genetic loci significantly linked with resistance. Three genes from two of these loci were targeted for deletion by CRISPR/Cas9 mutation in resistant accessions Est-1 and Uod-1. Deletion of Resistance to Plasmodiophora brassicae 1 (RPB1) rendered both lines susceptible to the P. brassicae pathotype P1+. Further analysis of rpb1 knock-out Est-1 and Uod-1 lines showed that the RPB1 protein is required for activation of downstream defence responses, such as the expression of phytoalexin biosynthesis gene CYP71A13. RPB1 has recently been shown to encode a cation channel localised in the endoplasmic reticulum. The clubroot susceptible Arabidopsis accession Col-0 lacks a functional RPB1 gene; when Col-0 is transformed with RPB1 expression driven by its native promoter it is capable of activating RPB1 transcription in response to infection, but this is not sufficient to confer resistance. Transient expression of RPB1 in Nicotiana tabacum induced programmed cell death in leaves. We conclude that RPB1 is a critical component of the defence response to P. brassicae infection in Arabidopsis, acting downstream of pathogen recognition but required for the elaboration of effective resistance.

The effect of DNA polymorphisms and natural variation on crossover hotspot activity in Arabidopsis hybrids

In hybrid organisms, genetically divergent homologous chromosomes pair and recombine during meiosis; however, the effect of specific types of polymorphisms on crossover is poorly understood. Here, to analyze this in Arabidopsis, we develop the seed-typing method that enables the massively parallel fine-mapping of crossovers by sequencing. We show that structural variants, observed in one of the generated intervals, do not change crossover frequency unless they are located directly within crossover hotspots. Both natural and Cas9-induced deletions result in lower hotspot activity but are not compensated by increases in immediately adjacent hotspots. To examine the effect of single nucleotide polymorphisms on crossover formation, we analyze hotspot activity in mismatch detection-deficient msh2 mutants. Surprisingly, polymorphic hotspots show reduced activity in msh2. In lines where only the hotspot-containing interval is heterozygous, crossover numbers increase above those in the inbred (homozygous). We conclude that MSH2 shapes crossover distribution by stimulating hotspot activity at polymorphic regions.

Role of pathogen's effectors in understanding host-pathogen interaction

Pathogens can pose challenges to plant growth and development at various stages of their life cycle. Two interconnected defense strategies prevent the growth of pathogens in plants, i.e., molecular patterns triggered immunity (PTI) and pathogenic effector-triggered immunity (ETI) that often provides resistance when PTI no longer functions as a result of pathogenic effectors. Plants may trigger an ETI defense response by directly or indirectly detecting pathogen effectors via their resistance proteins. A typical resistance protein is a nucleotide-binding receptor with leucine-rich sequences (NLRs) that undergo structural changes as they recognize their effectors and form associations with other NLRs. As a result of dimerization or oligomerization, downstream components activate "helper" NLRs, resulting in a response to ETI. It was thought that ETI is highly dependent on PTI. However, recent studies have found that ETI and PTI have symbiotic crosstalk, and both work together to create a robust system of plant defense. In this article, we have summarized the recent advances in understanding the plant's early immune response, its components, and how they cooperate in innate defense mechanisms. Moreover, we have provided the current perspective on engineering strategies for crop protection based on up-to-date knowledge.

DNA polymerase epsilon interacts with SUVH2/9 to repress the expression of genes associated with meiotic DSB hotspot in Arabidopsis

Meiotic recombination is initiated by the SPORULATION 11 (SPO11)-triggered formation of double-strand breaks (DSBs) that usually occur in open chromatin with active transcriptional features in many eukaryotes. However, gene transcription at DSB sites appears to be detrimental for repair, but the regulatory mechanisms governing transcription at meiotic DSB sites are largely undefined in plants. Here, we demonstrate that the largest DNA polymerase epsilon subunit POL2A interacts with SU(VAR)3 to 9 homologs SUVH2 and SUVH9. N-SIM (structured illumination microscopy) observation shows that the colocalization of SUVH2 with the meiotic DSB marker γ-H2AX is dependent on POL2A. RNA-seq of male meiocytes demonstrates that POL2A and SUVH2 jointly repress the expression of 865 genes, which have several known characteristics associated with meiotic DSB sites. Bisulfite-seq and small RNA-seq of male meiocytes support the idea that the silencing of these genes by POL2A and SUVH2/9 is likely independent of CHH methylation or 24-nt siRNA accumulation. Moreover, pol2a suvh2 suvh9 triple mutants have more severe defects in meiotic recombination and fertility compared with either pol2a or suvh2 suvh9. Our results not only identify a epigenetic regulatory mechanism for gene silencing in male meiocytes but also reveal roles for DNA polymerase and SUVH2/9 beyond their classic functions in mitosis.

Fast and Precise: How to Measure Meiotic Crossovers in Arabidopsis

During meiosis, homologous chromosomes (homologs) pair and undergo genetic recombination via assembly and disassembly of the synaptonemal complex. Meiotic recombination is initiated by excess formation of DNA double-strand breaks (DSBs), among which a subset are repaired by reciprocal genetic exchange, called crossovers (COs). COs generate genetic variations across generations, profoundly affecting genetic diversity and breeding. At least one CO between homologs is essential for the first meiotic chromosome segregation, but generally only one and fewer than three inter-homolog COs occur in plants. CO frequency and distribution are biased along chromosomes, suppressed in centromeres, and controlled by pro-CO, anti-CO, and epigenetic factors. Accurate and high-throughput detection of COs is important for our understanding of CO formation and chromosome behavior. Here, we review advanced approaches that enable precise measurement of the location, frequency, and genomic landscapes of COs in plants, with a focus on Arabidopsis thaliana.

Meiosis in crops: from genes to genomes

Meiosis is a key feature of sexual reproduction. During meiosis homologous chromosomes replicate, recombine, and randomly segregate, followed by the segregation of sister chromatids to produce haploid cells. The unique genotypes of recombinant gametes are an essential substrate for the selection of superior genotypes in natural populations and in plant breeding. In this review we summarize current knowledge on meiosis in diverse monocot and dicot crop species and provide a comprehensive resource of cloned meiotic mutants in six crop species (rice, maize, wheat, barley, tomato, and Brassica species). Generally, the functional roles of meiotic proteins are conserved between plant species, but we highlight notable differences in mutant phenotypes. The physical lengths of plant chromosomes vary greatly; for instance, wheat chromosomes are roughly one order of magnitude longer than those of rice. We explore how chromosomal distribution for crossover recombination can vary between species. We conclude that research on meiosis in crops will continue to complement that in Arabidopsis, and alongside possible applications in plant breeding will facilitate a better understanding of how the different stages of meiosis are controlled in plant species.

Natural variation identifies SNI1, the SMC5/6 component, as a modifier of meiotic crossover in Arabidopsis

Meiotic recombination plays a fundamental role in shaping genetic diversity in eukaryotes. Extensive variation in crossover rate exists between populations and species. The identity of modifier loci and their roles in genome evolution remain incompletely understood. We explored natural variation in Arabidopsis crossover and identified SNI1 as the causal gene underlying a major modifier locus. To date, SNI1 had no known role in crossover. SNI1 is a component of the SMC5/6 complex that is closely related to cohesin and condensin. Arabidopsis sni1 and other SMC5/6 mutants show similar effects on the interference-independent crossover pathway. Hence, our findings demonstrate that the SMC5/6 complex, which is known for its role in DNA damage repair, is also important for control of meiotic crossover.

The frequency and distribution of meiotic crossovers are tightly controlled; however, variation in this process can be observed both within and between species. Using crosses of two natural Arabidopsis thaliana accessions, Col and Ler, we mapped a crossover modifier locus to semidominant polymorphisms in SUPPRESSOR OF NPR1-1 INDUCIBLE 1 (SNI1), which encodes a component of the SMC5/6 complex. The sni1 mutant exhibits a modified pattern of recombination across the genome with crossovers elevated in chromosome distal regions but reduced in pericentromeres. Mutations in SNI1 result in reduced crossover interference and can partially restore the fertility of a Class I crossover pathway mutant, which suggests that the protein affects noninterfering crossover repair. Therefore, we tested genetic interactions between SNI1 and both RECQ4 and FANCM DNA helicases, which showed that additional Class II crossovers observed in the sni1 mutant are FANCM independent. Furthermore, genetic analysis of other SMC5/6 mutants confirms the observations of crossover redistribution made for SNI1. The study reveals the importance of the SMC5/6 complex in ensuring the proper progress of meiotic recombination in plants.

Selection for Plastic, Pathogen-Inducible Recombination in a Red Queen Model with Diploid Antagonists

Antagonistic interactions and co-evolution between a host and its parasite are known to cause oscillations in the population genetic structure of both species (Red Queen dynamics). Potentially, such oscillations may select for increased sex and recombination in the host, although theoretical models suggest that this happens under rather restricted values of selection intensity, epistasis, and other parameters. Here, we explore a model in which the diploid parasite succeeds to infect the diploid host only if their phenotypes at the interaction-mediating loci match. Whenever regular oscillations emerge in this system, we test whether plastic, pathogen-inducible recombination in the host can be favored over the optimal constant recombination. Two forms of the host recombination dependence on the parasite pressure were considered: either proportionally to the risk of infection (prevention strategy) or upon the fact of infection (remediation strategy). We show that both forms of plastic recombination can be favored, although relatively infrequently (up to 11% of all regimes with regular oscillations, and up to 20% of regimes with obligate parasitism). This happens under either strong overall selection and high recombination rate in the host, or weak overall selection and low recombination rate in the host. In the latter case, the system’s dynamics are considerably more complex. The prevention strategy is favored more often than the remediation one. It is noteworthy that plastic recombination can be favored even when any constant recombination is rejected, making plasticity an evolutionary mechanism for the rescue of host recombination.

Plant NLR diversity: the known unknowns of pan-NLRomes

Plants and pathogens constantly adapt to each other. As a consequence, many members of the plant immune system, and especially the intracellular nucleotide-binding site leucine-rich repeat receptors, also known as NOD-like receptors (NLRs), are highly diversified, both among family members in the same genome, and between individuals in the same species. While this diversity has long been appreciated, its true extent has remained unknown. With pan-genome and pan-NLRome studies becoming more and more comprehensive, our knowledge of NLR sequence diversity is growing rapidly, and pan-NLRomes provide powerful platforms for assigning function to NLRs. These efforts are an important step toward the goal of comprehensively predicting from sequence alone whether an NLR provides disease resistance, and if so, to which pathogens.

AFLAP: assembly-free linkage analysis pipeline using k-mers from genome sequencing data

Our assembly-free linkage analysis pipeline (AFLAP) identifies segregating markers as k-mers in the raw reads without using a reference genome assembly for calling variants and provides genotype tables for the construction of unbiased, high-density genetic maps without a genome assembly. AFLAP is validated and contrasted to a conventional workflow using simulated data. AFLAP is applied to whole genome sequencing and genotype-by-sequencing data of F1, F2, and recombinant inbred populations of two different plant species, producing genetic maps that are concordant with genome assemblies. The AFLAP-based genetic map for Bremia lactucae enables the production of a chromosome-scale genome assembly.

HIGH CROSSOVER RATE1 encodes PROTEIN PHOSPHATASE X1 and restricts meiotic crossovers in Arabidopsis

Meiotic crossovers are tightly restricted in most eukaryotes, despite an excess of initiating DNA double-strand breaks. The majority of plant crossovers are dependent on class I interfering repair, with a minority formed via the class II pathway. Class II repair is limited by anti-recombination pathways; however, similar pathways repressing class I crossovers have not been identified. Here, we performed a forward genetic screen in Arabidopsis using fluorescent crossover reporters to identify mutants with increased or decreased recombination frequency. We identified HIGH CROSSOVER RATE1 (HCR1) as repressing crossovers and encoding PROTEIN PHOSPHATASE X1. Genome-wide analysis showed that hcr1 crossovers are increased in the distal chromosome arms. MLH1 foci significantly increase in hcr1 and crossover interference decreases, demonstrating an effect on class I repair. Consistently, yeast two-hybrid and in planta assays show interaction between HCR1 and class I proteins, including HEI10, PTD, MSH5 and MLH1. We propose that HCR1 plays a major role in opposition to pro-recombination kinases to restrict crossovers in Arabidopsis.

A Truncated Singleton NLR Causes Hybrid Necrosis in Arabidopsis thaliana

Hybrid necrosis in plants arises from conflict between divergent alleles of immunity genes contributed by different parents, resulting in autoimmunity. We investigate a severe hybrid necrosis case in Arabidopsis thaliana, where the hybrid does not develop past the cotyledon stage and dies 3 weeks after sowing. Massive transcriptional changes take place in the hybrid, including the upregulation of most NLR (nucleotide-binding site leucine-rich repeat) disease-resistance genes. This is due to an incompatible interaction between the singleton TIR-NLR gene DANGEROUS MIX 10 (DM10), which was recently relocated from a larger NLR cluster, and an unlinked locus, DANGEROUS MIX 11 (DM11). There are multiple DM10 allelic variants in the global A. thaliana population, several of which have premature stop codons. One of these, which has a truncated LRR-PL (leucine-rich repeat [LRR]-post-LRR) region, corresponds to the DM10 risk allele. The DM10 locus and the adjacent genomic region in the risk allele carriers are highly differentiated from those in the nonrisk carriers in the global A. thaliana population, suggesting that this allele became geographically widespread only relatively recently. The DM11 risk allele is much rarer and found only in two accessions from southwestern Spain-a region from which the DM10 risk haplotype is absent-indicating that the ranges of DM10 and DM11 risk alleles may be nonoverlapping.

Genome-Wide Identification of Populus Malectin/Malectin-Like Domain-Containing Proteins and Expression Analyses Reveal Novel Candidates for Signaling and Regulation of Wood Development

Malectin domain (MD) is a ligand-binding protein motif of pro- and eukaryotes. It is particularly abundant in Viridiplantae, where it occurs as either a single (MD, PF11721) or tandemly duplicated domain (PF12819) called malectin-like domain (MLD). In herbaceous plants, MD- or MLD-containing proteins (MD proteins) are known to regulate development, reproduction, and resistance to various stresses. However, their functions in woody plants have not yet been studied. To unravel their potential role in wood development, we carried out genome-wide identification of MD proteins in the model tree species black cottonwood (Populus trichocarpa), and analyzed their expression and co-expression networks. P. trichocarpa had 146 MD genes assigned to 14 different clades, two of which were specific to the genus Populus. 87% of these genes were located on chromosomes, the rest being associated with scaffolds. Based on their protein domain organization, and in agreement with the exon-intron structures, the MD genes identified here could be classified into five superclades having the following domains: leucine-rich repeat (LRR)-MD-protein kinase (PK), MLD-LRR-PK, MLD-PK (CrRLK1L), MLD-LRR, and MD-Kinesin. Whereas the majority of MD genes were highly expressed in leaves, particularly under stress conditions, eighteen showed a peak of expression during secondary wall formation in the xylem and their co-expression networks suggested signaling functions in cell wall integrity, pathogen-associated molecular patterns, calcium, ROS, and hormone pathways. Thus, P. trichocarpa MD genes having different domain organizations comprise many genes with putative foliar defense functions, some of which could be specific to Populus and related species, as well as genes with potential involvement in signaling pathways in other tissues including developing wood.

MSH2 shapes the meiotic crossover landscape in relation to interhomolog polymorphism in Arabidopsis

During meiosis, DNA double-strand breaks undergo interhomolog repair to yield crossovers between homologous chromosomes. To investigate how interhomolog sequence polymorphism affects crossovers, we sequenced multiple recombinant populations of the model plant Arabidopsis thaliana. Crossovers were elevated in the diverse pericentromeric regions, showing a local preference for polymorphic regions. We provide evidence that crossover association with elevated diversity is mediated via the Class I crossover formation pathway, although very high levels of diversity suppress crossovers. Interhomolog polymorphism causes mismatches in recombining molecules, which can be detected by MutS homolog (MSH) mismatch repair protein heterodimers. Therefore, we mapped crossovers in a msh2 mutant, defective in mismatch recognition, using multiple hybrid backgrounds. Although total crossover numbers were unchanged in msh2 mutants, recombination was remodelled from the diverse pericentromeres towards the less-polymorphic sub-telomeric regions. Juxtaposition of megabase heterozygous and homozygous regions causes crossover remodelling towards the heterozygous regions in wild type Arabidopsis, but not in msh2 mutants. Immunostaining showed that MSH2 protein accumulates on meiotic chromosomes during prophase I, consistent with MSH2 regulating meiotic recombination. Our results reveal a pro-crossover role for MSH2 in regions of higher sequence diversity in A. thaliana.

Deciphering the Impact of a Bacterial Infection on Meiotic Recombination in Arabidopsis with Fluorescence Tagged Lines

Plants are under strong evolutionary pressure to maintain surveillance against pathogens. One major disease resistance mechanism is based on NB-LRR (NLR) proteins that specifically recognize pathogen effectors. The cluster organization of the NLR gene family could favor sequence exchange between NLR genes via recombination, favoring their evolutionary dynamics. Increasing data, based on progeny analysis, suggest the existence of a link between the perception of biotic stress and the production of genetic diversity in the offspring. This could be driven by an increased rate of meiotic recombination in infected plants, but this has never been strictly demonstrated. In order to test if pathogen infection can increase DNA recombination in pollen meiotic cells, we infected Arabidopsis Fluorescent Tagged Lines (FTL) with the virulent bacteria Pseudomonas syringae. We measured the meiotic recombination rate in two regions of chromosome 5, containing or not an NLR gene cluster. In all tested intervals, no significant difference in genetic recombination frequency between infected and control plants was observed. Although it has been reported that pathogen exposure can sometimes increase the frequency of recombinant progeny in plants, our findings suggest that meiotic recombination rate in Arabidopsis may be resilient to at least some pathogen attack. Alternative mechanisms are discussed.

ASY1 acts as a dosage-dependent antagonist of telomere-led recombination and mediates crossover interference in Arabidopsis

During meiosis, interhomolog recombination produces crossovers and noncrossovers to create genetic diversity. Meiotic recombination frequency varies at multiple scales, with high subtelomeric recombination and suppressed centromeric recombination typical in many eukaryotes. During recombination, sister chromatids are tethered as loops to a polymerized chromosome axis, which, in plants, includes the ASY1 HORMA domain protein and REC8-cohesin complexes. Using chromatin immunoprecipitation, we show an ascending telomere-to-centromere gradient of ASY1 enrichment, which correlates strongly with REC8-cohesin ChIP-seq data. We mapped crossovers genome-wide in the absence of ASY1 and observe that telomere-led recombination becomes dominant. Surprisingly, asy1/+ heterozygotes also remodel crossovers toward subtelomeric regions at the expense of the pericentromeres. Telomeric recombination increases in asy1/+ occur in distal regions where ASY1 and REC8 ChIP enrichment are lowest in wild type. In wild type, the majority of crossovers show interference, meaning that they are more widely spaced along the chromosomes than expected by chance. To measure interference, we analyzed double crossover distances, MLH1 foci, and fluorescent pollen tetrads. Interestingly, while crossover interference is normal in asy1/+, it is undetectable in asy1 mutants, indicating that ASY1 is required to mediate crossover interference. Together, this is consistent with ASY1 antagonizing telomere-led recombination and promoting spaced crossover formation along the chromosomes via interference. These findings provide insight into the role of the meiotic axis in patterning recombination frequency within plant genomes.

Structural evolution drives diversification of the large LRR-RLK gene family

●Cells are continuously exposed to chemical signals that they must discriminate between and respond to appropriately. In embryophytes, the leucine-rich repeat receptor-like kinases (LRR-RLKs) are signal receptors critical in development and defense. LRR-RLKs have diversified to hundreds of genes in many plant genomes. Although intensively studied, a well-resolved LRR-RLK gene tree has remained elusive. ●To resolve the LRR-RLK gene tree, we developed an improved gene discovery method based on iterative hidden Markov model searching and phylogenetic inference. We used this method to infer complete gene trees for each of the LRR-RLK subclades and reconstructed the deepest nodes of the full gene family. ●We discovered that the LRR-RLK gene family is even larger than previously thought, and that protein domain gains and losses are prevalent. These structural modifications, some of which likely predate embryophyte diversification, led to misclassification of some LRR-RLK variants as members of other gene families. Our work corrects this misclassification. ●Our results reveal ongoing structural evolution generating novel LRR-RLK genes. These new genes are raw material for the diversification of signaling in development and defense. Our methods also enable phylogenetic reconstruction in any large gene family.

Interacting Genomic Landscapes of REC8-Cohesin, Chromatin, and Meiotic Recombination in Arabidopsis

Meiosis recombines genetic variation and influences eukaryote genome evolution. During meiosis, DNA double-strand breaks (DSBs) enter interhomolog repair to yield crossovers and noncrossovers. DSB repair occurs as replicated sister chromatids are connected to a polymerized axis. Cohesin rings containing the REC8 kleisin subunit bind sister chromatids and anchor chromosomes to the axis. Here, we report the genomic landscape of REC8 using chromatin immunoprecipitation sequencing (ChIP-seq) in Arabidopsis (Arabidopsis thaliana). REC8 associates with regions of high nucleosome occupancy in multiple chromatin states, including histone methylation at H3K4 (expressed genes), H3K27 (silent genes), and H3K9 (silent transposons). REC8 enrichment is associated with suppression of meiotic DSBs and crossovers at the chromosome and fine scales. As REC8 enrichment is greatest in transposon-dense heterochromatin, we repeated ChIP-seq in kyp suvh5 suvh6 H3K9me2 mutants. Surprisingly, REC8 enrichment is maintained in kyp suvh5 suvh6 heterochromatin and no defects in centromeric cohesion were observed. REC8 occupancy within genes anti-correlates with transcription and is reduced in COPIA transposons that reactivate expression in kyp suvh5 suvh6 Abnormal axis structures form in rec8 that recruit DSB-associated protein foci and undergo synapsis, which is followed by chromosome fragmentation. Therefore, REC8 occupancy correlates with multiple chromatin states and is required to organize meiotic chromosome architecture and interhomolog recombination.

Architecture of Parallel Adaptation in Ten Lacustrine Threespine Stickleback Populations from the White Sea Area

Adaptation of threespine stickleback to freshwater involves parallel recruitment of freshwater alleles in clusters of closely linked sites, or divergence islands (DIs). However, it remains unclear to what extent the DIs and the alleles that constitute them coincide between populations that underwent adaptation to freshwater independently. We examine threespine sticklebacks from ten freshwater lakes that emerged 500–1500 years ago in the White Sea basin, with the emphasis on repeatability of genomic patterns of adaptation among the lake populations and the role of local recombination rate in the distribution and structure of DIs. The 65 detected DIs are clustered in the genome, forming 12 aggregations, and this clustering cannot be explained by the variation of the recombination rate. Only 21 of the DIs are present in all the freshwater populations, likely being indispensable for successful colonization of freshwater environment by the ancestral marine population. Within most DIs, the same set of single nucleotide polymorphisms (SNPs) distinguish marine and freshwater haplotypes in all the lake populations; however, in some DIs, freshwater alleles differ between populations, suggesting that they could have been established by recruitment of different haplotypes in different populations.

Chromosomal assembly and analyses of genome-wide recombination rates in the forest pathogenic fungus Armillaria ostoyae

Recombination shapes the evolutionary trajectory of populations and plays an important role in the faithful transmission of chromosomes during meiosis. Levels of sexual reproduction and recombination are important properties of host-pathogen interactions because the speed of antagonistic co-evolution depends on the ability of hosts and pathogens to generate genetic variation. However, our understanding of the importance of recombination is limited because large taxonomic groups remain poorly investigated. Here, we analyze recombination rate variation in the basidiomycete fungus Armillaria ostoyae, which is an aggressive pathogen on a broad range of conifers and other trees. We analyzed a previously constructed, dense genetic map based on 198 single basidiospore progeny from a cross. Progeny were genotyped at a genome-wide set of single-nucleotide polymorphism (SNP) markers using double digest restriction site associated DNA sequencing. Based on a linkage map of on 11,700 SNPs spanning 1007.5 cM, we assembled genomic scaffolds into 11 putative chromosomes of a total genome size of 56.6 Mb. We identified 1984 crossover events among all progeny and found that recombination rates were highly variable along chromosomes. Recombination hotspots tended to be in regions close to the telomeres and were more gene-poor than the genomic background. Genes in proximity to recombination hotspots were encoding on average shorter proteins and were enriched for pectin degrading enzymes. Our analyses enable more powerful population and genome-scale studies of a major tree pathogen.

An Ultra High-Density Arabidopsis thaliana Crossover Map That Refines the Influences of Structural Variation and Epigenetic Features

Many environmental, genetic, and epigenetic factors are known to affect the frequency and positioning of meiotic crossovers (COs). Suppression of COs by large, cytologically visible inversions and translocations has long been recognized, but relatively little is known about how smaller structural variants (SVs) affect COs. To examine fine-scale determinants of the CO landscape, including SVs, we used a rapid, cost-effective method for high-throughput sequencing to generate a precise map of >17,000 COs between the Col-0 and Ler-0 accessions of Arabidopsis thaliana COs were generally suppressed in regions with SVs, but this effect did not depend on the size of the variant region, and was only marginally affected by the variant type. CO suppression did not extend far beyond the SV borders and CO rates were slightly elevated in the flanking regions. Disease resistance gene clusters, which often exist as SVs, exhibited high CO rates at some loci, but there was a tendency toward depressed CO rates at loci where large structural differences exist between the two parents. Our high-density map also revealed in fine detail how CO positioning relates to genetic (DNA motifs) and epigenetic (chromatin structure) features of the genome. We conclude that suppression of COs occurs over a narrow region spanning large- and small-scale SVs, representing an influence on the CO landscape in addition to sequence and epigenetic variation along chromosomes.

The cytological basis of powdery mildew resistance in wild Chinese Vitis species

The wild Chinese grapevines (Vitis spp.) show varying levels of resistance to powdery mildew caused by Erysiphe necator that is an economically important disease of cultivated grapevines (Vitis vinifera). However, little information is available regarding the cytological mechanisms of powdery mildew resistance in these wild relatives. Here, we studied the cytological responses of three wild Chinese grapevine accessions after they were infected with E. necator (En) NAFU1 in comparison to the susceptible V. vinifera cv. 'Thompson Seedless' grape. The hyphal growth and sporulation of En NAFU1 were significantly restricted in wild species compared to 'Thompson Seedless', which appears to be associated with early cell wall deposition at the attempt sites, encasement of haustoria, and hypersensitive response-like cell death of penetrated epidermal cells. Moreover, endogenous free salicylic acid (SA) was more abundant in wild Chinese Vitis species than in 'Thompson Seedless' under pathogen-free condition. During En NAFU1 colonization, SA conjugates accumulated higher in wild grapevines than in 'Thompson Seedless'. In addition, the species-specific expression patterns of defense-associated genes during En NAFU1 colonization indicated that mechanisms underlying powdery mildew resistance are divergent among different wild Chinese Vitis species. These results contribute to understanding of mechanisms underlying defense responses of wild Chinese Vitis species against powdery mildew.

Linked-read sequencing of gametes allows efficient genome-wide analysis of meiotic recombination

Meiotic crossovers (COs) ensure proper chromosome segregation and redistribute the genetic variation that is transmitted to the next generation. Large populations and the demand for genome-wide, fine-scale resolution challenge existing methods for CO identification. Taking advantage of linked-read sequencing, we develop a highly efficient method for genome-wide identification of COs at kilobase resolution in pooled recombinants. We first test this method using a pool of Arabidopsis F₂ recombinants, and recapitulate results obtained from the same plants using individual whole-genome sequencing. By applying this method to a pool of pollen DNA from an F₁ plant, we establish a highly accurate CO landscape without generating or sequencing a single recombinant plant. The simplicity of this approach enables the simultaneous generation and analysis of multiple CO landscapes, accelerating the pace at which mechanisms for the regulation of recombination can be elucidated through efficient comparisons of genotypic and environmental effects on recombination.

Meiotic crossovers (COs) generate genetic variation and ensure proper chromosome segregation. Here, the authors develop a method for identifying COs at kilobase resolution in pooled recombinants using linked-read sequencing data, and apply it to investigate genome-wide CO landscapes of Arabidopsis thaliana.

Linked-read sequencing of gametes allows efficient genome-wide analysis of meiotic recombination

Meiotic crossovers (COs) ensure proper chromosome segregation and redistribute the genetic variation that is transmitted to the next generation. Large populations and the demand for genome-wide, fine-scale resolution challenge existing methods for CO identification. Taking advantage of linked-read sequencing, we develop a highly efficient method for genome-wide identification of COs at kilobase resolution in pooled recombinants. We first test this method using a pool of Arabidopsis F₂ recombinants, and recapitulate results obtained from the same plants using individual whole-genome sequencing. By applying this method to a pool of pollen DNA from an F₁ plant, we establish a highly accurate CO landscape without generating or sequencing a single recombinant plant. The simplicity of this approach enables the simultaneous generation and analysis of multiple CO landscapes, accelerating the pace at which mechanisms for the regulation of recombination can be elucidated through efficient comparisons of genotypic and environmental effects on recombination.

Natural Variation in TBP-ASSOCIATED FACTOR 4b Controls Meiotic Crossover and Germline Transcription in Arabidopsis

Meiotic crossover frequency varies within genomes, which influences genetic diversity and adaptation. In turn, genetic variation within populations can act to modify crossover frequency in cis and trans. To identify genetic variation that controls meiotic crossover frequency, we screened Arabidopsis accessions using fluorescent recombination reporters. We mapped a genetic modifier of crossover frequency in Col × Bur populations of Arabidopsis to a premature stop codon within TBP-ASSOCIATED FACTOR 4b (TAF4b), which encodes a subunit of the RNA polymerase II general transcription factor TFIID. The Arabidopsis taf4b mutation is a rare variant found in the British Isles, originating in South-West Ireland. Using genetics, genomics, and immunocytology, we demonstrate a genome-wide decrease in taf4b crossovers, with strongest reduction in the sub-telomeric regions. Using RNA sequencing (RNA-seq) from purified meiocytes, we show that TAF4b expression is meiocyte enriched, whereas its paralog TAF4 is broadly expressed. Consistent with the role of TFIID in promoting gene expression, RNA-seq of wild-type and taf4b meiocytes identified widespread transcriptional changes, including in genes that regulate the meiotic cell cycle and recombination. Therefore, TAF4b duplication is associated with acquisition of meiocyte-specific expression and promotion of germline transcription, which act directly or indirectly to elevate crossovers. This identifies a novel mode of meiotic recombination control via a general transcription factor.

Heterogeneous transposable elements as silencers, enhancers and targets of meiotic recombination

During meiosis, DNA double-strand breaks are initiated by the topoisomerase-like enzyme SPO11 and are repaired by inter-sister chromatid and inter-homologue DNA repair pathways. Genome-wide maps of initiating DNA double-strand breaks and inter-homologue repair events are now available for a number of mammalian, fungal and plant species. In mammals, PRDM9 specifies the location of meiotic recombination initiation via recognition of specific DNA sequence motifs by its C2H2 zinc finger array. In fungi and plants, meiotic recombination appears to be initiated less discriminately in accessible chromatin, including at gene promoters. Generally, meiotic crossover is suppressed in highly repetitive genomic regions that are made up of transposable elements (TEs), to prevent deleterious non-allelic homologous recombination events. However, recent and older studies have revealed intriguing relationships between meiotic recombination initiation and repair, and transposable elements. For instance, gene conversion events have been detected in maize centromeric retroelements, mouse MULE-MuDR DNA transposons undergo substantial meiotic recombination initiation, Arabidopsis Helitron TEs are among the hottest of recombination initiation hotspots, and human TE sequences can modify the crossover rate at adjacent PRDM9 motifs in cis. Here, we summarize the relationship between meiotic recombination and TEs, discuss recent insights from highly divergent eukaryotes and highlight outstanding questions in the field.

Ancient polymorphisms contribute to genome-wide variation by long-term balancing selection and divergent sorting in Boechera stricta

BACKGROUND

Genomic variation is widespread, and both neutral and selective processes can generate similar patterns in the genome. These processes are not mutually exclusive, so it is difficult to infer the evolutionary mechanisms that govern population and species divergence. Boechera stricta is a perennial relative of Arabidopsis thaliana native to largely undisturbed habitats with two geographic and ecologically divergent subspecies. Here, we delineate the evolutionary processes driving the genetic diversity and population differentiation in this species.

RESULTS

Using whole-genome re-sequencing data from 517 B. stricta accessions, we identify four genetic groups that diverged around 30-180 thousand years ago, with long-term small effective population sizes and recent population expansion after the Last Glacial Maximum. We find three genomic regions with elevated nucleotide diversity, totaling about 10% of the genome. These three regions of elevated nucleotide diversity show excess of intermediate-frequency alleles, higher absolute divergence (dXY), and lower relative divergence (FST) than genomic background, and significant enrichment in immune-related genes, reflecting long-term balancing selection. Scattered across the genome, we also find regions with both high FST and dXY among the groups, termed FST-islands. Population genetic signatures indicate that FST-islands with elevated divergence, which have experienced directional selection, are derived from divergent sorting of ancient polymorphisms.

CONCLUSIONS

Our results suggest that long-term balancing selection on disease resistance genes may have maintained ancestral haplotypes across different geographical lineages, and unequal sorting of balanced polymorphisms may have generated genomic regions with elevated divergence. This study highlights the importance of ancestral balanced polymorphisms as crucial components of genome-wide variation.

Complex Relationships between Chromatin Accessibility, Sequence Divergence, and Gene Expression in Arabidopsis thaliana

Variation in regulatory DNA is thought to drive phenotypic variation, evolution, and disease. Prior studies of regulatory DNA and transcription factors across animal species highlighted a fundamental conundrum: Transcription factor binding domains and cognate binding sites are conserved, while regulatory DNA sequences are not. It remains unclear how conserved transcription factors and dynamic regulatory sites produce conserved expression patterns across species. Here, we explore regulatory DNA variation and its functional consequences within Arabidopsis thaliana, using chromatin accessibility to delineate regulatory DNA genome-wide. Unlike in previous cross-species comparisons, the positional homology of regulatory DNA is maintained among A. thaliana ecotypes and less nucleotide divergence has occurred. Of the ∼50,000 regulatory sites in A. thaliana, we found that 15% varied in accessibility among ecotypes. Some of these accessibility differences were associated with extensive, previously unannotated sequence variation, encompassing many deletions and ancient hypervariable alleles. Unexpectedly, for the majority of such regulatory sites, nearby gene expression was unaffected. Nevertheless, regulatory sites with high levels of sequence variation and differential chromatin accessibility were the most likely to be associated with differential gene expression. Finally, and most surprising, we found that the vast majority of differentially accessible sites show no underlying sequence variation. We argue that these surprising results highlight the necessity to consider higher-order regulatory context in evaluating regulatory variation and predicting its phenotypic consequences.

Meiocyte-Specific and AtSPO11-1-Dependent Small RNAs and Their Association with Meiotic Gene Expression and Recombination

Meiotic recombination ensures accurate chromosome segregation and results in genetic diversity in sexually reproducing eukaryotes. Over the last few decades, the genetic regulation of meiotic recombination has been extensively studied in many organisms. However, the role of endogenous meiocyte-specific small RNAs (ms-sRNAs; 21-24 nucleotide [nt]) and their involvement in meiotic recombination are unclear. Here, we sequenced the total small RNA (sRNA) and messenger RNA populations from meiocytes and leaves of wild type Arabidopsis (Arabidopsis thaliana) and meiocytes of spo11-1, a mutant defective in double-strand break formation, and we discovered 2,409 ms-sRNA clusters, 1,660 of which areSPORULATION 11-1 (AtSPO11-1)-dependent. Unlike mitotic small interfering RNAs that are enriched in intergenic regions and associated with gene silencing, ms-sRNAs are significantly enriched in genic regions and exhibit a positive correlation with genes that are preferentially expressed in meiocytes (i.e. Arabidopsis SKP1-LIKE1 and RAD51), in a fashion unrelated to DNA methylation. We also found that AtSPO11-1-dependent sRNAs have distinct characteristics compared with ms-sRNAs and tend to be associated with two known types of meiotic recombination hotspot motifs (i.e. CTT-repeat and A-rich motifs). These results reveal different meiotic and mitotic sRNA landscapes and provide new insights into how sRNAs relate to gene expression in meiocytes and meiotic recombination.

Interhomolog polymorphism shapes meiotic crossover within the Arabidopsis RAC1 and RPP13 disease resistance genes

During meiosis, chromosomes undergo DNA double-strand breaks (DSBs), which can be repaired using a homologous chromosome to produce crossovers. Meiotic recombination frequency is variable along chromosomes and tends to concentrate in narrow hotspots. We mapped crossover hotspots located in the Arabidopsis thaliana RAC1 and RPP13 disease resistance genes, using varying haplotypic combinations. We observed a negative non-linear relationship between interhomolog divergence and crossover frequency within the hotspots, consistent with polymorphism locally suppressing crossover repair of DSBs. The fancm, recq4a recq4b, figl1 and msh2 mutants, or lines with increased HEI10 dosage, are known to show increased crossovers throughout the genome. Surprisingly, RAC1 crossovers were either unchanged or decreased in these genetic backgrounds, showing that chromosome location and local chromatin environment are important for regulation of crossover activity. We employed deep sequencing of crossovers to examine recombination topology within RAC1, in wild type, fancm, recq4a recq4b and fancm recq4a recq4b backgrounds. The RAC1 recombination landscape was broadly conserved in the anti-crossover mutants and showed a negative relationship with interhomolog divergence. However, crossovers at the RAC1 5'-end were relatively suppressed in recq4a recq4b backgrounds, further indicating that local context may influence recombination outcomes. Our results demonstrate the importance of interhomolog divergence in shaping recombination within plant disease resistance genes and crossover hotspots.

Where to Cross Over? Defining Crossover Sites in Plants

It is believed that recombination in meiosis serves to reshuffle genetic material from both parents to increase genetic variation in the progeny. At the same time, the number of crossovers is usually kept at a very low level. As a consequence, many organisms need to make the best possible use from the one or two crossovers that occur per chromosome in meiosis. From this perspective, the decision of where to allocate rare crossover events becomes an important issue, especially in self-pollinating plant species, which experience limited variation due to inbreeding. However, the freedom in crossover allocation is significantly limited by other, genetic and non-genetic factors, including chromatin structure. Here we summarize recent progress in our understanding of those processes with a special emphasis on plant genomes. First, we focus on factors which influence the distribution of recombination initiation sites and discuss their effects at both, the single hotspot level and at the chromosome scale. We also briefly explain the aspects of hotspot evolution and their regulation. Next, we analyze how recombination initiation sites translate into the development of crossovers and their location. Moreover, we provide an overview of the sequence polymorphism impact on crossover formation and chromosomal distribution.

A Critical Assessment of 60 Years of Maize Intragenic Recombination

Until the mid-1950s, it was believed that genetic crossovers did not occur within genes. Crossovers occurred between genes, the "beads on a string" model. Then in 1956, Seymour Benzer published his classic paper describing crossing over within a gene, intragenic recombination. This result from a bacteriophage gene prompted Oliver Nelson to study intragenic recombination in the maize Waxy locus. His studies along with subsequent work by others working with maize and other organisms described the outcomes of intragenic recombination and provided some of the earliest evidence that genes, not intergenic regions, were recombination hotspots. High-throughput genotyping approaches have since replaced single gene intragenic studies for characterizing the outcomes of recombination. These large-scale studies confirm that genes, or more generally genic regions, are the most active recombinogenic regions, and suggested a pattern of crossovers similar to the budding yeast Saccharomyces cerevisiae. In S. cerevisiae recombination is initiated by double-strand breaks (DSBs) near transcription start sites (TSSs) of genes producing a polarity gradient where crossovers preferentially resolve at the 5' end of genes. Intragenic studies in maize yielded less evidence for either polarity or for DSBs near TSSs initiating recombination and in certain respects resembled Schizosaccharomyces pombe or mouse. These different perspectives highlight the need to draw upon the strengths of different approaches and caution against relying on a single model system or approach for understanding recombination.

Low-coverage resequencing detects meiotic recombination pattern and features in tomato RILs

Traditional plant breeding relies on meiotic recombination for mixing of parental alleles to create novel allele combinations. Detailed analysis of recombination patterns in model organisms shows that recombination is tightly regulated within the genome, but frequencies vary extensively along chromosomes. Despite being a model organism for fruit developmental studies, high-resolution recombination patterns are lacking in tomato. In this study, we developed a novel methodology to use low-coverage resequencing to identify genome-wide recombination patterns and applied this methodology on 60 tomato Recombinant Inbred Lines (RILs). Our methodology identifies polymorphic markers from the low-coverage resequencing population data and utilizes the same data to locate the recombination breakpoints in individuals by using a variable sliding window. We identified 1,445 recombination sites comprising 112 recombination prone regions enriched for AT-rich DNA motifs. Furthermore, the recombination prone regions in tomato preferably occurred in gene promoters over intergenic regions, an observation consistent with Arabidopsis thaliana, Zea mays and Mimulus guttatus. Overall, our cost effective method and findings enhance the understanding of meiotic recombination in tomato and suggest evolutionarily conserved recombination associated genomic features.

High-resolution crossover mapping reveals similarities and differences of male and female recombination in maize

Meiotic crossovers (COs) are not uniformly distributed across the genome. Factors affecting this phenomenon are not well understood. Although many species exhibit large differences in CO numbers between sexes, sex-specific aspects of CO landscape are particularly poorly elucidated. Here, we conduct high-resolution CO mapping in maize. Our results show that CO numbers as well as their overall distribution are similar in male and female meioses. There are, nevertheless, dissimilarities at local scale. Male and female COs differ in their locations relative to transcription start sites in gene promoters and chromatin marks, including nucleosome occupancy and tri-methylation of lysine 4 of histone H3 (H3K4me3). Our data suggest that sex-specific factors not only affect male–female CO number disparities but also cause fine differences in CO positions. Differences between male and female CO landscapes indicate that recombination has distinct implications for population structure and gene evolution in male and in female meioses.

Sex-specific meiotic crossover (CO) landscapes have been identified in multiple species. Here, the authors show that male and female meioses in maize have similar CO landscapes, and differences between COs in the two sexes only exists in their location relative to transcription start sites and some chromatin marks.

Variation in recombination frequency and distribution across eukaryotes: patterns and processes

Recombination, the exchange of DNA between maternal and paternal chromosomes during meiosis, is an essential feature of sexual reproduction in nearly all multicellular organisms. While the role of recombination in the evolution of sex has received theoretical and empirical attention, less is known about how recombination rate itself evolves and what influence this has on evolutionary processes within sexually reproducing organisms. Here, we explore the patterns of, and processes governing recombination in eukaryotes. We summarize patterns of variation, integrating current knowledge with an analysis of linkage map data in 353 organisms. We then discuss proximate and ultimate processes governing recombination rate variation and consider how these influence evolutionary processes. Genome-wide recombination rates (cM/Mb) can vary more than tenfold across eukaryotes, and there is large variation in the distribution of recombination events across closely related taxa, populations and individuals. We discuss how variation in rate and distribution relates to genome architecture, genetic and epigenetic mechanisms, sex, environmental perturbations and variable selective pressures. There has been great progress in determining the molecular mechanisms governing recombination, and with the continued development of new modelling and empirical approaches, there is now also great opportunity to further our understanding of how and why recombination rate varies.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Nucleosomes and DNA methylation shape meiotic DSB frequency in Arabidopsis thaliana transposons and gene regulatory regions

Meiotic recombination initiates from DNA double-strand breaks (DSBs) generated by SPO11 topoisomerase-like complexes. Meiotic DSB frequency varies extensively along eukaryotic chromosomes, with hotspots controlled by chromatin and DNA sequence. To map meiotic DSBs throughout a plant genome, we purified and sequenced Arabidopsis thaliana SPO11-1-oligonucleotides. SPO11-1-oligos are elevated in gene promoters, terminators, and introns, which is driven by AT-sequence richness that excludes nucleosomes and allows SPO11-1 access. A positive relationship was observed between SPO11-1-oligos and crossovers genome-wide, although fine-scale correlations were weaker. This may reflect the influence of interhomolog polymorphism on crossover formation, downstream from DSB formation. Although H3K4me3 is enriched in proximity to SPO11-1-oligo hotspots at gene 5' ends, H3K4me3 levels do not correlate with DSBs. Repetitive transposons are thought to be recombination silenced during meiosis, to prevent nonallelic interactions and genome instability. Unexpectedly, we found high SPO11-1-oligo levels in nucleosome-depleted Helitron/Pogo/Tc1/Mariner DNA transposons, whereas retrotransposons were coldspots. High SPO11-1-oligo transposons are enriched within gene regulatory regions and in proximity to immunity genes, suggesting a role as recombination enhancers. As transposon mobility in plant genomes is restricted by DNA methylation, we used the met1 DNA methyltransferase mutant to investigate the role of heterochromatin in SPO11-1-oligo distributions. Epigenetic activation of meiotic DSBs in proximity to centromeres and transposons occurred in met1 mutants, coincident with reduced nucleosome occupancy, gain of transcription, and H3K4me3. Together, our work reveals a complex relationship between chromatin and meiotic DSBs within A. thaliana genes and transposons, with significance for the diversity and evolution of plant genomes.

Epigenetic activation of meiotic recombination near Arabidopsis thaliana centromeres via loss of H3K9me2 and non-CG DNA methylation

Eukaryotic centromeres contain the kinetochore, which connects chromosomes to the spindle allowing segregation. During meiosis, centromeres are suppressed for inter-homolog crossover, as recombination in these regions can cause chromosome missegregation and aneuploidy. Plant centromeres are surrounded by transposon-dense pericentromeric heterochromatin that is epigenetically silenced by histone 3 lysine 9 dimethylation (H3K9me2), and DNA methylation in CG and non-CG sequence contexts. However, the role of these chromatin modifications in control of meiotic recombination in the pericentromeres is not fully understood. Here, we show that disruption of Arabidopsis thaliana H3K9me2 and non-CG DNA methylation pathways, for example, via mutation of the H3K9 methyltransferase genes KYP/SUVH4 SUVH5 SUVH6, or the CHG DNA methyltransferase gene CMT3, increases meiotic recombination in proximity to the centromeres. Using immunocytological detection of MLH1 foci and genotyping by sequencing of recombinant plants, we observe that H3K9me2 and non-CG DNA methylation pathway mutants show increased pericentromeric crossovers. Increased pericentromeric recombination in H3K9me2/non-CG mutants occurs in hybrid and inbred backgrounds and likely involves contributions from both the interfering and noninterfering crossover repair pathways. We also show that meiotic DNA double-strand breaks (DSBs) increase in H3K9me2/non-CG mutants within the pericentromeres, via purification and sequencing of SPO11-1-oligonucleotides. Therefore, H3K9me2 and non-CG DNA methylation exert a repressive effect on both meiotic DSB and crossover formation in plant pericentromeric heterochromatin. Our results may account for selection of enhancer trap Dissociation (Ds) transposons into the CMT3 gene by recombination with proximal transposon launch-pads.

Nucleosomes and DNA methylation shape meiotic DSB frequency in Arabidopsis thaliana transposons and gene regulatory regions

Meiotic recombination initiates from DNA double-strand breaks (DSBs) generated by SPO11 topoisomerase-like complexes. Meiotic DSB frequency varies extensively along eukaryotic chromosomes, with hotspots controlled by chromatin and DNA sequence. To map meiotic DSBs throughout a plant genome, we purified and sequenced Arabidopsis thaliana SPO11-1-oligonucleotides. SPO11-1-oligos are elevated in gene promoters, terminators, and introns, which is driven by AT-sequence richness that excludes nucleosomes and allows SPO11-1 access. A positive relationship was observed between SPO11-1-oligos and crossovers genome-wide, although fine-scale correlations were weaker. This may reflect the influence of interhomolog polymorphism on crossover formation, downstream from DSB formation. Although H3K4me3 is enriched in proximity to SPO11-1-oligo hotspots at gene 5' ends, H3K4me3 levels do not correlate with DSBs. Repetitive transposons are thought to be recombination silenced during meiosis, to prevent nonallelic interactions and genome instability. Unexpectedly, we found high SPO11-1-oligo levels in nucleosome-depleted Helitron/Pogo/Tc1/Mariner DNA transposons, whereas retrotransposons were coldspots. High SPO11-1-oligo transposons are enriched within gene regulatory regions and in proximity to immunity genes, suggesting a role as recombination enhancers. As transposon mobility in plant genomes is restricted by DNA methylation, we used the met1 DNA methyltransferase mutant to investigate the role of heterochromatin in SPO11-1-oligo distributions. Epigenetic activation of meiotic DSBs in proximity to centromeres and transposons occurred in met1 mutants, coincident with reduced nucleosome occupancy, gain of transcription, and H3K4me3. Together, our work reveals a complex relationship between chromatin and meiotic DSBs within A. thaliana genes and transposons, with significance for the diversity and evolution of plant genomes.

Unleashing meiotic crossovers in hybrid plants

Meiotic crossovers shuffle parental genetic information, providing novel combinations of alleles on which natural or artificial selection can act. However, crossover events are relatively rare, typically one to three exchange points per chromosome pair. Recent work has identified three pathways limiting meiotic crossovers in Arabidopsis thaliana that rely on the activity of FANCM [Crismani W, et al. (2012) Science 336:1588-1590], RECQ4 [Séguéla-Arnaud M, et al. (2015) Proc Natl Acad Sci USA 112:4713-4718], and FIGL1 [Girard C, et al. (2015) PLoS Genet 11:e1005369]. Here we analyzed recombination in plants in which one, two, or all three of these pathways were disrupted in both pure line and hybrid contexts. The greatest effect was observed when combining recq4 and figl1 mutations, which increased the hybrid genetic map length from 389 to 3,037 cM. This corresponds to an unprecedented 7.8-fold increase in crossover frequency. Disrupting the three pathways did not further increase recombination, suggesting that some upper limit had been reached. The increase in crossovers is not uniform along chromosomes and rises from centromere to telomere. Finally, although in wild type recombination is much higher in male meiosis than in female meiosis (490 cM vs. 290 cM), female recombination is higher than male recombination in recq4 figl1 (3,200 cM vs. 2,720 cM), suggesting that the factors that make wild-type female meiosis less recombinogenic than male wild-type meiosis do not apply in the mutant context. The massive increase in recombination observed in recq4 figl1 hybrids opens the possibility of manipulating recombination to enhance plant breeding efficiency.

Unleashing meiotic crossovers in hybrid plants

Meiotic crossovers shuffle parental genetic information, providing novel combinations of alleles on which natural or artificial selection can act. However, crossover events are relatively rare, typically one to three exchange points per chromosome pair. Recent work has identified three pathways limiting meiotic crossovers in Arabidopsis thaliana that rely on the activity of FANCM [Crismani W, et al. (2012) Science 336:1588-1590], RECQ4 [Séguéla-Arnaud M, et al. (2015) Proc Natl Acad Sci USA 112:4713-4718], and FIGL1 [Girard C, et al. (2015) PLoS Genet 11:e1005369]. Here we analyzed recombination in plants in which one, two, or all three of these pathways were disrupted in both pure line and hybrid contexts. The greatest effect was observed when combining recq4 and figl1 mutations, which increased the hybrid genetic map length from 389 to 3,037 cM. This corresponds to an unprecedented 7.8-fold increase in crossover frequency. Disrupting the three pathways did not further increase recombination, suggesting that some upper limit had been reached. The increase in crossovers is not uniform along chromosomes and rises from centromere to telomere. Finally, although in wild type recombination is much higher in male meiosis than in female meiosis (490 cM vs. 290 cM), female recombination is higher than male recombination in recq4 figl1 (3,200 cM vs. 2,720 cM), suggesting that the factors that make wild-type female meiosis less recombinogenic than male wild-type meiosis do not apply in the mutant context. The massive increase in recombination observed in recq4 figl1 hybrids opens the possibility of manipulating recombination to enhance plant breeding efficiency.

Massive crossover elevation via combination of HEI10 and recq4a recq4b during Arabidopsis meiosis

During meiosis, homologous chromosomes undergo reciprocal crossovers, which generate genetic diversity and underpin classical crop improvement. Meiotic recombination initiates from DNA double-strand breaks (DSBs), which are processed into single-stranded DNA that can invade a homologous chromosome. The resulting joint molecules can ultimately be resolved as crossovers. In Arabidopsis, competing pathways balance the repair of ∼100-200 meiotic DSBs into ∼10 crossovers per meiosis, with the excess DSBs repaired as noncrossovers. To bias DSB repair toward crossovers, we simultaneously increased dosage of the procrossover E3 ligase gene HEI10 and introduced mutations in the anticrossovers helicase genes RECQ4A and RECQ4B As HEI10 and recq4a recq4b increase interfering and noninterfering crossover pathways, respectively, they combine additively to yield a massive meiotic recombination increase. Interestingly, we also show that increased HEI10 dosage increases crossover coincidence, which indicates an effect on interference. We also show that patterns of interhomolog polymorphism and heterochromatin drive recombination increases distally towards the subtelomeres in both HEI10 and recq4a recq4b backgrounds, while the centromeres remain crossover suppressed. These results provide a genetic framework for engineering meiotic recombination landscapes in plant genomes.

Tackling Plant Meiosis: From Model Research to Crop Improvement

Genetic engineering and traditional plant breeding, which harnesses the natural genetic variation that arises during meiosis, will have key roles to improve crop varieties and thus deliver Food Security in the future. Meiosis, a specialized cell division producing haploid gametes to maintain somatic diploidy following their fusion, assures genetic variation by regulated genetic exchange through homologous recombination. However, meiotic recombination events are restricted in their total number and their distribution along chromosomes limiting allelic variations in breeding programs. Thus, modifying the number and distribution of meiotic recombination events has great potential to improve and accelerate plant breeding. In recent years much progress has been made in understanding meiotic progression and recombination in plants. Many genes and factors involved in these processes have been identified primarily in Arabidopsis thaliana but also more recently in crops such as Brassica, rice, barley, maize, or wheat. These advances put researchers in the position to translate acquired knowledge to various crops likely improving and accelerating breeding programs. However, although fundamental aspects of meiotic progression and recombination are conserved between species, differences in genome size and organization (due to repetitive DNA content and ploidy level) exist, particularly among plants, that likely account for differences in meiotic progression and recombination patterns found between species. Thus, tools and approaches are needed to better understand differences and similarities in meiotic progression and recombination among plants, to study fundamental aspects of meiosis in a variety of plants including crops and non-model species, and to transfer knowledge into crop species. In this article, we provide an overview of tools and approaches available to study plant meiosis, highlight new techniques, give examples of areas of future research and review distinct aspects of meiosis in non-model species.

MeioSeed: a CellProfiler-based program to count fluorescent seeds for crossover frequency analysis in Arabidopsis thaliana

BACKGROUND

The formation of crossovers during meiosis is pivotal for the redistribution of traits among the progeny of sexually reproducing organisms. In plants the molecular mechanisms underlying the formation of crossovers have been well established, but relatively little is known about the factors that determine the exact location and the frequency of crossover events in the genome. In the model plant species Arabidopsis, research on these factors has been greatly facilitated by reporter lines containing linked fluorescence marker genes under control of promoters active in seeds or pollen, allowing for the visualization of crossover events by fluorescence microscopy. However, the usefulness of these reporter lines to screen for novel modulators of crossover frequency in a high throughput manner relies on the availability of programs that can accurately count fluorescent seeds. Such a program was previously not available in scientific literature.

RESULTS

Here we present MeioSeed, a novel CellProfiler-based program that accurately counts GFP and RFP fluorescent Arabidopsis seeds with adjustable detection thresholds for fluorescence intensity, making use of a robust seed classifier which was trained by machine learning in Ilastik. Using the previously published reporter line Col3-4/20 as an example, we explain the use of MeioSeed and the steps taken to optimize the thresholding settings of the program to fit the published model for recombination frequency and transgene segregation. The use of MeioSeed is illustrated by investigating salt stress as a novel abiotic trigger for changes in crossover frequency in Col3-4/20 (♂) × Ler-0 (♀) F₁ hybrids. Salt stress was found to trigger increases in crossover frequency between the marker genes of up to 70% compared to the control treatment without salt stress. Genotyping of control and salt treated populations revealed that the changes in crossover frequency were not limited to the region between the marker genes, but that fluctuations in crossover frequency are likely to occur genome-wide after treatment with high salt concentrations.

CONCLUSIONS

MeioSeed allows for the high throughput recognition and counting of fluorescent Arabidopsis seeds and can facilitate the screening for novel abiotic and biotic modulators of crossover frequency using reporter lines in Arabidopsis.

Advances towards Controlling Meiotic Recombination for Plant Breeding

Meiotic homologous recombination generates new combinations of preexisting genetic variation and is a crucial process in plant breeding. Within the last decade, our understanding of plant meiotic recombination and genome diversity has advanced considerably. Innovation in DNA sequencing technology has led to the exploration of high-resolution genetic and epigenetic information in plant genomes, which has helped to accelerate plant breeding practices via high-throughput genotyping, and linkage and association mapping. In addition, great advances toward understanding the genetic and epigenetic control mechanisms of meiotic recombination have enabled the expansion of breeding programs and the unlocking of genetic diversity that can be used for crop improvement. This review highlights the recent literature on plant meiotic recombination and discusses the translation of this knowledge to the manipulation of meiotic recombination frequency and location with regards to crop plant breeding.

Modification of meiotic recombination by natural variation in plants

Meiosis is a specialized cell division that produces haploid gametes required for sexual reproduction. During the first meiotic division, homologous chromosomes pair and undergo reciprocal crossing over, which recombines linked sequence variation. Meiotic recombination frequency varies extensively both within and between species. In this review, we will examine the molecular basis of meiotic recombination rate variation, with an emphasis on plant genomes. We first consider cis modification caused by polymorphisms at the site of recombination, or elsewhere on the same chromosome. We review cis effects caused by mismatches within recombining joint molecules, the effect of structural hemizygosity, and the role of specific DNA sequence motifs. In contrast, trans modification of recombination is exerted by polymorphic loci encoding diffusible molecules, which are able to modulate recombination on the same and/or other chromosomes. We consider trans modifiers that act to change total recombination levels, hotspot locations, or interactions between homologous and homeologous chromosomes in polyploid species. Finally, we consider the significance of genetic variation that modifies meiotic recombination for adaptation and evolution of plant species.

A transposon-directed epigenetic change in ZmCCT underlies quantitative resistance to Gibberella stalk rot in maize

A major resistance quantitative trait locus, qRfg1, significantly enhances maize resistance to Gibberella stalk rot, a devastating disease caused by Fusarium graminearum. However, the underlying molecular mechanism remains unknown. We adopted a map-based cloning approach to identify the resistance gene at qRfg1 and examined the dynamic epigenetic changes during qRfg1-mediated maize resistance to the disease. A CCT domain-containing gene, ZmCCT, is the causal gene at the qRfg1 locus and a polymorphic CACTA-like transposable element (TE1) c. 2.4 kb upstream of ZmCCT is the genetic determinant of allelic variation. The non-TE1 ZmCCT allele is in a poised state, with predictive bivalent chromatin enriched for both repressive (H3K27me3/H3K9me3) and active (H3K4me3) histone marks. Upon pathogen challenge, this non-TE1 ZmCCT allele was promptly induced by a rapid yet transient reduction in H3K27me3/H3K9me3 and a progressive decrease in H3K4me3, leading to disease resistance. However, TE1 insertion in ZmCCT caused selective depletion of H3K4me3 and enrichment of methylated GC to suppress the pathogen-induced ZmCCT expression, resulting in disease susceptibility. Moreover, ZmCCT-mediated resistance to Gibberella stalk rot is not affected by photoperiod sensitivity. This chromatin-based regulatory mechanism enables ZmCCT to be more precise and timely in defense against F. graminearum infection.

Evidence for Adaptive Introgression of Disease Resistance Genes Among Closely Related Arabidopsis Species

The generation and maintenance of functional variation in the pathogen defense system of plants is central to the constant evolutionary battle between hosts and parasites. If a species is susceptible to a given pathogen, hybridization and subsequent introgression of a resistance allele from a related species can potentially be an important source of new immunity and is therefore expected to be selected for in a process referred to as adaptive introgression. Here, we survey sequence variation in 10 resistance (R-) genes and compare them with 37 reference genes in natural populations of the two closely related and interfertile species: Arabidopsis lyrata and A. halleri. The R-genes are highly polymorphic in both species and show clear signs of trans-species polymorphisms. We show that A. lyrata and A. halleri have had a history of limited introgression for the reference genes. For the R-genes, the introgression rate has been significantly higher than for the reference genes, resulting in fewer fixed differences between species and a higher sharing of identical haplotypes. We conclude that R-genes likely cross the species boundaries at a higher rate than reference genes and therefore also that some of the increased diversity and trans-specific polymorphisms in R-genes is due to adaptive introgression.