Comparative genotyping of the Saccharomyces cerevisiae laboratory strains S288C and CEN.PK113-7D using oligonucleotide microarrays

Construction and evaluation of gRNA arrays for multiplex CRISPR-Cas9

Endonuclease system CRISPR-Cas9 represents a powerful toolbox for the budding yeast's Saccharomyces cerevisiae genome perturbation. The resulting double-strand breaks are preferentially repaired via highly efficient homologous recombination, which subsequently leads to marker-free genome editing. The goal of this study was to evaluate precise targeting of multiple loci simultaneously. To construct an array of independently expressing guide RNAs (gRNAs), the genes encoding them were assembled through a BioBrick construction procedure. We designed a multiplex CRISPR-Cas9 system for targeting 6 marker genes, whereby the gRNA array was expressed from a single plasmid. To evaluate the performance of the gRNA array, the activity of the designed system was assessed by the success rate of the introduction of perturbations within the target loci: successful gRNA expression, followed by target DNA double-strand breaks formation and their repair by homologous recombination led to premature termination of the coding sequence of the marker genes, resulting in the prevention of growth of the transformants on the corresponding selection media. In conclusion, we successfully introduced up to five simultaneous perturbations within single cells of yeast S. cerevisiae using the multiplex CRISPR-Cas9 system. While this has been done before, we here present an alternative sequential BioBrick assembly with the capability to accommodate many highly similar gRNA-expression cassettes, and an exhaustive evaluation of their performance.

Chromatin Regulators Ahc1p and Eaf3p Positively Influence Nitrogen Metabolism in Saccharomyces cerevisiae

There is a complex regulatory network of nitrogen metabolism in Saccharomyces cerevisiae, and many details of this regulatory network have not been revealed. This study explored the global regulation of nitrogen metabolism in S. cerevisiae from an epigenetic perspective. Comparative transcriptome analysis of S. cerevisiae S288C treated with 30 nitrogen sources identified nine chromatin regulators (CRs) that responded significantly to different nitrogen sources. Functional analysis showed that among the CRs identified, Ahc1p and Eaf3p promoted the utilization of non-preferred nitrogen sources through global regulation of nitrogen metabolism. Ahc1p regulated nitrogen metabolism through amino acid transport, nitrogen catabolism repression (NCR), and the Ssy1p-Ptr3p-Ssy5p signaling sensor system. Eaf3p regulated nitrogen metabolism via amino acid transport and NCR. The regulatory mechanisms of the effects of Ahc1p and Eaf3p on nitrogen metabolism depended on the function of their histone acetyltransferase complex ADA and NuA4. These epigenetic findings provided new insights for a deeper understanding of the nitrogen metabolism regulatory network in S. cerevisiae.

Gene co-expression network analysis reveals the positive impact of endocytosis and mitochondria-related genes over nitrogen metabolism in Saccharomyces cerevisiae

Nitrogen metabolism is essential for most cellular activities. Therefore, a deep understanding of its regulatory mechanisms is necessary for the efficient utilization of nitrogen sources for Saccharomyces cerevisiae. In this study, a gene co-expression network was constructed for S. cerevisiae S288C with different nitrogen sources. From this, a key gene co-expression module related to nitrogen source preference utilization was obtained, and 10 hub genes centrally located in the co-expression network were identified. Functional studies verified that the endocytosis-related genes CAP1 and END3 significantly increased the utilization of multiple non-preferred amino acids and reduced the accumulation of the harmful nitrogen metabolite precursor urea by regulating amino acid transporters and TOR pathway. The mitochondria-related gene ATP12, MRPL22, MRP1 and NAM9 significantly increased the utilization of multiple non-preferred amino acids and reduced accumulation of the urea by coordinately regulating nitrogen catabolism repression, Ssy1p-Ptr3p-Ssy5p signaling sensor system, amino acid transporters, TOR pathway and urea metabolism-related pathways. Furthermore, these data revealed the potential positive effects of endocytosis and mitochondrial ribosomes protein translation on nitrogen source preference. This study provides new analytical perspectives for complex regulatory networks involving nitrogen metabolism in S. cerevisiae.

Potassium and Sodium Salt Stress Characterization in the Yeasts Saccharomyces cerevisiae, Kluyveromyces marxianus, and Rhodotorula toruloides

Biotechnology requires efficient microbial cell factories. The budding yeast Saccharomyces cerevisiae is a vital cell factory, but more diverse cell factories are essential for the sustainable use of natural resources. Here, we benchmarked nonconventional yeasts Kluyveromyces marxianus and Rhodotorula toruloides against S. cerevisiae strains CEN.PK and W303 for their responses to potassium and sodium salt stress. We found an inverse relationship between the maximum growth rate and the median cell volume that was responsive to salt stress. The supplementation of K⁺ to CEN.PK cultures reduced Na⁺ toxicity and increased the specific growth rate 4-fold. The higher K⁺ and Na⁺ concentrations impaired ethanol and acetate metabolism in CEN.PK and acetate metabolism in W303. In R. toruloides cultures, these salt supplementations induced a trade-off between glucose utilization and cellular aggregate formation. Their combined use increased the beta-carotene yield by 60% compared with that of the reference. Neural network-based image analysis of exponential-phase cultures showed that the vacuole-to-cell volume ratio increased with increased cell volume for W303 and K. marxianus but not for CEN.PK and R. toruloides in response to salt stress. Our results provide insights into common salt stress responses in yeasts and will help design efficient bioprocesses. IMPORTANCE Characterization of microbial cell factories under industrially relevant conditions is crucial for designing efficient bioprocesses. Salt stress, typical in industrial bioprocesses, impinges upon cell volume and affects productivity. This study presents an open-source neural network-based analysis method to evaluate volumetric changes using yeast optical microscopy images. It allows quantification of cell and vacuole volumes relevant to cellular physiology. On applying salt stress in yeasts, we found that the combined use of K⁺ and Na⁺ improves the cellular fitness of Saccharomyces cerevisiae strain CEN.PK and increases the beta-carotene productivity in Rhodotorula toruloides, a commercially important antioxidant and a valuable additive in foods.

Population Analysis and Evolution of Saccharomyces cerevisiae Mitogenomes

The study of mitogenomes allows the unraveling of some paths of yeast evolution that are often not exposed when analyzing the nuclear genome. Although both nuclear and mitochondrial genomes are known to determine phenotypic diversity and fitness, no concordance has yet established between the two, mainly regarding strains' technological uses and/or geographical distribution. In the current work, we proposed a new method to align and analyze yeast mitogenomes, overcoming current difficulties that make it impossible to obtain comparable mitogenomes for a large number of isolates. To this end, 12,016 mitogenomes were considered, and we developed a novel approach consisting of the design of a reference sequence intended to be comparable between all mitogenomes. Subsequently, the population structure of 6646 Saccharomyces cerevisiae mitogenomes was assessed. Results revealed the existence of particular clusters associated with the technological use of the strains, in particular regarding clinical isolates, laboratory strains, and yeasts used for wine-associated activities. As far as we know, this is the first time that a positive concordance between nuclear and mitogenomes has been reported for S. cerevisiae, in terms of strains' technological applications. The results obtained highlighted the importance of including the mtDNA genome in evolutionary analysis, in order to clarify the origin and history of yeast species.

Nanopore sequencing enables near-complete de novo assembly of Saccharomyces cerevisiae reference strain CEN.PK113-7D

The haploid Saccharomyces cerevisiae strain CEN.PK113-7D is a popular model system for metabolic engineering and systems biology research. Current genome assemblies are based on short-read sequencing data scaffolded based on homology to strain S288C. However, these assemblies contain large sequence gaps, particularly in subtelomeric regions, and the assumption of perfect homology to S288C for scaffolding introduces bias. In this study, we obtained a near-complete genome assembly of CEN.PK113-7D using only Oxford Nanopore Technology's MinION sequencing platform. Fifteen of the 16 chromosomes, the mitochondrial genome and the 2-μm plasmid are assembled in single contigs and all but one chromosome starts or ends in a telomere repeat. This improved genome assembly contains 770 Kbp of added sequence containing 248 gene annotations in comparison to the previous assembly of CEN.PK113-7D. Many of these genes encode functions determining fitness in specific growth conditions and are therefore highly relevant for various industrial applications. Furthermore, we discovered a translocation between chromosomes III and VIII that caused misidentification of a MAL locus in the previous CEN.PK113-7D assembly. This study demonstrates the power of long-read sequencing by providing a high-quality reference assembly and annotation of CEN.PK113-7D and places a caveat on assumed genome stability of microorganisms.

Parallelised online biomass monitoring in shake flasks enables efficient strain and carbon source dependent growth characterisation of Saccharomyces cerevisiae

BACKGROUND

Baker's yeast, Saccharomyces cerevisiae, as one of the most often used workhorses in biotechnology has been developed into a huge family of application optimised strains in the last decades. Increasing numbers of strains render their characterisation highly challenging, even with the simple methods of growth-based analytics. Here we present a new sensor system for the automated, non-invasive and parallelisable monitoring of biomass in continuously shaken shake flask cultures, called CGQ ("cell growth quantifier"). The CGQ implements a dynamic approach of backscattered light measurement, allowing for efficient and accurate growth-based strain characterisation, as exemplarily demonstrated for the four most commonly used laboratory and industrial yeast strains, BY4741, W303-1A, CEN.PK2-1C and Ethanol Red.

RESULTS

Growth experiments revealed distinct carbon source utilisation differences between the investigated S. cerevisiae strains. Phenomena such as diauxic shifts, morphological changes and oxygen limitations were clearly observable in the growth curves. A strictly monotonic non-linear correlation of OD600 and the CGQ's backscattered light intensities was found, with strain-to-strain as well as growth-phase related differences. The CGQ measurements showed high resolution, sensitivity and smoothness even below an OD600 of 0.2 and were furthermore characterised by low background noise and signal drift in combination with high reproducibility.

CONCLUSIONS

With the CGQ, shake flask fermentations can be automatically monitored regarding biomass and growth rates with high resolution and parallelisation. This makes the CGQ a valuable tool for growth-based strain characterisation and development. The exceptionally high resolution allows for the identification of distinct metabolic differences and shifts as well as for morphologic changes. Applications that will benefit from that kind of automatized biomass monitoring include, amongst many others, the characterization of deregulated native or integrated heterologous pathways, the fast detection of co-fermentation as well as the realisation of rational and growth-data driven evolutionary engineering approaches.

The Saccharomyces cerevisiae poly(A) binding protein Pab1 as a target for eliciting stress tolerant phenotypes

When exploited as cell factories, Saccharomyces cerevisiae cells are exposed to harsh environmental stresses impairing titer, yield and productivity of the fermentative processes. The development of robust strains therefore represents a pivotal challenge for the implementation of cost-effective bioprocesses. Altering master regulators of general cellular rewiring represents a possible strategy to evoke shaded potential that may accomplish the desirable features. The poly(A) binding protein Pab1, as stress granules component, was here selected as the target for obtaining widespread alterations in mRNA metabolism, resulting in stress tolerant phenotypes. Firstly, we demonstrated that the modulation of Pab1 levels improves robustness against different stressors. Secondly, the mutagenesis of PAB1 and the application of a specific screening protocol on acetic acid enriched medium allowed the isolation of the further ameliorated mutant pab1 A60-9. These findings pave the way for a novel approach to unlock industrially promising phenotypes through the modulation of a post-transcriptional regulatory element.

The genome sequence of the popular hexose-transport-deficient Saccharomyces cerevisiae strain EBY.VW4000 reveals LoxP/Cre-induced translocations and gene loss

Saccharomyces cerevisiae harbours a large group of tightly controlled hexose transporters with different characteristics. Construction and characterization of S. cerevisiae EBY.VW4000, a strain devoid of glucose import, was a milestone in hexose-transporter research. This strain has become a widely used platform for discovery and characterization of transporters from a wide range of organisms. To abolish glucose uptake, 21 genes were knocked out, involving 16 successive deletion rounds with the LoxP/Cre system. Although such intensive modifications are known to increase the risk of genome alterations, the genome of EBY.VW4000 has hitherto not been characterized. Based on a combination of whole genome sequencing, karyotyping and molecular confirmation, the present study reveals that construction of EBY.VW4000 resulted in gene losses and chromosomal rearrangements. Recombinations between the LoxP scars have led to the assembly of four neo-chromosomes, truncation of two chromosomes and loss of two subtelomeric regions. Furthermore, sporulation and spore germination are severely impaired in EBY.VW4000. Karyotyping of the EBY.VW4000 lineage retraced its current chromosomal architecture to four translocations events occurred between the 6th and the 12th rounds of deletion. The presented data facilitate further studies on EBY.VW4000 and highlight the risks of genome alterations associated with repeated use of the LoxP/Cre system.

Benchmarking two commonly used Saccharomyces cerevisiae strains for heterologous vanillin-β-glucoside production

The yeast Saccharomyces cerevisiae is a widely used eukaryotic model organism and a key cell factory for production of biofuels and wide range of chemicals. From the broad palette of available yeast strains, the most popular are those derived from laboratory strain S288c and the industrially relevant CEN.PK strain series. Importantly, in recent years these two strains have been subjected to comparative "-omics" analyzes pointing out significant genotypic and phenotypic differences. It is therefore possible that the two strains differ significantly with respect to their potential as cell factories for production of specific compounds. To examine this possibility, we have reconstructed a de novo vanillin-β-glucoside pathway in an identical manner in S288c and CEN.PK strains. Characterization of the two resulting strains in two standard conditions revealed that the S288c background strain produced up to 10-fold higher amounts of vanillin-β-glucoside compared to CEN.PK. This study demonstrates that yeast strain background may play a major role in the outcome of newly developed cell factories for production of a given product.

CRISPR-Cas system enables fast and simple genome editing of industrial Saccharomyces cerevisiae strains

There is a demand to develop 3rd generation biorefineries that integrate energy production with the production of higher value chemicals from renewable feedstocks. Here, robust and stress-tolerant industrial strains of Saccharomyces cerevisiae will be suitable production organisms. However, their genetic manipulation is challenging, as they are usually diploid or polyploid. Therefore, there is a need to develop more efficient genetic engineering tools. We applied a CRISPR–Cas9 system for genome editing of different industrial strains, and show simultaneous disruption of two alleles of a gene in several unrelated strains with the efficiency ranging between 65% and 78%. We also achieved simultaneous disruption and knock-in of a reporter gene, and demonstrate the applicability of the method by designing lactic acid-producing strains in a single transformation event, where insertion of a heterologous gene and disruption of two endogenous genes occurred simultaneously. Our study provides a foundation for efficient engineering of industrial yeast cell factories.

•

We developed CRISPR–Cas9-based system for gene disruptions in industrial yeast.

•

We showed high rate of disruption efficiency in unrelated industrial strains.

•

Gene knock-in may be performed simultaneously with gene disruption.

•

Use of the described Cas9-based system results in marker-free stable genetic modifications.

•

The method was applied for single-step construction of lactic acid-producing strains.

De novo sequencing, assembly and analysis of the genome of the laboratory strain Saccharomyces cerevisiae CEN.PK113-7D, a model for modern industrial biotechnology

Saccharomyces cerevisiae CEN.PK 113-7D is widely used for metabolic engineering and systems biology research in industry and academia. We sequenced, assembled, annotated and analyzed its genome. Single-nucleotide variations (SNV), insertions/deletions (indels) and differences in genome organization compared to the reference strain S. cerevisiae S288C were analyzed. In addition to a few large deletions and duplications, nearly 3000 indels were identified in the CEN.PK113-7D genome relative to S288C. These differences were overrepresented in genes whose functions are related to transcriptional regulation and chromatin remodelling. Some of these variations were caused by unstable tandem repeats, suggesting an innate evolvability of the corresponding genes. Besides a previously characterized mutation in adenylate cyclase, the CEN.PK113-7D genome sequence revealed a significant enrichment of non-synonymous mutations in genes encoding for components of the cAMP signalling pathway. Some phenotypic characteristics of the CEN.PK113-7D strains were explained by the presence of additional specific metabolic genes relative to S288C. In particular, the presence of the BIO1 and BIO6 genes correlated with a biotin prototrophy of CEN.PK113-7D. Furthermore, the copy number, chromosomal location and sequences of the MAL loci were resolved. The assembled sequence reveals that CEN.PK113-7D has a mosaic genome that combines characteristics of laboratory strains and wild-industrial strains.

Laboratory evolution of new lactate transporter genes in a jen1Δ mutant of Saccharomyces cerevisiae and their identification as ADY2 alleles by whole-genome resequencing and transcriptome analysis

Laboratory evolution is a powerful approach in applied and fundamental yeast research, but complete elucidation of the molecular basis of evolved phenotypes remains a challenge. In this study, DNA microarray-based transcriptome analysis and whole-genome resequencing were used to investigate evolution of novel lactate transporters in Saccharomyces cerevisiae that can replace Jen1p, the only documented S. cerevisiae lactate transporter. To this end, a jen1Δ mutant was evolved for growth on lactate in serial batch cultures. Two independent evolution experiments yielded growth on lactate as sole carbon source (0.14 and 0.18 h(-1) , respectively). Transcriptome analysis did not provide leads, but whole-genome resequencing showed different single-nucleotide changes (C755G/Leu219Val and C655G/Ala252Gly) in the acetate transporter gene ADY2. Introduction of these ADY2 alleles in a jen1Δ ady2Δ strain enabled growth on lactate (0.14 h(-1) for Ady2p(Leu219Val) and 0.12 h(-1) for Ady2p(Ala252Gly) ), demonstrating that these alleles of ADY2 encode efficient lactate transporters. Depth of coverage of DNA sequencing, combined with karyotyping, gene deletions and diagnostic PCR, showed that an isochromosome III (c. 475 kb) with two additional copies of ADY2(C755G) had been formed via crossover between retrotransposons YCLWΔ15 and YCRCΔ6. The isochromosome formation shows how even short periods of selective pressure can cause substantial karyotype changes.

Genome-wide analytical approaches for reverse metabolic engineering of industrially relevant phenotypes in yeast

Successful reverse engineering of mutants that have been obtained by nontargeted strain improvement has long presented a major challenge in yeast biotechnology. This paper reviews the use of genome-wide approaches for analysis of Saccharomyces cerevisiae strains originating from evolutionary engineering or random mutagenesis. On the basis of an evaluation of the strengths and weaknesses of different methods, we conclude that for the initial identification of relevant genetic changes, whole genome sequencing is superior to other analytical techniques, such as transcriptome, metabolome, proteome, or array-based genome analysis. Key advantages of this technique over gene expression analysis include the independency of genome sequences on experimental context and the possibility to directly and precisely reproduce the identified changes in naive strains. The predictive value of genome-wide analysis of strains with industrially relevant characteristics can be further improved by classical genetics or simultaneous analysis of strains derived from parallel, independent strain improvement lineages.

Functional genomics for food fermentation processes

This review describes recent scientific and technological drivers of food fermentation research. In addition, a number of practical implications of the results of this development will be highlighted. The first part of the manuscript elaborates on the message that genome sequence information gives us an unprecedented view on the biodiversity of microbes in food fermentation. This information can be made applicable for tailoring relevant characteristics of food products through fermentation. The second part deals with the integration of genome sequence data into metabolic models and the use of these models for a number of topics that are relevant for food fermentation processes. The final part will be about metagenomics approaches to reveal the complexity and understand the functionality of undefined complex microbial consortia used in a diverse range of food fermentation processes.

High-frequency genetic contents variations in clinical Candida albicans isolates

Genome plasticity is a hallmark of Candida albicans and is believed to be an adaptation strategy. But the extent of such genomic variability is not well investigated. In this study, genetic contents of clinical C. albicans isolates were investigated at whole-genome level with array-based comparative genomic hybridization (array CGH) technology. It was revealed that C. albicans possessed variations of genetic contents, as well as aneuploidy. The variable genes were scattered across the chromosomes, as well clustered in particular regions, including sub-telomeric regions, retrotransposon-insertion sites and a variable region on chromosome 6.

Saccharomyces cerevisiae BY4741 and W303-1A laboratory strains differ in salt tolerance

Saccharomyces cerevisiae yeast cells serve as a model to elucidate the bases of salt tolerance and potassium homeostasis regulation in eukaryotic cells. In this study, we show that two widely used laboratory strains, BY4741 and W303-1A, differ not only in cell size and volume but also in their relative plasma-membrane potential (estimated with a potentiometric fluorescent dye diS-C3(3) and as Hygromycin B sensitivity) and tolerance to alkali-metal cations. W303-1A cells and their mutant derivatives lacking either uptake (trk1 trk2) or efflux (nha1) systems for alkali-metal cations are more tolerant to toxic sodium and lithium cations but also more sensitive to higher external concentrations of potassium than BY4741 cells and their mutants. Moreover, our results suggest that though the two strains do not differ in the total potassium content, the regulation of intracellular potassium homeostasis is probably not the same in BY4741 and W303-1A cells.

Defining species specific genome differences in malaria parasites

Background

In recent years a number of genome sequences for different plasmodium species have become available. This has allowed the identification of numerous conserved genes across the different species and has significantly enhanced our understanding of parasite biology. In contrast little is known about species specific differences between the different genomes partly due to the lower sequence coverage and therefore relatively poor annotation of some of the draft genomes particularly the rodent malarias parasite species.

Results

To improve the current annotation and gene identification status of the draft genomes of P. berghei, P. chabaudi and P. yoelii, we performed genome-wide comparisons between these three species. Through analyses via comparative genome hybridizations using a newly designed pan-rodent array as well as in depth bioinformatics analysis, we were able to improve on the coverage of the draft rodent parasite genomes by detecting orthologous genes between these related rodent parasite species. More than 1,000 orthologs for P. yoelii were now newly associated with a P. falciparum gene. In addition to extending the current core gene set for all plasmodium species this analysis also for the first time identifies a relatively small number of genes that are unique to the primate malaria parasites while a larger gene set is uniquely conserved amongst the rodent malaria parasites.

Conclusions

These findings allow a more thorough investigation of the genes that are important for host specificity in malaria.

Systematic identification of balanced transposition polymorphisms in Saccharomyces cerevisiae

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

Balanced structural polymorphisms are differences in the relative arrangement of genomic features within species that do not affect DNA copy number. Little is known about their prevalence or importance because they are difficult to observe. Here, we present a novel methodology for systematically identifying such polymorphisms based on the idea that single-copy DNA that occupies different genomic locations in two parents will segregate independently during meiosis and will therefore reveal itself as a copy number difference among a fraction of progeny. Comparative hybridization reveals multiple balanced structural polymorphisms that involve changes to gene order in two strains of yeast; the results are independently validated using paired-end whole genome shotgun sequencing. The longest transposed segment we identify comprises 13.5 kb near the telomere of chromosome XV in the S288C reference strain and contains several annotated genes. We map the location of this polymorphism in the non-reference strain using genome-wide genotypic data, which also reveals an appreciable frequency of ectopic recombination among transposed segment pairs. The breakpoints of the remaining polymorphisms are localized by the paired-end sequence data. Our work provides proof-of-principle for a very general approach to systematically identify all balanced genomic polymorphisms in two different genotypes and is a starting point for understanding the frequency, evolutionary origins, and functional consequences of this seldom-studied class of genomic structural variation in eukaryotes.

Comparative genomics of wild type yeast strains unveils important genome diversity

Background

Genome variability generates phenotypic heterogeneity and is of relevance for adaptation to environmental change, but the extent of such variability in natural populations is still poorly understood. For example, selected Saccharomyces cerevisiae strains are variable at the ploidy level, have gene amplifications, changes in chromosome copy number, and gross chromosomal rearrangements. This suggests that genome plasticity provides important genetic diversity upon which natural selection mechanisms can operate.

Results

In this study, we have used wild-type S. cerevisiae (yeast) strains to investigate genome variation in natural and artificial environments. We have used comparative genome hybridization on array (aCGH) to characterize the genome variability of 16 yeast strains, of laboratory and commercial origin, isolated from vineyards and wine cellars, and from opportunistic human infections. Interestingly, sub-telomeric instability was associated with the clinical phenotype, while Ty element insertion regions determined genomic differences of natural wine fermentation strains. Copy number depletion of ASP3 and YRF1 genes was found in all wild-type strains. Other gene families involved in transmembrane transport, sugar and alcohol metabolism or drug resistance had copy number changes, which also distinguished wine from clinical isolates.

Conclusion

We have isolated and genotyped more than 1000 yeast strains from natural environments and carried out an aCGH analysis of 16 strains representative of distinct genotype clusters. Important genomic variability was identified between these strains, in particular in sub-telomeric regions and in Ty-element insertion sites, suggesting that this type of genome variability is the main source of genetic diversity in natural populations of yeast. The data highlights the usefulness of yeast as a model system to unravel intraspecific natural genome diversity and to elucidate how natural selection shapes the yeast genome.

Phenotypic screening, transcriptional profiling, and comparative genomic analysis of an invasive and non-invasive strain of Candida albicans

Background

Invasion of host tissue by the human fungal pathogen Candida albicans is an important step during the development of candidosis. However, not all C. albicans strains possess the same invasive and virulence properties. For example, the two clinical isolates SC5314 and ATCC10231 differ in their ability to invade host tissue and cause experimental infections. Strain SC5314 is invasive whereas strain ATCC10231 is non-invasive and strongly attenuated in virulence compared to SC5314. In this study we compare the in vitro phenotypic, transcriptional and genomic profiles of these two widely used laboratory strains in order to determine the principal biological and genetic properties responsible for their differential virulence.

Results

In all media tested, the two strains showed the same metabolic flexibility, stress resistance, adhesion properties and hydrolytic enzyme secretion in vitro. However, differences were observed in response to cell-surface disturbing agents and alkaline pH. Furthermore, reduced hyphal formation in strain ATCC10231 under certain conditions correlated with reduced invasive properties in an in vitro invasion assay and a reduced ability to invade epithelial tissue. Despite these diverse phenotypic properties, no substantial genomic differences were detected by comparative genome hybridisation within the open reading frames. However, in vitro transcriptional profiling displayed major differences in the gene expression of these two strains, even under normal in vitro growth conditions.

Conclusion

Our data suggest that the reason for differential virulence of C. albicans strains is not due to the absence of specific genes, but rather due to differences in the expression, function or activity of common genes.

Growth temperature exerts differential physiological and transcriptional responses in laboratory and wine strains of Saccharomyces cerevisiae

Laboratory strains of Saccharomyces cerevisiae have been widely used as a model for studying eukaryotic cells and mapping the molecular mechanisms of many different human diseases. Industrial wine yeasts, on the other hand, have been selected on the basis of their adaptation to stringent environmental conditions and the organoleptic properties that they confer to wine. Here, we used a two-factor design to study the responses of a standard laboratory strain, CEN.PK113-7D, and an industrial wine yeast strain, EC1118, to growth temperatures of 15 degrees C and 30 degrees C in nitrogen-limited, anaerobic, steady-state chemostat cultures. Physiological characterization revealed that the growth temperature strongly impacted the biomass yield of both strains. Moreover, we found that the wine yeast was better adapted to mobilizing resources for biomass production and that the laboratory yeast exhibited higher fermentation rates. To elucidate mechanistic differences controlling the growth temperature response and underlying adaptive mechanisms between the strains, DNA microarrays and targeted metabolome analysis were used. We identified 1,007 temperature-dependent genes and 473 strain-dependent genes. The transcriptional response was used to identify highly correlated gene expression subnetworks within yeast metabolism. We showed that temperature differences most strongly affect nitrogen metabolism and the heat shock response. A lack of stress response element-mediated gene induction, coupled with reduced trehalose levels, indicated that there was a decreased general stress response at 15 degrees C compared to that at 30 degrees C. Differential responses among strains were centered on sugar uptake, nitrogen metabolism, and expression of genes related to organoleptic properties. Our study provides global insight into how growth temperature affects differential physiological and transcriptional responses in laboratory and wine strains of S. cerevisiae.

A systems biology perspective of wine fermentations

The yeast Saccharomyces cerevisiae is an important industrial microorganism. Nowadays, it is being used as a cell factory for the production of pharmaceuticals such as insulin, although this yeast has long been utilized in the bakery to raise dough, and in the production of alcoholic beverages, fermenting the sugars derived from rice, wheat, barley, corn and grape juice. S. cerevisiae has also been extensively used as a model eukaryotic system. In the last decade, genomic techniques have revealed important features of its molecular biology. For example, DNA array technologies are routinely used for determining gene expression levels in cells under different physiological conditions or environmental stimuli. Laboratory strains of S. cerevisiae are different from wine strains. For instance, laboratory yeasts are unable to completely transform all the sugar in the grape must into ethanol under winemaking conditions. In fact, standard culture conditions are usually very different from winemaking conditions, where multiple stresses occur simultaneously and sequentially throughout the fermentation. The response of wine yeasts to these stimuli differs in some aspects from laboratory strains, as suggested by the increasing number of studies in functional genomics being conducted on wine strains. In this paper we review the most recent applications of post-genomic techniques to understand yeast physiology in the wine industry. We also report recent advances in wine yeast strain improvement and propose a reference framework for integration of genomic information, bioinformatic tools and molecular biology techniques for cellular and metabolic engineering. Finally, we discuss the current state and future perspectives for using 'modern' biotechnology in the wine industry.

Molecular analysis of maltotriose active transport and fermentation by Saccharomyces cerevisiae reveals a determinant role for the AGT1 permease

Incomplete and/or sluggish maltotriose fermentation causes both quality and economic problems in the ale-brewing industry. Although it has been proposed previously that the sugar uptake must be responsible for these undesirable phenotypes, there have been conflicting reports on whether all the known alpha-glucoside transporters in Saccharomyces cerevisiae (MALx1, AGT1, and MPH2 and MPH3 transporters) allow efficient maltotriose utilization by yeast cells. We characterized the kinetics of yeast cell growth, sugar consumption, and ethanol production during maltose or maltotriose utilization by several S. cerevisiae yeast strains (both MAL constitutive and MAL inducible) and by their isogenic counterparts with specific deletions of the AGT1 gene. Our results clearly showed that yeast strains carrying functional permeases encoded by the MAL21, MAL31, and/or MAL41 gene in their plasma membranes were unable to utilize maltotriose. While both high- and low-affinity transport activities were responsible for maltose uptake from the medium, in the case of maltotriose, the only low-affinity (K(m), 36 +/- 2 mM) transport activity was mediated by the AGT1 permease. In conclusion, the AGT1 transporter is required for efficient maltotriose fermentation by S. cerevisiae yeasts, highlighting the importance of this permease for breeding and/or selection programs aimed at improving sluggish maltotriose fermentations.

Physiological and transcriptional responses of Saccharomyces cerevisiae to zinc limitation in chemostat cultures

Transcriptional responses of the yeast Saccharomyces cerevisiae to Zn availability were investigated at a fixed specific growth rate under limiting and abundant Zn concentrations in chemostat culture. To investigate the context dependency of this transcriptional response and eliminate growth rate-dependent variations in transcription, yeast was grown under several chemostat regimens, resulting in various carbon (glucose), nitrogen (ammonium), zinc, and oxygen supplies. A robust set of genes that responded consistently to Zn limitation was identified, and the set enabled the definition of the Zn-specific Zap1p regulon, comprised of 26 genes and characterized by a broader zinc-responsive element consensus (MHHAACCBYNMRGGT) than so far described. Most surprising was the Zn-dependent regulation of genes involved in storage carbohydrate metabolism. Their concerted down-regulation was physiologically relevant as revealed by a substantial decrease in glycogen and trehalose cellular content under Zn limitation. An unexpectedly large number of genes were synergistically or antagonistically regulated by oxygen and Zn availability. This combinatorial regulation suggested a more prominent involvement of Zn in mitochondrial biogenesis and function than hitherto identified.

Studying bacterial genome dynamics using microarray-based comparative genomic hybridization

Genome sequencing has revealed the remarkable amount of genetic diversity that can be encountered in bacterial genomes. In particular, the comparison of genome sequences from closely related strains has uncovered significant differences in gene content, hinting at the dynamic nature of bacterial genomes. The study of these genome dynamics is crucial to leveraging genomic information because the genome sequence of a single bacterial strain may not accurately represent the genome of the species. The dynamic nature of bacterial genome content has required us to apply the concepts of comparative genomics (CG) at the species level. Although direct genome sequence comparisons are an ideal method of performing CG, one current constraint is the limited availability of multiple genome sequences from a given bacterial species. DNA microarray-based comparative genomic hybridization (MCGH), which can be used to determine the presence or absence of thousands of genes in a single hybridization experiment, provides a powerful alternative for determining genome content and has been successfully used to investigate the genome dynamics of a wide number of bacterial species. Although MCGH-based studies have already provided a new vista on bacterial genome diversity, original methods for MCGH have been limited by the absence of novel gene sequences included in the microarray. New applications of the MCGH platform not only promise to accelerate the pace of novel gene discovery but will also help provide an integrated microarray-based approach to the study of bacterial CG.

A complex gene cluster for indole-diterpene biosynthesis in the grass endophyte Neotyphodium lolii

Lolitrems are a structurally diverse group of indole-diterpene mycotoxins synthesized by Epichloë/Neotyphodium endophytes in association with Pooid grasses. Using suppression subtractive hybridization combined with chromosome walking, two clusters of genes for lolitrem biosynthesis were isolated from Neotyphodium lolii, a mutualistic endophyte of perennial ryegrass. The first cluster contains five genes, ltmP, ltmQ, ltmF, ltmC, and ltmB, four of which appear to be orthologues of functionally characterized genes from Penicillium paxilli. The second cluster contains two genes, ltmE and ltmJ, that appear to be unique to lolitrem biosynthesis. The two clusters are separated by a 16 kb AT-rich sequence that includes two imperfect direct repeats. A previously isolated ltm cluster composed of ltmG, ltmM, and ltmK, is linked to these two new clusters by 35 kb of AT-rich retrotransposon relic sequence. All 10 genes at this complex LTM locus were highly expressed in planta but expression was very low or undetectable in mycelia. ltmM and ltmC were shown to be functional orthologues of P. paxilli paxM and paxC, respectively. This work provides a genetic foundation for elucidating the metabolic grid responsible for the diversity of indole-diterpenes synthesized by N. lolii.

Genome-Wide Detection of Polymorphisms at Nucleotide Resolution with a Single DNA Microarray

A central challenge of genomics is to detect, simply and inexpensively, all differences in sequence among the genomes of individual members of a species. We devised a system to detect all single-nucleotide differences between genomes with the use of data from a single hybridization to a whole-genome DNA microarray. This allowed us to detect a variety of spontaneous single-base pair substitutions, insertions, and deletions, and most (>90%) of the approximately 30,000 known single-nucleotide polymorphisms between two Saccharomyces cerevisiae strains. We applied this approach to elucidate the genetic basis of phenotypic variants and to identify the small number of single-base pair changes accumulated during experimental evolution of yeast.

Birth of a metabolic gene cluster in yeast by adaptive gene relocation

Although most eukaryotic genomes lack operons, they contain some physical clusters of genes that are related in function despite being unrelated in sequence. How these clusters are formed during evolution is unknown. The DAL cluster is the largest metabolic gene cluster in yeast and consists of six adjacent genes encoding proteins that enable Saccharomyces cerevisiae to use allantoin as a nitrogen source. We show here that the DAL cluster was assembled, quite recently in evolutionary terms, through a set of genomic rearrangements that happened almost simultaneously. Six of the eight genes involved in allantoin degradation, which were previously scattered around the genome, became relocated to a single subtelomeric site in an ancestor of S. cerevisiae and Saccharomyces castellii. These genomic rearrangements coincided with a biochemical reorganization of the purine degradation pathway, which switched to importing allantoin instead of urate. This change eliminated urate oxidase, one of several oxygen-consuming enzymes that were lost by yeasts that can grow vigorously in anaerobic conditions. The DAL cluster is located in a domain of modified chromatin involving both H2A.Z histone exchange and Hst1-Sum1-mediated histone deacetylation, and it may be a coadapted gene complex formed by epistatic selection.

A new approach for the analysis of bacterial microarray-based Comparative Genomic Hybridization: insights from an empirical study

BACKGROUND

Microarray-based Comparative Genomic Hybridization (M-CGH) has been used to characterize the extensive intraspecies genetic diversity found in bacteria at the whole-genome level. Although conventional microarray analytical procedures have proved adequate in handling M-CGH data, data interpretation using these methods is based on a continuous character model in which gene divergence and gene absence form a spectrum of decreasing gene conservation levels. However, whereas gene divergence may yet be accompanied by retention in gene function, gene absence invariably leads to loss of function. This distinction, if ignored, leads to a loss in the information to be gained from M-CGH data. We present here results from experiments in which two genome-sequenced strains of C. jejuni were compared against each other using M-CGH. Because the gene content of both strains was known a priori, we were able to closely examine the effects of sequence divergence and gene absence on M-CGH data in order to define analytical parameters for M-CGH data interpretation. This would facilitate the examination of the relative effects of sequence divergence or gene absence in comparative genomics analyses of multiple strains of any species for which genome sequence data and a DNA microarray are available.

RESULTS

As a first step towards improving the analysis of M-CGH data, we estimated the degree of experimental error in a series of experiments in which identical samples were compared against each other by M-CGH. This variance estimate was used to validate a Log Ratio-based methodology for identification of outliers in M-CGH data. We compared two genome strains by M-CGH to examine the effect of probe/target identity on the Log Ratios of signal intensities using prior knowledge of gene divergence and gene absence to establish Log Ratio thresholds for the identification of absent and conserved genes.

CONCLUSION

The results from this empirical study validate the Log Ratio thresholds that have been used in other studies to establish gene divergence/absence. Moreover, the analytical framework presented here enhances the information content derived from M-CGH data by shifting the focus from divergent/absent gene detection to accurate detection of conserved and absent genes. This approach closely aligns the technical limitations of M-CGH analysis with practical limitations on the biological interpretation of comparative genomics data.

Grr1p is required for transcriptional induction of amino acid permease genes and proper transcriptional regulation of genes in carbon metabolism of Saccharomyces cerevisiae

The F-box protein Grr1p is involved in cell cycle regulation, glucose repression and transcriptional induction of the amino acid permease (AAP) gene AGP1. We investigated the role of Grr1p in amino acid-mediated induction of AAP genes by performing batch cultivations with a wild-type strain and a grr1Delta strain and adding citrulline in the exponential phase. Whole-genome transcription analyses were performed on samples from each cultivation, both immediately before and 30 min after citrulline addition. Transcriptional induction of the AAP genes AGP1, BAP2, BAP3, DIP5, GNP1 and TAT1 is fully dependent on Grr1p. Comparison of the grr1Delta strain with the reference strain in the absence of citrulline revealed that GRR1 disruption leads to increased transcription of numerous genes. These encode enzymes in the tricarboxylic acid cycle, the pentose-phosphate pathway and both glucose and starch metabolism. Promoter analysis showed that many of the genes with increased transcription display Mig1p- and/or Msn2p/Msn4p-binding sites. Increased expression of glucose-repressed genes in the grr1Delta strain may be explained by the reduced expression of the hexose transporter genes HXT1, HXT2, HXT3 and HXT4 and a subsequent lowering of the glucose uptake; and the effect of GRR1 deletion on general carbon metabolism may therefore be indirect. Finally, none of the genes known to be primarily involved in cell cycle regulation displayed different expression levels in the grr1Delta cells as compared with the reference strain, suggesting that the role of Grr1p in cell cycle regulation does not include any transcriptional component.

Impact of ‘ome’ analyses on inverse metabolic engineering

Genome-wide or large-scale methodologies employed in functional genomics such as DNA sequencing, transcription profiling, proteomics, and metabolite profiling have become important tools in many metabolic engineering strategies. These techniques allow the identification of genetic differences and insight into their cellular effects. In the field of inverse metabolic engineering mapping of differences between strains with different degree of a certain desired phenotype and subsequent identification of factors conferring that phenotype are an essential part. Therefore, the tools of functional genomics in particular have the potential to promote and expand inverse metabolic engineering. Here, we review the use of functional genomics methods in inverse metabolic engineering, examples are presented, and we discuss the identification of targets for metabolic engineering with low fold changes using these techniques.

Role of transcriptional regulation in controlling fluxes in central carbon metabolism of Saccharomyces cerevisiae. A chemostat culture study

In contrast to batch cultivation, chemostat cultivation allows the identification of carbon source responses without interference by carbon-catabolite repression, accumulation of toxic products, and differences in specific growth rate. This study focuses on the yeast Saccharomyces cerevisiae, grown in aerobic, carbon-limited chemostat cultures. Genome-wide transcript levels and in vivo fluxes were compared for growth on two sugars, glucose and maltose, and for two C2-compounds, ethanol and acetate. In contrast to previous reports on batch cultures, few genes (180 genes) responded to changes of the carbon source by a changed transcript level. Very few transcript levels were changed when glucose as the growth-limiting nutrient was compared with maltose (33 transcripts), or when acetate was compared with ethanol (16 transcripts). Although metabolic flux analysis using a stoichiometric model revealed major changes in the central carbon metabolism, only 117 genes exhibited a significantly different transcript level when sugars and C2-compounds were provided as the growth-limiting nutrient. Despite the extensive knowledge on carbon source regulation in yeast, many of the carbon source-responsive genes encoded proteins with unknown or incompletely characterized biological functions. In silico promoter analysis of carbon source-responsive genes confirmed the involvement of several known transcriptional regulators and suggested the involvement of additional regulators. Transcripts involved in the glyoxylate cycle and gluconeogenesis showed a good correlation with in vivo fluxes. This correlation was, however, not observed for other important pathways, including the pentose-phosphate pathway, tricarboxylic acid cycle, and, in particular, glycolysis. These results indicate that in vivo fluxes in the central carbon metabolism of S. cerevisiae grown in steadystate, carbon-limited chemostat cultures are controlled to a large extent via post-transcriptional mechanisms.