Comparative genomics of protoploid Saccharomycetaceae

Quantitative phenotypic analysis of multistress response in Zygosaccharomyces rouxii complex

Zygosaccharomyces rouxii complex comprises three yeasts clusters sourced from sugar- and salt-rich environments: haploid Zygosaccharomyces rouxii, diploid Zygosaccharomyces sapae and allodiploid/aneuploid strains of uncertain taxonomic affiliations. These yeasts have been characterized with respect to gene copy number variation, karyotype variability and change in ploidy, but functional diversity in stress responses has not been explored yet. Here, we quantitatively analysed the stress response variation in seven strains of the Z. rouxii complex by modelling growth variables via model and model-free fitting methods. Based on the spline fit as most reliable modelling method, we resolved different interstrain responses to 15 environmental perturbations. Compared with Z. rouxii CBS 732(T) and Z. sapae strains ABT301(T) and ABT601, allodiploid strain ATCC 42981 and aneuploid strains CBS 4837 and CBS 4838 displayed higher multistress resistance and better performance in glycerol respiration even in the presence of copper. μ-based logarithmic phenotypic index highlighted that ABT601 is a slow-growing strain insensitive to stress, whereas ABT301(T) grows fast on rich medium and is sensitive to suboptimal conditions. Overall, the differences in stress response could imply different adaptation mechanisms to sugar- and salt-rich niches. The obtained phenotypic profiling contributes to provide quantitative insights for elucidating the adaptive mechanisms to stress in halo- and osmo-tolerant Zygosaccharomyces yeasts.

The genome sequence of the highly acetic acid-tolerant Zygosaccharomyces bailii-derived interspecies hybrid strain ISA1307, isolated from a sparkling wine plant

In this work, it is described the sequencing and annotation of the genome of the yeast strain ISA1307, isolated from a sparkling wine continuous production plant. This strain, formerly considered of the Zygosaccharomyces bailii species, has been used to study Z. bailii physiology, in particular, its extreme tolerance to acetic acid stress at low pH. The analysis of the genome sequence described in this work indicates that strain ISA1307 is an interspecies hybrid between Z. bailii and a closely related species. The genome sequence of ISA1307 is distributed through 154 scaffolds and has a size of around 21.2 Mb, corresponding to 96% of the genome size estimated by flow cytometry. Annotation of ISA1307 genome includes 4385 duplicated genes (∼90% of the total number of predicted genes) and 1155 predicted single-copy genes. The functional categories including a higher number of genes are ‘Metabolism and generation of energy’, ‘Protein folding, modification and targeting’ and ‘Biogenesis of cellular components’. The knowledge of the genome sequence of the ISA1307 strain is expected to contribute to accelerate systems-level understanding of stress resistance mechanisms in Z. bailii and to inspire and guide novel biotechnological applications of this yeast species/strain in fermentation processes, given its high resilience to acidic stress. The availability of the ISA1307 genome sequence also paves the way to a better understanding of the genetic mechanisms underlying the generation and selection of more robust hybrid yeast strains in the stressful environment of wine fermentations.

The complete genome of Blastobotrys (Arxula) adeninivorans LS3 - a yeast of biotechnological interest

BACKGROUND

The industrially important yeast Blastobotrys (Arxula) adeninivorans is an asexual hemiascomycete phylogenetically very distant from Saccharomyces cerevisiae. Its unusual metabolic flexibility allows it to use a wide range of carbon and nitrogen sources, while being thermotolerant, xerotolerant and osmotolerant.

RESULTS

The sequencing of strain LS3 revealed that the nuclear genome of A. adeninivorans is 11.8 Mb long and consists of four chromosomes with regional centromeres. Its closest sequenced relative is Yarrowia lipolytica, although mean conservation of orthologs is low. With 914 introns within 6116 genes, A. adeninivorans is one of the most intron-rich hemiascomycetes sequenced to date. Several large species-specific families appear to result from multiple rounds of segmental duplications of tandem gene arrays, a novel mechanism not yet described in yeasts. An analysis of the genome and its transcriptome revealed enzymes with biotechnological potential, such as two extracellular tannases (Atan1p and Atan2p) of the tannic-acid catabolic route, and a new pathway for the assimilation of n-butanol via butyric aldehyde and butyric acid.

CONCLUSIONS

The high-quality genome of this species that diverged early in Saccharomycotina will allow further fundamental studies on comparative genomics, evolution and phylogenetics. Protein components of different pathways for carbon and nitrogen source utilization were identified, which so far has remained unexplored in yeast, offering clues for further biotechnological developments. In the course of identifying alternative microorganisms for biotechnological interest, A. adeninivorans has already proved its strengthened competitiveness as a promising cell factory for many more applications.

Rapid identification of Zygosaccharomyces with genus-specific primers

There has been a recent and rapid increase in the number of species of the genus Zygosaccharomyces which now comprises Z. bailii, Z. bisporus, Z. gambellarensis, Z. kombuchaensis, Z. lentus, Z. machadoi, Z. mellis, Z. parabaillii, Z. pseudobailii, Z. pseudorouxii, Z. rouxii, Z. sapae, and Z. siamensis. Z. pseudorouxii is an unofficial name given to isolates closely related to the newly-described species Z. sapae. The Zygosaccharomyces genus contains species that are important as food and beverage spoilage organisms and others are associated with fermentations and sweet foodstuffs, such as honey. Their economic significance means that the ability to identify them rapidly is of significant importance. Although Z. rouxii and Z. bailii have been genome-sequenced the extent of sequence data for the others, especially the newly-discovered species, is sometimes extremely limited which makes identification slow. However, parts of the ITS1/5.8S/ITS2 rDNA region contain sequences of sufficient similarity within the genus and of sufficient difference with outgroups, to be potential regions for the design of genus-wide specific primers. We report here the development of genus-specific primers that can detect all the major Zygosaccharomyces species including all those associated with foods; the rare and localised species Z. machadoi and Z. gambellarensis are not detected. The size of the single amplicon produced varies between species and in some cases is sufficiently different to assign provisional species identification. Sequence data from rDNA regions are available for virtually all described yeast species in all genera, thus, prior to having sufficient sequence data from structural genes, rDNA regions may provide more generally suitable candidates for both genus-specific and species-specific primer design.

Genome sequence and analysis of methylotrophic yeast Hansenula polymorpha DL1

BACKGROUND

Hansenula polymorpha DL1 is a methylotrophic yeast, widely used in fundamental studies of methanol metabolism, peroxisome biogenesis and function, and also as a microbial cell factory for production of recombinant proteins and metabolic engineering towards the goal of high temperature ethanol production.

RESULTS

We have sequenced the 9 Mbp H. polymorpha DL1 genome and performed whole-genome analysis for the H. polymorpha transcriptome obtained from both methanol- and glucose-grown cells. RNA-seq analysis revealed the complex and dynamic character of the H. polymorpha transcriptome under the two studied conditions, identified abundant and highly unregulated expression of 40% of the genome in methanol grown cells, and revealed alternative splicing events. We have identified subtelomerically biased protein families in H. polymorpha, clusters of LTR elements at G + C-poor chromosomal loci in the middle of each of the seven H. polymorpha chromosomes, and established the evolutionary position of H. polymorpha DL1 within a separate yeast clade together with the methylotrophic yeast Pichia pastoris and the non-methylotrophic yeast Dekkera bruxellensis. Intergenome comparisons uncovered extensive gene order reshuffling between the three yeast genomes. Phylogenetic analyses enabled us to reveal patterns of evolution of methylotrophy in yeasts and filamentous fungi.

CONCLUSIONS

Our results open new opportunities for in-depth understanding of many aspects of H. polymorpha life cycle, physiology and metabolism as well as genome evolution in methylotrophic yeasts and may lead to novel improvements toward the application of H. polymorpha DL-1 as a microbial cell factory.

Global expression analysis of the yeast Lachancea (Saccharomyces) kluyveri reveals new URC genes involved in pyrimidine catabolism

Pyrimidines are important nucleic acid precursors which are constantly synthesized, degraded, and rebuilt in the cell. Four degradation pathways, two of which are found in eukaryotes, have been described. One of them, the URC pathway, has been initially discovered in our laboratory in the yeast Lachancea kluyveri. Here, we present the global changes in gene expression in L. kluyveri in response to different nitrogen sources, including uracil, uridine, dihydrouracil, and ammonia. The expression pattern of the known URC genes, URC1-6, helped to identify nine putative novel URC genes with a similar expression pattern. The microarray analysis provided evidence that both the URC and PYD genes are under nitrogen catabolite repression in L. kluyveri and are induced by uracil or dihydrouracil, respectively. We determined the function of URC8, which was found to catalyze the reduction of malonate semialdehyde to 3-hydroxypropionate, the final degradation product of the pathway. The other eight genes studied were all putative permeases. Our analysis of double deletion strains showed that the L. kluyveri Fui1p protein transported uridine, just like its homolog in Saccharomyces cerevisiae, but we demonstrated that is was not the only uridine transporter in L. kluyveri. We also showed that the L. kluyveri homologs of DUR3 and FUR4 do not have the same function that they have in S. cerevisiae, where they transport urea and uracil, respectively. In L. kluyveri, both of these deletion strains grew normally on uracil and urea.

Evolutionary mobility of the ribosomal DNA array in yeasts

The ribosomal DNA (rDNA) of eukaryotes is organized as large tandem arrays. Here, we compare the genomic locations of rDNA among yeast species and show that, despite its huge size (>1 Mb), the rDNA array has moved around the genome several times within the family Saccharomycetaceae. We identify an ancestral, nontelomeric, rDNA site that is conserved across many species including Saccharomyces cerevisiae. Within the genus Lachancea, however, the rDNA apparently transposed from the ancestral site to a new site internal to a different chromosome, becoming inserted into a short intergenic region beside a tRNA gene. In at least four other yeast lineages, the rDNA moved from the ancestral site to telomeric locations. Remarkably, both the ancestral rDNA site and the new site in Lachancea are adjacent to protein-coding genes whose products maintain the specialized chromatin structure of rDNA (HMO1 and CDC14, respectively). In almost every case where the rDNA was lost from the ancestral site, the entire array disappeared without any other rearrangements in the region, leaving just an intergenic spacer of less than 2 kb. The mechanism by which this large and complex locus moves around the genome is unknown, but we speculate that it may involve the formation of double-strand DNA breaks by Fob1 protein or the formation of extrachromosomal rDNA circles.

Genomes of Ashbya fungi isolated from insects reveal four mating-type loci, numerous translocations, lack of transposons, and distinct gene duplications

The filamentous fungus Ashbya gossypii is a cotton pathogen transmitted by insects. It is readily grown and manipulated in the laboratory and is commercially exploited as a natural overproducer of vitamin B2. Our previous genome analysis of A. gossypii isolate ATCC10895, collected in Trinidad nearly 100 years ago, revealed extensive synteny with the Saccharomyces cerevisiae genome, leading us to use it as a model organism to understand the evolution of filamentous growth. To further develop Ashbya as a model system, we have investigated the ecological niche of A. gossypii and isolated additional strains and a sibling species, both useful in comparative analysis. We isolated fungi morphologically similar to A. gossypii from different plant-feeding insects of the suborder Heteroptera, generated a phylogenetic tree based on rDNA-ITS sequences, and performed high coverage short read sequencing with one A. gossypii isolate from Florida, a new species, Ashbya aceri, isolated in North Carolina, and a genetically marked derivative of ATCC10895 intensively used for functional studies. In contrast to S. cerevisiae, all strains carry four not three mating type loci, adding a new puzzle in the evolution of Ashbya species. Another surprise was the genome identity of 99.9% between the Florida strain and ATCC10895, isolated in Trinidad. The A. aceri and A. gossypii genomes show conserved gene orders rearranged by eight translocations, 90% overall sequence identity, and fewer tandem duplications in the A. aceri genome. Both species lack transposable elements. Finally, our work identifies plant-feeding insects of the suborder Heteroptera as the most likely natural reservoir of Ashbya, and that infection of cotton and other plants may be incidental to the growth of the fungus in its insect host.

Diversity of mating-type chromosome structures in the yeast Zygosaccharomyces rouxii caused by ectopic exchanges between MAT-like loci

We investigated sex chromosome diversity in Zygosaccharomyces rouxii (Z. rouxii). In the current study, we show that the organization of the mating-type (MAT) locus is highly variable in the Z. rouxii population, indicating the MAT, HML, and HMR loci are translocation hotspots. Although NBRC1130 and CBS732 were originally two stocks of the type strain of the species, only NBRC1130 retains the original karyotype. A reciprocal translocation between the MAT and HMR loci appears to have occurred during the early passage culture of CBS732, which was used for genome sequencing. In NBRC1733, NBRC0686, NBRC0740 and NBRC1053, the terminal region of the chromosome containing the HMR locus was replaced with the chromosomal region to the left of the MAT or HML loci. The translocation events found in NBRC1733, NBRC0686, NBRC0740, and NBRC1053 were reconstructed under our experimental conditions using the DA2 background, and the reconstruction suggests that the frequency of this type of translocation is approximately 10(-7). These results suggest that the MAT and MAT-like loci were the susceptible regions in the genome, and the diversity of mating-type chromosome structures in Z. rouxii was caused by ectopic exchanges between MAT-like loci.

Detection and characterization of megasatellites in orthologous and nonorthologous genes of 21 fungal genomes

Megasatellites are large DNA tandem repeats, originally described in Candida glabrata, in protein-coding genes. Most of the genes in which megasatellites are found are of unknown function. In this work, we extended the search for megasatellites to 20 additional completely sequenced fungal genomes and extracted 216 megasatellites in 203 out of 142,121 genes, corresponding to the most exhaustive description of such genetic elements available today. We show that half of the megasatellites detected encode threonine-rich peptides predicted to be intrinsically disordered, suggesting that they may interact with several partners or serve as flexible linkers. Megasatellite motifs were clustered into several families. Their distribution in fungal genes shows that different motifs are found in orthologous genes and similar motifs are found in unrelated genes, suggesting that megasatellite formation or spreading does not necessarily track the evolution of their host genes. Altogether, these results suggest that megasatellites are created and lost during evolution of fungal genomes, probably sharing similar functions, although their primary sequences are not necessarily conserved.

Hsp90 regulates nongenetic variation in response to environmental stress

Nongenetic cell-to-cell variability often plays an important role for the survival of a clonal population in the face of fluctuating environments. However, the underlying mechanisms regulating such nongenetic heterogeneity remain elusive in most organisms. We report here that a clonal yeast population exhibits morphological heterogeneity when the level of Hsp90, a molecular chaperone, is reduced. The morphological heterogeneity is driven by the dosage of Cdc28 and Cla4, a key regulator of septin formation. Low Hsp90 levels reduce Cla4 protein stability and cause a subpopulation of cells to switch to a filamentous form that has been previously suggested to be beneficial under certain hostile environments. Moreover, Hsp90-dependent morphological heterogeneity can be induced by environmental stress and is conserved across diverse yeast species. Our results suggest that Hsp90 provides an evolutionarily conserved mechanism that links environmental stress to the induction of morphological diversity.

The spatiotemporal program of replication in the genome of Lachancea kluyveri

We generated a genome-wide replication profile in the genome of Lachancea kluyveri and assessed the relationship between replication and base composition. This species diverged from Saccharomyces cerevisiae before the ancestral whole genome duplication. The genome comprises eight chromosomes among which a chromosomal arm of 1 Mb has a G + C-content much higher than the rest of the genome. We identified 252 active replication origins in L. kluyveri and found considerable divergence in origin location with S. cerevisiae and with Lachancea waltii. Although some global features of S. cerevisiae replication are conserved: Centromeres replicate early, whereas telomeres replicate late, we found that replication origins both in L. kluyveri and L. waltii do not behave as evolutionary fragile sites. In L. kluyveri, replication timing along chromosomes alternates between regions of early and late activating origins, except for the 1 Mb GC-rich chromosomal arm. This chromosomal arm contains an origin consensus motif different from other chromosomes and is replicated early during S-phase. We showed that precocious replication results from the specific absence of late firing origins in this chromosomal arm. In addition, we found a correlation between GC-content and distance from replication origins as well as a lack of replication-associated compositional skew between leading and lagging strands specifically in this GC-rich chromosomal arm. These findings suggest that the unusual base composition in the genome of L. kluyveri could be linked to replication.

Evolutionary role of interspecies hybridization and genetic exchanges in yeasts

Forced interspecific hybridization has been used in yeasts for many years to study speciation or to construct artificial strains with novel fermentative and metabolic properties. Recent genome analyses indicate that natural hybrids are also generated spontaneously between yeasts belonging to distinct species, creating lineages with novel phenotypes, varied genetic stability, or altered virulence in the case of pathogens. Large segmental introgressions from evolutionarily distant species are also visible in some yeast genomes, suggesting that interspecific genetic exchanges occur during evolution. The origin of this phenomenon remains unclear, but it is likely based on weak prezygotic barriers, limited Dobzhansky-Muller (DM) incompatibilities, and rapid clonal expansions. Newly formed interspecies hybrids suffer rapid changes in the genetic contribution of each parent, including chromosome loss or aneuploidy, translocations, and loss of heterozygosity, that, except in a few recently studied cases, remain to be characterized more precisely at the genomic level by use of modern technologies. We review here known cases of natural or artificially formed interspecies hybrids between yeasts and discuss their potential importance in terms of genome evolution. Problems of meiotic fertility, ploidy constraint, gene and gene product compatibility, and nucleomitochondrial interactions are discussed and placed in the context of other known mechanisms of yeast genome evolution as a model for eukaryotes.

Comparative mitochondrial genomics within and among yeast species of the Lachancea genus

Yeasts are leading model organisms for mitochondrial genome studies. The explosion of complete sequence of yeast mitochondrial (mt) genomes revealed a wide diversity of organization and structure between species. Recently, genome-wide polymorphism survey on the mt genome of isolates of a single species, Lachancea kluyveri, was also performed. To compare the mitochondrial genome evolution at two hierarchical levels: within and among closely related species, we focused on five species of the Lachancea genus, which are close relatives of L. kluyveri. Hence, we sequenced the complete mt genome of L. dasiensis, L. nothofagi, L. mirantina, L. fantastica and L. meyersii. The phylogeny of the Lachancea genus was explored using these data. Analysis of intra- and interspecific variability across the whole Lachancea genus led to the same conclusions regarding the mitochondrial genome evolution. These genomes exhibit a similar architecture and are completely syntenic. Nevertheless, genome sizes vary considerably because of the variations of the intergenic regions and the intron content, contributing to mitochondrial genome plasticity. The high variability of the intergenic regions stands in contrast to the high level of similarity of protein sequences. Quantification of the selective constraints clearly revealed that most of the mitochondrial genes are under purifying selection in the whole genus.

Genome-wide metabolic (re-) annotation of Kluyveromyces lactis

BACKGROUND

Even before having its genome sequence published in 2004, Kluyveromyces lactis had long been considered a model organism for studies in genetics and physiology. Research on Kluyveromyces lactis is quite advanced and this yeast species is one of the few with which it is possible to perform formal genetic analysis. Nevertheless, until now, no complete metabolic functional annotation has been performed to the proteins encoded in the Kluyveromyces lactis genome.

RESULTS

In this work, a new metabolic genome-wide functional re-annotation of the proteins encoded in the Kluyveromyces lactis genome was performed, resulting in the annotation of 1759 genes with metabolic functions, and the development of a methodology supported by merlin (software developed in-house). The new annotation includes novelties, such as the assignment of transporter superfamily numbers to genes identified as transporter proteins. Thus, the genes annotated with metabolic functions could be exclusively enzymatic (1410 genes), transporter proteins encoding genes (301 genes) or have both metabolic activities (48 genes). The new annotation produced by this work largely surpassed the Kluyveromyces lactis currently available annotations. A comparison with KEGG's annotation revealed a match with 844 (~90%) of the genes annotated by KEGG, while adding 850 new gene annotations. Moreover, there are 32 genes with annotations different from KEGG.

CONCLUSIONS

The methodology developed throughout this work can be used to re-annotate any yeast or, with a little tweak of the reference organism, the proteins encoded in any sequenced genome. The new annotation provided by this study offers basic knowledge which might be useful for the scientific community working on this model yeast, because new functions have been identified for the so-called metabolic genes. Furthermore, it served as the basis for the reconstruction of a compartmentalized, genome-scale metabolic model of Kluyveromyces lactis, which is currently being finished.

A pipeline for automated annotation of yeast genome sequences by a conserved-synteny approach

BACKGROUND

Yeasts are a model system for exploring eukaryotic genome evolution. Next-generation sequencing technologies are poised to vastly increase the number of yeast genome sequences, both from resequencing projects (population studies) and from de novo sequencing projects (new species). However, the annotation of genomes presents a major bottleneck for de novo projects, because it still relies on a process that is largely manual.

RESULTS

Here we present the Yeast Genome Annotation Pipeline (YGAP), an automated system designed specifically for new yeast genome sequences lacking transcriptome data. YGAP does automatic de novo annotation, exploiting homology and synteny information from other yeast species stored in the Yeast Gene Order Browser (YGOB) database. The basic premises underlying YGAP's approach are that data from other species already tells us what genes we should expect to find in any particular genomic region and that we should also expect that orthologous genes are likely to have similar intron/exon structures. Additionally, it is able to detect probable frameshift sequencing errors and can propose corrections for them. YGAP searches intelligently for introns, and detects tRNA genes and Ty-like elements.

CONCLUSIONS

In tests on Saccharomyces cerevisiae and on the genomes of Naumovozyma castellii and Tetrapisispora blattae newly sequenced with Roche-454 technology, YGAP outperformed another popular annotation program (AUGUSTUS). For S. cerevisiae and N. castellii, 91-93% of YGAP's predicted gene structures were identical to those in previous manually curated gene sets. YGAP has been implemented as a webserver with a user-friendly interface at http://wolfe.gen.tcd.ie/annotation.

Mitochondrial genome evolution in a single protoploid yeast species

Mitochondria are organelles, which play a key role in some essential functions, including respiration, metabolite biosynthesis, ion homeostasis, and apoptosis. The vast numbers of mitochondrial DNA (mtDNA) sequences of various yeast species, which have recently been published, have also helped to elucidate the structural diversity of these genomes. Although a large corpus of data are now available on the diversity of yeast species, little is known so far about the mtDNA diversity in single yeast species. To study the genetic variations occurring in the mtDNA of wild yeast isolates, we performed a genome-wide polymorphism survey on the mtDNA of 18 Lachancea kluyveri (formerly Saccharomyces kluyveri) strains. We determined the complete mt genome sequences of strains isolated from various geographical locations (in North America, Asia, and Europe) and ecological niches (Drosophila, tree exudates, soil). The mt genome of the NCYC 543 reference strain is 51,525 bp long. It contains the same core of genes as Lachancea thermotolerans, the nearest relative to L. kluyveri. To explore the mt genome variations in a single yeast species, we compared the mtDNAs of the 18 isolates. The phylogeny and population structure of L. kluyveri provide clear-cut evidence for the existence of well-defined geographically isolated lineages. Although these genomes are completely syntenic, their size and the intron content were found to vary among the isolates studied. These genomes are highly polymorphic, showing an average diversity of 28.5 SNPs/kb and 6.6 indels/kb. Analysis of the SNP and indel patterns showed the existence of a particularly high overall level of polymorphism in the intergenic regions. The dN/dS ratios obtained are consistent with purifying selection in all these genes, with the noteworthy exception of the VAR1 gene, which gave a very high ratio. These data suggest that the intergenic regions have evolved very fast in yeast mitochondrial genomes.

Sos7, an essential component of the conserved Schizosaccharomyces pombe Ndc80-MIND-Spc7 complex, identifies a new family of fungal kinetochore proteins

Chromosome segregation is powered by the kinetochore, a large macromolecular structure assembled on centromeric chromatin. Attachment of sister chromatids to microtubules is mediated by the highly conserved tripartite KMN (acronym for KNL-1-Mis12-Ndc80) kinetochore network. In the fission yeast Schizosaccharomyces pombe, the equivalent complex is called NMS (Ndc80-MIND-Spc7). Here, we show that not all components of the NMS complex had been identified previously. A 10th NMS component exists, the essential Sos7 protein, which is a genetic and physical interaction partner of Spc7. The analysis of sos7 kinetochore-null mutant yeast strains demonstrated that Sos7 is central to NMS function. In particular, Sos7 is required for kinetochore targeting of Spc7 as well as components of the MIND complex. sos7 mutant strains show severe chromosome missegregation phenotypes and have compromised microtubule-kinetochore interactions. Sos7 is the founding member of a functionally conserved fungal kinetochore family not present in the point centromere carrying Saccharomycotina clusters, suggesting that the new Sos7 family might be a signature motif of fungi with regional centromeres.

Genome and transcriptome analysis of the food-yeast Candida utilis

The industrially important food-yeast Candida utilis is a Crabtree effect-negative yeast used to produce valuable chemicals and recombinant proteins. In the present study, we conducted whole genome sequencing and phylogenetic analysis of C. utilis, which showed that this yeast diverged long before the formation of the CUG and Saccharomyces/Kluyveromyces clades. In addition, we performed comparative genome and transcriptome analyses using next-generation sequencing, which resulted in the identification of genes important for characteristic phenotypes of C. utilis such as those involved in nitrate assimilation, in addition to the gene encoding the functional hexose transporter. We also found that an antisense transcript of the alcohol dehydrogenase gene, which in silico analysis did not predict to be a functional gene, was transcribed in the stationary-phase, suggesting a novel system of repression of ethanol production. These findings should facilitate the development of more sophisticated systems for the production of useful reagents using C. utilis.

A genome-scale metabolic model of the lipid-accumulating yeast Yarrowia lipolytica

BACKGROUND

Yarrowia lipolytica is an oleaginous yeast which has emerged as an important microorganism for several biotechnological processes, such as the production of organic acids, lipases and proteases. It is also considered a good candidate for single-cell oil production. Although some of its metabolic pathways are well studied, its metabolic engineering is hindered by the lack of a genome-scale model that integrates the current knowledge about its metabolism.

RESULTS

Combining in silico tools and expert manual curation, we have produced an accurate genome-scale metabolic model for Y. lipolytica. Using a scaffold derived from a functional metabolic model of the well-studied but phylogenetically distant yeast S. cerevisiae, we mapped conserved reactions, rewrote gene associations, added species-specific reactions and inserted specialized copies of scaffold reactions to account for species-specific expansion of protein families. We used physiological measures obtained under lab conditions to validate our predictions.

CONCLUSIONS

Y. lipolytica iNL895 represents the first well-annotated metabolic model of an oleaginous yeast, providing a base for future metabolic improvement, and a starting point for the metabolic reconstruction of other species in the Yarrowia clade and other oleaginous yeasts.

Fructanase and fructosyltransferase activity of non-Saccharomyces yeasts isolated from fermenting musts of Mezcal

Fructanase and fructosyltransferase are interesting for the tequila process and prebiotics production (functional food industry). In this study, one hundred thirty non-Saccharomyces yeasts isolated from "Mezcal de Oaxaca" were screened for fructanase and fructosyltransferase activity. On solid medium, fifty isolates grew on Agave tequilana fructans (ATF), inulin or levan. In liquid media, inulin and ATF induced fructanase activities of between 0.02 and 0.27U/ml depending of yeast isolate. High fructanase activity on sucrose was observed for Kluyveromyces marxianus and Torulaspora delbrueckii, while the highest fructanase activity on inulin and ATF was observed for Issatchenkia orientalis, Cryptococcus albidus, and Candida apicola. Zygosaccharomyces bisporus and Candida boidinii had a high hydrolytic activity on levan. Sixteen yeasts belonging to K. marxianus, T. delbrueckii and C. apicola species were positive for fructosyltransferase activity. Mezcal microbiota proved to showed to be a source for new fructanase and fructosyltransferases with potential application in the tequila and food industry.

Pichia sorbitophila, an Interspecies Yeast Hybrid, Reveals Early Steps of Genome Resolution After Polyploidization

Polyploidization is an important process in the evolution of eukaryotic genomes, but ensuing molecular mechanisms remain to be clarified. Autopolyploidization or whole-genome duplication events frequently are resolved in resulting lineages by the loss of single genes from most duplicated pairs, causing transient gene dosage imbalance and accelerating speciation through meiotic infertility. Allopolyploidization or formation of interspecies hybrids raises the problem of genetic incompatibility (Bateson-Dobzhansky-Muller effect) and may be resolved by the accumulation of mutational changes in resulting lineages. In this article, we show that an osmotolerant yeast species, Pichia sorbitophila, recently isolated in a concentrated sorbitol solution in industry, illustrates this last situation. Its genome is a mosaic of homologous and homeologous chromosomes, or parts thereof, that corresponds to a recently formed hybrid in the process of evolution. The respective parental contributions to this genome were characterized using existing variations in GC content. The genomic changes that occurred during the short period since hybrid formation were identified (e.g., loss of heterozygosity, unilateral loss of rDNA, reciprocal exchange) and distinguished from those undergone by the two parental genomes after separation from their common ancestor (i.e., NUMT (NUclear sequences of MiTochondrial origin) insertions, gene acquisitions, gene location movements, reciprocal translocation). We found that the physiological characteristics of this new yeast species are determined by specific but unequal contributions of its two parents, one of which could be identified as very closely related to an extant Pichia farinosa strain.

Genome evolution in the eremothecium clade of the Saccharomyces complex revealed by comparative genomics

We used comparative genomics to elucidate the genome evolution within the pre-whole-genome duplication genus Eremothecium. To this end, we sequenced and assembled the complete genome of Eremothecium cymbalariae, a filamentous ascomycete representing the Eremothecium type strain. Genome annotation indicated 4712 gene models and 143 tRNAs. We compared the E. cymbalariae genome with that of its relative, the riboflavin overproducer Ashbya (Eremothecium) gossypii, and the reconstructed yeast ancestor. Decisive changes in the Eremothecium lineage leading to the evolution of the A. gossypii genome include the reduction from eight to seven chromosomes, the downsizing of the genome by removal of 10% or 900 kb of DNA, mostly in intergenic regions, the loss of a TY3-Gypsy-type transposable element, the re-arrangement of mating-type loci, and a massive increase of its GC content. Key species-specific events are the loss of MNN1-family of mannosyltransferases required to add the terminal fourth and fifth α-1,3-linked mannose residue to O-linked glycans and genes of the Ehrlich pathway in E. cymbalariae and the loss of ZMM-family of meiosis-specific proteins and acquisition of riboflavin overproduction in A. gossypii. This reveals that within the Saccharomyces complex genome, evolution is not only based on genome duplication with subsequent gene deletions and chromosomal rearrangements but also on fungi associated with specific environments (e.g. involving fungal-insect interactions as in Eremothecium), which have encountered challenges that may be reflected both in genome streamlining and their biosynthetic potential.

Evolutionary erosion of yeast sex chromosomes by mating-type switching accidents

We investigate yeast sex chromosome evolution by comparing genome sequences from 16 species in the family Saccharomycetaceae, including data from genera Tetrapisispora, Kazachstania, Naumovozyma, and Torulaspora. We show that although most yeast species contain a mating-type (MAT) locus and silent HML and HMR loci structurally analogous to those of Saccharomyces cerevisiae, their detailed organization is highly variable and indicates that the MAT locus is a deletion hotspot. Over evolutionary time, chromosomal genes located immediately beside MAT have continually been deleted, truncated, or transposed to other places in the genome in a process that is gradually shortening the distance between MAT and HML. Each time a gene beside MAT is removed by deletion or transposition, the next gene on the chromosome is brought into proximity with MAT and is in turn put at risk for removal. This process has also continually replaced the triplicated sequence regions, called Z and X, that allow HML and HMR to be used as templates for DNA repair at MAT during mating-type switching. We propose that the deletion and transposition events are caused by evolutionary accidents during mating-type switching, combined with natural selection to keep MAT and HML on the same chromosome. The rate of deletion accelerated greatly after whole-genome duplication, probably because genes were redundant and could be deleted without requiring transposition. We suggest that, despite its mutational cost, switching confers an evolutionary benefit by providing a way for an isolated germinating spore to reform spores if the environment is too poor.

A consistent phylogenetic backbone for the fungi

The kingdom of fungi provides model organisms for biotechnology, cell biology, genetics, and life sciences in general. Only when their phylogenetic relationships are stably resolved, can individual results from fungal research be integrated into a holistic picture of biology. However, and despite recent progress, many deep relationships within the fungi remain unclear. Here, we present the first phylogenomic study of an entire eukaryotic kingdom that uses a consistency criterion to strengthen phylogenetic conclusions. We reason that branches (splits) recovered with independent data and different tree reconstruction methods are likely to reflect true evolutionary relationships. Two complementary phylogenomic data sets based on 99 fungal genomes and 109 fungal expressed sequence tag (EST) sets analyzed with four different tree reconstruction methods shed light from different angles on the fungal tree of life. Eleven additional data sets address specifically the phylogenetic position of Blastocladiomycota, Ustilaginomycotina, and Dothideomycetes, respectively. The combined evidence from the resulting trees supports the deep-level stability of the fungal groups toward a comprehensive natural system of the fungi. In addition, our analysis reveals methodologically interesting aspects. Enrichment for EST encoded data—a common practice in phylogenomic analyses—introduces a strong bias toward slowly evolving and functionally correlated genes. Consequently, the generalization of phylogenomic data sets as collections of randomly selected genes cannot be taken for granted. A thorough characterization of the data to assess possible influences on the tree reconstruction should therefore become a standard in phylogenomic analyses.

Identification of conserved gene clusters in multiple genomes based on synteny and homology

Background

Uncovering the relationship between the conserved chromosomal segments and the functional relatedness of elements within these segments is an important question in computational genomics. We build upon the series of works on gene teams and homology teams.

Results

Our primary contribution is a local sliding-window SYNS (SYNtenic teamS) algorithm that refines an existing family structure into orthologous sub-families by analyzing the neighborhoods around the members of a given family with a locally sliding window. The neighborhood analysis is done by computing conserved gene clusters. We evaluate our algorithm on the existing homologous families from the Genolevures database over five genomes of the Hemyascomycete phylum.

Conclusions

The result is an efficient algorithm that works on multiple genomes, considers paralogous copies of genes and is able to uncover orthologous clusters even in distant genomes. Resulting orthologous clusters are comparable to those obtained by manual curation.

IONS: Identification of Orthologs by Neighborhood and Similarity-an Automated Method to Identify Orthologs in Chromosomal Regions of Common Evolutionary Ancestry and its Application to Hemiascomycetous Yeasts

Comparative sequence analysis is widely used to infer gene function and study genome evolution and requires proper ortholog identification across different genomes. We have developed a program for the Identification of Orthologs in one-to-one relationship by Neighborhood and Similarity (IONS) between closely related species. The algorithm combines two levels of evidence to determine co-ancestrality at the genome scale: sequence similarity and shared neighborhood. The method was initially designed to provide anchor points for syntenic blocks within the Génolevures project concerning nine hemiascomycetous yeasts (about 50,000 genes) and is applicable to different input databases. Comparison based on use of a Rand index shows that the results are highly consistent with the pillars of the Yeast Gene Order Browser, a manually curated database. Compared with SYNERGY, another algorithm reporting homology relationships, our method’s main advantages are its automation and the absence of dataset-dependent parameters, facilitating consistent integration of newly released genomes.

The evolution of gene expression regulatory networks in yeasts

Gene regulation is a major source of phenotypic diversity between and within species. This aspect of evolution has long been addressed from the sole point of view of the genome sequence. The incredible development of transcriptomics approaches now allows one to actually study the topology and the properties of regulatory networks on an evolutionary perspective. This new discipline is called comparative functional genomics or comparative transcriptomics. This article reviews some of the main advances made in this field, using yeast species, and especially the species sequenced in the frame of the Genolevures program, as a model.

Systematic discovery of unannotated genes in 11 yeast species using a database of orthologous genomic segments

Background

In standard BLAST searches, no information other than the sequences of the query and the database entries is considered. However, in situations where two genes from different species have only borderline similarity in a BLAST search, the discovery that the genes are located within a region of conserved gene order (synteny) can provide additional evidence that they are orthologs. Thus, for interpreting borderline search results, it would be useful to know whether the syntenic context of a database hit is similar to that of the query. This principle has often been used in investigations of particular genes or genomic regions, but to our knowledge it has never been implemented systematically.

Results

We made use of the synteny information contained in the Yeast Gene Order Browser database for 11 yeast species to carry out a systematic search for protein-coding genes that were overlooked in the original annotations of one or more yeast genomes but which are syntenic with their orthologs. Such genes tend to have been overlooked because they are short, highly divergent, or contain introns. The key features of our software - called SearchDOGS - are that the database entries are classified into sets of genomic segments that are already known to be orthologous, and that very weak BLAST hits are retained for further analysis if their genomic location is similar to that of the query. Using SearchDOGS we identified 595 additional protein-coding genes among the 11 yeast species, including two new genes in Saccharomyces cerevisiae. We found additional genes for the mating pheromone a-factor in six species including Kluyveromyces lactis.

Conclusions

SearchDOGS has proven highly successful for identifying overlooked genes in the yeast genomes. We anticipate that our approach can be adapted for study of further groups of species, such as bacterial genomes. More generally, the concept of doing sequence similarity searches against databases to which external information has been added may prove useful in other settings.

Mechanisms of chromosome number evolution in yeast

The whole-genome duplication (WGD) that occurred during yeast evolution changed the basal number of chromosomes from 8 to 16. However, the number of chromosomes in post-WGD species now ranges between 10 and 16, and the number in non-WGD species (Zygosaccharomyces, Kluyveromyces, Lachancea, and Ashbya) ranges between 6 and 8. To study the mechanism by which chromosome number changes, we traced the ancestry of centromeres and telomeres in each species. We observe only two mechanisms by which the number of chromosomes has decreased, as indicated by the loss of a centromere. The most frequent mechanism, seen 8 times, is telomere-to-telomere fusion between two chromosomes with the concomitant death of one centromere. The other mechanism, seen once, involves the breakage of a chromosome at its centromere, followed by the fusion of the two arms to the telomeres of two other chromosomes. The only mechanism by which chromosome number has increased in these species is WGD. Translocations and inversions have cycled telomere locations, internalizing some previously telomeric genes and creating novel telomeric locations. Comparison of centromere structures shows that the length of the CDEII region is variable between species but uniform within species. We trace the complete rearrangement history of the Lachancea kluyveri genome since its common ancestor with Saccharomyces and propose that its exceptionally low level of rearrangement is a consequence of the loss of the non-homologous end joining (NHEJ) DNA repair pathway in this species.

The number of chromosomes in organisms often changes over evolutionary time. To study how the number changes, we compare several related species of yeast that share a common ancestor roughly 150 million years ago and have varying numbers of chromosomes. By inferring ancestral genome structures, we examine the changes in location of centromeres and telomeres, key elements that biologically define chromosomes. Their locations change over time by rearrangements of chromosome segments. By following these rearrangements, we trace an evolutionary path between existing centromeres and telomeres to those in the ancestral genomes, allowing us to identify the specific evolutionary events that caused changes in chromosome number. We show that, in these yeasts, chromosome number has generally decreased over time except for one notable exception: an event in an ancestor of several species where the whole genome was duplicated. Chromosome number reduction occurs by the simultaneous removal of a centromere from a chromosome and fusion of the rest of the chromosome to another that contains a working centromere. This process also results in telomere removal and the movement of genes from the ends of chromosomes to new locations in the middle of chromosomes.

Reinventing heterochromatin in budding yeasts: Sir2 and the origin recognition complex take center stage

The transcriptional silencing of the cryptic mating-type loci in Saccharomyces cerevisiae is one of the best-studied models of repressive heterochromatin. However, this type of heterochromatin, which is mediated by the Sir proteins, has a distinct molecular composition compared to the more ubiquitous type of heterochromatin found in Schizosaccharomyces pombe, other fungi, animals, and plants and characterized by the presence of HP1 (heterochromatin protein 1). This review discusses how the loss of important heterochromatin proteins, including HP1, in the budding yeast lineage presented an evolutionary opportunity for the development and diversification of alternative varieties of heterochromatin, in which the conserved deacetylase Sir2 and the replication protein Orc1 play key roles. In addition, we highlight how this diversification has been facilitated by gene duplications and has contributed to adaptations in lifestyle.

The intronome of budding yeasts

Whatever their abundance in genomes, spliceosomal introns are the signature of eukaryotic genes. The sequence of Saccharomyces cerevisiae, achieved fifteen years ago, revealed that this yeast has very few introns, but conserved intron boundaries typical for an intron definition mechanism. With the improvement and the development of new sequencing technologies, yeast genomes have been extensively sequenced during the last decade. We took advantage of this plethora of data to compile and assess the intron content of the protein-coding genes of 13 genomes representative of the evolution of hemiascomycetous yeasts. We first observed that intron paucity is a general rule and that the fastest evolving genomes tend to lose their introns more rapidly (e.g. S. cerevisiae versus Yarrowia lipolytica). Noticeable differences were also confirmed for 5' splice sites and branch point sites (BP) as well as for the relative position of the BP. These changes seemed to be correlated with the lineage specific evolution of splicing factors.

Phosphorylation network rewiring by gene duplication

In a comprehensive analysis of phosphoregulatory evolution in yeast, the authors observe that phosphorylation sites tend to be lost after gene duplication and protein network turnover reshuffles kinase–substrate relationships over time.

Elucidating how complex regulatory networks have assembled during evolution requires a detailed understanding of the evolutionary dynamics that follow gene duplication events, including changes in post-translational modifications. We compared the phosphorylation profiles of paralogous proteins in the budding yeast Saccharomyces cerevisiae to that of a species that diverged from the budding yeast before the duplication of those genes. We found that 100 million years of post-duplication divergence are sufficient for the majority of phosphorylation sites to be lost or gained in one paralog or the other, with a strong bias toward losses. However, some losses may be partly compensated for by the evolution of other phosphosites, as paralogous proteins tend to preserve similar numbers of phosphosites over time. We also found that up to 50% of kinase–substrate relationships may have been rewired during this period. Our results suggest that after gene duplication, proteins tend to subfunctionalize at the level of post-translational regulation and that even when phosphosites are preserved, there is a turnover of the kinases that phosphorylate them.

Transposable elements in yeasts

With the development of new sequencing technologies in the past decade, yeast genomes have been extensively sequenced and their structures investigated. Transposable elements (TEs) are ubiquitous in eukaryotes and constitute a limited part of yeast genomes. However, due to their ability to move in genomes and generate dispersed repeated sequences, they contribute to modeling yeast genomes and thereby induce plasticity. This review assesses the TE contents of yeast genomes investigated so far. Their diversity and abundance at the inter- and intraspecific levels are presented, and their effects on gene expression and genome stability is considered. Recent results concerning TE-host interactions are also analyzed.

Yeast cell morphology and sexual reproduction--a short overview and some considerations

Over the decades, basic research in life sciences has profited greatly from the study of the small unicellular fungal species Saccharomyces cerevisiae. This yeast turned out to be key for the identification and understanding of molecular mechanisms that underlay the basic functions of all eukaryotic cells. These include, but are not limited to, the regulatory mechanisms behind cellular reproduction (cell cycle control), cellular morphogenesis (cell polarity, cytoskeleton and membrane trafficking) and the management of cellular information (chromosome biology, transcription and translation). Rapid access to genomic information of many yeast species, combined with bioinformatics analyses, provide information on the evolutionary history of yeasts and the molecular ancestry of their constituents. The availability of a comprehensive list of experimental procedures for these organisms presents now a unique opportunity to learn about variations of molecular processes on an evolutionary scale. Yeast cell morphology is another interesting factor, since cellular shapes influence the interactions with the environment on all levels. In this overview article I provide a short summary of the relevant aspects of yeast cell morphology, in particular in relation to one of the most influencing processes, cellular reproduction by mating and meiosis.

Yeasty clocks: dating genomic changes in yeasts

Calibration of clocks to date evolutionary changes is of primary importance for comparative genomics. In the absence of fossil records, the dating of changes during yeast genome evolution can only rely on the properties of the genomes themselves, given the uncertainty of extrapolations using clocks from other organisms. In this work, we use the experimentally determined mutational rate of Saccharomyces cerevisiae to calculate the numbers of successive generations corresponding to observed sequence polymorphism between strains or species of other yeasts. We then examine synteny conservation across the entire subphylum of Saccharomycotina yeasts, and compare this second clock based on chromosomal rearrangements with the first one based on sequence divergence. A non-linear relationship is observed, that interestingly also applies to insects although, for equivalent sequence divergence, their rate of chromosomal rearrangements is higher than that of yeasts.

The GreenCut2 resource, a phylogenomically derived inventory of proteins specific to the plant lineage

The plastid is a defining structure of photosynthetic eukaryotes and houses many plant-specific processes, including the light reactions, carbon fixation, pigment synthesis, and other primary metabolic processes. Identifying proteins associated with catalytic, structural, and regulatory functions that are unique to plastid-containing organisms is necessary to fully define the scope of plant biochemistry. Here, we performed phylogenomics on 20 genomes to compile a new inventory of 597 nucleus-encoded proteins conserved in plants and green algae but not in non-photosynthetic organisms. 286 of these proteins are of known function, whereas 311 are not characterized. This inventory was validated as applicable and relevant to diverse photosynthetic eukaryotes using an additional eight genomes from distantly related plants (including Micromonas, Selaginella, and soybean). Manual curation of the known proteins in the inventory established its importance to plastid biochemistry. To predict functions for the 52% of proteins of unknown function, we used sequence motifs, subcellular localization, co-expression analysis, and RNA abundance data. We demonstrate that 18% of the proteins in the inventory have functions outside the plastid and/or beyond green tissues. Although 32% of proteins in the inventory have homologs in all cyanobacteria, unexpectedly, 30% are eukaryote-specific. Finally, 8% of the proteins of unknown function share no similarity to any characterized protein and are plant lineage-specific. We present this annotated inventory of 597 proteins as a resource for functional analyses of plant-specific biochemistry.

The Awesome Power of Yeast Evolutionary Genetics: New Genome Sequences and Strain Resources for the Saccharomyces sensu stricto Genus

High-quality, well-annotated genome sequences and standardized laboratory strains fuel experimental and evolutionary research. We present improved genome sequences of three species of Saccharomyces sensu stricto yeasts: S. bayanus var. uvarum (CBS 7001), S. kudriavzevii (IFO 1802^T and ZP 591), and S. mikatae (IFO 1815^T), and describe their comparison to the genomes of S. cerevisiae and S. paradoxus. The new sequences, derived by assembling millions of short DNA sequence reads together with previously published Sanger shotgun reads, have vastly greater long-range continuity and far fewer gaps than the previously available genome sequences. New gene predictions defined a set of 5261 protein-coding orthologs across the five most commonly studied Saccharomyces yeasts, enabling a re-examination of the tempo and mode of yeast gene evolution and improved inferences of species-specific gains and losses. To facilitate experimental investigations, we generated genetically marked, stable haploid strains for all three of these Saccharomyces species. These nearly complete genome sequences and the collection of genetically marked strains provide a valuable toolset for comparative studies of gene function, metabolism, and evolution, and render Saccharomyces sensu stricto the most experimentally tractable model genus. These resources are freely available and accessible through www.SaccharomycesSensuStricto.org.

Gene dispensability

Genome-wide mutagenesis studies indicate that up to about 90% of genes in bacteria and 80% in eukaryotes can be inactivated individually leaving an organism viable, often seemingly unaffected. Several strategies are used to learn what these apparently dispensable genes contribute to fitness. Assays of growth under hundreds of physical and chemical stresses are among the most effective experimental approaches. Comparative studies of genomic DNA sequences continue to be valuable in discriminating between the core bacterial genome and the more variable niche-specific genes. The concept of the core genome appears currently unfeasible for eukaryotes but progress has been made in understanding why they contain numerous gene duplicates.

Species-specific PCR primers for the rapid identification of yeasts of the genus Zygosaccharomyces

Species-specific primer pairs that produce a single band of known product size have been developed for members of the Zygosaccharomyces clade including Zygosaccharomyces bailii, Zygosaccharomyces bisporus, Zygosaccharomyces kombuchaensis, Zygosaccharomyces lentus, Zygosaccharomyces machadoi, Zygosaccharomyces mellis and Zygosaccharomyces rouxii. An existing primer pair for the provisional new species Zygosaccharomyces pseudorouxii has been confirmed as specific. The HIS3 gene, encoding imidazole-glycerolphosphate dehydratase, was used as the target gene. This housekeeping gene evolves slowly and is thus well conserved among different isolates, but shows a significant number of base pair changes between even closely related species, sufficient for species-specific primer design. The primers were tested on type and wild strains of the genus Zygosaccharomyces and on members of the Saccharomycetaceae. Sequencing of the D1/D2 region of rDNA was used to confirm the identification of all nonculture collection isolates. This approach used extracted genomic DNA, but in practice, it can be used efficiently with a rapid colony PCR protocol. The method also successfully detected known and new hybrid strains of Z. rouxii and Z. pseudorouxii. The method is rapid, robust and inexpensive. It requires little expertise by the user and is thus useful for preliminary, large-scale screens.

From sequence to function: Insights from natural variation in budding yeasts

Background

Natural variation offers a powerful approach for assigning function to DNA sequence—a pressing challenge in the age of high throughput sequencing technologies.

Scope of Review

Here we review comparative genomic approaches that are bridging the sequence–function and genotype–phenotype gaps. Reverse genomic approaches aim to analyse sequence to assign function, whereas forward genomic approaches start from a phenotype and aim to identify the underlying genotype responsible.

Major Conclusions

Comparative genomic approaches, pioneered in budding yeasts, have resulted in dramatic improvements in our understanding of the function of both genes and regulatory sequences. Analogous studies in other systems, including humans, demonstrate the ubiquity of comparative genomic approaches. Recently, forward genomic approaches, exploiting natural variation within yeast populations, have started to offer powerful insights into how genotype influences phenotype and even the ability to predict phenotypes.

General Significance

Comparative genomic experiments are defining the fundamental rules that govern complex traits in natural populations from yeast to humans.

This article is part of a Special Issue entitled Systems Biology of Microorganisms.