Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Saccharomyces cerevisiae: Population Divergence and Resistance to Oxidative Stress in Clinical, Domesticated and Wild Isolates

  • Stephanie Diezmann ,

    sd21@duke.edu

    Affiliation Department of Molecular Genetics & Microbiology, Duke University Medical Center, Durham, North Carolina, United States of America

  • Fred S. Dietrich

    Affiliation Department of Molecular Genetics & Microbiology, Duke University Medical Center, Durham, North Carolina, United States of America

Abstract

Background

Saccharomyces cerevisiae has been associated with human life for millennia in the brewery and bakery. Recently it has been recognized as an emerging opportunistic pathogen. To study the evolutionary history of S. cerevisiae, the origin of clinical isolates and the importance of a virulence-associated trait, population genetics and phenotypic assays have been applied to an ecologically diverse set of 103 strains isolated from clinics, breweries, vineyards, fruits, soil, commercial supplements and insect guts.

Methodology/Principal Findings

DNA sequence data from five nuclear DNA loci were analyzed for population structure and haplotype distribution. Additionally, all strains were tested for survival of oxidative stress, a trait associated with microbial pathogenicity. DNA sequence analyses identified three genetic subgroups within the recombining S. cerevisiae strains that are associated with ecology, geography and virulence. Shared alleles suggest that the clinical isolates contain genetic contribution from the fruit isolates. Clinical and fruit isolates exhibit high levels of recombination, unlike the genetically homogenous soil isolates in which no recombination was detected. However, clinical and soil isolates were more resistant to oxidative stress than any other population, suggesting a correlation between survival in oxidative stress and yeast pathogenicity.

Conclusions/Significance

Population genetic analyses of S. cerevisiae delineated three distinct groups, comprising primarily the (i) human-associated brewery and vineyard strains, (ii) clinical and fruit isolates (iii) and wild soil isolates from eastern U.S. The interactions between S. cerevisiae and humans potentiate yeast evolution and the development of genetically, ecologically and geographically divergent groups.

Introduction

S. cerevisiae has long been associated with humans as the fermentative agent in the production of bread, beer, and wine. Archeological evidence for the production of fermented beverages in China dates to 7,000 BCE, and molecular evidence demonstrating S. cerevisiae was the fermentative agent has been found in wine jars from ancient Egypt dating to 3,150 BCE [1], [2]. The close relationship between humans and yeast is further reflected in molecular signatures recovered from African artifacts that contained palm wine and from European wine and beer vessels that can be traced to Mesopotamia [3], [4]. Due to its close association with humankind, it has been speculated that yeast may have been the first living being domesticated [5]. Yet, the ecology of S. cerevisiae embraces a wider range than domesticated strains found in the vineyard and the brewery. Wild strains have been isolated from mushroom fruiting bodies as well as oak tree-associated soils and fluxes [6][11]. Wild isolates of S. cerevisiae are furthermore a major cause of spoilage of mango fruit and peach puree, and it has recently been identified in surveys of the fungal diversity in beetle guts [12][14].

The breadth of S. cerevisiae ecological diversity, encompassing domesticated and wild isolates, has spurred interest in the life history and population genetics of this species. The population biology of an organism so tightly associated with humans is challenging to study due to sampling bias, limited sampling sizes, human influence, and non-random sampling [5], [8], [15]. Consequently, questions concerning population structure and genetic diversity in S. cerevisiae have mainly been addressed using strains from grape berries, vineyards and other industrial applications. Fay and Benavides analyzed approximately 7 kb of coding and non-coding DNA sequences in 81 strains from vineyards, fermentation of sake, palm wine, ragi and cider, fruit sources, including lychee, fig, and mushrooms, oak tree and surrounding soil from New Jersey, and patients in the U.S. and Europe. This extensive analysis resulted in the recognition of domesticated, i.e. human-associated, and wild populations in S. cerevisiae [3]. Aa et al. examined 6.6 kb of coding and non-coding sequences from 27 strains including soil and oak-associated isolates from Pennsylvania, vineyard strains from two locations in Italy and strains from rotten figs from California [16]. The authors found signatures of distinct population structure, moderate levels of recombination and demonstrated that the oak isolates form a monophyletic group. In summary, previous studies indicated that ecology rather than geography coincided with population structure and that clinical isolates are only distantly related to isolates used in fermentation.

In addition to investigations into the origin and consequences of domestication, S. cerevisiae's has been reported as an emerging opportunistic pathogen [17], [18]. Since the late 1950's, there have been increasing case reports of S. cerevisiae causing infections [19]. S. cerevisiae and its commercially available preparations known as S. boulardii, that are used to treat antibiotic-related diarrhea, have been shown to cause a wide variety of infections, ranging from cutaneous infections and vaginitis to systemic infections of the bloodstream and vital organs in immunocompromised and immunocompetent individuals [20][33]. These infections are similar to those by the related yeast Candida albicans, the most common human fungal pathogen [20], [34], [35]. Because of its evolutionary kinship with C. albicans and its status as genetics model system, S. cerevisiae, the benevolent baker's yeast, has acquired enhanced scientific value as a model pathogen in the study of virulence-related traits [36], [37].

Microbial virulence-related traits promote host invasion, colonization and virulence. One of these is the ability to survive oxidative stress, exerted by radical oxygen species (ROS), an integral component of mammalian host defenses, is associated with virulence in various bacteria [38][42]. Other virulence traits that have been studied in S. cerevisiae include growth at high temperature [43], the formation of multiple colony phenotypes [44], pseudohyphal growth [45], and loss of mitochondrial genome function [46]. These studies showed that clinical isolates differed phenotypically from laboratory and wine strains, but could not exclude the possibility that the observed association between clinical origin and a virulence trait is due to shared ancestry rather than host adaptation. If clinical isolates share a common ancestor, the evolution of virulence could be attributed to an isolated event that imparted selective advantage to one or more progenitor, pathogenic strains. However, if clinical isolates exhibit multiple evolutionary histories, pathogenicity would more likely reflect an adaptive advantage conferred by the acquisition of multiple virulence traits in different strains. This scenario is illustrated by the example of growth at 37°C, which is required for pathogenicity. It is teleologically reasonable to hypothesize that this trait is likely to be acquired by yeasts growing on fruit, which attain high temperatures during decomposition. This is supported by “heat death point” studies, which showed that S. cerevisiae isolates from spoiled peach pure and grape juice were heat resistant [13], [47]. Hence, clinical isolates could have acquired the ability to survive high temperatures via recombination with isolates from decomposing organic matter.

Despite progress in understanding the population structure of domesticated and wild S. cerevisiae, little is known about the emergence of this versatile yeast as a pathogen and the role of selection and/or adaptation in this evolutionary process. This study aims to (i) test for genetic differentiation between strains of a broad range of origins, (ii) investigate the evolutionary origin of clinical isolates, and (iii) identify an association between a virulence-related trait and pathogenicity. The targeted virulence trait is survival of oxidative stress. To investigate population structure and evolution of pathogenicity, an ecologically diverse sample of 103 S. cerevisiae isolates, comprising seven populations, was analyzed. Clinical isolates with confirmed virulence were obtained from the California Institute for Medical Research in Stanford and Duke University Medical Center in Durham, NC [43]. The domesticated, fermentation-associated, isolates consisted of brewery strains from Europe, and vineyard isolates collected from small commercial vineyards in from North Carolina and Australia (AU) [9], [48]. The wild strains consisted of fruit isolates, collected from tropical monocultures and fermenting fruits from different locales around the world, soil isolates from an arboretum in Pennsylvania (PA) and parks in North Carolina (NC) and insect guts from Louisiana. Commercial S. boulardii strains from France, Italy and Germany, complete the sampling [12], [33].

The results confirm previous findings and generated several novel conclusions. (i) At least three divergent groups representing different evolutionary trajectories were detected in S. cerevisiae, expanding the previous findings. (ii) Clinical isolates are genetically similar to the isolates from fruit, which supports the hypothesis that monoculture or fermenting fruits may serve as a natural reservoir for the evolution of clinical isolates. (iii) Strains from the clinic and Pennsylvania soil tolerate oxidative stress better than any other group. However, the clinical isolates are genetically diverse while the soil isolates are identical at all loci. This strongly suggests that resistance to ROS is an adaptive property of the clinical strains.

Results

All strains are diploid isolates of S. cerevisiae

All isolates in this study belong to the species S. cerevisiae, as confirmed by ITS sequencing and all strains were determined by FACS analysis to be diploid (data not shown). To ensure correct estimates of population genetic parameters haplotype phase within each locus was determined. For the majority of isolates the correct locus haplotype phase could be determined by visual analysis of sequence alignments. For 22 strains the correct phase of one ore more loci was identified using PHASE 2.1.1. [49], [50]. For ten of those haplotype phase assignments were verified by cloning and sequencing. Consequently, co-dominant marker assignments and two haplotypes per strain were used in the subsequent analyses. For calculations of population genetic parameters, strains sharing the same origin were grouped into populations. Seven populations were recognized including clinic, fruit, brewery, NC and AU vineyard, NC and PA soil.

Within S. cerevisiae three genetically divergent groups are recognized

Principal component analysis (PCA), a transformation method that reduces multidimensional data sets to lower dimensions, was employed to detect putative structure among 87 isolates for which complete sequence data for five loci could be obtained (Table 1). Confidence envelope calculations of PCA results, revealed three major groups of 19, 19, and 30 isolates, each were identified and designated A, B, and C (Figure 1). The clustering was supported by UPGMA analysis (Unweighted Pair Group Method Arithmetic) of a genetic distance matrix derived from strains included in PCA (Figure S1). The composition of groups A, B, and C differs. The majority of clinical isolates (60%), a third of the fruit and brewery (33%) isolates are clustered in group A. All soil isolates from Pennsylvania and North Carolina comprise group B. All of the isolates from the Australian vineyards, half of the those from the North Carolina vineyards (50%), almost half the brewery isolates (42%), one-fourth of the fruit isolates (24%), one clinical isolate, and S. boulardii form group C. Note that two fruit isolates in group C stem from grapes, two fruit isolates from Holly and a papaya. Interestingly, all group C vineyard isolates have been collected from Vitis vinifera in Australia or North Carolina. The two identical NC vineyard isolates that are outside the major clusters were isolated from V. rotundifolia grapes, the Muscadine grapes that are native to the southeastern United States. Seven of the nine isolates collected in Adelaide (AU) are identical. As measure of genetic differentiation pairwise Fst values were calculated for groups A, B, and C. Fst values showed significant genetic differentiation between the three groups and analysis of molecular variance (AMOVA) indicated that 25% of the observed variance occurs within these groups and 62% between them. The index of association (IA), measuring the extent of linkage equilibrium, was significantly different from zero for group A but not for B and C.

thumbnail
Figure 1. Three divergent groups in S. cerevisiae.

PCA recognized three major groups (A, B, and C) in a combined analysis of complete sequence data for strains (N = 87) from all ecological backgrounds. The size (N) and index of association (IA), which is a measure of linkage equilibrium, of each group are given next to the group name. The IA of each group has been calculated using sequence data for all five loci of each strain included in that group. Strain origins are coded with symbols and colors (see legend to the left). The arrows inside group C point to fruit isolates that have been isolated from grapes but are not part of the Australian or North Carolinian sampling. Strains of particular interest due to their history as lab strains (RM11, S288c) or importance as a model fungal pathogen (YJM145) have their names attached. Numbers inside symbols indicate how many strains share this genotype. Pairwise Fst values and significance levels for comparisons by origin (left) and PCA group (right) are indicated.

https://doi.org/10.1371/journal.pone.0005317.g001

thumbnail
Table 1. Origin, survival in oxidative stress and haplotype configurations of 103 strains included in this study.

https://doi.org/10.1371/journal.pone.0005317.t001

Levels of recombination, genetic diversity and linkage differ among S. cerevisiae populations

Seven populations encompassed a total of 82 strains for which complete sequence data were obtained. Five population genetics parameters concerning genetic differentiation (Fst), recombination, linkage equilibrium (IA), nucleotide diversity (π), and deviation from Hardy-Weinberg equilibrium (HWE) were calculated (Figure 1, Table 2).

thumbnail
Table 2. Nucleotide diversity (π), minimum number of recombination events, index of association (IA), and Hardy-Weinberg Equilibrium (HWE) for seven defined populations.

https://doi.org/10.1371/journal.pone.0005317.t002

Genetic differentiation varied between populations from different origins (Figure 1). Clinical and fruit isolates were not significantly different from each other, but the clinical isolates were distinct from all other isolates. Fruit isolates differed from all except NC vineyard isolates. There was a large degree of genetic differentiation between AU vineyard and NC and PA soils (0.822, 0.845), and soil samples from NC and PA exhibited extensive genetic differentiation from each other (0.831). NC and AU Vineyard populations differed less from each other (0.282). The brewery population was most different from the soil populations (0.445, 0.551) and less different from vineyard and clinical populations (0.239, 0.101, 0.124).

Linkage equilibrium and the minimum number of recombination events calculated for each population differ between populations (Table 2). The IA was significant for the NC vineyard, brewery and fruit populations. Each of the clinic, fruit and brewery populations could be explained by a minimum of five recombination events. The soil isolates showed no evidence of recombination and the vineyard isolates of only a single recombination event among these five loci. IA and recombination results were reflected in the differences in nucleotide diversity between populations. Fruit, brewery, clinical, and NC vineyard populations exhibit high relative nucleotide diversity. Low π values were calculated for isolates from NC and PA soils and the AU vineyard population. π values are to a large degree in concordance with observed heterozygosity (Table 1). Between 33% and 88% of fruit, brewery, clinical and NC vineyard strains are heterozygous at one ore more loci. Indeed, four brewery isolates are heterozygous at every locus tested. NC and PA soil isolates are completely homozygous and 28% of AU vineyard isolates were heterozygous for at least one locus.

All populations were tested for deviations from HWE at each locus. A population that deviates from HWE does so because of non-random mating, mutation, and selection and other factors that affect population structure. Only the fruit population differs significantly at all loci from HWE (Table 2). The clinical population differs at three loci, NC soil, NC and AU vineyard, and brewery populations deviate at least at one locus. HWE calculations could not be conducted on the PA soil population.

Clinical and domesticated isolates of S. cerevisiae exhibit haplotype diversity, soil isolates do not

Haplotype networks were inferred from the DNA sequence data of five nuclear coding loci, totaling 2,527 bp [51] (Figure 2, Table 1). Although haplotype diversity varied between loci, from 18 for PHD1 to 7 for ACT1 and ADP1, a common pattern was identified. Each network is characterized by the presence of one to three dominant haplotypes. Those accumulated mutations and gave rise to six to nineteen minor haplotypes in different networks, some of which have not been sampled in this study and are marked accordingly (Figure 2).

thumbnail
Figure 2. Haplotype diversity at five nuclear loci.

Two haplotypes per strain were analyzed for five nuclear coding loci. Lengths of analyzed sequence data and number of sampled haplotypes are given for each locus. The size of each pie represents the number of identical haplotypes and the proportions indicate how many of those share a particular origin. If N>20, the number of haplotypes is indicated in the pie. The length of the connecting lines translates into nucleotide substitutions distinguishing one haplotype from another. Haplotypes collected from the clinic, the brewery and fruit sources are randomly distributed in each network, whereas soil isolates share one or two haplotypes. No correlation between strain origin and haplotype could be detected and no haplotype that unifies all strains from one origin.

https://doi.org/10.1371/journal.pone.0005317.g002

The haplotype diversity between populations differs widely, and three groups with different levels of diversity can be discerned. The soil populations are the least diverse. With two exceptions, all strains from NC and PA soils are represented by one of the dominant haplotypes in each network. In ADP1, NC soil isolates are represented by two haplotypes, and in PHD1, NC and PA soil isolates have different haplotypes. Both vineyard populations exhibit an intermediate amount of haplotype diversity – between two and four types per network. The clinical, brewery and fruit populations demonstrate the greatest diversity, in concordance with the population genetic parameters described above. They are represented by three to eleven different types in each network. Haplotypes from isolates of the same population do not cluster; for example, at the MLS1 locus, clinical isolates originate from all three common haplotypes, not just one (Figure 2).

Clinical and soil PA isolates of S. cerevisiae are highly resistant to oxidative stress

All isolates were tested for survival of oxidative stress exerted by tert-Butyl hydroperoxide, a stable organic analog of H2O2. Clinical isolates, PA soil isolates, and strains from insect guts exhibited the highest mean survival rates (Figure 3A, Table 1). NC soil, both vineyard populations and the S. boulardii isolates show decreased survival and intermediate levels of variation. Brewery and fruit populations display survival rates ranging from very low to very high. Analysis of variance (ANOVA) confirmed significant differences in oxidative stress between groups of N≥10 (P<0.0001). Clinical isolates differ significantly (P<0.001) from all groups but PA soil (Figure 3B). PA soil isolates differ significantly (P<0.001) from the vineyard and the brewery isolates and only slightly (P<0.05) from the NC soil and fruit isolates. The remaining pair-wise comparisons were either not or only marginally significant (P<0.05).

thumbnail
Figure 3. Survival of 103 S. cerevisiae strains in 20 mM TBHP.

Box and whiskers plot for survival in 20 mM TBHP as tested in the CFU assay. Strains are clustered and color-coded by origin (Table 1). For each group of strains the colored box entails the size of the 25th and 75th percentile. The horizontal line dividing the box is the median (50th percentile) and the whiskers represent the most extreme outliers with highest or lowest survival. Experiments were carried out in triplicates and all three values of each strain were used in this plot (A). The multiple pair-wise comparisons of all groups of N≥10 indicate that clinical isolates and strains from the soil in PA differ significantly from strains from other demographic groups. The legend denotes significance values obtained in ANOVA of all three experimental replicates of each strain (B).

https://doi.org/10.1371/journal.pone.0005317.g003

Discussion

Although S. cerevisiae is a species of inter-fertile isolates [9], [52], [53], PCA delineated three genetically distinct groups with different strain compositions. Group A, comprises isolates from genetically diverse and recombining populations collected in the clinic, the brewery and fruits. Group B includes the genetically homogenous but distinct PA and NC soil populations. Group C unites the distinct yet homogenous NC and AU vineyard isolates with genetically diverse brewery strains. These results are similar to the findings from the Fay and Benavides and Aa et al. studies where they reported that human association has shaped S. cerevisiae life history and led to domesticated isolates distinct from wild soil isolates [3], [16]. This study identified a third group that is dominated by clinical isolates and genetically different from domesticated and wild lineages. These three life history trajectories of S. cerevisiae are characterized by different levels of recombination and genetic differentiation that could explain their origin and maintenance.

Low levels of recombination and genetic diversity among soil populations (group B) suggest a predominantly clonal life style that contributes to differentiation from clinical and domesticated groups and maintenance of a wild lineage that also exhibits geographic structure because all of its isolates were collected in the eastern United States. Interestingly, the soil isolates share four of five haplotypes with beetle gut isolates, suggesting that insects acquire S. cerevisiae while foraging on the ground. In concordance with the Fay and Benavides results, ecology, rather than geography, appears to dominate evolution and origin of domestication in the vineyard and the brewery (group C) [3]. Although yeast isolates from NC and AU vineyards, collected from wine grapes, are genetically differentiated, they were recognized as most similar to each other and different from the Muscadine grape isolates in PCA. The importance of ecology is further emphasized by the observed distribution of brewery strains in PCA groups A and C, which suggests that brewery strains in wine-growing regions originated from grapes, and in areas lacking viniculture, the domesticated brewery yeasts were derived from wild fruit isolates.

The clinical and PA soil populations exhibited the highest mean resistance to ROS. Therefore ROS resistance in the clinical isolates arose either independently or was contributed by the soil isolates. The genetic diversity among the clinical isolates suggests multiple origins of the genetic material required for ROS resistance in these strains. Regardless, fungal resistance to ROS offers protection from oxidative host defenses and is undoubtedly an advantageous pathobiological property. Consistent with this conclusion, the probiotic strains of S. boulardii were not as resistant to ROS as the clinical strains. Indeed, strains of S. boulardii were less virulent compared to a clinical (YJM145) and a laboratory (Y55) strain [33]. The S. boulardii isolates used here were cultured from commercial products sold in France, Italy and Germany. Their close association with vineyard strains suggests that S. boulardii may have originated from vineyard strains.

Most recently, a study conducted by Kvitek, Will and Gasch combined investigations into genetic diversity with stress response and identified a distinct sake fermentation clade but no discrete clades for clinical or oak-associated isolates [54]. The discrepancy between the results reported by Kvitek et al. and our study can be explained by differences in strain sampling and number as well as the kind of assays conducted. The overlap in strain sampling is minimal. One PA soil and five clinical isolates were included in both studies. Furthermore, experimental procedures varied. While Kvitek et al. employed a plate-based assay, cells here were exposed to peroxide under defined conditions in liquid. Interestingly, strains isolated from similar environments share comparable stress response expression profiling patterns and oak isolates appear to have been selected for growth in this particular niche [54].

This investigation presents an exciting and provocative demonstration of the complex life history of S. cerevisiae beyond its prosaic service as a scientific tool in the laboratory or an agent of fermentation. As a species, S. cerevisiae entered two new life history trajectories while continuing its life in the soil and on decomposing fruit. Some strains of S. cerevisiae became domesticated in fermentation and brewing, and while others became pathogenic. With the accumulation of ecological and population genetic data, this versatile microbe becomes an invaluable model for evolutionary biology and population genetics. The current report offers a new paradigm for studying pathogenesis by identifying correlation(s) among the virulence traits of isolates and the ecology of their ancestral reservoirs. These correlations will identify the genotypes associated with pathogenic strains or the potential for pathogenicity and elucidate the evolution of pathogenicity.

Materials and Methods

Yeast strains, culturing and DNA extraction

S. cerevisiae strains were isolated from S. cerevisiae infected patients (N = 15), soil in Pennsylvania (PA, N = 10), soil in North Carolina (NC, N = 14), vineyards in NC (N = 10) and Australia (AUS, N = 14), various fruits (N = 18), brewery (N = 16), commercial S. boulardii preparations (N = 4), and the insect gut (N = 2) (Table 1). All strains were permanently stored at −80°C in 15% glycerol or cultivated on rich media plates (YPD, 1% yeast extract, 2% peptone, 2% dextrose, 2% agar). DNA was extracted using the CTAB buffer method [55]. NC soil samples were collected as part of this study using a previously established enrichment protocol and ITS sequencing for identification [9], [56].

DNA sequencing and sequence analysis

Genomic DNAs were diluted 1∶100 and 2 µl added to a 25 µl reaction from TaKaRa Ex Taq kit with a final primer concentration of 0.6 µM (Table 3). The PCR regime consisted of 5 min initial denaturation at 95°C, followed by 35 cycles of 30 sec at 95°C, 30 sec at the appropriate annealing temperature (Table 3), and 45 sec at 72°C, concluding with 10 min extension at 72°C. PCR products were purified using the Montáge-SEQ96Cleanup Kit (Millipore) or the ExoSAP protocol and their concentration determined electrophoretically [57]. Between 500 and 1000 ng PCR product were sequenced with BigDye chemistry version 3 according to the instructions supplied by Applied Biosystems. Chromatograms were assembled and analyzed in Sequencher 4.8 (Gene Codes Corp.) and sequences edited in MacClade 4.06 [58]. ITS sequence data were blasted against GenBank using the nr database and the blastn algorithm [59].

thumbnail
Table 3. Primers used for amplification of nuclear loci and ITS.

https://doi.org/10.1371/journal.pone.0005317.t003

The haplotype phase of nuclear loci in strains with more than one heterozygous site per locus was determined using PHASE 2.1.1. [49], [50]. Loci that received low confidence probabilities were cloned using the pCR2.1 TOPO TA Cloning Kit (Invitrogen) and sequenced as described above. Two haplotypes per strain were further analyzed. PCA, HWE, AMOVA, Fst, π, IA were calculated in GenAlEx 6, DnaSP 4.50.3, and MultiLocus 1.3 [60][62]. 98.9% confidence envelopes were calculated for three PCA groups based on three times the standard deviation of the PCA scores of strains included in each cluster [63]. Strains with missing data points were excluded from these analyses and statistical significance values were calculated in 999 permutations (Fst), and 10,000 (IA). Due to limited sampling size S. boulardii and beetle gut isolates were not included in HWE, Fst, π and IA calculations. An UPGMA tree based on a pair-wise genetic distances between isolates included in PCA was build in PHYLIP 3.68 [64]. Haplotype networks were inferred from MLS1, ACT1, ADP1, PHD1 and RPB1 in PAUP* 4.0b10 applying the parsimony criterion, conducting heuristic searches, and using tree-bisection-reconnection (TBR) as branch swapping algorithm. These five loci were chosen for analysis based on their prior employment in phylogenetic studies of the Saccharomycetales, population genetic analysis of S. cerevisiae and C. albicans and their differential expression patterns upon phagocytosis by macrophages (ACT1, RPB1 [37], ADP1 [65], MLS1 [66], PHD1 [16]).

Survival during oxidative stress

Two-day old cultures from YPD plates were inoculated into liquid synthetic defined media (SD, 2% glucose, 37.8 mM (NH4)2SO4, 1.7 g/l yeast nitrogen base (YNB) without amino acids and ammonium sulfate), incubated at 30°C while shaking (250 rpm), transferred once, and grown over night. In order to reduce variability in the cell suspension, the cells were washed twice with 0.9% NaCl solution and dissolved in 1× phosphate buffered saline. After adjusting the cell number to 2×103 cells/ml, 20 mM tert-Butyl hydroperoxide (TBHP, Sigma) were added and treated and control samples incubated for one hour at 30°C, while shaking. 100 µl cell suspensions from both samples were spread on YPD plates, incubated for 48 h at 30°C and colonies counted. Survival was calculated as the ratio of treated over untreated cells. Each strain was tested three times and results were clustered by origin and plotted in box-and-whisker plots. A one-way ANOVA with Bonferroni correction was carried out on experimental triplicates of each strain in every group to assess differences in variation for all groups with N≥10.

Fluorescence-Activated Cell Sorting (FACS)

The ploidy of each S. cerevisiae strain was assessed using a previously published procedure with two modifications [67]. Exponentially growing cells were used and the final cell number was adjusted to 5×107 cells/ml. The successful completion of the staining procedure was verified microscopically and cells were stored in the dark at 4°C until ready for FACS analysis. On average 9,860 gated events were measured per strain.

Supporting Information

Table S1.

Genbank accession numbers for haplotype sequences.

https://doi.org/10.1371/journal.pone.0005317.s001

(0.04 MB DOC)

Figure S1.

UPGMA tree. The tree was generated from a pair-wise genetic distance matrix based on haplotype data of the 87 strains that were included in PCA. The strains are color-coded by origin (legend) and groups recognized in PCA indicated at internodes in the tree. Strains marked with * denote isolates that are not included in PCA groups A, B or C confidence envelopes. The numbered bar below the tree indicates total genetic distance observed in the data set.

https://doi.org/10.1371/journal.pone.0005317.s002

(1.34 MB TIF)

Acknowledgments

John McCusker (clinic, S. boulardii), Cletus Kurtzman (fruit), Brad Nicholson (brewery), Paul Sniegowksi (soil, PA), Anne Rouse (vineyards, NC and AU), Wiley Schell, (clinic), and Meredith Blackwell (insect gut) contributed S. cerevisiae strains. Heath O'Brien and Suzanne Joneson assisted in collecting wild yeast in NC. Andria Allen and Lisa Bukovnik assisted with the generation of the DNA sequence data and the FACS facility at the Duke Comprehensive Cancer Center conducted FACS analyses. Thomas Mitchell, Rytas Vilgalys and Joseph Heitman provided insightful comments and suggestions on the manuscript.

Author Contributions

Conceived and designed the experiments: SD. Performed the experiments: SD. Analyzed the data: SD. Contributed reagents/materials/analysis tools: FSD. Wrote the paper: SD.

References

  1. 1. McGovern PE, Zhang J, Tang J, Zhang Z, Hall GR, et al. (2004) Fermented beverages of pre- and proto-historic China. Proceedings of the National Academy of Sciences 101: 17593–17598.
  2. 2. Cavalieri D, McGovern PE, Hartl DL, Mortimer R, Polsinelli M (2003) Evidence for S. cerevisiae fermentation in ancient wine. Journal of Molecular Evolution 57: S226–S232.
  3. 3. Fay JC, Benavides JA (2005) Evidence for domesticated and wild populations of Saccharomyces cerevisiae. PLoS Genetics 1: e5.
  4. 4. Legras J-L, Merdinoglu D, Cornuet J-M, Karst F (2007) Bread, beer and wine: Saccharomyces cerevisiae diversity reflects human history. Molecular Ecology 16: 2091–2102.
  5. 5. Vaughan-Martini A, Martini A (1995) Facts, myths and legends on the prime industrial microorganism. Journal of Industrial Microbiology 14: 514–522.
  6. 6. Capriotti A (1954) Yeasts in some Netherlands soils. Antonie van Leeuwenhoek 21: 145–156.
  7. 7. Capriotti A (1967) Yeasts from U.S.A. soils. Archiv für Mikrobiologie 57: 406–413.
  8. 8. Naumov GI, Naumova ES, Korhola M (1992) Genetic identification of natural Saccharomyces sensu stricto yeasts from Finland, Holland and Slovakia. Antonie van Leeuwenhoek 61: 237–243.
  9. 9. Sniegowski PD, Dombrowski PG, Fingerman E (2002) Saccharomyces cerevisiae and Saccharomyces paradoxus coexist in a natural woodland site in North America and display different levels of reproductive isolation from European conspecifics. FEMS Yeast Research 1: 299–306.
  10. 10. Naumov GI, Naumova ES, Sniegowski PD (1998) Saccharomyces paradoxus and Saccharomyces cerevisiae are associated with exudates of North American oaks. Canadian Journal of Microbiology 44: 1045–1050.
  11. 11. Goddard MR, Burt A (1999) Recurrent invasion and extinction of a selfish gene. Proceedings of the National Academy of Sciences 96: 13880–13885.
  12. 12. Suh S-O, McHugh JV, Pollock DD, Blackwell M (2005) The beetle gut: a hyperdiverse source of novel yeasts. Mycological Research 109: 261–265.
  13. 13. Garza S, Teixidó JA, Sanchis V, Viñas I, Condón S (1994) Heat resistance of Saccharomyces cerevisiae strains isolated from spoiled peach puree. International Journal of Food Microbiology 23: 209–213.
  14. 14. Toraskar MV, Modi VV (1989) Studies on spoilage of Alphonso Mangoes by Saccharomyces cerevisiae. Acta Horticulturae 231: 697–708.
  15. 15. Replansky T, Koufopanou V, Greig D, Bell G (2008) Saccharomyces sensu stricto as a model system for evolution and ecology. Trends in Ecology and Evolution 23: 494–501.
  16. 16. Aa E, Townsend JP, Adams RI, Nielsen KM, Taylor JW (2006) Population structure and gene evolution in Saccharomyces cerevisiae. FEMS Yeast Research 6: 702–715.
  17. 17. Skovgaard N (2007) New trends in emerging pathogens. International Journal of Food Microbiology 120: 217–224.
  18. 18. Hazen KC (1995) New and emerging yeast pathogens. Clinical Microbiology Reviews 8: 462–478.
  19. 19. Cimolai N, Gill MJ, Church D (1987) Saccharomyces cerevisiae fungemia: case report and review of the literature. Diagnostic Microbiology and Infectious Disease 8: 113–117.
  20. 20. Sobel JD, Vazquez J, Lynch M, Meriwether C, Zervos MJ (1993) Vaginitis due to Saccharomyces cerevisiae: epidemiology, clinical aspects, and therapy. Clinical Infectious Diseases 16: 93–99.
  21. 21. McCullough MJ, Clemons KV, Farina C, McCusker JH, Stevens DA (1998) Epidemiological investigation of vaginal Saccharomyces cerevisiae isolates by a genotypic method. Journal of Clinical Microbiology 36: 557–562.
  22. 22. Aucott JN, Fayen J, Grossnicklas H, Morrissey A, Lederman MM, et al. (1990) Invasive infection with Saccharomyces cerevisiae: report of three cases and review. Reviews of Infectious Diseases 12: 406–411.
  23. 23. Chertow GM, Marcantonio ER, Wells RG (1991) Saccharomyces cerevisiae empyema in a patient with esophago-pleural fistula complicating variceal sclerotherapy. Chest 99: 1518–1519.
  24. 24. Konecny P, Drummond FM, Tish KN, Tapsall JW (1999) Saccharomyces cerevisiae oesophagitis in an HIV-infected patient. International Journal of STD & AIDS 10: 821–822.
  25. 25. Tawfik OW, Papasian CJ, Dixon A, Potter LM (1989) Saccharomyces cerevisiae pneumonia in a patient with acquired immune deficiency syndrome. Journal of Clinical Microbiology 27: 1689–1691.
  26. 26. Smith D, Metzgar D, Wills C, Fierer J (2002) Fatal Saccharomyces cerevisiae aortic graft infection. Journal of Clinical Microbiology 40: 2691–2692.
  27. 27. Williams JS, Mufti GJ, Powell S, Salisbury JR, Higgins EM (2007) Saccharomyces cerevisiae emboli in an immunocompromised patient with relapsed acute myeloid leukaemia. Clinical and Experimental Dermatology 32: 395–397.
  28. 28. Cesaro S, Chinello P, Rossi L, Zanesco L (2000) Saccharomyces cerevisiae fungemia in a neutropenic patient treated with Saccharomyces boulardii. Support Care Cancer 8: 504–505.
  29. 29. Bassetti S, Frei R, Zimmerli W (1998) Fungemia with Saccharomyces cerevisiae after treatment with Saccharomyces boulardii. The American Journal of Medicine 105: 71–72.
  30. 30. Enache-Angoulvant A, Hennequin C (2005) Invasive Saccharomyces infection: a comprehensive review. Clinical Infectious Diseases 41: 1559–1568.
  31. 31. McFarland LV (1996) Saccharomyces boulardii is not Saccharomyces cerevisiae. Clinical Infectious Diseases 22: 200–201.
  32. 32. Chia JKS, Chan SM, Goldstein H (1995) Baker's yeast as adjunctive therapy for relapses of Clostridium difficile diarrhea. Clinical Infectious Diseases 20: 1581.
  33. 33. McCullough MJ, Clemons KV, McCusker JH, Stevens DA (1998) Species identification and virulence attributes of Saccharomyces boulardii (nom. inval.). Journal of Clinical Microbiology 36: 2613–2617.
  34. 34. McNeil MM, Nash SL, Hajjeh RA, Phelan MA, Conn LA, et al. (2001) Trends in mortality due to invasive mycotic diseases in the United States, 1980–1997. Clin Inf Dis 33: 641–647.
  35. 35. Zaoutis TE, Argon J, Chu J, Berlin JA, Walsh TJ, et al. (2005) The epidemiology and attributable outcomes of candidemia in adults and children hospitalized in the United States: a propensity analysis. Clinical Infectious Diseases 41: 1232–1239.
  36. 36. McCusker JH (2006) Saccharomyces cerevisiae - an emerging and model pathogenic fungus. In: Heitman J, Filler SG, Edwards JE Jr, Mitchell AP, editors. Molecular Principles of Fungal Pathogenesis. ASM Press. pp. 245–260.
  37. 37. Diezmann S, Cox CJ, Schönian G, Vilgalys RJ, Mitchell TG (2004) Phylogeny and evolution of medical species of Candida and related taxa: a multigenic analysis. Journal of Clinical Microbiology 42: 5624–5635.
  38. 38. Rea RB, Gahan CG, Hill C (2004) Disruption of putative regulatory loci in Listeria monocytogenes demonstrates a significant role for Fur and PerR in virulence. Infection and Immunity 72: 717–727.
  39. 39. Riboulet E, Verneuil N, La Carbona S, Sauvageot N, Auffray Y, et al. (2007) Relationships between oxidative stress response and virulence in Enterococcus faecalis. Journal of Molecular Microbiology and Biotechnology 13: 140–146.
  40. 40. Abshire KZ, Neidhardt FC (1993) Analysis of proteins synthesized by Salmonella thyphimurium during growth within a host macrophage. Journal of Bacteriology 175: 3734–3743.
  41. 41. Iyer GYN, Islam MF, Quastel JH (1961) Biochemical aspects of phagocytosis. Nature 192: 535–541.
  42. 42. Sbarra AJ, Karnovsky ML (1959) The biochemical basis of phagocytosis. I. metabolic changes during the ingestion of particles by polymorphnuclear leukocytes. Journal of Biological Chemistry 234: 1355–1362.
  43. 43. McCusker JH, Clemons KV, Stevens DA, Davis RW (1994) Genetic characterization of pathogenic Saccharomyces cerevisiae isolates. Genetics 136: 1261–1269.
  44. 44. Clemons KV, Hanson LC, Stevens DA (1996) Colony phenotype switching in clinical and non-clinical isolates of Saccharomyces cerevisiae. Journal of Medical and Veterinary Mycology 34: 259–264.
  45. 45. McCusker JH, Clemons KV, Stevens DA, Davis RW (1994) Saccharomyces cerevisiae virulence phenotype as determined with CD-1 mice is associated with the ability to grow at 42°C and form pseudohyphae. Infection and Immunity 62: 5447–5455.
  46. 46. Weger SD, Ganji A, Clemons KV, Byron JK, Minn Y, et al. (2002) Correlation of the frequency of petite formation by isolates of Saccharomyces cerevisiae with virulence. Medical Mycology 40: 161–168.
  47. 47. Aref H, Cruess WV (1934) An investigation of the thermal death point of Saccharomyces ellipsoideus. Journal of Bacteriology 27: 443–452.
  48. 48. Rouse AL (2007) Evolution of a yeast transcriptional regulatory network: genetic, phenotypic and fitness variation. Durham, NC: Duke University.
  49. 49. Stephens M, Donnelly P (2003) A comparison of Bayesian methods for haplotype reconstruction. American Journal of Human Genetics 73: 1162–1168.
  50. 50. Stephens M, Smith NJ, Donnelly P (2001) A new statistical method for haplotype reconstruction from population data. American Journal of Human Genetics 68: 978–989.
  51. 51. Templeton AR, Crandall KA, Sing CF (1992) A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. cladogram estimation. Genetics 132: 619–633.
  52. 52. Steinmetz LM, Sinha H, Richards DR, Spiegelman JI, Oefner PJ, et al. (2002) Dissecting the architecture of a quantitative trait locus. Nature 416: 326–330.
  53. 53. Perlstein EO, Ruderfer DM, Ramachandran G, Haggarty SJ, Kruglyak L, et al. (2006) Revealing complex traits with small molecules and naturally recombinant yeast strains. Chemistry & Biology 13: 319–327.
  54. 54. Kvitek DJ, Will JL, Gasch AP (2008) Variations in stress sensitivity and genomic expression in diverse S. cerevisiae isolates. PLoS Genetics 4: e1000223.
  55. 55. Gardes M, Bruns TD (1993) ITS primers with enhanced specificity for basidiomycetes - application to the identification of mycorrhizae and rusts. Molecular Ecology 2: 113–118.
  56. 56. White TJ, Bruns T, Lee S, Taylor JW (1990) Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ, editors. PCR Protocols: A Guide to Methods and Applications. New York: Academic Press, Inc. pp. 315–322.
  57. 57. Werle E, Schneider C, Renner M, Völker M, Fiehn W (1994) Convenient single-step, one tube purification of PCR products for direct sequencing. Nucleic Acids Research 22: 4354–4355.
  58. 58. Maddison WP, Maddison DR (1997) MacClade. Analysis of phylogeny and character evolution. Sunderland, Mass.: Sinauer Associates, Inc.
  59. 59. Altschul SF, Gish W, Miller WG, Myers EW, Lipman DJ (1990) Basic Local Alignment Search Tool. Journal of Molecular Biology 215: 403–410.
  60. 60. Peakall R, Smouse PE (2006) GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes 6: 288–295.
  61. 61. Rozas J, Sánchez-Delbarrio JC, Messeguer X, Rozas R (2003) DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19: 2496–2497.
  62. 62. Agapow P-M, Burt A (2001) Indices of multilocus linkage disequilibrium. Molecular Ecology Notes 1: 101–102.
  63. 63. Esseiva P, Anglada F, Dujourdy L, Taroni F, Margot P, et al. (2005) Chemical profiling and classification of illicit heroin by principal component analysis, calculation of inter sample correlation and artificial neural networks. Talanta 67: 360–367.
  64. 64. Felsenstein J (1989) PHYLIP - Phylogeny Inference Package (Version 3.2). Cladistics 5: 164–166.
  65. 65. Bougnoux ME, Tavanti A, Bouchier C, Gow NAR, Magnier A, et al. (2003) Collaborative Consensus for Optimized Multilocus Sequence Typing of Candida albicans. Journal of Clinical Microbiology 41: 5265–5266.
  66. 66. Lorenz MC, Fink GR (2001) The glyoxylate cycle is required for fungal virulence. Nature 412: 83–86.
  67. 67. Burke D, Dawson D, Stearns T (2000) Techniques & Protocols #18. Methods in Yeast Genetics: A Cold Spring Harbor Laboratory Course Manual. Cold Spring Harbor: Cold Spring Harbor Laboratory Press.
  68. 68. Bulmer GS, Marquez ML, Co-Barcelona L, Fromtling RA (1999) Yeasts and fluconazole susceptibility in the Philippines. Mycopathologia 146: 117–120.
  69. 69. Clemons KV, McCusker JH, Davis RW, Stevens DA (1994) Comparative pathogenesis of clinical and nonclinical isolates of Saccharomyces cerevisiae. The Journal of Infectious Diseases 169: 859–867.
  70. 70. Mortimer RK, Johnston JR (1986) Genealogy of principal strains of the yeast genetic stock center. Genetics 113: 35–43.
  71. 71. Mortimer RK, Romano P, Suzzi G, Polsinelli M (1994) Genome renewal: a new phenomenon revealed from a genetic study of 43 strains of Saccharomyces cerevisiae derived from natural fermentation of grape musts. Genetics 10: 1543–1552.
  72. 72. Matheny PB, Liu YJ, Ammirati JF, Hall BD (2002) Using RPB1 sequences to improve phylogenetic inference among mushrooms (Inocybe, Agaricales). American Journal of Botany 89: 688–698.