Search

Abstract Background Cost-effective production of biofuels from lignocellulose requires the fermentation of d-xylose. Many yeast species within and closely related to the genera Spathaspora and Scheffersomyces (both of the order Serinales) natively assimilate and ferment xylose. Other species consume xylose inefficiently, leading to extracellular accumulation of xylitol. Xylitol excretion is thought to be due to the different cofactor requirements of the first two steps of xylose metabolism. Xylose reductase (XR) generally uses NADPH to reduce xylose to xylitol, while xylitol dehydrogenase (XDH) generally uses NAD+ to oxidize xylitol to xylulose, creating an imbalanced redox pathway. This imbalance is thought to be particularly consequential in hypoxic or anoxic environments. Results We screened the growth of xylose-fermenting yeast species in high and moderate aeration and identified both ethanol producers and xylitol producers. Selected species were further characterized for their XR and XDH cofactor preferences by enzyme assays and gene expression patterns by RNA-Seq. Our data revealed that xylose metabolism is more redox balanced in some species, but it is strongly affected by oxygen levels. Under high aeration, most species switched from ethanol production to xylitol accumulation, despite the availability of ample oxygen to accept electrons from NADH. This switch was followed by decreases in enzyme activity and the expression of genes related to xylose metabolism, suggesting that bottlenecks in xylose fermentation are not always due to cofactor preferences. Finally, we expressed XYL genes from multiple Scheffersomyces species in a strain of Saccharomyces cerevisiae. Recombinant S. cerevisiae expressing XYL1 from Scheffersomyces xylosifermentans, which encodes an XR without a cofactor preference, showed improved anaerobic growth on xylose as the primary carbon source compared to S. cerevisiae strain expressing XYL genes from Scheffersomyces stipitis. Conclusion Collectively, our data do not support the hypothesis that xylitol accumulation occurs primarily due to differences in cofactor preferences between xylose reductase and xylitol dehydrogenase; instead, gene expression plays a major role in response to oxygen levels. We have also identified the yeast Sc. xylosifermentans as a potential source for genes that can be engineered into S. cerevisiae to improve xylose fermentation and biofuel production.

Summary: When Saccharomyces cerevisiae grows on mixtures of glucose and galactose, galactose utilization is repressed by glucose, and induction of the GAL gene network only occurs when glucose is exhausted. Contrary to reference GAL alleles, alternative alleles support faster growth on galactose, thus enabling distinct galactose utilization strategies maintained by balancing selection. Here, we report on new wild populations of Saccharomyces cerevisiae harboring alternative GAL versions and, for the first time, of Saccharomyces paradoxus alternative alleles. We also show that the non-functional GAL version found earlier in Saccharomyces kudriavzevii is phylogenetically related to the alternative versions, which constitutes a case of trans-specific maintenance of highly divergent alleles. Strains harboring the different GAL network variants show different levels of alleviation of glucose repression and growth proficiency on galactose. We propose that domestication involved specialization toward thriving in milk from a generalist ancestor partially adapted to galactose consumption in the plant niche.

IntroductionEukaryotic life depends on the functional elements encoded by both the nuclear genome and organellar genomes, such as those contained within the mitochondria. The content, size, and structure of the mitochondrial genome varies across organisms with potentially large implications for phenotypic variance and resulting evolutionary trajectories. Among yeasts in the subphylum Saccharomycotina, extensive differences have been observed in various species relative to the model yeast Saccharomyces cerevisiae, but mitochondrial genome sampling across many groups has been scarce, even as hundreds of nuclear genomes have become available.MethodsBy extracting mitochondrial assemblies from existing short-read genome sequence datasets, we have greatly expanded both the number of available genomes and the coverage across sparsely sampled clades.ResultsComparison of 353 yeast mitochondrial genomes revealed that, while size and GC content were fairly consistent across species, those in the genera Metschnikowia and Saccharomyces trended larger, while several species in the order Saccharomycetales, which includes S. cerevisiae, exhibited lower GC content. Extreme examples for both size and GC content were scattered throughout the subphylum. All mitochondrial genomes shared a core set of protein-coding genes for Complexes III, IV, and V, but they varied in the presence or absence of mitochondrially-encoded canonical Complex I genes. We traced the loss of Complex I genes to a major event in the ancestor of the orders Saccharomycetales and Saccharomycodales, but we also observed several independent losses in the orders Phaffomycetales, Pichiales, and Dipodascales. In contrast to prior hypotheses based on smaller-scale datasets, comparison of evolutionary rates in protein-coding genes showed no bias towards elevated rates among aerobically fermenting (Crabtree/Warburg-positive) yeasts. Mitochondrial introns were widely distributed, but they were highly enriched in some groups. The majority of mitochondrial introns were poorly conserved within groups, but several were shared within groups, between groups, and even across taxonomic orders, which is consistent with horizontal gene transfer, likely involving homing endonucleases acting as selfish elements.DiscussionAs the number of available fungal nuclear genomes continues to expand, the methods described here to retrieve mitochondrial genome sequences from these datasets will prove invaluable to ensuring that studies of fungal mitochondrial genomes keep pace with their nuclear counterparts.

Ploidy is an evolutionarily labile trait, and its variation across the tree of life has profound impacts on evolutionary trajectories and life histories. The immediate consequences and molecular causes of ploidy variation on organismal fitness are frequently less clear, although extreme mating type skews in some fungi hint at links between cell type and adaptive traits. Here, we report an unusual recurrent ploidy reduction in replicate populations of the budding yeast Saccharomyces eubayanus experimentally evolved for improvement of a key metabolic trait, the ability to use maltose as a carbon source. We find that haploids have a substantial, but conditional, fitness advantage in the absence of other genetic variation. Using engineered genotypes that decouple the effects of ploidy and cell type, we show that increased fitness is primarily due to the distinct transcriptional program deployed by haploid-like cell types, with a significant but smaller contribution from absolute ploidy. The link between cell-type specification and the carbon metabolism adaptation can be traced to the noncanonical regulation of a maltose transporter by a haploid-specific gene. This study provides novel mechanistic insight into the molecular basis of an environment-cell type fitness interaction and illustrates how selection on traits unexpectedly linked to ploidy states or cell types can drive karyotypic evolution in fungi.

Here, the authors describe the geographies, hosts, substrates, and phylogenetic relationships for 1,794 Saccharomyces strains. They provide insight into the genetic and phenotypic diversity in the genus, not seen through prior work focused on the model species Saccharomyces cerevisiae.

Abstract Background Lignocellulosic conversion residue (LCR) is the material remaining after deconstructed lignocellulosic biomass is subjected to microbial fermentation and treated to remove the biofuel. Technoeconomic analyses of biofuel refineries have shown that further microbial processing of this LCR into other bioproducts may help offset the costs of biofuel generation. Identifying organisms able to metabolize LCR is an important first step for harnessing the full chemical and economic potential of this material. In this study, we investigated the aerobic LCR utilization capabilities of 71 Streptomyces and 163 yeast species that could be engineered to produce valuable bioproducts. The LCR utilization by these individual microbes was compared to that of an aerobic mixed microbial consortium derived from a wastewater treatment plant as representative of a consortium with the highest potential for degrading the LCR components and a source of genetic material for future engineering efforts. Results We analyzed several batches of a model LCR by chemical oxygen demand (COD) and chromatography-based assays and determined that the major components of LCR were oligomeric and monomeric sugars and other organic compounds. Many of the Streptomyces and yeast species tested were able to grow in LCR, with some individual microbes capable of utilizing over 40% of the soluble COD. For comparison, the maximum total soluble COD utilized by the mixed microbial consortium was about 70%. This represents an upper limit on how much of the LCR could be valorized by engineered Streptomyces or yeasts into bioproducts. To investigate the utilization of specific components in LCR and have a defined media for future experiments, we developed a synthetic conversion residue (SynCR) to mimic our model LCR and used it to show lignocellulose-derived inhibitors (LDIs) had little effect on the ability of the Streptomyces species to metabolize SynCR. Conclusions We found that LCR is rich in carbon sources for microbial utilization and has vitamins, minerals, amino acids and other trace metabolites necessary to support growth. Testing diverse collections of Streptomyces and yeast species confirmed that these microorganisms were capable of growth on LCR and revealed a phylogenetic correlation between those able to best utilize LCR. Identification and quantification of the components of LCR enabled us to develop a synthetic LCR (SynCR) that will be a useful tool for examining how individual components of LCR contribute to microbial growth and as a substrate for future engineering efforts to use these microorganisms to generate valuable bioproducts.

Metabolic engineering strategies have been successfully implemented to improve the production of isobutanol, a next-generation biofuel, in Saccharomyces cerevisiae. Here, we explore how two of these strategies, pathway re-localization and redox cofactor-balancing, affect the performance and physiology of isobutanol producing strains. We equipped yeast with isobutanol cassettes which had either a mitochondrial or cytosolic localized isobutanol pathway and used either a redox-imbalanced (NADPH-dependent) or redox-balanced (NADH-dependent) ketol-acid reductoisomerase enzyme. We then conducted transcriptomic, proteomic and metabolomic analyses to elucidate molecular differences between the engineered strains. Pathway localization had a large effect on isobutanol production with the strain expressing the mitochondrial-localized enzymes producing 3.8-fold more isobutanol than strains expressing the cytosolic enzymes. Cofactor-balancing did not improve isobutanol titers and instead the strain with the redox-imbalanced pathway produced 1.5-fold more isobutanol than the balanced version, albeit at low overall pathway flux. Functional genomic analyses suggested that the poor performances of the cytosolic pathway strains were in part due to a shortage in cytosolic Fe–S clusters, which are required cofactors for the dihydroxyacid dehydratase enzyme. We then demonstrated that this cofactor limitation may be partially recovered by disrupting iron homeostasis with a fra2 mutation, thereby increasing cellular iron levels. The resulting isobutanol titer of the fra2 null strain harboring a cytosolic-localized isobutanol pathway outperformed the strain with the mitochondrial-localized pathway by 1.3-fold, demonstrating that both localizations can support flux to isobutanol.

ABSTRACT The unambiguous application of fungal names is important to communicate scientific findings. Names are critical for (clinical) diagnostics, legal compliance, and regulatory controls, such as biosafety, food security, quarantine regulations, and industrial applications. Consequently, the stability of the taxonomic system and the traceability of nomenclatural changes is crucial for a broad range of users and taxonomists. The unambiguous application of names is assured by the preservation of nomenclatural history and the physical organisms representing a name. Fungi are extremely diverse in terms of ecology, lifestyle, and methods of study. Predominantly unicellular fungi known as yeasts are usually investigated as living cultures. Methods to characterize yeasts include physiological (growth) tests and experiments to induce a sexual morph; both methods require viable cultures. Thus, the preservation and availability of viable reference cultures are important, and cultures representing reference material are cited in species descriptions. Historical surveys revealed drawbacks and inconsistencies between past practices and modern requirements as stated in the International Code of Nomenclature for Algae, Fungi, and Plants (ICNafp). Improper typification of yeasts is a common problem, resulting in a large number invalid yeast species names. With this opinion letter, we address the problem that culturable microorganisms, notably some fungi and algae, require specific provisions under the ICNafp. We use yeasts as a prominent example of fungi known from cultures. But viable type material is important not only for yeasts, but also for other cultivable Fungi that are characterized by particular morphological structures (a specific type of spores), growth properties, and secondary metabolites. We summarize potential proposals which, in our opinion, will improve the stability of fungal names, in particular by protecting those names for which the reference material can be traced back to the original isolate.

Replicate runs of maximum likelihood phylogenetic analyses can generate different tree topologies due to differences in parameters, such as random seeds. Here, Shen et al. demonstrate that replicate runs can generate substantially different tree topologies even with identical data and parameters.

Many industrial organisms are the result of recent or ancient allopolypoidy events. Here the authors iteratively combine the genomes of six yeast species to generate a viable hybrid.

[This corrects the article DOI: 10.1371/journal.pgen.1008304.].

During a survey of Nothofagus trees and their parasitic fungi in Andean Patagonia (Argentina), genetically distinct strains of Hanseniaspora were obtained from the sugar-containing stromata of parasitic Cyttaria spp. Phylogenetic analyses based on the single-gene sequences (encoding rRNA and actin) or on conserved, single-copy, orthologous genes from genome sequence assemblies revealed that these strains represent a new species closely related to Hanseniaspora valbyensis. Additionally, delimitation of this novel species was supported by genetic distance calculations using overall genome relatedness indices (OGRI) between the novel taxon and its closest relatives. To better understand the mode of speciation in Hanseniaspora, we examined genes that were retained or lost in the novel species in comparison to its closest relatives. These analyses show that, during diversification, this novel species and its closest relatives, H. valbyensis and Hanseniaspora jakobsenii, lost mitochondrial and other genes involved in the generation of precursor metabolites and energy, which could explain their slower growth and higher ethanol yields under aerobic conditions. Similarly, Hanseniaspora mollemarum lost the ability to sporulate, along with genes that are involved in meiosis and mating. Based on these findings, a formal description of the novel yeast species Hanseniaspora smithiae sp. nov. is proposed, with CRUB 1602H as the holotype.

Reverse ecology is the inference of ecological information from patterns of genomic variation. One rich, heretofore underutilized, source of ecologically relevant genomic information is codon optimality or adaptation. Bias toward codons that match the tRNA pool is robustly associated with high gene expression in diverse organisms, suggesting that codon optimization could be used in a reverse ecology framework to identify highly expressed, ecologically relevant genes. To test this hypothesis, we examined the relationship between optimal codon usage in the classic galactose metabolism (GAL) pathway and known ecological niches for 329 species of budding yeasts, a diverse subphylum of fungi. We find that optimal codon usage in the GAL pathway is positively correlated with quantitative growth on galactose, suggesting that GAL codon optimization reflects increased capacity to grow on galactose. Optimal codon usage in the GAL pathway is also positively correlated with human-associated ecological niches in yeasts of the CUG-Ser1 clade and with dairy-associated ecological niches in the family Saccharomycetaceae. For example, optimal codon usage of GAL genes is greater than 85% of all genes in the genome of the major human pathogen Candida albicans (CUG-Ser1 clade) and greater than 75% of genes in the genome of the dairy yeast Kluyveromyces lactis (family Saccharomycetaceae). We further find a correlation between optimization in the GALactose pathway genes and several genes associated with nutrient sensing and metabolism. This work suggests that codon optimization harbors information about the metabolic ecology of microbial eukaryotes. This information may be particularly useful for studying fungal dark matter-species that have yet to be cultured in the lab or have only been identified by genomic material.

Candida tropicalis is a human pathogen that primarily infects the immunocompromised. Whereas the genome of one isolate, C. tropicalis MYA-3404, was originally sequenced in 2009, there have been no large-scale, multi-isolate studies of the genetic and phenotypic diversity of this species. Here, we used whole genome sequencing and phenotyping to characterize 77 isolates of C. tropicalis from clinical and environmental sources from a variety of locations. We show that most C. tropicalis isolates are diploids with approximately 2-6 heterozygous variants per kilobase. The genomes are relatively stable, with few aneuploidies. However, we identified one highly homozygous isolate and six isolates of C. tropicalis with much higher heterozygosity levels ranging from 36-49 heterozygous variants per kilobase. Our analyses show that the heterozygous isolates represent two different hybrid lineages, where the hybrids share one parent (A) with most other C. tropicalis isolates, but the second parent (B or C) differs by at least 4% at the genome level. Four of the sequenced isolates descend from an AB hybridization, and two from an AC hybridization. The hybrids are MTLa/α heterozygotes. Hybridization, or mating, between different parents is therefore common in the evolutionary history of C. tropicalis. The new hybrids were predominantly found in environmental niches, including from soil. Hybridization is therefore unlikely to be associated with virulence. In addition, we used genotype-phenotype correlation and CRISPR-Cas9 editing to identify a genome variant that results in the inability of one isolate to utilize certain branched-chain amino acids as a sole nitrogen source.

The wild, cold-adapted parent of hybrid lager-brewing yeasts, Saccharomyces eubayanus, has a complex and understudied natural history. The exploration of this diversity can be used both to develop new brewing applications and to enlighten our understanding of the dynamics of yeast evolution in the wild. Here, we integrate whole genome sequence and phenotypic data of 200 S. eubayanus strains, the largest collection known to date. S. eubayanus has a multilayered population structure, consisting of two major populations that are further structured into six subpopulations. Four of these subpopulations are found exclusively in the Patagonian region of South America; one is found predominantly in Patagonia and sparsely in Oceania and North America; and one is specific to the Holarctic ecozone. Plant host associations differed between subpopulations and between S. eubayanus and its sister species, Saccharomyces uvarum. S. eubayanus is most abundant and genetically diverse in northern Patagonia, where some locations harbor more genetic diversity than is found outside of South America, suggesting that northern Patagonia east of the Andes was a glacial refugium for this species. All but one subpopulation shows isolation-by-distance, and gene flow between subpopulations is low. However, there are strong signals of ancient and recent outcrossing, including two admixed lineages, one that is sympatric with and one that is mostly isolated from its parental populations. Using our extensive biogeographical data, we build a robust model that predicts all known and a handful of additional regions of the globe that are climatically suitable for S. eubayanus, including Europe where host accessibility and competitive exclusion by other Saccharomyces species may explain its continued elusiveness. We conclude that this industrially relevant species has rich natural diversity with many factors contributing to its complex distribution and natural history.

The genetic code for amino acids is nearly universal, and among eukaryotic nuclear genomes the only known reassignments are of codon CUG in yeasts. Here, the authors identify a third independent CUG transition in budding yeasts that is still ongoing with alternative tRNAs present in the genome.

Abstract Background Associations between traits are prevalent in nature, occurring across a diverse range of taxa and traits. Individual traits may co-evolve with one other, and these correlations can be driven by factors intrinsic or extrinsic to an organism. However, few studies, especially in microbes, have simultaneously investigated both across a broad taxonomic range. Here we quantify pairwise associations among 48 traits across 784 diverse yeast species of the ancient budding yeast subphylum Saccharomycotina, assessing the effects of phylogenetic history, genetics, and ecology. Results We find extensive negative (traits that tend to not occur together) and positive (traits that tend to co-occur) pairwise associations among traits, as well as between traits and environments. These associations can largely be explained by the biological properties of the traits, such as overlapping biochemical pathways. The isolation environments of the yeasts explain a minor but significant component of the variance, while phylogeny (the retention of ancestral traits in descendant species) plays an even more limited role. Positive correlations are pervasive among carbon utilization traits and track with chemical structures (e.g., glucosides and sugar alcohols) and metabolic pathways, suggesting a molecular basis for the presence of suites of traits. In several cases, characterized genes from model organisms suggest that enzyme promiscuity and overlapping biochemical pathways are likely mechanisms to explain these macroevolutionary trends. Interestingly, fermentation traits are negatively correlated with the utilization of pentose sugars, which are major components of the plant biomass degraded by fungi and present major bottlenecks to the production of cellulosic biofuels. Finally, we show that mammalian pathogenic and commensal yeasts have a suite of traits that includes growth at high temperature and, surprisingly, the utilization of a narrowed panel of carbon sources. Conclusions These results demonstrate how both intrinsic physiological factors and extrinsic ecological factors drive the distribution of traits present in diverse organisms across macroevolutionary timescales.

Abstract Background Lignocellulosic biomass is a common resource across the globe, and its fermentation offers a promising option for generating renewable liquid transportation fuels. The deconstruction of lignocellulosic biomass releases sugars that can be fermented by microbes, but these processes also produce fermentation inhibitors, such as aromatic acids and aldehydes. Several research projects have investigated lignocellulosic biomass fermentation by the baker’s yeast Saccharomyces cerevisiae. Most projects have taken synthetic biological approaches or have explored naturally occurring diversity in S. cerevisiae to enhance stress tolerance, xylose consumption, or ethanol production. Despite these efforts, improved strains with new properties are needed. In other industrial processes, such as wine and beer fermentation, interspecies hybrids have combined important traits from multiple species, suggesting that interspecies hybridization may also offer potential for biofuel research. Results To investigate the efficacy of this approach for traits relevant to lignocellulosic biofuel production, we generated synthetic hybrids by crossing engineered xylose-fermenting strains of S. cerevisiae with wild strains from various Saccharomyces species. These interspecies hybrids retained important parental traits, such as xylose consumption and stress tolerance, while displaying intermediate kinetic parameters and, in some cases, heterosis (hybrid vigor). Next, we exposed them to adaptive evolution in ammonia fiber expansion-pretreated corn stover hydrolysate and recovered strains with improved fermentative traits. Genome sequencing showed that the genomes of these evolved synthetic hybrids underwent rearrangements, duplications, and deletions. To determine whether the genus Saccharomyces contains additional untapped potential, we screened a genetically diverse collection of more than 500 wild, non-engineered Saccharomyces isolates and uncovered a wide range of capabilities for traits relevant to cellulosic biofuel production. Notably, Saccharomyces mikatae strains have high innate tolerance to hydrolysate toxins, while some Saccharomyces species have a robust native capacity to consume xylose. Conclusions This research demonstrates that hybridization is a viable method to combine industrially relevant traits from diverse yeast species and that members of the genus Saccharomyces beyond S. cerevisiae may offer advantageous genes and traits of interest to the lignocellulosic biofuel industry.

Variation in synonymous codon usage is abundant across multiple levels of organization: between codons of an amino acid, between genes in a genome, and between genomes of different species. It is now well understood that variation in synonymous codon usage is influenced by mutational bias coupled with both natural selection for translational efficiency and genetic drift, but how these processes shape patterns of codon usage bias across entire lineages remains unexplored. To address this question, we used a rich genomic data set of 327 species that covers nearly one third of the known biodiversity of the budding yeast subphylum Saccharomycotina. We found that, while genome-wide relative synonymous codon usage (RSCU) for all codons was highly correlated with the GC content of the third codon position (GC3), the usage of codons for the amino acids proline, arginine, and glycine was inconsistent with the neutral expectation where mutational bias coupled with genetic drift drive codon usage. Examination between genes' effective numbers of codons and their GC3 contents in individual genomes revealed that nearly a quarter of genes (381,174/1,683,203; 23%), as well as most genomes (308/327; 94%), significantly deviate from the neutral expectation. Finally, by evaluating the imprint of translational selection on codon usage, measured as the degree to which genes' adaptiveness to the tRNA pool were correlated with selective pressure, we show that translational selection is widespread in budding yeast genomes (264/327; 81%). These results suggest that the contribution of translational selection and drift to patterns of synonymous codon usage across budding yeasts varies across codons, genes, and genomes; whereas drift is the primary driver of global codon usage across the subphylum, the codon bias of large numbers of genes in the majority of genomes is influenced by translational selection.

Cell-cycle checkpoints and DNA repair processes protect organisms from potentially lethal mutational damage. Compared to other budding yeasts in the subphylum Saccharomycotina, we noticed that a lineage in the genus Hanseniaspora exhibited very high evolutionary rates, low Guanine-Cytosine (GC) content, small genome sizes, and lower gene numbers. To better understand Hanseniaspora evolution, we analyzed 25 genomes, including 11 newly sequenced, representing 18/21 known species in the genus. Our phylogenomic analyses identify two Hanseniaspora lineages, a faster-evolving lineage (FEL), which began diversifying approximately 87 million years ago (mya), and a slower-evolving lineage (SEL), which began diversifying approximately 54 mya. Remarkably, both lineages lost genes associated with the cell cycle and genome integrity, but these losses were greater in the FEL. E.g., all species lost the cell-cycle regulator WHIskey 5 (WHI5), and the FEL lost components of the spindle checkpoint pathway (e.g., Mitotic Arrest-Deficient 1 [MAD1], Mitotic Arrest-Deficient 2 [MAD2]) and DNA-damage-checkpoint pathway (e.g., Mitosis Entry Checkpoint 3 [MEC3], RADiation sensitive 9 [RAD9]). Similarly, both lineages lost genes involved in DNA repair pathways, including the DNA glycosylase gene 3-MethylAdenine DNA Glycosylase 1 (MAG1), which is part of the base-excision repair pathway, and the DNA photolyase gene PHotoreactivation Repair deficient 1 (PHR1), which is involved in pyrimidine dimer repair. Strikingly, the FEL lost 33 additional genes, including polymerases (i.e., POLymerase 4 [POL4] and POL32) and telomere-associated genes (e.g., Repressor/activator site binding protein-Interacting Factor 1 [RIF1], Replication Factor A 3 [RFA3], Cell Division Cycle 13 [CDC13], Pbp1p Binding Protein [PBP2]). Echoing these losses, molecular evolutionary analyses reveal that, compared to the SEL, the FEL stem lineage underwent a burst of accelerated evolution, which resulted in greater mutational loads, homopolymer instabilities, and higher fractions of mutations associated with the common endogenously damaged base, 8-oxoguanine. We conclude that Hanseniaspora is an ancient lineage that has diversified and thrived, despite lacking many otherwise highly conserved cell-cycle and genome integrity genes and pathways, and may represent a novel, to our knowledge, system for studying cellular life without them.

At the molecular level, the evolution of new traits can be broadly divided between changes in gene expression and changes in protein-coding sequence. For proteins, the evolution of novel functions is generally thought to proceed through sequential point mutations or recombination of whole functional units. In Saccharomyces, the uptake of the sugar maltotriose into the cell is the primary limiting factor in its utilization, but maltotriose transporters are relatively rare, except in brewing strains. No known wild strains of Saccharomyces eubayanus, the cold-tolerant parent of hybrid lager-brewing yeasts (Saccharomyces cerevisiae x S. eubayanus), are able to consume maltotriose, which limits their ability to fully ferment malt extract. In one strain of S. eubayanus, we found a gene closely related to a known maltotriose transporter and were able to confer maltotriose consumption by overexpressing this gene or by passaging the strain on maltose. Even so, most wild strains of S. eubayanus lack native maltotriose transporters. To determine how this rare trait could evolve in naive genetic backgrounds, we performed an adaptive evolution experiment for maltotriose consumption, which yielded a single strain of S. eubayanus able to grow on maltotriose. We mapped the causative locus to a gene encoding a novel chimeric transporter that was formed by an ectopic recombination event between two genes encoding transporters that are unable to import maltotriose. In contrast to classic models of the evolution of novel protein functions, the recombination breakpoints occurred within a single functional domain. Thus, the ability of the new protein to carry maltotriose was likely acquired through epistatic interactions between independently evolved substitutions. By acquiring multiple mutations at once, the transporter rapidly gained a novel function, while bypassing potentially deleterious intermediate steps. This study provides an illuminating example of how recombination between paralogs can establish novel interactions among substitutions to create adaptive functions.

A novel yeast species was isolated from the sugar-rich stromata of Cyttaria hariotii collected from two different Nothofagus tree species in the Andean forests of Patagonia, Argentina. Phylogenetic analyses of the concatenated sequence of the rRNA gene sequences and the protein-coding genes for actin and translational elongation factor-1α indicated that the novel species belongs to the genus Hanseniaspora. De novo genome assembly of the strain CRUB 1928T yielded a 10.2-Mbp genome assembly predicted to encode 4452 protein-coding genes. The genome sequence data were compared to the genomes of other Hanseniaspora species using three different methods, an alignment-free distance measure, Kr, and two model-based estimations of DNA-DNA homology values, of which all provided indicative values to delineate species of Hanseniaspora. Given its potential role in a rare indigenous alcoholic beverage in which yeasts ferment sugars extracted from the stromata of Cytarria sp., we searched for the genes that may suggest adaptation of novel Hanseniaspora species to fermenting communities. The SSU1-like gene encoding a sulfite efflux pump, which, among Hanseniaspora, is present only in close relatives to the new species, was detected and analyzed, suggesting that this gene might be one factor that characterizes this novel species. We also discuss several candidate genes that likely underlie the physiological traits used for traditional taxonomic identification. Based on these results, a novel yeast species with the name Hanseniaspora gamundiae sp. nov. is proposed with CRUB 1928T (ex-types: ZIM 2545T = NRRL Y-63793T = PYCC 7262T; MycoBank number MB 824091) as the type strain. Furthermore, we propose the transfer of the Kloeckera species, K. hatyaiensis, K. lindneri and K. taiwanica to the genus Hanseniaspora as Hanseniaspora hatyaiensis comb. nov. (MB 828569), Hanseniaspora lindneri comb. nov. (MB 828566) and Hanseniaspora taiwanica comb. nov. (MB 828567).

Understanding the phylogenetic relationships among the yeasts of the subphylum Saccharomycotina is a prerequisite for understanding the evolution of their metabolisms and ecological lifestyles. In the last two decades, the use of rDNA and multilocus data sets has greatly advanced our understanding of the yeast phylogeny, but many deep relationships remain unsupported. In contrast, phylogenomic analyses have involved relatively few taxa and lineages that were often selected with limited considerations for covering the breadth of yeast biodiversity. Here we used genome sequence data from 86 publicly available yeast genomes representing nine of the 11 known major lineages and 10 nonyeast fungal outgroups to generate a 1233-gene, 96-taxon data matrix. Species phylogenies reconstructed using two different methods (concatenation and coalescence) and two data matrices (amino acids or the first two codon positions) yielded identical and highly supported relationships between the nine major lineages. Aside from the lineage comprised by the family Pichiaceae, all other lineages were monophyletic. Most interrelationships among yeast species were robust across the two methods and data matrices. However, eight of the 93 internodes conflicted between analyses or data sets, including the placements of: the clade defined by species that have reassigned the CUG codon to encode serine, instead of leucine; the clade defined by a whole genome duplication; and the species Ascoidea rubescens. These phylogenomic analyses provide a robust roadmap for future comparative work across the yeast subphylum in the disciplines of taxonomy, molecular genetics, evolutionary biology, ecology, and biotechnology. To further this end, we have also provided a BLAST server to query the 86 Saccharomycotina genomes, which can be found at http://y1000plus.org/blast.

The availability of genomes across the tree of life is highly biased toward vertebrates, pathogens, human disease models, and organisms with relatively small and simple genomes. Recent progress in genomics has enabled the de novo decoding of the genome of virtually any organism, greatly expanding its potential for understanding the biology and evolution of the full spectrum of biodiversity. The increasing diversity of sequencing technologies, assays, and de novo assembly algorithms have augmented the complexity of de novo genome sequencing projects in nonmodel organisms. To reduce the costs and challenges in de novo genome sequencing projects and streamline their experimental design and analysis, we developed iWGS (in silico Whole Genome Sequencer and Analyzer), an automated pipeline for guiding the choice of appropriate sequencing strategy and assembly protocols. iWGS seamlessly integrates the four key steps of a de novo genome sequencing project: data generation (through simulation), data quality control, de novo assembly, and assembly evaluation and validation. The last three steps can also be applied to the analysis of real data. iWGS is designed to enable the user to have great flexibility in testing the range of experimental designs available for genome sequencing projects, and supports all major sequencing technologies and popular assembly tools. Three case studies illustrate how iWGS can guide the design of de novo genome sequencing projects, and evaluate the performance of a wide variety of user-specified sequencing strategies and assembly protocols on genomes of differing architectures. iWGS, along with a detailed documentation, is freely available at https://github.com/zhouxiaofan1983/iWGS.

The genome sequences of more than 100 strains of the yeast Saccharomyces cerevisiae have been published. Unfortunately, most of these genome assemblies contain dozens to hundreds of gaps at repetitive sequences, including transposable elements, tRNAs, and subtelomeric regions, which is where novel genes generally reside. Relatively few strains have been chosen for genome sequencing based on their biofuel production potential, leaving an additional knowledge gap. Here, we describe the nearly complete genome sequence of GLBRCY22-3 (Y22-3), a strain of S. cerevisiae derived from the stress-tolerant wild strain NRRL YB-210 and subsequently engineered for xylose metabolism. After benchmarking several genome assembly approaches, we developed a pipeline to integrate Pacific Biosciences (PacBio) and Illumina sequencing data and achieved one of the highest quality genome assemblies for any S. cerevisiae strain. Specifically, the contig N50 is 693 kbp, and the sequences of most chromosomes, the mitochondrial genome, and the 2-micron plasmid are complete. Our annotation predicts 92 genes that are not present in the reference genome of the laboratory strain S288c, over 70% of which were expressed. We predicted functions for 43 of these genes, 28 of which were previously uncharacterized and unnamed. Remarkably, many of these genes are predicted to be involved in stress tolerance and carbon metabolism and are shared with a Brazilian bioethanol production strain, even though the strains differ dramatically at most genetic loci. The Y22-3 genome sequence provides an exceptionally high-quality resource for basic and applied research in bioenergy and genetics.

The pollen stores of bumble bees host diverse microbiota that influence overall colony fitness. Yet, the taxonomic identity of these symbiotic microbes is relatively unknown. In this descriptive study, we characterized the microbial community of pollen provisions within captive-bred bumble bee hives obtained from two commercial suppliers located in North America. Findings from 16S rRNA and ITS gene-based analyses revealed that pollen provisions from the captive-bred hives shared several microbial taxa that have been previously detected among wild populations. While diverse microbes across phyla Firmicutes, Proteobacteria, Bacteroidetes, Actinobacteria, and Ascomycota were detected in all commercial hives, significant differences were detected at finer-scale taxonomic resolution based on the supplier source. The causative agent of chalkbrood disease in honey bees, Ascosphaera apis, was detected in all hives obtained from one supplier source, although none of the hives showed symptoms of infection. The shared core microbiota across both commercial supplier sources consisted of two ubiquitous bee-associated groups, Lactobacillus and Wickerhamiella/Starmerella clade yeasts that potentially contribute to the beneficial function of the microbiome of bumble bee pollen provisions.

Fructophily is a rare trait that consists of the preference for fructose over other carbon sources. Here, we show that in a yeast lineage (the Wickerhamiella/Starmerella, W/S clade) comprised of fructophilic species thriving in the high-sugar floral niche, the acquisition of fructophily is concurrent with a wider remodeling of central carbon metabolism. Coupling comparative genomics with biochemical and genetic approaches, we gathered ample evidence for the loss of alcoholic fermentation in an ancestor of the W/S clade and subsequent reinstatement through either horizontal acquisition of homologous bacterial genes or modification of a pre-existing yeast gene. An enzyme required for sucrose assimilation was also acquired from bacteria, suggesting that the genetic novelties identified in the W/S clade may be related to adaptation to the high-sugar environment. This work shows how even central carbon metabolism can be remodeled by a surge of HGT events.

[This corrects the article DOI: 10.1371/journal.pgen.1006372.].

The inability of native Saccharomyces cerevisiae to convert xylose from plant biomass into biofuels remains a major challenge for the production of renewable bioenergy. Despite extensive knowledge of the regulatory networks controlling carbon metabolism in yeast, little is known about how to reprogram S. cerevisiae to ferment xylose at rates comparable to glucose. Here we combined genome sequencing, proteomic profiling, and metabolomic analyses to identify and characterize the responsible mutations in a series of evolved strains capable of metabolizing xylose aerobically or anaerobically. We report that rapid xylose conversion by engineered and evolved S. cerevisiae strains depends upon epistatic interactions among genes encoding a xylose reductase (GRE3), a component of MAP Kinase (MAPK) signaling (HOG1), a regulator of Protein Kinase A (PKA) signaling (IRA2), and a scaffolding protein for mitochondrial iron-sulfur (Fe-S) cluster biogenesis (ISU1). Interestingly, the mutation in IRA2 only impacted anaerobic xylose consumption and required the loss of ISU1 function, indicating a previously unknown connection between PKA signaling, Fe-S cluster biogenesis, and anaerobiosis. Proteomic and metabolomic comparisons revealed that the xylose-metabolizing mutant strains exhibit altered metabolic pathways relative to the parental strain when grown in xylose. Further analyses revealed that interacting mutations in HOG1 and ISU1 unexpectedly elevated mitochondrial respiratory proteins and enabled rapid aerobic respiration of xylose and other non-fermentable carbon substrates. Our findings suggest a surprising connection between Fe-S cluster biogenesis and signaling that facilitates aerobic respiration and anaerobic fermentation of xylose, underscoring how much remains unknown about the eukaryotic signaling systems that regulate carbon metabolism.

The evolutionary mechanisms leading to duplicate gene retention are well understood, but the long-term impacts of paralog differentiation on the regulation of metabolism remain underappreciated. Here we experimentally dissect the functions of two pairs of ancient paralogs of the GALactose sugar utilization network in two yeast species. We show that the Saccharomyces uvarum network is more active, even as over-induction is prevented by a second co-repressor that the model yeast Saccharomyces cerevisiae lacks. Surprisingly, removal of this repression system leads to a strong growth arrest, likely due to overly rapid galactose catabolism and metabolic overload. Alternative sugars, such as fructose, circumvent metabolic control systems and exacerbate this phenotype. We further show that S. cerevisiae experiences homologous metabolic constraints that are subtler due to how the paralogs have diversified. These results show how the functional differentiation of paralogs continues to shape regulatory network architectures and metabolic strategies long after initial preservation.

Lager-style beers constitute the vast majority of the beer market, and yet, the genetic origin of the yeast strains that brew them has been shrouded in mystery and controversy. Unlike ale-style beers, which are generally brewed with Saccharomyces cerevisiae, lagers are brewed at colder temperatures with allopolyploid hybrids of Saccharomyces eubayanus x S. cerevisiae. Since the discovery of S. eubayanus in 2011, additional strains have been isolated from South America, North America, Australasia, and Asia, but only interspecies hybrids have been isolated in Europe. Here, using genome sequence data, we examine the relationships of these wild S. eubayanus strains to each other and to domesticated lager strains. Our results support the existence of a relatively low-diversity (π = 0.00197) lineage of S. eubayanus whose distribution stretches across the Holarctic ecozone and includes wild isolates from Tibet, new wild isolates from North America, and the S. eubayanus parents of lager yeasts. This Holarctic lineage is closely related to a population with higher diversity (π = 0.00275) that has been found primarily in South America but includes some widely distributed isolates. A second diverse South American population (π = 0.00354) and two early-diverging Asian subspecies are more distantly related. We further show that no single wild strain from the Holarctic lineage is the sole closest relative of lager yeasts. Instead, different parts of the genome portray different phylogenetic signals and ancestry, likely due to outcrossing and incomplete lineage sorting. Indeed, standing genetic variation within this wild Holarctic lineage of S. eubayanus is responsible for genetic variation still segregating among modern lager-brewing hybrids. We conclude that the relationships among wild strains of S. eubayanus and their domesticated hybrids reflect complex biogeographical and genetic processes.