A computational pipeline to discover highly phylogenetically informative genes in sequenced genomes: application to Saccharomyces cerevisiae natural strains

Integrated genomics and phenotype microarray analysis of Saccharomyces cerevisiae industrial strains for rice wine fermentation and recombinant protein production

The industrial potential of Saccharomyces cerevisiae has extended beyond its traditional use in fermentation to various applications, including recombinant protein production. Herein, comparative genomics was performed with three industrial S. cerevisiae strains and revealed a heterozygous diploid genome for the 98-5 and KSD-YC strains (exploited for rice wine fermentation) and a haploid genome for strain Y2805 (used for recombinant protein production). Phylogenomic analysis indicated that Y2805 was closely associated with the reference strain S288C, whereas KSD-YC and 98-5 were grouped with Asian and European wine strains, respectively. Particularly, a single nucleotide polymorphism (SNP) in FDC1, involved in the biosynthesis of 4-vinylguaiacol (4-VG, a phenolic compound with a clove-like aroma), was found in KSD-YC, consistent with its lack of 4-VG production. Phenotype microarray (PM) analysis showed that KSD-YC and 98-5 displayed broader substrate utilization than S288C and Y2805. The SNPs detected by genome comparison were mapped to the genes responsible for the observed phenotypic differences. In addition, detailed information on the structural organization of Y2805 selection markers was validated by Sanger sequencing. Integrated genomics and PM analysis elucidated the evolutionary history and genetic diversity of industrial S. cerevisiae strains, providing a platform to improve fermentation processes and genetic manipulation.

MinION Sequencing of Yeast Mock Communities To Assess the Effect of Databases and ITS-LSU Markers on the Reliability of Metabarcoding Analysis

Microbial communities play key roles both for humans and the environment. They are involved in ecosystem functions, maintaining their stability, and provide important services, such as carbon cycle and nitrogen cycle. Acting both as symbionts and as pathogens, description of the structure and composition of these communities is important. Metabarcoding uses ribosomal DNA (rDNA) (eukaryotic) or rRNA gene (prokaryotic) sequences for identification of species present in a site and measuring their abundance. This procedure requires several technical steps that could be source of bias producing a distorted view of the real community composition. In this work, we took advantage of an innovative "long-read" next-generation sequencing (NGS) technology (MinION) amplifying the DNA spanning from the internal transcribed spacer (ITS) to large subunit (LSU) that can be read simultaneously in this platform, providing more information than "short-read" systems. The experimental system consisted of six fungal mock communities composed of species present at various relative amounts to mimic natural situations characterized by predominant and low-frequency species. The influence of the sequencing platform (MinION and Illumina MiSeq) and the effect of different reference databases and marker sequences on metagenomic identification of species were evaluated. The results showed that the ITS-based database provided more accurate species identification than LSU. Furthermore, a procedure based on a preliminary identification with standard reference databases followed by the production of custom databases, including only the best outputs of the first step, is proposed. This additional step improved the estimate of species proportion of the mock communities and reduced the number of ghost species not really present in the simulated communities. IMPORTANCE Metagenomic analyses are fundamental in many research areas; therefore, improvement of methods and protocols for the description of microbial communities becomes more and more necessary. Long-read sequencing could be used for reducing biases due to the multicopy nature of rDNA sequences and short-read limitations. However, these novel technologies need to be assessed and standardized with controlled experiments, such as mock communities. The interest behind this work was to evaluate how long reads performed identification and quantification of species mixed in precise proportions and how the choice of database affects such analyses. Development of a pipeline that mitigates the effect of the barcoding sequences and the impact of the reference database on metagenomic analyses can help microbiome studies go one step further.

Next-Generation Sequencing Applications for the Study of Fungal Pathogens

Next-generation sequencing (NGS) has become a widely used technology in biological research. NGS applications for clinical pathogen detection have become vital technologies. It is increasingly common to perform fast, accurate, and specific detection of clinical specimens using NGS. Pathogenic fungi with high virulence and drug resistance cause life-threatening clinical infections. NGS has had a significant biotechnological impact on detecting bacteria and viruses but is not equally applicable to fungi. There is a particularly urgent clinical need to use NGS to help identify fungi causing infections and prevent negative impacts. This review summarizes current research on NGS applications for fungi and offers a visual method of fungal detection. With the development of NGS and solutions for overcoming sequencing limitations, we suggest clinicians test specimens as soon as possible when encountering infections of unknown cause, suspected infections in vital organs, or rapidly progressive disease.

Comparison of Phylogenetic Tree Topologies for Nitrogen Associated Genes Partially Reconstruct the Evolutionary History of Saccharomyces cerevisiae

Massive sequencing projects executed in Saccharomyces cerevisiae have revealed in detail its population structure. The recent "1002 yeast genomes project" has become the most complete catalogue of yeast genetic diversity and a powerful resource to analyse the evolutionary history of genes affecting specific phenotypes. In this work, we selected 22 nitrogen associated genes and analysed the sequence information from the 1011 strains of the "1002 yeast genomes project". We constructed a total evidence (TE) phylogenetic tree using concatenated information, which showed a 27% topology similarity with the reference (REF) tree of the "1002 yeast genomes project". We also generated individual phylogenetic trees for each gene and compared their topologies, identifying genes with similar topologies (suggesting a shared evolutionary history). Furthermore, we pruned the constructed phylogenetic trees to compare the REF tree topology versus the TE tree and the individual genes trees, considering each phylogenetic cluster/subcluster within the population, observing genes with cluster/subcluster topologies of high similarity to the REF tree. Finally, we used the pruned versions of the phylogenetic trees to compare four strains considered as representatives of S. cerevisiae clean lineages, observing for 15 genes that its cluster topologies match 100% the REF tree, supporting that these strains represent main lineages of yeast population. Altogether, our results showed the potential of tree topologies comparison for exploring the evolutionary history of a specific group of genes.

Population genomics reveals evolution and variation of Saccharomyces cerevisiae in the human and insects gut

The quest to discover the variety of ecological niches inhabited by Saccharomyces cerevisiae has led to research in areas as diverse as wineries, oak trees and insect guts. The discovery of fungal communities in the human gastrointestinal tract suggested the host's gut as a potential reservoir for yeast adaptation. Here, we report the existence of yeast populations associated with the human gut (HG) that differ from those isolated from other human body sites. Phylogenetic analysis on 12 microsatellite loci and 1715 combined CDSs from whole-genome sequencing revealed three subclusters of HG strains with further evidence of clonal colonization within the host's gut. The presence of such subclusters was supported by other genomic features, such as copy number variation, absence/introgressions of CDSs and relative polymorphism frequency. Functional analysis of CDSs specific of the different subclusters suggested possible alterations in cell wall composition and sporulation features. The phenotypic analysis combined with immunological profiling of these strains further showed that sporulation was related with strain-specific genomic characteristics in the immune recognition pattern. We conclude that both genetic and environmental factors involved in cell wall remodelling and sporulation are the main drivers of adaptation in S. cerevisiae populations in the human gut.

Social wasp intestines host the local phenotypic variability of Saccharomyces cerevisiae strains

Nowadays, the presence of Saccharomyces cerevisiae has been assessed in both wild and human-related environments. Social wasps have been shown to maintain and vector S. cerevisiae among different environments. The availability of strains isolated from wasp intestines represents a striking opportunity to assess whether the strains found in wasp intestines are characterized by peculiar traits. We analysed strains isolated from the intestines of social wasps and compared them with strains isolated from other sources, all collected in a restricted geographic area. We evaluated the production of volatile metabolites during grape must fermentation, the resistance to different stresses and the ability to exploit various carbon sources. Wasp strains, in addition to representing a wide range of S. cerevisiae genotypes, also represent large part of the phenotypes characterizing the sympatric set of yeast strains; their higher production of acetic acid and ethyl acetate could reflect improved ability to attract insects. Our findings suggest that the relationship between yeasts and wasps should be preserved, to safeguard not only the natural variance of this microorganism but also the interests of wine-makers, who could take advantage from the exploitation of their phenotypic variability. Copyright © 2016 John Wiley & Sons, Ltd.

Multi-gene phylogenetic analysis reveals that shochu-fermenting Saccharomyces cerevisiae strains form a distinct sub-clade of the Japanese sake cluster

Shochu is a traditional Japanese distilled spirit. The formation of the distinguishing flavour of shochu produced in individual distilleries is attributed to putative indigenous yeast strains. In this study, we performed the first (to our knowledge) phylogenetic classification of shochu strains based on nucleotide gene sequences. We performed phylogenetic classification of 21 putative indigenous shochu yeast strains isolated from 11 distilleries. All of these strains were shown or confirmed to be Saccharomyces cerevisiae, sharing species identification with 34 known S. cerevisiae strains (including commonly used shochu, sake, ale, whisky, bakery, bioethanol and laboratory yeast strains and clinical isolate) that were tested in parallel. Our analysis used five genes that reflect genome-level phylogeny for the strain-level classification. In a first step, we demonstrated that partial regions of the ZAP1, THI7, PXL1, YRR1 and GLG1 genes were sufficient to reproduce previous sub-species classifications. In a second step, these five analysed regions from each of 25 strains (four commonly used shochu strains and the 21 putative indigenous shochu strains) were concatenated and used to generate a phylogenetic tree. Further analysis revealed that the putative indigenous shochu yeast strains form a monophyletic group that includes both the shochu yeasts and a subset of the sake group strains; this cluster is a sister group to other sake yeast strains, together comprising a sake-shochu group. Differences among shochu strains were small, suggesting that it may be possible to correlate subtle phenotypic differences among shochu flavours with specific differences in genome sequences. Copyright © 2017 John Wiley & Sons, Ltd.

The genetic architecture of low-temperature adaptation in the wine yeast Saccharomyces cerevisiae

Background

Low-temperature growth and fermentation of wine yeast can enhance wine aroma and make them highly desirable traits for the industry. Elucidating response to cold in Saccharomyces cerevisiae is, therefore, of paramount importance to select or genetically improve new wine strains. As most enological traits of industrial importance in yeasts, adaptation to low temperature is a polygenic trait regulated by many interacting loci.

Results

In order to unravel the genetic determinants of low-temperature fermentation, we mapped quantitative trait loci (QTLs) by bulk segregant analyses in the F13 offspring of two Saccharomyces cerevisiae industrial strains with divergent performance at low temperature. We detected four genomic regions involved in the adaptation at low temperature, three of them located in the subtelomeric regions (chromosomes XIII, XV and XVI) and one in the chromosome XIV. The QTL analysis revealed that subtelomeric regions play a key role in defining individual variation, which emphasizes the importance of these regions’ adaptive nature.

Conclusions

The reciprocal hemizygosity analysis (RHA), run to validate the genes involved in low-temperature fermentation, showed that genetic variation in mitochondrial proteins, maintenance of correct asymmetry and distribution of phospholipid in the plasma membrane are key determinants of low-temperature adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3572-2) contains supplementary material, which is available to authorized users.

A New Orthology Assessment Method for Phylogenomic Data: Unrooted Phylogenetic Orthology

Current sequencing technologies are making available unprecedented amounts of genetic data for a large variety of species including nonmodel organisms. Although many phylogenomic surveys spend considerable time finding orthologs from the wealth of sequence data, these results do not transcend the original study and after being processed for specific phylogenetic purposes these orthologs do not become stable orthology hypotheses. We describe a procedure to detect and document the phylogenetic distribution of orthologs allowing researchers to use this information to guide selection of loci best suited to test specific evolutionary questions. At the core of this pipeline is a new phylogenetic orthology method that is neither affected by the position of the root nor requires explicit assignment of outgroups. We discuss the properties of this new orthology assessment method and exemplify its utility for phylogenomics using a small insects dataset. In addition, we exemplify the pipeline to identify and document stable orthologs for the group of orb-weaving spiders (Araneoidea) using RNAseq data. The scripts used in this study, along with sample files and additional documentation, are available at https://github.com/ballesterus/UPhO.

Social wasps are a Saccharomyces mating nest

The reproductive ecology of Saccharomyces cerevisiae is still largely unknown. Recent evidence of interspecific hybridization, high levels of strain heterozygosity, and prion transmission suggest that outbreeding occurs frequently in yeasts. Nevertheless, the place where yeasts mate and recombine in the wild has not been identified. We found that the intestine of social wasps hosts highly outbred S. cerevisiae strains as well as a rare S. cerevisiae×S. paradoxus hybrid. We show that the intestine of Polistes dominula social wasps favors the mating of S. cerevisiae strains among themselves and with S. paradoxus cells by providing a succession of environmental conditions prompting cell sporulation and spores germination. In addition, we prove that heterospecific mating is the only option for European S. paradoxus strains to survive in the gut. Taken together, these findings unveil the best hidden secret of yeast ecology, introducing the insect gut as an environmental alcove in which crosses occur, maintaining and generating the diversity of the ascomycetes.

Fast method for identifying inter- and intra-species Saccharomyces hybrids in extensive genetic improvement programs based on yeast breeding

AIMS

The present work proposes a two-step molecular strategy to select inter- and intra-species Saccharomyces hybrids obtained by spore-to-spore mating, one of the most used methods for generating improved hybrids from homothallic wine yeasts.

METHODS AND RESULTS

As low spore viability and haplo-selfing are the main causes of failed mating, at first, we used colony screening PCR (csPCR) of discriminative gene markers to select hybrids directly on dissection plate and discard homozygous diploid colonies arisen from one auto-diploidized progenitor. Then, pre-selected candidates were submitted to recursive streaking and conventional PCR in order to discriminate between the hybrids with stable genomic background and the false-positive admixtures of progenitor cells both undergone haplo-selfing. csPCRs of internal transcribed spacer (ITS) 1 or 2, and the subsequent digestion with diagnostic endonucleases HaeIII and RsaI, respectively, were efficient to select six new Saccharomyces cerevisiae × Saccharomyces uvarum hybrids from 64 crosses. Intragenic minisatellite regions in PIR3, HSP150, and DAN4 genes showed high inter-strain size variation detectable by cost-effective agarose gel electrophoresis and were successful to validate six new intra-species S. cerevisiae hybrids from 34 crosses.

CONCLUSIONS

Both protocols reduce significantly the number of massive DNA extractions, prevent misinterpretations caused by one or both progenitors undergone haplo-selfing, and can be easily implemented in yeast labs without any specific instrumentation.

SIGNIFICANCE AND IMPACT OF THE STUDY

The study provides a method for the marker-assisted selection of several inter- and intra-species yeast hybrids in a cost-effective, rapid and reproducible manner.

Genomic evolution of Saccharomyces cerevisiae under Chinese rice wine fermentation

Rice wine fermentation represents a unique environment for the evolution of the budding yeast, Saccharomyces cerevisiae. To understand how the selection pressure shaped the yeast genome and gene regulation, we determined the genome sequence and transcriptome of a S. cerevisiae strain YHJ7 isolated from Chinese rice wine (Huangjiu), a popular traditional alcoholic beverage in China. By comparing the genome of YHJ7 to the lab strain S288c, a Japanese sake strain K7, and a Chinese industrial bioethanol strain YJSH1, we identified many genomic sequence and structural variations in YHJ7, which are mainly located in subtelomeric regions, suggesting that these regions play an important role in genomic evolution between strains. In addition, our comparative transcriptome analysis between YHJ7 and S288c revealed a set of differentially expressed genes, including those involved in glucose transport (e.g., HXT2, HXT7) and oxidoredutase activity (e.g., AAD10, ADH7). Interestingly, many of these genomic and transcriptional variations are directly or indirectly associated with the adaptation of YHJ7 strain to its specific niches. Our molecular evolution analysis suggested that Japanese sake strains (K7/UC5) were derived from Chinese rice wine strains (YHJ7) at least approximately 2,300 years ago, providing the first molecular evidence elucidating the origin of Japanese sake strains. Our results depict interesting insights regarding the evolution of yeast during rice wine fermentation, and provided a valuable resource for genetic engineering to improve industrial wine-making strains.

Ribosomal DNA sequence heterogeneity reflects intraspecies phylogenies and predicts genome structure in two contrasting yeast species

The ribosomal RNA encapsulates a wealth of evolutionary information, including genetic variation that can be used to discriminate between organisms at a wide range of taxonomic levels. For example, the prokaryotic 16S rDNA sequence is very widely used both in phylogenetic studies and as a marker in metagenomic surveys and the internal transcribed spacer region, frequently used in plant phylogenetics, is now recognized as a fungal DNA barcode. However, this widespread use does not escape criticism, principally due to issues such as difficulties in classification of paralogous versus orthologous rDNA units and intragenomic variation, both of which may be significant barriers to accurate phylogenetic inference. We recently analyzed data sets from the Saccharomyces Genome Resequencing Project, characterizing rDNA sequence variation within multiple strains of the baker's yeast Saccharomyces cerevisiae and its nearest wild relative Saccharomyces paradoxus in unprecedented detail. Notably, both species possess single locus rDNA systems. Here, we use these new variation datasets to assess whether a more detailed characterization of the rDNA locus can alleviate the second of these phylogenetic issues, sequence heterogeneity, while controlling for the first. We demonstrate that a strong phylogenetic signal exists within both datasets and illustrate how they can be used, with existing methodology, to estimate intraspecies phylogenies of yeast strains consistent with those derived from whole-genome approaches. We also describe the use of partial Single Nucleotide Polymorphisms, a type of sequence variation found only in repetitive genomic regions, in identifying key evolutionary features such as genome hybridization events and show their consistency with whole-genome Structure analyses. We conclude that our approach can transform rDNA sequence heterogeneity from a problem to a useful source of evolutionary information, enabling the estimation of highly accurate phylogenies of closely related organisms, and discuss how it could be extended to future studies of multilocus rDNA systems. [concerted evolution; genome hydridisation; phylogenetic analysis; ribosomal DNA; whole genome sequencing; yeast]

Genomic sequence diversity and population structure of Saccharomyces cerevisiae assessed by RAD-seq

The budding yeast Saccharomyces cerevisiae is important for human food production and as a model organism for biological research. The genetic diversity contained in the global population of yeast strains represents a valuable resource for a number of fields, including genetics, bioengineering, and studies of evolution and population structure. Here, we apply a multiplexed, reduced genome sequencing strategy (restriction site-associated sequencing or RAD-seq) to genotype a large collection of S. cerevisiae strains isolated from a wide range of geographical locations and environmental niches. The method permits the sequencing of the same 1% of all genomes, producing a multiple sequence alignment of 116,880 bases across 262 strains. We find diversity among these strains is principally organized by geography, with European, North American, Asian, and African/S. E. Asian populations defining the major axes of genetic variation. At a finer scale, small groups of strains from cacao, olives, and sake are defined by unique variants not present in other strains. One population, containing strains from a variety of fermentations, exhibits high levels of heterozygosity and a mixture of alleles from European and Asian populations, indicating an admixed origin for this group. We propose a model of geographic differentiation followed by human-associated admixture, primarily between European and Asian populations and more recently between European and North American populations. The large collection of genotyped yeast strains characterized here will provide a useful resource for the broad community of yeast researchers.

Role of social wasps in Saccharomyces cerevisiae ecology and evolution

Saccharomyces cerevisiae is one of the most important model organisms and has been a valuable asset to human civilization. However, despite its extensive use in the last 9,000 y, the existence of a seasonal cycle outside human-made environments has not yet been described. We demonstrate the role of social wasps as vector and natural reservoir of S. cerevisiae during all seasons. We provide experimental evidence that queens of social wasps overwintering as adults (Vespa crabro and Polistes spp.) can harbor yeast cells from autumn to spring and transmit them to their progeny. This result is mirrored by field surveys of the genetic variability of natural strains of yeast. Microsatellites and sequences of a selected set of loci able to recapitulate the yeast strain's evolutionary history were used to compare 17 environmental wasp isolates with a collection of strains from grapes from the same region and more than 230 strains representing worldwide yeast variation. The wasp isolates fall into subclusters representing the overall ecological and industrial yeast diversity of their geographic origin. Our findings indicate that wasps are a key environmental niche for the evolution of natural S. cerevisiae populations, the dispersion of yeast cells in the environment, and the maintenance of their diversity. The close relatedness of several wasp isolates with grape and wine isolates reflects the crucial role of human activities on yeast population structure, through clonal expansion and selection of specific strains during the biotransformation of fermented foods, followed by dispersal mediated by insects and other animals.