Population genomics of the fission yeast Schizosaccharomyces pombe

Pheromone Response and Mating Behavior in Fission Yeast

Most ascomycete fungi, including the fission yeast Schizosaccharomyces pombe, secrete two peptidyl mating pheromones: C-terminally modified and unmodified peptides. S. pombe has two mating types, plus and minus, which secrete two different pheromones, P-factor (unmodified) and M-factor (modified), respectively. These pheromones are specifically recognized by receptors on the cell surface of cells of opposite mating types, which trigger a pheromone response. Recognition between pheromones and their corresponding receptors is important for mate discrimination; therefore, genetic changes in pheromone or receptor genes affect mate recognition and cause reproductive isolation that limits gene flow between populations. Such genetic variation in recognition via the pheromone/receptor system may drive speciation. Our recent studies reported that two pheromone receptors in S. pombe might have different stringencies in pheromone recognition. In this review, we focus on the molecular mechanism of pheromone response and mating behavior, emphasizing pheromone diversification and its impact on reproductive isolation in S. pombe and closely related fission yeast species. We speculate that the "asymmetric" system might allow flexible adaptation to pheromone mutational changes while maintaining stringent recognition of mating partners. The loss of pheromone activity results in the extinction of an organism's lineage. Therefore, genetic changes in pheromones and their receptors may occur gradually and/or coincidently before speciation. Our findings suggest that the M-factor plays an important role in partner discrimination, whereas P-factor communication allows flexible adaptation to create variations in S. pombe. Our inferences provide new insights into the evolutionary mechanisms underlying pheromone diversification.

Enforcement of Postzygotic Species Boundaries in the Fungal Kingdom

Understanding the molecular basis of speciation is a primary goal in evolutionary biology. The formation of the postzygotic reproductive isolation that causes hybrid dysfunction, thereby reducing gene flow between diverging populations, is crucial for speciation. Using various advanced approaches, including chromosome replacement, hybrid introgression and transcriptomics, population genomics, and experimental evolution, scientists have revealed multiple mechanisms involved in postzygotic barriers in the fungal kingdom. These results illuminate both unique and general features of fungal speciation. Our review summarizes experiments on fungi exploring how Dobzhansky-Muller incompatibility, killer meiotic drive, chromosome rearrangements, and antirecombination contribute to postzygotic reproductive isolation. We also discuss possible evolutionary forces underlying different reproductive isolation mechanisms and the potential roles of the evolutionary arms race under the Red Queen hypothesis and epigenetic divergence in speciation.

Segregating Complete Tf2 Elements Are Largely Neutral in Fission Yeast

Transposable elements (TEs) comprise a large proportion of the eukaryote genomes. Yet it remains poorly understood how TEs influence the fitness of the hosts carrying them. Here, we empirically test the impact of TEs on the host fitness in the fission yeast Schizosaccharomyces pombe. We find that two families of TEs (Tf1 and Tf2 elements), both of which belong to long terminal repeat retrotransposons, are highly polymorphic among individual S. pombe strains. Only 13 complete Tf2 elements are identified in S. pombe laboratory strain 972. These 13 Tf2 elements integrated into host genomes in very recent time and are segregating within the S. pombe population. Through knocking out each of the 13 Tf2 elements in S. pombe strain 972, we find Tf2 knockout does not affect the host fitness, and Tf2 elements do not alter the expression of nearby genes. Challenged by diverse forms of stress, the Tf2 knockout strains do not exhibit different growth rates from wild-type strain. Together, we conclude that segregating complete Tf2 elements insertions are largely neutral to host fitness in the fission yeast. Our study provides genome-wide empirical support for the selfish nature of TEs in fission yeast.

Genome Comparisons of the Fission Yeasts Reveal Ancient Collinear Loci Maintained by Natural Selection

Fission yeasts have a unique life history and exhibit distinct evolutionary patterns from other yeasts. Besides, the species demonstrate stable genome structures despite the relatively fast evolution of their genomic sequences. To reveal what could be the reason for that, comparative genomic analyses were carried out. Our results provided evidence that the structural and sequence evolution of the fission yeasts were correlated. Moreover, we revealed ancestral locally collinear blocks (aLCBs), which could have been inherited from their last common ancestor. These aLCBs proved to be the most conserved regions of the genomes as the aLCBs contain almost eight genes/blocks on average in the same orientation and order across the species. Gene order of the aLCBs is mainly fission-yeast-specific but supports the idea of filamentous ancestors. Nevertheless, the sequences and gene structures within the aLCBs are as mutable as any sequences in other parts of the genomes. Although genes of certain Gene Ontology (GO) categories tend to cluster at the aLCBs, those GO enrichments are not related to biological functions or high co-expression rates, they are, rather, determined by the density of essential genes and Rec12 cleavage sites. These data and our simulations indicated that aLCBs might not only be remnants of ancestral gene order but are also maintained by natural selection.

Variables Influencing Differences in Sequence Conservation in the Fission Yeast Schizosaccharomyces pombe

Which variables determine the constraints on gene sequence evolution is one of the most central questions in molecular evolution. In the fission yeast Schizosaccharomyces pombe, an important model organism, the variables influencing the rate of sequence evolution have yet to be determined. Previous studies in other single celled organisms have generally found gene expression levels to be most significant, with numerous other variables such as gene length and functional importance identified as having a smaller impact. Using publicly available data, we used partial least squares regression, principal components regression, and partial correlations to determine the variables most strongly associated with sequence evolution constraints. We identify centrality in the protein–protein interactions network, amino acid composition, and cellular location as the most important determinants of sequence conservation. However, each factor only explains a small amount of variance, and there are numerous variables having a significant or heterogeneous influence. Our models explain more than half of the variance in dN, raising the possibility that future refined models could quantify the role of stochastics in evolutionary rate variation.

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.

High Complexity and Degree of Genetic Variation in Brettanomyces bruxellensis Population

Genome-wide characterization of genetic variants of a large population of individuals within the same species is essential to have a deeper insight into its evolutionary history as well as the genotype–phenotype relationship. Population genomic surveys have been performed in multiple yeast species, including the two model organisms, Saccharomyces cerevisiae and Schizosaccharomyces pombe. In this context, we sought to characterize at the population level the Brettanomyces bruxellensis yeast species, which is a major cause of wine spoilage and can contribute to the specific flavor profile of some Belgium beers. We have completely sequenced the genome of 53 B. bruxellensis strains isolated worldwide. The annotation of the reference genome allowed us to define the gene content of this species. As previously suggested, our genomic data clearly highlighted that genetic diversity variation is related to ploidy level, which is variable in the B. bruxellensis species. Genomes are punctuated by multiple loss-of-heterozygosity regions, whereas aneuploidies as well as segmental duplications are uncommon. Interestingly, triploid genomes are more prone to gene copy number variation than diploids. Finally, the pangenome of the species was reconstructed and was found to be small with few accessory genes compared with S. cerevisiae. The pangenome is composed of 5,409 ORFs (open reading frames) among which 5,106 core ORFs and 303 ORFs that are variable within the population. All these results highlight the different trajectories of species evolution and consequently the interest of establishing population genomic surveys in more populations.

The rate and molecular spectrum of mutation are selectively maintained in yeast

What determines the rate (μ) and molecular spectrum of mutation is a fundamental question. The prevailing hypothesis asserts that natural selection against deleterious mutations has pushed μ to the minimum achievable in the presence of genetic drift, or the drift barrier. Here we show that, contrasting this hypothesis, μ substantially exceeds the drift barrier in diverse organisms. Random mutation accumulation (MA) in yeast frequently reduces μ, and deleting the newly discovered mutator gene PSP2 nearly halves μ. These results, along with a comparison between the MA and natural yeast strains, demonstrate that μ is maintained above the drift barrier by stabilizing selection. Similar comparisons show that the mutation spectrum such as the universal AT mutational bias is not intrinsic but has been selectively preserved. These findings blur the separation of mutation from selection as distinct evolutionary forces but open the door to alleviating mutagenesis in various organisms by genome editing.

How natural selection shapes the rate and molecular spectrum of mutations is debated. Yeast mutation accumulation experiments identify a gene promoting mutagenesis and show stabilizing selection maintaining the mutation rate above the drift barrier. Selection also preserves the mutation spectrum.

Ascorbic acid supplementation suppresses cadmium-derived alterations in the fission yeast Schizosaccharomyces pombe

Cadmium (Cd) a highly toxic environmental pollutant, that does not have any physiological function in the organism, represents a great concern for human health as it can be easily transported from its environmental sources to the food chain. Food, water, and air are the major sources of Cd exposure to the population. Cd-mediated impairments of the basic cellular properties largely depend on its ability to enhance the formation of reactive oxygen species (ROS) and thus triggers oxidative stress to the cell. With the use of fission yeast Schizosaccharomyces pombe (S. pombe) as the model organism, we have analyzed the impact of Cd on the cell growth intensity, as it represents the fundamental feature of all living organisms. Cells were incubated with different Cd concentrations for 3, 6, and 9 hours to investigate the effect of Cd on cell growth in a time and dose-dependent manner. Further possible Cd-derived alterations, as the peroxidation of membrane lipids or the functional impairment of the enzymatic antioxidant protection mechanisms, were investigated by determination of the MDA content and via catalase (CAT) activity detection. Moreover, ascorbic acid (AsA) pre-treatment was subjected to investigate the assumed positive effect of AsA against Cd toxicity. We show here on one hand that cells suffer under the influence of Cd, but on the other hand, they substantially profit from AsA supplementation. Because S. pombe is known to shares many molecular, and biochemical similarities with higher organisms, the effect of AsA in cadmium toxicity elimination might be expected to a similar extent also in other cell types.

Large-scale sequencing and comparative analysis of oenological Saccharomyces cerevisiae strains supported by nanopore refinement of key genomes

Saccharomyces cerevisiae has long been part of human activities related to the production of food and wine. The industrial demand for fermented beverages with well-defined and stable characteristics boosted the isolation and selection of strains conferring a distinctive aroma profile to the final product. To uncover variants characterizing oenological strains, the sequencing of 65 new S. cerevisiae isolates, and the comparison with other 503 publicly available genomes were performed. A hybrid approach based on short Illumina and long Oxford Nanopore reads allowed the in-depth investigation of eleven genomes and the identification of putative laterally transferred regions and structural variants. A comparative analysis between clusters of strains belonging to different datasets allowed the identification of novel relevant genetic features including single nucleotide polymorphisms, insertions and structural variants. Detection of oenological single nucleotide variants shed light on the existence of different levels of modulation for the mevalonate pathway relevant for the biosynthesis of aromatic compounds.

A model of autowave self-organization as a hierarchy of active media in the biological evolution

Within the framework of the active media concept, we develop a biophysical model of autowave self-organization which is treated as a hierarchy of active media in the evolution of the biosphere. We also propose a mathematical model of the autowave process of speciation in a flow of mutations for the three main taxonometric groups (prokaryotes, unicellular and multicellular eukaryotes) with a naturally determined lower boundary of living matter (the appearance of prokaryotes) and an open upper boundary for the formation of new species. It is shown that the fluctuation-bifurcation description of the evolution for the formation of new taxonometric groups as a trajectory of transformation of small fluctuations into giant ones adequately reflects the process of self-organization during the formation of taxa. The major concepts of biological evolution, conditions of hierarchy formation as a fundamental manifestation of self-organization and complexity in the evolution of biological systems are considered.

Low Base-Substitution Mutation Rate but High Rate of Slippage Mutations in the Sequence Repeat-Rich Genome of Dictyostelium discoideum

We describe the rate and spectrum of spontaneous mutations for the social amoeba Dictyostelium discoideum, a key model organism in molecular, cellular, evolutionary and developmental biology. Whole-genome sequencing of 37 mutation accumulation lines of D. discoideum after an average of 1,500 cell divisions yields a base-substitution mutation rate of 2.47 × 10⁻¹¹ per site per generation, substantially lower than that of most eukaryotic and prokaryotic organisms, and of the same order of magnitude as in the ciliates Paramecium tetraurelia and Tetrahymena thermophila. Known for its high genomic AT content and abundance of simple sequence repeats, we observe that base-substitution mutations in D. discoideum are highly A/T biased. This bias likely contributes both to the high genomic AT content and to the formation of simple sequence repeats in the AT-rich genome of Dictyostelium discoideum. In contrast to the situation in other surveyed unicellular eukaryotes, indel rates far exceed the base-substitution mutation rate in this organism with a high proportion of 3n indels, particularly in regions without simple sequence repeats. Like ciliates, D. discoideum has a large effective population size, reducing the power of random genetic drift, magnifying the effect of selection on replication fidelity, in principle allowing D. discoideum to evolve an extremely low base-substitution mutation rate.

Atypical meiosis can be adaptive in outcrossed Schizosaccharomyces pombe due to wtf meiotic drivers

Killer meiotic drivers are genetic parasites that destroy ‘sibling’ gametes lacking the driver allele. The fitness costs of drive can lead to selection of unlinked suppressors. This suppression could involve evolutionary tradeoffs that compromise gametogenesis and contribute to infertility. Schizosaccharomyces pombe, an organism containing numerous gamete (spore)-killing wtf drivers, offers a tractable system to test this hypothesis. Here, we demonstrate that in scenarios analogous to outcrossing, wtf drivers generate a fitness landscape in which atypical spores, such as aneuploids and diploids, are advantageous. In this context, wtf drivers can decrease the fitness costs of mutations that disrupt meiotic fidelity and, in some circumstances, can even make such mutations beneficial. Moreover, we find that S. pombe isolates vary greatly in their ability to make haploid spores, with some isolates generating up to 46% aneuploid or diploid spores. This work empirically demonstrates the potential for meiotic drivers to shape the evolution of gametogenesis.

Adaptation to Industrial Stressors Through Genomic and Transcriptional Plasticity in a Bioethanol Producing Fission Yeast Isolate

Schizosaccharomyces pombe is a model unicellular eukaryote with ties to the basic research, oenology and industrial biotechnology sectors. While most investigations into S. pombe cell biology utilize Leupold’s 972h^- laboratory strain background, recent studies have described a wealth of genetic and phenotypic diversity within wild populations of S. pombe including stress resistance phenotypes which may be of interest to industry. Here we describe the genomic and transcriptomic characterization of Wilmar-P, an S. pombe isolate used for bioethanol production from sugarcane molasses at industrial scale. Novel sequences present in Wilmar-P but not in the laboratory S. pombe genome included multiple coding sequences with near-perfect nucleotide identity to Schizosaccharomyces octosporus sequences. Wilmar-P also contained a ∼100kb duplication in the right arm of chromosome III, a region harboring ght5⁺, the predominant hexose transporter encoding gene. Transcriptomic analysis of Wilmar-P grown in molasses revealed strong downregulation of core environmental stress response genes and upregulation of hexose transporters and drug efflux pumps compared to laboratory S. pombe. Finally, examination of the regulatory network of Scr1, which is involved in the regulation of several genes differentially expressed on molasses, revealed expanded binding of this transcription factor in Wilmar-P compared to laboratory S. pombe in the molasses condition. Together our results point to both genomic plasticity and transcriptomic adaptation as mechanisms driving phenotypic adaptation of Wilmar-P to the molasses environment and therefore adds to our understanding of genetic diversity within industrial fission yeast strains and the capacity of this strain for commercial scale bioethanol production.

Intraspecific Diversity of Fission Yeast Mitochondrial Genomes

The fission yeast Schizosaccharomyces pombe is an important model organism, but its natural diversity and evolutionary history remain under-studied. In particular, the population genomics of the S. pombe mitochondrial genome (mitogenome) has not been thoroughly investigated. Here, we assembled the complete circular-mapping mitogenomes of 192 S. pombe isolates de novo, and found that these mitogenomes belong to 69 nonidentical sequence types ranging from 17,618 to 26,910 bp in length. Using the assembled mitogenomes, we identified 20 errors in the reference mitogenome and discovered two previously unknown mitochondrial introns. Analyzing sequence diversity of these 69 types of mitogenomes revealed two highly distinct clades, with only three mitogenomes exhibiting signs of inter-clade recombination. This diversity pattern suggests that currently available S. pombe isolates descend from two long-separated ancestral lineages. This conclusion is corroborated by the diversity pattern of the recombination-repressed K-region located between donor mating-type loci mat2 and mat3 in the nuclear genome. We estimated that the two ancestral S. pombe lineages diverged about 31 million generations ago. These findings shed new light on the evolution of S. pombe and the data sets generated in this study will facilitate future research on genome evolution.

Ancestral Admixture Is the Main Determinant of Global Biodiversity in Fission Yeast

Mutation and recombination are key evolutionary processes governing phenotypic variation and reproductive isolation. We here demonstrate that biodiversity within all globally known strains of Schizosaccharomyces pombe arose through admixture between two divergent ancestral lineages. Initial hybridization was inferred to have occurred ∼20-60 sexual outcrossing generations ago consistent with recent, human-induced migration at the onset of intensified transcontinental trade. Species-wide heritable phenotypic variation was explained near-exclusively by strain-specific arrangements of alternating ancestry components with evidence for transgressive segregation. Reproductive compatibility between strains was likewise predicted by the degree of shared ancestry. To assess the genetic determinants of ancestry block distribution across the genome, we characterized the type, frequency, and position of structural genomic variation using nanopore and single-molecule real-time sequencing. Despite being associated with double-strand break initiation points, over 800 segregating structural variants exerted overall little influence on the introgression landscape or on reproductive compatibility between strains. In contrast, we found strong ancestry disequilibrium consistent with negative epistatic selection shaping genomic ancestry combinations during the course of hybridization. This study provides a detailed, experimentally tractable example that genomes of natural populations are mosaics reflecting different evolutionary histories. Exploiting genome-wide heterogeneity in the history of ancestral recombination and lineage-specific mutations sheds new light on the population history of S. pombe and highlights the importance of hybridization as a creative force in generating biodiversity.

Fitness Landscape of the Fission Yeast Genome

The relationship between DNA sequence, biochemical function, and molecular evolution is relatively well-described for protein-coding regions of genomes, but far less clear in noncoding regions, particularly, in eukaryote genomes. In part, this is because we lack a complete description of the essential noncoding elements in a eukaryote genome. To contribute to this challenge, we used saturating transposon mutagenesis to interrogate the Schizosaccharomyces pombe genome. We generated 31 million transposon insertions, a theoretical coverage of 2.4 insertions per genomic site. We applied a five-state hidden Markov model (HMM) to distinguish insertion-depleted regions from insertion biases. Both raw insertion-density and HMM-defined fitness estimates showed significant quantitative relationships to gene knockout fitness, genetic diversity, divergence, and expected functional regions based on transcription and gene annotations. Through several analyses, we conclude that transposon insertions produced fitness effects in 66-90% of the genome, including substantial portions of the noncoding regions. Based on the HMM, we estimate that 10% of the insertion depleted sites in the genome showed no signal of conservation between species and were weakly transcribed, demonstrating limitations of comparative genomics and transcriptomics to detect functional units. In this species, 3'- and 5'-untranslated regions were the most prominent insertion-depleted regions that were not represented in measures of constraint from comparative genomics. We conclude that the combination of transposon mutagenesis, evolutionary, and biochemical data can provide new insights into the relationship between genome function and molecular evolution.

Killer Meiotic Drive and Dynamic Evolution of the wtf Gene Family

Natural selection works best when the two alleles in a diploid organism are transmitted to offspring at equal frequencies. Despite this, selfish loci known as meiotic drivers that bias their own transmission into gametes are found throughout eukaryotes. Drive is thought to be a powerful evolutionary force, but empirical evolutionary analyses of drive systems are limited by low numbers of identified meiotic drive genes. Here, we analyze the evolution of the wtf gene family of Schizosaccharomyces pombe that contains both killer meiotic drive genes and suppressors of drive. We completed assemblies of all wtf genes for two S. pombe isolates, as well as a subset of wtf genes from over 50 isolates. We find that wtf copy number can vary greatly between isolates and that amino acid substitutions, expansions and contractions of DNA sequence repeats, and nonallelic gene conversion between family members all contribute to dynamic wtf gene evolution. This work demonstrates the power of meiotic drive to foster rapid evolution and identifies a recombination mechanism through which transposons can indirectly mobilize meiotic drivers.

Asymmetric diversification of mating pheromones in fission yeast

In fungi, mating between partners depends on the molecular recognition of two peptidyl mating pheromones by their respective receptors. The fission yeast Schizosaccharomyces pombe (Sp) has two mating types, Plus (P) and Minus (M). The mating pheromones P-factor and M-factor, secreted by P and M cells, are recognized by the receptors mating type auxiliary minus 2 (Mam2) and mating type auxiliary plus 3 (Map3), respectively. Our recent study demonstrated that a few mutations in both M-factor and Map3 can trigger reproductive isolation in S. pombe. Here, we explored the mechanism underlying reproductive isolation through genetic changes of pheromones/receptors in nature. We investigated the diversity of genes encoding the pheromones and their receptor in 150 wild S. pombe strains. Whereas the amino acid sequences of M-factor and Map3 were completely conserved, those of P-factor and Mam2 were very diverse. In addition, the P-factor gene contained varying numbers of tandem repeats of P-factor (4–8 repeats). By exploring the recognition specificity of pheromones between S. pombe and its close relative Schizosaccharomyces octosporus (So), we found that So-M-factor did not have an effect on S. pombe P cells, but So-P-factor had a partial effect on S. pombe M cells. Thus, recognition of M-factor seems to be stringent, whereas that of P-factor is relatively relaxed. We speculate that asymmetric diversification of the two pheromones might be facilitated by the distinctly different specificities of the two receptors. Our findings suggest that M-factor communication plays an important role in defining the species, whereas P-factor communication is able to undergo a certain degree of flexible adaptation–perhaps as a first step toward prezygotic isolation in S. pombe.

An asymmetric pheromone/receptor system in the fission yeast Schizosaccharomyces pombe might allow flexible adaptation of pheromones to mutational changes while maintaining stringent recognition for mating partners, perhaps as a first step toward prezygotic mating isolation.

The emergence of a new species might occur when two groups can no longer mate. Although such reproductive isolation is considered a key evolutionary process, the mechanisms by which it actually occurs have been confined to conjecture. The two sexes (Plus [P] and Minus [M]) of S. pombe each secrete a pheromone (P-factor and M-factor), which binds to a corresponding receptor (mating type auxiliary minus 2 [Mam2] and mating type auxiliary plus 3 [Map3]) on cells of the opposite sex. The interaction between a pheromone and its receptor is essential for successful mating. Here, we explored conservation of the mating pheromone communication system among 150 wild S. pombe strains of different geographical origins and the closely related species S. octosporus. We found that 1) the M-factor/Map3 interaction was completely conserved, whereas the P-factor/Mam2 interaction was very diverse in the strains investigated, and 2) most of the P-factor variants were functional across species. Thus, we have revealed an asymmetric pheromone/receptor system in fungal mating: namely, whereas M-factor communication operates extremely stringently, P-factor communication has the flexibility to create variations, perhaps facilitating prezygotic isolation in S. pombe.

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.

Uncovering Natural Longevity Alleles from Intercrossed Pools of Aging Fission Yeast Cells

Chronological lifespan of non-dividing yeast cells is a quantitative trait that reflects cellular aging. By monitoring allele frequencies in aging segregant pools, Ellis et al. uncover regulatory variants in the 5'-untranslated regions of two genes...

Quantitative traits often show large variation caused by multiple genetic factors . One such trait is the chronological lifespan of non-dividing yeast cells, serving as a model for cellular aging. Screens for genetic factors involved in aging typically assay mutants of protein-coding genes. To identify natural genetic variants contributing to cellular aging, we exploited two strains of the fission yeast, Schizosaccharomyces pombe, that differ in chronological lifespan. We generated segregant pools from these strains and subjected them to advanced intercrossing over multiple generations to break up linkage groups. We chronologically aged the intercrossed segregant pool, followed by genome sequencing at different times to detect genetic variants that became reproducibly enriched as a function of age. A region on Chromosome II showed strong positive selection during aging. Based on expected functions, two candidate variants from this region in the long-lived strain were most promising to be causal: small insertions and deletions in the 5′-untranslated regions of ppk31 and SPBC409.08. Ppk31 is an ortholog of Rim15, a conserved kinase controlling cell proliferation in response to nutrients, while SPBC409.08 is a predicted spermine transmembrane transporter. Both Rim15 and the spermine-precursor, spermidine, are implicated in aging as they are involved in autophagy-dependent lifespan extension. Single and double allele replacement suggests that both variants, alone or combined, have subtle effects on cellular longevity. Furthermore, deletion mutants of both ppk31 and SPBC409.08 rescued growth defects caused by spermidine. We propose that Ppk31 and SPBC409.08 may function together to modulate lifespan, thus linking Rim15/Ppk31 with spermidine metabolism.

Repeated evolution of self-compatibility for reproductive assurance

Sexual reproduction in eukaryotes requires the fusion of two compatible gametes of opposite sexes or mating types. To meet the challenge of finding a mating partner with compatible gametes, evolutionary mechanisms such as hermaphroditism and self-fertilization have repeatedly evolved. Here, by combining the insights from comparative genomics, computer simulations and experimental evolution in fission yeast, we shed light on the conditions promoting separate mating types or self-compatibility by mating-type switching. Analogous to multiple independent transitions between switchers and non-switchers in natural populations mediated by structural genomic changes, novel switching genotypes readily evolved under selection in the experimental populations. Detailed fitness measurements accompanied by computer simulations show the benefits and costs of switching during sexual and asexual reproduction, governing the occurrence of both strategies in nature. Our findings illuminate the trade-off between the benefits of reproductive assurance and its fitness costs under benign conditions facilitating the evolution of self-compatibility.

Mating-type switching enables self-compatible reproduction in fungi, but switching ability is variable even within species. Here, the authors find de novo evolution of switching genotypes in experimentally evolved fission yeast populations and show a trade-off between mating success and growth.

The natural diversity and ecology of fission yeast

While the fission yeast is a powerful model of eukaryote biology, there have been few studies of quantitative genetics, phenotypic or genetic diversity. Here I survey the small collection of fission yeast diversity research. I discuss what we can infer about the ecology and origins of Schizosaccharomyces pombe from microbiology field studies and the few strains that have been collected.

High-Quality de Novo Genome Assembly of the Dekkera bruxellensis Yeast Using Nanopore MinION Sequencing

Genetic variation in natural populations represents the raw material for phenotypic diversity. Species-wide characterization of genetic variants is crucial to have a deeper insight into the genotype-phenotype relationship. With the advent of new sequencing strategies and more recently the release of long-read sequencing platforms, it is now possible to explore the genetic diversity of any nonmodel organisms, representing a fundamental resource for biological research. In the frame of population genomic surveys, a first step is to obtain the complete sequence and high-quality assembly of a reference genome. Here, we sequenced and assembled a reference genome of the nonconventional Dekkera bruxellensis yeast. While this species is a major cause of wine spoilage, it paradoxically contributes to the specific flavor profile of some Belgium beers. In addition, an extreme karyotype variability is observed across natural isolates, highlighting that D. bruxellensis genome is very dynamic. The whole genome of the D. bruxellensis UMY321 isolate was sequenced using a combination of Nanopore long-read and Illumina short-read sequencing data. We generated the most complete and contiguous de novo assembly of D. bruxellensis to date and obtained a first glimpse into the genomic variability within this species by comparing the sequences of several isolates. This genome sequence is therefore of high value for population genomic surveys and represents a reference to study genome dynamic in this yeast species.

Population Genomics of Paramecium Species

Population-genomic analyses are essential to understanding factors shaping genomic variation and lineage-specific sequence constraints. The dearth of such analyses for unicellular eukaryotes prompted us to assess genomic variation in Paramecium, one of the most well-studied ciliate genera. The Paramecium aurelia complex consists of ∼15 morphologically indistinguishable species that diverged subsequent to two rounds of whole-genome duplications (WGDs, as long as 320 MYA) and possess extremely streamlined genomes. We examine patterns of both nuclear and mitochondrial polymorphism, by sequencing whole genomes of 10-13 worldwide isolates of each of three species belonging to the P. aurelia complex: P. tetraurelia, P. biaurelia, P. sexaurelia, as well as two outgroup species that do not share the WGDs: P. caudatum and P. multimicronucleatum. An apparent absence of global geographic population structure suggests continuous or recent dispersal of Paramecium over long distances. Intergenic regions are highly constrained relative to coding sequences, especially in P. caudatum and P. multimicronucleatum that have shorter intergenic distances. Sequence diversity and divergence are reduced up to ∼100-150 bp both upstream and downstream of genes, suggesting strong constraints imposed by the presence of densely packed regulatory modules. In addition, comparison of sequence variation at non-synonymous and synonymous sites suggests similar recent selective pressures on paralogs within and orthologs across the deeply diverging species. This study presents the first genome-wide population-genomic analysis in ciliates and provides a valuable resource for future studies in evolutionary and functional genetics in Paramecium.

Budding Yeast Strains and Genotype-Phenotype Mapping

A small number of well-studied laboratory strains of Saccharomyces cerevisiae, mostly derived from S288C, are used in yeast research. Although powerful, studies for understanding S288C do not always capture the phenotypic essence or the genetic complexity of S. cerevisiae biology. This is particularly problematic for multilocus phenotypes identified in laboratory strains because these loci have never been jointly exposed to natural selection and the corresponding phenotypes do not represent optimization for any particular purpose or environment. Isolation and sequencing of new natural yeast strains also reveal that the total sequence diversity of the S. cerevisiae global population is poorly sampled in common laboratory strains. Here we discuss methodologies required for using the natural genetic variation in yeast to complete a genotype-phenotype map.

Molecular Population Genetics

Molecular population genetics aims to explain genetic variation and molecular evolution from population genetics principles. The field was born 50 years ago with the first measures of genetic variation in allozyme loci, continued with the nucleotide sequencing era, and is currently in the era of population genomics. During this period, molecular population genetics has been revolutionized by progress in data acquisition and theoretical developments. The conceptual elegance of the neutral theory of molecular evolution or the footprint carved by natural selection on the patterns of genetic variation are two examples of the vast number of inspiring findings of population genetics research. Since the inception of the field, Drosophila has been the prominent model species: molecular variation in populations was first described in Drosophila and most of the population genetics hypotheses were tested in Drosophila species. In this review, we describe the main concepts, methods, and landmarks of molecular population genetics, using the Drosophila model as a reference. We describe the different genetic data sets made available by advances in molecular technologies, and the theoretical developments fostered by these data. Finally, we review the results and new insights provided by the population genomics approach, and conclude by enumerating challenges and new lines of inquiry posed by increasingly large population scale sequence data.

Transient structural variations have strong effects on quantitative traits and reproductive isolation in fission yeast

Large structural variations (SVs) within genomes are more challenging to identify than smaller genetic variants but may substantially contribute to phenotypic diversity and evolution. We analyse the effects of SVs on gene expression, quantitative traits and intrinsic reproductive isolation in the yeast Schizosaccharomyces pombe. We establish a high-quality curated catalogue of SVs in the genomes of a worldwide library of S. pombe strains, including duplications, deletions, inversions and translocations. We show that copy number variants (CNVs) show a variety of genetic signals consistent with rapid turnover. These transient CNVs produce stoichiometric effects on gene expression both within and outside the duplicated regions. CNVs make substantial contributions to quantitative traits, most notably intracellular amino acid concentrations, growth under stress and sugar utilization in winemaking, whereas rearrangements are strongly associated with reproductive isolation. Collectively, these findings have broad implications for evolution and for our understanding of quantitative traits including complex human diseases.

Negative epistasis: a route to intraspecific reproductive isolation in yeast?

Exploring the molecular bases of intraspecific reproductive isolation captures the ongoing phenotypic consequences of genetic divergence and provides insights into the early onset of speciation. Recent species-wide surveys using natural populations of yeasts demonstrated that intrinsic post-zygotic reproductive isolation segregates readily within the same species, and revealed the multiplicity of the genetic mechanisms underlying such processes. These advances deepened our current understandings and opened further perspectives regarding the complete picture of molecular and evolutionary origins driving the onset of intraspecific reproductive isolation in yeasts.

Population genomics of yeasts: towards a comprehensive view across a broad evolutionary scale

With the advent of high-throughput technologies for sequencing, the complete description of the genetic variation that occurs in populations, also known as population genomics, is foreseeable but far from being reached. Explaining the forces that govern patterns of genetic variation is essential to elucidate the evolutionary history of species. Genetic variation results from a wide assortment of evolutionary forces, among which mutation, selection, recombination and drift play major roles in shaping genomes. In addition, exploring the genetic variation within a population also corresponds to the first step towards dissecting the genotype-phenotype relationship. In this context, yeast species are of particular interest because they represent a unique resource for studying the evolution of intraspecific genetic diversity in a phylum spanning a broad evolutionary scale. Here, we briefly review recent progress in yeast population genomics and provide some perspective on this rapidly evolving field. In fact, we truly believe that it is of interest to supplement comparative and early population genomic studies with the deep sequencing of more extensive sets of individuals from the same species. In parallel, it would be more than valuable to uncover the intraspecific variation of a large number of unexplored species, including those that are closely and more distantly related. Altogether, these data would enable substantially more powerful genomic scans for functional dissection.

Bulk Segregant Analysis Reveals the Genetic Basis of a Natural Trait Variation in Fission Yeast

Although the fission yeast Schizosaccharomyces pombe is a well-established model organism, studies of natural trait variations in this species remain limited. To assess the feasibility of segregant-pool-based mapping of phenotype-causing genes in natural strains of fission yeast, we investigated the cause of a maltose utilization defect (Mal(-)) of the S. pombe strain CBS5557 (originally known as Schizosaccharomyces malidevorans). Analyzing the genome sequence of CBS5557 revealed 955 nonconservative missense substitutions, and 61 potential loss-of-function variants including 47 frameshift indels, 13 early stop codons, and 1 splice site mutation. As a side benefit, our analysis confirmed 146 sequence errors in the reference genome and improved annotations of 27 genes. We applied bulk segregant analysis to map the causal locus of the Mal(-) phenotype. Through sequencing the segregant pools derived from a cross between CBS5557 and the laboratory strain, we located the locus to within a 2.23-Mb chromosome I inversion found in most S. pombe isolates including CBS5557. To map genes within the inversion region that occupies 18% of the genome, we created a laboratory strain containing the same inversion. Analyzing segregants from a cross between CBS5557 and the inversion-containing laboratory strain narrowed down the locus to a 200-kb interval and led us to identify agl1, which suffers a 5-bp deletion in CBS5557, as the causal gene. Interestingly, loss of agl1 through a 34-kb deletion underlies the Mal(-) phenotype of another S. pombe strain CGMCC2.1628. This work adapts and validates the bulk segregant analysis method for uncovering trait-gene relationship in natural fission yeast strains.

The Spontaneous Mutation Rate in the Fission Yeast Schizosaccharomyces pombe

The rate at which new mutations arise in the genome is a key factor in the evolution and adaptation of species. Here we describe the rate and spectrum of spontaneous mutations for the fission yeast Schizosaccharomyces pombe, a key model organism with many similarities to higher eukaryotes. We undertook an ∼1700-generation mutation accumulation (MA) experiment with a haploid S. pombe, generating 422 single-base substitutions and 119 insertion-deletion mutations (indels) across the 96 replicates. This equates to a base-substitution mutation rate of 2.00 × 10(-10) mutations per site per generation, similar to that reported for the distantly related budding yeast Saccharomyces cerevisiae. However, these two yeast species differ dramatically in their spectrum of base substitutions, the types of indels (S. pombe is more prone to insertions), and the pattern of selection required to counteract a strong AT-biased mutation rate. Overall, our results indicate that GC-biased gene conversion does not play a major role in shaping the nucleotide composition of the S. pombe genome and suggest that the mechanisms of DNA maintenance may have diverged significantly between fission and budding yeasts. Unexpectedly, CpG sites appear to be excessively liable to mutation in both species despite the likely absence of DNA methylation.

The genomic and phenotypic diversity of Schizosaccharomyces pombe

Natural variation within species reveals aspects of genome evolution and function. The fission yeast Schizosaccharomyces pombe is an important model for eukaryotic biology, but researchers typically use one standard laboratory strain. To extend the usefulness of this model, we surveyed the genomic and phenotypic variation in 161 natural isolates. We sequenced the genomes of all strains, finding moderate genetic diversity (π = 3 × 10(-3) substitutions/site) and weak global population structure. We estimate that dispersal of S. pombe began during human antiquity (∼340 BCE), and ancestors of these strains reached the Americas at ∼1623 CE. We quantified 74 traits, finding substantial heritable phenotypic diversity. We conducted 223 genome-wide association studies, with 89 traits showing at least one association. The most significant variant for each trait explained 22% of the phenotypic variance on average, with indels having larger effects than SNPs. This analysis represents a rich resource to examine genotype-phenotype relationships in a tractable model.