HGT in the human and skin commensal Malassezia: A bacterially derived flavohemoglobin is required for NO resistance and host interaction.

The skin of humans and animals is colonized by commensal and pathogenic fungi and bacteria that share this ecological niche and have established microbial interactions. Malassezia are the most abundant fungal skin inhabitant of warm-blooded animals and have been implicated in skin diseases and systemic disorders, including Crohn's disease and pancreatic cancer. Flavohemoglobin is a key enzyme involved in microbial nitrosative stress resistance and nitric oxide degradation. Comparative genomics and phylogenetic analyses within the Malassezia genus revealed that flavohemoglobin-encoding genes were acquired through independent horizontal gene transfer events from different donor bacteria that are part of the mammalian microbiome. Through targeted gene deletion and functional complementation in Malassezia sympodialis, we demonstrated that bacterially derived flavohemoglobins are cytoplasmic proteins required for nitric oxide detoxification and nitrosative stress resistance under aerobic conditions. RNA-sequencing analysis revealed that endogenous accumulation of nitric oxide resulted in up-regulation of genes involved in stress response and down-regulation of the MalaS7 allergen-encoding genes. Solution of the high-resolution X-ray crystal structure of Malassezia flavohemoglobin revealed features conserved with both bacterial and fungal flavohemoglobins. In vivo pathogenesis is independent of Malassezia flavohemoglobin. Lastly, we identified an additional 30 genus- and species-specific horizontal gene transfer candidates that might have contributed to the evolution of this genus as the most common inhabitants of animal skin.

The lichen symbiosis re-viewed through the genomes of Cladonia grayi and its algal partner Asterochloris glomerata.

BACKGROUND: Lichens, encompassing 20,000 known species, are symbioses between specialized fungi (mycobionts), mostly ascomycetes, and unicellular green algae or cyanobacteria (photobionts). Here we describe the first parallel genomic analysis of the mycobiont Cladonia grayi and of its green algal photobiont Asterochloris glomerata. We focus on genes/predicted proteins of potential symbiotic significance, sought by surveying proteins differentially activated during early stages of mycobiont and photobiont interaction in coculture, expanded or contracted protein families, and proteins with differential rates of evolution. RESULTS: A) In coculture, the fungus upregulated small secreted proteins, membrane transport proteins, signal transduction components, extracellular hydrolases and, notably, a ribitol transporter and an ammonium transporter, and the alga activated DNA metabolism, signal transduction, and expression of flagellar components. B) Expanded fungal protein families include heterokaryon incompatibility proteins, polyketide synthases, and a unique set of G-protein α subunit paralogs. Expanded algal protein families include carbohydrate active enzymes and a specific subclass of cytoplasmic carbonic anhydrases. The alga also appears to have acquired by horizontal gene transfer from prokaryotes novel archaeal ATPases and Desiccation-Related Proteins. Expanded in both symbionts are signal transduction components, ankyrin domain proteins and transcription factors involved in chromatin remodeling and stress responses. The fungal transportome is contracted, as are algal nitrate assimilation genes. C) In the mycobiont, slow-evolving proteins were enriched for components involved in protein translation, translocation and sorting. CONCLUSIONS: The surveyed genes affect stress resistance, signaling, genome reprogramming, nutritional and structural interactions. The alga carries many genes likely transferred horizontally through viruses, yet we found no evidence of inter-symbiont gene transfer. The presence in the photobiont of meiosis-specific genes supports the notion that sexual reproduction occurs in Asterochloris while they are free-living, a phenomenon with implications for the adaptability of lichens and the persistent autonomy of the symbionts. The diversity of the genes affecting the symbiosis suggests that lichens evolved by accretion of many scattered regulatory and structural changes rather than through introduction of a few key innovations. This predicts that paths to lichenization were variable in different phyla, which is consistent with the emerging consensus that ascolichens could have had a few independent origins.

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.

The deacetylase Sir2 from the yeast Clavispora lusitaniae lacks the evolutionarily conserved capacity to generate subtelomeric heterochromatin.

Deacetylases of the Sir2 or sirtuin family are thought to regulate life cycle progression and life span in response to nutrient availability. This family has undergone successive rounds of duplication and diversification, enabling the enzymes to perform a wide variety of biological functions. Two evolutionarily conserved functions of yeast Sir2 proteins are the generation of repressive chromatin in subtelomeric domains and the suppression of unbalanced recombination within the tandem rDNA array. Here, we describe the function of the Sir2 ortholog ClHst1 in the yeast Clavispora lusitaniae, an occasional opportunistic pathogen. ClHst1 was localized to the non-transcribed spacer regions of the rDNA repeats and deacetylated histones at these loci, indicating that, like other Sir2 proteins, ClHst1 modulates chromatin structure at the rDNA repeats. However, we found no evidence that ClHst1 associates with subtelomeric regions or impacts gene expression directly. This surprising observation highlights the plasticity of sirtuin function. Related yeast species, including Candida albicans, possess an additional Sir2 family member. Thus, it is likely that the ancestral Candida SIR2/HST1 gene was duplicated and subfunctionalized, such that HST1 retained the capacity to regulate rDNA whereas SIR2 had other functions, perhaps including the generation of subtelomeric chromatin. After subsequent species diversification, the SIR2 paralog was apparently lost in the C. lusitaniae lineage. Thus, C. lusitaniae presents an opportunity to discover how subtelomeric chromatin can be reconfigured.

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.

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.

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.

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.

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.