Sporulation environment drives phenotypic variation in the pathogen Aspergillus fumigatus.

Aspergillus fumigatus causes more than 300,000 life-threatening infections annually and is widespread across varied environments with a single colony producing thousands of conidia, genetically identical dormant spores. Conidia are easily wind-dispersed to new environments where they can germinate and, if inhaled by susceptible hosts, cause disease. Using high-throughput single-cell analysis via flow cytometry we analyzed conidia produced and germinated in nine environmentally and medically relevant conditions (complete medium, minimal medium, high temperature, excess copper, excess iron, limited iron, excess salt, excess reactive oxygen species, and limited zinc). We found that germination phenotypes vary among genetically identical individuals, that the environment of spore production determines the size of spores and the degree of germination heterogeneity, and that the environment of spore production impacts virulence in a Galleria mellonella host.

Genetic variation in aneuploidy prevalence and tolerance across Saccharomyces cerevisiae lineages.

Individuals carrying an aberrant number of chromosomes can vary widely in their expression of aneuploidy phenotypes. A major unanswered question is the degree to which an individual's genetic makeup influences its tolerance of karyotypic imbalance. Here we investigated within-species variation in aneuploidy prevalence and tolerance, using Saccharomyces cerevisiae as a model for eukaryotic biology. We analyzed genotypic and phenotypic variation recently published for over 1,000 S. cerevisiae strains spanning dozens of genetically defined clades and ecological associations. Our results show that the prevalence of chromosome gain and loss varies by clade and can be better explained by differences in genetic background than ecology. The relationships between lineages with high aneuploidy frequencies suggest that increased aneuploidy prevalence emerged multiple times in S. cerevisiae evolution. Separate from aneuploidy prevalence, analyzing growth phenotypes revealed that some genetic backgrounds-such as the European Wine lineage-show fitness costs in aneuploids compared to euploids, whereas other clades with high aneuploidy frequencies show little evidence of major deleterious effects. Our analysis confirms that chromosome gain can produce phenotypic benefits, which could influence evolutionary trajectories. These results have important implications for understanding genetic variation in aneuploidy prevalence in health, disease, and evolution.

Evolution of Ty1 copy number control in yeast by horizontal transfer and recombination.

Transposable elements constitute a large fraction of most eukaryotic genomes. Insertion of mobile DNA sequences typically has deleterious effects on host fitness, and thus diverse mechanisms have evolved to control mobile element proliferation. Mobility of the Ty1 retrotransposon in Saccharomyces yeasts is regulated by copy number control (CNC) mediated by a self-encoded restriction factor derived from the Ty1 gag capsid gene that inhibits virus-like particle function. Here, we survey a panel of wild and human-associated strains of S. cerevisiae and S. paradoxus to investigate how genomic Ty1 content influences variation in Ty1 mobility. We observe high levels of mobility for a tester element with a gag sequence from the canonical Ty1 subfamily in permissive strains that either lack full-length Ty1 elements or only contain full-length copies of the Ty1' subfamily that have a divergent gag sequence. In contrast, low levels of canonical Ty1 mobility are observed in restrictive strains carrying full-length Ty1 elements containing a canonical gag sequence. Phylogenomic analysis of full-length Ty1 elements revealed that Ty1' is the ancestral subfamily present in wild strains of S. cerevisiae, and that canonical Ty1 in S. cerevisiae is a derived subfamily that acquired gag from S. paradoxus by horizontal transfer and recombination. Our results provide evidence that variation in the ability of S. cerevisiae and S. paradoxus strains to repress canonical Ty1 transposition via CNC is regulated by the genomic content of different Ty1 subfamilies, and that self-encoded forms of transposon control can spread across species boundaries by horizontal transfer.

Phased Diploid Genome Assemblies for Three Strains of Candida albicans from Oak Trees.

Although normally a harmless commensal, Candida albicans, it is also one of the most common causes of bloodstream infections in the U.S. Candida albicans has long been considered an obligate commensal, however, recent studies suggest it can live outside animal hosts. Here, we have generated PacBio sequences and phased genome assemblies for three C. albicans strains from oak trees (NCYC 4144, NCYC 4145, and NCYC 4146). PacBio datasets are high depth (over 400 fold coverage) and more than half of the sequencing data are contained in reads longer than 15 kb. Primary assemblies showed high contiguity with several chromosomes for each strain recovered as single contigs, and greater than half of the alternative haplotype sequence was assembled in haplotigs at least 174 kb long. Using these assemblies we were able to identify structural polymorphisms, including a polymorphic inversion over 100 kb in length. These results show that phased de novo diploid assemblies for C. albicans can enable the study of genomic variation within and among strains of an important fungal pathogen.

mSphere of Influence: the Wild Genetic Diversity of Our Closest Yeast Companions.

Douda Bensasson uses the population genomics of model yeast species to understand how wild yeast colonize new environments, such as humans or their food. In this mSphere of Influence article, she reflects on how the discovery of "Surprisingly diverged populations of Saccharomyces cerevisiae in natural environments remote from human activity" (Q.-M. Wang, W.-Q. Liu, G. Liti, S.-A. Wang, and F.-Y. Bai, Mol Ecol 21:5404-5417, 2012, https://doi.org/10.1111/j.1365-294X.2012.05732.x) showed that a field survey and population genetic analysis of old growth forests could "unveil the hidden part of the iceberg" of natural variation in S. cerevisiae that went unnoticed for over a hundred years of yeast research.

Diverse Lineages of Candida albicans Live on Old Oaks.

The human pathogen Candida albicans is considered an obligate commensal of animals, yet it is occasionally isolated from trees, shrubs, and grass. We generated genome sequence data for three strains of C. albicans that we isolated from oak trees in an ancient wood pasture, and compared these to the genomes of over 200 clinical strains. C. albicans strains from oak are similar to clinical C. albicans in that they are predominantly diploid and can become homozygous at the mating locus through whole-chromosome loss of heterozygosity. Oak strains differed from clinical strains in showing slightly higher levels of heterozygosity genome-wide. Using phylogenomic analyses and in silico chromosome painting, we show that each oak strain is more closely related to strains from humans and other animals than to strains from other oaks. The high genetic diversity of C. albicans from old oaks shows that they can live in this environment for extended periods of time.

Habitat Predicts Levels of Genetic Admixture in Saccharomyces cerevisiae.

Genetic admixture can provide material for populations to adapt to local environments, and this process has played a crucial role in the domestication of plants and animals. The model yeast, Saccharomyces cerevisiae, has been domesticated multiple times for the production of wine, sake, beer, and bread, but the high rate of admixture between yeast lineages has so far been treated as a complication for population genomic analysis. Here, we make use of the low recombination rate at centromeres to investigate admixture in yeast using a classic Bayesian approach and a locus-by-locus phylogenetic approach. Using both approaches, we find that S. cerevisiae from stable oak woodland habitats are less likely to show recent genetic admixture compared with those isolated from transient habitats such as fruits, wine, or human infections. When woodland yeast strains do show recent genetic admixture, the degree of admixture is lower than in strains from other habitats. Furthermore, S. cerevisiae populations from oak woodlands are genetically isolated from each other, with only occasional migration between woodlands and local fruit habitats. Application of the phylogenetic approach suggests that there is a previously undetected population in North Africa that is the closest outgroup to the European S. cerevisiae, including the domesticated Wine population. Careful testing for admixture in S. cerevisiae leads to a better understanding of the underlying population structure of the species and will be important for understanding the selective processes underlying domestication in this economically important species.

Adaptive divergence in wine yeasts and their wild relatives suggests a prominent role for introgressions and rapid evolution at noncoding sites.

In Saccharomyces cerevisiae, the main yeast in wine fermentation, the opportunity to examine divergence at the molecular level between a domesticated lineage and its wild counterpart arose recently due to the identification of the closest relatives of wine strains, a wild population associated with Mediterranean oaks. As genomic data are available for a considerable number of representatives belonging to both groups, we used population genomics to estimate the degree and distribution of nucleotide variation between wine yeasts and their closest wild relatives. We found widespread genomewide divergence, particularly at noncoding sites, which, together with above average divergence in trans-acting DNA binding proteins, may suggest an important role for divergence at the level of transcriptional regulation. Nine outlier regions putatively under strong divergent selection were highlighted by a genomewide scan under stringent conditions. Several cases of introgressions, originating in the sibling species Saccharomyces paradoxus, were also identified in the Mediterranean oak population. FZF1 and SSU1, mostly known for conferring sulphite resistance in wine yeasts, were among the introgressed genes, although not fixed. Because the introgressions detected in our study are not found in wine strains, we hypothesize that ongoing divergent ecological selection segregates the two forms between the different niches. Together, our results provide a first insight into the extent and kind of divergence between wine yeasts and their closest wild relatives.

Summer temperature can predict the distribution of wild yeast populations.

The wine yeast, Saccharomyces cerevisiae, is the best understood microbial eukaryote at the molecular and cellular level, yet its natural geographic distribution is unknown. Here we report the results of a field survey for S. cerevisiae,S. paradoxus and other budding yeast on oak trees in Europe. We show that yeast species differ in their geographic distributions, and investigated which ecological variables can predict the isolation rate of S. paradoxus, the most abundant species. We find a positive association between trunk girth and S. paradoxus abundance suggesting that older trees harbor more yeast. S. paradoxus isolation frequency is also associated with summer temperature, showing highest isolation rates at intermediate temperatures. Using our statistical model, we estimated a range of summer temperatures at which we expect high S. paradoxus isolation rates, and show that the geographic distribution predicted by this optimum temperature range is consistent with the worldwide distribution of sites where S. paradoxus has been isolated. Using laboratory estimates of optimal growth temperatures for S. cerevisiae relative to S. paradoxus, we also estimated an optimum range of summer temperatures for S. cerevisiae. The geographic distribution of these optimum temperatures is consistent with the locations where wild S. cerevisiae have been reported, and can explain why only human-associated S. cerevisiae strains are isolated at northernmost latitudes. Our results provide a starting point for targeted isolation of S. cerevisiae from natural habitats, which could lead to a better understanding of climate associations and natural history in this important model microbe.

A population genomics insight into the Mediterranean origins of wine yeast domestication.

The domestication of the wine yeast Saccharomyces cerevisiae is thought to be contemporary with the development and expansion of viticulture along the Mediterranean basin. Until now, the unavailability of wild lineages prevented the identification of the closest wild relatives of wine yeasts. Here, we enlarge the collection of natural lineages and employ whole-genome data of oak-associated wild isolates to study a balanced number of anthropic and natural S. cerevisiae strains. We identified industrial variants and new geographically delimited populations, including a novel Mediterranean oak population. This population is the closest relative of the wine lineage as shown by a weak population structure and further supported by genomewide population analyses. A coalescent model considering partial isolation with asymmetrical migration, mostly from the wild group into the Wine group, and population growth, was found to be best supported by the data. Importantly, divergence time estimates between the two populations agree with historical evidence for winemaking. We show that three horizontally transmitted regions, previously described to contain genes relevant to wine fermentation, are present in the Wine group but not in the Mediterranean oak group. This represents a major discontinuity between the two populations and is likely to denote a domestication fingerprint in wine yeasts. Taken together, these results indicate that Mediterranean oaks harbour the wild genetic stock of domesticated wine yeasts.

Evolutionary genomics of transposable elements in Saccharomyces cerevisiae.

Saccharomyces cerevisiae is one of the premier model systems for studying the genomics and evolution of transposable elements. The availability of the S. cerevisiae genome led to unprecedented insights into its five known transposable element families (the LTR retrotransposons Ty1-Ty5) in the years shortly after its completion. However, subsequent advances in bioinformatics tools for analysing transposable elements and the recent availability of genome sequences for multiple strains and species of yeast motivates new investigations into Ty evolution in S. cerevisiae. Here we provide a comprehensive phylogenetic and population genetic analysis of all Ty families in S. cerevisiae based on a systematic re-annotation of Ty elements in the S288c reference genome. We show that previous annotation efforts have underestimated the total copy number of Ty elements for all known families. In addition, we identify a new family of Ty3-like elements related to the S. paradoxus Ty3p which is composed entirely of degenerate solo LTRs. Phylogenetic analyses of LTR sequences identified three families with short-branch, recently active clades nested among long branch, inactive insertions (Ty1, Ty3, Ty4), one family with essentially all recently active elements (Ty2) and two families with only inactive elements (Ty3p and Ty5). Population genomic data from 38 additional strains of S. cerevisiae show that the majority of Ty insertions in the S288c reference genome are fixed in the species, with insertions in active clades being predominantly polymorphic and insertions in inactive clades being predominantly fixed. Finally, we use comparative genomic data to provide evidence that the Ty2 and Ty3p families have arisen in the S. cerevisiae genome by horizontal transfer. Our results demonstrate that the genome of a single individual contains important information about the state of TE population dynamics within a species and suggest that horizontal transfer may play an important role in shaping the genomic diversity of transposable elements in unicellular eukaryotes.

Evidence for a high mutation rate at rapidly evolving yeast centromeres.

BACKGROUND: Although their role in cell division is essential, centromeres evolve rapidly in animals, plants and yeasts. Unlike the complex centromeres of plants and aminals, the point centromeres of Saccharomcyes yeasts can be readily sequenced to distinguish amongst the possible explanations for fast centromere evolution. RESULTS: Using DNA sequences of all 16 centromeres from 34 strains of Saccharomyces cerevisiae and population genomic data from Saccharomyces paradoxus, I show that centromeres in both species evolve 3 times more rapidly even than selectively unconstrained DNA. Exceptionally high levels of polymorphism seen in multiple yeast populations suggest that rapid centromere evolution does not result from the repeated selective sweeps expected under meiotic drive. I further show that there is little evidence for crossing-over or gene conversion within centromeres, although there is clear evidence for recombination in their immediate vicinity. Finally I show that the mutation spectrum at centromeres is consistent with the pattern of spontaneous mutation elsewhere in the genome. CONCLUSIONS: These results indicate that rapid centromere evolution is a common phenomenon in yeast species. Furthermore, these results suggest that rapid centromere evolution does not result from the mutagenic effect of gene conversion, but from a generalised increase in the mutation rate, perhaps arising from the unusual chromatin structure at centromeres in yeast and other eukaryotes.

Population genomics of domestic and wild yeasts.

Since the completion of the genome sequence of Saccharomyces cerevisiae in 1996 (refs 1, 2), there has been a large increase in complete genome sequences, accompanied by great advances in our understanding of genome evolution. Although little is known about the natural and life histories of yeasts in the wild, there are an increasing number of studies looking at ecological and geographic distributions, population structure and sexual versus asexual reproduction. Less well understood at the whole genome level are the evolutionary processes acting within populations and species that lead to adaptation to different environments, phenotypic differences and reproductive isolation. Here we present one- to fourfold or more coverage of the genome sequences of over seventy isolates of the baker's yeast S. cerevisiae and its closest relative, Saccharomyces paradoxus. We examine variation in gene content, single nucleotide polymorphisms, nucleotide insertions and deletions, copy numbers and transposable elements. We find that phenotypic variation broadly correlates with global genome-wide phylogenetic relationships. S. paradoxus populations are well delineated along geographic boundaries, whereas the variation among worldwide S. cerevisiae isolates shows less differentiation and is comparable to a single S. paradoxus population. Rather than one or two domestication events leading to the extant baker's yeasts, the population structure of S. cerevisiae consists of a few well-defined, geographically isolated lineages and many different mosaics of these lineages, supporting the idea that human influence provided the opportunity for cross-breeding and production of new combinations of pre-existing variations.

Rapid evolution of yeast centromeres in the absence of drive.

To find the most rapidly evolving regions in the yeast genome we compared most of chromosome III from three closely related lineages of the wild yeast Saccharomyces paradoxus. Unexpectedly, the centromere appears to be the fastest-evolving part of the chromosome, evolving even faster than DNA sequences unlikely to be under selective constraint (i.e., synonymous sites after correcting for codon usage bias and remnant transposable elements). Centromeres on other chromosomes also show an elevated rate of nucleotide substitution. Rapid centromere evolution has also been reported for some plants and animals and has been attributed to selection for inclusion in the egg or the ovule at female meiosis. But Saccharomyces yeasts have symmetrical meioses with all four products surviving, thus providing no opportunity for meiotic drive. In addition, yeast centromeres show the high levels of polymorphism expected under a neutral model of molecular evolution. We suggest that yeast centromeres suffer an elevated rate of mutation relative to other chromosomal regions and they change through a process of "centromere drift," not drive.

Population genomics of the wild yeast Saccharomyces paradoxus: Quantifying the life cycle.

Most microbes have complex life cycles with multiple modes of reproduction that differ in their effects on DNA sequence variation. Population genomic analyses can therefore be used to estimate the relative frequencies of these different modes in nature. The life cycle of the wild yeast Saccharomyces paradoxus is complex, including clonal reproduction, outcrossing, and two different modes of inbreeding. To quantify these different aspects we analyzed DNA sequence variation in the third chromosome among 20 isolates from two populations. Measures of mutational and recombinational diversity were used to make two independent estimates of the population size. In an obligately sexual population these values should be approximately equal. Instead there is a discrepancy of about three orders of magnitude between our two estimates of population size, indicating that S. paradoxus goes through a sexual cycle approximately once in every 1,000 asexual generations. Chromosome III also contains the mating type locus (MAT), which is the most outbred part in the entire genome, and by comparing recombinational diversity as a function of distance from MAT we estimate the frequency of matings to be approximately 94% from within the same tetrad, 5% with a clonemate after switching the mating type, and 1% outcrossed. Our study illustrates the utility of population genomic data in quantifying life cycles.

Release and persistence of extracellular DNA in the environment.

The introduction of genetically modified organisms (GMOs) has called for an improved understanding of the fate of DNA in various environments, because extracellular DNA may also be important for transferring genetic information between individuals and species. Accumulating nucleotide sequence data suggest that acquisition of foreign DNA by horizontal gene transfer (HGT) is of considerable importance in bacterial evolution. The uptake of extracellular DNA by natural transformation is one of several ways bacteria can acquire new genetic information given sufficient size, concentration and integrity of the DNA. We review studies on the release, breakdown and persistence of bacterial and plant DNA in soil, sediment and water, with a focus on the accessibility of the extracellular nucleic acids as substrate for competent bacteria. DNA fragments often persist over time in many environments, thereby facilitating their detection and characterization. Nevertheless, the long-term physical persistence of DNA fragments of limited size observed by PCR and Southern hybridization often contrasts with the short-term availability of extracellular DNA to competent bacteria studied in microcosms. The main factors leading to breakdown of extracellular DNA are presented. There is a need for improved methods for accurately determining the degradation routes and the persistence, integrity and potential for horizontal transfer of DNA released from various organisms throughout their lifecycles.

Recent LTR retrotransposon insertion contrasts with waves of non-LTR insertion since speciation in Drosophila melanogaster.

LTR and non-LTR retrotransposons exhibit distinct patterns of abundance within the Drosophila melanogaster genome, yet the causes of these differences remain unknown. Here we investigate whether genomic differences between LTR and non-LTR retrotransposons reflect systematic differences in their insertion history. We find that for 17 LTR and 10 non-LTR retrotransposon families that evolve under a pseudogene-like mode of evolution, most elements from LTR families have integrated in the very recent past since colonization of non-African habitats ( approximately 16,000 years ago), whereas elements from non-LTR families have been accumulating in overlapping waves since the divergence of D. melanogaster from its sister species, Drosophila simulans ( approximately 5.4 Mya). LTR elements are significantly younger than non-LTR elements, individually and by family, in regions of high and low recombination, and in genic and intergenic regions. We show that analysis of transposable element (TE) nesting provides a method to calculate transposition rates from genome sequences, which we estimate to be one to two orders of magnitude lower than those that are based on mutation accumulation studies. Recent LTR integration provides a nonequilibrium alternative for the low population frequency of LTR elements in this species, a pattern that is classically interpreted as evidence for selection against the transpositional increase of TEs. Our results call for a new class of population genetic models that incorporate TE copy number, allele frequency, and the age of insertions to provide more powerful and robust inferences about the forces that control the evolution of TEs in natural populations.

Mitochondrial genome sequences and comparative genomics of Phytophthora ramorum and P. sojae.

The sequences of the mitochondrial genomes of the oomycetes Phytophthora ramorum and P. sojae were determined during the course of complete nuclear genome sequencing (Tyler et al., Science, 313:1261,2006). Both mitochondrial genomes are circular mapping, with sizes of 39,314 bp for P. ramorum and 42,977 bp for P. sojae. Each contains a total of 37 recognizable protein-encoding genes, 26 or 25 tRNAs (P. ramorum and P. sojae, respectively) specifying 19 amino acids, six more open reading frames (ORFs) that are conserved, presumably due to functional constraint, across Phytophthora species (P. sojae, P. ramorum, and P. infestans), six ORFs that are unique for P. sojae and one that is unique for P. ramorum. Non-coding regions comprise about 11.5 and 18.4% of the genomes of P. ramorum and P. sojae, respectively. Relative to P. sojae, there is an inverted repeat of 1,150 bp in P. ramorum that includes an unassigned unique ORF, a tRNA gene, and adjacent non-coding sequences, but otherwise the gene order in both species is identical. Comparisons of these genomes with published sequences of the P. infestans mitochondrial genome reveals a number of similarities, but the gene order in P. infestans differed in two adjacent locations due to inversions and specific regions of the genomes exhibited greater divergence than others. For example, the breakpoints for the inversions observed in P. infestans corresponded to regions of high sequence divergence in comparisons between P. ramorum and P. sojae and the location of a hypervariable microsatellite sequence (eight repeats of 24 bp) in the P. sojae genome corresponds to a site of major length variation in P. infestans. Although the overwhelming majority of each genome is conserved (81-92%), there are a number of genes that evolve more rapidly than others. Some of these rapidly evolving genes appear specific to Phytophthora, arose recently, and future evaluation of their function and the effects of their loss could prove fruitful for understanding the phylogeny of these devastating plant pathogens.

Transition-transversion bias is not universal: a counter example from grasshopper pseudogenes.

Comparisons of the DNA sequences of metazoa show an excess of transitional over transversional substitutions. Part of this bias is due to the relatively high rate of mutation of methylated cytosines to thymine. Postmutation processes also introduce a bias, particularly selection for codon-usage bias in coding regions. It is generally assumed, however, that there is a universal bias in favour of transitions over transversions, possibly as a result of the underlying chemistry of mutation. Surprisingly, this underlying trend has been evaluated only in two types of metazoan, namely Drosophila and the Mammalia. Here, we investigate a third group, and find no such bias. We characterize the point substitution spectrum in Podisma pedestris, a grasshopper species with a very large genome. The accumulation of mutations was surveyed in two pseudogene families, nuclear mitochondrial and ribosomal DNA sequences. The cytosine-guanine (CpG) dinucleotides exhibit the high transition frequencies expected of methylated sites. The transition rate at other cytosine residues is significantly lower. After accounting for this methylation effect, there is no significant difference between transition and transversion rates. These results contrast with reports from other taxa and lead us to reject the hypothesis of a universal transition/transversion bias. Instead we suggest fundamental interspecific differences in point substitution processes.

Phytophthora genome sequences uncover evolutionary origins and mechanisms of pathogenesis.

Draft genome sequences have been determined for the soybean pathogen Phytophthora sojae and the sudden oak death pathogen Phytophthora ramorum. Oömycetes such as these Phytophthora species share the kingdom Stramenopila with photosynthetic algae such as diatoms, and the presence of many Phytophthora genes of probable phototroph origin supports a photosynthetic ancestry for the stramenopiles. Comparison of the two species' genomes reveals a rapid expansion and diversification of many protein families associated with plant infection such as hydrolases, ABC transporters, protein toxins, proteinase inhibitors, and, in particular, a superfamily of 700 proteins with similarity to known oömycete avirulence genes.