Multiple independent losses of crossover interference during yeast evolutionary history.

Meiotic recombination is essential for the accurate chromosome segregation and the generation of genetic diversity through crossover and gene conversion events. Although this process has been studied extensively in a few selected model species, understanding how its properties vary across species remains limited. For instance, the ancestral ZMM pathway that generates interference-dependent crossovers has undergone multiple losses throughout evolution, suggesting variations in the regulation of crossover formation. In this context, we first characterized the meiotic recombination landscape and properties of the Kluyveromyces lactis budding yeast. We then conducted a comprehensive analysis of 29,151 recombination events (19, 212 COs and 9, 939 NCOs) spanning 577 meioses in the five budding yeast species Saccharomyces cerevisiae, Saccharomyces paradoxus, Lachancea kluyveri, Lachancea waltii and K. lactis. Eventually, we found that the Saccharomyces yeasts displayed higher recombination rates compared to the non-Saccharomyces yeasts. In addition, bona fide crossover interference and associated crossover homeostasis were detected in the Saccharomyces species only, adding L. kluyveri and K. lactis to the list of budding yeast species that lost crossover interference. Finally, recombination hotspots, although highly conserved within the Saccharomyces yeasts are not conserved beyond the Saccharomyces genus. Overall, these results highlight great variability in the recombination landscape and properties through budding yeasts evolution.

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.

Loss of Heterozygosity Spectrum Depends on Ploidy Level in Natural Yeast Populations.

The appearance of genomic variations such as loss of heterozygosity (LOH) has a significant impact on phenotypic diversity observed in a population. Recent large-scale yeast population genomic surveys have shown a high frequency of these events in natural isolates and more particularly in polyploids. However, the frequency, extent, and spectrum of LOH in polyploid organisms have never been explored and are poorly characterized to date. Here, we accumulated 5,163 LOH events over 1,875 generations in 76 mutation accumulation (MA) lines comprising nine natural heterozygous diploid, triploid, and tetraploid natural S. cerevisiae isolates from different ecological and geographical origins. We found that the rate and spectrum of LOH are variable across ploidy levels. Of the total accumulated LOH events, 8.5%, 21%, and 70.5% of them were found in diploid, triploid, and tetraploid MA lines, respectively. Our results clearly show that the frequency of generated LOH events increases with ploidy level. In fact, the cumulative LOH rates were estimated to be 9.3 × 10-3, 2.2 × 10-2, and 8.4 × 10-2 events per division for diploids, triploids, and tetraploids, respectively. In addition, a clear bias toward the accumulation of interstitial and short LOH tracts is observed in triploids and tetraploids compared with diploids. The variation of the frequency and spectrum of LOH events across ploidy level could be related to the genomic instability, characterizing higher ploidy isolates.

ZMIZ1 Variants Cause a Syndromic Neurodevelopmental Disorder.

ZMIZ1 is a coactivator of several transcription factors, including p53, the androgen receptor, and NOTCH1. Here, we report 19 subjects with intellectual disability and developmental delay carrying variants in ZMIZ1. The associated features include growth failure, feeding difficulties, microcephaly, facial dysmorphism, and various other congenital malformations. Of these 19, 14 unrelated subjects carried de novo heterozygous single-nucleotide variants (SNVs) or single-base insertions/deletions, 3 siblings harbored a heterozygous single-base insertion, and 2 subjects had a balanced translocation disrupting ZMIZ1 or involving a regulatory region of ZMIZ1. In total, we identified 13 point mutations that affect key protein regions, including a SUMO acceptor site, a central disordered alanine-rich motif, a proline-rich domain, and a transactivation domain. All identified variants were absent from all available exome and genome databases. In vitro, ZMIZ1 showed impaired coactivation of the androgen receptor. In vivo, overexpression of ZMIZ1 mutant alleles in developing mouse brains using in utero electroporation resulted in abnormal pyramidal neuron morphology, polarization, and positioning, underscoring the importance of ZMIZ1 in neural development and supporting mutations in ZMIZ1 as the cause of a rare neurodevelopmental syndrome.

Pan-transcriptome reveals a large accessory genome contribution to gene expression variation in yeast.

Gene expression is an essential step in the translation of genotypes into phenotypes. However, little is known about the transcriptome architecture and the underlying genetic effects at the species level. Here we generated and analyzed the pan-transcriptome of ~1,000 yeast natural isolates across 4,977 core and 1,468 accessory genes. We found that the accessory genome is an underappreciated driver of transcriptome divergence. Global gene expression patterns combined with population structure showed that variation in heritable expression mainly lies within subpopulation-specific signatures, for which accessory genes are overrepresented. Genome-wide association analyses consistently highlighted that accessory genes are associated with proportionally more variants with larger effect sizes, illustrating the critical role of the accessory genome on the transcriptional landscape within and between populations.

Telomere-to-telomere assemblies of 142 strains characterize the genome structural landscape in Saccharomyces cerevisiae.

Pangenomes provide access to an accurate representation of the genetic diversity of species, both in terms of sequence polymorphisms and structural variants (SVs). Here we generated the Saccharomyces cerevisiae Reference Assembly Panel (ScRAP) comprising reference-quality genomes for 142 strains representing the species' phylogenetic and ecological diversity. The ScRAP includes phased haplotype assemblies for several heterozygous diploid and polyploid isolates. We identified circa (ca.) 4,800 nonredundant SVs that provide a broad view of the genomic diversity, including the dynamics of telomere length and transposable elements. We uncovered frequent cases of complex aneuploidies where large chromosomes underwent large deletions and translocations. We found that SVs can impact gene expression near the breakpoints and substantially contribute to gene repertoire evolution. We also discovered that horizontally acquired regions insert at chromosome ends and can generate new telomeres. Overall, the ScRAP demonstrates the benefit of a pangenome in understanding genome evolution at population scale.

Loss of heterozygosity results in rapid but variable genome homogenization across yeast genetic backgrounds.

The dynamics and diversity of the appearance of genetic variants play an essential role in the evolution of the genome and the shaping of biodiversity. Recent population-wide genome sequencing surveys have highlighted the importance of loss of heterozygosity (LOH) events and have shown that they are a neglected part of the genetic diversity landscape. To assess the extent, variability, and spectrum, we explored the accumulation of LOH events in 169 heterozygous diploid Saccharomyces cerevisiae mutation accumulation lines across nine genetic backgrounds. In total, we detected a large set of 22,828 LOH events across distinct genetic backgrounds with a heterozygous level ranging from 0.1% to 1%. LOH events are very frequent with a rate consistently much higher than the mutation rate, showing their importance for genome evolution. We observed that the interstitial LOH (I-LOH) events, resulting in internal short LOH tracts, were much frequent (n = 19,660) than the terminal LOH (T-LOH) events, that is, tracts extending to the end of the chromosome (n = 3168). However, the spectrum, the rate, and the fraction of the genome under LOH vary across genetic backgrounds. Interestingly, we observed that the more the ancestors were heterozygous, the more they accumulated T-LOH events. In addition, frequent short I-LOH tracts are a signature of the lines derived from hybrids with low spore fertility. Finally, we found lines showing almost complete homozygotization during vegetative progression. Overall, our results highlight that the variable dynamics of the LOH accumulation across distinct genetic backgrounds might lead to rapid differential genome evolution during vegetative growth.

Dissection of quantitative trait loci in the Lachancea waltii yeast species highlights major hotspots.

Dissecting the genetic basis of complex trait remains a real challenge. The budding yeast Saccharomyces cerevisiae has become a model organism for studying quantitative traits, successfully increasing our knowledge in many aspects. However, the exploration of the genotype-phenotype relationship in non-model yeast species could provide a deeper insight into the genetic basis of complex traits. Here, we have studied this relationship in the Lachancea waltii species which diverged from the S. cerevisiae lineage prior to the whole-genome duplication. By performing linkage mapping analyses in this species, we identified 86 quantitative trait loci (QTL) impacting the growth in a large number of conditions. The distribution of these loci across the genome has revealed two major QTL hotspots. A first hotspot corresponds to a general growth QTL, impacting a wide range of conditions. By contrast, the second hotspot highlighted a trade-off with a disadvantageous allele for drug-free conditions which proved to be advantageous in the presence of several drugs. Finally, a comparison of the detected QTL in L. waltii with those which had been previously identified for the same trait in a closely related species, Lachancea kluyveri was performed. This analysis clearly showed the absence of shared QTL across these species. Altogether, our results represent a first step toward the exploration of the genetic architecture of quantitative trait across different yeast species.

Sex without crossing over in the yeast Saccharomycodes ludwigii.

BACKGROUND: Intermixing of genomes through meiotic reassortment and recombination of homologous chromosomes is a unifying theme of sexual reproduction in eukaryotic organisms and is considered crucial for their adaptive evolution. Previous studies of the budding yeast species Saccharomycodes ludwigii suggested that meiotic crossing over might be absent from its sexual life cycle, which is predominated by fertilization within the meiotic tetrad. RESULTS: We demonstrate that recombination is extremely suppressed during meiosis in Sd. ludwigii. DNA double-strand break formation by the conserved transesterase Spo11, processing and repair involving interhomolog interactions are required for normal meiosis but do not lead to crossing over. Although the species has retained an intact meiotic gene repertoire, genetic and population analyses suggest the exceptionally rare occurrence of meiotic crossovers in its genome. A strong AT bias of spontaneous mutations and the absence of recombination are likely responsible for its unusually low genomic GC level. CONCLUSIONS: Sd. ludwigii has followed a unique evolutionary trajectory that possibly derives fitness benefits from the combination of frequent mating between products of the same meiotic event with the extreme suppression of meiotic recombination. This life style ensures preservation of heterozygosity throughout its genome and may enable the species to adapt to its environment and survive with only minimal levels of rare meiotic recombination. We propose Sd. ludwigii as an excellent natural forum for the study of genome evolution and recombination rates.

De novo assembly and annotation of three Leptosphaeria genomes using Oxford Nanopore MinION sequencing.

Leptosphaeria maculans and Leptosphaeria biglobosa are ascomycete phytopathogens of Brassica napus (oilseed rape, canola). Here we report the complete sequence of three Leptosphaeria genomes (L. maculans JN3, L. maculans Nz-T4 and L. biglobosa G12-14). Nz-T4 and G12-14 genome assemblies were generated de novo and the reference JN3 genome assembly was improved using Oxford Nanopore MinION reads. The new assembly of L. biglobosa showed the existence of AT rich regions and pointed to a genome compartmentalization previously unsuspected following Illumina sequencing. Moreover nanopore sequencing allowed us to generate a chromosome-level assembly for the L. maculans reference isolate, JN3. The genome annotation was supported by integrating conserved proteins and RNA sequencing from Leptosphaeria-infected samples. The newly produced high-quality assemblies and annotations of those three Leptosphaeria genomes will allow further studies, notably focused on the tripartite interaction between L. maculans, L. biglobosa and oilseed rape. The discovery of as yet unknown effectors will notably allow progress in B. napus breeding towards L. maculans resistance.