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

Meiotic recombination has been deeply characterized in a few model species only, notably in the budding yeast Saccharomyces cerevisiae. Interestingly, most members of the ZMM pathway that implements meiotic crossover interference in S. cerevisiae have been lost in Lachancea yeast species after the divergence of Lachancea kluyveri from the rest of the clade. This suggests major differences in the control of crossover distribution. After investigating meiosis in L. kluyveri, we determined the meiotic recombination landscape of Lachancea waltii and identified several characteristics that should help understand better the underlying mechanisms. Such characteristics include systematic regions of loss of heterozygosity (LOH) in L. waltii hybrids, compatible with dysregulated Spo11-mediated DNA double strand breaks (DSB) independently of meiosis. They include a higher recombination rate in L. waltii than in L. kluyveri despite the lack of multiple ZMM pro-crossover factors. L. waltii exhibits an elevated frequency of zero-crossover bivalents as L. kluyveri but opposite to S. cerevisiae. L. waltii gene conversion tracts lengths are comparable to those observed in S. cerevisiae and shorter than in L. kluyveri despite the lack of Mlh2, a factor limiting conversion tracts size in S. cerevisiae. L. waltii recombination hotspots are 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, in line with the loss of several ZMM genes, we found only residual crossover interference in L. waltii likely coming from the modest interference existing between recombination precursors.

Genetic structure and evolutionary diversity of mating-type (MAT) loci in Hypsizygus marmoreus

The mating compatibility in fungi is generally governed by genes located within a single or two unlinked mating type (MAT) loci. Hypsizygus marmoreus is an edible mushroom in the order Agaricales with a tetrapolar system, which contains two unlinked MAT loci-homeodomain (HD) transcription factor genes and pheromone/pheromone receptor genes (P/R). In this study, we analyzed the genetic structure and diversity of MAT loci in tetrapolar system of H. marmoreus through sequencing of 54 heterokaryon and 8 homokaryon strains. Although within the HD loci, the gene order was conserved, the gene contents were variable, and the HD loci haplotypes were further classified into four types. By analyzing the structure, phylogeny, and the HD transmissibility based on the progeny of these four HD mating-type loci types, we found that they were heritable and tightly linked at the HD loci. The P/R loci genes were found to comprise three pheromone receptors, three pheromones, and two pheromone receptor-like genes. Intra- and inter-specific phylogenetic analyses of pheromone receptors revealed that the STE3 genes were divided into three groups, and we thus theorize that they diverged before speciation. Comparative analysis of the MAT regions among 73 Basidiomycete species indicated that the diversity of HD and P/R loci in Agaricales and Boletales may contribute to mating compatibility. The number of HD genes were not correlated with the tetrapolar or bipolar systems. In H. marmoreus, the expression levels of these genes at HD and P/R loci of compatible strains were found higher than in those of homonuclear/homokaryotic strains, indicating that these mating genes acted as switches for mating processes. Further collinear analysis of HD loci in interspecific species found that HD loci contains conserved recombination hotspots showing major rearrangements in Coprinopsis cinerea and Schizophyllum commune, suggesting different mechanisms for evolution of physically linked MAT loci in these groups. It seems likely that gene rearrangements are common in Agaricales fungi around HD loci. Together, our study provides insights into the genomic basis of mating compatibility in H. marmoreus.

Aborting meiosis allows recombination in sterile diploid yeast hybrids

Hybrids between diverged lineages contain novel genetic combinations but an impaired meiosis often makes them evolutionary dead ends. Here, we explore to what extent an aborted meiosis followed by a return-to-growth (RTG) promotes recombination across a panel of 20 Saccharomyces cerevisiae and S. paradoxus diploid hybrids with different genomic structures and levels of sterility. Genome analyses of 275 clones reveal that RTG promotes recombination and generates extensive regions of loss-of-heterozygosity in sterile hybrids with either a defective meiosis or a heavily rearranged karyotype, whereas RTG recombination is reduced by high sequence divergence between parental subgenomes. The RTG recombination preferentially arises in regions with low local heterozygosity and near meiotic recombination hotspots. The loss-of-heterozygosity has a profound impact on sexual and asexual fitness, and enables genetic mapping of phenotypic differences in sterile lineages where linkage analysis would fail. We propose that RTG gives sterile yeast hybrids access to a natural route for genome recombination and adaptation.

Deep learning identifies and quantifies recombination hotspot determinants

Recombination is one of the essential genetic processes for sexually reproducing organisms, which can happen more frequently in some regions, called recombination hotspots. Although several factors, such as PRDM9 binding motifs, are known to be related to the hotspots, their contributions to the recombination hotspots have not been quantified, and other determinants are yet to be elucidated. Here, we develop a computational method, RHSNet, based on deep learning and signal processing, to identify and quantify the hotspot determinants in a purely data-driven manner, utilizing datasets from various studies, populations, sexes, and species. In addition to being able to identify hotspot regions and the well-known determinants accurately, RHSNet is sensitive to the difference between different PRDM9 alleles and different sexes, and can generalize to PRDM9-lacking species. The cross-sex, cross-population, and cross-species studies suggest that the proposed method has the potential to identify and quantify the evolutionary determinant motifs.

Epigenetic Marks and Variation of Sequence-Based Information Along Genomic Regions Are Predictive of Recombination Hot/Cold Spots in Saccharomyces cerevisiae

Characterization and identification of recombination hotspots provide important insights into the mechanism of recombination and genome evolution. In contrast with existing sequence-based models for predicting recombination hotspots which were defined in a ORF-based manner, here, we first defined recombination hot/cold spots based on public high-resolution Spo11-oligo-seq data, then characterized them in terms of DNA sequence and epigenetic marks, and finally presented classifiers to identify hotspots. We found that, in addition to some previously discovered DNA-based features like GC-skew, recombination hotspots in yeast can also be characterized by some remarkable features associated with DNA physical properties and shape. More importantly, by using DNA-based features and several epigenetic marks, we built several classifiers to discriminate hotspots from coldspots, and found that SVM classifier performs the best with an accuracy of ∼92%, which is also the highest among the models in comparison. Feature importance analysis combined with prediction results show that epigenetic marks and variation of sequence-based features along the hotspots contribute dominantly to hotspot identification. By using incremental feature selection method, an optimal feature subset that consists of much less features was obtained without sacrificing prediction accuracy.

PRDM9 directs recombination hotspots to prevent competition with meiotic transcription

PRDM9 drives recombination hotspots in some mammals, including mice and apes. Non-functional orthologs of PRDM9 are present in a wide variety of vertebrates, but why it is functionally maintained in some lineages is not clear. One possible explanation is that PRDM9 plays a role in ensuring that meiosis is successful. During meiosis, available DNA may be a limiting resource given the tight packaging of chromosomes and could lead to competition between two key processes: meiotic transcription and recombination. Here we explore this potential competition and the role that PRDM9 might play in their interaction. Leveraging existing mouse genomic data, we use resampling schemes that simulate shuffled features along the genome and models that account for the rarity of features in the genome, to test if PRDM9 influences interactions between recombination hotspots and meiotic transcription in a whole genome framework. We also explored possible DNA sequence motifs associated to clusters of hotspots not tied to transcription or PRDM9. We find evidence of competition between meiotic transcription and recombination, with PRDM9 appearing to relocate recombination to avoid said conflict. We also find that retrotransposons may be playing a role in directing hotspots in the absence of other factors.

Evolutionary implications of recombination differences across diverging populations of Anopheles

Recombination is one of the main evolutionary mechanisms responsible for changing the genomic architecture of populations; and in essence, it is the main mechanism by which novel combinations of alleles, haplotypes, are formed. A clear picture that has emerged across study systems is that recombination is highly variable, even among closely related species. However, it is only until very recently that we have started to understand how recombination variation between populations of the same species impact genetic diversity and divergence. Here, we used whole-genome sequence data to build fine-scale recombination maps for nine populations within two species of Anopheles, Anopheles gambiae and Anopheles coluzzii. The genome-wide recombination averages were on the same order of magnitude for all populations except one. Yet, we identified significant differences in fine-scale recombination rates among all population comparisons. We report that effective population sizes, and presence of a chromosomal inversion has major contribution to recombination rate variation along the genome and across populations. We identified over 400 highly variable recombination hotspots across all populations, where only 9.6% are shared between two or more populations. Additionally, our results are consistent with recombination hotspots contributing to both genetic diversity and absolute divergence (dxy) between populations and species of Anopheles. However, we also show that recombination has a small impact on population genetic differentiation as estimated with FST. The minimal impact that recombination has on genetic differentiation across populations represents the first empirical evidence against recent theoretical work suggesting that variation in recombination along the genome can mask or impair our ability to detect signatures of selection. Our findings add new understanding to how recombination rates vary within species, and how this major evolutionary mechanism can maintain and contribute to genetic variation and divergence within a prominent malaria vector.