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.

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.

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.

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.

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.

Morphogenesis in skin is governed by discrete sets of differentially expressed microRNAs.

During embryogenesis, multipotent progenitors within the single-layered surface epithelium differentiate to form the epidermis and its appendages. Here, we show that microRNAs (miRNAs) have an essential role in orchestrating these events. We cloned more than 100 miRNAs from skin and show that epidermis and hair follicles differentially express discrete miRNA families. To explore the functional significance of this finding, we conditionally targeted Dicer1 gene ablation in embryonic skin progenitors. Within the first week after loss of miRNA expression, cell fate specification and differentiation were not markedly impaired, and in the interfollicular epidermis, apoptosis was not markedly increased. Notably, however, developing hair germs evaginate rather than invaginate, thereby perturbing the epidermal organization. Here we characterize miRNAs in skin, the existence of which was hitherto unappreciated, and demonstrate their differential expression and importance in the morphogenesis of epithelial tissues within this vital organ.

Distinct evolutionary trajectories following loss of RNA interference in Cryptococcus neoformans.

While increased mutation rates typically have negative consequences in multicellular organisms, hypermutation can be advantageous for microbes adapting to the environment. Previously, we identified two hypermutator Cryptococcus neoformans clinical isolates that rapidly develop drug resistance due to transposition of a retrotransposon, Cnl1. Cnl1-mediated hypermutation is caused by a nonsense mutation in the gene encoding a novel RNAi component, Znf3, combined with a tremendous transposon burden. To elucidate adaptative mechanisms following RNAi loss, two bioinformatic pipelines were developed to identify RNAi loss-of-function mutations in a collection of 387 sequenced C. neoformans isolates. Remarkably, several RNAi-loss isolates were identified that are not hypermutators and have not accumulated transposons. To test if these RNAi loss-of-function mutations can cause hypermutation, the mutations were introduced into a non-hypermutator strain with a high transposon burden, which resulted in a hypermutator phenotype. To further investigate if RNAi-loss isolates can become hypermutators, in vitro passaging was performed. Although no hypermutators were found in two C. neoformans RNAi-loss strains after short-term passage, hypermutation was observed in a passaged C. deneoformans strain with increased transposon burden. Additionally, when an RNAi-loss isolate was crossed with an isolate containing a high Cnl1 burden, F1 hypermutator progeny were identified with distinct mutational spectra. In addition to Cnl1 transpositions, insertions of a novel gigantic DNA transposon KDZ1 (~11 kb), contributed to hypermutation in the progeny. Our results suggest that RNAi loss is relatively common (7/387, ~1.8%) and enables distinct evolutionary trajectories: hypermutation following transposon accumulation or survival without hypermutation.

Intrinsically disordered sequences can tune fungal growth and the cell cycle for specific temperatures.

Temperature can impact every reaction essential to a cell. For organisms that cannot regulate their own temperature, adapting to temperatures that fluctuate unpredictably and on variable timescales is a major challenge. Extremes in the magnitude and frequency of temperature changes are increasing across the planet, raising questions as to how the biosphere will respond. To examine mechanisms of adaptation to temperature, we collected wild isolates from different climates of the fungus Ashbya gossypii, which has a compact genome of only ∼4,600 genes. We found control of the nuclear division cycle and polarized morphogenesis, both critical processes for fungal growth, were temperature sensitive and varied among the isolates. The phenotypes were associated with naturally varying sequences within the glutamine-rich region (QRR) IDR of an RNA-binding protein called Whi3. This protein regulates both nuclear division and polarized growth via its ability to form biomolecular condensates. In cells and in cell-free reconstitution assays, we found that temperature tunes the properties of Whi3-based condensates. Exchanging Whi3 sequences between isolates was sufficient to rescue temperature-sensitive phenotypes, and specifically, a heptad repeat sequence within the QRR confers temperature-sensitive behavior. Together, these data demonstrate that sequence variation in the size and composition of an IDR can promote cell adaptation to growth at specific temperature ranges. These data demonstrate the power of IDRs as tuning knobs for rapid adaptation to environmental fluctuations.

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.

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.

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.

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.

Genome structure of a Saccharomyces cerevisiae strain widely used in bioethanol production.

Bioethanol is a biofuel produced mainly from the fermentation of carbohydrates derived from agricultural feedstocks by the yeast Saccharomyces cerevisiae. One of the most widely adopted strains is PE-2, a heterothallic diploid naturally adapted to the sugar cane fermentation process used in Brazil. Here we report the molecular genetic analysis of a PE-2 derived diploid (JAY270), and the complete genome sequence of a haploid derivative (JAY291). The JAY270 genome is highly heterozygous (approximately 2 SNPs/kb) and has several structural polymorphisms between homologous chromosomes. These chromosomal rearrangements are confined to the peripheral regions of the chromosomes, with breakpoints within repetitive DNA sequences. Despite its complex karyotype, this diploid, when sporulated, had a high frequency of viable spores. Hybrid diploids formed by outcrossing with the laboratory strain S288c also displayed good spore viability. Thus, the rearrangements that exist near the ends of chromosomes do not impair meiosis, as they do not span regions that contain essential genes. This observation is consistent with a model in which the peripheral regions of chromosomes represent plastic domains of the genome that are free to recombine ectopically and experiment with alternative structures. We also explored features of the JAY270 and JAY291 genomes that help explain their high adaptation to industrial environments, exhibiting desirable phenotypes such as high ethanol and cell mass production and high temperature and oxidative stress tolerance. The genomic manipulation of such strains could enable the creation of a new generation of industrial organisms, ideally suited for use as delivery vehicles for future bioenergy technologies.

Non-coding RNA prediction and verification in Saccharomyces cerevisiae.

Non-coding RNA (ncRNA) play an important and varied role in cellular function. A significant amount of research has been devoted to computational prediction of these genes from genomic sequence, but the ability to do so has remained elusive due to a lack of apparent genomic features. In this work, thermodynamic stability of ncRNA structural elements, as summarized in a Z-score, is used to predict ncRNA in the yeast Saccharomyces cerevisiae. This analysis was coupled with comparative genomics to search for ncRNA genes on chromosome six of S. cerevisiae and S. bayanus. Sets of positive and negative control genes were evaluated to determine the efficacy of thermodynamic stability for discriminating ncRNA from background sequence. The effect of window sizes and step sizes on the sensitivity of ncRNA identification was also explored. Non-coding RNA gene candidates, common to both S. cerevisiae and S. bayanus, were verified using northern blot analysis, rapid amplification of cDNA ends (RACE), and publicly available cDNA library data. Four ncRNA transcripts are well supported by experimental data (RUF10, RUF11, RUF12, RUF13), while one additional putative ncRNA transcript is well supported but the data are not entirely conclusive. Six candidates appear to be structural elements in 5' or 3' untranslated regions of annotated protein-coding genes. This work shows that thermodynamic stability, coupled with comparative genomics, can be used to predict ncRNA with significant structural elements.

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.

Microsporidia evolved from ancestral sexual fungi.

Microsporidia are obligate, intracellular eukaryotic pathogens that infect animal cells, including humans [1]. Previous studies suggested microsporidia share a common ancestor with fungi [2-7]. However, the exact nature of this phylogenetic relationship is unclear because of unusual features of microsporidial genomes, which are compact with fewer and highly divergent genes [8]. As a consequence, it is unclear whether microsporidia evolved from a specific fungal lineage, or whether microsporidia are a sister group to all fungi. Here, we present evidence addressing this controversial question that is independent of sequence-based phylogenetic reconstruction, but rather based on genome structure. In the zygomycete basal fungal lineage, the sex locus is a syntenic gene cluster governing sexual reproduction in which a high mobility group (HMG) transcription-factor gene is flanked by triose-phosphate transporter (TPT) and RNA helicase genes [9]. Strikingly, microsporidian genomes harbor a sex-related locus with the same genes in the same order. Genome-wide synteny analysis reveals multiple other loci conserved between microsporidia and zygomycetes to the exclusion of all other fungal lineages with sequenced genomes. These findings support the hypothesis that microsporidia are true fungi that descended from a zygomycete ancestor and suggest microsporidia may have an extant sexual cycle.

Same-sex mating and the origin of the Vancouver Island Cryptococcus gattii outbreak.

Genealogy can illuminate the evolutionary path of important human pathogens. In some microbes, strict clonal reproduction predominates, as with the worldwide dissemination of Mycobacterium leprae, the cause of leprosy. In other pathogens, sexual reproduction yields clones with novel attributes, for example, enabling the efficient, oral transmission of the parasite Toxoplasma gondii. However, the roles of clonal or sexual propagation in the origins of many other microbial pathogen outbreaks remain unknown, like the recent fungal meningoencephalitis outbreak on Vancouver Island, Canada, caused by Cryptococcus gattii. Here we show that the C. gattii outbreak isolates comprise two distinct genotypes. The majority of isolates are hypervirulent and have an identical genotype that is unique to the Pacific Northwest. A minority of the isolates are significantly less virulent and share an identical genotype with fertile isolates from an Australian recombining population. Genotypic analysis reveals evidence of sexual reproduction, in which the majority genotype is the predicted offspring. However, instead of the classic a-alpha sexual cycle, the majority outbreak clone appears to have descended from two alpha mating-type parents. Analysis of nuclear content revealed a diploid environmental isolate homozygous for the major genotype, an intermediate produced during same-sex mating. These studies demonstrate how cryptic same-sex reproduction can enable expansion of a human pathogen to a new geographical niche and contribute to the ongoing production of infectious spores. This has implications for the emergence of other microbial pathogens and inbreeding in host range expansion in the fungal and other kingdoms.

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.