Meiotic recombination initiation in and around retrotransposable elements in Saccharomyces cerevisiae

Meiotic prophase length modulates Tel1-dependent DNA double-strand break interference

During meiosis, genetic recombination is initiated by the formation of many DNA double-strand breaks (DSBs) catalysed by the evolutionarily conserved topoisomerase-like enzyme, Spo11, in preferred genomic sites known as hotspots. DSB formation activates the Tel1/ATM DNA damage responsive (DDR) kinase, locally inhibiting Spo11 activity in adjacent hotspots via a process known as DSB interference. Intriguingly, in S. cerevisiae, over short genomic distances (<15 kb), Spo11 activity displays characteristics of concerted activity or clustering, wherein the frequency of DSB formation in adjacent hotspots is greater than expected by chance. We have proposed that clustering is caused by a limited number of sub-chromosomal domains becoming primed for DSB formation. Here, we provide evidence that DSB clustering is abolished when meiotic prophase timing is extended via deletion of the NDT80 transcription factor. We propose that extension of meiotic prophase enables most cells, and therefore most chromosomal domains within them, to reach an equilibrium state of similar Spo11-DSB potential, reducing the impact that priming has on estimates of coincident DSB formation. Consistent with this view, when Tel1 is absent but Ndt80 is present and thus cells are able to rapidly exit meiotic prophase, genome-wide maps of Spo11-DSB formation are skewed towards pericentromeric regions and regions that load pro-DSB factors early-revealing regions of preferential priming-but this effect is abolished when NDT80 is deleted. Our work highlights how the stochastic nature of Spo11-DSB formation in individual cells within the limited temporal window of meiotic prophase can cause localised DSB clustering-a phenomenon that is exacerbated in tel1Δ cells due to the dual roles that Tel1 has in DSB interference and meiotic prophase checkpoint control.

Dna2 removes toxic ssDNA-RPA filaments generated from meiotic recombination-associated DNA synthesis

During the repair of DNA double-strand breaks (DSBs), de novo synthesized DNA strands can displace the parental strand to generate single-strand DNAs (ssDNAs). Many programmed DSBs and thus many ssDNAs occur during meiosis. However, it is unclear how these ssDNAs are removed for the complete repair of meiotic DSBs. Here, we show that meiosis-specific depletion of Dna2 (dna2-md) results in an abundant accumulation of RPA and an expansion of RPA from DSBs to broader regions in Saccharomyces cerevisiae. As a result, DSB repair is defective and spores are inviable, although the levels of crossovers/non-crossovers seem to be unaffected. Furthermore, Dna2 induction at pachytene is highly effective in removing accumulated RPA and restoring spore viability. Moreover, the depletion of Pif1, an activator of polymerase δ required for meiotic recombination-associated DNA synthesis, and Pif1 inhibitor Mlh2 decreases and increases RPA accumulation in dna2-md, respectively. In addition, blocking DNA synthesis during meiotic recombination dramatically decreases RPA accumulation in dna2-md. Together, our findings show that meiotic DSB repair requires Dna2 to remove ssDNA-RPA filaments generated from meiotic recombination-associated DNA synthesis. Additionally, we showed that Dna2 also regulates DSB-independent RPA distribution.

UGDR: a generic pipeline to detect recombined regions in polyploid and complex hybrid yeast genomes

MOTIVATION

In eukaryotes, homologous recombination between the parental genomes frequently occurs during the evolutionary conserved process of meiosis, generating the genetic diversity transmitted by the gametes. The genome-wide determination of the frequency and location of the recombination events can now be efficiently performed by genotyping the offspring's polymorphic markers. However, genotyping recombination in complex hybrid genomes with existing methods remains challenging because of their strain and ploidy specificity and the degree of diversity and complexity of the parental genomes, especially in [Formula: see text] polyploids.

RESULTS

We present UGDR, a pipeline to genotype the polymorphisms of complex hybrid yeast genomes. It is based on optimal mapping strategies of NGS reads, comparative analyses of the allelic ratio variation and read depth coverage. We tested the UGDR pipeline with sequencing reads from recombined hybrid diploid yeast strains and various clinical strains exhibiting different degrees of ploidy. UGDR allows to plot the markers distribution and recombination profile per chromosome.

CONCLUSION

UGDR detects and plots recombination events in haploids and polyploid yeasts, which facilitates the discovery and understanding of the yeast genetic recombination map and identify new out-performing recombinants.

A widespread inversion polymorphism conserved among Saccharomyces species is caused by recurrent homogenization of a sporulation gene family

Saccharomyces genomes are highly collinear and show relatively little structural variation, both within and between species of this yeast genus. We investigated the only common inversion polymorphism known in S. cerevisiae, which affects a 24-kb 'flip/flop' region containing 15 genes near the centromere of chromosome XIV. The region exists in two orientations, called reference (REF) and inverted (INV). Meiotic recombination in this region is suppressed in crosses between REF and INV orientation strains such as the BY x RM cross. We find that the inversion polymorphism is at least 17 million years old because it is conserved across the genus Saccharomyces. However, the REF and INV isomers are not ancient alleles but are continually being re-created by re-inversion of the region within each species. Inversion occurs due to continual homogenization of two almost identical 4-kb sequences that form an inverted repeat (IR) at the ends of the flip/flop region. The IR consists of two pairs of genes that are specifically and strongly expressed during the late stages of sporulation. We show that one of these gene pairs, YNL018C/YNL034W, codes for a protein that is essential for spore formation. YNL018C and YNL034W are the founder members of a gene family, Centroid, whose members in other Saccharomycetaceae species evolve fast, duplicate frequently, and are preferentially located close to centromeres. We tested the hypothesis that Centroid genes are a meiotic drive system, but found no support for this idea.

Genomic diversity, chromosomal rearrangements, and interspecies hybridization in the Ogataea polymorpha species complex

The methylotrophic yeast Ogataea polymorpha has long been a useful system for recombinant protein production, as well as a model system for methanol metabolism, peroxisome biogenesis, thermotolerance, and nitrate assimilation. It has more recently become an important model for the evolution of mating-type switching. Here, we present a population genomics analysis of 47 isolates within the O. polymorpha species complex, including representatives of the species O. polymorpha, Ogataea parapolymorpha, Ogataea haglerorum, and Ogataea angusta. We found low levels of nucleotide sequence diversity within the O. polymorpha species complex and identified chromosomal rearrangements both within and between species. In addition, we found that one isolate is an interspecies hybrid between O. polymorpha and O. parapolymorpha and present evidence for loss of heterozygosity following hybridization.

Meiosis and beyond - understanding the mechanistic and evolutionary processes shaping the germline genome

The separation of germ cell populations from the soma is part of the evolutionary transition to multicellularity. Only genetic information present in the germ cells will be inherited by future generations, and any molecular processes affecting the germline genome are therefore likely to be passed on. Despite its prevalence across taxonomic kingdoms, we are only starting to understand details of the underlying micro-evolutionary processes occurring at the germline genome level. These include segregation, recombination, mutation and selection and can occur at any stage during germline differentiation and mitotic germline proliferation to meiosis and post-meiotic gamete maturation. Selection acting on germ cells at any stage from the diploid germ cell to the haploid gametes may cause significant deviations from Mendelian inheritance and may be more widespread than previously assumed. The mechanisms that affect and potentially alter the genomic sequence and allele frequencies in the germline are pivotal to our understanding of heritability. With the rise of new sequencing technologies, we are now able to address some of these unanswered questions. In this review, we comment on the most recent developments in this field and identify current gaps in our knowledge.

Genome assembly of the basket willow, Salix viminalis, reveals earliest stages of sex chromosome expansion

BACKGROUND

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

RESULTS

Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

CONCLUSIONS

Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

Genome assembly of the basket willow, Salix viminalis, reveals earliest stages of sex chromosome expansion

BACKGROUND

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

RESULTS

Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

CONCLUSIONS

Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

Advances Towards How Meiotic Recombination Is Initiated: A Comparative View and Perspectives for Plant Meiosis Research

Meiosis is an essential cell-division process for ensuring genetic diversity across generations. Meiotic recombination ensures the accuracy of genetic interchange between homolous chromosomes and segregation of parental alleles. Programmed DNA double-strand breaks (DSBs), catalyzed by the evolutionarily conserved topoisomerase VIA (a subunit of the archaeal type II DNA topoisomerase)-like enzyme Spo11 and several other factors, is a distinctive feature of meiotic recombination initiation. The meiotic DSB formation and its regulatory mechanisms are similar among species, but certain aspects are distinct. In this review, we introduced the cumulative knowledge of the plant proteins crucial for meiotic DSB formation and technical advances in DSB detection. We also summarized the genome-wide DSB hotspot profiles for different model organisms. Moreover, we highlighted the classical views and recent advances in our knowledge of the regulatory mechanisms that ensure the fidelity of DSB formation, such as multifaceted kinase-mediated phosphorylation and the consequent high-dimensional changes in chromosome structure. We provided an overview of recent findings concerning DSB formation, distribution and regulation, all of which will help us to determine whether meiotic DSB formation is evolutionarily conserved or varies between plants and other organisms.

Integrative Meta-Assembly Pipeline (IMAP): Chromosome-level genome assembler combining multiple de novo assemblies

Background

Genomic data have become major resources to understand complex mechanisms at fine-scale temporal and spatial resolution in functional and evolutionary genetic studies, including human diseases, such as cancers. Recently, a large number of whole genomes of evolving populations of yeast (Saccharomyces cerevisiae W303 strain) were sequenced in a time-dependent manner to identify temporal evolutionary patterns. For this type of study, a chromosome-level sequence assembly of the strain or population at time zero is required to compare with the genomes derived later. However, there is no fully automated computational approach in experimental evolution studies to establish the chromosome-level genome assembly using unique features of sequencing data.

Methods and results

In this study, we developed a new software pipeline, the integrative meta-assembly pipeline (IMAP), to build chromosome-level genome sequence assemblies by generating and combining multiple initial assemblies using three de novo assemblers from short-read sequencing data. We significantly improved the continuity and accuracy of the genome assembly using a large collection of sequencing data and hybrid assembly approaches. We validated our pipeline by generating chromosome-level assemblies of yeast strains W303 and SK1, and compared our results with assemblies built using long-read sequencing and various assembly evaluation metrics. We also constructed chromosome-level sequence assemblies of S. cerevisiae strain Sigma1278b, and three commonly used fungal strains: Aspergillus nidulans A713, Neurospora crassa 73, and Thielavia terrestris CBS 492.74, for which long-read sequencing data are not yet available. Finally, we examined the effect of IMAP parameters, such as reference and resolution, on the quality of the final assembly of the yeast strains W303 and SK1.

Conclusions

We developed a cost-effective pipeline to generate chromosome-level sequence assemblies using only short-read sequencing data. Our pipeline combines the strengths of reference-guided and meta-assembly approaches. Our pipeline is available online at http://github.com/jkimlab/IMAP including a Docker image, as well as a Perl script, to help users install the IMAP package, including several prerequisite programs. Users can use IMAP to easily build the chromosome-level assembly for the genome of their interest.

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

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

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

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

Hotspots for Initiation of Meiotic Recombination

Homologous chromosomes must pair and recombine to ensure faithful chromosome segregation during meiosis, a specialized type of cell division that occurs in sexually reproducing eukaryotes. Meiotic recombination initiates by programmed induction of DNA double-strand breaks (DSBs) by the conserved type II topoisomerase-like enzyme SPO11. A subset of meiotic DSBs are resolved as crossovers, whereby reciprocal exchange of DNA occurs between homologous chromosomes. Importantly, DSBs are non-randomly distributed along eukaryotic chromosomes, forming preferentially in permissive regions known as hotspots. In many species, including plants, DSB hotspots are located within nucleosome-depleted regions. DSB localization is governed by interconnected factors, including cis-regulatory elements, transcription factor binding, and chromatin accessibility, as well as by higher-order chromosome architecture. The spatiotemporal control of DSB formation occurs within a specialized chromosomal structure characterized by sister chromatids organized into linear arrays of chromatin loops that are anchored to a proteinaceous axis. Although SPO11 and its partner proteins required for DSB formation are bound to the axis, DSBs occur preferentially within the chromatin loops, which supports the "tethered-loop/axis model" for meiotic recombination. In this mini review, we discuss insights gained from recent efforts to define and profile DSB hotspots at high resolution in eukaryotic genomes. These advances are deepening our understanding of how meiotic recombination shapes genetic diversity and genome evolution in diverse species.

Coevolution between transposable elements and recombination

One of the most striking patterns of genome structure is the tight, typically negative, association between transposable elements (TEs) and meiotic recombination rates. While this is a highly recurring feature of eukaryotic genomes, the mechanisms driving correlations between TEs and recombination remain poorly understood, and distinguishing cause versus effect is challenging. Here, we review the evidence for a relation between TEs and recombination, and discuss the underlying evolutionary forces. Evidence to date suggests that overall TE densities correlate negatively with recombination, but the strength of this correlation varies across element types, and the pattern can be reversed. Results suggest that heterogeneity in the strength of selection against ectopic recombination and gene disruption can drive TE accumulation in regions of low recombination, but there is also strong evidence that the regulation of TEs can influence local recombination rates. We hypothesize that TE insertion polymorphism may be important in driving within-species variation in recombination rates in surrounding genomic regions. Furthermore, the interaction between TEs and recombination may create positive feedback, whereby TE accumulation in non-recombining regions contributes to the spread of recombination suppression. Further investigation of the coevolution between recombination and TEs has important implications for our understanding of the evolution of recombination rates and genome structure.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

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

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

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

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

Sequencing Spo11 Oligonucleotides for Mapping Meiotic DNA Double-Strand Breaks in Yeast

Meiosis is a specialized form of cell division resulting in reproductive cells with a reduced, usually haploid, genome complement. A key step after premeiotic DNA replication is the occurrence of homologous recombination at multiple places throughout the genome, initiated with the formation of DNA double-strand breaks (DSBs) catalyzed by the topoisomerase-like protein Spo11. DSBs are distributed non-randomly in genomes, and understanding the mechanisms that shape this distribution is important for understanding how meiotic recombination influences heredity and genome evolution. Several methods exist for mapping where Spo11 acts. Of these, sequencing of Spo11-associated oligonucleotides (Spo11 oligos) is the most precise, specifying the locations of DNA breaks to the base pair. In this chapter we detail the steps involved in Spo11-oligo mapping in the SK1 strain of budding yeast Saccharomyces cerevisiae, from harvesting cells of highly synchronous meiotic cultures, through preparation of sequencing libraries, to the mapping pipeline used for processing the data.

Programming sites of meiotic crossovers using Spo11 fusion proteins

Meiotic recombination shapes the genetic diversity transmitted upon sexual reproduction. However, its non-random distribution along the chromosomes constrains the landscape of potential genetic combinations. For a variety of purposes, it is desirable to expand the natural repertoire by changing the distribution of crossovers in a wide range of eukaryotes. Toward this end, we report the local stimulation of meiotic recombination at a number of chromosomal sites by tethering the natural Spo11 protein to various DNA-binding modules: full-length DNA binding proteins, zinc fingers (ZFs), transcription activator-like effector (TALE) modules, and the CRISPR-Cas9 system. In the yeast Saccharomyces cerevisiae, each strategy is able to stimulate crossover frequencies in naturally recombination-cold regions. The binding and cleavage efficiency of the targeting Spo11 fusions (TSF) are variable, being dependent on the chromosomal regions and potential competition with endogenous factors. TSF-mediated genome interrogation distinguishes naturally recombination-cold regions that are flexible and can be warmed-up (gene promoters and coding sequences), from those that remain refractory (gene terminators and centromeres). These results describe new generic experimental strategies to increase the genetic diversity of gametes, which should prove useful in plant breeding and other applications.

Genomic features shaping the landscape of meiotic double-strand-break hotspots in maize

Meiotic recombination is the most important source of genetic variation in higher eukaryotes. It is initiated by formation of double-strand breaks (DSBs) in chromosomal DNA in early meiotic prophase. The DSBs are subsequently repaired, resulting in crossovers (COs) and noncrossovers (NCOs). Recombination events are not distributed evenly along chromosomes but cluster at recombination hotspots. How specific sites become hotspots is poorly understood. Studies in yeast and mammals linked initiation of meiotic recombination to active chromatin features present upstream from genes, such as absence of nucleosomes and presence of trimethylation of lysine 4 in histone H3 (H3K4me3). Core recombination components are conserved among eukaryotes, but it is unclear whether this conservation results in universal characteristics of recombination landscapes shared by a wide range of species. To address this question, we mapped meiotic DSBs in maize, a higher eukaryote with a large genome that is rich in repetitive DNA. We found DSBs in maize to be frequent in all chromosome regions, including sites lacking COs, such as centromeres and pericentromeric regions. Furthermore, most DSBs are formed in repetitive DNA, predominantly Gypsy retrotransposons, and only one-quarter of DSB hotspots are near genes. Genic and nongenic hotspots differ in several characteristics, and only genic DSBs contribute to crossover formation. Maize hotspots overlap regions of low nucleosome occupancy but show only limited association with H3K4me3 sites. Overall, maize DSB hotspots exhibit distribution patterns and characteristics not reported previously in other species. Understanding recombination patterns in maize will shed light on mechanisms affecting dynamics of the plant genome.

Genomic and chromatin features shaping meiotic double-strand break formation and repair in mice

The SPO11-generated DNA double-strand breaks (DSBs) that initiate meiotic recombination occur non-randomly across genomes, but mechanisms shaping their distribution and repair remain incompletely understood. Here, we expand on recent studies of nucleotide-resolution DSB maps in mouse spermatocytes. We find that trimethylation of histone H3 lysine 36 around DSB hotspots is highly correlated, both spatially and quantitatively, with trimethylation of H3 lysine 4, consistent with coordinated formation and action of both PRDM9-dependent histone modifications. In contrast, the DSB-responsive kinase ATM contributes independently of PRDM9 to controlling hotspot activity, and combined action of ATM and PRDM9 can explain nearly two-thirds of the variation in DSB frequency between hotspots. DSBs were modestly underrepresented in most repetitive sequences such as segmental duplications and transposons. Nonetheless, numerous DSBs form within repetitive sequences in each meiosis and some classes of repeats are preferentially targeted. Implications of these findings are discussed for evolution of PRDM9 and its role in hybrid strain sterility in mice. Finally, we document the relationship between mouse strain-specific DNA sequence variants within PRDM9 recognition motifs and attendant differences in recombination outcomes. Our results provide further insights into the complex web of factors that influence meiotic recombination patterns.

Variation of the meiotic recombination landscape and properties over a broad evolutionary distance in yeasts

Meiotic recombination is a major factor of genome evolution, deeply characterized in only a few model species, notably the yeast Saccharomyces cerevisiae. Consequently, little is known about variations of its properties across species. In this respect, we explored the recombination landscape of Lachancea kluyveri, a protoploid yeast species that diverged from the Saccharomyces genus more than 100 million years ago and we found striking differences with S. cerevisiae. These variations include a lower recombination rate, a higher frequency of chromosomes segregating without any crossover and the absence of recombination on the chromosome arm containing the sex locus. In addition, although well conserved within the Saccharomyces clade, the S. cerevisiae recombination hotspots are not conserved over a broader evolutionary distance. Finally and strikingly, we found evidence of frequent reversal of commitment to meiosis, resulting in return to mitotic growth after allele shuffling. Identification of this major but underestimated evolutionary phenomenon illustrates the relevance of exploring non-model species.

Meiotic recombination promotes accurate chromosome segregation and genetic diversity. To date, the mechanisms and rules lying behind recombination were dissected using model organisms such as the budding yeast Saccharomyces cerevisiae. To assess the conservation and variation of this process over a broad evolutionary distance, we explored the meiotic recombination landscape in Lachancea kluyveri, a budding yeast species that diverged from S. cerevisiae more than 100 million years ago. The meiotic recombination map we generated revealed that the meiotic recombination landscape and properties significantly vary across distantly related yeast species, raising the yet to confirm possibility that recombination hotspots conservation across yeast species depends on synteny conservation. Finally, the frequent meiotic reversions we observed led us to re-evaluate their evolutionary importance.

Recombination rate variation in mice from an isolated island

Recombination rate is a heritable trait that varies among individuals. Despite the major impact of recombination rate on patterns of genetic diversity and the efficacy of selection, natural variation in this phenotype remains poorly characterized. We present a comparison of genetic maps, sampling 1212 meioses, from a unique population of wild house mice (Mus musculus domesticus) that recently colonized remote Gough Island. Crosses to a mainland reference strain (WSB/EiJ) reveal pervasive variation in recombination rate among Gough Island mice, including subchromosomal intervals spanning up to 28% of the genome. In spite of this high level of polymorphism, the genomewide recombination rate does not significantly vary. In general, we find that recombination rate varies more when measured in smaller genomic intervals. Using the current standard genetic map of the laboratory mouse to polarize intervals with divergent recombination rates, we infer that the majority of evolutionary change occurred in one of the two tested lines of Gough Island mice. Our results confirm that natural populations harbour a high level of recombination rate polymorphism and highlight the disparities in recombination rate evolution across genomic scales.

Extensive Recombination of a Yeast Diploid Hybrid through Meiotic Reversion

In somatic cells, recombination between the homologous chromosomes followed by equational segregation leads to loss of heterozygosity events (LOH), allowing the expression of recessive alleles and the production of novel allele combinations that are potentially beneficial upon Darwinian selection. However, inter-homolog recombination in somatic cells is rare, thus reducing potential genetic variation. Here, we explored the property of S. cerevisiae to enter the meiotic developmental program, induce meiotic Spo11-dependent double-strand breaks genome-wide and return to mitotic growth, a process known as Return To Growth (RTG). Whole genome sequencing of 36 RTG strains derived from the hybrid S288c/SK1 diploid strain demonstrates that the RTGs are bona fide diploids with mosaic recombined genome, derived from either parental origin. Individual RTG genome-wide genotypes are comprised of 5 to 87 homozygous regions due to the loss of heterozygous (LOH) events of various lengths, varying between a few nucleotides up to several hundred kilobases. Furthermore, we show that reiteration of the RTG process shows incremental increases of homozygosity. Phenotype/genotype analysis of the RTG strains for the auxotrophic and arsenate resistance traits validates the potential of this procedure of genome diversification to rapidly map complex traits loci (QTLs) in diploid strains without undergoing sexual reproduction.

High-Resolution Global Analysis of the Influences of Bas1 and Ino4 Transcription Factors on Meiotic DNA Break Distributions in Saccharomyces cerevisiae

Meiotic recombination initiates with DNA double-strand breaks (DSBs) made by Spo11. In Saccharomyces cerevisiae, many DSBs occur in "hotspots" coinciding with nucleosome-depleted gene promoters. Transcription factors (TFs) stimulate DSB formation in some hotspots, but TF roles are complex and variable between locations. Until now, available data for TF effects on global DSB patterns were of low spatial resolution and confined to a single TF. Here, we examine at high resolution the contributions of two TFs to genome-wide DSB distributions: Bas1, which was known to regulate DSB activity at some loci, and Ino4, for which some binding sites were known to be within strong DSB hotspots. We examined fine-scale DSB distributions in TF mutant strains by deep sequencing oligonucleotides that remain covalently bound to Spo11 as a byproduct of DSB formation, mapped Bas1 and Ino4 binding sites in meiotic cells, evaluated chromatin structure around DSB hotspots, and measured changes in global messenger RNA levels. Our findings show that binding of these TFs has essentially no predictive power for DSB hotspot activity and definitively support the hypothesis that TF control of DSB numbers is context dependent and frequently indirect. TFs often affected the fine-scale distributions of DSBs within hotspots, and when seen, these effects paralleled effects on local chromatin structure. In contrast, changes in DSB frequencies in hotspots did not correlate with quantitative measures of chromatin accessibility, histone H3 lysine 4 trimethylation, or transcript levels. We also ruled out hotspot competition as a major source of indirect TF effects on DSB distributions. Thus, counter to prevailing models, roles of these TFs on DSB hotspot strength cannot be simply explained via chromatin "openness," histone modification, or compensatory interactions between adjacent hotspots.

A high-definition view of functional genetic variation from natural yeast genomes

The question of how genetic variation in a population influences phenotypic variation and evolution is of major importance in modern biology. Yet much is still unknown about the relative functional importance of different forms of genome variation and how they are shaped by evolutionary processes. Here we address these questions by population level sequencing of 42 strains from the budding yeast Saccharomyces cerevisiae and its closest relative S. paradoxus. We find that genome content variation, in the form of presence or absence as well as copy number of genetic material, is higher within S. cerevisiae than within S. paradoxus, despite genetic distances as measured in single-nucleotide polymorphisms being vastly smaller within the former species. This genome content variation, as well as loss-of-function variation in the form of premature stop codons and frameshifting indels, is heavily enriched in the subtelomeres, strongly reinforcing the relevance of these regions to functional evolution. Genes affected by these likely functional forms of variation are enriched for functions mediating interaction with the external environment (sugar transport and metabolism, flocculation, metal transport, and metabolism). Our results and analyses provide a comprehensive view of genomic diversity in budding yeast and expose surprising and pronounced differences between the variation within S. cerevisiae and that within S. paradoxus. We also believe that the sequence data and de novo assemblies will constitute a useful resource for further evolutionary and population genomics studies.