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

WHO elements - A new category of selfish genetic elements at the borderline between homing elements and transposable elements

Homing genetic elements are a form of selfish DNA that inserts into a specific target site in the genome and spreads through the population by a process of biased inheritance. Two well-known types of homing element, called inteins and homing introns, were discovered decades ago. In this review we describe WHO elements, a newly discovered type of homing element that constitutes a distinct third category but is rare, having been found only in a few yeast species so far. WHO elements are inferred to spread using the same molecular homing mechanism as inteins and introns: they encode a site-specific endonuclease that cleaves the genome at the target site, making a DNA break that is subsequently repaired by copying the element. For most WHO elements, the target site is in the glycolytic gene FBA1. WHO elements differ from inteins and homing introns in two fundamental ways: they do not interrupt their host gene (FBA1), and they occur in clusters. The clusters were formed by successive integrations of different WHO elements into the FBA1 locus, the result of an 'arms race' between the endonuclease and its target site. We also describe one family of WHO elements (WHO10) that is no longer specifically associated with the FBA1 locus and instead appears to have become transposable, inserting at random genomic sites in Torulaspora globosa with up to 26 copies per strain. The WHO family of elements is therefore at the borderline between homing genetic elements and transposable elements.

What are the 100 most cited fungal genera?

The global diversity of fungi has been estimated between 2 to 11 million species, of which only about 155 000 have been named. Most fungi are invisible to the unaided eye, but they represent a major component of biodiversity on our planet, and play essential ecological roles, supporting life as we know it. Although approximately 20 000 fungal genera are presently recognised, the ecology of most remains undetermined. Despite all this diversity, the mycological community actively researches some fungal genera more commonly than others. This poses an interesting question: why have some fungal genera impacted mycology and related fields more than others? To address this issue, we conducted a bibliometric analysis to identify the top 100 most cited fungal genera. A thorough database search of the Web of Science, Google Scholar, and PubMed was performed to establish which genera are most cited. The most cited 10 genera are Saccharomyces, Candida, Aspergillus, Fusarium, Penicillium, Trichoderma, Botrytis, Pichia, Cryptococcus and Alternaria. Case studies are presented for the 100 most cited genera with general background, notes on their ecology and economic significance and important research advances. This paper provides a historic overview of scientific research of these genera and the prospect for further research. Citation: Bhunjun CS, Chen YJ, Phukhamsakda C, Boekhout T, Groenewald JZ, McKenzie EHC, Francisco EC, Frisvad JC, Groenewald M, Hurdeal VG, Luangsa-ard J, Perrone G, Visagie CM, Bai FY, Błaszkowski J, Braun U, de Souza FA, de Queiroz MB, Dutta AK, Gonkhom D, Goto BT, Guarnaccia V, Hagen F, Houbraken J, Lachance MA, Li JJ, Luo KY, Magurno F, Mongkolsamrit S, Robert V, Roy N, Tibpromma S, Wanasinghe DN, Wang DQ, Wei DP, Zhao CL, Aiphuk W, Ajayi-Oyetunde O, Arantes TD, Araujo JC, Begerow D, Bakhshi M, Barbosa RN, Behrens FH, Bensch K, Bezerra JDP, Bilański P, Bradley CA, Bubner B, Burgess TI, Buyck B, Čadež N, Cai L, Calaça FJS, Campbell LJ, Chaverri P, Chen YY, Chethana KWT, Coetzee B, Costa MM, Chen Q, Custódio FA, Dai YC, Damm U, de Azevedo Santiago ALCM, De Miccolis Angelini RM, Dijksterhuis J, Dissanayake AJ, Doilom M, Dong W, Alvarez-Duarte E, Fischer M, Gajanayake AJ, Gené J, Gomdola D, Gomes AAM, Hausner G, He MQ, Hou L, Iturrieta-González I, Jami F, Jankowiak R, Jayawardena RS, Kandemir H, Kiss L, Kobmoo N, Kowalski T, Landi L, Lin CG, Liu JK, Liu XB, Loizides M, Luangharn T, Maharachchikumbura SSN, Makhathini Mkhwanazi GJ, Manawasinghe IS, Marin-Felix Y, McTaggart AR, Moreau PA, Morozova OV, Mostert L, Osiewacz HD, Pem D, Phookamsak R, Pollastro S, Pordel A, Poyntner C, Phillips AJL, Phonemany M, Promputtha I, Rathnayaka AR, Rodrigues AM, Romanazzi G, Rothmann L, Salgado-Salazar C, Sandoval-Denis M, Saupe SJ, Scholler M, Scott P, Shivas RG, Silar P, Souza-Motta CM, Silva-Filho AGS, Spies CFJ, Stchigel AM, Sterflinger K, Summerbell RC, Svetasheva TY, Takamatsu S, Theelen B, Theodoro RC, Thines M, Thongklang N, Torres R, Turchetti B, van den Brule T, Wang XW, Wartchow F, Welti S, Wijesinghe SN, Wu F, Xu R, Yang ZL, Yilmaz N, Yurkov A, Zhao L, Zhao RL, Zhou N, Hyde KD, Crous PW (2024). What are the 100 most cited fungal genera? Studies in Mycology 108: 1-411. doi: 10.3114/sim.2024.108.01.

Genomic factors shape carbon and nitrogen metabolic niche breadth across Saccharomycotina yeasts

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Two general paradigms have been proposed to explain this variation: (i) trade-offs between performance efficiency and breadth and (ii) the joint influence of extrinsic (environmental) and intrinsic (genomic) factors. We assembled genomic, metabolic, and ecological data from nearly all known species of the ancient fungal subphylum Saccharomycotina (1154 yeast strains from 1051 species), grown in 24 different environmental conditions, to examine niche breadth evolution. We found that large differences in the breadth of carbon utilization traits between yeasts stem from intrinsic differences in genes encoding specific metabolic pathways, but we found limited evidence for trade-offs. These comprehensive data argue that intrinsic factors shape niche breadth variation in microbes.

Unlocking the genome of the non-sourdough Kazachstania humilis MAW1: insights into inhibitory factors and phenotypic properties

BACKGROUND

Ascomycetous budding yeasts are ubiquitous environmental microorganisms important in food production and medicine. Due to recent intensive genomic research, the taxonomy of yeast is becoming more organized based on the identification of monophyletic taxa. This includes genera important to humans, such as Kazachstania. Until now, Kazachstania humilis (previously Candida humilis) was regarded as a sourdough-specific yeast. In addition, any antibacterial activity has not been associated with this species.

RESULTS

Previously, we isolated a yeast strain that impaired bio-hydrogen production in a dark fermentation bioreactor and inhibited the growth of Gram-positive and Gram-negative bacteria. Here, using next generation sequencing technologies, we sequenced the genome of this strain named K. humilis MAW1. This is the first genome of a K. humilis isolate not originating from a fermented food. We used novel phylogenetic approach employing the 18 S-ITS-D1-D2 region to show the placement of the K. humilis MAW1 among other members of the Kazachstania genus. This strain was examined by global phenotypic profiling, including carbon sources utilized and the influence of stress conditions on growth. Using the well-recognized bacterial model Escherichia coli AB1157, we show that K. humilis MAW1 cultivated in an acidic medium inhibits bacterial growth by the disturbance of cell division, manifested by filament formation. To gain a greater understanding of the inhibitory effect of K. humilis MAW1, we selected 23 yeast proteins with recognized toxic activity against bacteria and used them for Blast searches of the K. humilis MAW1 genome assembly. The resulting panel of genes present in the K. humilis MAW1 genome included those encoding the 1,3-β-glucan glycosidase and the 1,3-β-glucan synthesis inhibitor that might disturb the bacterial cell envelope structures.

CONCLUSIONS

We characterized a non-sourdough-derived strain of K. humilis, including its genome sequence and physiological aspects. The MAW1, together with other K. humilis strains, shows the new organization of the mating-type locus. The revealed here pH-dependent ability to inhibit bacterial growth has not been previously recognized in this species. Our study contributes to the building of genome sequence-based classification systems; better understanding of K.humilis as a cell factory in fermentation processes and exploring bacteria-yeast interactions in microbial communities.

The metabolic domestication syndrome of budding yeast

Cellular metabolism evolves through changes in the structure and quantitative states of metabolic networks. Here, we explore the evolutionary dynamics of metabolic states by focusing on the collection of metabolite levels, the metabolome, which captures key aspects of cellular physiology. Using a phylogenetic framework, we profiled metabolites in 27 populations of nine budding yeast species, providing a graduated view of metabolic variation across multiple evolutionary time scales. Metabolite levels evolve more rapidly and independently of changes in the metabolic network's structure, providing complementary information to enzyme repertoire. Although metabolome variation accumulates mainly gradually over time, it is profoundly affected by domestication. We found pervasive signatures of convergent evolution in the metabolomes of independently domesticated clades of Saccharomyces cerevisiae. Such recurring metabolite differences between wild and domesticated populations affect a substantial part of the metabolome, including rewiring of the TCA cycle and several amino acids that influence aroma production, likely reflecting adaptation to human niches. Overall, our work reveals previously unrecognized diversity in central metabolism and the pervasive influence of human-driven selection on metabolite levels in yeasts.

High quality de novo genome assembly of the non-conventional yeast Kazachstania bulderi describes a potential low pH production host for biorefineries

Kazachstania bulderi is a non-conventional yeast species able to grow efficiently on glucose and δ-gluconolactone at low pH. These unique traits make K. bulderi an ideal candidate for use in sustainable biotechnology processes including low pH fermentations and the production of green chemicals including organic acids. To accelerate strain development with this species, detailed information of its genetics is needed. Here, by employing long read sequencing we report a high-quality phased genome assembly for three strains of K. bulderi species, including the type strain. The sequences were assembled into 12 chromosomes with a total length of 14 Mb, and the genome was fully annotated at structural and functional levels, including allelic and structural variants, ribosomal array and mating type locus. This high-quality reference genome provides a resource to advance our fundamental knowledge of biotechnologically relevant non-conventional yeasts and to support the development of genetic tools for manipulating such strains towards their use as production hosts in biotechnological processes.

POInTbrowse: orthology prediction and synteny exploration for paleopolyploid genomes

We describe POInT_browse, a web portal that gives access to the orthology inferences made for polyploid genomes with POInT, the Polyploidy Orthology Inference Tool. Ancient, or paleo-, polyploidy events are widely distributed across the eukaryotic phylogeny, and the combination of duplicated and lost duplicated genes that these polyploidies produce can confound the identification of orthologous genes between genomes. POInT uses conserved synteny and phylogenetic models to infer orthologous genes between genomes with a shared polyploidy. It also gives confidence estimates for those orthology inferences. POInT_browse gives both graphical and query-based access to these inferences from 12 different polyploidy events, allowing users to visualize genomic regions produced by polyploidies and perform batch queries for each polyploidy event, downloading genes trees and coding sequences for orthologous genes meeting user-specified criteria. POInT_browse and the associated data are online at https://wgd.statgen.ncsu.edu .

Parallel expansion and divergence of an adhesin family in pathogenic yeasts

Opportunistic yeast pathogens arose multiple times in the Saccharomycetes class, including the recently emerged, multidrug-resistant (MDR) Candida auris. We show that homologs of a known yeast adhesin family in Candida albicans, the Hyr/Iff-like (Hil) family, are enriched in distinct clades of Candida species as a result of multiple, independent expansions. Following gene duplication, the tandem repeat-rich region in these proteins diverged extremely rapidly and generated large variations in length and β-aggregation potential, both of which are known to directly affect adhesion. The conserved N-terminal effector domain was predicted to adopt a β-helical fold followed by an α-crystallin domain, making it structurally similar to a group of unrelated bacterial adhesins. Evolutionary analyses of the effector domain in C. auris revealed relaxed selective constraint combined with signatures of positive selection, suggesting functional diversification after gene duplication. Lastly, we found the Hil family genes to be enriched at chromosomal ends, which likely contributed to their expansion via ectopic recombination and break-induced replication. Combined, these results suggest that the expansion and diversification of adhesin families generate variation in adhesion and virulence within and between species and are a key step toward the emergence of fungal pathogens.

Isolation and characterization of Zygosaccharomyces sp. yeast strains from miso

There is currently great interest in the salt-tolerant yeast strains used to produce miso and soy sauce. Since the isolation of Zygosaccharomyces sp. strain from Japanese miso more than 60 years, several hybrid strains have been identified in fermented foods. Studies have shown that the active mating-type locus of the original Zygosaccharomyces sp. yeast strain is located between the T-subgenome sequence and the P-subgenome sequence. In this study, 32 salt-tolerant Zygosaccharomyces sp. yeast strains were isolated from five miso factories in Hiroshima Prefecture, Japan. Analysis by flow cytometry revealed that 27 strains were diploid and five strains were haploid. PCR analysis indicated that the 27 diploid strains had the same chromosomal structure of the active mating-type (MAT) locus as the original yeast strain isolated from miso 60 years ago. In addition, the 27 diploid strains were allodiploid, namely, natural hybrids of Z. rouxii and a related species, while the five haploid strains were all Z. rouxii. We found that cells of yeast strains isolated from miso changed morphologically when co-cultured with a yeast strain of opposite mating-type under nitrogen starvation conditions. The DNA sequence of the active mating-type locus and the results of cell morphology changes by co-culture were consistent with the mating type of each strain shown in the mating experiments. These findings will be useful for the future production of miso and soy sauce.

The DNA replication protein Orc1 from the yeast Torulaspora delbrueckii is required for heterochromatin formation but not as a silencer-binding protein

To understand the process by which new protein functions emerge, we examined how the yeast heterochromatin protein Sir3 arose through gene duplication from the conserved DNA replication protein Orc1. Orc1 is a subunit of the origin recognition complex (ORC), which marks origins of DNA replication. In Saccharomyces cerevisiae, Orc1 also promotes heterochromatin assembly by recruiting the structural proteins Sir1-4 to silencer DNA. In contrast, the paralog of Orc1, Sir3, is a nucleosome-binding protein that spreads across heterochromatic loci in conjunction with other Sir proteins. We previously found that a non-duplicated Orc1 from the yeast Kluyveromyces lactis behaved like ScSir3 but did not have a silencer-binding function like ScOrc1. Moreover, K. lactis lacks Sir1, the protein that interacts directly with ScOrc1 at the silencer. Here, we examined whether the emergence of Sir1 coincided with Orc1 acting as a silencer-binding protein. In the non-duplicated species Torulaspora delbrueckii, which has an ortholog of Sir1 (TdKos3), we found that TdOrc1 spreads across heterochromatic loci independently of ORC, as ScSir3 and KlOrc1 do. This spreading is dependent on the nucleosome binding BAH domain of Orc1 and on Sir2 and Kos3. However, TdOrc1 does not have a silencer-binding function: T. delbrueckii silencers do not require ORC binding sites to function, and Orc1 and Kos3 do not appear to interact. Instead, Orc1 and Kos3 both spread across heterochromatic loci with other Sir proteins. Thus, Orc1 and Sir1/Kos3 originally had different roles in heterochromatin formation than they do now in S. cerevisiae.

Mating-Type Switching in Budding Yeasts, from Flip/Flop Inversion to Cassette Mechanisms

Mating-type switching is a natural but unusual genetic control process that regulates cell identity in ascomycete yeasts. It involves physically replacing one small piece of genomic DNA by another, resulting in replacement of the master regulatory genes in the mating pathway and hence a switch of cell type and mating behavior. In this review, we concentrate on recent progress that has been made on understanding the origins and evolution of mating-type switching systems in budding yeasts (subphylum Saccharomycotina). Because of the unusual nature and the complexity of the mechanism in Saccharomyces cerevisiae, mating-type switching was assumed until recently to have originated only once or twice during yeast evolution. However, comparative genomics analysis now shows that switching mechanisms arose many times independently-at least 11 times in budding yeasts and once in fission yeasts-a dramatic example of convergent evolution. Most of these lineages switch mating types by a flip/flop mechanism that inverts a section of a chromosome and is simpler than the well-characterized 3-locus cassette mechanism (MAT/HML/HMR) used by S. cerevisiae. Mating-type switching (secondary homothallism) is one of the two possible mechanisms by which a yeast species can become self-fertile. The other mechanism (primary homothallism) has also emerged independently in multiple evolutionary lineages of budding yeasts, indicating that homothallism has been favored strongly by natural selection. Recent work shows that HO endonuclease, which makes the double-strand DNA break that initiates switching at the S. cerevisiae MAT locus, evolved from an unusual mobile genetic element that originally targeted a glycolytic gene, FBA1.

Functional divergence in the proteins encoded by ARO80 from S. uvarum, S. kudriavzevii and S. cerevisiae explain differences in the aroma production during wine fermentation

Phenylethanol (PE) and phenylethyl acetate (PEA) are commonly desired compounds in wine because of their rose‐like aroma. The yeast S. cerevisiae produces the PE either through de novo biosynthesis by shikimate pathway followed by the Ehrlich pathway or the direct phenylalanine catabolism via Ehrlich pathway, and then converted into PEA. Previous work demonstrated that, compared to S. cerevisiae, other Saccharomyces species, such as S. kudriavzevii and S. uvarum, produce higher concentrations of PE and PEA from the precursor phenylalanine, which indicates differential activities of the biosynthetic‐involved enzymes. A previous in‐silico analysis suggested that the transcriptional activator Aro80p is one of the best candidates to explain these differences. An improved functional analysis identified significant radical amino acid changes in the S. uvarum and S. kudriavzevii Aro80p that could impact the expression of the catabolic genes ARO9 and ARO10, and hence, the production of PE from phenylalanine. Indeed, wine S. cerevisiae strains carrying the S. uvarum and S. kudriavzevii ARO80 alleles increased the production of both compounds in the presence of phenylalanine by increasing the expression of ARO9 and ARO10. This study provides novel insights of the unidentified Aro80p regulatory region and the potential usage of alternatives ARO80 alleles to enhance the PE and PEA concentration in wine.

OUP accepted manuscript

Recent studies have suggested that species of the Kazachstania genus may be interesting models of yeast domestication. Among these, Kazachstania barnettii has been isolated from various microbially transformed foodstuffs such as sourdough bread and kefir. In the present work, we sequence, assemble, and annotate the complete genomes of two K. barnettii strains: CLIB 433, being one of the two reference strains for K. barnettii that was isolated as a spoilage organism in soft drink, and CLIB 1767, recently isolated from artisan bread-making sourdough. Both assemblies are of high quality with N50 statistics greater than 1.3 Mb and BUSCO score greater than 99%. An extensive comparison of the two obtained genomes revealed very few differences between the two K. barnettii strains, considering both genome structure and gene content. The proposed genome assemblies will constitute valuable references for future comparative genomic, population genomic, or transcriptomic studies of the K. barnettii species.

Insights on life cycle and cell identity regulatory circuits for unlocking genetic improvement in Zygosaccharomyces and Kluyveromyces yeasts

Evolution has provided a vast diversity of yeasts that play fundamental roles in nature and society. This diversity is not limited to genotypically homogeneous species with natural interspecies hybrids and allodiploids that blur species boundaries frequently isolated. Thus, life cycle and the nature of breeding systems have profound effects on genome variation, shaping heterozygosity, genotype diversity and ploidy level. The apparent enrichment of hybrids in industry-related environments suggests that hybridization provides an adaptive route against stressors and creates interest in developing new hybrids for biotechnological uses. For example, in the Saccharomyces genus where regulatory circuits controlling cell identity, mating competence and meiosis commitment have been extensively studied, this body of knowledge is being used to combine interesting traits into synthetic F1 hybrids, to bypass F1 hybrid sterility and to dissect complex phenotypes by bulk segregant analysis. Although these aspects are less known in other industrially promising yeasts, advances in whole-genome sequencing and analysis are changing this and new insights are being gained, especially in the food-associated genera Zygosaccharomyces and Kluyveromyces. We discuss this new knowledge and highlight how deciphering cell identity circuits in these lineages will contribute significantly to identify the genetic determinants underpinning complex phenotypes and open new avenues for breeding programmes.

Genomic analysis revealed conserved acid tolerance mechanisms from native micro-organisms in fermented feed

AIMS

Fermented feed is an agricultural practice used in many regions of the world to improve the growth performance of farm animals. This study aimed to identify and evaluate the lactic acid bacteria and yeast involved in the production of fermented feed.

METHODS AND RESULTS

We isolated and described two micro-organisms from autochthonous microbiota origin present in a regional feed product, Lactobacillus paracasei IBR07 (Lacticaseibacillus paracasei) and Kazachstania unispora IBR014 (Saccharomyces unisporum). Genome sequence analyses were performed to characterize both micro-organisms. Potential pathways involved in the acid response, tolerance and persistence were predicted in both genomes. Although L. paracasei and K. unispora are considered safe for animal feed, we analysed the presence of virulence factors, antibiotic resistance and pathogenicity islands. Furthermore, the Galleria mellonella model was used to support the safety of both isolates.

CONCLUSIONS

We conclude that IBR07 and IBR014 strains are good candidates to be used as starter cultures for feed fermentation.

SIGNIFICANCE AND IMPACT OF THE STUDY

The data presented here will be helpful to explore other biotechnological aspects and constitute a starting point for further studies to establish the consumption benefit of fermented feed in farm animal production.

Biotechnological Importance of Torulaspora delbrueckii: From the Obscurity to the Spotlight

Torulaspora delbrueckii has attracted interest in recent years, especially due to its biotechnological potential, arising from its flavor- and aroma-enhancing properties when used in wine, beer or bread dough fermentation, as well as from its remarkable resistance to osmotic and freezing stresses. In the present review, genomic, biochemical, and phenotypic features of T. delbrueckii are described, comparing them with other species, particularly with the biotechnologically well-established yeast, Saccharomyces cerevisiae. We conclude about the aspects that make this yeast a promising biotechnological model to be exploited in a wide range of industries, particularly in wine and bakery. A phylogenetic analysis was also performed, using the core proteome of T. delbrueckii, to compare the number of homologous proteins relative to the most closely related species, understanding the phylogenetic placement of this species with robust support. Lastly, the genetic tools available for T. delbrueckii improvement are discussed, focusing on adaptive laboratorial evolution and its potential.

Conservation of a DNA Replication Motif among Phylogenetically Distant Budding Yeast Species

Eukaryotic DNA replication begins at genomic loci termed origins, which are bound by the origin recognition complex (ORC). Although ORC is conserved across species, the sequence composition of origins is more varied. In the budding yeast Saccharomyces cerevisiae, the ORC-binding motif consists of an A/T-rich 17 bp “extended ACS” sequence adjacent to a B1 element composed of two 3-bp motifs. Similar sequences occur at origins in closely related species, but it is not clear when this type of replication origin arose and whether it predated a whole-genome duplication that occurred around 100 Ma in the budding yeast lineage. To address these questions, we identified the ORC-binding sequences in the nonduplicated species Torulaspora delbrueckii. We used chromatin immunoprecipitation followed by sequencing and identified 190 ORC-binding sites distributed across the eight T. delbrueckii chromosomes. Using these sites, we identified an ORC-binding motif that is nearly identical to the known motif in S. cerevisiae. We also found that the T. delbrueckii ORC-binding sites function as origins in T. delbrueckii when cloned onto a plasmid and that the motif is required for plasmid replication. Finally, we compared an S. cerevisiae origin with two T. delbrueckii ORC-binding sites and found that they conferred similar stabilities to a plasmid. These results reveal that the ORC-binding motif arose prior to the whole-genome duplication and has been maintained for over 100 Myr.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

Hexose transport in Torulaspora delbrueckii: identification of Igt1, a new dual-affinity transporter

Torulaspora delbrueckii is a yeast species receiving increasing attention from the biotechnology industry, with particular relevance in the wine, beer and baking sectors. However, little is known about its sugar transporters and sugar transport capacity, frequently a rate-limiting step of sugar metabolism and efficient fermentation. Actually, only one glucose transporter, Lgt1, has been characterized so far. Here we report the identification and characterization of a second glucose transporter gene, IGT1, located in a cluster, upstream of LGT1 and downstream of two other putative hexose transporters. Functional characterization of IGT1 in a Saccharomyces cerevisiae hxt-null strain revealed that it encodes a transporter able to mediate uptake of glucose, fructose and mannose and established that its affinity, as measured by Km, could be modulated by glucose concentration in the medium. In fact, IGT1-transformed S. cerevisiae hxt-null cells, grown in 0.1% glucose displayed biphasic glucose uptake kinetics with an intermediate- (Km = 6.5 ± 2.0 mM) and a high-affinity (Km = 0.10 ± 0.01 mM) component, whereas cells grown in 2% glucose displayed monophasic kinetics with an intermediate-affinity (Km of 11.5 ± 1.5 mM). This work contributes to a better characterization of glucose transport in T. delbrueckii, with relevant implications for its exploitation in the food industry.

Improvement of Torulaspora delbrueckii Genome Annotation: Towards the Exploitation of Genomic Features of a Biotechnologically Relevant Yeast

Saccharomyces cerevisiae is the most commonly used yeast in wine, beer, and bread fermentations. However, Torulaspora delbrueckii has attracted interest in recent years due to its properties, ranging from its ability to produce flavor- and aroma-enhanced wine to its ability to survive longer in frozen dough. In this work, publicly available genomes of T. delbrueckii were explored and their annotation was improved. A total of 32 proteins were additionally annotated for the first time in the type strain CBS1146, in comparison with the previous annotation available. In addition, the annotation of the remaining three T. delbrueckii strains was performed for the first time. eggNOG-mapper was used to perform the functional annotation of the deduced T. delbrueckii coding genes, offering insights into its biological significance, and revealing 24 clusters of orthologous groups (COGs), which were gathered in three main functional categories: information storage and processing (28% of the proteins), cellular processing and signaling (27%), and metabolism (23%). Small intraspecies variability was found when considering the functional annotation of the four available T. delbrueckii genomes. A comparative study was also conducted between the T. delbrueckii genome and those from 386 fungal species, revealing a high number of homologous genes with species from the Zygotorulaspora and Zygosaccharomyces genera, but also with Lachancea and S. cerevisiae. Lastly, the phylogenetic placement of T. delbrueckii was clarified using the core homologs that were found across 204 common protein sequences of 386 fungal species and strains.

Regulation of mating and mating-type-specific genes in Zygosaccharomyces sp. yeast

Zygosaccharomyces sp. is an industrially important yeast for the production traditional fermented foods in Japan. At present, however, there is no easy method for mating Zygosaccharomyces sp. strains in the laboratory; furthermore, little is known about the expression of mating-type-specific genes in this yeast. Here, mating was observed when Zygosaccharomyces sp. was subjected to nitrogen-starvation conditions. The expression of mating-type-specific genes, Zygo STE6 and Zygo MFα1, was induced under nitrogen-starvation conditions, as confirmed by lacZ reporter assay. This expression was mating-type-specific: Zygo STE6 was expressed specifically for mating-type a, whereas and Zygo MFα1 was expressed specifically for mating-type α. Yeast strains Zygosaccharomyces rouxii DL2 and DA2, derived from type strain Z. rouxii CBS732, did not show mating even under nitrogen-starvation conditions. Gene sequencing revealed that the Zygo STE12 in Z. rouxii CBS732 has a frameshift mutation. Under nitrogen starvation, mating was observed in both DL2 and DA2 transformed with the wild-type Zygo STE12. The expression of Zygo STE6 in Z. rouxii DL2 transformed with wild-type Zygo STE12 under nitrogen-starvation conditions was confirmed by lacZ reporter assay. Collectively, these results revealed that, under nitrogen-starvation conditions, Zygosaccharomyces sp. can mate and mating-type-specific genes are expressed. Furthermore, Zygo Ste12 is essential for both mating and the expression of mating-type-specific genes in Zygosaccharomyces sp.

Identification of a novel gene controlling homothallism in the yeast Kazachstania naganishii isolated in Japan

A novel gene controlling homothallic life cycle was identified in the yeast Kazachstania naganishii isolated in Japan. This gene was isolated by means of complementing a mutation, mti1, which had led to heterothallism from original homothallism in the yeast. The configuration of original mutation in MTI1 gene revealed that a truncated product is formed due to occurrence of a stop codon by a nucleotide insertion. When the gene was disrupted with a marker, the disruptant spore clone was haploid and stably heterothallic. Disfunction of the gene caused inability to self-diploidize due to defect of mating-type interconversion. The gene MTI1 (for Mating Type Interconversion) is a weak homolog of the Saccharomyces cerevisiae VID22/ENV11, which has been reported to function in vacuolar protein processing. K. naganishii has a gene representing significant homology with the HO gene of S. cerevisiae on chromosome V, which has not been clarified to be involved in regulation of life cycle in K. naganishii. The MTI1 gene defined in this study is located on K. naganishii chromosome IV and does not represent significant homology to the above ScHO-like gene and any other genes concerning life cycles of yeasts. From the viewpoint of gene evolution, it is extremely interesting that the MTI1 gene is a new type of gene controlling homothallism in addition to an HO-type gene, leading to discovery of an unknown mechanism regulating life cycles in yeasts.

Genome editing in the yeast Nakaseomyces delphensis and description of its complete sexual cycle

The environmental yeast Nakaseomyces delphensis is, phylogenetically, the closest known species to Candida glabrata, a major fungal pathogen of humans. C. glabrata is haploid and described as asexual, while N. delphensis is also haploid, but has been described as competent for mating and meiosis. Both genomes contain homologues of all the genes necessary for sexual reproduction and also the genes for Ho-dependent mating-type switching, like Saccharomyces cerevisiae. We first report the construction of genetically engineered strains of N. delphensis, including by CRISPR-Cas 9 gene editing. We also report the description of the sexual cycle of N. delphensis. We show that it undergoes Ho-dependent mating-type switching in culture and that deletion of the HO gene prevents such switching and allows maintenance of stable, separate, MATa and MATalpha haploid strains. Rare, genetically selected diploids can be obtained through mating of haploid strains, mutated or not for the HO gene. In contrast to HO/HO diploids, which behave as expected, Δho/Δho diploids exhibit unusual profiles in flow cytometry. Both types of diploids can produce recombined haploid cells, which grow like the original haploid-type strain. Our experiments thus allow the genetic manipulation of N. delphensis and the reconstruction, in the laboratory, of its entire life cycle.

Kluyveromyces osmophilus is not a synonym of Zygosaccharomyces mellis; reinstatement as Zygosaccharomyces osmophilus comb. nov

Kluyveromyces osmophilus, a single-strain species isolated from Mozambique sugar, has been treated a synonym of Zygosaccharomyces mellis. Analyses of D1/D2 LSU rRNA gene sequences confirmed that the species belongs to the genus Zygosaccharomyces but showed it to be distinct from strains of Z. mellis. During studies of yeasts associated with stingless bees in Brazil, nine additional isolates of the species were obtained from unripe and ripe honey and pollen of Scaptotrigona cfr. bipunctata, as well as ripe honey of Tetragonisca angustula. The D1/D2 sequences of the Brazilian isolates were identical to those of the type strain of K. osmophilus CBS 5499 (=ATCC 22027), indicating that they represent the same species. Phylogenomic analyses using 4038 orthologous genes support the reinstatement of K. osmophilus as a member of the genus Zygosaccharomyces. We, therefore, propose the name Zygosaccharomyces osmophilus comb. nov. (lectotype ATCC 22027; MycoBank no. MB 833739).

The yeast mating-type switching endonuclease HO is a domesticated member of an unorthodox homing genetic element family

The mating-type switching endonuclease HO plays a central role in the natural life cycle of Saccharomyces cerevisiae, but its evolutionary origin is unknown. HO is a recent addition to yeast genomes, present in only a few genera close to Saccharomyces. Here we show that HO is structurally and phylogenetically related to a family of unorthodox homing genetic elements found in Torulaspora and Lachancea yeasts. These WHO elements home into the aldolase gene FBA1, replacing its 3' end each time they integrate. They resemble inteins but they operate by a different mechanism that does not require protein splicing. We show that a WHO protein cleaves Torulaspora delbrueckii FBA1 efficiently and in an allele-specific manner, leading to DNA repair by gene conversion or NHEJ. The DNA rearrangement steps during WHO element homing are very similar to those during mating-type switching, and indicate that HO is a domesticated WHO-like element.

In the same way as a sperm from a male and an egg from a female join together to form an embryo in most animals, yeast cells have two sexes that coordinate how they reproduce. These are called “mating types” and, rather than male or female, an individual yeast cell can either be mating type “a” or “alpha”. Every yeast cell contains the genes for both mating types, and each cell’s mating type is determined by which of those genes it has active.

Only one mating type gene can be ‘on’ at a time, but some yeast species can swap mating type on demand by switching the corresponding genes ‘on’ or ‘off’. This switch is unusual. Rather than simply activate one of the genes it already has, the yeast cell keeps an inactive version of each mating type gene tucked away, makes a copy of the gene it wants to be active and pastes that copy into a different location in its genome. To do all of this yeast need another gene called HO. This gene codes for an enzyme that cuts the DNA at the location of the active mating type gene. This makes an opening that allows the cell to replace the ‘a’ gene with the ‘alpha’ gene, or vice versa. This system allows yeast cells to continue mating even if all the cells in a colony start off as the same mating type. But, cutting into the DNA is risky, and can damage the health of the cell. So, why did yeast cells evolve a system that could cause them harm?

To find out where the HO gene came from, Coughlan et al. searched through all the available genomes from yeast species for other genes with similar sequences and identified a cluster which they nicknamed “weird HO” genes, or WHO genes for short. Testing these genes revealed that they also code for enzymes that make cuts in the yeast genome, but the way the cell repairs the cuts is different. The WHO genes are jumping genes. When the enzyme encoded by a WHO gene makes a cut in the genome, the yeast cell copies the gene into the gap, allowing the gene to ‘jump’ from one part of the genome to another. It is possible that this was the starting point for the evolution of the HO gene. Changes to a WHO gene could have allowed it to cut into the mating type region of the yeast genome, giving the yeast an opportunity to ‘domesticate’ it. Over time, the yeast cell stopped the WHO gene from jumping into the gap and started using the cut to change its mating type.

Understanding how cells adapt genes for different purposes is a key question in evolutionary biology. There are many other examples of domesticated jumping genes in other organisms, including in the human immune system. Understanding the evolution of HO not only sheds light on how yeast mating type switching evolved, but on how other species might harness and adapt their genes.

Alternative Lengthening of Telomeres in the Budding Yeast Naumovozyma castellii

The enzyme telomerase ensures the integrity of linear chromosomes by maintaining telomere length. As a hallmark of cancer, cell immortalization and unlimited proliferation is gained by reactivation of telomerase. However, a significant fraction of cancer cells instead uses alternative telomere lengthening mechanisms to ensure telomere function, collectively known as Alternative Lengthening of Telomeres (ALT). Although the budding yeast Naumovozyma castellii (Saccharomyces castellii) has a proficient telomerase activity, we demonstrate here that telomeres in N. castellii are efficiently maintained by a novel ALT mechanism after telomerase knockout. Remarkably, telomerase-negative cells proliferate indefinitely without any major growth crisis and display wild-type colony morphology. Moreover, ALT cells maintain linear chromosomes and preserve a wild-type DNA organization at the chromosome termini, including a short stretch of terminal telomeric sequence. Notably, ALT telomeres are elongated by the addition of ∼275 bp repeats containing a short telomeric sequence and the subtelomeric DNA located just internally (TelKO element). Although telomeres may be elongated by several TelKO repeats, no dramatic genome-wide amplification occurs, thus indicating that the repeat addition may be regulated. Intriguingly, a short interstitial telomeric sequence (ITS) functions as the initiation point for the addition of the TelKO element. This implies that N. castellii telomeres are structurally predisposed to efficiently switch to the ALT mechanism as a response to telomerase dysfunction.

Multiple Reinventions of Mating-type Switching during Budding Yeast Evolution

Cell type in budding yeasts is determined by the genotype at the mating-type (MAT) locus, but yeast species differ widely in their mating compatibility systems and life cycles. Among sexual yeasts, heterothallic species are those in which haploid strains fall into two distinct and stable mating types (MATa and MATα), whereas homothallic species are those that can switch mating types or that appear not to have distinct mating types [1, 2]. The evolutionary history of these mating compatibility systems is uncertain, particularly regarding the number and direction of transitions between homothallism and heterothallism, and regarding whether the process of mating-type switching had a single origin [3, 4, 5]. Here, we inferred the mating compatibility systems of 332 budding yeast species from their genome sequences. By reference to a robust phylogenomic tree [6], we detected evolutionary transitions between heterothallism and homothallism, and among different forms of homothallism. We find that mating-type switching has arisen independently at least 11 times during yeast evolution and that transitions from heterothallism to homothallism greatly outnumber transitions in the opposite direction (31 versus 3). Although the 3-locus MAT-HML-HMR mechanism of mating-type switching as seen in Saccharomyces cerevisiae had a single evolutionary origin in budding yeasts, simpler “flip/flop” mechanisms of switching evolved separately in at least 10 other groups of yeasts. These results point to the adaptive value of homothallism and mating-type switching to unicellular fungi.

•

Mating-type switching by flip-flopping arose at least 10 separate times

•

Mating-type switching by copy-and-paste arose only once in budding yeasts

•

Transitions toward homothallism outnumber transitions away from homothallism

•

Heterothallic species can become homothallic by DNA introgression

Positive Selection Evidence in Xylose-Related Genes Suggests Methylglyoxal Reductase as a Target for the Improvement of Yeasts' Fermentation in Industry

Xylose assimilation and fermentation are important traits for second generation ethanol production. However, some genomic features associated with this pentose sugar's metabolism remain unknown in yeasts. Comparative genomics studies have led to important insights in this field, but we are still far from completely understanding endogenous yeasts' xylose metabolism. In this work, we carried out a deep evolutionary analysis suited for comparative genomics of xylose-consuming yeasts, searching for of positive selection on genes associated with glucose and xylose metabolism in the xylose-fermenters' clade. Our investigation detected positive selection fingerprints at this clade not only among sequences of important genes for xylose metabolism, such as xylose reductase and xylitol dehydrogenase, but also in genes expected to undergo neutral evolution, such as the glycolytic gene phosphoglycerate mutase. In addition, we present expansion, positive selection marks, and convergence as evidence supporting the hypothesis that natural selection is shaping the evolution of the little studied methylglyoxal reductases. We propose a metabolic model suggesting that selected codons among these proteins caused a putative change in cofactor preference from NADPH to NADH that alleviates cellular redox imbalance. These findings provide a wider look into pentose metabolism of yeasts and add this previously overlooked piece into the intricate puzzle of oxidative imbalance. Although being extensively discussed in evolutionary works the awareness of selection patterns is recent in biotechnology researches, rendering insights to surpass the reached status quo in many of its subareas.

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.

Discovery and Evolution of New Domains in Yeast Heterochromatin Factor Sir4 and Its Partner Esc1

Sir4 is a core component of heterochromatin found in yeasts of the Saccharomycetaceae family, whose general hallmark is to harbor a three-loci mating-type system with two silent loci. However, a large part of the Sir4 amino acid sequences has remained unexplored, belonging to the dark proteome. Here, we analyzed the phylogenetic profile of yet undescribed foldable regions present in Sir4 as well as in Esc1, an Sir4-interacting perinuclear anchoring protein. Within Sir4, we identified a new conserved motif (TOC) adjacent to the N-terminal KU-binding motif. We also found that the Esc1-interacting region of Sir4 is a Dbf4-related H-BRCT domain, only present in species possessing the HO endonuclease and in Kluveryomyces lactis. In addition, we found new motifs within Esc1 including a motif (Esc1-F) that is unique to species where Sir4 possesses an H-BRCT domain. Mutagenesis of conserved amino acids of the Sir4 H-BRCT domain, known to play a critical role in the Dbf4 function, shows that the function of this domain is separable from the essential role of Sir4 in transcriptional silencing and the protection from HO-induced cutting in Saccharomyces cerevisiae. In the more distant methylotrophic clade of yeasts, which often harbor a two-loci mating-type system with one silent locus, we also found a yet undescribed H-BRCT domain in a distinct protein, the ISWI2 chromatin-remodeling factor subunit Itc1. This study provides new insights on yeast heterochromatin evolution and emphasizes the interest of using sensitive methods of sequence analysis for identifying hitherto ignored functional regions within the dark proteome.

Comparative Genomics Between Saccharomyces kudriavzevii and S. cerevisiae Applied to Identify Mechanisms Involved in Adaptation

Yeasts belonging to the Saccharomyces genus play an important role in human-driven fermentations. The species S. cerevisiae has been widely studied because it is the dominant yeast in most fermentations and it has been widely used as a model eukaryotic organism. Recently, other species of the Saccharomyces genus are gaining interest to solve the new challenges that the fermentation industry are facing. One of these species is S. kudriavzevii, which exhibits interesting physiological properties compared to S. cerevisiae, such as a better adaptation to grow at low temperatures, a higher glycerol synthesis and lower ethanol production. The aim of this study is to understand the molecular basis behind these phenotypic differences of biotechnological interest by using a species-based comparative genomics approach. In this work, we sequenced, assembled and annotated two new genomes of S. kudriavzevii. We used a combination of different statistical methods to identify functional divergence, signatures of positive selection and acceleration of substitution rates at specific amino acid sites of proteins in S. kudriavzevii when compared to S. cerevisiae, and vice versa. We provide a list of candidate genes in which positive selection could be acting during the evolution of both S. cerevisiae and S. kudriavzevii clades. Some of them could be related to certain important differences in metabolism previously reported by other authors such us DAL3 and ARO4, involved in nitrogen assimilation and amino acid biosynthesis. In addition, three of those genes (FBA1, ZIP1, and RQC2) showed accelerated evolutionary rates in Sk branch. Finally, genes of the riboflavin biosynthesis were also among those genes with a significant higher rate of nucleotide substitution and those proteins have amino acid positions contributing to functional divergence.

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.

Interplay of Chimeric Mating-Type Loci Impairs Fertility Rescue and Accounts for Intra-Strain Variability in Zygosaccharomyces rouxii Interspecies Hybrid ATCC42981

The pre-whole genome duplication (WGD) Zygosaccharomyces clade comprises several allodiploid strain/species with industrially interesting traits. The salt-tolerant yeast ATCC42981 is a sterile and allodiploid strain which contains two subgenomes, one of them resembling the haploid parental species Z. rouxii. Recently, different mating-type-like (MTL) loci repertoires were reported for ATCC42981 and the Japanese strain JCM22060, which are considered two stocks of the same strain. MTL reconstruction by direct sequencing approach is challenging due to gene redundancy, structure complexities, and allodiploid nature of ATCC42981. Here, DBG2OLC and MaSuRCA hybrid de novo assemblies of ONT and Illumina reads were combined with in vitro long PCR to definitively solve these incongruences. ATCC42981 exhibits several chimeric MTL loci resulting from reciprocal translocation between parental haplotypes and retains two MATa/MATα expression loci, in contrast to MATα in JCM22060. Consistently to these reconstructions, JCM22060, but not ATCC42981, undergoes mating and meiosis. To ascertain whether the damage of one allele at the MAT locus regains the complete sexual cycle in ATCC42981, we removed the MATα expressed locus by gene deletion. The resulting MATa/- hemizygous mutants did not show any evidence of sporulation, as well as of self- and out-crossing fertility, probably because incomplete silencing at the chimeric HMLα cassette masks the loss of heterozygosity at the MAT locus. We also found that MATα deletion switched off a2 transcription, an activator of a-specific genes in pre-WGD species. These findings suggest that regulatory scheme of cell identity needs to be further investigated in Z. rouxii protoploid yeast.

Linking gene expression and oenological traits: Comparison between Torulaspora delbrueckii and Saccharomyces cerevisiae strains

Wine fermentations typically involve the yeast Saccharomyces cerevisiae. However, many other yeast species participate to the fermentation process, some with interesting oenological traits. In this study the species Torulaspora delbrueckii, used occasionally in mixed or sequential fermentation with S. cerevisiae to improve wine sensory profile, was investigated to understand the physiological differences between the two. Next generation sequencing was used to characterize the transcriptome of T. delbrueckii and highlight the different genomic response of these yeasts during growth under wine-like conditions. Of particular interest were the basic differences in the glucose fermentation pathway and the formation of aromatic and flavour compounds such as glycerol, esters and acetic acid. Paralog genes were missing in glycolysis and glycerol biosynthesis in T. delbrueckii. Results indicate the tendency of T. delbrueckii to produce less acetic acid relied on a higher expression of alcoholic fermentation related genes, whereas acetate esters were influenced by the absence of esterases, ATF1-2. Additionally, in the Δbap2 S. cerevisiae strain, the final concentration of short branched chain ethyl esters (SBCEEs) was related to branched chain amino acid (BCAA) uptake. In conclusion, different adaption strategies are apparent for T. delbrueckii and S. cerevisiae yeasts, an understanding of which will allow winemakers to make better use of such microbial tools to achieve a desired wine sensory outcome.

Translatome and transcriptome analysis of TMA20 (MCT-1) and TMA64 (eIF2D) knockout yeast strains

TMA20 (MCT-1), TMA22 (DENR) and TMA64 (eIF2D) are eukaryotic translation factors involved in ribosome recycling and re-initiation. They operate with P-site bound tRNA in post-termination or (re-)initiation translation complexes, thus participating in the removal of 40S ribosomal subunit from mRNA stop codons after termination and controlling translation re-initiation on mRNAs with upstream open reading frames (uORFs), as well as de novo initiation on some specific mRNAs. Here we report ribosomal profiling data of S.cerevisiae strains with individual deletions of TMA20, TMA64 or both TMA20 and TMA64 genes. We provide RNA-Seq and Ribo-Seq data from yeast strains grown in the rich YPD or minimal SD medium. We illustrate our data by plotting differential distribution of ribosomal-bound mRNA fragments throughout uORFs in 5'-untranslated region (5' UTR) of GCN4 mRNA and on mRNA transcripts encoded in MAT locus in the mutant and wild-type strains, thus providing a basis for investigation of the role of these factors in the stress response, mating and sporulation. We also document a shift of transcription start site of the APC4 gene which occurs when the neighboring TMA64 gene is replaced by the standard G418-resistance cassette used for the creation of the Yeast Deletion Library. This shift results in dramatic deregulation of the APC4 gene expression, as revealed by our Ribo-Seq data, which can be probably used to explain strong genetic interactions of TMA64 with genes involved in the cell cycle and mitotic checkpoints. Raw RNA-Seq and Ribo-Seq data as well as all gene counts are available in NCBI Gene Expression Omnibus (GEO) repository under GEO accession GSE122039 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE122039).

Intrinsic cooperativity potentiates parallel cis-regulatory evolution

Convergent evolutionary events in independent lineages provide an opportunity to understand why evolution favors certain outcomes over others. We studied such a case where a large set of genes-those coding for the ribosomal proteins-gained cis-regulatory sequences for a particular transcription regulator (Mcm1) in independent fungal lineages. We present evidence that these gains occurred because Mcm1 shares a mechanism of transcriptional activation with an ancestral regulator of the ribosomal protein genes, Rap1. Specifically, we show that Mcm1 and Rap1 have the inherent ability to cooperatively activate transcription through contacts with the general transcription factor TFIID. Because the two regulatory proteins share a common interaction partner, the presence of one ancestral cis-regulatory sequence can 'channel' random mutations into functional sites for the second regulator. At a genomic scale, this type of intrinsic cooperativity can account for a pattern of parallel evolution involving the fixation of hundreds of substitutions.

Comparative whole-genome analysis reveals artificial selection effects on Ustilago esculenta genome

Ustilago esculenta, infects Zizania latifolia, and induced host stem swollen to be a popular vegetable called Jiaobai in China. It is the long-standing artificial selection that maximizes the occurrence of favourable Jiaobai, and thus maintaining the plant-fungi interaction and modulating the fungus evolving from plant pathogen to entophyte. In this study, whole genome of U. esculenta was sequenced and transcriptomes of the fungi and its host were analysed. The 20.2 Mb U. esculenta draft genome of 6,654 predicted genes including mating, primary metabolism, secreted proteins, shared a high similarity to related Smut fungi. But U. esculenta prefers RNA silencing not repeat-induced point in defence and has more introns per gene, indicating relatively slow evolution rate. The fungus also lacks some genes in amino acid biosynthesis pathway which were filled by up-regulated host genes and developed distinct amino acid response mechanism to balance the infection-resistance interaction. Besides, U. esculenta lost some surface sensors, important virulence factors and host range-related effectors to maintain the economic endophytic life. The elucidation of the U. esculenta genomic information as well as expression profiles can not only contribute to more comprehensive insights into the molecular mechanism underlying artificial selection but also into smut fungi-host interactions.

Genome Mining of Non-Conventional Yeasts: Search and Analysis of MAL Clusters and Proteins

Genomic clustering of functionally related genes is rare in yeasts and other eukaryotes with only few examples available. Here, we summarize our data on a nontelomeric MAL cluster of a non-conventional methylotrophic yeast Ogataea (Hansenula) polymorpha containing genes for α-glucosidase MAL1, α-glucoside permease MAL2 and two hypothetical transcriptional activators. Using genome mining, we detected MAL clusters of varied number, position and composition in many other maltose-assimilating non-conventional yeasts from different phylogenetic groups. The highest number of MAL clusters was detected in Lipomyces starkeyi while no MAL clusters were found in Schizosaccharomyces pombe and Blastobotrys adeninivorans. Phylograms of α-glucosidases and α-glucoside transporters of yeasts agreed with phylogenesis of the respective yeast species. Substrate specificity of unstudied α-glucosidases was predicted from protein sequence analysis. Specific activities of Scheffersomycesstipitis α-glucosidases MAL7, MAL8, and MAL9 heterologously expressed in Escherichia coli confirmed the correctness of the prediction-these proteins were verified promiscuous maltase-isomaltases. α-Glucosidases of earlier diverged yeasts L. starkeyi, B. adeninivorans and S. pombe showed sequence relatedness with α-glucosidases of filamentous fungi and bacilli.

Genome Sequence of Torulaspora microellipsoides CLIB 830T

We report here the genome sequence of the ascomycetous yeast Torulaspora microellipsoides CLIB 830^T. 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.

We report here the genome sequence of the ascomycetous yeast Torulaspora microellipsoides CLIB 830^T. 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.

The complete mitochondrial DNA sequence from Kazachstania sinensis reveals a general +1C frameshift mechanism in CTGY codons

The complete mitochondrial DNA (mtDNA) sequence from Kazachstania sinensis was analysed and compared to mtDNA from related yeasts. It contained the same set of genes; however, it only contained 23 tRNAs, as the trnR2 gene was absent. Most of the 12 introns within cox1, cob and rnl genes were inserted in the same sites as in other yeasts; however, two introns in rnl were in unusual positions. Traits such as gene order and GC cluster number were more related to Saccharomyces than to the other Kazachstania or linked clades. The most exceptional feature was the +1 frameshift in cox3, atp6 and cob open reading frames that was also found in other Kazachstania, Nakaseomyces delphensis and Candida glabrata. Comparison of DNA and protein sequences revealed the universal sites of +1C frameshifts were either CTGT or CTGC sequences. Moreover, an A→G substitution was found at position 37 in the anticodon stem loop tRNA gene for cysteine in all species with frameshifts but not in other sibling yeasts. This substitution allowed strong Watson-Crick base-pairing between an unmodified G (ACG) and the skipped C in the CTGY, leading to this quadruplet being read as cysteine.

Genome Sequence of Australian Indigenous Wine Yeast Torulaspora delbrueckii COFT1 Using Nanopore Sequencing

Here, we report the first sequenced genome of an indigenous Australian wine isolate of Torulaspora delbrueckii using the Oxford Nanopore MinION and Illumina HiSeq sequencing platforms. The genome size is 9.4 Mb and contains 4,831 genes.

Misidentification of genome assemblies in public databases: The case of Naumovozyma dairenensis and proposal of a protocol to correct misidentifications

Online sequence databases such as NCBI GenBank serve as a tremendously useful platform for researchers to share and reuse published data. However, submission systems lack control for errors such as organism misidentification, which once entered in the database can be propagated and mislead downstream analyses. Here we present an illustrating case of misidentification of Candida albicans from a clinical sample as Naumovozyma dairenensis based on whole-genome shotgun data. Analyses of phylogenetic markers, read mapping and single nucleotide polymorphisms served to correct the identification. We propose that the routine use of such analyses could help to detect misidentifications arising from unsupervised analyses and correct them before they enter the databases. Finally, we discuss broader implications of such misidentifications and the difficulty of correcting them once they are in the records.

High-oleate yeast oil without polyunsaturated fatty acids

BACKGROUND

Oleate-enriched triacylglycerides are well-suited for lubricant applications that require high oxidative stability. Fatty acid carbon chain length and degree of desaturation are key determinants of triacylglyceride properties and the ability to manipulate fatty acid composition in living organisms is critical to developing a source of bio-based oil tailored to meet specific application requirements.

RESULTS

We sought to engineer the oleaginous yeast Yarrowia lipolytica for production of high-oleate triacylglyceride oil. We studied the effect of deletions and overexpressions in the fatty acid and triacylglyceride synthesis pathways to identify modifications that increase oleate levels. Oleic acid accumulation in triacylglycerides was promoted by exchanging the native ∆9 fatty acid desaturase and glycerol-3-phosphate acyltransferase with heterologous enzymes, as well as deletion of the Δ12 fatty acid desaturase and expression of a fatty acid elongase. By combining these engineering steps, we eliminated polyunsaturated fatty acids and created a Y. lipolytica strain that accumulates triglycerides with > 90% oleate content.

CONCLUSIONS

High-oleate content and lack of polyunsaturates distinguish this triacylglyceride oil from plant and algal derived oils. Its composition renders the oil suitable for applications that require high oxidative stability and further demonstrates the potential of Y. lipolytica as a producer of tailored lipid profiles.

Reconstructing Yeasts Phylogenies and Ancestors from Whole Genome Data

Phylogenetic studies aim to discover evolutionary relationships and histories. These studies are based on similarities of morphological characters and molecular sequences. Currently, widely accepted phylogenetic approaches are based on multiple sequence alignments, which analyze shared gene datasets and concatenate/coalesce these results to a final phylogeny with maximum support. However, these approaches still have limitations, and often have conflicting results with each other. Reconstructing ancestral genomes helps us understand mechanisms and corresponding consequences of evolution. Most existing genome level phylogeny and ancestor reconstruction methods can only process simplified real genome datasets or simulated datasets with identical genome content, unique genome markers, and limited types of evolutionary events. Here, we provide an alternative way to resolve phylogenetic problems based on analyses of real genome data. We use phylogenetic signals from all types of genome level evolutionary events, and overcome the conflicting issues existing in traditional phylogenetic approaches. Further, we build an automated computational pipeline to reconstruct phylogenies and ancestral genomes for two high-resolution real yeast genome datasets. Comparison results with recent studies and publications show that we reconstruct very accurate and robust phylogenies and ancestors. Finally, we identify and analyze the conserved syntenic blocks among reconstructed ancestral genomes and present yeast species.

Mechanism for Restoration of Fertility in Hybrid Zygosaccharomyces rouxii Generated by Interspecies Hybridization

The mechanism of whole-genome duplication (WGD) in yeast has been intensively studied because it has a large impact on yeast evolution. WGD has shaped the genomic architecture of modern Saccharomyces cerevisiae; however, the mechanism for restoring fertility after interspecies hybridization, which would be involved in the process of WGD, has not been thoroughly elucidated. In this study, we obtained a draft genome sequence of the salt-tolerant yeast Zygosaccharomyces rouxii NBRC110957 and revealed that it is a hybrid lineage of Z. rouxii (allodiploid) with two subgenomes equivalent to NBRC1876. Because this allodiploid yeast can mate with other allodiploid strains and form spores, it can be a good model of restoring fertility after interspecies hybridization. We observed that NBRC110957 and NBRC1876 contain six mating-type-like (MTL) loci. There are no large deletions or deleterious mutations in MTL loci, except for several-base-pair deletions in the X region in certain MTL loci. We also assigned only one mating-type (MAT) locus that exclusively determines mating types from six MTL loci. These results suggest that it is possible to recover mating competence regardless of whether cells lose one MAT locus through random gene loss by mitotically dividing after interspecies hybridization. Moreover, we propose that perturbation of gene expression and substantial breakdown of MAT heterozygosity caused by chromosomal rearrangement at MTL loci play roles in restoring the mating competence of allodiploids. This scenario can provide a mechanism for restoring fertility after interspecies hybridization that is compatible with random gene loss models and suggests genomic plasticity during WGD in yeast.IMPORTANCE A whole-genome duplication occurred in an ancestor of the baker's yeast Saccharomyces cerevisiae The origins of this complex and multifaceted process, which requires intra- or interspecies hybridization followed by dysfunction of one mating-type (MAT) locus to regain mating competence, has not been thoroughly elucidated. In this study, we provide a mechanism for regaining fertility in an interspecies hybrid, Zygosaccharomyces rouxii The draft genome sequence analysis and mating test showed that the Z. rouxii strain used in this study is an intact interspecies hybrid, suggesting that it is possible to recover fertility regardless of whether cells lose one MAT locus.

Self-fertility in Chromocrea spinulosa is a consequence of direct repeat-mediated loss of MAT1-2, subsequent imbalance of nuclei differing in mating type, and recognition between unlike nuclei in a common cytoplasm

The filamentous fungus Chromocrea spinulosa (Trichoderma spinulosum) exhibits both self-fertile (homothallic) and self-sterile (heterothallic) sexual reproductive behavior. Self-fertile strains produce progeny cohorts that are 50% homothallic, 50% heterothallic. Heterothallic progeny can mate only with homothallic strains, and progeny also segregate 50% homothallic, 50% heterothallic. Sequencing of the mating type (MAT) region of homothallic and heterothallic strains revealed that both carry an intact MAT1-1 locus with three MAT1-1 genes (MAT1-1-1, MAT1-1-2, MAT1-1-3), as previously described for the Sordariomycete group of filamentous fungi. Homothallic strains, however, have a second version of MAT with the MAT1-2 locus genetically linked to MAT1-1. In this version, the MAT1-1-1 open reading frame is split into a large and small fragment and the truncated ends are bordered by 115bp direct repeats (DR). The MAT1-2-1 gene and additional sequences are inserted between the repeats. To understand the mechanism whereby C. spinulosa can exhibit both homothallic and heterothallic behavior, we utilized molecular manipulation to delete one of the DRs from a homothallic strain and insert MAT1-2 into a heterothallic strain. Mating assays indicated that: i) the DRs are key to homothallic behavior, ii) looping out of MAT1-2-1 via intra-molecular homologous recombination between the DRs in self-fertile strains results in two nuclear types in an individual (one carrying both MAT1-1 and MAT1-2 and one carrying MAT1-1 only), iii) self-fertility is achieved by inter-nuclear recognition between these two nuclear types before meiosis, iv) the two types of nuclei are in unequal proportion, v) having both an intact MAT1-1-1 and MAT1-2-1 gene in a single nucleus is not sufficient for self-fertility, and vi) the large truncated MAT1-1-1 fragment is expressed. Comparisons with MAT regions of Trichoderma reesei and Trichoderma virens suggest that several crossovers between misaligned parental MAT chromosomes may have led to the MAT architecture of homothallic C. spinulosa.

Fungi employ one of two mating tactics for sexual reproduction: self-sterile/heterothallic species can mate only with a genetically distinct partner while self-fertile/homothallic species do not require a partner. In ascomycetes, sexual reproduction is controlled by master regulators encoded by the mating-type (MAT) locus. The architecture of MAT differs in heterothallic versus homothallic species; heterothallics carry one of two forms (MAT1-1 or MAT1-2) per nucleus, whereas most homothallics carry both MAT forms in a single nucleus. There are intriguing exceptions. For example, the yeast models, Saccharomyces cerevisiae, and Schizosaccharomyces pombe undergo reversible MAT switching, not demonstrated in filamentous fungi. Here, we describe the mating mechanism in Chromocrea spinulosa (Trichoderma spinulosum), a filamentous ascomycete that exhibits both homothallic and heterothallic behavior. Self-fertile strains produce progeny cohorts that are 50% homothallic, 50% heterothallic. Self-sterile strains can mate only with homothallic strains, and when this occurs, homothallic and heterothallic progeny are also produced in a 1:1 ratio. By MAT sequencing and manipulation, we discovered unique MAT architecture and determined that self-fertility is achieved by deletion of MAT1-2 from most homothallic nuclei and subsequent inter-nuclear recognition between the resulting two, unevenly present, nuclear types in a common cytoplasm.