Domestication signatures in the non-conventional yeast Lachancea cidri.

Evaluating domestication signatures beyond model organisms is essential for a thorough understanding of the genotype-phenotype relationship in wild and human-related environments. Structural variations (SVs) can significantly impact phenotypes playing an important role in the physiological adaptation of species to different niches, including during domestication. A detailed characterization of the fitness consequences of these genomic rearrangements, however, is still limited in non-model systems, largely due to the paucity of direct comparisons between domesticated and wild isolates. Here, we used a combination of sequencing strategies to explore major genomic rearrangements in a Lachancea cidri yeast strain isolated from cider (CBS2950) and compared them to those in eight wild isolates from primary forests. Genomic analysis revealed dozens of SVs, including a large reciprocal translocation (~16 kb and 500 kb) present in the cider strain, but absent from all wild strains. Interestingly, the number of SVs was higher relative to single-nucleotide polymorphisms in the cider strain, suggesting a significant role in the strain's phenotypic variation. The set of SVs identified directly impacts dozens of genes and likely underpins the greater fermentation performance in the L. cidri CBS2950. In addition, the large reciprocal translocation affects a proline permease (PUT4) regulatory region, resulting in higher PUT4 transcript levels, which agrees with higher ethanol tolerance, improved cell growth when using proline, and higher amino acid consumption during fermentation. These results suggest that SVs are responsible for the rapid physiological adaptation of yeast to a human-related environment and demonstrate the key contribution of SVs in adaptive fermentative traits in non-model species.IMPORTANCEThe exploration of domestication signatures associated with human-related environments has predominantly focused on studies conducted on model organisms, such as Saccharomyces cerevisiae, overlooking the potential for comparisons across other non-Saccharomyces species. In our research, employing a combination of long- and short-read data, we found domestication signatures in Lachancea cidri, a non-model species recently isolated from fermentative environments in cider in France. The significance of our study lies in the identification of large array of major genomic rearrangements in a cider strain compared to wild isolates, which underly several fermentative traits. These domestication signatures result from structural variants, which are likely responsible for the phenotypic differences between strains, providing a rapid path of adaptation to human-related environments.

Telomerase-independent survival leads to a mosaic of complex subtelomere rearrangements in Chlamydomonas reinhardtii.

Telomeres and subtelomeres, the genomic regions located at chromosome extremities, are essential for genome stability in eukaryotes. In the absence of the canonical maintenance mechanism provided by telomerase, telomere shortening induces genome instability. The landscape of the ensuing genome rearrangements is not accessible by short-read sequencing. Here, we leverage Oxford Nanopore Technologies long-read sequencing to survey the extensive repertoire of genome rearrangements in telomerase mutants of the model green microalga Chlamydomonas reinhardtii In telomerase-mutant strains grown for hundreds of generations, most chromosome extremities were capped by short telomere sequences that were either recruited de novo from other loci or maintained in a telomerase-independent manner. Other extremities did not end with telomeres but only with repeated subtelomeric sequences. The subtelomeric elements, including rDNA, were massively rearranged and involved in breakage-fusion-bridge cycles, translocations, recombinations, and chromosome circularization. These events were established progressively over time and displayed heterogeneity at the subpopulation level. New telomere-capped extremities composed of sequences originating from more internal genomic regions were associated with high DNA methylation, suggesting that de novo heterochromatin formation contributes to the restoration of chromosome end stability in C. reinhardtii The diversity of alternative strategies present in the same organism to maintain chromosome integrity and the variety of rearrangements found in telomerase mutants are remarkable, and illustrate genome plasticity at short timescales.

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.

Late Pleistocene-dated divergence between South Hemisphere populations of the non-conventional yeast L. cidri.

Most organisms belonging to the Saccharomycotina subphylum have high genetic diversity and a vast repertoire of metabolisms and lifestyles. Lachancea cidri is an ideal yeast model for exploring the interplay between genetics, ecological function and evolution. Lachancea cidri diverged from the Saccharomyces lineage before the whole-genome duplication and is distributed across the South Hemisphere, displaying an important ecological success. We applied phylogenomics to investigate the genetic variation of L. cidri isolates obtained from Australia and South America. Our approach revealed the presence of two main lineages according to their geographic distribution (Aus and SoAm). Estimation of the divergence time suggests that SoAm and Aus lineages diverged near the last glacial maximum event during the Pleistocene (64-8 KYA). Interestingly, we found that the French reference strain is closely related to the Australian strains, with a recent divergence (405-51 YA), likely associated to human movements. Additionally, we identified different lineages within the South American population, revealing that Patagonia contains a similar genetic diversity comparable to that of other lineages in S. cerevisiae. These findings support the idea of a Pleistocene-dated divergence between South Hemisphere lineages, where the Nothofagus and Araucaria ecological niches likely favoured the extensive distribution of L. cidri in Patagonia.

An evolutionary model identifies the main evolutionary biases for the evolution of genome-replication profiles.

Recent results comparing the temporal program of genome replication of yeast species belonging to the Lachancea clade support the scenario that the evolution of the replication timing program could be mainly driven by correlated acquisition and loss events of active replication origins. Using these results as a benchmark, we develop an evolutionary model defined as birth-death process for replication origins and use it to identify the evolutionary biases that shape the replication timing profiles. Comparing different evolutionary models with data, we find that replication origin birth and death events are mainly driven by two evolutionary pressures, the first imposes that events leading to higher double-stall probability of replication forks are penalized, while the second makes less efficient origins more prone to evolutionary loss. This analysis provides an empirically grounded predictive framework for quantitative evolutionary studies of the replication timing program.

A Versatile Protocol to Generate Translocations in Yeast Genomes Using CRISPR/Cas9.

Genomic engineering methods represent powerful tools to examine chromosomal modifications and to subsequently study their impacts on cellular phenotypes. However, quantifying the fitness impact of translocations, independently from base substitutions or the insertion of genetic markers, remains a challenge. Here we report a rapid and straightforward protocol for engineering either targeted reciprocal translocations at the base pair level of resolution between two chromosomes or multiple simultaneous rearrangements in the yeast genome, without inserting any marker sequence in the chromosomes. Our CRISPR/Cas9-based method consists of inducing either (1) two double-strand breaks (DSBs) in two different chromosomes with two distinct guide RNAs (gRNAs) while providing specifically designed homologous donor DNA forcing the trans-repair of chromosomal extremities to generate a targeted reciprocal translocation or (2) multiple DSBs with a single gRNA targeting dispersed repeated sequences and leaving endogenous uncut copies of the repeat to be used as donor DNA, thereby generating multiple translocations, often associated with large segmental duplications (Fleiss, et al. PLoS Genet 15:e1008332, 2019).

The budding yeast life cycle: More complex than anticipated?

The budding yeast, Saccharomyces cerevisiae, has served as a model for nearly a century to understand the principles of the eukaryotic life cycle. The canonical life cycle of S. cerevisiae comprises a regular alternation between haploid and diploid phases. Haploid gametes generated by sporulation are expected to quickly restore the diploid phase mainly through inbreeding via intratetrad mating or haploselfing, thereby promoting genome homozygotization. However, recent large population genomics data unveiled that heterozygosity and polyploidy are unexpectedly common. This raises the interesting paradox of a haplo-diplobiontic species being well-adapted to inbreeding and able to maintain high levels of heterozygosity and polyploidy, thereby suggesting an unanticipated complexity of the yeast life cycle. Here, we propose that unprogrammed mating type switching, heterothallism, reduced spore formation and viability, cell-cell fusion and dioecy could play key and uncharted contributions to generate and maintain heterozygosity through polyploidization.

Highly Contiguous Nanopore Genome Assembly of Chlamydomonas reinhardtii CC-1690.

The current Chlamydomonas reinhardtii reference genome remains fragmented due to gaps stemming from large repetitive regions. To overcome the vast majority of these gaps, publicly available Oxford Nanopore Technology data were used to create a new reference-quality de novo genome assembly containing only 21 contigs, 30/34 telomeric ends, and a genome size of 111 Mb.

Phylogenetic Reconstruction Based on Synteny Block and Gene Adjacencies.

Gene order can be used as an informative character to reconstruct phylogenetic relationships between species independently from the local information present in gene/protein sequences. PhyChro is a reconstruction method based on chromosomal rearrangements, applicable to a wide range of eukaryotic genomes with different gene contents and levels of synteny conservation. For each synteny breakpoint issued from pairwise genome comparisons, the algorithm defines two disjoint sets of genomes, named partial splits, respectively, supporting the two block adjacencies defining the breakpoint. Considering all partial splits issued from all pairwise comparisons, a distance between two genomes is computed from the number of partial splits separating them. Tree reconstruction is achieved through a bottom-up approach by iteratively grouping sister genomes minimizing genome distances. PhyChro estimates branch lengths based on the number of synteny breakpoints and provides confidence scores for the branches. PhyChro performance is evaluated on two data sets of 13 vertebrates and 21 yeast genomes by using up to 130,000 and 179,000 breakpoints, respectively, a scale of genomic markers that has been out of reach until now. PhyChro reconstructs very accurate tree topologies even at known problematic branching positions. Its robustness has been benchmarked for different synteny block reconstruction methods. On simulated data PhyChro reconstructs phylogenies perfectly in almost all cases, and shows the highest accuracy compared with other existing tools. PhyChro is very fast, reconstructing the vertebrate and yeast phylogenies in <15 min.

MUM&Co: accurate detection of all SV types through whole-genome alignment.

SUMMARY: MUM&Co is a single bash script to detect structural variations (SVs) utilizing whole-genome alignment (WGA). Using MUMmer's nucmer alignment, MUM&Co can detect insertions, deletions, tandem duplications, inversions and translocations greater than 50 bp. Its versatility depends upon the WGA and therefore benefits from contiguous de-novo assemblies generated by third generation sequencing technologies. Benchmarked against five WGA SV-calling tools, MUM&Co outperforms all tools on simulated SVs in yeast, plant and human genomes and performs similarly in two real human datasets. Additionally, MUM&Co is particularly unique in its ability to find inversions in both simulated and real datasets. Lastly, MUM&Co's primary output is an intuitive tabulated file containing a list of SVs with only necessary genomic details. AVAILABILITY AND IMPLEMENTATION: https://github.com/SAMtoBAM/MUMandCo. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Molecular profiling of beer wort fermentation diversity across natural Saccharomyces eubayanus isolates.

The utilization of S. eubayanus has recently become a topic of interest due to the novel organoleptic properties imparted to beer. However, the utilization of S. eubayanus in brewing requires the comprehension of the mechanisms that underlie fermentative differences generated from its natural genetic variability. Here, we evaluated fermentation performance and volatile compound production in ten genetically distinct S. eubayanus strains in a brewing fermentative context. The evaluated strains showed a broad phenotypic spectrum, some of them exhibiting a high fermentation capacity and high levels of volatile esters and/or higher alcohols. Subsequently, we obtained molecular profiles by generating 'end-to-end' genome assemblies, as well as metabolome and transcriptome profiling of two Patagonian isolates exhibiting significant differences in beer aroma profiles. These strains showed clear differences in concentrations of intracellular metabolites, including amino acids, such as valine, leucine and isoleucine, likely impacting the production of 2-methylpropanol and 3-methylbutanol. These differences in the production of volatile compounds are attributed to gene expression variation, where the most profound differentiation is attributed to genes involved in assimilatory sulfate reduction, which in turn validates phenotypic differences in H2 S production. This study lays a solid foundation for future research to improve fermentation performance and select strains for new lager styles based on aroma and metabolic profiles.

Reshuffling yeast chromosomes with CRISPR/Cas9.

Genome engineering is a powerful approach to study how chromosomal architecture impacts phenotypes. However, quantifying the fitness impact of translocations independently from the confounding effect of base substitutions has so far remained challenging. We report a novel application of the CRISPR/Cas9 technology allowing to generate with high efficiency both uniquely targeted and multiple concomitant reciprocal translocations in the yeast genome. Targeted translocations are constructed by inducing two double-strand breaks on different chromosomes and forcing the trans-chromosomal repair through homologous recombination by chimerical donor DNAs. Multiple translocations are generated from the induction of several DSBs in LTR repeated sequences and promoting repair using endogenous uncut LTR copies as template. All engineered translocations are markerless and scarless. Targeted translocations are produced at base pair resolution and can be sequentially generated one after the other. Multiple translocations result in a large diversity of karyotypes and are associated in many instances with the formation of unanticipated segmental duplications. To test the phenotypic impact of translocations, we first recapitulated in a lab strain the SSU1/ECM34 translocation providing increased sulphite resistance to wine isolates. Surprisingly, the same translocation in a laboratory strain resulted in decreased sulphite resistance. However, adding the repeated sequences that are present in the SSU1 promoter of the resistant wine strain induced sulphite resistance in the lab strain, yet to a lower level than that of the wine isolate, implying that additional polymorphisms also contribute to the phenotype. These findings illustrate the advantage brought by our technique to untangle the phenotypic impacts of structural variations from confounding effects of base substitutions. Secondly, we showed that strains with multiple translocations, even those devoid of unanticipated segmental duplications, display large phenotypic diversity in a wide range of environmental conditions, showing that simply reconfiguring chromosome architecture is sufficient to provide fitness advantages in stressful growth conditions.

André Goffeau's imprinting on second generation yeast "genomologists".

All authors of the present paper have worked in labs that participated to the sequencing effort of the Saccharomyces cerevisiae reference genome, and we owe to this the fact that we have all chosen to work on genomics of yeasts. S. cerevisiae has been a popular model species for genetics since the 20th century as well as being a model for general eukaryotic cellular processes. Although it has also been used empirically in fermentation for millennia, there was until recently, a lack of knowledge about the natural and evolutionary history of this yeast. The achievement of the international effort to sequence its genome was the foundation for understanding many eukaryotic biological processes but also represented the first step towards the study of the genome and ecological diversity of yeast populations worldwide. We will describe recent advances in yeast comparative and population genomics that find their origins in the S. cerevisiae genome project initiated and pursued by André Goffeau.

Characterizing meiotic chromosomes' structure and pairing using a designer sequence optimized for Hi-C.

In chromosome conformation capture experiments (Hi-C), the accuracy with which contacts are detected varies due to the uneven distribution of restriction sites along genomes. In addition, repeated sequences or homologous regions remain indistinguishable because of the ambiguities they introduce during the alignment of the sequencing reads. We addressed both limitations by designing and engineering 144 kb of a yeast chromosome with regularly spaced restriction sites (Syn-HiC design). In the Syn-HiC region, Hi-C signal-to-noise ratio is enhanced and can be used to measure the shape of an unbiased distribution of contact frequencies, allowing to propose a robust definition of a Hi-C experiment resolution. The redesigned region is also distinguishable from its native homologous counterpart in an otherwise isogenic diploid strain. As a proof of principle, we tracked homologous chromosomes during meiotic prophase in synchronized and pachytene-arrested cells and captured important features of their spatial reorganization, such as chromatin restructuration into arrays of Rec8-delimited loops, centromere declustering, individualization, and pairing. Overall, we illustrate the promises held by redesigning genomic regions to explore complex biological questions.

The evolution of the temporal program of genome replication.

Genome replication is highly regulated in time and space, but the rules governing the remodeling of these programs during evolution remain largely unknown. We generated genome-wide replication timing profiles for ten Lachancea yeasts, covering a continuous evolutionary range from closely related to more divergent species. We show that replication programs primarily evolve through a highly dynamic evolutionary renewal of the cohort of active replication origins. We found that gained origins appear with low activity yet become more efficient and fire earlier as they evolutionarily age. By contrast, origins that are lost comprise the complete range of firing strength. Additionally, they preferentially occur in close vicinity to strong origins. Interestingly, despite high evolutionary turnover, active replication origins remain regularly spaced along chromosomes in all species, suggesting that origin distribution is optimized to limit large inter-origin intervals. We propose a model on the evolutionary birth, death, and conservation of active replication origins.

A Molecular Portrait of De Novo Genes in Yeasts.

New genes, with novel protein functions, can evolve "from scratch" out of intergenic sequences. These de novo genes can integrate the cell's genetic network and drive important phenotypic innovations. Therefore, identifying de novo genes and understanding how the transition from noncoding to coding occurs are key problems in evolutionary biology. However, identifying de novo genes is a difficult task, hampered by the presence of remote homologs, fast evolving sequences and erroneously annotated protein coding genes. To overcome these limitations, we developed a procedure that handles the usual pitfalls in de novo gene identification and predicted the emergence of 703 de novo gene candidates in 15 yeast species from 2 genera whose phylogeny spans at least 100 million years of evolution. We validated 85 candidates by proteomic data, providing new translation evidence for 25 of them through mass spectrometry experiments. We also unambiguously identified the mutations that enabled the transition from noncoding to coding for 30 Saccharomyces de novo genes. We established that de novo gene origination is a widespread phenomenon in yeasts, only a few being ultimately maintained by selection. We also found that de novo genes preferentially emerge next to divergent promoters in GC-rich intergenic regions where the probability of finding a fortuitous and transcribed ORF is the highest. Finally, we found a more than 3-fold enrichment of de novo genes at recombination hot spots, which are GC-rich and nucleosome-free regions, suggesting that meiotic recombination contributes to de novo gene emergence in yeasts.

Contrasting evolutionary genome dynamics between domesticated and wild yeasts.

Structural rearrangements have long been recognized as an important source of genetic variation, with implications in phenotypic diversity and disease, yet their detailed evolutionary dynamics remain elusive. Here we use long-read sequencing to generate end-to-end genome assemblies for 12 strains representing major subpopulations of the partially domesticated yeast Saccharomyces cerevisiae and its wild relative Saccharomyces paradoxus. These population-level high-quality genomes with comprehensive annotation enable precise definition of chromosomal boundaries between cores and subtelomeres and a high-resolution view of evolutionary genome dynamics. In chromosomal cores, S. paradoxus shows faster accumulation of balanced rearrangements (inversions, reciprocal translocations and transpositions), whereas S. cerevisiae accumulates unbalanced rearrangements (novel insertions, deletions and duplications) more rapidly. In subtelomeres, both species show extensive interchromosomal reshuffling, with a higher tempo in S. cerevisiae. Such striking contrasts between wild and domesticated yeasts are likely to reflect the influence of human activities on structural genome evolution.

Dissection of quantitative traits by bulk segregant mapping in a protoploid yeast species.

Since more than a decade ago, Saccharomyces cerevisiae has been used as a model to dissect complex traits, revealing the genetic basis of a large number of traits in fine detail. However, to have a more global view of the genetic architecture of traits across species, the examination of the molecular basis of phenotypes within non-conventional species would undoubtedly be valuable. In this respect, the Saccharomycotina yeasts represent ideal and potential non-model organisms. Here we sought to assess the feasibility of genetic mapping by bulk segregant analysis in the protoploid Lachancea kluyveri (formerly S. kluyveri) yeast species, a distantly related species to S. cerevisiae For this purpose, we designed a fluorescent mating-type marker, compatible with any mating-competent strains representative of this species, to rapidly create a large population of haploid segregants (>10(5) cells). Quantitative trait loci can be mapped by selecting and sequencing an enriched pool of progeny with extreme phenotypic values. As a test bed, we applied this strategy and mapped the causal loci underlying halotolerance phenotypes in L. kluyveri Overall, this study demonstrates that bulk segregant mapping is a powerful way for investigating the genetic basis of natural variations in non-model yeast organisms and more precisely in L. kluyveri.

Reconstruction of ancestral chromosome architecture and gene repertoire reveals principles of genome evolution in a model yeast genus.

Reconstructing genome history is complex but necessary to reveal quantitative principles governing genome evolution. Such reconstruction requires recapitulating into a single evolutionary framework the evolution of genome architecture and gene repertoire. Here, we reconstructed the genome history of the genus Lachancea that appeared to cover a continuous evolutionary range from closely related to more diverged yeast species. Our approach integrated the generation of a high-quality genome data set; the development of AnChro, a new algorithm for reconstructing ancestral genome architecture; and a comprehensive analysis of gene repertoire evolution. We found that the ancestral genome of the genus Lachancea contained eight chromosomes and about 5173 protein-coding genes. Moreover, we characterized 24 horizontal gene transfers and 159 putative gene creation events that punctuated species diversification. We retraced all chromosomal rearrangements, including gene losses, gene duplications, chromosomal inversions and translocations at single gene resolution. Gene duplications outnumbered losses and balanced rearrangements with 1503, 929, and 423 events, respectively. Gene content variations between extant species are mainly driven by differential gene losses, while gene duplications remained globally constant in all lineages. Remarkably, we discovered that balanced chromosomal rearrangements could be responsible for up to 14% of all gene losses by disrupting genes at their breakpoints. Finally, we found that nonsynonymous substitutions reached fixation at a coordinated pace with chromosomal inversions, translocations, and duplications, but not deletions. Overall, we provide a granular view of genome evolution within an entire eukaryotic genus, linking gene content, chromosome rearrangements, and protein divergence into a single evolutionary framework.