Comparative Genomics of Saccharomyces cerevisiae Natural Isolates for Bioenergy Production

Lignocellulosic plant material is a viable source of biomass to produce alternative energy including ethanol and other biofuels. However, several factors—including toxic byproducts from biomass pretreatment and poor fermentation of xylose and other pentose sugars—currently limit the efficiency of microbial biofuel production. To begin to understand the genetic basis of desirable traits, we characterized three strains of Saccharomyces cerevisiae with robust growth in a pretreated lignocellulosic hydrolysate or tolerance to stress conditions relevant to industrial biofuel production, through genome and transcriptome sequencing analysis. All stress resistant strains were highly mosaic, suggesting that genetic admixture may contribute to novel allele combinations underlying these phenotypes. Strain-specific gene sets not found in the lab strain were functionally linked to the tolerances of particular strains. Furthermore, genes with signatures of evolutionary selection were enriched for functional categories important for stress resistance and included stress-responsive signaling factors. Comparison of the strains’ transcriptomic responses to heat and ethanol treatment—two stresses relevant to industrial bioethanol production—pointed to physiological processes that were related to particular stress resistance profiles. Many of the genotype-by-environment expression responses occurred at targets of transcription factors with signatures of positive selection, suggesting that these strains have undergone positive selection for stress tolerance. Our results generate new insights into potential mechanisms of tolerance to stresses relevant to biofuel production, including ethanol and heat, present a backdrop for further engineering, and provide glimpses into the natural variation of stress tolerance in wild yeast strains.

[1]  Akira Nishimura,et al.  Whole-Genome Sequencing of Sake Yeast Saccharomyces cerevisiae Kyokai no. 7 , 2011, DNA research : an international journal for rapid publication of reports on genes and genomes.

[2]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[3]  Edith D. Wong,et al.  Saccharomyces Genome Database: the genomics resource of budding yeast , 2011, Nucleic Acids Res..

[4]  Julie C. Mitchell,et al.  chipD: a web tool to design oligonucleotide probes for high-density tiling arrays , 2010, Nucleic Acids Res..

[5]  M. Johnston,et al.  Remarkably ancient balanced polymorphisms in a multi-locus gene network , 2010, Nature.

[6]  Audrey P Gasch,et al.  Stress-activated genomic expression changes serve a preparative role for impending stress in yeast. , 2008, Molecular biology of the cell.

[7]  S. Casaregola,et al.  Amplification of a Zygosaccharomyces bailii DNA Segment in Wine Yeast Genomes by Extrachromosomal Circular DNA Formation , 2011, PloS one.

[8]  J. Plotkin,et al.  The Population Genetics of dN/dS , 2008, PLoS genetics.

[9]  Ash A. Alizadeh,et al.  Genome-wide analysis of DNA copy number variation in breast cancer using DNA microarrays , 1999, Nature Genetics.

[10]  Ronald W. Davis,et al.  Functional profiling of the Saccharomyces cerevisiae genome , 2002, Nature.

[11]  R. Durbin,et al.  Efficient de novo assembly of large genomes using compressed data structures. , 2012, Genome research.

[12]  Genome sequencing and comparative analysis of Saccharomyces cerevisiae strain YJM 789 , 2007 .

[13]  J. Townsend,et al.  Population structure and gene evolution in Saccharomyces cerevisiae. , 2006, FEMS yeast research.

[14]  J. Warringer,et al.  Genetic Basis of Variations in Nitrogen Source Utilization in Four Wine Commercial Yeast Strains , 2013, PloS one.

[15]  L. Carey,et al.  Coevolution Trumps Pleiotropy: Carbon Assimilation Traits Are Independent of Metabolic Network Structure in Budding Yeast , 2013, PloS one.

[16]  Saeed Tavazoie,et al.  Molecular Systems Biology 6; Article number 378; doi:10.1038/msb.2010.33 Citation: Molecular Systems Biology 6:378 , 2022 .

[17]  Gavin Sherlock,et al.  Bulk Segregant Analysis by High-Throughput Sequencing Reveals a Novel Xylose Utilization Gene from Saccharomyces cerevisiae , 2010, PLoS genetics.

[18]  Gianni Liti,et al.  Advances in Quantitative Trait Analysis in Yeast , 2012, PLoS genetics.

[19]  Ronald W. Davis,et al.  Whole-genome sequencing of the efficient industrial fuel-ethanol fermentative Saccharomyces cerevisiae strain CAT-1 , 2012, Molecular Genetics and Genomics.

[20]  Brigitte Cambon,et al.  Eukaryote-to-eukaryote gene transfer events revealed by the genome sequence of the wine yeast Saccharomyces cerevisiae EC1118 , 2009, Proceedings of the National Academy of Sciences.

[21]  S. Hara,et al.  Cloning and nucleotide sequence of the KHS killer gene of Saccharomyces cerevisiae. , 1990, Agricultural and biological chemistry.

[22]  Andrew W. Murray,et al.  Rapid Expansion and Functional Divergence of Subtelomeric Gene Families in Yeasts , 2010, Current Biology.

[23]  Justin C. Fay,et al.  Genomic Sequence Diversity and Population Structure of Saccharomyces cerevisiae Assessed by RAD-seq , 2013, G3: Genes, Genomes, Genetics.

[24]  Bruce E Dale,et al.  Multifaceted characterization of cell wall decomposition products formed during ammonia fiber expansion (AFEX) and dilute acid based pretreatments. , 2010, Bioresource technology.

[25]  Alok J. Saldanha,et al.  Java Treeview - extensible visualization of microarray data , 2004, Bioinform..

[26]  F. Dietrich,et al.  The Reacquisition of Biotin Prototrophy in Saccharomyces cerevisiae Involved Horizontal Gene Transfer, Gene Duplication and Gene Clustering , 2007, Genetics.

[27]  Audrey P. Gasch,et al.  Natural Variation in the Yeast Glucose-Signaling Network Reveals a New Role for the Mig3p Transcription Factor , 2012, G3: Genes | Genomes | Genetics.

[28]  P. Donnelly,et al.  Inference of population structure using multilocus genotype data. , 2000, Genetics.

[29]  S. Salzberg,et al.  Versatile and open software for comparing large genomes , 2004, Genome Biology.

[30]  David Botstein,et al.  GO: : TermFinder--open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes , 2004, Bioinform..

[31]  Chikara Furusawa,et al.  Comprehensive phenotypic analysis for identification of genes affecting growth under ethanol stress in Saccharomyces cerevisiae. , 2009, FEMS yeast research.

[32]  Yan Zhou,et al.  Enhancement of the Initial Rate of Ethanol Fermentation Due to Dysfunction of Yeast Stress Response Components Msn2p and/or Msn4p , 2010, Applied and Environmental Microbiology.

[33]  Steven Salzberg,et al.  TigrScan and GlimmerHMM: two open source ab initio eukaryotic gene-finders , 2004, Bioinform..

[34]  J. Nielsen,et al.  Adaptively evolved yeast mutants on galactose show trade-offs in carbon utilization on glucose. , 2013, Metabolic engineering.

[35]  Miguel C. Teixeira,et al.  Genome-Wide Identification of Saccharomyces cerevisiae Genes Required for Maximal Tolerance to Ethanol , 2009, Applied and Environmental Microbiology.

[36]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[37]  A. Gasch,et al.  Exploiting Natural Variation in Saccharomyces cerevisiae to Identify Genes for Increased Ethanol Resistance , 2010, Genetics.

[38]  A. Gasch Yeast genomic expression studies using DNA microarrays. , 2002, Methods in enzymology.

[39]  Robin D Dowell,et al.  Genotype to Phenotype: A Complex Problem , 2010, Science.

[40]  Anneli Petersson,et al.  Increased tolerance and conversion of inhibitors in lignocellulosic hydrolysates by Saccharomyces cerevisiae , 2007 .

[41]  G. Liti,et al.  Surprisingly diverged populations of Saccharomyces cerevisiae in natural environments remote from human activity , 2012, Molecular ecology.

[42]  M. Takemura,et al.  Saccharomyces cerevisiae Σ1278b Has Novel Genes of the N-Acetyltransferase Gene Superfamily Required for l-Proline Analogue Resistance , 2000, Journal of bacteriology.

[43]  N. Friedman,et al.  Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data , 2011, Nature Biotechnology.

[44]  Mark D. Robinson,et al.  FunSpec: a web-based cluster interpreter for yeast , 2002, BMC Bioinformatics.

[45]  Daniel J. Kvitek,et al.  Variations in Stress Sensitivity and Genomic Expression in Diverse S. cerevisiae Isolates , 2008, PLoS genetics.

[46]  Justin C. Fay,et al.  Incipient Balancing Selection through Adaptive Loss of Aquaporins in Natural Saccharomyces cerevisiae Populations , 2010, PLoS genetics.

[47]  Yoshio Katakura,et al.  Identification of target genes conferring ethanol stress tolerance to Saccharomyces cerevisiae based on DNA microarray data analysis. , 2007, Journal of biotechnology.

[48]  Justin C. Fay,et al.  A Catalog of Neutral and Deleterious Polymorphism in Yeast , 2008, PLoS genetics.

[49]  Peilin Jia,et al.  Genome sequencing and comparative analysis of Saccharomyces cerevisiae strain YJM789 , 2007, Proceedings of the National Academy of Sciences.

[50]  Adrian E. Raftery,et al.  mclust Version 4 for R : Normal Mixture Modeling for Model-Based Clustering , Classification , and Density Estimation , 2012 .

[51]  Anders Blomberg,et al.  Trait Variation in Yeast Is Defined by Population History , 2011, PLoS genetics.

[52]  Venkatesh Balan,et al.  Phenotypic selection of a wild Saccharomyces cerevisiae strain for simultaneous saccharification and co-fermentation of AFEX™ pretreated corn stover , 2013, Biotechnology for Biofuels.

[53]  B. Dale,et al.  The impacts of pretreatment on the fermentability of pretreated lignocellulosic biomass: a comparative evaluation between ammonia fiber expansion and dilute acid pretreatment , 2009, Biotechnology for biofuels.

[54]  B. Barrell,et al.  Life with 6000 Genes , 1996, Science.

[55]  Leonid Kruglyak,et al.  Comprehensive polymorphism survey elucidates population structure of Saccharomyces cerevisiae , 2009, Nature.

[56]  Fabiana M. Duarte,et al.  Genome structure of a Saccharomyces cerevisiae strain widely used in bioethanol production. , 2009, Genome research.

[57]  Benjamin M. Bolstad,et al.  affy - analysis of Affymetrix GeneChip data at the probe level , 2004, Bioinform..

[58]  Jun-Yi Leu,et al.  A Tradeoff Drives the Evolution of Reduced Metal Resistance in Natural Populations of Yeast , 2011, PLoS genetics.

[59]  David Botstein,et al.  Yeast metabolic and signaling genes are required for heat-shock survival and have little overlap with the heat-induced genes , 2013, Proceedings of the National Academy of Sciences.

[60]  Anthony R. Borneman,et al.  Whole-Genome Comparison Reveals Novel Genetic Elements That Characterize the Genome of Industrial Strains of Saccharomyces cerevisiae , 2011, PLoS genetics.

[61]  Ash A. Alizadeh,et al.  Genome-wide analysis of DNA copy-number changes using cDNA microarrays , 1999, Nature Genetics.

[62]  D. Hartl,et al.  Population genetic variation in genome-wide gene expression. , 2003, Molecular biology and evolution.

[63]  L. Jønson,et al.  Genome‐wide identification of genes required for growth of Saccharomyces cerevisiae under ethanol stress , 2006, Yeast.

[64]  Martin Dragosits,et al.  Adaptive laboratory evolution – principles and applications for biotechnology , 2013, Microbial Cell Factories.

[65]  Sandrine Dudoit,et al.  Polygenic and directional regulatory evolution across pathways in Saccharomyces , 2010, Proceedings of the National Academy of Sciences.

[66]  A. Gasch,et al.  Harnessing Genetic Diversity in Saccharomyces cerevisiae for Fermentation of Xylose in Hydrolysates of Alkaline Hydrogen Peroxide-Pretreated Biomass , 2013, Applied and Environmental Microbiology.

[67]  T. Borodina,et al.  Transcriptome analysis by strand-specific sequencing of complementary DNA , 2009, Nucleic acids research.

[68]  K. Iwashita,et al.  The AWA1 Gene Is Required for the Foam-Forming Phenotype and Cell Surface Hydrophobicity of Sake Yeast , 2002, Applied and Environmental Microbiology.

[69]  A. F. Bennett,et al.  An experimental test of evolutionary trade-offs during temperature adaptation , 2007, Proceedings of the National Academy of Sciences.

[70]  Angus H. Forgan,et al.  Comparative genome analysis of a Saccharomyces cerevisiae wine strain. , 2008, FEMS yeast research.

[71]  Robert P. Davey,et al.  Population genomics of domestic and wild yeasts , 2008, Nature.

[72]  D. Botstein,et al.  Genomic expression programs in the response of yeast cells to environmental changes. , 2000, Molecular biology of the cell.

[73]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[74]  Jun Yu,et al.  KaKs_Calculator 2.0: A Toolkit Incorporating Gamma-Series Methods and Sliding Window Strategies , 2010, Genom. Proteom. Bioinform..

[75]  Jean YH Yang,et al.  Bioconductor: open software development for computational biology and bioinformatics , 2004, Genome Biology.

[76]  M. Aigle,et al.  RTM1: a member of a new family of telomeric repeated genes in yeast. , 1995, Genetics.

[77]  Paul M. Magwene,et al.  Outcrossing, mitotic recombination, and life-history trade-offs shape genome evolution in Saccharomyces cerevisiae , 2011, Proceedings of the National Academy of Sciences.