The genomic landscape of meiotic crossovers and gene conversions in Arabidopsis thaliana

A Pedigree-Based Map of Crossovers and Noncrossovers in Aye-Ayes (Daubentonia madagascariensis)

Gaining a better understanding of the rates and patterns of meiotic recombination is crucial for improving evolutionary genomic modeling, with applications ranging from demographic to selective inference. Although previous research has provided important insights into the landscape of crossovers in humans and other haplorrhines, our understanding of both the considerably more common outcome of recombination (i.e. noncrossovers) as well as the landscapes in more distantly related primates (i.e. strepsirrhines) remains limited owing to difficulties associated with both the identification of noncrossover tracts as well as species sampling. Thus, in order to elucidate recombination patterns in this understudied branch of the primate clade, we here characterize crossover and noncrossover landscapes in aye-ayes utilizing whole-genome sequencing data from six three-generation pedigrees and three two-generation multi-sibling families, and in so doing provide novel insights into this important evolutionary process shaping genomic diversity in one of the world's most critically endangered primate species.

Scalable eQTL mapping using single-nucleus RNA-sequencing of recombined gametes from a small number of individuals

Phenotypic differences between individuals of a species are often caused by differences in gene expression, which are in turn caused by genetic variation. Expression quantitative trait locus (eQTL) analysis is a methodology by which we can identify such causal variants. Scaling eQTL analysis is costly due to the expense of generating mapping populations, and the collection of matched transcriptomic and genomic information. We developed a rapid eQTL analysis approach using single-cell/nucleus RNA sequencing of gametes from a small number of heterozygous individuals. Patterns of inherited polymorphisms are used to infer the recombinant genomes of thousands of individual gametes and identify how different haplotypes correlate with variation in gene expression. Applied to Arabidopsis pollen nuclei, our approach uncovers both cis- and trans-eQTLs, ultimately mapping variation in a master regulator of sperm cell development that affects the expression of hundreds of genes. This establishes snRNA-sequencing as a powerful, cost-effective method for the mapping of meiotic recombination, addressing the scalability challenges of eQTL analysis and enabling eQTL mapping in specific cell-types.

Mixed Outcomes in Recombination Rates After Domestication: Revisiting Theory and Data

The process of domestication has altered many phenotypes. Selection on these phenotypes has long been hypothesised to indirectly select for increases in the genome-wide recombination rate. This hypothesis is potentially consistent with theory on the evolution of the recombination rate, but empirical support has been unclear. We review relevant theory, lab-based experiments, and data comparing recombination rates in wild progenitors and their domesticated counterparts. We utilise population sequencing data and a deep learning method to infer genome-wide recombination rates for new comparisons of chicken/red junglefowl, sheep/mouflon, and goat/bezoar. We find evidence of increased recombination in domesticated goats compared to bezoars but more mixed results in chicken and generally decreased recombination in domesticated sheep compared to mouflon. Our results add to a growing body of literature in plants and animals that finds no consistent evidence of an increase in genome-wide recombination with domestication.

The recombination landscape of introgression in yeast

Meiotic recombination is an evolutionary force that acts by breaking up genomic linkage, increasing the efficacy of selection. Recombination is initiated with a double-strand break which is resolved via a crossover, which involves the reciprocal exchange of genetic material between homologous chromosomes, or a non-crossover, which results in small tracts of non-reciprocal exchange of genetic material. Crossover and non-crossover rates vary between species, populations, individuals, and across the genome. In recent years, recombination rate has been associated with the distribution of ancestry derived from past interspecific hybridization (introgression) in a variety of species. We explore this interaction of recombination and introgression by sequencing spores and detecting crossovers and non-crossovers from two crosses of the yeast Saccharomyces uvarum. One cross is between strains which each contain introgression from their sister species, S. eubayanus, while the other cross has no introgression present. We find that the recombination landscape is significantly different between S. uvarum crosses, and that some of these differences can be explained by the presence of introgression in one cross. Crossovers are significantly reduced in heterozygous introgression compared to syntenic regions in the cross without introgression. This translates to reduced allele shuffling within introgressed regions, and an overall reduction of shuffling on most chromosomes with introgression compared to the syntenic regions and chromosomes without introgression. Our results suggest that hybridization can significantly influence the recombination landscape, and that the reduction in allele shuffling contributes to the initial purging of introgression in the generations following a hybridization event.

The identification and analysis of meristematic mutations within the apple tree that developed the RubyMac sport mutation

BACKGROUND

Understanding the molecular basis of sport mutations in fruit trees has the potential to accelerate generation of improved cultivars.

RESULTS

For this, we analyzed the genome of the apple tree that developed the RubyMac phenotype through a sport mutation that led to the characteristic fruit coloring of this variety. Overall, we found 46 somatic mutations that distinguished the mutant and wild-type branches of the tree. In addition, we found 54 somatic gene conversions (i.e., loss-of-heterozygosity mutations) that also distinguished the two parts of the tree. Approximately 20% of the mutations were specific to individual cell lineages, suggesting that they originated from the corresponding meristematic layers. Interestingly, the de novo mutations were enriched for GC =  > AT transitions while the gene conversions showed the opposite bias for AT =  > GC transitions, suggesting that GC-biased gene conversions have the potential to counteract the AT-bias of de novo mutations. By comparing the gene expression patterns in fruit skins from mutant and wild-type branches, we found 56 differentially expressed genes including 18 involved in anthocyanin biosynthesis. While none of the differently expressed genes harbored a somatic mutation, we found that some of them in regions of the genome that were recently associated with natural variation in fruit coloration.

CONCLUSION

Our analysis revealed insights in the characteristics of somatic change, which not only included de novo mutations but also gene conversions. Some of these somatic changes displayed strong candidate mutations for the change in fruit coloration in RubyMac.

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.

Diversity in Recombination Hotspot Characteristics and Gene Structure Shape Fine-Scale Recombination Patterns in Plant Genomes

During the meiosis of many eukaryote species, crossovers tend to occur within narrow regions called recombination hotspots. In plants, it is generally thought that gene regulatory sequences, especially promoters and 5' to 3' untranslated regions, are enriched in hotspots, but this has been characterized in a handful of species only. We also lack a clear description of fine-scale variation in recombination rates within genic regions and little is known about hotspot position and intensity in plants. To address this question, we constructed fine-scale recombination maps from genetic polymorphism data and inferred recombination hotspots in 11 plant species. We detected gradients of recombination in genic regions in most species, yet gradients varied in intensity and shape depending on specific hotspot locations and gene structure. To further characterize recombination gradients, we decomposed them according to gene structure by rank and number of exons. We generalized the previously observed pattern that recombination hotspots are organized around the boundaries of coding sequences, especially 5' promoters. However, our results also provided new insight into the relative importance of the 3' end of genes in some species and the possible location of hotspots away from genic regions in some species. Variation among species seemed driven more by hotspot location among and within genes than by differences in size or intensity among species. Our results shed light on the variation in recombination rates at a very fine scale, revealing the diversity and complexity of genic recombination gradients emerging from the interaction between hotspot location and gene structure.

Male-biased recombination at chromosome ends in a songbird revealed by precisely mapping crossover positions

Recombination plays a crucial role in evolution by generating novel haplotypes and disrupting linkage between genes, thereby enhancing the efficiency of selection. Here, we analyze the genomes of 12 great reed warblers (Acrocephalus arundinaceus) in a 3-generation pedigree to identify precise crossover positions along the chromosomes. We located more than 200 crossovers and found that these were highly concentrated toward the telomeric ends of the chromosomes. Apart from this major pattern in the recombination landscape, we found significantly higher frequencies of crossovers in genic compared with intergenic regions, and in exons compared with introns. Moreover, while the number of recombination events was similar between the sexes, the crossovers were located significantly closer to the ends of paternal compared with maternal chromosomes. In conclusion, our study of the great reed warbler revealed substantial variation in crossover frequencies within chromosomes, with a distinct bias toward the sub-telomeric regions, particularly on the paternal side. These findings emphasize the importance of thoroughly screening the entire length of chromosomes to characterize the recombination landscape and uncover potential sex-biases in recombination.

The recombination landscape of introgression in yeast

Meiotic recombination is an evolutionary force that acts by breaking up genomic linkage, increasing the efficacy of selection. Recombination is initiated with a double-strand break which is resolved via a crossover, which involves the reciprocal exchange of genetic material between homologous chromosomes, or a non-crossover, which results in small tracts of non-reciprocal exchange of genetic material. Crossover and non-crossover rates vary between species, populations, individuals, and across the genome. In recent years, recombination rate has been associated with the distribution of ancestry derived from past interspecific hybridization (introgression) in a variety of species. We explore this interaction of recombination and introgression by sequencing spores and detecting crossovers and non-crossovers from two crosses of the yeast Saccharomyces uvarum. One cross is between strains which each contain introgression from their sister species, S. eubayanus, while the other cross has no introgression present. We find that the recombination landscape is significantly different between S. uvarum crosses, and that some of these differences can be explained by the presence of introgression in one cross. Crossovers are reduced and non-crossovers are increased in heterozygous introgression compared to syntenic regions in the cross without introgression. This translates to reduced allele shuffling within introgressed regions, and an overall reduction of shuffling on most chromosomes with introgression compared to the syntenic regions and chromosomes without introgression. Our results suggest that hybridization can significantly influence the recombination landscape, and that the reduction in allele shuffling contributes to the initial purging of introgression in the generations following a hybridization event.

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.

MCPtaggR: R package for accurate genotype calling in reduced representation sequencing data by eliminating error-prone markers based on genome comparison

Reduced representation sequencing (RRS) offers cost-effective, high-throughput genotyping platforms such as genotyping-by-sequencing (GBS). RRS reads are typically mapped onto a reference genome. However, mapping reads harbouring mismatches against the reference can potentially result in mismapping and biased mapping, leading to the detection of error-prone markers that provide incorrect genotype information. We established a genotype-calling pipeline named mappable collinear polymorphic tag genotyping (MCPtagg) to achieve accurate genotyping by eliminating error-prone markers. MCPtagg was designed for the RRS-based genotyping of a population derived from a biparental cross. The MCPtagg pipeline filters out error-prone markers prior to genotype calling based on marker collinearity information obtained by comparing the genome sequences of the parents of a population to be genotyped. A performance evaluation on real GBS data from a rice F2 population confirmed its effectiveness. Furthermore, our performance test using a genome assembly that was obtained by genome sequence polishing on an available genome assembly suggests that our pipeline performs well with converted genomes, rather than necessitating de novo assembly. This demonstrates its flexibility and scalability. The R package, MCPtaggR, was developed to provide functions for the pipeline and is available at https://github.com/tomoyukif/MCPtaggR.

Allelic gene conversion softens selective sweeps

The prominence of positive selection, in which beneficial mutations are favored by natural selection and rapidly increase in frequency, is a subject of intense debate. Positive selection can result in selective sweeps, in which the haplotype(s) bearing the adaptive allele "sweep" through the population, thereby removing much of the genetic diversity from the region surrounding the target of selection. Two models of selective sweeps have been proposed: classical sweeps, or "hard sweeps", in which a single copy of the adaptive allele sweeps to fixation, and "soft sweeps", in which multiple distinct copies of the adaptive allele leave descendants after the sweep. Soft sweeps can be the outcome of recurrent mutation to the adaptive allele, or the presence of standing genetic variation consisting of multiple copies of the adaptive allele prior to the onset of selection. Importantly, soft sweeps will be common when populations can rapidly adapt to novel selective pressures, either because of a high mutation rate or because adaptive alleles are already present. The prevalence of soft sweeps is especially controversial, and it has been noted that selection on standing variation or recurrent mutations may not always produce soft sweeps. Here, we show that the inverse is true: selection on single-origin de novo mutations may often result in an outcome that is indistinguishable from a soft sweep. This is made possible by allelic gene conversion, which "softens" hard sweeps by copying the adaptive allele onto multiple genetic backgrounds, a process we refer to as a "pseudo-soft" sweep. We carried out a simulation study examining the impact of gene conversion on sweeps from a single de novo variant in models of human, Drosophila, and Arabidopsis populations. The fraction of simulations in which gene conversion had produced multiple haplotypes with the adaptive allele upon fixation was appreciable. Indeed, under realistic demographic histories and gene conversion rates, even if selection always acts on a single-origin mutation, sweeps involving multiple haplotypes are more likely than hard sweeps in large populations, especially when selection is not extremely strong. Thus, even when the mutation rate is low or there is no standing variation, hard sweeps are expected to be the exception rather than the rule in large populations. These results also imply that the presence of signatures of soft sweeps does not necessarily mean that adaptation has been especially rapid or is not mutation limited.

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.

The human fungal pathogen Aspergillus fumigatus can produce the highest known number of meiotic crossovers

Sexual reproduction involving meiosis is essential in most eukaryotes. This produces offspring with novel genotypes, both by segregation of parental chromosomes as well as crossovers between homologous chromosomes. A sexual cycle for the opportunistic human pathogenic fungus Aspergillus fumigatus is known, but the genetic consequences of meiosis have remained unknown. Among other Aspergilli, it is known that A. flavus has a moderately high recombination rate with an average of 4.2 crossovers per chromosome pair, whereas A. nidulans has in contrast a higher rate with 9.3 crossovers per chromosome pair. Here, we show in a cross between A. fumigatus strains that they produce an average of 29.9 crossovers per chromosome pair and large variation in total map length across additional strain crosses. This rate of crossovers per chromosome is more than twice that seen for any known organism, which we discuss in relation to other genetic model systems. We validate this high rate of crossovers through mapping of resistance to the laboratory antifungal acriflavine by using standing variation in an undescribed ABC efflux transporter. We then demonstrate that this rate of crossovers is sufficient to produce one of the common multidrug resistant haplotypes found in the cyp51A gene (TR34/L98H) in crosses among parents harboring either of 2 nearby genetic variants, possibly explaining the early spread of such haplotypes. Our results suggest that genomic studies in this species should reassess common assumptions about linkage between genetic regions. The finding of an unparalleled crossover rate in A. fumigatus provides opportunities to understand why these rates are not generally higher in other eukaryotes.

Expanding the stdpopsim species catalog, and lessons learned for realistic genome simulations

Simulation is a key tool in population genetics for both methods development and empirical research, but producing simulations that recapitulate the main features of genomic datasets remains a major obstacle. Today, more realistic simulations are possible thanks to large increases in the quantity and quality of available genetic data, and the sophistication of inference and simulation software. However, implementing these simulations still requires substantial time and specialized knowledge. These challenges are especially pronounced for simulating genomes for species that are not well-studied, since it is not always clear what information is required to produce simulations with a level of realism sufficient to confidently answer a given question. The community-developed framework stdpopsim seeks to lower this barrier by facilitating the simulation of complex population genetic models using up-to-date information. The initial version of stdpopsim focused on establishing this framework using six well-characterized model species (Adrion et al., 2020). Here, we report on major improvements made in the new release of stdpopsim (version 0.2), which includes a significant expansion of the species catalog and substantial additions to simulation capabilities. Features added to improve the realism of the simulated genomes include non-crossover recombination and provision of species-specific genomic annotations. Through community-driven efforts, we expanded the number of species in the catalog more than threefold and broadened coverage across the tree of life. During the process of expanding the catalog, we have identified common sticking points and developed the best practices for setting up genome-scale simulations. We describe the input data required for generating a realistic simulation, suggest good practices for obtaining the relevant information from the literature, and discuss common pitfalls and major considerations. These improvements to stdpopsim aim to further promote the use of realistic whole-genome population genetic simulations, especially in non-model organisms, making them available, transparent, and accessible to everyone.

Patterns, mechanisms, and consequences of homoeologous exchange in allopolyploid angiosperms: a genomic and epigenomic perspective

Allopolyploids result from hybridization between different evolutionary lineages coupled with genome doubling. Homoeologous chromosomes (chromosomes with common shared ancestry) may undergo recombination immediately after allopolyploid formation and continue over successive generations. The outcome of this meiotic pairing behavior is dynamic and complex. Homoeologous exchanges (HEs) may lead to the formation of unbalanced gametes, reduced fertility, and selective disadvantage. By contrast, HEs could act as sources of novel evolutionary substrates, shifting the relative dosage of parental gene copies, generating novel phenotypic diversity, and helping the establishment of neo-allopolyploids. However, HE patterns vary among lineages, across generations, and even within individual genomes and chromosomes. The causes and consequences of this variation are not fully understood, though interest in this evolutionary phenomenon has increased in the last decade. Recent technological advances show promise in uncovering the mechanistic basis of HEs. Here, we describe recent observations of the common patterns among allopolyploid angiosperm lineages, underlying genomic and epigenomic features, and consequences of HEs. We identify critical research gaps and discuss future directions with far-reaching implications in understanding allopolyploid evolution and applying them to the development of important phenotypic traits of polyploid crops.

Reliable genotyping of recombinant genomes using a robust hidden Markov model

Meiotic recombination is an essential mechanism during sexual reproduction and includes the exchange of chromosome segments between homologous chromosomes. New allelic combinations are transmitted to the new generation, introducing novel genetic variation in the offspring genomes. With the improvement of high-throughput whole-genome sequencing technologies, large numbers of recombinant individuals can now be sequenced with low sequencing depth at low costs, necessitating computational methods for reconstructing their haplotypes. The main challenge is the uncertainty in haplotype calling that arises from the low information content of a single genomic position. Straightforward sliding window-based approaches are difficult to tune and fail to place recombination breakpoints precisely. Hidden Markov model (HMM)-based approaches, on the other hand, tend to over-segment the genome. Here, we present RTIGER, an HMM-based model that exploits in a mathematically precise way the fact that true chromosome segments typically have a certain minimum length. We further separate the task of identifying the correct haplotype sequence from the accurate placement of haplotype borders, thereby maximizing the accuracy of border positions. By comparing segmentations based on simulated data with known underlying haplotypes, we highlight the reasons for RTIGER outperforming traditional segmentation approaches. We then analyze the meiotic recombination pattern of segregants of 2 Arabidopsis (Arabidopsis thaliana) accessions and a previously described hyper-recombining mutant. RTIGER is available as an R package with an efficient Julia implementation of the core algorithm.

Cycles of satellite and transposon evolution in Arabidopsis centromeres

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation^1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox^3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

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.

Control of meiotic crossing over in plant breeding

Meiotic crossing over is the main mechanism for constructing a new allelic composition of individual chromosomes and is necessary for the proper distribution of homologous chromosomes between gametes. The parameters of meiotic crossing over that have developed in the course of evolution are determined by natural selection and do not fully suit the tasks of selective breeding research. This review summarizes the results of experimental studies aimed at increasing the frequency of crossovers and redistributing their positions along chromosomes using genetic manipulations at different stages of meiotic recombination. The consequences of inactivation and/or overexpression of the SPO11 genes, the products of which generate meiotic double-strand breaks in DNA, for the redistribution of crossover positions in the genome of various organisms are discussed. The results of studies concerning the effect of inactivation or overexpression of genes encoding RecA-like recombinases on meiotic crossing over, including those in cultivated tomato (Solanum lycopersicum L.) and its interspecific hybrids, are summarized. The consequences of inactivation of key genes of the mismatch repair system are discussed. Their suppression made it possible to significantly increase the frequency of meiotic recombination between homeologues in the interspecific hybrid yeast Saccharomyces cerevisiae × S. paradoxus and between homologues in arabidopsis plants (Arabidopsis thaliana L.). Also discussed are attempts to extrapolate these results to other plant species, in which a decrease in reproductive properties and microsatellite instability in the genome have been noted. The most significant results on the meiotic recombination frequency increase upon inactivation of the FANCM, TOP3α, RECQ4, FIGL1 crossover repressor genes and upon overexpression of the HEI10 crossover enhancer gene are separately described. In some experiments, the increase of meiotic recombination frequency by almost an order of magnitude and partial redistribution of the crossover positions along chromosomes were achieved in arabidopsis while fully preserving fecundity. Similar results have been obtained for some crops.

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.

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.

Dissecting the Meiotic Recombination Patterns in a Brassica napus Double Haploid Population Using 60K SNP Array

Meiotic recombination not only maintains the stability of the chromosome structure but also creates genetic variations for adapting to changeable environments. A better understanding of the mechanism of crossover (CO) patterns at the population level is useful for crop improvement. However, there are limited cost-effective and universal methods to detect the recombination frequency at the population level in Brassica napus. Here, the Brassica 60K Illumina Infinium SNP array (Brassica 60K array) was used to systematically study the recombination landscape in a double haploid (DH) population of B. napus. It was found that COs were unevenly distributed across the whole genome, and a higher frequency of COs existed at the distal ends of each chromosome. A considerable number of genes (more than 30%) in the CO hot regions were associated with plant defense and regulation. In most tissues, the average gene expression level in the hot regions (CO frequency of greater than 2 cM/Mb) was significantly higher than that in the regions with a CO frequency of less than 1 cM/Mb. In addition, a bin map was constructed with 1995 recombination bins. For seed oil content, Bin 1131 to 1134, Bin 1308 to 1311, Bin 1864 to 1869, and Bin 2184 to 2230 were identified on chromosomes A08, A09, C03, and C06, respectively, which could explain 8.5%, 17.3%, 8.6%, and 3.9% of the phenotypic variation. These results could not only deepen our understanding of meiotic recombination in B. napus at the population level, and provide useful information for rapeseed breeding in the future, but also provided a reference for studying CO frequency in other species.

Genomic targets of HOP2 are enriched for features found at recombination hotspots

HOP2 is a conserved protein that plays a positive role in homologous chromosome pairing and a separable role in preventing illegitimate connections between nonhomologous chromosome regions during meiosis. We employed ChIP-seq to discover that Arabidopsis HOP2 binds along the length of all chromosomes, except for centromeric and nucleolar organizer regions, and no binding sites were detected in the organelle genomes. A large number of reads were assigned to the HOP2 locus itself, yet TAIL-PCR and SNP analysis of the aligned sequences indicate that many of these reads originate from the transforming T-DNA, supporting the role of HOP2 in preventing nonhomologous exchanges. The 292 ChIP-seq peaks are largely found in promoter regions and downstream from genes, paralleling the distribution of recombination hotspots, and motif analysis revealed that there are several conserved sequences that are also enriched at crossover sites. We conducted coimmunoprecipitation of HOP2 followed by LC-MS/MS and found enrichment for several proteins, including some histone variants and modifications that are also known to be associated with recombination hotspots. We propose that HOP2 may be directed to chromatin motifs near double strand breaks, where homology checks are proposed to occur.

Lessons from the meiotic recombination landscape of the ZMM deficient budding yeast Lachancea waltii

Meiotic recombination is a driving force for genome evolution, deeply characterized in a few model species, notably in the budding yeast Saccharomyces cerevisiae. Interestingly, Zip2, Zip3, Zip4, Spo16, Msh4, and Msh5, members of the so-called ZMM pathway that implements the interfering meiotic crossover pathway in S. cerevisiae, have been lost in Lachancea yeast species after the divergence of Lachancea kluyveri from the rest of the clade. In this context, after investigating meiosis in L. kluyveri, we determined the meiotic recombination landscape of Lachancea waltii. Attempts to generate diploid strains with fully hybrid genomes invariably resulted in strains with frequent whole-chromosome aneuploidy and multiple extended regions of loss of heterozygosity (LOH), which mechanistic origin is so far unclear. Despite the lack of multiple ZMM pro-crossover factors in L. waltii, numbers of crossovers and noncrossovers per meiosis were higher than in L. kluyveri but lower than in S. cerevisiae, for comparable genome sizes. Similar to L. kluyveri but opposite to S. cerevisiae, L. waltii exhibits an elevated frequency of zero-crossover bivalents. Lengths of gene conversion tracts for both crossovers and non-crossovers in L. waltii were comparable to those observed in S. cerevisiae and shorter than in L. kluyveri despite the lack of Mlh2, a factor limiting conversion tract size in S. cerevisiae. L. waltii recombination hotspots were not shared with either S. cerevisiae or L. kluyveri, showing that meiotic recombination hotspots can evolve at a rather limited evolutionary scale within budding yeasts. Finally, L. waltii crossover interference was reduced relative to S. cerevisiae, with interference being detected only in the 25 kb distance range. Detection of positive inference only at short distance scales in the absence of multiple ZMM factors required for interference-sensitive crossovers in other systems likely reflects interference between early recombination precursors such as DSBs.

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.

The Linkage-Disequilibrium and Recombinational Landscape in Daphnia pulex

By revealing the influence of recombinational activity beyond what can be achieved with controlled crosses, measures of linkage disequilibrium (LD) in natural populations provide a powerful means of defining the recombinational landscape within which genes evolve. In one of the most comprehensive studies of this sort ever performed, involving whole-genome analyses on nearly 1,000 individuals of the cyclically parthenogenetic microcrustacean Daphnia pulex, the data suggest a relatively uniform pattern of recombination across the genome. Patterns of LD are quite consistent among populations; average rates of recombination are quite similar for all chromosomes; and although some chromosomal regions have elevated recombination rates, the degree of inflation is not large, and the overall spatial pattern of recombination is close to the random expectation. Contrary to expectations for models in which crossing-over is the primary mechanism of recombination, and consistent with data for other species, the distance-dependent pattern of LD indicates excessively high levels at both short and long distances and unexpectedly low levels of decay at long distances, suggesting significant roles for factors such as nonindependent mutation, population subdivision, and recombination mechanisms unassociated with crossing over. These observations raise issues regarding the classical LD equilibrium model widely applied in population genetics to infer recombination rates across various length scales on chromosomes.

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.

Epigenetic Dynamics and Regulation of Plant Male Reproduction

Flowering plant male germlines develop within anthers and undergo epigenetic reprogramming with dynamic changes in DNA methylation, chromatin modifications, and small RNAs. Profiling the epigenetic status using different technologies has substantially accumulated information on specific types of cells at different stages of male reproduction. Many epigenetically related genes involved in plant gametophyte development have been identified, and the mutation of these genes often leads to male sterility. Here, we review the recent progress on dynamic epigenetic changes during pollen mother cell differentiation, microsporogenesis, microgametogenesis, and tapetal cell development. The reported epigenetic variations between male fertile and sterile lines are summarized. We also summarize the epigenetic regulation-associated male sterility genes and discuss how epigenetic mechanisms in plant male reproduction can be further revealed.

Natural and artificial sources of genetic variation used in crop breeding: A baseline comparator for genome editing

Traditional breeding has successfully selected beneficial traits for food, feed, and fibre crops over the last several thousand years. The last century has seen significant technological advancements particularly in marker assisted selection and the generation of induced genetic variation, including over the last few decades, through mutation breeding, genetic modification, and genome editing. While regulatory frameworks for traditional varietal development and for genetic modification with transgenes are broadly established, those for genome editing are lacking or are still evolving in many regions. In particular, the lack of "foreign" recombinant DNA in genome edited plants and that the resulting SNPs or INDELs are indistinguishable from those seen in traditional breeding has challenged development of new legislation. Where products of genome editing and other novel breeding technologies possess no transgenes and could have been generated via traditional methods, we argue that it is logical and proportionate to apply equivalent legislative oversight that already exists for traditional breeding and novel foods. This review analyses the types and the scale of spontaneous and induced genetic variation that can be selected during traditional plant breeding activities. It provides a base line from which to judge whether genetic changes brought about by techniques of genome editing or other reverse genetic methods are indeed comparable to those routinely found using traditional methods of plant breeding.

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.

Knock-down of gene expression throughout meiosis and pollen formation by virus-induced gene silencing in Arabidopsis thaliana

Through the inactivation of genes that act during meiosis it is possible to direct the genetic make-up of plants in subsequent generations and optimize breeding schemes. Offspring may show higher recombination of parental alleles resulting from elevated crossover (CO) incidence, or by omission of meiotic divisions, offspring may become polyploid. However, stable mutations in genes essential for recombination, or for either one of the two meiotic divisions, can have pleiotropic effects on plant morphology and line stability, for instance by causing lower fertility. Therefore, it is often favorable to temporarily change gene expression during meiosis rather than relying on stable null mutants. It was previously shown that virus-induced gene silencing (VIGS) can be used to transiently reduce CO frequencies. We asked if VIGS could also be used to modify other processes throughout meiosis and during pollen formation in Arabidopsis thaliana. Here, we show that VIGS-mediated knock-down of FIGL1, RECQ4A/B, OSD1 and QRT2 can induce (i) an increase in chiasma numbers, (ii) unreduced gametes and (iii) pollen tetrads. We further show that VIGS can target both sexes and different genetic backgrounds and can simultaneously silence different gene copies. The successful knock-down of these genes in A. thaliana suggests that VIGS can be exploited to manipulate any process during or shortly after meiosis. Hence, the transient induction of changes in inheritance patterns can be used as a powerful tool for applied research and biotechnological applications.

inGAP-family: Accurate Detection of Meiotic Recombination Loci and Causal Mutations by Filtering Out Artificial Variants due to Genome Complexities

Accurately identifying DNA polymorphisms can bridge the gap between phenotypes and genotypes and is essential for molecular marker assisted genetic studies. Genome complexities, including large-scale structural variations, bring great challenges to bioinformatic analysis for obtaining high-confidence genomic variants, as sequence differences between non-allelic loci of two or more genomes can be misinterpreted as polymorphisms. It is important to correctly filter out artificial variants to avoid false genotyping or estimation of allele frequencies. Here, we present an efficient and effective framework, inGAP-family, to discover, filter, and visualize DNA polymorphisms and structural variants (SVs) from alignment of short reads. Applying this method to polymorphism detection on real datasets shows that elimination of artificial variants greatly facilitates the precise identification of meiotic recombination points as well as causal mutations in mutant genomes or quantitative trait loci. In addition, inGAP-family provides a user-friendly graphical interface for detecting polymorphisms and SVs, further evaluating predicted variants and identifying mutations related to genotypes. It is accessible at https://sourceforge.net/projects/ingap-family/.

RecombineX: A generalized computational framework for automatic high-throughput gamete genotyping and tetrad-based recombination analysis

Meiotic recombination is an essential biological process that ensures faithful chromosome segregation and promotes parental allele shuffling. Tetrad analysis is a powerful approach to quantify the genetic makeups and recombination landscapes of meiotic products. Here we present RecombineX (https://github.com/yjx1217/RecombineX), a generalized computational framework that automates the full workflow of marker identification, gamete genotyping, and tetrad-based recombination profiling based on any organism or genetic background with batch processing capability. Aside from conventional reference-based analysis, RecombineX can also perform analysis based on parental genome assemblies, which facilitates analyzing meiotic recombination landscapes in their native genomic contexts. Additional features such as copy number variation profiling and missing genotype inference further enhance downstream analysis. RecombineX also includes a dedicate module for simulating the genomes and reads of recombinant tetrads, which enables fine-tuned simulation-based hypothesis testing. This simulation module revealed the power and accuracy of RecombineX even when analyzing tetrads with very low sequencing depths (e.g., 1-2X). Tetrad sequencing data from the budding yeast Saccharomyces cerevisiae and green alga Chlamydomonas reinhardtii were further used to demonstrate the accuracy and robustness of RecombineX for organisms with both small and large genomes, manifesting RecombineX as an all-around one stop solution for future tetrad analysis. Interestingly, our re-analysis of the budding yeast tetrad sequencing data with RecombineX and Oxford Nanopore sequencing revealed two unusual structural rearrangement events that were not noticed before, which exemplify the occasional genome instability triggered by meiosis.

High-Resolution Estimates of Crossover and Noncrossover Recombination from a Captive Baboon Colony

Homologous recombination has been extensively studied in humans and a handful of model organisms. Much less is known about recombination in other species, including nonhuman primates. Here, we present a study of crossovers (COs) and noncrossover (NCO) recombination in olive baboons (Papio anubis) from two pedigrees containing a total of 20 paternal and 17 maternal meioses, and compare these results to linkage disequilibrium (LD) based recombination estimates from 36 unrelated olive baboons. We demonstrate how COs, combined with LD-based recombination estimates, can be used to identify genome assembly errors. We also quantify sex-specific differences in recombination rates, including elevated male CO and reduced female CO rates near telomeres. Finally, we add to the increasing body of evidence suggesting that while most NCO recombination tracts in mammals are short (e.g., <500 bp), there is a non-negligible fraction of longer (e.g., >1 kb) NCO tracts. For NCO tracts shorter than 10 kb, we fit a mixture of two (truncated) geometric distributions model to the NCO tract length distribution and estimate that >99% of all NCO tracts are very short (mean 24 bp), but the remaining tracts can be quite long (mean 4.3 kb). A single geometric distribution model for NCO tract lengths is incompatible with the data, suggesting that LD-based methods for estimating NCO recombination rates that make this assumption may need to be modified.

Gene conversion: a non-Mendelian process integral to meiotic recombination

Meiosis is undoubtedly the mechanism that underpins Mendelian genetics. Meiosis is a specialised, reductional cell division which generates haploid gametes (reproductive cells) carrying a single chromosome complement from diploid progenitor cells harbouring two chromosome sets. Through this process, the hereditary material is shuffled and distributed into haploid gametes such that upon fertilisation, when two haploid gametes fuse, diploidy is restored in the zygote. During meiosis the transient physical connection of two homologous chromosomes (one originally inherited from each parent) each consisting of two sister chromatids and their subsequent segregation into four meiotic products (gametes), is what enables genetic marker assortment forming the core of Mendelian laws. The initiating events of meiotic recombination are DNA double-strand breaks (DSBs) which need to be repaired in a certain way to enable the homologous chromosomes to find each other. This is achieved by DSB ends searching for homologous repair templates and invading them. Ultimately, the repair of meiotic DSBs by homologous recombination physically connects homologous chromosomes through crossovers. These physical connections provided by crossovers enable faithful chromosome segregation. That being said, the DSB repair mechanism integral to meiotic recombination also produces genetic transmission distortions which manifest as postmeiotic segregation events and gene conversions. These processes are non-reciprocal genetic exchanges and thus non-Mendelian.

Manipulation of Meiotic Recombination to Hasten Crop Improvement

Reciprocal (cross-overs = COs) and non-reciprocal (gene conversion) DNA exchanges between the parental chromosomes (the homologs) during meiotic recombination are, together with mutation, the drivers for the evolution and adaptation of species. In plant breeding, recombination combines alleles from genetically diverse accessions to generate new haplotypes on which selection can act. In recent years, a spectacular progress has been accomplished in the understanding of the mechanisms underlying meiotic recombination in both model and crop plants as well as in the modulation of meiotic recombination using different strategies. The latter includes the stimulation and redistribution of COs by either modifying environmental conditions (e.g., T°), harnessing particular genomic situations (e.g., triploidy in Brassicaceae), or inactivating/over-expressing meiotic genes, notably some involved in the DNA double-strand break (DSB) repair pathways. These tools could be particularly useful for shuffling diversity in pre-breeding generations. Furthermore, thanks to the site-specific properties of genome editing technologies the targeting of meiotic recombination at specific chromosomal regions nowadays appears an attainable goal. Directing COs at desired chromosomal positions would allow breaking linkage situations existing between favorable and unfavorable alleles, the so-called linkage drag, and accelerate genetic gain. This review surveys the recent achievements in the manipulation of meiotic recombination in plants that could be integrated into breeding schemes to meet the challenges of deploying crops that are more resilient to climate instability, resistant to pathogens and pests, and sparing in their input requirements.

GC content of plant genes is linked to past gene duplications

The frequency of G and C nucleotides in genomes varies from species to species, and sometimes even between different genes in the same genome. The monocot grasses have a bimodal distribution of genic GC content absent in dicots. We categorized plant genes from 5 dicots and 4 monocot grasses by synteny to related species and determined that syntenic genes have significantly higher GC content than non-syntenic genes at their 5`-end in the third position within codons for all 9 species. Lower GC content is correlated with gene duplication, as lack of synteny to distantly related genomes is associated with past interspersed gene duplications. Two mutation types can account for biased GC content, mutation of methylated C to T and gene conversion from A to G. Gene conversion involves non-reciprocal exchanges between homologous alleles and is not detectable when the alleles are identical or heterozygous for presence-absence variation, both likely situations for genes duplicated to new loci. Gene duplication can cause production of siRNA which can induce targeted methylation, elevating mC→T mutations. Recently duplicated plant genes are more frequently methylated and less likely to undergo gene conversion, each of these factors synergistically creating a mutational environment favoring AT nucleotides. The syntenic genes with high GC content in the grasses compose a subset that have undergone few duplications, or for which duplicate copies were purged by selection. We propose a “biased gene duplication / biased mutation” (BDBM) model that may explain the origin and trajectory of the observed link between duplication and genic GC bias. The BDBM model is supported by empirical data based on joint analyses of 9 angiosperm species with their genes categorized by duplication status, GC content, methylation levels and functional classes.

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.

Chromosome-Wide Characterization of Intragenic Crossover in Shiitake Mushroom, Lentinula edodes

Meiotic crossover plays a critical role in generating genetic variations and is a central component of breeding. However, our understanding of crossover in mushroom-forming fungi is limited. Here, in Lentinula edodes, we characterized the chromosome-wide intragenic crossovers, by utilizing the single-nucleotide polymorphisms (SNPs) datasets of an F₁ haploid progeny. A total of 884 intragenic crossovers were identified in 110 single-spore isolates, the majority of which were closer to transcript start sites. About 71.5% of the intragenic crossovers were clustered into 65 crossover hotspots. A 10 bp motif (GCTCTCGAAA) was significantly enriched in the hotspot regions. Crossover frequencies around mating-type A (MAT-A) loci were enhanced and formed a hotspot in L. edodes. Genome-wide quantitative trait loci (QTLs) mapping identified sixteen crossover-QTLs, contributing 8.5–29.1% of variations. Most of the detected crossover-QTLs were co-located with crossover hotspots. Both cis- and trans-QTLs contributed to the nonuniformity of crossover along chromosomes. On chr2, we identified a QTL hotspot that regulated local, global crossover variation and crossover hotspot in L. edodes. These findings and observations provide a comprehensive view of the crossover landscape in L. edodes, and advance our understandings of conservation and diversity of meiotic recombination in mushroom-forming fungi.

The genetic and epigenetic landscape of the Arabidopsis centromeres

Centromeres attach chromosomes to spindle microtubules during cell division and, despite this conserved role, show paradoxically rapid evolution and are typified by complex repeats. We used long-read sequencing to generate the Col-CEN Arabidopsis thaliana genome assembly that resolves all five centromeres. The centromeres consist of megabase-scale tandemly repeated satellite arrays, which support CENTROMERE SPECIFIC HISTONE H3 (CENH3) occupancy and are densely DNA methylated, with satellite variants private to each chromosome. CENH3 preferentially occupies satellites that show the least amount of divergence and occur in higher-order repeats. The centromeres are invaded by ATHILA retrotransposons, which disrupt genetic and epigenetic organization. Centromeric crossover recombination is suppressed, yet low levels of meiotic DNA double-strand breaks occur that are regulated by DNA methylation. We propose that Arabidopsis centromeres are evolving through cycles of satellite homogenization and retrotransposon-driven diversification.

Evolution of crossover interference enables stable autopolyploidy by ensuring pairwise partner connections in Arabidopsis arenosa

Polyploidy is a major driver of evolutionary change. Autopolyploids, which arise by within-species whole-genome duplication, carry multiple nearly identical copies of each chromosome. This presents an existential challenge to sexual reproduction. Meiotic chromosome segregation requires formation of DNA crossovers (COs) between two homologous chromosomes. How can this outcome be achieved when more than two essentially equivalent partners are available? We addressed this question by comparing diploid, neo-autotetraploid, and established autotetraploid Arabidopsis arenosa using new approaches for analysis of meiotic CO patterns in polyploids. We discover that crossover interference, the classical process responsible for patterning of COs in diploid meiosis, is defective in the neo-autotetraploid but robust in the established autotetraploid. The presented findings suggest that, initially, diploid-like interference fails to act effectively on multivalent pairing and accompanying pre-CO recombination interactions and that stable autopolyploid meiosis can emerge by evolution of a “supercharged” interference process, which can now act effectively on such configurations. Thus, the basic interference mechanism responsible for simplifying CO patterns along chromosomes in diploid meiosis has evolved the capability to also simplify CO patterns among chromosomes in autopolyploids, thereby promoting bivalent formation. We further show that evolution of stable autotetraploidy preadapts meiosis to higher ploidy, which in turn has interesting mechanistic and evolutionary implications.

•

In a neo-autotetraploid, aberrant crossover interference confers aberrant meiosis

•

In a stable autotetraploid, regular crossover interference confers regular meiosis

•

Crossover and synaptic patterns point to evolution of “supercharged” interference

•

Accordingly, evolution of stable autotetraploidy preadapts to higher ploidies

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.

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.

High-Frequency Homologous Recombination Occurred Preferentially in Populus

Homologous recombination (HR), the most significant event in meiosis, has important implications for genetic diversity and evolution in organisms. Heteroduplex DNA (hDNA), the product of HR, can be captured by artificially induced chromosome doubling during the development of the embryo sac to inhibit postmeiotic segregation, subsequently, and hDNAs are directly detected using codominant simple sequence repeat (SSR) markers. In the present study, two hybrid triploid populations derived from doubling the chromosomes of the embryo sac induced by high temperature in Populus tomentosa served as starting materials. Eighty-seven, 62, and 79 SSR markers on chromosomes 01, 04, and 19, respectively, that were heterozygous in the maternal parent and different from the paternal parent were screened to detect and characterize the hDNA in P. tomentosa. The results showed that the hDNA frequency patterns on chromosomes changed slightly when the number of SSR primers increased. The highest hDNA frequency occurred at the adjacent terminal on chromosomes, which was slightly higher than those at the terminals in the two genotypic individuals, and the hDNA frequency gradually decreased as the locus-centromere distance decreased. With the increase in the number of SSR markers employed for detection, the number of recombination events (REs) detected significantly increased. In regions with high methylation or long terminal repeat (LTR) retrotransposon enrichment, the frequency of hDNA was low, and high frequencies were observed in regions with low sequence complexity and high gene density. High-frequency recombination occurring at high gene density regions strongly affected the association between molecular markers and quantitative trait loci (QTLs), which was an important factor contributing to the difficulty encountered by MAS in achieving the expected breeding results.

Genome-wide variability in recombination activity is associated with meiotic chromatin organization

Recombination enables reciprocal exchange of genomic information between parental chromosomes and successful segregation of homologous chromosomes during meiosis. Errors in this process lead to negative health outcomes, whereas variability in recombination rate affects genome evolution. In mammals, most crossovers occur in hotspots defined by PRDM9 motifs, although PRDM9 binding peaks are not all equally hot. We hypothesize that dynamic patterns of meiotic genome folding are linked to recombination activity. We apply an integrative bioinformatics approach to analyze how three-dimensional (3D) chromosomal organization during meiosis relates to rates of double-strand-break (DSB) and crossover (CO) formation at PRDM9 binding peaks. We show that active, spatially accessible genomic regions during meiotic prophase are associated with DSB-favored loci, which further adopt a transient locally active configuration in early prophase. Conversely, crossover formation is depleted among DSBs in spatially accessible regions during meiotic prophase, particularly within gene bodies. We also find evidence that active chromatin regions have smaller average loop sizes in mammalian meiosis. Collectively, these findings establish that differences in chromatin architecture along chromosomal axes are associated with variable recombination activity. We propose an updated framework describing how 3D organization of brush-loop chromosomes during meiosis may modulate recombination.

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.

Recombination Rate Variation and Infrequent Sex Influence Genetic Diversity in Chlamydomonas reinhardtii

Recombination confers a major evolutionary advantage by breaking up linkage disequilibrium between harmful and beneficial mutations, thereby facilitating selection. However, in species that are only periodically sexual, such as many microbial eukaryotes, the realized rate of recombination is also affected by the frequency of sex, meaning that infrequent sex can increase the effects of selection at linked sites despite high recombination rates. Despite this, the rate of sex of most facultatively sexual species is unknown. Here, we use genomewide patterns of linkage disequilibrium to infer fine-scale recombination rate variation in the genome of the facultatively sexual green alga Chlamydomonas reinhardtii. We observe recombination rate variation of up to two orders of magnitude and find evidence of recombination hotspots across the genome. Recombination rate is highest flanking genes, consistent with trends observed in other nonmammalian organisms, though intergenic recombination rates vary by intergenic tract length. We also find a positive relationship between nucleotide diversity and physical recombination rate, suggesting a widespread influence of selection at linked sites in the genome. Finally, we use estimates of the effective rate of recombination to calculate the rate of sex that occurs in natural populations, estimating a sexual cycle roughly every 840 generations. We argue that the relatively infrequent rate of sex and large effective population size creates a population genetic environment that increases the influence of selection on linked sites across the genome.

Chasing breeding footprints through structural variations in Cucumis melo and wild relatives

Cucumis melo (melon or muskmelon) is an important crop in the family of the Cucurbitaceae. Melon is cross pollinated and domesticated at several locations throughout the breeding history, resulting in highly diverse genetic structure in the germplasm. Yet, the relations among the groups and cultivars are still incomplete. We shed light on the melonbreeding history, analyzing structural variations ranging from 50 bp up to 100 kb, identified from whole genome sequences of 100 selected melon accessions and wild relatives. Phylogenetic trees based on SV types completely resolve cultivars and wild accessions into two monophyletic groups and clustering of cultivars largely correlates with their geographic origin. Taking into account morphology, we found six mis-categorized cultivars. Unique inversions are more often shared between cultivars, carrying advantageous genes and do not directly originate from wild species. Approximately 60% of the inversion breaks carry a long poly A/T motif, and following observations in other plant species, suggest that inversions in melon likely resulted from meiotic recombination events. We show that resistance genes in the linkage V region are expanded in the cultivar genomes compared to wild relatives. Furthermore, particular agronomic traits such as fruit ripening, fragrance, and stress response are specifically selected for in the melon subspecies. These results represent distinctive footprints of selective breeding that shaped today's melon. The sequences and genomic relations between land races, wild relatives, and cultivars will serve the community to identify genetic diversity, optimize experimental designs, and enhance crop development.