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

Homoeologous crossovers are distally biased and underlie genomic instability in first-generation neo-allopolyploid Arabidopsis suecica

First-generation polyploids often suffer from more meiotic errors and lower fertility than established wild polyploid populations. One such example is the allopolyploid model species Arabidopsis suecica which originated c. 16 000 generations ago. We present here a comparison of meiosis and its outcomes in naturally evolved and first-generation 'synthetic' A. suecica using a combination of cytological and genomic approaches. We show that while meiosis in natural lines is largely diploid-like, synthetic lines have high levels of meiotic errors including incomplete synapsis and nonhomologous crossover formation. Whole-genome re-sequencing of progeny revealed 20-fold higher levels of homoeologous exchange and eightfold higher aneuploidy originating from synthetic parents. Homoeologous exchanges showed a strong distal bias and occurred predominantly in genes, regularly generating novel protein variants. We also observed that homoeologous exchanges can generate megabase scale INDELs when occurring in regions of inverted synteny. Finally, we observed evidence of sex-specific differences in adaptation to polyploidy with higher success in reciprocal crosses to natural lines when synthetic plants were used as the female parent. Our results directly link cytological phenotypes in A. suecica with their genomic outcomes, demonstrating that homoeologous crossovers underlie genomic instability in neo-allopolyploids and are more distally biased than homologous crossovers.

Alternating between even and odd ploidy levels switches on and off the recombination control, even near the centromeres

Meiotic recombination is a key biological process in plant evolution and breeding, as it generates genetic diversity in each generation through the formation of crossovers (COs). However, due to their importance in genome stability, COs are highly regulated in frequency and distribution. We previously demonstrated that this strict regulation of COs can be modified, both in terms of CO frequency and distribution, in allotriploid Brassica hybrids (2n = 3x = 29; AAC) resulting from a cross between Brassica napus (2n = 4x = 38; AACC) and Brassica rapa (2n = 2x = 20; AA). Using the recently updated B. napus genome now including pericentromeres, we demonstrated that COs occur in these cold regions in allotriploids, as close as 375 kb from the centromere. Reverse transcription quantitative PCR (RT-qPCR) of various meiotic genes indicated that Class I COs are likely involved in the increased recombination frequency observed in allotriploids. We also demonstrated that this modified recombination landscape can be maintained via successive generations of allotriploidy (odd ploidy level). This deregulated meiotic behavior reverts to strict regulation in allotetraploid (even ploidy level) progeny in the second generation. Overall, we provide an easy way to manipulate tight recombination control in a polyploid crop.

Enhanced recombination empowers the detection and mapping of Quantitative Trait Loci

Modern plant breeding, such as genomic selection and gene editing, is based on the knowledge of the genetic architecture of desired traits. Quantitative trait loci (QTL) analysis, which combines high throughput phenotyping and genotyping of segregating populations, is a powerful tool to identify these genetic determinants and to decipher the underlying mechanisms. However, meiotic recombination, which shuffles genetic information between generations, is limited: Typically only one to two exchange points, called crossovers, occur between a pair of homologous chromosomes. Here we test the effect on QTL analysis of boosting recombination, by mutating the anti-crossover factors RECQ4 and FIGL1 in Arabidopsis thaliana full hybrids and lines in which a single chromosome is hybrid. We show that increasing recombination ~6-fold empowers the detection and resolution of QTLs, reaching the gene scale with only a few hundred plants. Further, enhanced recombination unmasks some secondary QTLs undetected under normal recombination. These results show the benefits of enhanced recombination to decipher the genetic bases of traits.

A catalogue of recombination coldspots in interspecific tomato hybrids

Increasing natural resistance and resilience in plants is key for ensuring food security within a changing climate. Breeders improve these traits by crossing cultivars with their wild relatives and introgressing specific alleles through meiotic recombination. However, some genomic regions are devoid of recombination especially in crosses between divergent genomes, limiting the combinations of desirable alleles. Here, we used pooled-pollen sequencing to build a map of recombinant and non-recombinant regions between tomato and five wild relatives commonly used for introgressive tomato breeding. We detected hybrid-specific recombination coldspots that underscore the role of structural variations in modifying recombination patterns and maintaining genetic linkage in interspecific crosses. Crossover regions and coldspots show strong association with specific TE superfamilies exhibiting differentially accessible chromatin between somatic and meiotic cells. About two-thirds of the genome are conserved coldspots, located mostly in the pericentromeres and enriched with retrotransposons. The coldspots also harbor genes associated with agronomic traits and stress resistance, revealing undesired consequences of linkage drag and possible barriers to breeding. We presented examples of linkage drag that can potentially be resolved by pairing tomato with other wild species. Overall, this catalogue will help breeders better understand crossover localization and make informed decisions on generating new tomato varieties.

Pushed waves, trailing edges, and extreme events: Eco-evolutionary dynamics of a geographic range shift in the owl limpet, Lottia gigantea

As climatic variation re-shapes global biodiversity, understanding eco-evolutionary feedbacks during species range shifts is of increasing importance. Theory on range expansions distinguishes between two different forms: "pulled" and "pushed" waves. Pulled waves occur when the source of the expansion comes from low-density peripheral populations, while pushed waves occur when recruitment to the expanding edge is supplied by high-density populations closer to the species' core. How extreme events shape pushed/pulled wave expansion events, as well as trailing-edge declines/contractions, remains largely unexplored. We examined eco-evolutionary responses of a marine invertebrate (the owl limpet, Lottia gigantea) that increased in abundance during the 2014-2016 marine heatwaves near the poleward edge of its geographic range in the northeastern Pacific. We used whole-genome sequencing from 19 populations across >11 degrees of latitude to characterize genomic variation, gene flow, and demographic histories across the species' range. We estimated present-day dispersal potential and past climatic stability to identify how contemporary and historical seascape features shape genomic characteristics. Consistent with expectations of a pushed wave, we found little genomic differentiation between core and leading-edge populations, and higher genomic diversity at range edges. A large and well-mixed population in the northern edge of the species' range is likely a result of ocean current anomalies increasing larval settlement and high-dispersal potential across biogeographic boundaries. Trailing-edge populations have higher differentiation from core populations, possibly driven by local selection and limited gene flow, as well as high genomic diversity likely as a result of climatic stability during the Last Glacial Maximum. Our findings suggest that extreme events can drive poleward range expansions that carry the adaptive potential of core populations, while also cautioning that trailing-edge extirpations may threaten unique evolutionary variation. This work highlights the importance of understanding how both trailing and leading edges respond to global change and extreme events.

Meiotic double-strand break repair DNA synthesis tracts in Arabidopsis thaliana

We report here the successful labelling of meiotic prophase I DNA synthesis in the flowering plant, Arabidopsis thaliana. Incorporation of the thymidine analogue, EdU, enables visualisation of the footprints of recombinational repair of programmed meiotic DNA double-strand breaks (DSB), with ~400 discrete, SPO11-dependent, EdU-labelled chromosomal foci clearly visible at pachytene and later stages of meiosis. This number equates well with previous estimations of 200-300 DNA double-strand breaks per meiosis in Arabidopsis, confirming the power of this approach to detect the repair of most or all SPO11-dependent meiotic DSB repair events. The chromosomal distribution of these DNA-synthesis foci accords with that of early recombination markers and MLH1, which marks Class I crossover sites. Approximately 10 inter-homologue cross-overs (CO) have been shown to occur in each Arabidopsis male meiosis and, athough very probably under-estimated, an equivalent number of inter-homologue gene conversions (GC) have been described. Thus, at least 90% of meiotic recombination events, and very probably more, have not previously been accessible for analysis. Visual examination of the patterns of the foci on the synapsed pachytene chromosomes corresponds well with expectations from the different mechanisms of meiotic recombination and notably, no evidence for long Break-Induced Replication DNA synthesis tracts was found. Labelling of meiotic prophase I, SPO11-dependent DNA synthesis holds great promise for further understanding of the molecular mechanisms of meiotic recombination, at the heart of reproduction and evolution of eukaryotes.

HEI10 coarsening, chromatin and sequence polymorphism shape the plant meiotic recombination landscape

Meiosis is a conserved eukaryotic cell division that produces spores required for sexual reproduction. During meiosis, chromosomes pair and undergo programmed DNA double-strand breaks, followed by homologous repair that can result in reciprocal crossovers. Crossover formation is highly regulated with typically few events per homolog pair. Crossovers additionally show wider spacing than expected from uniformly random placement - defining the phenomenon of interference. In plants, the conserved HEI10 E3 ligase is initially loaded along meiotic chromosomes, before maturing into a small number of foci, corresponding to crossover locations. We review the coarsening model that explains these dynamics as a diffusion and aggregation process, resulting in approximately evenly spaced HEI10 foci. We review how underlying chromatin states, and the presence of interhomolog polymorphisms, shape the meiotic recombination landscape, in light of the coarsening model. Finally, we consider future directions to understand the control of meiotic recombination in plant genomes.

The structure, function, and evolution of plant centromeres

Centromeres are essential regions of eukaryotic chromosomes responsible for the formation of kinetochore complexes, which connect to spindle microtubules during cell division. Notably, although centromeres maintain a conserved function in chromosome segregation, the underlying DNA sequences are diverse both within and between species and are predominantly repetitive in nature. The repeat content of centromeres includes high-copy tandem repeats (satellites), and/or specific families of transposons. The functional region of the centromere is defined by loading of a specific histone 3 variant (CENH3), which nucleates the kinetochore and shows dynamic regulation. In many plants, the centromeres are composed of satellite repeat arrays that are densely DNA methylated and invaded by centrophilic retrotransposons. In some cases, the retrotransposons become the sites of CENH3 loading. We review the structure of plant centromeres, including monocentric, holocentric, and metapolycentric architectures, which vary in the number and distribution of kinetochore attachment sites along chromosomes. We discuss how variation in CENH3 loading can drive genome elimination during early cell divisions of plant embryogenesis. We review how epigenetic state may influence centromere identity and discuss evolutionary models that seek to explain the paradoxically rapid change of centromere sequences observed across species, including the potential roles of recombination. We outline putative modes of selection that could act within the centromeres, as well as the role of repeats in driving cycles of centromere evolution. Although our primary focus is on plant genomes, we draw comparisons with animal and fungal centromeres to derive a eukaryote-wide perspective of centromere structure and function.

Methylomes as key features for predicting recombination in some plant species

Knowing how chromosome recombination works is essential for plant breeding. It enables the design of crosses between different varieties to combine desirable traits and create new ones. This is because the meiotic crossovers between homologous chromatids are not purely random, and various strategies have been developed to describe and predict such exchange events. Recent studies have used methylation data to predict chromosomal recombination in rice using machine learning models. This approach proved successful due to the presence of a positive correlation between the CHH context cytosine methylation and recombination rates in rice chromosomes. This paper assesses the question if methylation can be used to predict recombination in four plant species: Arabidopsis, maize, sorghum, and tomato. The results indicate a positive association between CHH context methylation and recombination rates in certain plant species, with varying degrees of strength in their relationships. The CG and CHG methylation contexts show negative correlation with recombination. Methylation data was key effectively in predicting recombination in sorghum and tomato, with a mean determination coefficient of 0.65 ± 0.11 and 0.76 ± 0.05, respectively. In addition, the mean correlation values between predicted and experimental recombination rates were 0.83 ± 0.06 for sorghum and 0.90 ± 0.05 for tomato, confirming the significance of methylomes in both monocotyledonous and dicotyledonous species. The predictions for Arabidopsis and maize were not as accurate, likely due to the comparatively weaker relationships between methylation contexts and recombination, in contrast to sorghum and tomato, where stronger associations were observed. To enhance the accuracy of predictions, further evaluations using data sets closely related to each other might prove beneficial. In general, this methylome-based method holds great potential as a reliable strategy for predicting recombination rates in various plant species, offering valuable insights to breeders in their quest to develop novel and improved varieties.

Investigating structural variant, indel and single nucleotide polymorphism differentiation between locally adapted Atlantic salmon populations

Genomic structural variants (SVs) are now recognized as an integral component of intraspecific polymorphism and are known to contribute to evolutionary processes in various organisms. However, they are inherently difficult to detect and genotype from readily available short-read sequencing data, and therefore remain poorly documented in wild populations. Salmonid species displaying strong interpopulation variability in both life history traits and habitat characteristics, such as Atlantic salmon (Salmo salar), offer a prime context for studying adaptive polymorphism, but the contribution of SVs to fine-scale local adaptation has yet to be explored. Here, we performed a comparative analysis of SVs, single nucleotide polymorphisms (SNPs) and small indels (<50 bp) segregating in the Romaine and Puyjalon salmon, two putatively locally adapted populations inhabiting neighboring rivers (Québec, Canada) and showing pronounced variation in life history traits, namely growth, fecundity, and age at maturity and smoltification. We first catalogued polymorphism using a hybrid SV characterization approach pairing both short- (16X) and long-read sequencing (20X) for variant discovery with graph-based genotyping of SVs across 60 salmon genomes, along with characterization of SNPs and small indels from short reads. We thus identified 115,907 SVs, 8,777,832 SNPs and 1,089,321 short indels, with SVs covering 4.8 times more base pairs than SNPs. All three variant types revealed a highly congruent population structure and similar patterns of F ST and density variation along the genome. Finally, we performed outlier detection and redundancy analysis (RDA) to identify variants of interest in the putative local adaptation of Romaine and Puyjalon salmon. Genes located near these variants were enriched for biological processes related to nervous system function, suggesting that observed variation in traits such as age at smoltification could arise from differences in neural development. This study therefore demonstrates the feasibility of large-scale SV characterization and highlights its relevance for salmonid population genomics.

Structural Variants and Speciation: Multiple Processes at Play

Research on the genomic architecture of speciation has increasingly revealed the importance of structural variants (SVs) that affect the presence, abundance, position, and/or direction of a nucleotide sequence. SVs include large chromosomal rearrangements such as fusion/fissions and inversions and translocations, as well as smaller variants such as duplications, insertions, and deletions (CNVs). Although we have ample evidence that SVs play a key role in speciation, the underlying mechanisms differ depending on the type and length of the SV, as well as the ecological, demographic, and historical context. We review predictions and empirical evidence for classic processes such as underdominance due to meiotic aberrations and the coupling effect of recombination suppression before exploring how recent sequencing methodologies illuminate the prevalence and diversity of SVs. We discuss specific properties of SVs and their impact throughout the genome, highlighting that multiple processes are at play, and possibly interacting, in the relationship between SVs and speciation.

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.

Epigenetic regulation during meiosis and crossover

Meiosis is a distinctive type of cell division that reorganizes genetic material between generations. The initial stages of meiosis consist of several crucial steps which include double strand break, homologous chromosome pairing, break repair and crossover. Crossover frequency varies depending on the position on the chromosome, higher at euchromatin region and rare at heterochromatin, centromeres, telomeres and ribosomal DNA. Crossover positioning is dependent on various factors, especially epigenetic modifications. DNA methylation, histone post-translational modifications, histone variants and non-coding RNAs are most probably playing an important role in positioning of crossovers on a chromosomal level as well as hotspot level. DNA methylation negatively regulates crossover frequency and its effect is visible in centromeres, pericentromeres and heterochromatin regions. Pericentromeric chromatin and heterochromatin mark studies have been a centre of attraction in meiosis. Crossover hotspots are associated with euchromatin regions having specific chromatin modifications such as H3K4me3, H2A.Z. and H3 acetylation. This review will provide the current understanding of the epigenetic role in plants during meiotic recombination, chromosome synapsis, double strand break and hotspots with special attention to euchromatin and heterochromatin marks. Further, the role of epigenetic modifications in regulating meiosis and crossover in other organisms is also discussed.

Comprehensive dissection of meiotic DNA double-strand breaks and crossovers in cucumber

Meiotic recombination drives genetic diversity and crop genome optimization. In plant breeding, parents with favorable traits are crossed to create elite varieties. Different hybridizations produce diverse types of segment reshuffling between homologous chromosomes. However, little is known about the factors that cause hybrid-specific changes in crossovers (COs). Here, we constructed 2 F2 populations from crosses between a semiwild and 2 domesticated cucumber (Cucumis sativus) accessions and examined CO events. COs mainly occurred around genes and differed unevenly along chromosomes between the 2 hybrids. Fine-scale CO distributions were suppressed in regions of heterozygous structural variations (SVs) and were accelerated by high sequence polymorphism. C. sativus RADiation sensitive 51A (CsRAD51A) binding, histone H3 lysine 4 trimethylation (H3K4me3) modification, chromatin accessibility, and hypomethylation were positively associated with global CO landscapes and in local DNA double-strand break (DSB) hotspots and genes. The frequency and suppression of COs could be roughly predicted based on multiomic information. Differences in CO events between hybrids could be partially traced to distinct genetic and epigenetic features and were significantly associated with specific DSB hotspots and heterozygous SVs. Our findings identify the genomic and epigenetic features that contribute to CO formation and hybrid-specific divergence in cucumber and provide theoretical support for selecting parental combinations and manipulating recombination events at target genomic regions during plant breeding.

MSH2 stimulates interfering and inhibits non-interfering crossovers in response to genetic polymorphism

Meiotic crossovers can be formed through the interfering pathway, in which one crossover prevents another from forming nearby, or by an independent non-interfering pathway. In Arabidopsis, local sequence polymorphism between homologs can stimulate interfering crossovers in a MSH2-dependent manner. To understand how MSH2 regulates crossovers formed by the two pathways, we combined Arabidopsis mutants that elevate non-interfering crossovers with msh2 mutants. We demonstrate that MSH2 blocks non-interfering crossovers at polymorphic loci, which is the opposite effect to interfering crossovers. We also observe MSH2-independent crossover inhibition at highly polymorphic sites. We measure recombination along the chromosome arms in lines differing in patterns of heterozygosity and observe a MSH2-dependent crossover increase at the boundaries between heterozygous and homozygous regions. Here, we show that MSH2 is a master regulator of meiotic DSB repair in Arabidopsis, with antagonistic effects on interfering and non-interfering crossovers, which shapes the crossover landscape in relation to interhomolog polymorphism.

Meiotic recombination is confirmed to be unusually high in the fission yeast Schizosaccharomyces pombe

In most eukaryotes, meiotic crossovers (COs) are limited to 1-3 per chromosome, and are prevented from occurring close to one another by CO interference. The fission yeast Schizosaccharomyces pombe, an exception to these general rules, was reported to have the highest CO number per chromosome and no or weak interference. However, global CO frequency was indirectly estimated, calling for confirmation. Here, we used an innovative strategy to determine COs genome-wide in S. pombe. We confirmed weak CO interference, acting at physical distances compatible with the patterning of recombination precursors. We revealed a slight co-variation in CO number between chromosomes, suggesting that a limiting pro-CO factor varies between meiocytes. CO number per chromosome varies proportionally with chromosome size, with the three chromosomes having, on average, 15.9, 12.5, and 7.0 COs, respectively. This reinforces S. pombe's status as the eukaryote with the highest CO number per chromosome described to date.

Transition to Self-compatibility Associated With Dominant S-allele in a Diploid Siberian Progenitor of Allotetraploid Arabidopsis kamchatica Revealed by Arabidopsis lyrata Genomes

A transition to selfing can be beneficial when mating partners are scarce, for example, due to ploidy changes or at species range edges. Here, we explain how self-compatibility evolved in diploid Siberian Arabidopsis lyrata, and how it contributed to the establishment of allotetraploid Arabidopsis kamchatica. First, we provide chromosome-level genome assemblies for two self-fertilizing diploid A. lyrata accessions, one from North America and one from Siberia, including a fully assembled S-locus for the latter. We then propose a sequence of events leading to the loss of self-incompatibility in Siberian A. lyrata, date this independent transition to ∼90 Kya, and infer evolutionary relationships between Siberian and North American A. lyrata, showing an independent transition to selfing in Siberia. Finally, we provide evidence that this selfing Siberian A. lyrata lineage contributed to the formation of the allotetraploid A. kamchatica and propose that the selfing of the latter is mediated by the loss-of-function mutation in a dominant S-allele inherited from A. lyrata.

Genomic features of meiotic crossovers in diploid potato

Meiotic recombination plays an important role in genome evolution and crop improvement. Potato (Solanum tuberosum L.) is the most important tuber crop in the world, but research about meiotic recombination in potato is limited. Here, we resequenced 2163 F₂ clones derived from five different genetic backgrounds and identified 41 945 meiotic crossovers. Some recombination suppression in euchromatin regions was associated with large structural variants. We also detected five shared crossover hotspots. The number of crossovers in each F₂ individual from the accession Upotato 1 varied from 9 to 27, with an average of 15.5, 78.25% of which were mapped within 5 kb of their presumed location. We show that 57.1% of the crossovers occurred in gene regions, with poly-A/T, poly-AG, AT-rich, and CCN repeats enriched in the crossover intervals. The recombination rate is positively related with gene density, SNP density, Class II transposon, and negatively related with GC density, repeat sequence density and Class I transposon. This study deepens our understanding of meiotic crossovers in potato and provides useful information for diploid potato breeding.

The regulation of meiotic crossover distribution: a coarse solution to a century-old mystery?

Meiotic crossovers, which are exchanges of genetic material between homologous chromosomes, are more evenly and distantly spaced along chromosomes than expected by chance. This is because the occurrence of one crossover reduces the likelihood of nearby crossover events - a conserved and intriguing phenomenon called crossover interference. Although crossover interference was first described over a century ago, the mechanism allowing coordination of the fate of potential crossover sites half a chromosome away remains elusive. In this review, we discuss the recently published evidence supporting a new model for crossover patterning, coined the coarsening model, and point out the missing pieces that are still needed to complete this fascinating puzzle.

Heterozygous inversion breakpoints suppress meiotic crossovers by altering recombination repair outcomes

Heterozygous chromosome inversions suppress meiotic crossover (CO) formation within an inversion, potentially because they lead to gross chromosome rearrangements that produce inviable gametes. Interestingly, COs are also severely reduced in regions nearby but outside of inversion breakpoints even though COs in these regions do not result in rearrangements. Our mechanistic understanding of why COs are suppressed outside of inversion breakpoints is limited by a lack of data on the frequency of noncrossover gene conversions (NCOGCs) in these regions. To address this critical gap, we mapped the location and frequency of rare CO and NCOGC events that occurred outside of the dl-49 chrX inversion in D. melanogaster. We created full-sibling wildtype and inversion stocks and recovered COs and NCOGCs in the syntenic regions of both stocks, allowing us to directly compare rates and distributions of recombination events. We show that COs outside of the proximal inversion breakpoint are distributed in a distance-dependent manner, with strongest suppression near the inversion breakpoint. We find that NCOGCs occur evenly throughout the chromosome and, importantly, are not suppressed near inversion breakpoints. We propose a model in which COs are suppressed by inversion breakpoints in a distance-dependent manner through mechanisms that influence DNA double-strand break repair outcome but not double-strand break formation. We suggest that subtle changes in the synaptonemal complex and chromosome pairing might lead to unstable interhomolog interactions during recombination that permits NCOGC formation but not CO formation.

Technology-driven approaches for meiosis research in tomato and wild relatives

Meiosis is a specialized cell division during reproduction where one round of chromosomal replication is followed by genetic recombination and two rounds of segregation to generate recombined, ploidy-reduced spores. Meiosis is crucial to the generation of new allelic combinations in natural populations and artificial breeding programs. Several plant species are used in meiosis research including the cultivated tomato (Solanum lycopersicum) which is a globally important crop species. Here we outline the unique combination of attributes that make tomato a powerful model system for meiosis research. These include the well-characterized behavior of chromosomes during tomato meiosis, readily available genomics resources, capacity for genome editing, clonal propagation techniques, lack of recent polyploidy and the possibility to generate hybrids with twelve related wild species. We propose that further exploitation of genome bioinformatics, genome editing and artificial intelligence in tomato will help advance the field of plant meiosis research. Ultimately this will help address emerging themes including the evolution of meiosis, how recombination landscapes are determined, and the effect of temperature on meiosis.

Fine mapping of meiotic crossovers in Brassica oleracea reveals patterns and variations depending on direction and combination of crosses

Meiotic recombination is crucial for assuring proper segregation of parental chromosomes and generation of novel allelic combinations. As this process is tightly regulated, identifying factors influencing rate, and distribution of meiotic crossovers (COs) is of major importance, notably for plant breeding programs. However, high-resolution recombination maps are sparse in most crops including the Brassica genus and knowledge about intraspecific variation and sex differences is lacking. Here, we report fine-scale resolution recombination landscapes for 10 female and 10 male crosses in Brassica oleracea, by analyzing progenies of five large four-way-cross populations from two reciprocally crossed F1s per population. Parents are highly diverse inbred lines representing major crops, including broccoli, cauliflower, cabbage, kohlrabi, and kale. We produced approximately 4.56T Illumina data from 1248 progenies and identified 15 353 CO across the 10 reciprocal crosses, 51.13% of which being mapped to <10 kb. We revealed fairly similar Mb-scale recombination landscapes among all cross combinations and between the sexes, and provided evidence that these landscapes are largely independent of sequence divergence. We evidenced strong influence of gene density and large structural variations on CO formation in B. oleracea. Moreover, we found extensive variations in CO number depending on the direction and combination of the initial parents crossed with, for the first time, a striking interdependency between these factors. These data improve our current knowledge on meiotic recombination and are important for Brassica breeders.

Global dissection of the recombination landscape in soybean using a high-density 600K SoySNP array

Recombination is crucial for crop breeding because it can break linkage drag and generate novel allele combinations. However, the high-resolution recombination landscape and its driving forces in soybean are largely unknown. Here, we constructed eight recombinant inbred line (RIL) populations and genotyped individual lines using the high-density 600K SoySNP array, which yielded a high-resolution recombination map with 5636 recombination sites at a resolution of 1.37 kb. The recombination rate was negatively correlated with transposable element density and GC content but positively correlated with gene density. Interestingly, we found that meiotic recombination was enriched at the promoters of active genes. Further investigations revealed that chromatin accessibility and active epigenetic modifications promoted recombination. Our findings provide important insights into the control of homologous recombination and thus will increase our ability to accelerate soybean breeding by manipulating meiotic recombination rate.

Crossover patterning in plants

KEY MESSAGE

Chromatin state, and dynamic loading of pro-crossover protein HEI10 at recombination intermediates shape meiotic chromosome patterning in plants. Meiosis is the basis of sexual reproduction, and its basic progression is conserved across eukaryote kingdoms. A key feature of meiosis is the formation of crossovers which result in the reciprocal exchange of segments of maternal and paternal chromosomes. This exchange generates chromosomes with new combinations of alleles, increasing the efficiency of both natural and artificial selection. Crossovers also form a physical link between homologous chromosomes at metaphase I which is critical for accurate chromosome segregation and fertility. The patterning of crossovers along the length of chromosomes is a highly regulated process, and our current understanding of its regulation forms the focus of this review. At the global scale, crossover patterning in plants is largely governed by the classically observed phenomena of crossover interference, crossover homeostasis and the obligatory crossover which regulate the total number of crossovers and their relative spacing. The molecular actors behind these phenomena have long remained obscure, but recent studies in plants implicate HEI10 and ZYP1 as key players in their coordination. In addition to these broad forces, a wealth of recent studies has highlighted how genomic and epigenomic features shape crossover formation at both chromosomal and local scales, revealing that crossovers are primarily located in open chromatin associated with gene promoters and terminators with low nucleosome occupancy.

High-resolution mapping reveals hotspots and sex-biased recombination in Populus trichocarpa

Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.

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.

Gamete Binning to Achieve Haplotype-Resolved Genome Assembly

Haplotype-resolved genome assemblies remain a challenge in practice. Here, we provide a step-by-step guide on gamete binning, a method to generate haplotype-resolved genome assemblies for diploid species. The protocol starts by phasing heterozygous variants to individual haplotypes of specific chromosomes using the genome information of individual haploid gametes of the focal individual. Using phased variants, the whole-genome sequencing reads from the diploid genome can be genotyped and assigned into groups, which represent the individual haplotypes of each of the chromosomes. Finally, haplotype-specific chromosomes can be assembled independently using standard assembly tools. First applications of gamete binning revealed a haplotyping accuracy over 99%, which outperformed sequence-only or Hi-C-based haplotype-resolved genome assemblies.Availability: github.com/schneebergerlab/GameteBinning_prac .

Three problems in the genetics of speciation by selection

Speciation is the process by which barriers to gene flow evolve between populations. Although we now know that speciation is largely driven by natural selection, knowledge of the agents of selection and the genetic and genomic mechanisms that facilitate divergence is required for a satisfactory theory of speciation. In this essay, we highlight three advances/problems in our understanding of speciation that have arisen from studies of the genes and genomic regions that underlie the evolution of reproductive isolation. First, we describe how the identification of “speciation” genes makes it possible to identify the agents of selection causing the evolution of reproductive isolation, while also noting that the link between the genetics of phenotypic divergence and intrinsic postzygotic reproductive barriers remains tenuous. Second, we discuss the important role of recombination suppressors in facilitating speciation with gene flow, but point out that the means and timing by which reproductive barriers become associated with recombination cold spots remains uncertain. Third, we establish the importance of ancient genetic variation in speciation, although we argue that the focus of speciation studies on evolutionarily young groups may bias conclusions in favor of ancient variation relative to new mutations.

Coexpression of MEIOTIC-TOPOISOMERASE VIB-dCas9 with guide RNAs specific to a recombination hotspot is insufficient to increase crossover frequency in Arabidopsis

During meiosis, homologous chromosomes pair and recombine, which can result in reciprocal crossovers that increase genetic diversity. Crossovers are unevenly distributed along eukaryote chromosomes and show repression in heterochromatin and the centromeres. Within the chromosome arms, crossovers are often concentrated in hotspots, which are typically in the kilobase range. The uneven distribution of crossovers along chromosomes, together with their low number per meiosis, creates a limitation during crop breeding, where recombination can be beneficial. Therefore, targeting crossovers to specific genome locations has the potential to accelerate crop improvement. In plants, meiotic crossovers are initiated by DNA double-strand breaks that are catalyzed by SPO11 complexes, which consist of 2 catalytic (SPO11-1 and SPO11-2) and 2 noncatalytic subunits (MTOPVIB). We used the model plant Arabidopsis thaliana to coexpress an MTOPVIB-dCas9 fusion protein with guide RNAs specific to the 3a crossover hotspot. We observed that this was insufficient to significantly change meiotic crossover frequency or pattern within 3a. We discuss the implications of our findings for targeting meiotic recombination within plant genomes.

The megabase-scale crossover landscape is largely independent of sequence divergence

Meiotic recombination frequency varies along chromosomes and strongly correlates with sequence divergence. However, the causal relationship between recombination landscapes and polymorphisms is unclear. Here, we characterize the genome-wide recombination landscape in the quasi-absence of polymorphisms, using Arabidopsis thaliana homozygous inbred lines in which a few hundred genetic markers were introduced through mutagenesis. We find that megabase-scale recombination landscapes in inbred lines are strikingly similar to the recombination landscapes in hybrids, with the notable exception of heterozygous large rearrangements where recombination is prevented locally. In addition, the megabase-scale recombination landscape can be largely explained by chromatin features. Our results show that polymorphisms are not a major determinant of the shape of the megabase-scale recombination landscape but rather favour alternative models in which recombination and chromatin shape sequence divergence across the genome.

The frequency of recombination varies along chromosomes and highly correlates with sequence divergence. Here, the authors show that polymorphisms are not a major determinant of the megabase-scale recombination landscape in Arabidopsis, which is rather determined by chromatin accessibility and DNA methylation.

Recombination landscape divergence between populations is marked by larger low-recombining regions in domesticated rye

The genomic landscape of recombination plays an essential role in evolution. Patterns of recombination are highly variable along chromosomes, between sexes, individuals, populations, and species. In many eukaryotes, recombination rates are elevated in sub-telomeric regions and drastically reduced near centromeres, resulting in large low-recombining (LR) regions. The processes of recombination are influenced by genetic factors, such as different alleles of genes involved in meiosis and chromatin structure, as well as external environmental stimuli like temperature and overall stress. In this work, we focused on the genomic landscapes of recombination in a collection of 916 rye (Secale cereale) individuals. By analysing population structure among individuals of different domestication status and geographic origin, we detected high levels of admixture, reflecting the reproductive biology of a self-incompatible, wind-pollinating grass species. We then analysed patterns of recombination in overlapping subpopulations, which revealed substantial variation in the physical size of LR regions, with a tendency for larger LR regions in domesticated subpopulations. Genome-wide association scans (GWAS) for LR region size revealed a major quantitative-trait-locus (QTL) at which, among 18 annotated genes, an ortholog of histone H4 acetyltransferase ESA1 was located. Rye individuals belonging to domesticated subpopulations showed increased synaptonemal complex length, but no difference in crossover frequency, indicating that only the recombination landscape is different. Furthermore, the genomic region harbouring rye ScESA1 showed moderate patterns of selection in domesticated subpopulations, suggesting that larger LR regions were indirectly selected for during domestication to achieve more homogeneous populations for agricultural use.

Allele segregation analysis of F1 hybrids between independent Brassica allohexaploid lineages

In the Brassica genus, we find both diploid species (one genome) and allotetraploid species (two different genomes) but no naturally occurring hexaploid species (three different genomes, AABBCC). Although hexaploids can be produced via human intervention, these neo-polyploids have quite unstable genomes and usually suffer from severe genome reshuffling. Whether these genome rearrangements continue in later generations and whether genomic arrangements follow similar, reproducible patterns between different lineages is still unknown. We crossed Brassica hexaploids resulting from different species combinations to produce five F₁ hybrids and analyzed the karyotypes of the parents and the F₁ hybrids, as well as allele segregation in a resulting test-cross population via molecular karyotyping using SNP array genotyping. Although some genomic regions were found to be more likely to be duplicated, deleted, or rearranged, a consensus pattern was not shared between genotypes. Brassica hexaploids had a high tolerance for fixed structural rearrangements, but which rearrangements occur and become fixed over many generations does not seem to show either strong reproducibility or to indicate selection for stability. On average, we observed 10 de novo chromosome rearrangements contributed almost equally from both parents to the F₁ hybrids. At the same time, the F₁ hybrid meiosis produced on average 8.6 new rearrangements. Hence, the increased heterozygosity in the F₁ hybrid did not significantly improve genome stability in our hexaploid hybrids and might have had the opposite effect. However, hybridization between lineages was readily achieved and may be exploited for future genetics and breeding purposes.

Epigenomic and structural events preclude recombination in Brassica napus

Meiotic recombination is a major evolutionary process generating genetic diversity at each generation in sexual organisms. However, this process is highly regulated, with the majority of crossovers lying in the distal chromosomal regions that harbor low DNA methylation levels. Even in these regions, some islands without recombination remain, for which we investigated the underlying causes. Genetic maps were established in two Brassica napus hybrids to detect the presence of such large nonrecombinant islands. The role played by DNA methylation and structural variations in this local absence of recombination was determined by performing bisulfite sequencing and whole genome comparisons. Inferred structural variations were validated using either optical mapping or oligo fluorescence in situ hybridization. Hypermethylated or inverted regions between Brassica genomes were associated with the absence of recombination. Pairwise comparisons of nine B. napus genome assemblies revealed that such inversions occur frequently and may contain key agronomic genes such as resistance to biotic stresses. We conclude that such islands without recombination can have different origins, such as DNA methylation or structural variations in B. napus. It is thus essential to take into account these features in breeding programs as they may hamper the efficient combination of favorable alleles in elite varieties.

Genomic prediction of the recombination rate variation in barley - A route to highly recombinogenic genotypes

Meiotic recombination is not only fundamental to the adaptation of sexually reproducing eukaryotes in nature but increased recombination rates facilitate the combination of favourable alleles into a single haplotype in breeding programmes. The main objectives of this study were to (i) assess the extent and distribution of the recombination rate variation in cultivated barley (Hordeum vulgare L.), (ii) quantify the importance of the general and specific recombination effects, and (iii) evaluate a genomic selection approach's ability to predict the recombination rate variation. Genetic maps were created for the 45 segregating populations that were derived from crosses among 23 spring barley inbreds with origins across the world. The genome-wide recombination rate among populations ranged from 0.31 to 0.73 cM/Mbp. The crossing design used in this study allowed to separate the general recombination effects (GRE) of individual parental inbreds from the specific recombination effects (SRE) caused by the combinations of parental inbreds. The variance of the genome-wide GRE was found to be about eight times the variance of the SRE. This finding indicated that parental inbreds differ in the efficiency of their recombination machinery. The ability to predict the chromosome or genome-wide recombination rate of an inbred ranged from 0.80 to 0.85. These results suggest that a reliable screening of large genetic materials for their potential to cause a high extent of genetic recombination in their progeny is possible, allowing to systematically manipulate the recombination rate using natural variation.

Male and female recombination landscapes of diploid Arabidopsis arenosa

The number and placement of meiotic crossover events during meiosis have important implications for the fidelity of chromosome segregation as well as patterns of inheritance. Despite the functional importance of recombination, recombination landscapes vary widely among and within species, and this can have a strong impact on evolutionary processes. A good knowledge of recombination landscapes is important for model systems in evolutionary and ecological genetics, since it can improve interpretation of genomic patterns of differentiation and genome evolution, and provides an important starting point for understanding the causes and consequences of recombination rate variation. Arabidopsis arenosa is a powerful evolutionary genetic model for studying the molecular basis of adaptation and recombination rate evolution. Here, we generate genetic maps for 2 diploid A. arenosa individuals from distinct genetic lineages where we have prior knowledge that meiotic genes show evidence of selection. We complement the genetic maps with cytological approaches to map and quantify recombination rates, and test the idea that these populations might have distinct patterns of recombination. We explore how recombination differs at the level of populations, individuals, sexes and genomic regions. We show that the positioning of crossovers along a chromosome correlates with their number, presumably a consequence of crossover interference, and discuss how this effect can cause differences in recombination landscape among sexes or species. We identify several instances of female segregation distortion. We found that averaged genome-wide recombination rate is lower and sex differences subtler in A. arenosa than in Arabidopsis thaliana.

Sexual dimorphism and the effect of wild introgressions on recombination in cassava (Manihot esculenta Crantz) breeding germplasm

Recombination has essential functions in meiosis, evolution, and breeding. The frequency and distribution of crossovers dictate the generation of new allele combinations and can vary across species and between sexes. Here, we examine recombination landscapes across the 18 chromosomes of cassava (Manihot esculenta Crantz) with respect to male and female meioses and known introgressions from the wild relative Manihot glaziovii. We used SHAPEIT2 and duoHMM to infer crossovers from genotyping-by-sequencing data and a validated multigenerational pedigree from the International Institute of Tropical Agriculture cassava breeding germplasm consisting of 7020 informative meioses. We then constructed new genetic maps and compared them to an existing map previously constructed by the International Cassava Genetic Map Consortium. We observed higher recombination rates in females compared to males, and lower recombination rates in M. glaziovii introgression segments on chromosomes 1 and 4, with suppressed recombination along the entire length of the chromosome in the case of the chromosome 4 introgression. Finally, we discuss hypothesized mechanisms underlying our observations of heterochiasmy and crossover suppression and discuss the broader implications for plant breeding.

High-Throughput Fluorescent Pollen Tetrad Analysis Using DeepTetrad

Meiotic recombination initiates from ~100-200 s of programmed DNA double stranded breaks (DSBs) in plants. Meiotic DSBs can be repaired using homologous chromosomes to generate a crossover . Meiotic crossover is critical for chromosomal segregation and increasing genetic variation. The number of crossovers is limited to one and three per chromosome pair in most plant species. Genetic, epigenetic, and environmental factors control crossover frequency and distribution. Due to the limited number of crossovers it is challenging to measure crossover frequency along chromosomes. We adapted fluorescence-tagged lines (FTLs ) that contain quartet1 mutations and linked transgenes expressing dsRed, eYFP, and eCFP in pollen tetrads into the deep learning-based image analysis tool, DeepTetrad. DeepTetrad enables the measurement of crossover frequency and interference by classifying 12 types of tetrads from three-color FTLs in a high-throughput manner, using conventional microscope instruments and a Linux machine. Here, we provide detailed procedures for preparing tetrad samples, tetrad imaging, running DeepTetrad, and analysis of DeepTetrad outputs. DeepTetrad-based measurements of crossover frequency and interference ratio will accelerate the genetic dissection of meiotic crossover control.

CTAB DNA Extraction and Genotyping-by-Sequencing to Map Meiotic Crossovers in Plants

Reciprocal DNA crossovers between chromosomes form new allelic combinations and contribute to the formation of novel genetic diversity. Crossovers are formed during meiosis of germ cells and these recombination events have influenced plant genome evolution and are used during breeding to create improved plant varieties. Meiotic crossovers are not uniformly formed across the genome but instead occur in regions with low nucleosome density. The recombination landscape differs between the model plant organism Arabidopsis thaliana and crops such as rice and maize. Genotyping-by-sequencing is a technique that can detect crossover location and provide information on the recombination landscape genome-wide. This technique can be used to compare crossover position between ecotypes, species, and mutant lines to gain information on factors controlling meiotic recombination. In this protocol, we describe the steps to purify DNA from plant tissue, prepare 96 DNA libraries in parallel and perform quality control before next-generation sequencing.

Quantitative modelling of fine-scale variations in the Arabidopsis thaliana crossover landscape

In, essentially, all species where meiotic crossovers (COs) have been studied, they occur preferentially in open chromatin, typically near gene promoters and to a lesser extent, at the end of genes. Here, in the case of Arabidopsis thaliana, we unveil further trends arising when one considers contextual information, namely summarised epigenetic status, gene or intergenic region size, and degree of divergence between homologs. For instance, we find that intergenic recombination rate is reduced if those regions are less than 1.5 kb in size. Furthermore, we propose that the presence of single nucleotide polymorphisms enhances the rate of CO formation compared to when homologous sequences are identical, in agreement with previous works comparing rates in adjacent homozygous and heterozygous blocks. Lastly, by integrating these different effects, we produce a quantitative and predictive model of the recombination landscape that reproduces much of the experimental variation.

Pan-genome of Raphanus highlights genetic variation and introgression among domesticated, wild, and weedy radishes

Post-polyploid diploidization associated with descending dysploidy and interspecific introgression drives plant genome evolution by unclear mechanisms. Raphanus is an economically and ecologically important Brassiceae genus and model system for studying post-polyploidization genome evolution and introgression. Here, we report the de novo sequence assemblies for 11 genomes covering most of the typical sub-species and varieties of domesticated, wild and weedy radishes from East Asia, South Asia, Europe, and America. Divergence among the species, sub-species, and South/East Asian types coincided with Quaternary glaciations. A genus-level pan-genome was constructed with family-based, locus-based, and graph-based methods, and whole-genome comparisons revealed genetic variations ranging from single-nucleotide polymorphisms (SNPs) to inversions and translocations of whole ancestral karyotype (AK) blocks. Extensive gene flow occurred between wild, weedy, and domesticated radishes. High frequencies of genome reshuffling, biased retention, and large-fragment translocation have shaped the genomic diversity. Most variety-specific gene-rich blocks showed large structural variations. Extensive translocation and tandem duplication of dispensable genes were revealed in two large rearrangement-rich islands. Disease resistance genes mostly resided on specific and dispensable loci. Variations causing the loss of function of enzymes modulating gibberellin deactivation were identified and could play an important role in phenotype divergence and adaptive evolution. This study provides new insights into the genomic evolution underlying post-polyploid diploidization and lays the foundation for genetic improvement of radish crops, biological control of weeds, and protection of wild species' germplasms.

Domestication shapes recombination patterns in tomato

Meiotic recombination is a biological process of key importance in breeding, to generate genetic diversity and develop novel or agronomically relevant haplotypes. In crop tomato, recombination is curtailed as manifested by linkage disequilibrium decay over a longer distance and reduced diversity compared with wild relatives. Here, we compared domesticated and wild populations of tomato and found an overall conserved recombination landscape, with local changes in effective recombination rate in specific genomic regions. We also studied the dynamics of recombination hotspots resulting from domestication and found that loss of such hotspots is associated with selective sweeps, most notably in the pericentromeric heterochromatin. We detected footprints of genetic changes and structural variants, among them associated with transposable elements, linked with hotspot divergence during domestication, likely causing fine-scale alterations to recombination patterns and resulting in linkage drag.

Selective sorting of ancestral introgression in maize and teosinte along an elevational cline

While often deleterious, hybridization can also be a key source of genetic variation and pre-adapted haplotypes, enabling rapid evolution and niche expansion. Here we evaluate these opposing selection forces on introgressed ancestry between maize (Zea mays ssp. mays) and its wild teosinte relative, mexicana (Zea mays ssp. mexicana). Introgression from ecologically diverse teosinte may have facilitated maize's global range expansion, in particular to challenging high elevation regions (> 1500 m). We generated low-coverage genome sequencing data for 348 maize and mexicana individuals to evaluate patterns of introgression in 14 sympatric population pairs, spanning the elevational range of mexicana, a teosinte endemic to the mountains of Mexico. While recent hybrids are commonly observed in sympatric populations and mexicana demonstrates fine-scale local adaptation, we find that the majority of mexicana ancestry tracts introgressed into maize over 1000 generations ago. This mexicana ancestry seems to have maintained much of its diversity and likely came from a common ancestral source, rather than contemporary sympatric populations, resulting in relatively low FST between mexicana ancestry tracts sampled from geographically distant maize populations. Introgressed mexicana ancestry in maize is reduced in lower-recombination rate quintiles of the genome and around domestication genes, consistent with pervasive selection against introgression. However, we also find mexicana ancestry increases across the sampled elevational gradient and that high introgression peaks are most commonly shared among high-elevation maize populations, consistent with introgression from mexicana facilitating adaptation to the highland environment. In the other direction, we find patterns consistent with adaptive and clinal introgression of maize ancestry into sympatric mexicana at many loci across the genome, suggesting that maize also contributes to adaptation in mexicana, especially at the lower end of its elevational range. In sympatric maize, in addition to high introgression regions we find many genomic regions where selection for local adaptation maintains steep gradients in introgressed mexicana ancestry across elevation, including at least two inversions: the well-characterized 14 Mb Inv4m on chromosome 4 and a novel 3 Mb inversion Inv9f surrounding the macrohairless1 locus on chromosome 9. Most outlier loci with high mexicana introgression show no signals of sweeps or local sourcing from sympatric populations and so likely represent ancestral introgression sorted by selection, resulting in correlated but distinct outcomes of introgression in different contemporary maize populations.

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.

The Formation of Bivalents and the Control of Plant Meiotic Recombination

During the first meiotic division, the segregation of homologous chromosomes depends on the physical association of the recombined homologous DNA molecules. The physical tension due to the sites of crossing-overs (COs) is essential for the meiotic spindle to segregate the connected homologous chromosomes to the opposite poles of the cell. This equilibrated partition of homologous chromosomes allows the first meiotic reductional division. Thus, the segregation of homologous chromosomes is dependent on their recombination. In this review, we will detail the recent advances in the knowledge of the mechanisms of recombination and bivalent formation in plants. In plants, the absence of meiotic checkpoints allows observation of subsequent meiotic events in absence of meiotic recombination or defective meiotic chromosomal axis formation such as univalent formation instead of bivalents. Recent discoveries, mainly made in Arabidopsis, rice, and maize, have highlighted the link between the machinery of double-strand break (DSB) formation and elements of the chromosomal axis. We will also discuss the implications of what we know about the mechanisms regulating the number and spacing of COs (obligate CO, CO homeostasis, and interference) in model and crop plants.

Genomic structural variants constrain and facilitate adaptation in natural populations of Theobroma cacao, the chocolate tree

Genomic structural variants (SVs) can play important roles in adaptation and speciation. Yet the overall fitness effects of SVs are poorly understood, partly because accurate population-level identification of SVs requires multiple high-quality genome assemblies. Here, we use 31 chromosome-scale, haplotype-resolved genome assemblies of Theobroma cacao-an outcrossing, long-lived tree species that is the source of chocolate-to investigate the fitness consequences of SVs in natural populations. Among the 31 accessions, we find over 160,000 SVs, which together cover eight times more of the genome than single-nucleotide polymorphisms and short indels (125 versus 15 Mb). Our results indicate that a vast majority of these SVs are deleterious: they segregate at low frequencies and are depleted from functional regions of the genome. We show that SVs influence gene expression, which likely impairs gene function and contributes to the detrimental effects of SVs. We also provide empirical support for a theoretical prediction that SVs, particularly inversions, increase genetic load through the accumulation of deleterious nucleotide variants as a result of suppressed recombination. Despite the overall detrimental effects, we identify individual SVs bearing signatures of local adaptation, several of which are associated with genes differentially expressed between populations. Genes involved in pathogen resistance are strongly enriched among these candidates, highlighting the contribution of SVs to this important local adaptation trait. Beyond revealing empirical evidence for the evolutionary importance of SVs, these 31 de novo assemblies provide a valuable resource for genetic and breeding studies in T cacao.

Crossover-active regions of the wheat genome are distinguished by DMC1, the chromosome axis, H3K27me3, and signatures of adaptation

The hexaploid bread wheat genome comprises over 16 gigabases of sequence across 21 chromosomes. Meiotic crossovers are highly polarized along the chromosomes, with elevation in the gene-dense distal regions and suppression in the Gypsy retrotransposon-dense centromere-proximal regions. We profiled the genomic landscapes of the meiotic recombinase DMC1 and the chromosome axis protein ASY1 in wheat and investigated their relationships with crossovers, chromatin state, and genetic diversity. DMC1 and ASY1 chromatin immunoprecipitation followed by sequencing (ChIP-seq) revealed strong co-enrichment in the distal, crossover-active regions of the wheat chromosomes. Distal ChIP-seq enrichment is consistent with spatiotemporally biased cytological immunolocalization of DMC1 and ASY1 close to the telomeres during meiotic prophase I. DMC1 and ASY1 ChIP-seq peaks show significant overlap with genes and transposable elements in the Mariner and Mutator superfamilies. However, DMC1 and ASY1 ChIP-seq peaks were detected along the length of each chromosome, including in low-crossover regions. At the fine scale, crossover elevation at DMC1 and ASY1 peaks and genes correlates with enrichment of the Polycomb histone modification H3K27me3. This indicates a role for facultative heterochromatin, coincident with high DMC1 and ASY1, in promoting crossovers in wheat and is reflected in distalized H3K27me3 enrichment observed via ChIP-seq and immunocytology. Genes with elevated crossover rates and high DMC1 and ASY1 ChIP-seq signals are overrepresented for defense-response and immunity annotations, have higher sequence polymorphism, and exhibit signatures of selection. Our findings are consistent with meiotic recombination promoting genetic diversity, shaping host–pathogen co-evolution, and accelerating adaptation by increasing the efficiency of selection.

A Modified Meiotic Recombination in Brassica napus Largely Improves Its Breeding Efficiency

Simple Summary

The selection of varieties more resilient to disease and climate change requires generating new genetic diversity for breeding. The main mechanism for reshuffling genetic information is through the recombination of chromosomes during meiosis. We showed in oilseed rape (Brassica napus, AACC, 2n = 4x = 38), which is a natural hybrid formed from a cross between turnip (B. rapa, AA, 2n = 2x = 20) and cabbage (B. oleracea, CC, 2n = 2x = 18), that there is significantly more crossovers occurring along the entire A chromosomes in allotriploid AAC (crossbetween B. napus and B. rapa) than in diploid AA or allotetraploid AACC hybrids. We demonstrated that these allotriploid AAC hybrids are highly efficient to introduce new variability within oilseed rape varieties, notably by enabling the introduction of small genomic regions carrying genes controlling agronomically interesting traits.

Abstract

Meiotic recombination is the main tool used by breeders to generate biodiversity, allowing genetic reshuffling at each generation. It enables the accumulation of favorable alleles while purging deleterious mutations. However, this mechanism is highly regulated with the formation of one to rarely more than three crossovers, which are not randomly distributed. In this study, we showed that it is possible to modify these controls in oilseed rape (Brassica napus, AACC, 2n = 4x = 38) and that it is linked to AAC allotriploidy and not to polyploidy per se. To that purpose, we compared the frequency and the distribution of crossovers along A chromosomes from hybrids carrying exactly the same A nucleotide sequence, but presenting three different ploidy levels: AA, AAC and AACC. Genetic maps established with 202 SNPs anchored on reference genomes revealed that the crossover rate is 3.6-fold higher in the AAC allotriploid hybrids compared to AA and AACC hybrids. Using a higher SNP density, we demonstrated that smaller and numerous introgressions of B. rapa were present in AAC hybrids compared to AACC allotetraploid hybrids, with 7.6 Mb vs. 16.9 Mb on average and 21 B. rapa regions per plant vs. nine regions, respectively. Therefore, this boost of recombination is highly efficient to reduce the size of QTL carried in cold regions of the oilseed rape genome, as exemplified here for a QTL conferring blackleg resistance.

BREC: an R package/Shiny app for automatically identifying heterochromatin boundaries and estimating local recombination rates along chromosomes

Background

Meiotic recombination is a vital biological process playing an essential role in genome's structural and functional dynamics. Genomes exhibit highly various recombination profiles along chromosomes associated with several chromatin states. However, eu-heterochromatin boundaries are not available nor easily provided for non-model organisms, especially for newly sequenced ones. Hence, we miss accurate local recombination rates necessary to address evolutionary questions.

Results

Here, we propose an automated computational tool, based on the Marey maps method, allowing to identify heterochromatin boundaries along chromosomes and estimating local recombination rates. Our method, called BREC (heterochromatin Boundaries and RECombination rate estimates) is non-genome-specific, running even on non-model genomes as long as genetic and physical maps are available. BREC is based on pure statistics and is data-driven, implying that good input data quality remains a strong requirement. Therefore, a data pre-processing module (data quality control and cleaning) is provided. Experiments show that BREC handles different markers' density and distribution issues.

Conclusions

BREC's heterochromatin boundaries have been validated with cytological equivalents experimentally generated on the fruit fly Drosophila melanogaster genome, for which BREC returns congruent corresponding values. Also, BREC's recombination rates have been compared with previously reported estimates. Based on the promising results, we believe our tool has the potential to help bring data science into the service of genome biology and evolution. We introduce BREC within an R-package and a Shiny web-based user-friendly application yielding a fast, easy-to-use, and broadly accessible resource. The BREC R-package is available at the GitHub repository https://github.com/GenomeStructureOrganization.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04233-1.

Evolution of Recombination Landscapes in Diverging Populations of Bread Wheat

Reciprocal exchanges of DNA (crossovers) that occur during meiosis are mandatory to ensure the production of fertile gametes in sexually reproducing species. They also contribute to shuffle parental alleles into new combinations thereby fueling genetic variation and evolution. However, due to biological constraints, the recombination landscape is highly heterogeneous along the genome which limits the range of allelic combinations and the adaptability of populations. An approach to better understand the constraints on the recombination process is to study how it evolved in the past. In this work, we tackled this question by constructing recombination profiles in four diverging bread wheat (Triticum aestivum L.) populations established from 371 landraces genotyped at 200,062 SNPs. We used linkage disequilibrium (LD) patterns to estimate in each population the past distribution of recombination along the genome and characterize its fine-scale heterogeneity. At the megabase scale, recombination rates derived from LD patterns were consistent with family-based estimates obtained from a population of 406 recombinant inbred lines. Among the four populations, recombination landscapes were positively correlated between each other and shared a statistically significant proportion of highly recombinant intervals. However, this comparison also highlighted that the similarity in recombination landscapes between populations was significantly decreasing with their genetic differentiation in most regions of the genome. This observation was found to be robust to SNPs ascertainment and demography and suggests a relatively rapid evolution of factors determining the fine-scale localization of recombination in bread wheat.

Rewiring Meiosis for Crop Improvement

Meiosis is a specialized cell division that contributes to halve the genome content and reshuffle allelic combinations between generations in sexually reproducing eukaryotes. During meiosis, a large number of programmed DNA double-strand breaks (DSBs) are formed throughout the genome. Repair of meiotic DSBs facilitates the pairing of homologs and forms crossovers which are the reciprocal exchange of genetic information between chromosomes. Meiotic recombination also influences centromere organization and is essential for proper chromosome segregation. Accordingly, meiotic recombination drives genome evolution and is a powerful tool for breeders to create new varieties important to food security. Modifying meiotic recombination has the potential to accelerate plant breeding but it can also have detrimental effects on plant performance by breaking beneficial genetic linkages. Therefore, it is essential to gain a better understanding of these processes in order to develop novel strategies to facilitate plant breeding. Recent progress in targeted recombination technologies, chromosome engineering, and an increasing knowledge in the control of meiotic chromosome segregation has significantly increased our ability to manipulate meiosis. In this review, we summarize the latest findings and technologies on meiosis in plants. We also highlight recent attempts and future directions to manipulate crossover events and control the meiotic division process in a breeding perspective.