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.

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.

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 environmentally-acquired human pathogen C. gattii is responsible for ongoing and expanding outbreaks in the Western United States and Canada. C. gattii comprises four distinct molecular types: VGI, VGII, VGIII, and VGIV. Molecular types VGI, VGII, and VGIII have been isolated from patients and animals throughout the Western US. The Pacific Northwest and Canadian outbreak is primarily caused by C. gattii VGII. VGIII is responsible for ongoing infections in HIV/AIDS patients in Southern California. However, only two environmental C. gattii isolates have ever been identified from the Californian environment: CBS7750 (VGII) and WM161 (VGIII). We sought to collect environmental samples from areas that had confirmed reports of clinical or veterinary infections. Here we report the isolation of C. gattii VGI and VGIII from environmental soil and tree samples. C. gattii isolates were obtained from three novel tree species: Canary Island pine, American sweetgum, and a Pohutukawa tree. Genetic analysis provides robust evidence that these environmental isolates are the source of human infections.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Cryptococcus neoformans var. grubii is a major human pathogen responsible for deadly meningoencephalitis in immunocompromised patients. Here, we report the sequencing and annotation of its genome. Evidence for extensive intron splicing, antisense transcription, non-coding RNAs, and alternative polyadenylation indicates the potential for highly intricate regulation of gene expression in this opportunistic pathogen. In addition, detailed molecular, genetic, and genomic studies were performed to characterize structural features of the genome, including centromeres and origins of replication. Finally, the phenotypic and genome re-sequencing analysis of a collection of isolates of the reference H99 strain resulting from laboratory passage revealed that microevolutionary processes during in vitro culturing of pathogenic fungi can impact virulence.

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.

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.

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.

Coupled-channel calculation of nonelastic cross sections using a density-functional structure model

A microscopic calculation of reaction cross sections for nucleon-nucleus scattering was performed by coupling the elastic channel to all particle-hole excitations in the target and one-nucleon pickup channels. The particle-hole states may be regarded as doorway states through which the flux flows to more complicated configurations, and subsequently to long-lived compound nucleus resonances. Target excitations for (40,48)Ca, 58Ni, 90Zr, and 144Sm were described in a random-phase framework using a Skyrme functional. Reaction cross sections obtained agreed very well with experimental data and predictions of a fitted optical potential. Couplings between inelastic states were found to be negligible, while the pickup channels contribute significantly. For the first time observed absorptions are completely accounted for by explicit channel coupling, for incident energies between 10 and 40 MeV.

Evolution of the sex-related locus and genomic features shared in microsporidia and fungi

Background

Microsporidia are obligate intracellular, eukaryotic pathogens that infect a wide range of animals from nematodes to humans, and in some cases, protists. The preponderance of evidence as to the origin of the microsporidia reveals a close relationship with the fungi, either within the kingdom or as a sister group to it. Recent phylogenetic studies and gene order analysis suggest that microsporidia share a particularly close evolutionary relationship with the zygomycetes.

Methodology/Principal Findings

Here we expanded this analysis and also examined a putative sex-locus for variability between microsporidian populations. Whole genome inspection reveals a unique syntenic gene pair (RPS9-RPL21) present in the vast majority of fungi and the microsporidians but not in other eukaryotic lineages. Two other unique gene fusions (glutamyl-prolyl tRNA synthetase and ubiquitin-ribosomal subunit S30) that are present in metazoans, choanoflagellates, and filasterean opisthokonts are unfused in the fungi and microsporidians. One locus previously found to be conserved in many microsporidian genomes is similar to the sex locus of zygomycetes in gene order and architecture. Both sex-related and sex loci harbor TPT, HMG, and RNA helicase genes forming a syntenic gene cluster. We sequenced and analyzed the sex-related locus in 11 different Encephalitozoon cuniculi isolates and the sibling species E. intestinalis (3 isolates) and E. hellem (1 isolate). There was no evidence for an idiomorphic sex-related locus in this Encephalitozoon species sample. According to sequence-based phylogenetic analyses, the TPT and RNA helicase genes flanking the HMG genes are paralogous rather than orthologous between zygomycetes and microsporidians.

Conclusion/Significance

The unique genomic hallmarks between microsporidia and fungi are independent of sequence based phylogenetic comparisons and further contribute to define the borders of the fungal kingdom and support the classification of microsporidia as unusual derived fungi. And the sex/sex-related loci appear to have been subject to frequent gene conversion and translocations in microsporidia and zygomycetes.

Morphological and genomic characterization of Filobasidiella depauperata: a homothallic sibling species of the pathogenic cryptococcus species complex

The fungal species Cryptococcus neoformans and Cryptococcus gattii cause respiratory and neurological disease in animals and humans following inhalation of basidiospores or desiccated yeast cells from the environment. Sexual reproduction in C. neoformans and C. gattii is controlled by a bipolar system in which a single mating type locus (MAT) specifies compatibility. These two species are dimorphic, growing as yeast in the asexual stage, and producing hyphae, basidia, and basidiospores during the sexual stage. In contrast, Filobasidiella depauperata, one of the closest related species, grows exclusively as hyphae and it is found in association with decaying insects. Examination of two available strains of F. depauperata showed that the life cycle of this fungal species shares features associated with the unisexual or same-sex mating cycle in C. neoformans. Therefore, F. depauperata may represent a homothallic and possibly an obligately sexual fungal species. RAPD genotyping of 39 randomly isolated progeny from isolate CBS7855 revealed a new genotype pattern in one of the isolated basidiospores progeny, therefore suggesting that the homothallic cycle in F. depauperata could lead to the emergence of new genotypes. Phylogenetic analyses of genes linked to MAT in C. neoformans indicated that two of these genes in F. depauperata, MYO2 and STE20, appear to form a monophyletic clade with the MATa alleles of C. neoformans and C. gattii, and thus these genes may have been recruited to the MAT locus before F. depauperata diverged. Furthermore, the ancestral MATa locus may have undergone accelerated evolution prior to the divergence of the pathogenic Cryptococcus species since several of the genes linked to the MATa locus appear to have a higher number of changes and substitutions than their MATalpha counterparts. Synteny analyses between C. neoformans and F. depauperata showed that genomic regions on other chromosomes displayed conserved gene order. In contrast, the genes linked to the MAT locus of C. neoformans showed a higher number of chromosomal translocations in the genome of F. depauperata. We therefore propose that chromosomal rearrangements appear to be a major force driving speciation and sexual divergence in these closely related pathogenic and saprobic species.

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.

Saccharomyces cerevisiae: population divergence and resistance to oxidative stress in clinical, domesticated and wild isolates

BACKGROUND

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

METHODOLOGY/PRINCIPAL FINDINGS

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

CONCLUSIONS/SIGNIFICANCE

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

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.

Recent advances in DNA sequence technology have made it possible to sequence entire genomes. Once a genome is sequenced, it becomes necessary to identify the set of genes and other functional elements within the genome. This is particularly challenging as much of the genomic sequence does not appear to perform any function and is loosely referred to as “junk.” Identifying functional elements among the “junk” is difficult. Experimental methods have been developed for this purpose but they are time-consuming, expensive, and often provide an incomplete picture. Thus, it is important to develop the ability to identify these functional elements using computational methods. Protein-coding genes are relatively easy to identify computationally, but other categories of functional elements present a significantly greater challenge. In this work, we used a computational approach to identify genes that do not encode for a protein but rather function as an RNA molecule. We then used experimental methods to verify our predictions and thereby validate the computational method.

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.

Comparative genomic analysis of fungal genomes reveals intron-rich ancestors

Analysis of intron gain and loss in fungal genomes provides support for an intron-rich fungus-animal ancestor.

Background

Eukaryotic protein-coding genes are interrupted by spliceosomal introns, which are removed from transcripts before protein translation. Many facets of spliceosomal intron evolution, including age, mechanisms of origins, the role of natural selection, and the causes of the vast differences in intron number between eukaryotic species, remain debated. Genome sequencing and comparative analysis has made possible whole genome analysis of intron evolution to address these questions.

Results

We analyzed intron positions in 1,161 sets of orthologous genes across 25 eukaryotic species. We find strong support for an intron-rich fungus-animal ancestor, with more than four introns per kilobase, comparable to the highest known modern intron densities. Indeed, the fungus-animal ancestor is estimated to have had more introns than any of the extant fungi in this study. Thus, subsequent fungal evolution has been characterized by widespread and recurrent intron loss occurring in all fungal clades. These results reconcile three previously proposed methods for estimation of ancestral intron number, which previously gave very different estimates of ancestral intron number for eight eukaryotic species, as well as a fourth more recent method. We do not find a clear inverse correspondence between rates of intron loss and gain, contrary to the predictions of selection-based proposals for interspecific differences in intron number.

Conclusion

Our results underscore the high intron density of eukaryotic ancestors and the widespread importance of intron loss through eukaryotic evolution.

Comparative hybridization reveals extensive genome variation in the AIDS-associated pathogen Cryptococcus neoformans

Extensive genome variation in the AIDS-associated pathogen Cryptococcus neoformans is revealed through comparative genome hybridization between strains of different mating type, molecular subtype and ploidy.

Background

Genome variability can have a profound influence on the virulence of pathogenic microbes. The availability of genome sequences for two strains of the AIDS-associated fungal pathogen Cryptococcus neoformans presented an opportunity to use comparative genome hybridization (CGH) to examine genome variability between strains of different mating type, molecular subtype, and ploidy.

Results

Initially, CGH was used to compare the approximately 100 kilobase MATa and MATα mating-type regions in serotype A and D strains to establish the relationship between the Log2 ratios of hybridization signals and sequence identity. Subsequently, we compared the genomes of the environmental isolate NIH433 (MATa) and the clinical isolate NIH12 (MATα) with a tiling array of the genome of the laboratory strain JEC21 derived from these strains. In this case, CGH identified putative recombination sites and the origins of specific segments of the JEC21 genome. Similarly, CGH analysis revealed marked variability in the genomes of strains representing the VNI, VNII, and VNB molecular subtypes of the A serotype, including disomy for chromosome 13 in two strains. Additionally, CGH identified differences in chromosome content between three strains with the hybrid AD serotype and revealed that chromosome 1 from the serotype A genome is preferentially retained in all three strains.

Conclusion

The genomes of serotypes A, D, and AD strains exhibit extensive variation that spans the range from small differences (such as regions of divergence, deletion, or amplification) to the unexpected disomy for chromosome 13 in haploid strains and preferential retention of specific chromosomes in naturally occurring diploids.

The reacquisition of biotin prototrophy in Saccharomyces cerevisiae involved horizontal gene transfer, gene duplication and gene clustering

The synthesis of biotin, a vitamin required for many carboxylation reactions, is a variable trait in Saccharomyces cerevisiae. Many S. cerevisiae strains, including common laboratory strains, contain only a partial biotin synthesis pathway. We here report the identification of the first step necessary for the biotin synthesis pathway in S. cerevisiae. The biotin auxotroph strain S288c was able to grow on media lacking biotin when BIO1 and the known biotin synthesis gene BIO6 were introduced together on a plasmid vector. BIO1 is a paralog of YJR154W, a gene of unknown function and adjacent to BIO6. The nature of BIO1 illuminates the remarkable evolutionary history of the biotin biosynthesis pathway in S. cerevisiae. This pathway appears to have been lost in an ancestor of S. cerevisiae and subsequently rebuilt by a combination of horizontal gene transfer and gene duplication followed by neofunctionalization. Unusually, for S. cerevisiae, most of the genes required for biotin synthesis in S. cerevisiae are grouped in two subtelomeric gene clusters. The BIO1-BIO6 functional cluster is an example of a cluster of genes of "dispensable function," one of the few categories of genes in S. cerevisiae that are positionally clustered.

Computational and experimental identification of novel human imprinted genes

Imprinted genes are essential in embryonic development, and imprinting dysregulation contributes to human disease. We report two new human imprinted genes: KCNK9 is predominantly expressed in the brain, is a known oncogene, and may be involved in bipolar disorder and epilepsy, while DLGAP2 is a candidate bladder cancer tumor suppressor. Both genes lie on chromosome 8, not previously suspected to contain imprinted genes. We identified these genes, along with 154 others, based on the predictions of multiple classification algorithms using DNA sequence characteristics as features. Our findings demonstrate that DNA sequence characteristics, including recombination hot spots, are sufficient to accurately predict the imprinting status of individual genes in the human genome.

Genome sequencing and comparative analysis of Saccharomyces cerevisiae strain YJM789

We sequenced the genome of Saccharomyces cerevisiae strain YJM789, which was derived from a yeast isolated from the lung of an AIDS patient with pneumonia. The strain is used for studies of fungal infections and quantitative genetics because of its extensive phenotypic differences to the laboratory reference strain, including growth at high temperature and deadly virulence in mouse models. Here we show that the approximately 12-Mb genome of YJM789 contains approximately 60,000 SNPs and approximately 6,000 indels with respect to the reference S288c genome, leading to protein polymorphisms with a few known cases of phenotypic changes. Several ORFs are found to be unique to YJM789, some of which might have been acquired through horizontal transfer. Localized regions of high polymorphism density are scattered over the genome, in some cases spanning multiple ORFs and in others concentrated within single genes. The sequence of YJM789 contains clues to pathogenicity and spurs the development of more powerful approaches to dissecting the genetic basis of complex hereditary traits.

Evolution of Budding Yeast Prion-determinant Sequences Across Diverse Fungi

Prions are transmissible self-replicating alternative states of proteins. Four prions ([PSI+], [URE3], [RNQ+] and [NU+]) can be inherited cytoplasmically in Saccharomyces cerevisiae laboratory strains. In the case of [PSI+], there is increasing evidence that prion formation may engender mechanisms to uncover hidden genetic variation. Here, we have analysed the evolution of the prion-determinant (PD) domains across 21 fungi, focusing on compositional biases, repeats and substitution rates. We find evidence for constraint on all four PD domains, but each domain has its own evolutionary dynamics. For [PSI+], the Q/N bias is maintained in fungal clades that diverged one billion years ago, with purifying selection observed within the Saccharomyces species. The degree of Q/N bias is correlated with the degree of local homology to prion-associated repeats, which occur rarely in other proteins (<1% of sequences for the proteomes studied). The evolutionary conservation of Q/N bias in Sup35p is unusual, with only eight other S. cerevisiae proteins showing similar, phylogenetically deep patterns of bias conservation. The [URE3] PD domain is unique to Hemiascomycota; part of the PD domain shows purifying selection, whereas another part engenders bias changes between clades. Also, like for Sup35p, the [RNQ+] and [NU+] PD domains show purifying selection in Saccharomyces species. Additionally, in each proteome, we observe on average several hundred yeast-prion-like domains, with fewest in fission yeast. Our findings on yeast prion evolution provide further support for the functional significance of these molecules.

Recent evolution of the human pathogen Cryptococcus neoformans by intervarietal transfer of a 14-gene fragment

The availability of the whole-genome sequence from the 2 known varieties of the human pathogenic fungus Cryptococcus neoformans provides an opportunity to study the relative contribution of divergence and introgression during the process of speciation in a genetically tractable organism. At the genomic level, these varieties are nearly completely syntenic, share approximately 85-90% nucleotide identity, and are believed to have diverged approximately 18 MYA. Via a comparative genomic approach, we identified a 14-gene region (approximately 40 kb) that is nearly identical between the 2 varieties that resulted from a nonreciprocal transfer event from var. grubii to var. neoformans approximately 2 MYA. The majority of clinical and environmental var. neoformans strains from around the world contain this sequence obtained from var. grubii. This introgression event likely occurred via an incomplete intervarietal sexual cycle, creating a hybrid intermediate where mobile elements common to both lineages mediated the exchange. The subsequent duplication in laboratory strains of a fragment of this same genomic region supports evolutionary theories that instabilities in subtelomeric regions promote adaptive evolution through gene amplification and subsequent adaptation. Along with a more ancient predicted transfer event in C. neoformans and a recently reported example from Saccharomyces cerevisiae, these data indicate that DNA exchange between closely related sympatric varieties or species may be a recurrent theme in the evolution of fungal species. It further suggests that although evolutionary divergence is the primary force driving speciation, rare introgression events also play a potentially important role.

Phylogenomic analysis of non-ribosomal peptide synthetases in the genus Aspergillus

Fungi from the genus Aspergillus are important saprophytes and opportunistic human fungal pathogens that contribute in these and other diverse ways to human well-being. Part of their impact on human well-being stems from the production of small molecular weight secondary metabolites, which may contribute to the ability of these fungi to cause invasive fungal infections and allergic diseases. In this study, we identified one group of enzymes responsible for secondary metabolite production in five Aspergillus species, the non-ribosomal peptide synthetases (NRPS). Hidden Markov models were used to search the genome databases of A. fumigatus, A. flavus, A. terreus, A. nidulans, and A. oryzae for domains conserved in NRPS proteins. A genealogy of adenylation domains was utilized to identify orthologous and unique NRPS among the Aspergillus species examined, as well as gain an understanding of the potential evolution of Aspergillus NRPS. mRNA abundance of the 14 NRPS identified in the A. fumigatus genome was analyzed using real-time reverse transcriptase PCR in different environmental conditions to gain a preliminary understanding of the possible functions of the NRPSs' peptide products. Our results suggest that Aspergillus species contain conserved and unique NRPS genes with a complex evolutionary history. This result suggests that the genus Aspergillus produces a substantial diversity of non-ribosomally synthesized peptides. Further analysis of these genes and their peptide products may identify important roles for secondary metabolites produced by NRPS in Aspergillus physiology, ecology, and fungal pathogenicity.

Evidence of mRNA-mediated intron loss in the human-pathogenic fungus Cryptococcus neoformans

Introns are a defining feature of eukaryotic genomes, though the mechanism of intron gain or loss is not well understood. Reverse transcription of mRNA followed by homologous recombination with the genome has been posited as a mechanism of intron loss, though little direct evidence of recent loss events has been described to support this model. We find supporting evidence for an mRNA-mediated mechanism of loss through comparative genome analyses that revealed a recent loss of 10 adjacent introns in a 22-exon gene in the human-pathogenic fungus Cryptococcus neoformans. We surveyed the gene structures of the entire genomes of Cryptococcus gattii, which diverged from the C. neoformans lineage 37 million years ago (Mya), and C. neoformans var. grubii and var. neoformans, which diverged 18 Mya. Our comparison revealed greater than 99.9% intron conservation, with evidence from 20 genes showing evidence of intron loss, but no convincing evidence of intron gain. Our findings confirm that Cryptococcus introns have been quite stable over recent evolutionary time, with occasional mRNA-mediated intron loss events.

Yeast diversity sampling on the San Juan Islands reveals no evidence for the spread of the Vancouver IslandCryptococcus gattiioutbreak to this locale

Biological diversity has been estimated for various phyla of life, such as insects and mammals, but in the microbe world is has been difficult to determine species richness and abundance. Here we describe a study of species diversity of fungi with a yeast-like colony morphology from the San Juan Islands, a group of islands that lies southeast of Vancouver Island, Canada. Our sampling revealed that the San Juan archipelago biosphere contains a diverse range of such fungi predominantly belonging to the Basidiomycota, particularly of the order Tremellales. One member of this group, Cryptococcus gattii, is the etiological agent of a current and ongoing outbreak of cryptococcosis on nearby Vancouver Island. Our sampling did not, however, reveal this species. While the lack of recovery of C. gattii does not preclude its presence on the San Juan Islands, our results suggest that the Strait of Juan de Fuca may be serving as a geographical barrier to restrict the dispersal of this primary human fungal pathogen into the United States.

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.

Clinical and environmental isolates of Cryptococcus gattii from Australia that retain sexual fecundity

Cryptococcus gattii is a primary pathogenic yeast that causes disease in both animals and humans. It is closely related to Cryptococcus neoformans and diverged from a common ancestor approximately 40 million years ago. While C. gattii has a characterized sexual cycle dependent upon a dimorphic region of the genome known as the MAT locus, mating has rarely been observed in this species. In this study, we identify for the first time clinical (both human and veterinary) and environmental isolates from Australia that retain sexual fecundity. A collection of 120 isolates from a variety of geographic locations was analyzed for molecular type, mating type, and the ability to develop mating structures when cocultured with fertile tester strains. Nine isolates produced dikaryotic filaments with paired nuclei, fused clamp connections, and basidiospores. DNA sequence analysis of three genes (URA5, the MATalpha-specific SXI1alpha gene, and the MATa-specific SXI2a gene) revealed little or no variability in URA5 and SXI2a, respectively. However across the 108 MATalpha strains sequenced, the SXI1alpha gene was found to exist as 11 different alleles. Phylogenetic analysis found most variation to occur in the more fertile genotypes. Although some lineages of Australian C. gattii have retained the ability to mate, the majority of isolates were sterile, suggesting that asexuality is the dominant mode of propagation in these populations.

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.

Identification and characterization of upstream open reading frames (uORF) in the 5' untranslated regions (UTR) of genes in Saccharomyces cerevisiae

We have taken advantage of recently sequenced hemiascomycete fungal genomes to computationally identify additional genes potentially regulated by upstream open reading frames (uORFs). Our approach is based on the observation that the structure, including the uORFs, of the post-transcriptionally uORF regulated Saccharomyces cerevisiae genes GCN4 and CPA1 is conserved in related species. Thirty-eight candidate genes for which uORFs were found in multiple species were identified and tested. We determined by 5' RACE that 15 of these 38 genes are transcribed. Most of these 15 genes have only a single uORF in their 5' UTR, and the length of these uORFs range from 3 to 24 codons. We cloned seven full-length UTR sequences into a luciferase (LUC) reporter system. Luciferase activity and mRNA level were compared between the wild-type UTR construct and a construct where the uORF start codon was mutated. The translational efficiency index (TEI) of each construct was calculated to test the possible regulatory function on translational level. We hypothesize that uORFs in the UTR of RPC11, TPK1, FOL1, WSC3, and MKK1 may have translational regulatory roles while uORFs in the 5' UTR of ECM7 and IMD4 have little effect on translation under the conditions tested.

Contribution of horizontal gene transfer to the evolution of Saccharomyces cerevisiae

The genomes of the hemiascomycetes Saccharomyces cerevisiae and Ashbya gossypii have been completely sequenced, allowing a comparative analysis of these two genomes, which reveals that a small number of genes appear to have entered these genomes as a result of horizontal gene transfer from bacterial sources. One potential case of horizontal gene transfer in A. gossypii and 10 potential cases in S. cerevisiae were identified, of which two were investigated further. One gene, encoding the enzyme dihydroorotate dehydrogenase (DHOD), is potentially a case of horizontal gene transfer, as shown by sequencing of this gene from additional bacterial and fungal species to generate sufficient data to construct a well-supported phylogeny. The DHOD-encoding gene found in S. cerevisiae, URA1 (YKL216W), appears to have entered the Saccharomycetaceae after the divergence of the S. cerevisiae lineage from the Candida albicans lineage and possibly since the divergence from the A. gossypii lineage. This gene appears to have come from the Lactobacillales, and following its acquisition the endogenous eukaryotic DHOD gene was lost. It was also shown that the bacterially derived horizontally transferred DHOD is required for anaerobic synthesis of uracil in S. cerevisiae. The other gene discussed in detail is BDS1, an aryl- and alkyl-sulfatase gene of bacterial origin that we have shown allows utilization of sulfate from several organic sources. Among the eukaryotes, this gene is found in S. cerevisiae and Saccharomyces bayanus and appears to derive from the alpha-proteobacteria.

Mapping of transcription start sites in Saccharomyces cerevisiae using 5' SAGE

A minimally addressed area in Saccharomyces cerevisiae research is the mapping of transcription start sites (TSS). Mapping of TSS in S.cerevisiae has the potential to contribute to our understanding of gene regulation, transcription, mRNA stability and aspects of RNA biology. Here, we use 5' SAGE to map 5' TSS in S.cerevisiae. Tags identifying the first 15-17 bases of the transcripts are created, ligated to form ditags, amplified, concatemerized and ligated into a vector to create a library. Each clone sequenced from this library identifies 10-20 TSS. We have identified 13,746 unique, unambiguous sequence tags from 2231 S.cerevisiae genes. TSS identified in this study are consistent with published results, with primer extension results described here, and are consistent with expectations based on previous work on transcription initiation. We have aligned the sequence flanking 4637 TSS to identify the consensus sequence A(A(rich))5NPyA(A/T)NN(A(rich))6, which confirms and expands the previous reported PyA(A/T)Pu consensus pattern. The TSS data allowed the identification of a previously unrecognized gene, uncovered errors in previous annotation, and identified potential regulatory RNAs and upstream open reading frames in 5'-untranslated region.

Identification of Cryptococcus neoformans temperature-regulated genes with a genomic-DNA microarray

The ability to survive and proliferate at 37 degrees C is an essential virulence attribute of pathogenic microorganisms. A partial-genome microarray was used to profile gene expression in the human-pathogenic fungus Cryptococcus neoformans during growth at 37 degrees C. Genes with orthologs involved in stress responses were induced during growth at 37 degrees C, suggesting that a conserved transcriptional program is used by C. neoformans to alter gene expression during stressful conditions. A gene encoding the transcription factor homolog Mga2 was induced at 37 degrees C and found to be important for high-temperature growth. Genes encoding fatty acid biosynthetic enzymes were identified as potential targets of Mga2, suggesting that membrane remodeling is an important component of adaptation to high growth temperatures. mga2Delta mutants were extremely sensitive to the ergosterol synthesis inhibitor fluconazole, indicating a coordination of the synthesis of membrane component precursors. Unexpectedly, genes involved in amino acid and pyrimidine biosynthesis were repressed at 37 degrees C, but components of these pathways were found to be required for high-temperature growth. Our findings demonstrate the utility of even partial-genome microarrays for delineating regulatory cascades that contribute to microbial pathogenesis.

Convergent evolution of chromosomal sex-determining regions in the animal and fungal kingdoms

Sexual identity is governed by sex chromosomes in plants and animals, and by mating type (MAT) loci in fungi. Comparative analysis of the MAT locus from a species cluster of the human fungal pathogen Cryptococcus revealed sequential evolutionary events that fashioned this large, highly unusual region. We hypothesize that MAT evolved via four main steps, beginning with acquisition of genes into two unlinked sex-determining regions, forming independent gene clusters that then fused via chromosomal translocation. A transitional tripolar intermediate state then converted to a bipolar system via gene conversion or recombination between the linked and unlinked sex-determining regions. MAT was subsequently subjected to intra- and interallelic gene conversion and inversions that suppress recombination. These events resemble those that shaped mammalian sex chromosomes, illustrating convergent evolution in sex-determining structures in the animal and fungal kingdoms.

A comparative genomic analysis of the sex determining region in fungi reveals a remarkable similarity between its evolution and the events which shaped mammalian sex chromosomes

Phylogenetic footprint analysis of IGF2 in extant mammals

Genomic imprinting results in monoallelic gene transcription that is directed by cis-acting regulatory elements epigenetically marked in a parent-of-origin-dependent manner. We performed phylogenetic sequence and epigenetic comparisons of IGF2 between the nonimprinted platypus (Ornithorhynchus anatinus) and imprinted opossum (Didelphis virginiana), mouse (Mus musculus), and human (Homo sapiens) to determine if their divergent imprint status would reflect differences in the conservation of genomic elements important in the regulation of imprinting. We report herein that IGF2 imprinting does not correlate evolutionarily with differential intragenic methylation, nor is it associated with motif 13, a reported IGF2-specific "imprint signature" located in the coding region. Instead, IGF2 imprinting is strongly associated with both the lack of short interspersed transposable elements (SINEs) and an intragenic conserved inverted repeat that contains candidate CTCF-binding sites, a role not previously ascribed to this particular sequence element. Our results are the first to demonstrate that comparative footprint analysis of species from evolutionarily distant mammalian clades, and exhibiting divergent imprint status is a powerful bioinformatics-based approach for identifying cis-acting elements potentially involved not only in the origins of genomic imprinting, but also in its maintenance in humans.

Complete genomic nucleotide sequence of the temperate bacteriophage Aa Phi 23 of Actinobacillus actinomycetemcomitans

The entire double-stranded DNA genome of the Actinobacillus actinomycetemcomitans bacteriophage Aa Phi 23 was sequenced. Linear DNA contained in the phage particles is circularly permuted and terminally redundant. Therefore, the physical map of the phage genome is circular. Its size is 43,033 bp with an overall molar G+C content of 42.5 mol%. Sixty-six potential open reading frames (ORFs) were identified, including an ORF resulting from a translational frameshift. A putative function could be assigned to 23 of them. Twenty-three other ORFs share homologies only with hypothetical proteins present in several bacteria or bacteriophages, and 20 ORFs seem to be specific for phage Aa Phi 23. The organization of the phage genome and several genetic functions share extensive similarities to that of the lambdoid phages. However, Aa Phi 23 encodes a DNA adenine methylase, and the DNA packaging strategy is more closely related to the P22 system. The attachment sites of Aa Phi 23 (attP) and several A. actinomycetemcomitans hosts (attB) are 49 bp long.

The Ashbya gossypii Genome as a Tool for Mapping the Ancient Saccharomyces cerevisiae Genome

We have sequenced and annotated the genome of the filamentous ascomycete Ashbya gossypii. With a size of only 9.2 megabases, encoding 4718 protein-coding genes, it is the smallest genome of a free-living eukaryote yet characterized. More than 90% of A. gossypii genes show both homology and a particular pattern of synteny with Saccharomyces cerevisiae. Analysis of this pattern revealed 300 inversions and translocations that have occurred since divergence of these two species. It also provided compelling evidence that the evolution of S. cerevisiae included a whole genome duplication or fusion of two related species and showed, through inferred ancient gene orders, which of the duplicated genes lost one copy and which retained both copies.

Reinvestigation of the Saccharomyces cerevisiae genome annotation by comparison to the genome of a related fungus: Ashbya gossypii

BACKGROUND

The recently sequenced genome of the filamentous fungus Ashbya gossypii revealed remarkable similarities to that of the budding yeast Saccharomyces cerevisiae both at the level of homology and synteny (conservation of gene order). Thus, it became possible to reinvestigate the S. cerevisiae genome in the syntenic regions leading to an improved annotation.

RESULTS

We have identified 23 novel S. cerevisiae open reading frames (ORFs) as syntenic homologs of A. gossypii genes; for all but one, homologs are present in other eukaryotes including humans. Other comparisons identified 13 overlooked introns and suggested 69 potential sequence corrections resulting in ORF extensions or ORF fusions with improved homology to the syntenic A. gossypii homologs. Of the proposed corrections, 25 were tested and confirmed by resequencing. In addition, homologs of nearly 1,000 S. cerevisiae ORFs, presently annotated as hypothetical, were found in A. gossypii at syntenic positions and can therefore be considered as authentic genes. Finally, we suggest that over 400 S. cerevisiae ORFs that overlap other ORFs in S. cerevisiae and for which no homolog can be detected in A. gossypii should be regarded as spurious.

CONCLUSIONS

Although, the S. cerevisiae genome is rightly considered as one of the most accurately sequenced and annotated eukaryotic genomes, we have shown that it still benefits substantially from comparison to the completed sequence and syntenic gene map of A. gossypii, an evolutionarily related fungus. This type of approach will strongly support the annotation of more complex genomes such as the human and murine genomes.

Verification of a new gene on Saccharomyces cerevisiae chromosome III

We present verification of the existence of a previously unannotated, intron-containing gene of 134 amino acids (predicted molecular weight approximately 15.5 kDa) located on chromosome III of Saccharomyces cerevisiae. Strains carrying a deletion of this gene, which we have called YCL012C, reveal no apparent phenotype. Orthologues of YCL012C are present in related species S. bayanus, S. paradoxus and Ashbya gossypii. Comparison with other fungal sequences reveals that orthologues of this gene are likely present in Schizosaccharomyces pombe, Neurospora crassa and Cryptococcus neoformans as well, indicating that YCL012C is a widely conserved fungal gene.

Mating-type locus of Cryptococcus neoformans: a step in the evolution of sex chromosomes

The sexual development and virulence of the fungal pathogen Cryptococcus neoformans is controlled by a bipolar mating system determined by a single locus that exists in two alleles, alpha and a. The alpha and a mating-type alleles from two divergent varieties were cloned and sequenced. The C. neoformans mating-type locus is unique, spans >100 kb, and contains more than 20 genes. MAT-encoded products include homologs of regulators of sexual development in other fungi, pheromone and pheromone receptors, divergent components of a MAP kinase cascade, and other proteins with no obvious function in mating. The alpha and a alleles of the mating-type locus have extensively rearranged during evolution and strain divergence but are stable during genetic crosses and in the population. The C. neoformans mating-type locus is strikingly different from the other known fungal mating-type loci, sharing features with the self-incompatibility systems and sex chromosomes of algae, plants, and animals. Our study establishes a new paradigm for mating-type loci in fungi with implications for the evolution of cell identity and self/nonself recognition.

The nucleotide sequence of Saccharomyces cerevisiae chromosome V

Here we report the sequence of 569,202 base pairs of Saccharomyces cerevisiae chromosome V. Analysis of the sequence revealed a centromere, two telomeres and 271 open reading frames (ORFs) plus 13 tRNAs and four small nuclear RNAs. There are two Tyl transposable elements, each of which contains an ORF (included in the count of 271). Of the ORFs, 78 (29%) are new, 81 (30%) have potential homologues in the public databases, and 112 (41%) are previously characterized yeast genes.

The nucleotide sequence of Saccharomyces cerevisiae chromosome XVI

The nucleotide sequence of the 948,061 base pairs of chromosome XVI has been determined, completing the sequence of the yeast genome. Chromosome XVI was the last yeast chromosome identified, and some of the genes mapped early to it, such as GAL4, PEP4 and RAD1 (ref. 2) have played important roles in the development of yeast biology. The architecture of this final chromosome seems to be typical of the large yeast chromosomes, and shows large duplications with other yeast chromosomes. Chromosome XVI contains 487 potential protein-encoding genes, 17 tRNA genes and two small nuclear RNA genes; 27% of the genes have significant similarities to human gene products, and 48% are new and of unknown biological function. Systematic efforts to explore gene function have begun.