Biomolecular condensates in fungi are tuned to function at specific temperatures.

Temperature can impact every reaction and molecular interaction essential to a cell. For organisms that cannot regulate their own temperature, a major challenge is how to adapt to temperatures that fluctuate unpredictability and on variable timescales. Biomolecular condensation offers a possible mechanism for encoding temperature-responsiveness and robustness into cell biochemistry and organization. To explore this idea, we examined temperature adaptation in a filamentous-growing fungus called Ashbya gossypii that engages biomolecular condensates containing the RNA-binding protein Whi3 to regulate mitosis and morphogenesis. We collected wild isolates of Ashbya that originate in different climates and found that mitotic asynchrony and polarized growth, which are known to be controlled by the condensation of Whi3, are temperature sensitive. Sequence analysis in the wild strains revealed changes to specific domains within Whi3 known to be important in condensate formation. Using an in vitro condensate reconstitution assay we found that temperature impacts the relative abundance of protein to RNA within condensates and that this directly impacts the material properties of the droplets. Finally, we found that exchanging Whi3 genes between warm and cold isolates was sufficient to rescue some, but not all, condensate-related phenotypes. Together these data demonstrate that material properties of Whi3 condensates are temperature sensitive, that these properties are important for function, and that sequence optimizes properties for a given climate.

RNA viruses, M satellites, chromosomal killer genes, and killer/nonkiller phenotypes in the 100-genomes S. cerevisiae strains.

We characterized previously identified RNA viruses (L-A, L-BC, 20S, and 23S), L-A-dependent M satellites (M1, M2, M28, and Mlus), and M satellite-dependent killer phenotypes in the Saccharomyces cerevisiae 100-genomes genetic resource population. L-BC was present in all strains, albeit in 2 distinct levels, L-BChi and L-BClo; the L-BC level is associated with the L-BC genotype. L-BChi, L-A, 20S, 23S, M1, M2, and Mlus (M28 was absent) were in fewer strains than the similarly inherited 2µ plasmid. Novel L-A-dependent phenotypes were identified. Ten M+ strains exhibited M satellite-dependent killing (K+) of at least 1 of the naturally M0 and cured M0 derivatives of the 100-genomes strains; in these M0 strains, sensitivities to K1+, K2+, and K28+ strains varied. Finally, to complement our M satellite-encoded killer toxin analysis, we assembled the chromosomal KHS1 and KHR1 killer genes and used naturally M0 and cured M0 derivatives of the 100-genomes strains to assess and characterize the chromosomal killer phenotypes.

Comparative analysis of RNA enrichment methods for preparation of Cryptococcus neoformans RNA sequencing libraries.

RNA sequencing (RNA-Seq) experiments focused on gene expression involve removal of ribosomal RNA (rRNA) because it is the major RNA constituent of cells. This process, called RNA enrichment, is done primarily to reduce cost: without rRNA removal, deeper sequencing must be performed to compensate for the sequencing reads wasted on rRNA. The ideal RNA enrichment method removes all rRNA without affecting other RNA in the sample. We tested the performance of three RNA enrichment methods on RNA isolated from Cryptococcus neoformans, a fungal pathogen of humans. We find that the RNase H depletion method is more efficient in depleting rRNA and more specific in recapitulating non-rRNA levels present in unenriched controls than the commonly-used Poly(A) isolation method. The RNase H depletion method is also more effective than the Ribo-Zero depletion method as measured by rRNA depletion efficiency and recapitulation of protein-coding RNA levels present in unenriched controls, while the Ribo-Zero depletion method more closely recapitulates annotated non-coding RNA (ncRNA) levels. Finally, we leverage these data to accurately map the C. neoformans mitochondrial rRNA genes, and also demonstrate that RNA-Seq data generated with the RNase H and Ribo-Zero depletion methods can be used to explore novel C. neoformans long non-coding RNA genes.

Mitochondrial Genome Variation Affects Multiple Respiration and Nonrespiration Phenotypes in Saccharomyces cerevisiae.

Mitochondrial genome variation and its effects on phenotypes have been widely analyzed in higher eukaryotes but less so in the model eukaryote Saccharomyces cerevisiae Here, we describe mitochondrial genome variation in 96 diverse S. cerevisiae strains and assess associations between mitochondrial genotype and phenotypes as well as nuclear-mitochondrial epistasis. We associate sensitivity to the ATP synthase inhibitor oligomycin with SNPs in the mitochondrially encoded ATP6 gene. We describe the use of iso-nuclear F1 pairs, the mitochondrial genome equivalent of reciprocal hemizygosity analysis, to identify and analyze mitochondrial genotype-dependent phenotypes. Using iso-nuclear F1 pairs, we analyze the oligomycin phenotype-ATP6 association and find extensive nuclear-mitochondrial epistasis. Similarly, in iso-nuclear F1 pairs, we identify many additional mitochondrial genotype-dependent respiration phenotypes, for which there was no association in the 96 strains, and again find extensive nuclear-mitochondrial epistasis that likely contributes to the lack of association in the 96 strains. Finally, in iso-nuclear F1 pairs, we identify novel mitochondrial genotype-dependent nonrespiration phenotypes: resistance to cycloheximide, ketoconazole, and copper. We discuss potential mechanisms and the implications of mitochondrial genotype and of nuclear-mitochondrial epistasis effects on respiratory and nonrespiratory quantitative traits.

2µ plasmid in Saccharomyces species and in Saccharomyces cerevisiae.

We determined that extrachromosomal 2µ plasmid was present in 67 of the Saccharomyces cerevisiae 100-genome strains; in addition to variation in the size and copy number of 2µ, we identified three distinct classes of 2µ. We identified 2µ presence/absence and class associations with populations, clinical origin and nuclear genotypes. We also screened genome sequences of S. paradoxus, S. kudriavzevii, S. uvarum, S. eubayanus, S. mikatae, S. arboricolus and S. bayanus strains for both integrated and extrachromosomal 2µ. Similar to S. cerevisiae, we found no integrated 2µ sequences in any S. paradoxus strains. However, we identified part of 2µ integrated into the genomes of some S. uvarum, S. kudriavzevii, S. mikatae and S. bayanus strains, which were distinct from each other and from all extrachromosomal 2µ. We identified extrachromosomal 2µ in one S. paradoxus, one S. eubayanus, two S. bayanus and 13 S. uvarum strains. The extrachromosomal 2µ in S. paradoxus, S. eubayanus and S. cerevisiae were distinct from each other. In contrast, the extrachromosomal 2µ in S. bayanus and S. uvarum strains were identical with each other and with one of the three classes of S. cerevisiae 2µ, consistent with interspecific transfer.

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.

Structures of naturally evolved CUP1 tandem arrays in yeast indicate that these arrays are generated by unequal nonhomologous recombination.

An important issue in genome evolution is the mechanism by which tandem duplications are generated from single-copy genes. In the yeast Saccharomyces cerevisiae, most strains contain tandemly duplicated copies of CUP1, a gene that encodes a copper-binding metallothionein. By screening 101 natural isolates of S. cerevisiae, we identified five different types of CUP1-containing repeats, as well as strains that only had one copy of CUP1. A comparison of the DNA sequences of these strains indicates that the CUP1 tandem arrays were generated by unequal nonhomologous recombination events from strains that had one CUP1 gene.

Cryptococcus gattii VGIII isolates causing infections in HIV/AIDS patients in Southern California: identification of the local environmental source as arboreal.

Ongoing Cryptococcus gattii outbreaks in the Western United States and Canada illustrate the impact of environmental reservoirs and both clonal and recombining propagation in driving emergence and expansion of microbial pathogens. C. gattii comprises four distinct molecular types: VGI, VGII, VGIII, and VGIV, with no evidence of nuclear genetic exchange, indicating these represent distinct species. C. gattii VGII isolates are causing the Pacific Northwest outbreak, whereas VGIII isolates frequently infect HIV/AIDS patients in Southern California. VGI, VGII, and VGIII have been isolated from patients and animals in the Western US, suggesting these molecular types occur in the environment. However, only two environmental isolates of C. gattii have ever been reported from California: CBS7750 (VGII) and WM161 (VGIII). The incongruence of frequent clinical presence and uncommon environmental isolation suggests an unknown C. gattii reservoir in California. Here we report frequent isolation of C. gattii VGIII MATα and MATa isolates and infrequent isolation of VGI MATα from environmental sources in Southern California. VGIII isolates were obtained from soil debris associated with tree species not previously reported as hosts from sites near residences of infected patients. These isolates are fertile under laboratory conditions, produce abundant spores, and are part of both locally and more distantly recombining populations. MLST and whole genome sequence analysis provide compelling evidence that these environmental isolates are the source of human infections. Isolates displayed wide-ranging virulence in macrophage and animal models. When clinical and environmental isolates with indistinguishable MLST profiles were compared, environmental isolates were less virulent. Taken together, our studies reveal an environmental source and risk of C. gattii to HIV/AIDS patients with implications for the >1,000,000 cryptococcal infections occurring annually for which the causative isolate is rarely assigned species status. Thus, the C. gattii global health burden could be more substantial than currently appreciated.

The reference genome sequence of Saccharomyces cerevisiae: then and now.

The genome of the budding yeast Saccharomyces cerevisiae was the first completely sequenced from a eukaryote. It was released in 1996 as the work of a worldwide effort of hundreds of researchers. In the time since, the yeast genome has been intensively studied by geneticists, molecular biologists, and computational scientists all over the world. Maintenance and annotation of the genome sequence have long been provided by the Saccharomyces Genome Database, one of the original model organism databases. To deepen our understanding of the eukaryotic genome, the S. cerevisiae strain S288C reference genome sequence was updated recently in its first major update since 1996. The new version, called "S288C 2010," was determined from a single yeast colony using modern sequencing technologies and serves as the anchor for further innovations in yeast genomic science.

Unisexual and heterosexual meiotic reproduction generate aneuploidy and phenotypic diversity de novo in the yeast Cryptococcus neoformans.

Aneuploidy is known to be deleterious and underlies several common human diseases, including cancer and genetic disorders such as trisomy 21 in Down's syndrome. In contrast, aneuploidy can also be advantageous and in fungi confers antifungal drug resistance and enables rapid adaptive evolution. We report here that sexual reproduction generates phenotypic and genotypic diversity in the human pathogenic yeast Cryptococcus neoformans, which is globally distributed and commonly infects individuals with compromised immunity, such as HIV/AIDS patients, causing life-threatening meningoencephalitis. C. neoformans has a defined a-α opposite sexual cycle; however, >99% of isolates are of the α mating type. Interestingly, α cells can undergo α-α unisexual reproduction, even involving genotypically identical cells. A central question is: Why would cells mate with themselves given that sex is costly and typically serves to admix preexisting genetic diversity from genetically divergent parents? In this study, we demonstrate that α-α unisexual reproduction frequently generates phenotypic diversity, and the majority of these variant progeny are aneuploid. Aneuploidy is responsible for the observed phenotypic changes, as chromosome loss restoring euploidy results in a wild-type phenotype. Other genetic changes, including diploidization, chromosome length polymorphisms, SNPs, and indels, were also generated. Phenotypic/genotypic changes were not observed following asexual mitotic reproduction. Aneuploidy was also detected in progeny from a-α opposite-sex congenic mating; thus, both homothallic and heterothallic sexual reproduction can generate phenotypic diversity de novo. Our study suggests that the ability to undergo unisexual reproduction may be an evolutionary strategy for eukaryotic microbial pathogens, enabling de novo genotypic and phenotypic plasticity and facilitating rapid adaptation to novel environments.

Reconstructing genome evolution in historic samples of the Irish potato famine pathogen.

Responsible for the Irish potato famine of 1845-49, the oomycete pathogen Phytophthora infestans caused persistent, devastating outbreaks of potato late blight across Europe in the 19th century. Despite continued interest in the history and spread of the pathogen, the genome of the famine-era strain remains entirely unknown. Here we characterize temporal genomic changes in introduced P. infestans. We shotgun sequence five 19th-century European strains from archival herbarium samples--including the oldest known European specimen, collected in 1845 from the first reported source of introduction. We then compare their genomes to those of extant isolates. We report multiple distinct genotypes in historical Europe and a suite of infection-related genes different from modern strains. At virulence-related loci, several now-ubiquitous genotypes were absent from the historical gene pool. At least one of these genotypes encodes a virulent phenotype in modern strains, which helps explain the 20th century's episodic replacements of European P. infestans lineages.

Genomes of Ashbya fungi isolated from insects reveal four mating-type loci, numerous translocations, lack of transposons, and distinct gene duplications.

The filamentous fungus Ashbya gossypii is a cotton pathogen transmitted by insects. It is readily grown and manipulated in the laboratory and is commercially exploited as a natural overproducer of vitamin B2. Our previous genome analysis of A. gossypii isolate ATCC10895, collected in Trinidad nearly 100 years ago, revealed extensive synteny with the Saccharomyces cerevisiae genome, leading us to use it as a model organism to understand the evolution of filamentous growth. To further develop Ashbya as a model system, we have investigated the ecological niche of A. gossypii and isolated additional strains and a sibling species, both useful in comparative analysis. We isolated fungi morphologically similar to A. gossypii from different plant-feeding insects of the suborder Heteroptera, generated a phylogenetic tree based on rDNA-ITS sequences, and performed high coverage short read sequencing with one A. gossypii isolate from Florida, a new species, Ashbya aceri, isolated in North Carolina, and a genetically marked derivative of ATCC10895 intensively used for functional studies. In contrast to S. cerevisiae, all strains carry four not three mating type loci, adding a new puzzle in the evolution of Ashbya species. Another surprise was the genome identity of 99.9% between the Florida strain and ATCC10895, isolated in Trinidad. The A. aceri and A. gossypii genomes show conserved gene orders rearranged by eight translocations, 90% overall sequence identity, and fewer tandem duplications in the A. aceri genome. Both species lack transposable elements. Finally, our work identifies plant-feeding insects of the suborder Heteroptera as the most likely natural reservoir of Ashbya, and that infection of cotton and other plants may be incidental to the growth of the fungus in its insect host.

High-throughput genome sequencing of lichenizing fungi to assess gene loss in the ammonium transporter/ammonia permease gene family.

BACKGROUND: Horizontal gene transfer has shaped the evolution of the ammonium transporter/ammonia permease gene family. Horizontal transfers of ammonium transporter/ammonia permease genes into the fungi include one transfer from archaea to the filamentous ascomycetes associated with the adaptive radiation of the leotiomyceta. The horizontally transferred gene has subsequently been lost in most of the group but has been selectively retained in lichenizing fungi. However, some groups of lichens appear to have secondarily lost the archaeal ammonium transporter. Definitive assessment of gene loss can only be made via whole genome sequencing. RESULTS: Ammonium transporter/ammonia permease gene sequences were recovered from the assembled genomes of eight lichenizing fungi in key clades including the Caliciales, the Peltigerales, the Ostropomycetidae, the Acarosporomycetidae, the Verrucariales, the Arthoniomycetidae and the Lichinales. The genes recovered were included in a refined phylogenetic analysis. The hypothesis that lichens symbiotic with a nitrogen-fixing cyanobacterium as a primary photobiont or lichens living in high nitrogen environments lose the plant-like ammonium transporters was upheld, but did not account for additional losses of ammonium transporters/ammonia permeases in the lichens from the Acarosporomycetidae, Chaetotheriomycetes and Arthoniomycetes. In addition, the four ammonium transporter/ammonia permease genes from Cladonia grayi were shown to be functional by expressing the lichen genes in a strain of Saccharomyces cerevisiae in which all three native ammonium transporters were deleted, and assaying for growth on limiting ammonia as a sole nitrogen source. CONCLUSIONS: Given sufficient coverage, next-generation sequencing technology can definitively address the loss of a gene in a genome when using environmental DNA isolated from lichen thalli collected from their natural habitats. Lichen-forming fungi have been losing ammonium transporters/ammonia permease genes at a slower rate than the most closely related non-lichenized lineages. These horizontally transferred genes in the Cladonia grayi genome encode functional ammonium transporters/ammonia permeases.

Discovery of a modified tetrapolar sexual cycle in Cryptococcus amylolentus and the evolution of MAT in the Cryptococcus species complex.

Sexual reproduction in fungi is governed by a specialized genomic region called the mating-type locus (MAT). The human fungal pathogenic and basidiomycetous yeast Cryptococcus neoformans has evolved a bipolar mating system (a, α) in which the MAT locus is unusually large (>100 kb) and encodes >20 genes including homeodomain (HD) and pheromone/receptor (P/R) genes. To understand how this unique bipolar mating system evolved, we investigated MAT in the closely related species Tsuchiyaea wingfieldii and Cryptococcus amylolentus and discovered two physically unlinked loci encoding the HD and P/R genes. Interestingly, the HD (B) locus sex-specific region is restricted (∼2 kb) and encodes two linked and divergently oriented homeodomain genes in contrast to the solo HD genes (SXI1α, SXI2a) of C. neoformans and Cryptococcus gattii. The P/R (A) locus contains the pheromone and pheromone receptor genes but has expanded considerably compared to other outgroup species (Cryptococcus heveanensis) and is linked to many of the genes also found in the MAT locus of the pathogenic Cryptococcus species. Our discovery of a heterothallic sexual cycle for C. amylolentus allowed us to establish the biological roles of the sex-determining regions. Matings between two strains of opposite mating-types (A1B1×A2B2) produced dikaryotic hyphae with fused clamp connections, basidia, and basidiospores. Genotyping progeny using markers linked and unlinked to MAT revealed that meiosis and uniparental mitochondrial inheritance occur during the sexual cycle of C. amylolentus. The sexual cycle is tetrapolar and produces fertile progeny of four mating-types (A1B1, A1B2, A2B1, and A2B2), but a high proportion of progeny are infertile, and fertility is biased towards one parental mating-type (A1B1). Our studies reveal insights into the plasticity and transitions in both mechanisms of sex determination (bipolar versus tetrapolar) and sexual reproduction (outcrossing versus inbreeding) with implications for similar evolutionary transitions and processes in fungi, plants, and animals.

High-resolution genome-wide analysis of irradiated (UV and γ-rays) diploid yeast cells reveals a high frequency of genomic loss of heterozygosity (LOH) events.

In diploid eukaryotes, repair of double-stranded DNA breaks by homologous recombination often leads to loss of heterozygosity (LOH). Most previous studies of mitotic recombination in Saccharomyces cerevisiae have focused on a single chromosome or a single region of one chromosome at which LOH events can be selected. In this study, we used two techniques (single-nucleotide polymorphism microarrays and high-throughput DNA sequencing) to examine genome-wide LOH in a diploid yeast strain at a resolution averaging 1 kb. We examined both selected LOH events on chromosome V and unselected events throughout the genome in untreated cells and in cells treated with either γ-radiation or ultraviolet (UV) radiation. Our analysis shows the following: (1) spontaneous and damage-induced mitotic gene conversion tracts are more than three times larger than meiotic conversion tracts, and conversion tracts associated with crossovers are usually longer and more complex than those unassociated with crossovers; (2) most of the crossovers and conversions reflect the repair of two sister chromatids broken at the same position; and (3) both UV and γ-radiation efficiently induce LOH at doses of radiation that cause no significant loss of viability. Using high-throughput DNA sequencing, we also detected new mutations induced by γ-rays and UV. To our knowledge, our study represents the first high-resolution genome-wide analysis of DNA damage-induced LOH events performed in any eukaryote.

Multiple horizontal gene transfers of ammonium transporters/ammonia permeases from prokaryotes to eukaryotes: toward a new functional and evolutionary classification.

The proteins of the ammonium transporter/methylammonium permease/Rhesus factor family (AMT/MEP/Rh family) are responsible for the movement of ammonia or ammonium ions across the cell membrane. Although it has been established that the Rh proteins are distantly related to the other members of the family, the evolutionary history of the AMT/MEP/Rh family remains unclear. Here, we use phylogenetic analysis to infer the evolutionary history of this family of proteins across 191 genomes representing all main lineages of life and to provide a new classification of the proteins in this family. Our phylogenetic analysis suggests that what has heretofore been conceived of as a protein family with two clades (AMT/MEP and Rh) is instead a protein family with three clades (AMT, MEP, and Rh). We show that the AMT/MEP/Rh family illustrates two contrasting modes of gene transmission: The AMT family as defined here exhibits vertical gene transfer (i.e., standard parent-to-offspring inheritance), whereas the MEP family as defined here is characterized by several ancient independent horizontal gene transfers (HGTs). These ancient HGT events include a gene replacement during the early evolution of the fungi, which could be a defining trait for the kingdom Fungi, a gene gain from hyperthermophilic chemoautolithotrophic prokaryotes during the early evolution of land plants (Embryophyta), and an independent gain of this same gene in the filamentous ascomycetes (Pezizomycotina) that was subsequently lost in most lineages but retained in even distantly related lichenized fungi. This recircumscription of the ammonium transporters/ammonia permeases family into MEP and AMT families informs the debate on the mechanism of transport in these proteins and on the nature of the transported molecule because published crystal structures of proteins from the MEP and Rh clades may not be representative of the AMT clade. The clades as depicted in this phylogenetic study appear to correspond to functionally different groups, with AMTs and ammonia permeases forming two distinct and possibly monophyletic groups.

A diverse population of Cryptococcus gattii molecular type VGIII in southern Californian HIV/AIDS patients.

Cryptococcus gattii infections in southern California have been reported in patients with HIV/AIDS. In this study, we examined the molecular epidemiology, population structure, and virulence attributes of isolates collected from HIV/AIDS patients in Los Angeles County, California. We show that these isolates consist almost exclusively of VGIII molecular type, in contrast to the VGII molecular type isolates causing the North American Pacific Northwest outbreak. The global VGIII population structure can be divided into two molecular groups, VGIIIa and VGIIIb. Isolates from the Californian patients are virulent in murine and macrophage models of infection, with VGIIIa significantly more virulent than VGIIIb. Several VGIII isolates are highly fertile and produce abundant sexual spores that may serve as infectious propagules. The a and α VGIII MAT locus alleles are largely syntenic with limited rearrangements compared to the known VGI (a/α) and VGII (α) MAT loci, but each has unique characteristics including a distinct deletion flanking the 5' VGIII MATa alleles and the α allele is more heterogeneous than the a allele. Our studies indicate that C. gattii VGIII is endemic in southern California, with other isolates originating from the neighboring regions of Mexico, and in rarer cases from Oregon and Washington state. Given that >1,000,000 cases of cryptococcal infection and >620,000 attributable mortalities occur annually in the context of the global AIDS pandemic, our findings suggest a significant burden of C. gattii may be unrecognized, with potential prognostic and therapeutic implications. These results signify the need to classify pathogenic Cryptococcus cases and highlight possible host differences among the C. gattii molecular types influencing infection of immunocompetent (VGI/VGII) vs. immunocompromised (VGIII/VGIV) hosts.

A whole-genome analysis of premature termination codons.

We sequenced the genomes of ten unrelated individuals and identified heterozygous stop codon-gain variants in protein-coding genes: we then sequenced their transcriptomes and assessed the expression levels of the stop codon-gain alleles. An ANOVA showed statistically significant differences between their expression levels (p=4×10(-16)). This difference was almost entirely accounted for by whether the stop codon-gain variant had a second, non-protein-truncating function in or near an alternate transcript: stop codon-gains without alternate functions were generally not found in the cDNA (p=3×10(-5)). Additionally, stop codon-gain variants in two intronless genes were not expressed, an unexpected outcome given previous studies. In this study, stop codon-gain variants were either well expressed in all individuals or were never expressed. Our finding that stop codon-gain variants were generally expressed only when they had an alternate function suggests that most naturally occurring stop codon-gain variants in protein-coding genes are either not transcribed or have their transcripts destroyed.

Oxidative stress survival in a clinical Saccharomyces cerevisiae isolate is influenced by a major quantitative trait nucleotide.

One of the major challenges in characterizing eukaryotic genetic diversity is the mapping of phenotypes that are the cumulative effect of multiple alleles. We have investigated tolerance of oxidative stress in the yeast Saccharomyces cerevisiae, a trait showing phenotypic variation in the population. Initial crosses identified that this is a quantitative trait. Microorganisms experience oxidative stress in many environments, including during infection of higher eukaryotes. Natural variation in oxidative stress tolerance is an important aspect of response to oxidative stress exerted by the human immune system and an important trait in microbial pathogens. A clinical isolate of the usually benign yeast S. cerevisiae was found to survive oxidative stress significantly better than the laboratory strain. We investigated the genetic basis of increased peroxide survival by crossing those strains, phenotyping 1500 segregants, and genotyping of high-survival segregants by hybridization of bulk and single segregant DNA to microarrays. This effort has led to the identification of an allele of the transcription factor Rds2 as contributing to stress response. Rds2 has not previously been associated with the survival of oxidative stress. The identification of its role in the oxidative stress response here is an example of a specific trait that appears to be beneficial to Saccharomyces cerevisiae when growing as a pathogen. Understanding the role of this fungal-specific transcription factor in pathogenicity will be important in deciphering how fungi infect and colonize the human host and could eventually lead to a novel drug target.

Insights into evolution of multicellular fungi from the assembled chromosomes of the mushroom Coprinopsis cinerea (Coprinus cinereus).

The mushroom Coprinopsis cinerea is a classic experimental model for multicellular development in fungi because it grows on defined media, completes its life cycle in 2 weeks, produces some 10(8) synchronized meiocytes, and can be manipulated at all stages in development by mutation and transformation. The 37-megabase genome of C. cinerea was sequenced and assembled into 13 chromosomes. Meiotic recombination rates vary greatly along the chromosomes, and retrotransposons are absent in large regions of the genome with low levels of meiotic recombination. Single-copy genes with identifiable orthologs in other basidiomycetes are predominant in low-recombination regions of the chromosome. In contrast, paralogous multicopy genes are found in the highly recombining regions, including a large family of protein kinases (FunK1) unique to multicellular fungi. Analyses of P450 and hydrophobin gene families confirmed that local gene duplications drive the expansions of paralogous copies and the expansions occur in independent lineages of Agaricomycotina fungi. Gene-expression patterns from microarrays were used to dissect the transcriptional program of dikaryon formation (mating). Several members of the FunK1 kinase family are differentially regulated during sexual morphogenesis, and coordinate regulation of adjacent duplications is rare. The genomes of C. cinerea and Laccaria bicolor, a symbiotic basidiomycete, share extensive regions of synteny. The largest syntenic blocks occur in regions with low meiotic recombination rates, no transposable elements, and tight gene spacing, where orthologous single-copy genes are overrepresented. The chromosome assembly of C. cinerea is an essential resource in understanding the evolution of multicellularity in the fungi.