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

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.

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.

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.

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.

Understanding and exploiting uniparental genome elimination in plants: insights from Arabidopsis thaliana

Uniparental genome elimination (UGE) refers to the preferential exclusion of one set of the parental chromosome complement during embryogenesis following successful fertilization, giving rise to uniparental haploid progeny. This artificially induced phenomenon was documented as one of the consequences of distant (wide) hybridization in plants. Ten decades since its discovery, attempts to unravel the molecular mechanism behind this process remained elusive due to a lack of genetic tools and genomic resources in the species exhibiting UGE. Hence, its successful adoption in agronomic crops for in planta (in vivo) haploid production remains implausible. Recently, Arabidopsis thaliana has emerged as a model system to unravel the molecular basis of UGE. It is now possible to simulate the genetic consequences of distant crosses in an A. thaliana intraspecific cross by a simple modification of centromeres, via the manipulation of the centromere-specific histone H3 variant gene, CENH3. Thus, the experimental advantages conferred by A. thaliana have been used to elucidate and exploit the benefits of UGE in crop breeding. In this review, we discuss developments and prospects of CENH3 gene-mediated UGE and other in planta haploid induction strategies to illustrate its potential in expediting plant breeding and genetics in A. thaliana and other model plants.

A Central Role for Genetics in Plant Biology

This article describes my involvement in the development of genetics as an essential tool in the integrated study of plant biology. My research comes from a strong background in plant genetics based on my education as a plant breeder at Wageningen University and collaborations with plant physiologists and molecular geneticists in Wageningen and the wider scientific community. It initially involved the isolation and physiological characterization of mutants defective in biosynthesis or mode of action of plant hormones, photoreceptors and traits such as flowering time in both Arabidopsis and tomato. I also generated a genetic map of Arabidopsis. Subsequently, the exploitation of natural variation became a main area of interest, including the molecular identification of underlying genetic differences. The integration of various disciplines and the adoption of Arabidopsis as a main model species contributed strongly to the impressive progress in our knowledge of plant biology over the past 40 years.

Brassica carinata genome characterization clarifies U's triangle model of evolution and polyploidy in Brassica

Ethiopian mustard (Brassica carinata) in the Brassicaceae family possesses many excellent agronomic traits. Here, the high-quality genome sequence of B. carinata is reported. Characterization revealed a genome anchored to 17 chromosomes with a total length of 1.087 Gb and an N50 scaffold length of 60 Mb. Repetitive sequences account for approximately 634 Mb or 58.34% of the B. carinata genome. Notably, 51.91% of 97,149 genes are confined to the terminal 20% of chromosomes as a result of the expansion of repeats in pericentromeric regions. Brassica carinata shares one whole-genome triplication event with the five other species in U's triangle, a classic model of evolution and polyploidy in Brassica. Brassica carinata was deduced to have formed ∼0.047 Mya, which is slightly earlier than B. napus but later than B. juncea. Our analysis indicated that the relationship between the two subgenomes (BcaB and BcaC) is greater than that between other two tetraploid subgenomes (BjuB and BnaC) and their respective diploid parents. RNA-seq datasets and comparative genomic analysis were used to identify several key genes in pathways regulating disease resistance and glucosinolate metabolism. Further analyses revealed that genome triplication and tandem duplication played important roles in the expansion of those genes in Brassica species. With the genome sequencing of B. carinata completed, the genomes of all six Brassica species in U's triangle are now resolved. The data obtained from genome sequencing, transcriptome analysis, and comparative genomic efforts in this study provide valuable insights into the genome evolution of the six Brassica species in U's triangle.

Comparative transcriptomic analysis of thermally stressed Arabidopsis thaliana meiotic recombination mutants

BACKGROUND

Meiosis is a specialized cell division that underpins sexual reproduction in most eukaryotes. During meiosis, interhomolog meiotic recombination facilitates accurate chromosome segregation and generates genetic diversity by shuffling parental alleles in the gametes. The frequency of meiotic recombination in Arabidopsis has a U-shaped curve in response to environmental temperature, and is dependent on the Type I, crossover (CO) interference-sensitive pathway. The mechanisms that modulate recombination frequency in response to temperature are not yet known.

RESULTS

In this study, we compare the transcriptomes of thermally-stressed meiotic-stage anthers from msh4 and mus81 mutants that mediate the Type I and Type II meiotic recombination pathways, respectively. We show that heat stress reduces the number of expressed genes regardless of genotype. In addition, msh4 mutants have a distinct gene expression pattern compared to mus81 and wild type controls. Interestingly, ASY1, which encodes a HORMA domain protein that is a component of meiotic chromosome axes, is up-regulated in wild type and mus81 but not in msh4. In addition, SDS the meiosis-specific cyclin-like gene, DMC1 the meiosis-specific recombinase, SYN1/REC8 the meiosis-specific cohesion complex component, and SWI1 which functions in meiotic sister chromatid cohesion are up-regulated in all three genotypes. We also characterize 51 novel, previously unannotated transcripts, and show that their promoter regions are associated with A-rich meiotic recombination hotspot motifs.

CONCLUSIONS

Our transcriptomic analysis of msh4 and mus81 mutants enhances our understanding of how the Type I and Type II meiotic CO pathway respond to environmental temperature stress and might provide a strategy to manipulate recombination levels in plants.

A phased genome based on single sperm sequencing reveals crossover pattern and complex relatedness in tea plants

For diploid organisms that are highly heterozygous, a phased haploid genome can greatly aid in functional genomic, population genetic and breeding studies. Based on the genome sequencing of 135 single sperm cells of the elite tea cultivar 'Fudingdabai', we herein phased the genome of Camellia sinensis, one of the most popular beverage crops worldwide. High-resolution genetic and recombination maps of Fudingdabai were constructed, which revealed that crossover (CO) positions were frequently located in the 5' and 3' ends of annotated genes, while CO distributions across the genome were random. The low CO frequency in tea can be explained by strong CO interference, and CO simulation revealed the proportion of interference insensitive CO ranged from 5.2% to 11.7%. We furthermore developed a method to infer the relatedness between tea accessions and detected complex kinship and genetic signatures of 106 tea accessions. Among them, 59 accessions were closely related with Fudingdabai and 31 of them were first-degree relatives. We additionally identified genes displaying allele specific expression patterns between the two haplotypes of Fudingdabai and genes displaying significantly differential expression levels between Fudingdabai and other haplotypes. These results lay the foundation for further investigation of genetic and epigenetic factors underpinning the regulation of gene expression and provide insights into the evolution of tea plants as well as a valuable genetic resource for future breeding efforts.

Sequence analysis of wheat subtelomeres reveals a high polymorphism among homoeologous chromosomes

Bread wheat, Triticum aestivum L., is one of the most important crops in the world. Understanding its genome organization (allohexaploid; AABBDD; 2n = 6x = 42) is essential for geneticists and plant breeders. Particularly, the knowledge of how homologous chromosomes (equivalent chromosomes from the same genome) specifically recognize each other to pair at the beginning of meiosis, the cellular process to generate gametes in sexually reproducing organisms, is fundamental for plant breeding and has a big influence on the fertility of wheat plants. Initial homologous chromosome interactions contribute to specific recognition and pairing between homologues at the onset of meiosis. Understanding the molecular basis of these critical processes can help to develop genetic tools in a breeding context to promote interspecific chromosome associations in hybrids or interspecific genetic crosses to facilitate the transfer of desirable agronomic traits from related species into a crop like wheat. The terminal regions of chromosomes, which include telomeres and subtelomeres, participate in chromosome recognition and pairing. We present a detailed molecular analysis of subtelomeres of wheat chromosome arms 1AS, 4AS, 7AS, 7BS and 7DS. Results showed a high polymorphism in the subtelomeric region among homoeologues (equivalent chromosomes from related genomes) for all the features analyzed, including genes, transposable elements, repeats, GC content, predicted CpG islands, recombination hotspots and targeted sequence motifs for relevant DNA-binding proteins. These polymorphisms might be the molecular basis for the specificity of homologous recognition and pairing in initial chromosome interactions at the beginning of meiosis in wheat.

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

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

Extended TAQing system for large-scale plant genome reorganization

We previously developed a large-scale genome restructuring technology called the TAQing system. It can induce genomic rearrangements by introducing transient and conditional formation of DNA double-strand breaks (DSBs) via heat activation of a restriction enzyme TaqI, which can cleave DNA at 5'-TCGA-3' sequences in the genome at higher temperatures (37-42°C). Such heat treatment sometimes confers lethal damage in certain plant species and TaqI cannot induce rearrangements in AT-rich regions. To overcome such problems we developed an extended TAQing (Ex-TAQing) system, which enables the use of a wider range of restriction enzymes active at standard plant-growing temperatures. We established the Ex-TAQing system using MseI that can efficiently cleave DNA at room temperature (at temperatures ranging from 22 to 25°C) and the 5'-TTAA-3' sequence which is highly abundant in the Arabidopsis genome. A synthetic intron-spanning MseI gene, which was placed downstream of a heat-shock-inducible promoter, was conditionally expressed upon milder heat treatment (33°C) to generate DSBs in Arabidopsis chromosomes. Genome resequencing revealed various types of genomic rearrangements, including copy number variations, translocation and direct end-joining at MseI cleavage sites. The Ex-TAQing system could induce large-scale rearrangements in diploids more frequently (17.4%, n = 23) than the standard TAQing system. The application of this system to tetraploids generated several strains with chromosomal rearrangements associated with beneficial phenotypes, such as high salinity stress tolerance and hypersensitivity to abscisic acid. We have developed the Ex-TAQing system, allowing more diverse patterns of genomic rearrangements, by employing various types of endonucleases and have opened a way to expand the capacity for artificial genome reorganization.