Genome-wide association analysis of clinical vs. nonclinical origin provides insights into Saccharomyces cerevisiae pathogenesis

From beer to breadboards: yeast as a force for biological innovation

The history of yeast Saccharomyces cerevisiae, aka brewer's or baker's yeast, is intertwined with our own. Initially domesticated 8,000 years ago to provide sustenance to our ancestors, for the past 150 years, yeast has served as a model research subject and a platform for technology. In this review, we highlight many ways in which yeast has served to catalyze the fields of functional genomics, genome editing, gene-environment interaction investigation, proteomics, and bioinformatics-emphasizing how yeast has served as a catalyst for innovation. Several possible futures for this model organism in synthetic biology, drug personalization, and multi-omics research are also presented.

Lipid saturation controls nuclear envelope function

The nuclear envelope (NE) is a spherical double membrane with elastic properties. How NE shape and elasticity are regulated by lipid chemistry is unknown. Here we discover lipid acyl chain unsaturation as essential for NE and nuclear pore complex (NPC) architecture and function. Increased lipid saturation rigidifies the NE and the endoplasmic reticulum into planar, polygonal membranes, which are fracture prone. These membranes exhibit a micron-scale segregation of lipids into ordered and disordered phases, excluding NPCs from the ordered phase. Balanced lipid saturation is required for NPC integrity, pore membrane curvature and nucleocytoplasmic transport. Oxygen deprivation amplifies the impact of saturated lipids, causing NE rigidification and rupture. Conversely, lipid droplets buffer saturated lipids to preserve NE architecture. Our study uncovers a fundamental link between lipid acyl chain structure and the integrity of the cell nucleus with implications for nuclear membrane malfunction in ischaemic tissues.

Genomic and Molecular Identification of Genes Contributing to the Caspofungin Paradoxical Effect in Aspergillus fumigatus

Aspergillus fumigatus is a deadly opportunistic fungal pathogen responsible for ~100,000 annual deaths. Azoles are the first line antifungal agent used against A. fumigatus, but azole resistance has rapidly evolved making treatment challenging. Caspofungin is an important second-line therapy against invasive pulmonary aspergillosis, a severe A. fumigatus infection. Caspofungin functions by inhibiting β-1,3-glucan synthesis, a primary and essential component of the fungal cell wall. A phenomenon termed the caspofungin paradoxical effect (CPE) has been observed in several fungal species where at higher concentrations of caspofungin, chitin replaces β-1,3-glucan, morphology returns to normal, and growth rate increases. CPE appears to occur in vivo, and it is therefore clinically important to better understand the genetic contributors to CPE. We applied genomewide association (GWA) analysis and molecular genetics to identify and validate candidate genes involved in CPE. We quantified CPE across 67 clinical isolates and conducted three independent GWA analyses to identify genetic variants associated with CPE. We identified 48 single nucleotide polymorphisms (SNPs) associated with CPE. We used a CRISPR/Cas9 approach to generate gene deletion mutants for seven genes harboring candidate SNPs. Two null mutants, ΔAfu3g13230 and ΔAfu4g07080 (dscP), resulted in reduced basal growth rate and a loss of CPE. We further characterized the dscP phosphatase-null mutant and observed a significant reduction in conidia production and extremely high sensitivity to caspofungin at both low and high concentrations. Collectively, our work reveals the contribution of Afu3g13230 and dscP in CPE and sheds new light on the complex genetic interactions governing this phenotype. IMPORTANCE This is one of the first studies to apply genomewide association (GWA) analysis to identify genes involved in an Aspergillus fumigatus phenotype. A. fumigatus is an opportunistic fungal pathogen that causes hundreds of thousands of infections and ~100,000 deaths each year, and antifungal resistance has rapidly evolved in this species. A phenomenon called the caspofungin paradoxical effect (CPE) occurs in some isolates, where high concentrations of the drug lead to increased growth rate. There is clinical relevance in understanding the genetic basis of this phenotype, since caspofungin concentrations could lead to unintended adverse clinical outcomes in certain cases. Using GWA analysis, we identified several interesting candidate polymorphisms and genes and then generated gene deletion mutants to determine whether these genes were important for CPE. Two of these mutant strains (ΔAfu3g13230 and ΔAfu4g07080/ΔdscP) displayed a loss of the CPE. This study sheds light on the genes involved in clinically important phenotype CPE.

Extensive simulations assess the performance of genome-wide association mapping in various Saccharomyces cerevisiae subpopulations

With the advent of high throughput sequencing technologies, genome-wide association studies (GWAS) have become a powerful paradigm for dissecting the genetic origins of the observed phenotypic variation. We recently completely sequenced the genome of 1011 Saccharomyces cerevisiae isolates, laying a strong foundation for GWAS. To assess the feasibility and the limits of this approach, we performed extensive simulations using five selected subpopulations as well as the total set of 1011 genomes. We measured the ability to detect the causal genetic variants involved in Mendelian and more complex traits using a linear mixed model approach. The results showed that population structure is well accounted for and is not the main problem when the sample size is high enough. While the genetic determinant of a Mendelian trait is easily mapped in all studied subpopulations, discrepancies are seen between datasets when performing GWAS on a complex trait in terms of detection, false positive and false negative rate. Finally, we performed GWAS on the different defined subpopulations using a real quantitative trait (resistance to copper sulfate) and showed the feasibility of this approach. The performance of each dataset depends simultaneously on several factors such as sample size, relatedness and population evolutionary history.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Population genomics of the pathogenic yeast Candida tropicalis identifies hybrid isolates in environmental samples

Candida tropicalis is a human pathogen that primarily infects the immunocompromised. Whereas the genome of one isolate, C. tropicalis MYA-3404, was originally sequenced in 2009, there have been no large-scale, multi-isolate studies of the genetic and phenotypic diversity of this species. Here, we used whole genome sequencing and phenotyping to characterize 77 isolates of C. tropicalis from clinical and environmental sources from a variety of locations. We show that most C. tropicalis isolates are diploids with approximately 2-6 heterozygous variants per kilobase. The genomes are relatively stable, with few aneuploidies. However, we identified one highly homozygous isolate and six isolates of C. tropicalis with much higher heterozygosity levels ranging from 36-49 heterozygous variants per kilobase. Our analyses show that the heterozygous isolates represent two different hybrid lineages, where the hybrids share one parent (A) with most other C. tropicalis isolates, but the second parent (B or C) differs by at least 4% at the genome level. Four of the sequenced isolates descend from an AB hybridization, and two from an AC hybridization. The hybrids are MTLa/α heterozygotes. Hybridization, or mating, between different parents is therefore common in the evolutionary history of C. tropicalis. The new hybrids were predominantly found in environmental niches, including from soil. Hybridization is therefore unlikely to be associated with virulence. In addition, we used genotype-phenotype correlation and CRISPR-Cas9 editing to identify a genome variant that results in the inability of one isolate to utilize certain branched-chain amino acids as a sole nitrogen source.

Genome-Wide Association for Itraconazole Sensitivity in Non-resistant Clinical Isolates of Aspergillus fumigatus

Aspergillus fumigatus is a potentially lethal opportunistic pathogen that infects over ~200,000 people and causes ~100,000 deaths per year globally. Treating A. fumigatus infections is particularly challenging because of the recent emergence of azole-resistance. The majority of studies focusing on the molecular mechanisms underlying azole resistance have examined azole-resistant isolates. However, isolates that are susceptible to azoles also display variation in their sensitivity, presenting a unique opportunity to identify genes contributing to azole sensitivity. Here, we used genome-wide association (GWA) analysis to identify loci involved in azole sensitivity by analyzing the association between 68,853 SNPs and itraconazole (ITCZ) minimum inhibitory concentration (MIC) in 76 clinical isolates of A. fumigatus from Japan. Population structure analysis suggests the presence of four distinct populations, with ITCZ MICs distributed relatively evenly across populations. We independently conducted GWA when treating ITCZ MIC as a quantitative trait and a binary trait, and identified two SNPs with strong associations in both analyses. These SNPs fell within the coding regions of Afu2g02220 and Afu2g02140. We functionally validated Afu2g02220 by knocking it out using a CRISPR/Cas9 approach, because orthologs of this gene are involved in sterol modification and ITCZ targets the ergosterol biosynthesis pathway. Knockout strains displayed no difference in growth compared to the parent strain in minimal media, yet a minor but consistent inhibition of growth in the presence of 0.15 μg/ml ITCZ. Our results suggest that GWA paired with efficient gene deletion is a powerful and unbiased strategy for identifying the genetic basis of complex traits in A. fumigatus.

Baker's Yeast Clinical Isolates Provide a Model for How Pathogenic Yeasts Adapt to Stress

Global outbreaks of drug-resistant fungi such as Candida auris are thought to be due at least in part to excessive use of antifungal drugs. Baker's yeast Saccharomyces cerevisiae has gained importance as an emerging opportunistic fungal pathogen that can cause infections in immunocompromised patients. Analyses of over 1000 S. cerevisiae isolates are providing rich resources to better understand how fungi can grow in human environments. A large percentage of clinical S. cerevisiae isolates are heterozygous across many nucleotide sites, and a significant proportion are of mixed ancestry and/or are aneuploid or polyploid. Such features potentially facilitate adaptation to new environments. These observations provide strong impetus for expanding genomic and molecular studies on clinical and wild isolates to understand the prevalence of genetic diversity and instability-generating mechanisms, and how they are selected for and maintained. Such work can also lead to the identification of new targets for antifungal drugs.

Gene-gene and gene-environment interactions in complex traits in yeast

One of the fundamental questions in biology is how the genotype regulates the phenotype. An increasing number of studies indicate that, in most cases, the effect of a genetic locus on the phenotype is context-dependent, i.e. it is influenced by the genetic background and the environment in which the phenotype is measured. Still, the majority of the studies, in both model organisms and humans, that map the genetic regulation of phenotypic variation in complex traits primarily identify additive loci with independent effects. This does not reflect an absence of the contribution of genetic interactions to phenotypic variation, but instead is a consequence of the technical limitations in mapping gene-gene interactions (GGI) and gene-environment interactions (GEI). Yeast, with its detailed molecular understanding, diverse population genomics and ease of genetic manipulation, is a unique and powerful resource to study the contributions of GGI and GEI in the regulation of phenotypic variation. Here we review studies in yeast that have identified GGI and GEI that regulate phenotypic variation, and discuss the contribution of these findings in explaining missing heritability of complex traits, and how observations from these GGI and GEI studies enhance our understanding of the mechanisms underlying genetic robustness and adaptability that shape the architecture of the genotype-phenotype map.

Use of Population Genetics to Assess the Ecology, Evolution, and Population Structure of Coccidioides

Although Coccidioides genotypes are highly genetically variable, they cluster into discrete populations, which has implications for human infections.

During the past 20 years, a general picture of the genetic diversity and population structure of Coccidioides, the causal agent of coccidioidomycosis (Valley fever), has emerged. The genus consists of 2 genetically diverse species, C. immitis and C. posadasii, each of which contains 1 or more distinct populations with limited gene flow. Genotypic data indicate that C. immitis is divided into 2 subpopulations (central and southern California populations) and C. posadasii is divided into 3 subpopulations (Arizona, Mexico, and Texas/South America populations). However, admixture within and among these populations and the current paucity of environmental isolates limit our understanding of the population genetics of Coccidioides. We assessed population structure of Coccidioides in Arizona by analyzing 495 clinical and environmental isolates. Our findings confirm the population structure as previously described and indicate a finer scale population structure in Arizona. Environmental isolates appear to have higher genetic diversity than isolates from human patients.

Genome-Wide Screen for Saccharomyces cerevisiae Genes Contributing to Opportunistic Pathogenicity in an Invertebrate Model Host

Environmental opportunistic pathogens can exploit vulnerable hosts through expression of traits selected for in their natural environments. Pathogenicity is itself a complicated trait underpinned by multiple complex traits, such as thermotolerance, morphology, and stress response. The baker's yeast, Saccharomyces cerevisiae, is a species with broad environmental tolerance that has been increasingly reported as an opportunistic pathogen of humans. Here we leveraged the genetic resources available in yeast and a model insect species, the greater waxmoth Galleria mellonella, to provide a genome-wide analysis of pathogenicity factors. Using serial passaging experiments of genetically marked wild-type strains, a hybrid strain was identified as the most fit genotype across all replicates. To dissect the genetic basis for pathogenicity in the hybrid isolate, bulk segregant analysis was performed which revealed eight quantitative trait loci significantly differing between the two bulks with alleles from both parents contributing to pathogenicity. A second passaging experiment with a library of deletion mutants for most yeast genes identified a large number of mutations whose relative fitness differed in vivovs.in vitro, including mutations in genes controlling cell wall integrity, mitochondrial function, and tyrosine metabolism. Yeast is presumably subjected to a massive assault by the innate insect immune system that leads to melanization of the host and to a large bottleneck in yeast population size. Our data support that resistance to the innate immune response of the insect is key to survival in the host and identifies shared genetic mechanisms between S. cerevisiae and other opportunistic fungal pathogens.

High-resolution mapping of cis-regulatory variation in budding yeast

Genetic variants affecting gene-expression levels are a major source of phenotypic variation. The approximate locations of these variants can be mapped as expression quantitative trait loci (eQTLs); however, a major limitation of eQTLs is their low resolution, which precludes investigation of the causal variants and their molecular mechanisms. Here we report RNA-seq and full genome sequences for 85 diverse isolates of the yeast Saccharomyces cerevisiae-including wild, domesticated, and human clinical strains-which allowed us to perform eQTL mapping with 50-fold higher resolution than previously possible. In addition to variants in promoters, we uncovered an important role for variants in 3'UTRs, especially those affecting binding of the PUF family of RNA-binding proteins. The eQTLs are predominantly under negative selection, particularly those affecting essential genes and conserved genes. However, applying the sign test for lineage-specific selection revealed the polygenic up-regulation of dozens of biofilm suppressor genes in strains isolated from human patients, consistent with the key role of biofilms in fungal pathogenicity. In addition, a single variant in the promoter of a biofilm suppressor, NIT3, showed the strongest genome-wide association with clinical origin. Altogether, our results demonstrate the power of high-resolution eQTL mapping in understanding the molecular mechanisms of regulatory variation, as well as the natural selection acting on this variation that drives adaptation to environments, ranging from laboratories to vineyards to the human body.

Bacterial Genetic Architecture of Ecological Interactions in Co-culture by GWAS-Taking Escherichia coli and Staphylococcus aureus as an Example

How a species responds to such a biotic environment in the community, ultimately leading to its evolution, has been a topic of intense interest to ecological evolutionary biologists. Until recently, limited knowledge was available regarding the genotypic changes that underlie phenotypic changes. Our study implemented GWAS (Genome-Wide Association Studies) to illustrate the genetic architecture of ecological interactions that take place in microbial populations. By choosing 45 such interspecific pairs of Escherichia coli and Staphylococcus aureus strains that were all genotyped throughout the entire genome, we employed Q-ROADTRIPS to analyze the association between single SNPs and microbial abundance measured at each time point for bacterial populations reared in monoculture and co-culture, respectively. We identified a large number of SNPs and indels across the genomes (35.69 G clean data of E. coli and 50.41 G of S. aureus). We reported 66 and 111 SNPs that were associated with interaction in E. coli and S. aureus, respectively. 23 out of 66 polymorphic changes resulted in amino acid alterations.12 significant genes, such as murE, treA, argS, and relA, which were also identified in previous evolutionary studies. In S. aureus, 111 SNPs detected in coding sequences could be divided into 35 non-synonymous and 76 synonymous SNPs. Our study illustrated the potential of genome-wide association methods for studying rapidly evolving traits in bacteria. Genetic association study methods will facilitate the identification of genetic elements likely to cause phenotypes of interest and provide targets for further laboratory investigation.

Budding Yeast Strains and Genotype-Phenotype Mapping

A small number of well-studied laboratory strains of Saccharomyces cerevisiae, mostly derived from S288C, are used in yeast research. Although powerful, studies for understanding S288C do not always capture the phenotypic essence or the genetic complexity of S. cerevisiae biology. This is particularly problematic for multilocus phenotypes identified in laboratory strains because these loci have never been jointly exposed to natural selection and the corresponding phenotypes do not represent optimization for any particular purpose or environment. Isolation and sequencing of new natural yeast strains also reveal that the total sequence diversity of the S. cerevisiae global population is poorly sampled in common laboratory strains. Here we discuss methodologies required for using the natural genetic variation in yeast to complete a genotype-phenotype map.

Genome Sequences of Cyberlindnera fabianii 65, Pichia kudriavzevii 129, and Saccharomyces cerevisiae 131 Isolated from Fermented Masau Fruits in Zimbabwe

Cyberlindnera fabianii 65, Pichia kudriavzevii 129, and Saccharomyces cerevisiae 131 have been isolated from the microbiota of fermented masau fruits. C. fabianii and P. kudriavzevii especially harbor promising features for biotechnology and food applications. Here, we present the draft annotated genome sequences of these isolates.

A small secreted protein in Zymoseptoria tritici is responsible for avirulence on wheat cultivars carrying the Stb6 resistance gene

Zymoseptoria tritici is the causal agent of Septoria tritici blotch, a major pathogen of wheat globally and the most damaging pathogen of wheat in Europe. A gene-for-gene (GFG) interaction between Z. tritici and wheat cultivars carrying the Stb6 resistance gene has been postulated for many years, but the genes have not been identified. We identified AvrStb6 by combining quantitative trait locus mapping in a cross between two Swiss strains with a genome-wide association study using a natural population of c. 100 strains from France. We functionally validated AvrStb6 using ectopic transformations. AvrStb6 encodes a small, cysteine-rich, secreted protein that produces an avirulence phenotype on wheat cultivars carrying the Stb6 resistance gene. We found 16 nonsynonymous single nucleotide polymorphisms among the tested strains, indicating that AvrStb6 is evolving very rapidly. AvrStb6 is located in a highly polymorphic subtelomeric region and is surrounded by transposable elements, which may facilitate its rapid evolution to overcome Stb6 resistance. AvrStb6 is the first avirulence gene to be functionally validated in Z. tritici, contributing to our understanding of avirulence in apoplastic pathogens and the mechanisms underlying GFG interactions between Z. tritici and wheat.

Using Population and Comparative Genomics to Understand the Genetic Basis of Effector-Driven Fungal Pathogen Evolution

Epidemics caused by fungal plant pathogens pose a major threat to agro-ecosystems and impact global food security. High-throughput sequencing enabled major advances in understanding how pathogens cause disease on crops. Hundreds of fungal genomes are now available and analyzing these genomes highlighted the key role of effector genes in disease. Effectors are small secreted proteins that enhance infection by manipulating host metabolism. Fungal genomes carry 100s of putative effector genes, but the lack of homology among effector genes, even for closely related species, challenges evolutionary and functional analyses. Furthermore, effector genes are often found in rapidly evolving chromosome compartments which are difficult to assemble. We review how population and comparative genomics toolsets can be combined to address these challenges. We highlight studies that associated genome-scale polymorphisms with pathogen lifestyles and adaptation to different environments. We show how genome-wide association studies can be used to identify effectors and other pathogenicity-related genes underlying rapid adaptation. We also discuss how the compartmentalization of fungal genomes into core and accessory regions shapes the evolution of effector genes. We argue that an understanding of genome evolution provides important insight into the trajectory of host-pathogen co-evolution.

Validation of Genome-Wide Association Studies as a Tool to Identify Virulence Factors in Parastagonospora nodorum

Parastagonospora nodorum is a necrotrophic fungal pathogen causing Septoria nodorum blotch on wheat. We have identified nine necrotrophic effector-host dominant sensitivity gene interactions, and we have cloned three of the necrotrophic effector genes, including SnToxA, SnTox1, and SnTox3. Because sexual populations of P. nodorum are difficult to develop under lab conditions, genome-wide association study (GWAS) is the best population genomic approach to identify genomic regions associated with traits using natural populations. In this article, we used a global collection of 191 P. nodorum isolates from which we identified 2,983 single-nucleotide polymorphism (SNP) markers and gene markers for SnToxA and SnTox3 to evaluate the power of GWAS on two popular wheat breeding lines that were sensitive to SnToxA and SnTox3. Strong marker trait associations (MTA) with P. nodorum virulence that mapped to SnTox3 and SnToxA were first identified using the marker set described above. A novel locus in the P. nodorum genome associated with virulence was also identified as a result of this analysis. To evaluate whether a sufficient level of marker saturation was available, we designed a set of primers every 1 kb in the genomic regions containing SnToxA and SnTox3. Polymerase chain reaction amplification was performed across the 191 isolates and the presence/absence polymorphism was scored and used as the genotype. The marker proximity necessary to identify MTA flanking SnToxA and SnTox3 ranged from 4 to 5 and 1 to 7 kb, respectively. Similar analysis was performed on the novel locus. Using a 45% missing data threshold, two more SNP were identified spanning a 4.6-kb genomic region at the novel locus. These results showed that the rate of linkage disequilibrium (LD) decay in P. nodorum and, likely, other fungi is high compared with plants and animals. The fast LD decay in P. nodorum is an advantage only if sufficient marker density is attained. Based on our results with the SnToxA and SnTox3 regions, markers are needed every 9 or 8 kb, respectively, or in every gene, to guarantee that genes associated with a quantitative trait such as virulence are not missed.

Whole Genome Analysis of 132 Clinical Saccharomyces cerevisiae Strains Reveals Extensive Ploidy Variation

Budding yeast has undergone several independent transitions from commercial to clinical lifestyles. The frequency of such transitions suggests that clinical yeast strains are derived from environmentally available yeast populations, including commercial sources. However, despite their important role in adaptive evolution, the prevalence of polyploidy and aneuploidy has not been extensively analyzed in clinical strains. In this study, we have looked for patterns governing the transition to clinical invasion in the largest screen of clinical yeast isolates to date. In particular, we have focused on the hypothesis that ploidy changes have influenced adaptive processes. We sequenced 144 yeast strains, 132 of which are clinical isolates. We found pervasive large-scale genomic variation in both overall ploidy (34% of strains identified as 3n/4n) and individual chromosomal copy numbers (36% of strains identified as aneuploid). We also found evidence for the highly dynamic nature of yeast genomes, with 35 strains showing partial chromosomal copy number changes and eight strains showing multiple independent chromosomal events. Intriguingly, a lineage identified to be baker's/commercial derived with a unique damaging mutation in NDC80 was particularly prone to polyploidy, with 83% of its members being triploid or tetraploid. Polyploidy was in turn associated with a >2× increase in aneuploidy rates as compared to other lineages. This dataset provides a rich source of information on the genomics of clinical yeast strains and highlights the potential importance of large-scale genomic copy variation in yeast adaptation.

Identification of Diverse Mycoviruses through Metatranscriptomics Characterization of the Viromes of Five Major Fungal Plant Pathogens

Mycoviruses can have a marked effect on natural fungal communities and influence plant health and productivity. However, a comprehensive picture of mycoviral diversity is still lacking. To characterize the viromes of five widely dispersed plant-pathogenic fungi, Colletotrichum truncatum, Macrophomina phaseolina, Diaporthe longicolla, Rhizoctonia solani, and Sclerotinia sclerotiorum, a high-throughput sequencing-based metatranscriptomic approach was used to detect viral sequences. Total RNA and double-stranded RNA (dsRNA) from mycelia and RNA from samples enriched for virus particles were sequenced. Sequence data were assembled de novo, and contigs with predicted amino acid sequence similarities to viruses in the nonredundant protein database were selected. The analysis identified 72 partial or complete genome segments representing 66 previously undescribed mycoviruses. Using primers specific for each viral contig, at least one fungal isolate was identified that contained each virus. The novel mycoviruses showed affinity with 15 distinct lineages: Barnaviridae, Benyviridae, Chrysoviridae, Endornaviridae, Fusariviridae, Hypoviridae, Mononegavirales, Narnaviridae, Ophioviridae, Ourmiavirus, Partitiviridae, Tombusviridae, Totiviridae, Tymoviridae, and Virgaviridae More than half of the viral sequences were predicted to be members of the Mitovirus genus in the family Narnaviridae, which replicate within mitochondria. Five viral sequences showed strong affinity with three families (Benyviridae, Ophioviridae, and Virgaviridae) that previously contained no mycovirus species. The genomic information provides insight into the diversity and taxonomy of mycoviruses and coevolution of mycoviruses and their fungal hosts.

IMPORTANCE

Plant-pathogenic fungi reduce crop yields, which affects food security worldwide. Plant host resistance is considered a sustainable disease management option but may often be incomplete or lacking for some crops to certain fungal pathogens or strains. In addition, the rising issues of fungicide resistance demand alternative strategies to reduce the negative impacts of fungal pathogens. Those fungus-infecting viruses (mycoviruses) that attenuate fungal virulence may be welcome additions for mitigation of plant diseases. By high-throughput sequencing of the RNAs from 275 isolates of five fungal plant pathogens, 66 previously undescribed mycoviruses were identified. In addition to identifying new potential biological control agents, these results expand the grand view of the diversity of mycoviruses and provide possible insights into the importance of intracellular and extracellular transmission in fungus-virus coevolution.

Genome-Wide Association Study Identifies Novel Candidate Genes for Aggressiveness, Deoxynivalenol Production, and Azole Sensitivity in Natural Field Populations of Fusarium graminearum

Genome-wide association studies can identify novel genomic regions and genes that affect quantitative traits. Fusarium head blight is a destructive disease caused by Fusarium graminearum that exhibits several quantitative traits, including aggressiveness, mycotoxin production, and fungicide resistance. Restriction site-associated DNA sequencing was performed for 220 isolates of F. graminearum. A total of 119 isolates were phenotyped for aggressiveness and deoxynivalenol (DON) production under natural field conditions across four environments. The effective concentration of propiconazole that inhibits isolate growth in vitro by 50% was calculated for 220 strains. Approximately 29,000 single nucleotide polymorphism markers were associated to each trait, resulting in 50, 29, and 74 quantitative trait nucleotides (QTNs) that were significantly associated to aggressiveness, DON production, and propiconazole sensitivity, respectively. Approximately 41% of these QTNs caused nonsynonymous substitutions in predicted exons, while the remainder were synonymous substitutions or located in intergenic regions. Three QTNs associated with propiconazole sensitivity were significant after Bonferroni correction. These QTNs were located in genes not previously associated with azole sensitivity. The majority of the detected QTNs were located in genes with predicted regulatory functions, suggesting that nucleotide variation in regulatory genes plays a major role in the corresponding quantitative trait variation.

Population genomics of yeasts: towards a comprehensive view across a broad evolutionary scale

With the advent of high-throughput technologies for sequencing, the complete description of the genetic variation that occurs in populations, also known as population genomics, is foreseeable but far from being reached. Explaining the forces that govern patterns of genetic variation is essential to elucidate the evolutionary history of species. Genetic variation results from a wide assortment of evolutionary forces, among which mutation, selection, recombination and drift play major roles in shaping genomes. In addition, exploring the genetic variation within a population also corresponds to the first step towards dissecting the genotype-phenotype relationship. In this context, yeast species are of particular interest because they represent a unique resource for studying the evolution of intraspecific genetic diversity in a phylum spanning a broad evolutionary scale. Here, we briefly review recent progress in yeast population genomics and provide some perspective on this rapidly evolving field. In fact, we truly believe that it is of interest to supplement comparative and early population genomic studies with the deep sequencing of more extensive sets of individuals from the same species. In parallel, it would be more than valuable to uncover the intraspecific variation of a large number of unexplored species, including those that are closely and more distantly related. Altogether, these data would enable substantially more powerful genomic scans for functional dissection.

Sucrose and Saccharomyces cerevisiae: a relationship most sweet

Sucrose is an abundant, readily available and inexpensive substrate for industrial biotechnology processes and its use is demonstrated with much success in the production of fuel ethanol in Brazil. Saccharomyces cerevisiae, which naturally evolved to efficiently consume sugars such as sucrose, is one of the most important cell factories due to its robustness, stress tolerance, genetic accessibility, simple nutrient requirements and long history as an industrial workhorse. This minireview is focused on sucrose metabolism in S. cerevisiae, a rather unexplored subject in the scientific literature. An analysis of sucrose availability in nature and yeast sugar metabolism was performed, in order to understand the molecular background that makes S. cerevisiae consume this sugar efficiently. A historical overview on the use of sucrose and S. cerevisiae by humans is also presented considering sugarcane and sugarbeet as the main sources of this carbohydrate. Physiological aspects of sucrose consumption are compared with those concerning other economically relevant sugars. Also, metabolic engineering efforts to alter sucrose catabolism are presented in a chronological manner. In spite of its extensive use in yeast-based industries, a lot of basic and applied research on sucrose metabolism is imperative, mainly in fields such as genetics, physiology and metabolic engineering.

On the Mapping of Epistatic Genetic Interactions in Natural Isolates: Combining Classical Genetics and Genomics

Genetic variation within species is the substrate of evolution. Epistasis, which designates the non-additive interaction between loci affecting a specific phenotype, could be one of the possible outcomes of genetic diversity. Dissecting the basis of such interactions is of current interest in different fields of biology, from exploring the gene regulatory network, to complex disease genetics, to the onset of reproductive isolation and speciation. We present here a general workflow to identify epistatic interactions between independently evolving loci in natural populations of the yeast Saccharomyces cerevisiae. The idea is to exploit the genetic diversity present in the species by evaluating a large number of crosses and analyzing the phenotypic distribution in the offspring. For a cross of interest, both parental strains would have a similar phenotypic value, whereas the resulting offspring would have a bimodal distribution of the phenotype, possibly indicating the presence of epistasis. Classical segregation analysis of the tetrads uncovers the penetrance and complexity of the interaction. In addition, this segregation could serve as the guidelines for choosing appropriate mapping strategies to narrow down the genomic regions involved. Depending on the segregation patterns observed, we propose different mapping strategies based on bulk segregant analysis or consecutive backcrosses followed by high-throughput genome sequencing. Our method is generally applicable to all systems with a haplodiplobiontic life cycle and allows high resolution mapping of interacting loci that govern various DNA polymorphisms from single nucleotide mutations to large-scale structural variations.

Characterization of Cell Wall Proteins in Saccharomyces cerevisiae Clinical Isolates Elucidates Hsp150p in Virulence

The budding yeast Saccharomyces cerevisiae has recently been described as an emerging opportunistic fungal pathogen. Fungal cell wall mannoproteins have been demonstrated to be involved in adhesion to inert surfaces and might be engaged in virulence. In this study, we observed four clinical isolates of S. cerevisiae with relatively hydrophobic cell surfaces. Yeast cell wall subproteome was evaluated quantitatively by liquid chromatography/tandem mass spectrometry. We identified totally 25 cell wall proteins (CWPs) from log-phase cells, within which 15 CWPs were quantified. The abundance of Scw10p, Pst1p, and Hsp150p/Pir2p were at least 2 folds higher in the clinical isolates than in S288c lab strain. Hsp150p is one of the members in Pir family conserved in pathogenic fungi Candida glabrata and Candida albicans. Overexpression of Hsp150p in lab strain increased cell wall integrity and potentially enhanced the virulence of yeast. Altogether, these results demonstrated that quantitative cell wall subproteome was analyzed in clinical isolates of S. cerevisiae, and several CWPs, especially Hsp150p, were found to be expressed at higher levels which presumably contribute to strain virulence and fungal pathogenicity.

Review of current methods for characterizing virulence and pathogenicity potential of industrial Saccharomyces cerevisiae strains towards humans

Most industrial Saccharomyces cerevisiae strains used in food or biotechnology processes are benign. However, reports of S. cerevisiae infections have emerged and novel strains continue to be developed. In order to develop recommendations for the human health risk assessment of S. cerevisiae strains, we conducted a literature review of current methods used to characterize their pathogenic potential and evaluated their relevance towards risk assessment. These studies revealed that expression of virulence traits in S. cerevisiae is complex and depends on many factors. Given the opportunistic nature of this organism, an approach using multiple lines of evidence is likely necessary for the reasonable prediction of the pathogenic potential of a particular strain. Risk assessment of S. cerevisiae strains would benefit from more research towards the comparison of virulent and non-virulent strains in order to better understand those genotypic and phenotypic traits most likely to be associated with pathogenicity.

The 100-genomes strains, an S. cerevisiae resource that illuminates its natural phenotypic and genotypic variation and emergence as an opportunistic pathogen

Saccharomyces cerevisiae, a well-established model for species as diverse as humans and pathogenic fungi, is more recently a model for population and quantitative genetics. S. cerevisiae is found in multiple environments—one of which is the human body—as an opportunistic pathogen. To aid in the understanding of the S. cerevisiae population and quantitative genetics, as well as its emergence as an opportunistic pathogen, we sequenced, de novo assembled, and extensively manually edited and annotated the genomes of 93 S. cerevisiae strains from multiple geographic and environmental origins, including many clinical origin strains. These 93 S. cerevisiae strains, the genomes of which are near-reference quality, together with seven previously sequenced strains, constitute a novel genetic resource, the “100-genomes” strains. Our sequencing coverage, high-quality assemblies, and annotation provide unprecedented opportunities for detailed interrogation of complex genomic loci, examples of which we demonstrate. We found most phenotypic variation to be quantitative and identified population, genotype, and phenotype associations. Importantly, we identified clinical origin associations. For example, we found that an introgressed PDR5 was present exclusively in clinical origin mosaic group strains; that the mosaic group was significantly enriched for clinical origin strains; and that clinical origin strains were much more copper resistant, suggesting that copper resistance contributes to fitness in the human host. The 100-genomes strains are a novel, multipurpose resource to advance the study of S. cerevisiae population genetics, quantitative genetics, and the emergence of an opportunistic pathogen.

The fascinating and secret wild life of the budding yeast S. cerevisiae

The budding yeast Saccharomyces cerevisiae has been used in laboratory experiments for over a century and has been instrumental in understanding virtually every aspect of molecular biology and genetics. However, it wasn't until a decade ago that the scientific community started to realise how little was known about this yeast's ecology and natural history, and how this information was vitally important for interpreting its biology. Recent large-scale population genomics studies coupled with intensive field surveys have revealed a previously unappreciated wild lifestyle of S. cerevisiae outside the restrictions of human environments and laboratories. The recent discovery that Chinese isolates harbour almost twice as much genetic variation as isolates from the rest of the world combined suggests that Asia is the likely origin of the modern budding yeast.

DOI:http://dx.doi.org/10.7554/eLife.05835.001

Saccharomyces cerevisiae: a nomadic yeast with no niche?

Different species are usually thought to have specific adaptations, which allow them to occupy different ecological niches. But recent neutral ecology theory suggests that species diversity can simply be the result of random sampling, due to finite population sizes and limited dispersal. Neutral models predict that species are not necessarily adapted to specific niches, but are functionally equivalent across a range of habitats. Here, we evaluate the ecology of Saccharomyces cerevisiae, one of the most important microbial species in human history. The artificial collection, concentration and fermentation of large volumes of fruit for alcohol production produce an environment in which S. cerevisiae thrives, and therefore it is assumed that fruit is the ecological niche that S. cerevisiae inhabits and has adapted to. We find very little direct evidence that S. cerevisiae is adapted to fruit, or indeed to any other specific niche. We propose instead a neutral nomad model for S. cerevisiae, which we believe should be used as the starting hypothesis in attempting to unravel the ecology of this important microbe.

It is assumed that Saccharomyces cerevisiae is adapted to inhabit fruits; however, we find very little evidence for adaptation to any niche. Instead, we propose a neutral nomad model for S. cerevisiae.

Of parasites and men

The living world has evolved and is evolving through interspecific relationships between organisms. The diversity of these interactions is enormous going from mutualism to parasitism. Humans live with a multitude of microorganisms, essential for their biology. However, interactions are not always advantageous. Indeed, many organisms might become pathogens, such as the Plasmodium species, the causative agents of malaria. Like many other microorganisms, they are «Machiavellian» in their capacity to elaborate a range of reproduction strategies, giving them a huge advantage in terms of adaptation. Here, we discuss the role played by parasites in the ecology and evolution of living organisms and particularly of humans. In the study of infectious diseases, humans are legitimately the focal point, although they represent only one ecosystem among many others and not taking this into account certainly biases our global view of the system. Indeed, we know only a minimal fraction of the microorganisms we live with. However, parasites have shaped and are still shaping the human genome. Several genetic signatures are the proofs of the selection pressures by parasites that humankind has endured during its evolution. But, ultimately, what are the solutionsfor living with pathogens? Should we eradicate them or should we learn how to control and manage them?

A genome-wide association study identifies genomic regions for virulence in the non-model organism Heterobasidion annosum s.s

The dense single nucleotide polymorphisms (SNP) panels needed for genome wide association (GWA) studies have hitherto been expensive to establish and use on non-model organisms. To overcome this, we used a next generation sequencing approach to both establish SNPs and to determine genotypes. We conducted a GWA study on a fungal species, analysing the virulence of Heterobasidion annosum s.s., a necrotrophic pathogen, on its hosts Picea abies and Pinus sylvestris. From a set of 33,018 single nucleotide polymorphisms (SNP) in 23 haploid isolates, twelve SNP markers distributed on seven contigs were associated with virulence (P<0.0001). Four of the contigs harbour known virulence genes from other fungal pathogens and the remaining three harbour novel candidate genes. Two contigs link closely to virulence regions recognized previously by QTL mapping in the congeneric hybrid H. irregulare × H. occidentale. Our study demonstrates the efficiency of GWA studies for dissecting important complex traits of small populations of non-model haploid organisms with small genomes.

Tapping into yeast diversity

Domesticated organisms demonstrate our capacity to influence wild species but also provide us with the opportunity to understand rapid evolution in the context of substantially altered environments and novel selective pressures. Recent advances in genetics and genomics have brought unprecedented insights into the domestication of many organisms and have opened new avenues for further improvements to be made. Yet, our ability to engineer biological systems is not without limits; genetic manipulation is often quite difficult. The budding yeast, Saccharomyces cerevisiae, is not only one of the most powerful model organisms, but is also the premier producer of fermented foods and beverages around the globe. As a model system, it entertains a hefty workforce dedicated to deciphering its genome and the function it encodes at a rich mechanistic level. As a producer, it is used to make leavened bread, and dozens of different alcoholic beverages, such as beer and wine. Yet, applying the awesome power of yeast genetics to understanding its origins and evolution requires some knowledge of its wild ancestors and the environments from which they were derived. A number of surprisingly diverse lineages of S. cerevisiae from both primeval and secondary forests in China have been discovered by Wang and his colleagues. These lineages substantially expand our knowledge of wild yeast diversity and will be a boon to elucidating the ecology, evolution and domestication of this academic and industrial workhorse.

DNA mutations mediate microevolution between host-adapted forms of the pathogenic fungus Cryptococcus neoformans

The disease cryptococcosis, caused by the fungus Cryptococcus neoformans, is acquired directly from environmental exposure rather than transmitted person-to-person. One explanation for the pathogenicity of this species is that interactions with environmental predators select for virulence. However, co-incubation of C. neoformans with amoeba can cause a “switch” from the normal yeast morphology to a pseudohyphal form, enabling fungi to survive exposure to amoeba, yet conversely reducing virulence in mammalian models of cryptococcosis. Like other human pathogenic fungi, C. neoformans is capable of microevolutionary changes that influence the biology of the organism and outcome of the host-pathogen interaction. A yeast-pseudohyphal phenotypic switch also happens under in vitro conditions. Here, we demonstrate that this morphological switch, rather than being under epigenetic control, is controlled by DNA mutation since all pseudohyphal strains bear mutations within genes encoding components of the RAM pathway. High rates of isolation of pseudohyphal strains can be explained by the physical size of RAM pathway genes and a hypermutator phenotype of the strain used in phenotypic switching studies. Reversion to wild type yeast morphology in vitro or within a mammalian host can occur through different mechanisms, with one being counter-acting mutations. Infection of mice with RAM mutants reveals several outcomes: clearance of the infection, asymptomatic maintenance of the strains, or reversion to wild type forms and progression of disease. These findings demonstrate a key role of mutation events in microevolution to modulate the ability of a fungal pathogen to cause disease.

Many diseases are contracted from the environment, rather than from sick people. It is unclear why those species are able to cause disease, since the selective pressures in the environment are presumed to be very different from those found within the host. Cryptococcus neoformans is a fungus that causes life-threatening lung and central nervous system disease in approximately one million people each year. The fungus is inhaled from environmental sources. One hypothesis to account for C. neoformans virulence is that amoeba are predators for this fungus, and surviving strains are pre-selected to be virulent in the human host. On the other hand, experiments have found that amoeba eat C. neoformans. A pseudohyphal cell type can survive, and while protecting against amoeba these cells are unable to cause disease in mouse models. We predicted that the pseudohyphal morphology reflected a change in function of a pathway of genes, and found that all pseudohyphal isolates contain mutations within genes for this pathway. The pseudohyphal trait is unstable, with reversion to normal yeast growth by counter-acting mutations. These mutations can occur during the course of mammalian infection. Our results show that mutation events account for a microevolution system currently described as phenotypic switching, and that mutations, at least under experimental conditions, can regulate pathogen adaptation and influence its host range.

On the prospects of whole-genome association mapping in Saccharomyces cerevisiae

Advances in sequencing technology have enabled whole-genome sequences to be obtained from multiple individuals within species, particularly in model organisms with compact genomes. For example, 36 genome sequences of Saccharomyces cerevisiae are now publicly available, and SNP data are available for even larger collections of strains. One potential use of these resources is mapping the genetic basis of phenotypic variation through genome-wide association (GWA) studies, with the benefit that associated variants can be studied experimentally with greater ease than in outbred populations such as humans. Here, we evaluate the prospects of GWA studies in S. cerevisiae strains through extensive simulations and a GWA study of mitochondrial copy number. We demonstrate that the complex and heterogeneous patterns of population structure present in yeast populations can lead to a high type I error rate in GWA studies of quantitative traits, and that methods typically used to control for population stratification do not provide adequate control of the type I error rate. Moreover, we show that while GWA studies of quantitative traits in S. cerevisiae may be difficult depending on the particular set of strains studied, association studies to map cis-acting quantitative trait loci (QTL) and Mendelian phenotypes are more feasible. We also discuss sampling strategies that could enable GWA studies in yeast and illustrate the utility of this approach in Saccharomyces paradoxus. Thus, our results provide important practical insights into the design and interpretation of GWA studies in yeast, and other model organisms that possess complex patterns of population structure.

Polygenic cis-regulatory adaptation in the evolution of yeast pathogenicity

The acquisition of new genes, via horizontal transfer or gene duplication/diversification, has been the dominant mechanism thus far implicated in the evolution of microbial pathogenicity. In contrast, the role of many other modes of evolution—such as changes in gene expression regulation—remains unknown. A transition to a pathogenic lifestyle has recently taken place in some lineages of the budding yeast Saccharomyces cerevisiae. Here we identify a module of physically interacting proteins involved in endocytosis that has experienced selective sweeps for multiple cis-regulatory mutations that down-regulate gene expression levels in a pathogenic yeast. To test if these adaptations affect virulence, we created a panel of single-allele knockout strains whose hemizygous state mimics the genes' adaptive down-regulations, and measured their virulence in a mammalian host. Despite having no growth advantage in standard laboratory conditions, nearly all of the strains were more virulent than their wild-type progenitor, suggesting that these adaptations likely played a role in the evolution of pathogenicity. Furthermore, genetic variants at these loci were associated with clinical origin across 88 diverse yeast strains, suggesting the adaptations may have contributed to the virulence of a wide range of clinical isolates. We also detected pleiotropic effects of these adaptations on a wide range of morphological traits, which appear to have been mitigated by compensatory mutations at other loci. These results suggest that cis-regulatory adaptation can occur at the level of physically interacting modules and that one such polygenic adaptation led to increased virulence during the evolution of a pathogenic yeast.