Population genomic analysis reveals highly conserved mitochondrial genomes in the yeast species Lachancea thermotolerans

The increasing availability of mitochondrial (mt) sequence data from various yeasts provides a tool to study genomic evolution within and between different species. While the genomes from a range of lineages are available, there is a lack of information concerning intraspecific mtDNA diversity. Here, we analyzed the mt genomes of 50 strains from Lachancea thermotolerans, a protoploid yeast species that has been isolated from several locations (Europe, Asia, Australia, South Africa, and North / South America) and ecological sources (fruit, tree exudate, plant material, and grape and agave fermentations). Protein-coding genes from the mtDNA were used to construct a phylogeny, which reflected a similar, yet less resolved topology than the phylogenetic tree of 50 nuclear genes. In comparison to its sister species Lachancea kluyveri, L. thermotolerans has a smaller mt genome. This is due to shorter intergenic regions and fewer introns, of which the latter are only found in COX1. We revealed that L. kluyveri and L. thermotolerans share similar levels of intraspecific divergence concerning the nuclear genomes. However, L. thermotolerans has a more highly conserved mt genome with the coding regions characterized by low rates of nonsynonymous substitution. Thus, in the mt genomes of L. thermotolerans, stronger purifying selection and lower mutation rates potentially shape genome diversity in contract to what was found for L. kluyveri, demonstrating that the factors driving mt genome evolution are different even between closely related species

Comprehensive survey of condition-specific reproductive isolation reveals genetic incompatibility in yeast

Genetic variation within a species could cause negative epistasis leading to reduced hybrid fitness and post-zygotic reproductive isolation. Recent studies in yeasts revealed chromosomal rearrangements as a major mechanism dampening intraspecific hybrid fertility on rich media. Here, by analysing a large number of Saccharomyces cerevisiae crosses on different culture conditions, we show environment-specific genetic incompatibility segregates readily within yeast and contributes to reproductive isolation. Over 24% (117 out of 481) of cases tested show potential epistasis, among which 6.7% (32 out of 481) are severe, with at least 20% of progeny loss on tested conditions. Based on the segregation patterns, we further characterize a two-locus Dobzhansky-Muller incompatibility case leading to offspring respiratory deficiency caused by nonsense mutation in a nuclear-encoding mitochondrial gene and tRNA suppressor. We provide evidence that this precise configuration could be adaptive in fluctuating environments, highlighting the role of ecological selection in the onset of genetic incompatibility and reproductive isolation in yeast

Differences in environmental stress response between yeasts is consistent with species-specific lifestyles

Defining how organisms respond to environmental change has always been an important step toward understating their adaptive capacity and physiology. Variation in transcription during stress has been widely described in model species, especially in the yeastSaccharomyces cerevisiae, which helped to shape general rules regarding how cells cope with environmental constraints as well as decipher the functions of many genes. Now, comparison of the environmental stress response (ESR) across species is essential to obtain a better insight into the common and species-specific features of stress defense. In this context, we explored the transcriptional landscape of the yeastLachancea kluyveri(formerlySaccharomyces kluyveri) in response to diverse stresses, using RNA-seq. We investigated variation in gene expression and observed a link between genetic plasticity and environmental sensitivity. We identified the ESR genes in this species and compared them to those already found inS. cerevisiae We observed common features between the two species as well as divergence in the regulatory networks involved. Interestingly, some changes were related to differences in species lifestyle. Thus, we were able to decipher how adaptation to stress has evolved among different yeast species. Finally, by analyzing patterns of coexpression, we were able to propose potential biological functions for 42% of genes and furthermore annotate 301 genes for which no function could be assigned by homology. This large dataset allowed for the characterization of the evolution of gene regulation and provides an efficient tool to assess gene function

Computational Models for Prediction of Yeast Strain Potential for Winemaking from Phenotypic Profiles

Saccharomyces cerevisiae strains from diverse natural habitats harbour a vast amount of phenotypic diversity, driven by interactions between yeast and the respective environment. In grape juice fermentations, strains are exposed to a wide array of biotic and abiotic stressors, which may lead to strain selection and generate naturally arising strain diversity. Certain phenotypes are of particular interest for the winemaking industry and could be identified by screening of large number of different strains. The objective of the present work was to use data mining approaches to identify those phenotypic tests that are most useful to predict a strain's potential for winemaking. We have constituted a S. cerevisiae collection comprising 172 strains of worldwide geographical origins or technological applications. Their phenotype was screened by considering 30 physiological traits that are important from an oenological point of view. Growth in the presence of potassium bisulphite, growth at 40 degrees C, and resistance to ethanol were mostly contributing to strain variability, as shown by the principal component analysis. In the hierarchical clustering of phenotypic profiles the strains isolated from the same wines and vineyards were scattered throughout all clusters, whereas commercial winemaking strains tended to co-cluster. Mann-Whitney test revealed significant associations between phenotypic results and strain's technological application or origin. Naive Bayesian classifier identified 3 of the 30 phenotypic tests of growth in iprodion (0.05 mg/mL), cycloheximide (0.1 mu g/mL) and potassium bisulphite (150 mg/mL) that provided most information for the assignment of a strain to the group of commercial strains. The probability of a strain to be assigned to this group was 27% using the entire phenotypic profile and increased to 95%, when only results from the three tests were considered. Results show the usefulness of computational approaches to simplify strain selection procedures.Ines Mendes and Ricardo Franco-Duarte are recipients of a fellowship from the Portuguese Science Foundation, FCT (SFRH/BD/74798/2010, SFRH/BD/48591/2008, respectively) and Joao Drumonde-Neves is recipient of a fellowship from the Azores government (M3.1.2/F/006/2008 (DRCT)). Financial support was obtained from FEDER funds through the program COMPETE and by national funds through FCT by the projects FCOMP-01-0124-008775 (PTDC/AGR-ALI/103392/2008) and PTDC/AGR-ALI/121062/2010. Lan Umek and Blaz Zupan acknowledge financial support from Slovene Research Agency (P2-0209). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.info:eu-repo/semantics/publishedVersio

Single-cell phenomics reveals intra-species variation of phenotypic noise in yeast

BACKGROUND: Most quantitative measures of phenotypic traits represent macroscopic contributions of large numbers of cells. Yet, cells of a tissue do not behave similarly, and molecular studies on several organisms have shown that regulations can be highly stochastic, sometimes generating diversified cellular phenotypes within tissues. Phenotypic noise, defined here as trait variability among isogenic cells of the same type and sharing a common environment, has therefore received a lot of attention. Given the potential fitness advantage provided by phenotypic noise in fluctuating environments, the possibility that it is directly subjected to evolutionary selection is being considered. For selection to act, phenotypic noise must differ between contemporary genotypes. Whether this is the case or not remains, however, unclear because phenotypic noise has very rarely been quantified in natural populations. RESULTS: Using automated image analysis, we describe here the phenotypic diversity of S. cerevisiae morphology at single-cell resolution. We profiled hundreds of quantitative traits in more than 1,000 cells of 37 natural strains, which represent various geographical and ecological origins of the species. We observed abundant trait variation between strains, with no correlation with their ecological origin or population history. Phenotypic noise strongly depended on the strain background. Noise variation was largely trait-specific (specific strains showing elevated noise for subset of traits) but also global (a few strains displaying elevated noise for many unrelated traits). CONCLUSIONS: Our results demonstrate that phenotypic noise does differ quantitatively between natural populations. This supports the possibility that, if noise is adaptive, microevolution may tune it in the wild. This tuning may happen on specific traits or by varying the degree of global phenotypic buffering

Whole-Genome Sequencing and Intraspecific Analysis of the Yeast Species Lachancea quebecensis

The gold standard in yeast population genomics has been the model organism Saccharomyces cerevisiae. However, the exploration of yeast species outside the Saccharomyces genus is essential to broaden the understanding of genome evolution. Here, we report the analyses of whole-genome sequences of nineisolates from the recently described yeast species Lachancea quebecensis. The genome of one isolate was assembled and annotated, and the intraspecific variability within L. quebecensis was surveyed by comparing the sequences from the eight other isolates to this reference sequence. Our study revealed that these strains harbor genomes with an average nucleotide diversity of pi = 2 x 10(-3) which is slightly lower, although on the same order of magnitude, as that previously determined for S. cerevisiae (pi = 4 x 10(-3)). Our results show that even though these isolates were all obtained from a relatively isolated geographic location, the same ecological source, and represent a smaller sample size than is available for S. cerevisiae, the levels of divergence are similar to those observed in this model species. This divergence is essentially linked to the presence of two distinct clusters delineated according to geographic location. However, even with relatively similar ranges of genome divergence, L. quebecensis has an extremely low global phenotypic variance of 0.062 compared with 0.59 previously determined in S. cerevisiae

Heterologous Stop Codon Readthrough of Metazoan Readthrough Candidates in Yeast

Recent analysis of genomic signatures in mammals, flies, and worms indicates that functional translational stop codon readthrough is considerably more abundant in metazoa than previously recognized, but this analysis provides only limited clues about the function or mechanism of readthrough. If an mRNA known to be read through in one species is also read through in another, perhaps these questions can be studied in a simpler setting. With this end in mind, we have investigated whether some of the readthrough genes in human, fly, and worm also exhibit readthrough when expressed in S. cerevisiae. We found that readthrough was highest in a gene with a post-stop hexamer known to trigger readthrough, while other metazoan readthrough genes exhibit borderline readthrough in S. cerevisiae.National Institutes of Health (U.S.) (5U54HG004555-03

Correction: AGAPE (Automated Genome Analysis PipelinE) for Pan-Genome Analysis of Saccharomyces cerevisiae

The characterization and public release of genome sequences from thousands of organisms is expanding the scope for genetic variation studies. However, understanding the phenotypic consequences of genetic variation remains a challenge in eukaryotes due to the complexity of the genotype-phenotype map. One approach to this is the intensive study of model systems for which diverse sources of information can be accumulated and integrated. Saccharomyces cerevisiae is an extensively studied model organism, with well-known protein functions and thoroughly curated phenotype data. To develop and expand the available resources linking genomic variation with function in yeast, we aim to model the pan-genome of S. cerevisiae. To initiate the yeast pan-genome, we newly sequenced or re-sequenced the genomes of 25 strains that are commonly used in the yeast research community using advanced sequencing technology at high quality. We also developed a pipeline for automated pan-genome analysis, which integrates the steps of assembly, annotation, and variation calling. To assign strain-specific functional annotations, we identified genes that were not present in the reference genome. We classified these according to their presence or absence across strains and characterized each group of genes with known functional and phenotypic features. The functional roles of novel genes not found in the reference genome and associated with strains or groups of strains appear to be consistent with anticipated adaptations in specific lineages. As more S. cerevisiae strain genomes are released, our analysis can be used to collate genome data and relate it to lineage-specific patterns of genome evolution. Our new tool set will enhance our understanding of genomic and functional evolution in S. cerevisiae, and will be available to the yeast genetics and molecular biology community

On the complexity of the Saccharomyces bayanus taxon: hybridization and potential hybrid speciation

Although the genus Saccharomyces has been thoroughly studied, some species in the genus has not yet been accurately resolved; an example is S. bayanus, a taxon that includes genetically diverse lineages of pure and hybrid strains. This diversity makes the assignation and classification of strains belonging to this species unclear and controversial. They have been subdivided by some authors into two varieties (bayanus and uvarum), which have been raised to the species level by others. In this work, we evaluate the complexity of 46 different strains included in the S. bayanus taxon by means of PCR-RFLP analysis and by sequencing of 34 gene regions and one mitochondrial gene. Using the sequence data, and based on the S. bayanus var. bayanus reference strain NBRC 1948, a hypothetical pure S. bayanus was reconstructed for these genes that showed alleles with similarity values lower than 97% with the S. bayanus var. uvarum strain CBS 7001, and of 99¿100% with the non S. cerevisiae portion in S. pastorianus Weihenstephan 34/70 and with the new species S. eubayanus. Among the S. bayanus strains under study, different levels of homozygosity, hybridization and introgression were found; however, no pure S. bayanus var. bayanus strain was identified. These S. bayanus hybrids can be classified into two types: homozygous (type I) and heterozygous hybrids (type II), indicating that they have been originated by different hybridization processes. Therefore, a putative evolutionary scenario involving two different hybridization events between a S. bayanus var. uvarum and unknown European S. eubayanus-like strains can be postulated to explain the genomic diversity observed in our S. bayanus var. bayanus strains

Systematic Identification of Balanced Transposition Polymorphisms in Saccharomyces cerevisiae

High-throughput techniques for detecting DNA polymorphisms generally do not identify changes in which the genomic position of a sequence, but not its copy number, varies among individuals. To explore such balanced structural polymorphisms, we used array-based Comparative Genomic Hybridization (aCGH) to conduct a genome-wide screen for single-copy genomic segments that occupy different genomic positions in the standard laboratory strain of Saccharomyces cerevisiae (S90) and a polymorphic wild isolate (Y101) through analysis of six tetrads from a cross of these two strains. Paired-end high-throughput sequencing of Y101 validated four of the predicted rearrangements. The transposed segments contained one to four annotated genes each, yet crosses between S90 and Y101 yielded mostly viable tetrads. The longest segment comprised 13.5 kb near the telomere of chromosome XV in the S288C reference strain and Southern blotting confirmed its predicted location on chromosome IX in Y101. Interestingly, inter-locus crossover events between copies of this segment occurred at a detectable rate. The presence of low-copy repetitive sequences at the junctions of this segment suggests that it may have arisen through ectopic recombination. Our methodology and findings provide a starting point for exploring the origins, phenotypic consequences, and evolutionary fate of this largely unexplored form of genomic polymorphism