Draft genome sequence of the Saccharomyces cerevisiae SpyCas9 expressing strain IMX2600, a laboratory and platform strain from the CEN.PK lineage for cell-factory research.

The biobased-economy aims to create a circular biotechnology ecosystem to transition from a fossil fuel-based to a sustainable industry based on biomass. For this, new microbial cell-factories are essential. We present the draft genome of the CEN.PK-derived Saccharomyces cerevisiae SpyCas9 expressing strain (IMX2600), that serve as chassis of new cell-factories.

Pantothenate Auxotrophy in a Naturally Occurring Biocontrol Yeast.

The genus Hanseniaspora is characterized by some of the smallest genomes among budding yeasts. These fungi are primarily found on plant surfaces and in fermented products and represent promising biocontrol agents against notorious fungal plant pathogens. In this work, we identify pantothenate auxotrophy of a Hanseniaspora meyeri isolate that shows strong antagonism against the plant pathogen Fusarium oxysporum. Furthermore, strong biocontrol activity in vitro required both pantothenate and biotin in the growth medium. We show that the H. meyeri isolate APC 12.1 can obtain the vitamin from plants and other fungi. The underlying reason for the auxotrophy is the lack of two key pantothenate biosynthesis genes, but six genes encoding putative pantothenate transporters are present in the genome. By constructing and using a Saccharomyces cerevisiae reporter strain, we identified one Hanseniaspora transporter that conferred pantothenate uptake activity to S. cerevisiae. Pantothenate auxotrophy is rare and has been described in only a few bacteria and in S. cerevisiae strains that were isolated from sake. Such auxotrophic strains may seem an unexpected and unlikely choice as potential biocontrol agents, but they may be particularly competitive in their ecological niche and their specific growth requirements are an inherent biocontainment strategy preventing uncontrolled growth in the environment. Auxotrophic strains, such as the H. meyeri isolate APC 12.1, may thus represent a promising strategy for developing biocontrol agents that will be easier to register than prototrophic strains, which are normally used for such applications. IMPORTANCE As a precursor of the essential coenzyme A (CoA), pantothenate is present in all organisms. Plants, bacteria, and fungi are known to synthesize this vitamin, while animals must obtain it through their diet. Pantothenate auxotrophy has not been described in naturally occurring, environmental fungi and is an unexpected property for an antagonistic yeast. Here, we report that yeasts from the genus Hanseniaspora lack key enzymes for pantothenate biosynthesis and identify a transporter responsible for the acquisition of pantothenate from the environment. Hanseniaspora isolates are strong antagonists of fungal plant pathogens. Their pantothenate auxotrophy is a natural biocontainment feature that could make such isolates interesting candidates for new biocontrol approaches and allow easier registration as plant protection agents than prototrophic strains.

Identification of fungal dihydrouracil-oxidase genes by expression in Saccharomyces cerevisiae.

Analysis of predicted fungal proteomes revealed a large family of sequences that showed similarity to the Saccharomyces cerevisiae Class-I dihydroorotate dehydrogenase Ura1, which supports synthesis of pyrimidines under aerobic and anaerobic conditions. However, expression of codon-optimised representatives of this gene family, from the ascomycete Alternaria alternata and the basidiomycete Schizophyllum commune, only supported growth of an S. cerevisiae ura1Δ mutant when synthetic media were supplemented with dihydrouracil. A hypothesis that these genes encode NAD(P)+-dependent dihydrouracil dehydrogenases (EC 1.3.1.1 or 1.3.1.2) was rejected based on absence of complementation in anaerobic cultures. Uracil- and thymine-dependent oxygen consumption and hydrogen-peroxide production by cell extracts of S. cerevisiae strains expressing the A. alternata and S. commune genes showed that, instead, they encode active dihydrouracil oxidases (DHO, EC1.3.3.7). DHO catalyses the reaction dihydrouracil + O2 → uracil + H2O2 and was only reported in the yeast Rhodotorula glutinis (Owaki in J Ferment Technol 64:205-210, 1986). No structural gene for DHO was previously identified. DHO-expressing strains were highly sensitive to 5-fluorodihydrouracil (5F-dhu) and plasmids bearing expression cassettes for DHO were readily lost during growth on 5F-dhu-containing media. These results show the potential applicability of fungal DHO genes as counter-selectable marker genes for genetic modification of S. cerevisiae and other organisms that lack a native DHO. Further research should explore the physiological significance of this enigmatic and apparently widespread fungal enzyme.

Respiratory reoxidation of NADH is a key contributor to high oxygen requirements of oxygen-limited cultures of Ogataea parapolymorpha.

While thermotolerance is an attractive trait for yeasts used in industrial ethanol production, oxygen requirements of known thermotolerant species are incompatible with process requirements. Analysis of oxygen-sufficient and oxygen-limited chemostat cultures of the facultatively fermentative, thermotolerant species Ogataea parapolymorpha showed its minimum oxygen requirements to be an order of magnitude larger than those reported for the thermotolerant yeast Kluyveromyces marxianus. High oxygen requirements of O. parapolymorpha coincided with a near absence of glycerol, a key NADH/NAD+ redox-cofactor-balancing product in many other yeasts, in oxygen-limited cultures. Genome analysis indicated absence of orthologs of the Saccharomyces cerevisiae glycerol-3-phosphate-phosphatase genes GPP1 and GPP2. Co-feeding of acetoin, whose conversion to 2,3-butanediol enables reoxidation of cytosolic NADH, supported a 2.5-fold increase of the biomass concentration in oxygen-limited cultures. An O. parapolymorpha strain in which key genes involved in mitochondrial reoxidation of NADH were inactivated did produce glycerol, but transcriptome analysis did not reveal a clear candidate for a responsible phosphatase. Expression of S. cerevisiae GPD2, which encodes NAD+-dependent glycerol-3-phosphate dehydrogenase, and GPP1 supported increased glycerol production by oxygen-limited chemostat cultures of O. parapolymorpha. These results identify dependence on respiration for NADH reoxidation as a key contributor to unexpectedly high oxygen requirements of O. parapolymorpha.

Genome sequence data of the antagonistic soil-borne yeast Cyberlindnera sargentensis (SHA 17.2).

Cyberlindnera sargentensis strain SHA 17.2, isolated from a Swiss soil sample, exhibited strong antagonistic activity against several plant pathogenic fungi in vitro and was highly competitive against other yeasts in soil. As a basis for identifying the mechanisms underlying its strong antagonistic activity, we have sequenced the genome of C. sargentensis (SHA 17.2) by long- and short read sequencing, de novo assembled them into seven contigs/chromosomes and a mitogenome (total genome size 11.4 Mbp), and annotated 5455 genes. This high-quality genome is the reference for transcriptome and proteome analyses aiming at elucidating the mode of action of C. sargentensis against fungal plant pathogens. It will thus serve as a resource for identifying potential biocontrol genes and performing comparative genomics analyses of yeast genomes.

Class-II dihydroorotate dehydrogenases from three phylogenetically distant fungi support anaerobic pyrimidine biosynthesis.

BACKGROUND: In most fungi, quinone-dependent Class-II dihydroorotate dehydrogenases (DHODs) are essential for pyrimidine biosynthesis. Coupling of these Class-II DHODHs to mitochondrial respiration makes their in vivo activity dependent on oxygen availability. Saccharomyces cerevisiae and closely related yeast species harbor a cytosolic Class-I DHOD (Ura1) that uses fumarate as electron acceptor and thereby enables anaerobic pyrimidine synthesis. Here, we investigate DHODs from three fungi (the Neocallimastigomycete Anaeromyces robustus and the yeasts Schizosaccharomyces japonicus and Dekkera bruxellensis) that can grow anaerobically but, based on genome analysis, only harbor a Class-II DHOD. RESULTS: Heterologous expression of putative Class-II DHOD-encoding genes from fungi capable of anaerobic, pyrimidine-prototrophic growth (Arura9, SjURA9, DbURA9) in an S. cerevisiae ura1Δ strain supported aerobic as well as anaerobic pyrimidine prototrophy. A strain expressing DbURA9 showed delayed anaerobic growth without pyrimidine supplementation. Adapted faster growing DbURA9-expressing strains showed mutations in FUM1, which encodes fumarase. GFP-tagged SjUra9 and DbUra9 were localized to S. cerevisiae mitochondria, while ArUra9, whose sequence lacked a mitochondrial targeting sequence, was localized to the yeast cytosol. Experiments with cell extracts showed that ArUra9 used free FAD and FMN as electron acceptors. Expression of SjURA9 in S. cerevisiae reproducibly led to loss of respiratory competence and mitochondrial DNA. A cysteine residue (C265 in SjUra9) in the active sites of all three anaerobically active Ura9 orthologs was shown to be essential for anaerobic activity of SjUra9 but not of ArUra9. CONCLUSIONS: Activity of fungal Class-II DHODs was long thought to be dependent on an active respiratory chain, which in most fungi requires the presence of oxygen. By heterologous expression experiments in S. cerevisiae, this study shows that phylogenetically distant fungi independently evolved Class-II dihydroorotate dehydrogenases that enable anaerobic pyrimidine biosynthesis. Further structure-function studies are required to understand the mechanistic basis for the anaerobic activity of Class-II DHODs and an observed loss of respiratory competence in S. cerevisiae strains expressing an anaerobically active DHOD from Sch. japonicus.

A squalene-hopene cyclase in Schizosaccharomyces japonicus represents a eukaryotic adaptation to sterol-limited anaerobic environments.

Biosynthesis of sterols, which are key constituents of canonical eukaryotic membranes, requires molecular oxygen. Anaerobic protists and deep-branching anaerobic fungi are the only eukaryotes in which a mechanism for sterol-independent growth has been elucidated. In these organisms, tetrahymanol, formed through oxygen-independent cyclization of squalene by a squalene-tetrahymanol cyclase, acts as a sterol surrogate. This study confirms an early report [C. J. E. A. Bulder, Antonie Van Leeuwenhoek, 37, 353-358 (1971)] that Schizosaccharomyces japonicus is exceptional among yeasts in growing anaerobically on synthetic media lacking sterols and unsaturated fatty acids. Mass spectrometry of lipid fractions of anaerobically grown Sch. japonicus showed the presence of hopanoids, a class of cyclic triterpenoids not previously detected in yeasts, including hop-22(29)-ene, hop-17(21)-ene, hop-21(22)-ene, and hopan-22-ol. A putative gene in Sch. japonicus showed high similarity to bacterial squalene-hopene cyclase (SHC) genes and in particular to those of Acetobacter species. No orthologs of the putative Sch. japonicus SHC were found in other yeast species. Expression of the Sch. japonicus SHC gene (Sjshc1) in Saccharomyces cerevisiae enabled hopanoid synthesis and stimulated anaerobic growth in sterol-free media, thus indicating that one or more of the hopanoids produced by SjShc1 could at least partially replace sterols. Use of hopanoids as sterol surrogates represents a previously unknown adaptation of eukaryotic cells to anaerobic growth. The fast anaerobic growth of Sch. japonicus in sterol-free media is an interesting trait for developing robust fungal cell factories for application in anaerobic industrial processes.

Engineering the thermotolerant industrial yeast Kluyveromyces marxianus for anaerobic growth.

Current large-scale, anaerobic industrial processes for ethanol production from renewable carbohydrates predominantly rely on the mesophilic yeast Saccharomyces cerevisiae. Use of thermotolerant, facultatively fermentative yeasts such as Kluyveromyces marxianus could confer significant economic benefits. However, in contrast to S. cerevisiae, these yeasts cannot grow in the absence of oxygen. Responses of K. marxianus and S. cerevisiae to different oxygen-limitation regimes were analyzed in chemostats. Genome and transcriptome analysis, physiological responses to sterol supplementation and sterol-uptake measurements identified absence of a functional sterol-uptake mechanism as a key factor underlying the oxygen requirement of K. marxianus. Heterologous expression of a squalene-tetrahymanol cyclase enabled oxygen-independent synthesis of the sterol surrogate tetrahymanol in K. marxianus. After a brief adaptation under oxygen-limited conditions, tetrahymanol-expressing K. marxianus strains grew anaerobically on glucose at temperatures of up to 45 °C. These results open up new directions in the development of thermotolerant yeast strains for anaerobic industrial applications.

Critical parameters and procedures for anaerobic cultivation of yeasts in bioreactors and anaerobic chambers.

All known facultatively fermentative yeasts require molecular oxygen for growth. Only in a small number of yeast species, these requirements can be circumvented by supplementation of known anaerobic growth factors such as nicotinate, sterols and unsaturated fatty acids. Biosynthetic oxygen requirements of yeasts are typically small and, unless extensive precautions are taken to minimize inadvertent entry of trace amounts of oxygen, easily go unnoticed in small-scale laboratory cultivation systems. This paper discusses critical points in the design of anaerobic yeast cultivation experiments in anaerobic chambers and laboratory bioreactors. Serial transfer or continuous cultivation to dilute growth factors present in anaerobically pre-grown inocula, systematic inclusion of control strains and minimizing the impact of oxygen diffusion through tubing are identified as key elements in experimental design. Basic protocols are presented for anaerobic-chamber and bioreactor experiments.

Identification of Oxygen-Independent Pathways for Pyridine Nucleotide and Coenzyme A Synthesis in Anaerobic Fungi by Expression of Candidate Genes in Yeast.

Neocallimastigomycetes are unique examples of strictly anaerobic eukaryotes. This study investigates how these anaerobic fungi bypass reactions involved in synthesis of pyridine nucleotide cofactors and coenzyme A that, in canonical fungal pathways, require molecular oxygen. Analysis of Neocallimastigomycetes proteomes identified a candidate l-aspartate-decarboxylase (AdcA) and l-aspartate oxidase (NadB) and quinolinate synthase (NadA), constituting putative oxygen-independent bypasses for coenzyme A synthesis and pyridine nucleotide cofactor synthesis. The corresponding gene sequences indicated acquisition by ancient horizontal gene transfer (HGT) events involving bacterial donors. To test whether these enzymes suffice to bypass corresponding oxygen-requiring reactions, they were introduced into fms1Δ and bna2Δ Saccharomyces cerevisiae strains. Expression of nadA and nadB from Piromyces finnis and adcA from Neocallimastix californiae conferred cofactor prototrophy under aerobic and anaerobic conditions. This study simulates how HGT can drive eukaryotic adaptation to anaerobiosis and provides a basis for elimination of auxotrophic requirements in anaerobic industrial applications of yeasts and fungi. IMPORTANCE NAD (NAD+) and coenzyme A (CoA) are central metabolic cofactors whose canonical biosynthesis pathways in fungi require oxygen. Anaerobic gut fungi of the Neocallimastigomycota phylum are unique eukaryotic organisms that adapted to anoxic environments. Analysis of Neocallimastigomycota genomes revealed that these fungi might have developed oxygen-independent biosynthetic pathways for NAD+ and CoA biosynthesis, likely acquired through horizontal gene transfer (HGT) from prokaryotic donors. We confirmed functionality of these putative pathways under anaerobic conditions by heterologous expression in the yeast Saccharomyces cerevisiae. This approach, combined with sequence comparison, offers experimental insight on whether HGT events were required and/or sufficient for acquiring new traits. Moreover, our results demonstrate an engineering strategy for enabling S. cerevisiae to grow anaerobically in the absence of the precursor molecules pantothenate and nicotinate, thereby contributing to alleviate oxygen requirements and to move closer to prototrophic anaerobic growth of this industrially relevant yeast.

Vitamin requirements and biosynthesis in Saccharomyces cerevisiae.

Chemically defined media for yeast cultivation (CDMY) were developed to support fast growth, experimental reproducibility, and quantitative analysis of growth rates and biomass yields. In addition to mineral salts and a carbon substrate, popular CDMYs contain seven to nine B-group vitamins, which are either enzyme cofactors or precursors for their synthesis. Despite the widespread use of CDMY in fundamental and applied yeast research, the relation of their design and composition to the actual vitamin requirements of yeasts has not been subjected to critical review since their first development in the 1940s. Vitamins are formally defined as essential organic molecules that cannot be synthesized by an organism. In yeast physiology, use of the term "vitamin" is primarily based on essentiality for humans, but the genome of the Saccharomyces cerevisiae reference strain S288C harbours most of the structural genes required for synthesis of the vitamins included in popular CDMY. Here, we review the biochemistry and genetics of the biosynthesis of these compounds by S. cerevisiae and, based on a comparative genomics analysis, assess the diversity within the Saccharomyces genus with respect to vitamin prototrophy.

Origin of Lactose Fermentation in Kluyveromyces lactis by Interspecies Transfer of a Neo-functionalized Gene Cluster during Domestication.

Humans have used yeasts to make cheese and kefir for millennia, but the ability to ferment the milk sugar lactose is found in only a few yeast species, of which the foremost is Kluyveromyces lactis [1]. Two genes, LAC12 (lactose permease) and LAC4 (lactase), are sufficient for lactose uptake and hydrolysis to glucose and galactose [2]. Here, we show that these genes have a complex evolutionary history in the genus Kluyveromyces that is likely the result of human activity during domestication. We show that the ancestral Lac12 was bifunctional, able to import both lactose and cellobiose into the cell. These disaccharides were then hydrolyzed by Lac4 in the case of lactose or Cel2 in the case of cellobiose. A second cellobiose transporter, Cel1, was also present ancestrally. In the K. lactis lineage, the ancestral LAC12 and LAC4 were lost and a separate upheaval in the sister species K. marxianus resulted in loss of CEL1 and quadruplication of LAC12. One of these LAC12 genes became neofunctionalized to encode an efficient lactose transporter capable of supporting fermentation, specifically in dairy strains of K. marxianus, where it formed a LAC4-LAC12-CEL2 gene cluster, although another remained a cellobiose transporter. Then, the ability to ferment lactose was acquired very recently by K. lactis var. lactis by introgression of LAC12 and LAC4 on a 15-kb subtelomeric region from a dairy strain of K. marxianus. The genomic history of the LAC genes shows that strong selective pressures were imposed on yeasts by early dairy farmers.

Snf2 controls pulcherriminic acid biosynthesis and antifungal activity of the biocontrol yeast Metschnikowia pulcherrima.

Metschnikowia pulcherrima synthesises the pigment pulcherrimin, from cyclodileucine (cyclo(Leu-Leu)) as a precursor, and exhibits strong antifungal activity against notorious plant pathogenic fungi. This yeast therefore has great potential for biocontrol applications against fungal diseases; particularly in the phyllosphere where this species is frequently found. To elucidate the molecular basis of the antifungal activity of M. pulcherrima, we compared a wild-type strain with a spontaneously occurring, pigmentless, weakly antagonistic mutant derivative. Whole genome sequencing of the wild-type and mutant strains identified a point mutation that creates a premature stop codon in the transcriptional regulator gene SNF2 in the mutant. Complementation of the mutant strain with the wild-type SNF2 gene restored pigmentation and recovered the strong antifungal activity. Mass spectrometry (UPLC HR HESI-MS) proved the presence of the pulcherrimin precursors cyclo(Leu-Leu) and pulcherriminic acid and identified new precursor and degradation products of pulcherriminic acid and/or pulcherrimin. All of these compounds were identified in the wild-type and complemented strain, but were undetectable in the pigmentless snf2 mutant strain. These results thus identify Snf2 as a regulator of antifungal activity and pulcherriminic acid biosynthesis in M. pulcherrima and provide a starting point for deciphering the molecular functions underlying the antagonistic activity of this yeast.

Coverage-Versus-Length Plots, a Simple Quality Control Step for de Novo Yeast Genome Sequence Assemblies.

Illumina sequencing has revolutionized yeast genomics, with prices for commercial draft genome sequencing now below $200. The popular SPAdes assembler makes it simple to generate a de novo genome assembly for any yeast species. However, whereas making genome assemblies has become routine, understanding what they contain is still challenging. Here, we show how graphing the information that SPAdes provides about the length and coverage of each scaffold can be used to investigate the nature of an assembly, and to diagnose possible problems. Scaffolds derived from mitochondrial DNA, ribosomal DNA, and yeast plasmids can be identified by their high coverage. Contaminating data, such as cross-contamination from other samples in a multiplex sequencing run, can be identified by its low coverage. Scaffolds derived from the bacteriophage PhiX174 and Lambda DNAs that are frequently used as molecular standards in Illumina protocols can also be detected. Assemblies of yeast genomes with high heterozygosity, such as interspecies hybrids, often contain two types of scaffold: regions of the genome where the two alleles assembled into two separate scaffolds and each has a coverage level C, and regions where the two alleles co-assembled (collapsed) into a single scaffold that has a coverage level 2C Visualizing the data with Coverage-vs.-Length (CVL) plots, which can be done using Microsoft Excel or Google Sheets, provides a simple method to understand the structure of a genome assembly and detect aberrant scaffolds or contigs. We provide a Python script that allows assemblies to be filtered to remove contaminants identified in CVL plots.

Population genomics shows no distinction between pathogenic Candida krusei and environmental Pichia kudriavzevii: One species, four names.

We investigated genomic diversity of a yeast species that is both an opportunistic pathogen and an important industrial yeast. Under the name Candida krusei, it is responsible for about 2% of yeast infections caused by Candida species in humans. Bloodstream infections with C. krusei are problematic because most isolates are fluconazole-resistant. Under the names Pichia kudriavzevii, Issatchenkia orientalis and Candida glycerinogenes, the same yeast, including genetically modified strains, is used for industrial-scale production of glycerol and succinate. It is also used to make some fermented foods. Here, we sequenced the type strains of C. krusei (CBS573T) and P. kudriavzevii (CBS5147T), as well as 30 other clinical and environmental isolates. Our results show conclusively that they are the same species, with collinear genomes 99.6% identical in DNA sequence. Phylogenetic analysis of SNPs does not segregate clinical and environmental isolates into separate clades, suggesting that C. krusei infections are frequently acquired from the environment. Reduced resistance of strains to fluconazole correlates with the presence of one gene instead of two at the ABC11-ABC1 tandem locus. Most isolates are diploid, but one-quarter are triploid. Loss of heterozygosity is common, including at the mating-type locus. Our PacBio/Illumina assembly of the 10.8 Mb CBS573T genome is resolved into 5 complete chromosomes, and was annotated using RNAseq support. Each of the 5 centromeres is a 35 kb gene desert containing a large inverted repeat. This species is a member of the genus Pichia and family Pichiaceae (the methylotrophic yeasts clade), and so is only distantly related to other pathogenic Candida species.

Zygosaccharomyces pseudobailii, another yeast interspecies hybrid that regained fertility by damaging one of its MAT loci.

Interspecies hybridization is an important evolutionary mechanism in yeasts. The genus Zygosaccharomyces in particular contains numerous hybrid strains and/or species. Here, we investigated the genome of Zygosaccharomyces strain MT15, an isolate from Maotai-flavor Chinese liquor fermentation. We found that it is an interspecies hybrid and identified it as Zygosaccharomyces pseudobailii. The Z. bailii species complex consists of three species: Z. bailii, which is not a hybrid and whose 10 Mb genome is designated 'A', and two hybrid species Z. parabailii ('AB' genome, 20 Mb) and Z. pseudobailii ('AC' genome, 20 Mb). The A, B and C subgenomes are all approximately 7%-10% different from one another in nucleotide sequence, and are derived from three different parental species. Despite being hybrids, Z. pseudobailii and Z. parabailii are capable of mating and sporulating. We previously showed that Z. parabailii regained fertility when one copy of its MAT locus became broken into two parts, causing the allodiploid hybrid to behave as a haploid gamete. In Z. pseudobailii, we find that a very similar process occurred after hybridization, when a deletion of 1.5 kb inactivated one of the two copies of its MAT locus. The half-sibling species Z. parabailii and Z. pseudobailii therefore went through remarkably parallel but independent steps to regain fertility after they were formed by separate interspecies hybridizations.

Genome Sequence of Torulaspora microellipsoides CLIB 830T.

We report here the genome sequence of the ascomycetous yeast Torulaspora microellipsoides CLIB 830T A reference genome for this species, which has been found as a donor of genetic material in wine strains of Saccharomyces cerevisiae, will undoubtedly give clues to our understanding of horizontal transfer mechanisms between species in the wine environment.

Ploidy Variation in Kluyveromyces marxianus Separates Dairy and Non-dairy Isolates.

Kluyveromyces marxianus is traditionally associated with fermented dairy products, but can also be isolated from diverse non-dairy environments. Because of thermotolerance, rapid growth and other traits, many different strains are being developed for food and industrial applications but there is, as yet, little understanding of the genetic diversity or population genetics of this species. K. marxianus shows a high level of phenotypic variation but the only phenotype that has been clearly linked to a genetic polymorphism is lactose utilisation, which is controlled by variation in the LAC12 gene. The genomes of several strains have been sequenced in recent years and, in this study, we sequenced a further nine strains from different origins. Analysis of the Single Nucleotide Polymorphisms (SNPs) in 14 strains was carried out to examine genome structure and genetic diversity. SNP diversity in K. marxianus is relatively high, with up to 3% DNA sequence divergence between alleles. It was found that the isolates include haploid, diploid, and triploid strains, as shown by both SNP analysis and flow cytometry. Diploids and triploids contain long genomic tracts showing loss of heterozygosity (LOH). All six isolates from dairy environments were diploid or triploid, whereas 6 out 7 isolates from non-dairy environment were haploid. This also correlated with the presence of functional LAC12 alleles only in dairy haplotypes. The diploids were hybrids between a non-dairy and a dairy haplotype, whereas triploids included three copies of a dairy haplotype.

Transcriptional Response to Lactic Acid Stress in the Hybrid Yeast Zygosaccharomyces parabailii.

Lactic acid has a wide range of applications starting from its undissociated form, and its production using cell factories requires stress-tolerant microbial hosts. The interspecies hybrid yeast Zygosaccharomyces parabailii has great potential to be exploited as a novel host for lactic acid production, due to high organic acid tolerance at low pH and a fermentative metabolism with a high growth rate. Here we used mRNA sequencing (RNA-seq) to analyze Z. parabailii's transcriptional response to lactic acid added exogenously, and we explore the biological mechanisms involved in tolerance. Z. parabailii contains two homeologous copies of most genes. Under lactic acid stress, the two genes in each homeolog pair tend to diverge in expression to a significantly greater extent than under control conditions, indicating that stress tolerance is facilitated by interactions between the two gene sets in the hybrid. Lactic acid induces downregulation of genes related to cell wall and plasma membrane functions, possibly altering the rate of diffusion of lactic acid into cells. Genes related to iron transport and redox processes were upregulated, suggesting an important role for respiratory functions and oxidative stress defense. We found differences in the expression profiles of genes putatively regulated by Haa1 and Aft1/Aft2, previously described as lactic acid responsive in Saccharomyces cerevisiae Furthermore, formate dehydrogenase (FDH) genes form a lactic acid-responsive gene family that has been specifically amplified in Z. parabailii in comparison to other closely related species. Our study provides a useful starting point for the engineering of Z. parabailii as a host for lactic acid production.IMPORTANCE Hybrid yeasts are important in biotechnology because of their tolerance to harsh industrial conditions. The molecular mechanisms of tolerance can be studied by analyzing differential gene expression under conditions of interest and relating gene expression patterns to protein functions. However, hybrid organisms present a challenge to the standard use of mRNA sequencing (RNA-seq) to study transcriptional responses to stress, because their genomes contain two similar copies of almost every gene. Here we used stringent mapping methods and a high-quality genome sequence to study the transcriptional response to lactic acid stress in Zygosaccharomyces parabailii ATCC 60483, a natural interspecies hybrid yeast that contains two complete subgenomes that are approximately 7% divergent in sequence. Beyond the insights we gained into lactic acid tolerance in this study, the methods we developed will be broadly applicable to other yeast hybrid strains.