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.

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.

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.

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.

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.

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.

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.

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.

Population genomics reveals chromosome-scale heterogeneous evolution in a protoploid yeast.

Yeast species represent an ideal model system for population genomic studies but large-scale polymorphism surveys have only been reported for species of the Saccharomyces genus so far. Hence, little is known about intraspecific diversity and evolution in yeast. To obtain a new insight into the evolutionary forces shaping natural populations, we sequenced the genomes of an expansive worldwide collection of isolates from a species distantly related to Saccharomyces cerevisiae: Lachancea kluyveri (formerly S. kluyveri). We identified 6.5 million single nucleotide polymorphisms and showed that a large introgression event of 1 Mb of GC-rich sequence in the chromosomal arm probably occurred in the last common ancestor of all L. kluyveri strains. Our population genomic data clearly revealed that this 1-Mb region underwent a molecular evolution pattern very different from the rest of the genome. It is characterized by a higher recombination rate, with a dramatically elevated A:T → G:C substitution rate, which is the signature of an increased GC-biased gene conversion. In addition, the predicted base composition at equilibrium demonstrates that the chromosome-scale compositional heterogeneity will persist after the genome has reached mutational equilibrium. Altogether, the data presented herein clearly show that distinct recombination and substitution regimes can coexist and lead to different evolutionary patterns within a single genome.

Ulysses: accurate detection of low-frequency structural variations in large insert-size sequencing libraries.

MOTIVATION: The detection of structural variations (SVs) in short-range Paired-End (PE) libraries remains challenging because SV breakpoints can involve large dispersed repeated sequences, or carry inherent complexity, hardly resolvable with classical PE sequencing data. In contrast, large insert-size sequencing libraries (Mate-Pair libraries) provide higher physical coverage of the genome and give access to repeat-containing regions. They can thus theoretically overcome previous limitations as they are becoming routinely accessible. Nevertheless, broad insert size distributions and high rates of chimerical sequences are usually associated to this type of libraries, which makes the accurate annotation of SV challenging. RESULTS: Here, we present Ulysses, a tool that achieves drastically higher detection accuracy than existing tools, both on simulated and real mate-pair sequencing datasets from the 1000 Human Genome project. Ulysses achieves high specificity over the complete spectrum of variants by assessing, in a principled manner, the statistical significance of each possible variant (duplications, deletions, translocations, insertions and inversions) against an explicit model for the generation of experimental noise. This statistical model proves particularly useful for the detection of low frequency variants. SV detection performed on a large insert Mate-Pair library from a breast cancer sample revealed a high level of somatic duplications in the tumor and, to a lesser extent, in the blood sample as well. Altogether, these results show that Ulysses is a valuable tool for the characterization of somatic mosaicism in human tissues and in cancer genomes.

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.

Filling annotation gaps in yeast genomes using genome-wide contact maps.

MOTIVATIONS: De novo sequencing of genomes is followed by annotation analyses aiming at identifying functional genomic features such as genes, non-coding RNAs or regulatory sequences, taking advantage of diverse datasets. These steps sometimes fail at detecting non-coding functional sequences: for example, origins of replication, centromeres and rDNA positions have proven difficult to annotate with high confidence. Here, we demonstrate an unconventional application of Chromosome Conformation Capture (3C) technique, which typically aims at deciphering the average 3D organization of genomes, by showing how functional information about the sequence can be extracted solely from the chromosome contact map. RESULTS: Specifically, we describe a combined experimental and bioinformatic procedure that determines the genomic positions of centromeres and ribosomal DNA clusters in yeasts, including species where classical computational approaches fail. For instance, we determined the centromere positions in Naumovozyma castellii, where these coordinates could not be obtained previously. Although computed centromere positions were characterized by conserved synteny with neighboring species, no consensus sequences could be found, suggesting that centromeric binding proteins or mechanisms have significantly diverged. We also used our approach to refine centromere positions in Kuraishia capsulata and to identify rDNA positions in Debaryomyces hansenii. Our study demonstrates how 3C data can be used to complete the functional annotation of eukaryotic genomes. AVAILABILITY AND IMPLEMENTATION: The source code is provided in the Supplementary Material. This includes a zipped file with the Python code and a contact matrix of Saccharomyces cerevisiae. CONTACT: romain.koszul@pasteur.fr SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A set of genetically diverged Saccharomyces cerevisiae strains with markerless deletions of multiple auxotrophic genes.

Genome analysis of over 70 Saccharomyces strains revealed the existence of five groups of genetically diverged S. cerevisiae wild-type isolates, which feature distinct genetic backgrounds and reflect the natural diversity existing among the species. The strains originated from different geographical and ecological niches (Malaysian, West African, North American, Wine/European and Sake) and represent clean, non-mosaic lineages of S. cerevisiae, meaning that their genomes differ essentially by monomorphic and private SNPs. In this study, one representative strain for each of the five S. cerevisiae clean lineages was selected and mutated for several auxotroph genes by clean markerless deletions, so that all dominant markers remained available for further genetic manipulations. A set of 50 strains was assembled, including eight haploid and two diploid strains for each lineage. These strains carry different combinations of leu2∆0, lys2∆0, met15∆0, ura3∆0 and/or ura3∆::KanMX-barcoded deletions with marker configurations resembling that of the BY series, which will allow large-scale crossing with existing deletion collections. This new set of genetically tractable strains provides a powerful tool kit to explore the impact of natural variation on complex biological processes.

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.

The mutational profile of the yeast genome is shaped by replication.

Despite the scrutiny that has been directed for years at the yeast genome, relatively little is known about the impact of replication on the substitution dynamics in Saccharomyces cerevisiae. Here, we show that the mutation rate increases with the replication timing by more than 30% between the earliest and the latest replicating regions. In addition, we found a mutational asymmetry associated with the polarity of replication resulting in higher rates of substitutions toward C and A than toward G and T in leading strands (reciprocally more substitutions toward G and T in lagging strands). Such mutational asymmetries applied over long evolutionary periods should generate compositional skews between the two DNA strands. Thus, we show that the leading replicating strands present an excess of C over G and of A over T in the genome of S. cerevisiae (reciprocally an excess of G + T over C + A in lagging strands). We also show that the nucleotide frequencies at mutational equilibrium predict a compositional skew at equilibrium very close to the observed skew between leading and lagging strands, suggesting that compositional equilibrium has been nearly attained in the present day genome of S. cerevisiae. Surprisingly, the direction of this skew is inverted compared with the one in the human genome.

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.

Unusual composition of a yeast chromosome arm is associated with its delayed replication.

The 11.3-Mb genome of the yeast Lachancea (Saccharomyces) kluyveri displays an intriguing compositional heterogeneity: a region of approximately 1 Mb, covering almost the whole left arm of chromosome C (C-left), has an average GC content of 52.9%, which is significantly higher than the 40.4% global GC content of the rest of the genome. This region contains the MAT locus, which remains normal in composition. The excess of GC base pairs affects both coding and noncoding sequences, and thus is not due to selective pressure acting on protein sequences. It leads to a strong codon usage bias and alters the amino acid composition of the 457 proteins encoded on C-left that do not show obvious bias for functional categories, or the presence of paralogs or orthologs of essential genes of Saccharomyces cerevisiae. They share significant synteny conservation with other species of the Saccharomycetaceae, and phylogenetic analysis indicates that C-left originates from a Lachancea species. In contrast, there is a complete absence of transposable elements in C-left, whereas 18 elements per megabase are distributed across the rest of the genome. Comparative hybridization of synchronized cells using high-density genome arrays reveals that C-left is replicated later during S phase than the rest of the genome. Two possible primary causes of this major compositional heterogeneity are discussed: an ancient hybridization of two related species with very distinct GC composition, or an intrinsic mechanism, possibly associated with the loss of the silent cassettes from C-left that progressively increased the GC content and generated the delayed replication of this chromosomal arm.