Evidence from comparative genomics for a complete sexual cycle in the 'asexual' pathogenic yeast Candida glabrata

Tn-seq of the Candida glabrata reference strain CBS138 reveals epigenetic plasticity, structural variation, and intrinsic mechanisms of resistance to micafungin

C. glabrata is an opportunistic pathogen that can resist common antifungals and rapidly acquire multidrug resistance. A large amount of genetic variation exists between isolates, which complicates generalizations. Portable Tn-seq methods can efficiently provide genome-wide information on strain differences and genetic mechanisms. Using the Hermes transposon, the CBS138 reference strain and a commonly studied derivative termed 2001 were subjected to Tn-seq in control conditions and after exposure to varying doses of the clinical antifungal micafungin. The approach revealed large differences between these strains, including a 131 kb tandem duplication and a variety of fitness differences. Additionally, both strains exhibited up to 1000-fold increased transposon accessibility in subtelomeric regions relative to the BG2 strain, indicative of open subtelomeric chromatin in these isolates and large epigenetic variation within the species. Unexpectedly, the Pdr1 transcription factor conferred resistance to micafungin through targets other than CDR1 . Other micafungin resistance pathways were also revealed including mannosyltransferase activity and biosynthesis of the lipid precursor sphingosine, the drugging of which by SDZ 90-215 or myriocin enhanced the potency of micafungin in vitro . These findings provide insights into complexity of the C. glabrata species as well as strategies for improving antifungal efficacy.

Summary

Candida glabrata is an emerging pathogen with large genetic diversity and genome plasticity. The type strain CBS138 and a laboratory derivative were mutagenized with the Hermes transposon and profiled using Tn-seq. Numerous genes that regulate innate and acquired resistance to an important clinical antifungal were uncovered, including a pleiotropic drug resistance gene (PDR1) and a duplication of part of one chromosome. Compounds that target PDR1 and other genes may augment the potency of existing antifungals.

Is There a Relationship Between Mating and Pathogenesis in Two Human Fungal Pathogens, Candida albicans and Candida glabrata?

Purpose of Review

Human fungal pathogens are rapidly increasing in incidence and readily able to evade the host immune responses. Our ability to study the genetic behind this has been limited due to the apparent lack of a sexual cycle and forward genetic tools. In this review, we discuss the evolution of mating, meiosis, and pathogenesis and if these processes are advantageous to pathogens.

Recent Findings

This review summarises what is currently known about the sexual cycles of two important human fungal pathogens, Candida albicans and Candida glabrata. This includes the identification of parasexual cycle in C. albicans and the observed low levels of recombination in C. glabrata populations.

Summary

In this review, we present what is currently known about the mating types and mating/sexual cycles of two clinically important human fungal pathogens, Candida albicans and Candida glabrata. We discuss the evolution of meiosis using the knowledge that has been amassed from the decades of studying Saccharomyces cerevisiae and how this can be applied to fungal pathogens. We further discuss how the evolution of pathogenesis has played a role in influencing mating processes in human fungal pathogens and compare sexual cycles between C. albicans and C. glabrata, highlighting knowledge gaps and suggesting how these two fungi have evolved distinct mating niches to allow the development of disease in a human host.

Ploidy Variation and Spontaneous Haploid-Diploid Switching of Candida glabrata Clinical Isolates

The human fungal pathogen Candida glabrata is phylogenetically closely related to Saccharomyces cerevisiae, a model eukaryotic organism. Unlike S. cerevisiae, which has both haploid and diploid forms and a complete sexual cycle, C. glabrata has long been considered a haploid and asexual species. In this study, we analyzed the ploidy states of 500 clinical isolates of C. glabrata from four Chinese hospitals and found that approximately 4% of the isolates were in or able to spontaneously switch to an aneuploid (genomic DNA, 1N-2N), diploid (2N), or hyperdiploid (>2N) form under in vivo or in vitro conditions. Stable diploid cells were identified in 3% of the isolates (15/500). Of particular interest, one clinical strain existed only in the diploid form. Multilocus sequence typing (MLST) assays revealed two major genetic clusters (A and B) of C. glabrata isolates. Most of the isolates (70%) from China belonged to the A cluster, whereas most of the isolates from other countries (such as Iran, Japan, United States, and European countries) belonged to the B cluster. Further investigation indicated that C. glabrata cells of different ploidy forms differed in a number of respects, including morphologies, antifungal susceptibility, virulence, and global gene expression profiles. Additionally, C. glabrata could undergo spontaneous switching between the diploid and haploid forms under both in vitro and in vivo conditions. Given the absence of an apparent sexual phase, one would expect that the ploidy shifts could function as an alternative strategy that promotes genetic diversity and benefits the ability of the fungus to rapidly adapt to the changing environment. IMPORTANCE The human fungal pathogen Candida glabrata has long been thought to be a haploid organism. Here, we report the population structure and ploidy states of 500 clinical isolates of C. glabrata from China. To our surprise, we found that the ploidy of a subset of clinical isolates varied dramatically. Some isolates were in or able to switch to an aneuploid, diploid, or hyperdiploid form. C. glabrata cells with different ploidy differed in a number of biological respects, including morphology, antifungal susceptibility, virulence, and global gene expression profile. Given the absence of an apparent sexual phase in this fungus, we propose that ploidy switching could be a strategy for rapid adaptation to environmental changes and could function as an alternative to sexual reproduction.

Population genetics and microevolution of clinical Candida glabrata reveals recombinant sequence types and hyper-variation within mitochondrial genomes, virulence genes, and drug targets

Candida glabrata is the second most common etiological cause of worldwide systemic candidiasis in adult patients. Genome analysis of 68 isolates from 8 hospitals across Scotland, together with 83 global isolates, revealed insights into the population genetics and evolution of C. glabrata. Clinical isolates of C. glabrata from across Scotland are highly genetically diverse, including at least 19 separate sequence types that have been recovered previously in globally diverse locations, and 1 newly discovered sequence type. Several sequence types had evidence for ancestral recombination, suggesting transmission between distinct geographical regions has coincided with genetic exchange arising in new clades. Three isolates were missing MATα1, potentially representing a second mating type. Signatures of positive selection were identified in every sequence type including enrichment for epithelial adhesins thought to facilitate fungal adhesin to human epithelial cells. In patent microevolution was identified from 7 sets of recurrent cases of candidiasis, revealing an enrichment for nonsynonymous and frameshift indels in cell surface proteins. Microevolution within patients also affected epithelial adhesins genes, and several genes involved in drug resistance including the ergosterol synthesis gene ERG4 and the echinocandin target FKS1/2, the latter coinciding with a marked drop in fluconazole minimum inhibitory concentration. In addition to nuclear genome diversity, the C. glabrata mitochondrial genome was particularly diverse, with reduced conserved sequence and conserved protein-encoding genes in all nonreference ST15 isolates. Together, this study highlights the genetic diversity within the C. glabrata population that may impact virulence and drug resistance, and 2 major mechanisms generating this diversity: microevolution and genetic exchange/recombination.

Genome plasticity in Candida albicans: A cutting-edge strategy for evolution, adaptation, and survival

Candida albicans is the most implicated fungal species that grows as a commensal or opportunistic pathogen in the human host. It is associated with many life-threatening infections, especially in immunocompromised persons. The genome of Candida albicans is very flexible and can withstand a wide assortment of variations in a continuously changing environment. Thus, genome plasticity is central to its adaptation and has long been of considerable interest. C. albicans has a diploid heterozygous genome that is highly dynamic and can display variation from small to large scale chromosomal rearrangement and aneuploidy, which have implications in drug resistance, virulence, and pathogenicity. This review presents an up-to-date overview of recent genomic studies involving C. albicans. It discusses the accumulating evidence that shows how mitotic recombination events, ploidy dynamics, aneuploidy, and loss of heterozygosity (LOH) influence evolution, adaptation, and survival in C. albicans. Understanding the factors that affect the genome is crucial for a proper understanding of species and rapid development and adjustment of therapeutic strategies to mitigate their spread.

Using Synthetic Genetic Interactions in Candida glabrata as a Novel Method to Detect Genes with Roles in Antifungal Drug Resistance

Synthetic genetic interaction analysis is a powerful genetic strategy that analyzes the fitness and phenotypes of single- and double-gene mutant cells in order to dissect the interactions between genes, categorize into biological pathways, and characterize genes of unknown function. It has been extensively employed in model organisms for fundamental, systems-level assessment of the interactions between genes. However, more recently, genetic interaction mapping has been applied to fungal pathogens and has been instrumental for the study of clinically important infectious organisms. This protocol herein explains in the detail the methodology and analysis that can be employed to develop interaction maps in microbial pathogens. Such techniques can aid in bridging our understanding of complex genetic networks, with applications to diverse microbial pathogens to further our understanding of virulence, the use of antimicrobial therapies, and host-pathogen interactions.

DNA damage response of major fungal pathogen Candida glabrata offers clues to explain its genetic diversity

How cells respond to DNA damage is key to maintaining genome integrity or facilitating genetic change. In fungi, DNA damage responses have been extensively characterized in the model budding yeast Saccharomyces cerevisiae, which is generally not pathogenic. However, it is not clear how closely these responses resemble those in fungal pathogens, in which genetic change plays an important role in the evolutionary arms race between pathogen and host and the evolution of antifungal drug resistance. A close relative of S. cerevisiae, Candida glabrata, is an opportunistic pathogen that displays high variability in chromosome structure among clinical isolates and rapidly evolves antifungal drug resistance. The mechanisms facilitating such genomic flexibility and evolvability in this organism are unknown. Recently we characterized the DNA damage response of C. glabrata and identified several features that distinguish it from the well characterized DNA damage response of S. cerevisiae. First, we discovered that, in contrast to the established paradigm, C. glabrata effector kinase Rad53 is not hyperphosphorylated upon DNA damage. We also uncovered evidence of an attenuated DNA damage checkpoint response, wherein in the presence of DNA damage C. glabrata cells did not accumulate in S-phase and proceeded with cell division, leading to aberrant mitoses and cell death. Finally, we identified evidence of transcriptional rewiring of the DNA damage response of C. glabrata relative to S. cerevisiae, including an upregulation of genes involved in mating and meiosis-processes that have not been reported in C. glabrata. Together, these results open new possibilities and raise tantalizing questions of how this major fungal pathogen facilitates genetic change.

The genome of Geosiphon pyriformis reveals ancestral traits linked to the emergence of the arbuscular mycorrhizal symbiosis

Arbuscular mycorrhizal fungi (AMF) (subphylum Glomeromycotina)¹ are among the most prominent symbionts and form the Arbuscular Mycorrhizal symbiosis (AMS) with over 70% of known land plants.^2,3 AMS allows plants to efficiently acquire poorly soluble soil nutrients⁴ and AMF to receive photosynthetically fixed carbohydrates. This plant-fungus symbiosis dates back more than 400 million years⁵ and is thought to be one of the key innovations that allowed the colonization of lands by plants.⁶ Genomic and genetic analyses of diverse plant species started to reveal the molecular mechanisms that allowed the evolution of this symbiosis on the host side, but how and when AMS abilities emerged in AMF remain elusive. Comparative phylogenomics could be used to understand the evolution of AMS.^7,8 However, the availability of genome data covering basal AMF phylogenetic nodes (Archaeosporales, Paraglomerales) is presently based on fragmentary protein coding datasets.⁹Geosiphon pyriformis (Archaeosporales) is the only fungus known to produce endosymbiosis with nitrogen-fixing cyanobacteria (Nostoc punctiforme) presumably representing the ancestral AMF state.^10-12 Unlike other AMF, it forms long fungal cells ("bladders") that enclose cyanobacteria. Once in the bladder, the cyanobacteria are photosynthetically active and fix nitrogen, receiving inorganic nutrients and water from the fungus. Arguably, G. pyriformis represents an ideal candidate to investigate the origin of AMS and the emergence of a unique endosymbiosis. Here, we aimed to advance knowledge in these questions by sequencing the genome of G. pyriformis, using a re-discovered isolate.

A Noncanonical DNA Damage Checkpoint Response in a Major Fungal Pathogen

In order to preserve genome integrity, all cells must mount appropriate responses to DNA damage, including slowing down or arresting the cell cycle to give the cells time to repair the damage and changing gene expression, for example to induce genes involved in DNA repair. The Rad53 protein kinase is a conserved central mediator of these responses in eukaryotic cells, and its extensive phosphorylation upon DNA damage is necessary for its activation and subsequent activity.

DNA damage checkpoints are key guardians of genome integrity. Eukaryotic cells respond to DNA damage by triggering extensive phosphorylation of Rad53/CHK2 effector kinase, whereupon activated Rad53/CHK2 mediates further aspects of checkpoint activation, including cell cycle arrest and transcriptional changes. Budding yeast Candida glabrata, closely related to model eukaryote Saccharomyces cerevisiae, is an opportunistic pathogen characterized by high genetic diversity and rapid emergence of drug-resistant mutants. However, the mechanisms underlying this genetic variability are unclear. We used Western blotting and mass spectrometry to show that, unlike S. cerevisiae, C. glabrata cells exposed to DNA damage did not induce C. glabrata Rad53 (CgRad53) phosphorylation. Furthermore, flow cytometry analysis showed that, unlike S. cerevisiae, C. glabrata cells did not accumulate in S phase upon DNA damage. Consistent with these observations, time-lapse microscopy showed C. glabrata cells continuing to divide in the presence of DNA damage, resulting in mitotic errors and cell death. Finally, transcriptome sequencing (RNAseq) analysis revealed transcriptional rewiring of the DNA damage response in C. glabrata and identified several key protectors of genome stability upregulated by DNA damage in S. cerevisiae but downregulated in C. glabrata, including proliferating cell nuclear antigen (PCNA). Together, our results reveal a noncanonical fungal DNA damage response in C. glabrata, which may contribute to rapidly generating genetic change and drug resistance.

Genome editing in the yeast Nakaseomyces delphensis and description of its complete sexual cycle

The environmental yeast Nakaseomyces delphensis is, phylogenetically, the closest known species to Candida glabrata, a major fungal pathogen of humans. C. glabrata is haploid and described as asexual, while N. delphensis is also haploid, but has been described as competent for mating and meiosis. Both genomes contain homologues of all the genes necessary for sexual reproduction and also the genes for Ho-dependent mating-type switching, like Saccharomyces cerevisiae. We first report the construction of genetically engineered strains of N. delphensis, including by CRISPR-Cas 9 gene editing. We also report the description of the sexual cycle of N. delphensis. We show that it undergoes Ho-dependent mating-type switching in culture and that deletion of the HO gene prevents such switching and allows maintenance of stable, separate, MATa and MATalpha haploid strains. Rare, genetically selected diploids can be obtained through mating of haploid strains, mutated or not for the HO gene. In contrast to HO/HO diploids, which behave as expected, Δho/Δho diploids exhibit unusual profiles in flow cytometry. Both types of diploids can produce recombined haploid cells, which grow like the original haploid-type strain. Our experiments thus allow the genetic manipulation of N. delphensis and the reconstruction, in the laboratory, of its entire life cycle.

Towards yeast taxogenomics: lessons from novel species descriptions based on complete genome sequences

In recent years, ‘multi-omic’ sciences have affected all aspects of fundamental and applied biological research. Yeast taxonomists, though somewhat timidly, have begun to incorporate complete genomic sequences into the description of novel taxa, taking advantage of these powerful data to calculate more reliable genetic distances, construct more robust phylogenies, correlate genotype with phenotype and even reveal cryptic sexual behaviors. However, the use of genomic data in formal yeast species descriptions is far from widespread. The present review examines published examples of genome-based species descriptions of yeasts, highlights relevant bioinformatic approaches, provides recommendations for new users and discusses some of the challenges facing the genome-based systematics of yeasts.

The Mechanisms of Mating in Pathogenic Fungi-A Plastic Trait

The impact of fungi on human and plant health is an ever-increasing issue. Recent studies have estimated that human fungal infections result in an excess of one million deaths per year and plant fungal infections resulting in the loss of crop yields worth approximately 200 million per annum. Sexual reproduction in these economically important fungi has evolved in response to the environmental stresses encountered by the pathogens as a method to target DNA damage. Meiosis is integral to this process, through increasing diversity through recombination. Mating and meiosis have been extensively studied in the model yeast Saccharomyces cerevisiae, highlighting that these mechanisms have diverged even between apparently closely related species. To further examine this, this review will inspect these mechanisms in emerging important fungal pathogens, such as Candida, Aspergillus, and Cryptococcus. It shows that both sexual and asexual reproduction in these fungi demonstrate a high degree of plasticity.

Molecular Typing of Candida glabrata

The yeast Candida glabrata has emerged, second only to Candida albicans, to be one of the most frequently isolated fungi in clinical specimen from human. Its frequent resistance towards azole antifungal drugs and the high capacity to form biofilms on indwelling catheters of individual isolates render it an often difficult to treat pathogen. Hence, there is a notably increasing scientific and clinical interest in this species. This has led to the development of a variety of molecular tools for genetic modification, strain collections, and last but not least different approaches to analyse the population structure among isolates of different geographical and clinical contexts. Often, these are used to study correlations (or the absence thereof) with different pathogenicity, virulence, or drug resistance traits. Three molecular methods have been used to type within the C. glabrata population on a genetic level by multiple studies: multi-locus sequence typing, microsatellite length polymorphisms, and clustering of whole-genome sequencing data, and these are subject of this review.

Mating genotypes and susceptibility profiles of clinical isolates of Candida glabrata from Turkey

The sexual cycle of Candida glabrata is not known; however, genomic evidence is indicative of recombination among subpopulations and the genome harbours genes necessary for undergoing mating and meiosis, which may increase fitness. The relationship between specific mating type-like (MTL) loci and antifungal susceptibility is not well understood in C. glabrata. We investigated different combinations of clinical C. glabrata isolate mating types and their antifungal susceptibility profiles. Allele profiles of the mating genes of 103 clinical C. glabrata isolates were identified, and their antifungal susceptibility to azoles, echinocandins and amphotericin B were compared. The majority (88.3%) of screened isolates harboured the a allele in the locus. The MTL1, MTL2 and MTL3 loci harboured a (88.3%), a (95.1%), and α (71.8%) alleles, respectively. The C. glabrata isolates were susceptible to echinocandins but displayed high minimal inhibitory concentrations (MICs) for azoles. The MIC ranges and MIC90 values of all isolates were 1.0 to ≥64 and 8.0 μg mL^-1 for fluconazole, 0.06 to ≥16.0 and 0.5 μg mL^-1 for voriconazole, 0.06 to ≥16.0 and 1.0 μg mL^-1 for posaconazole, ≤0.015 to 0.06, and 0.03 μg mL^-1 for caspofungin, ≤0.015 to 0.06 and 0.015 μg mL^-1 for anidulafungin and 0.5-2 and 2.0 μg mL^-1 for amphotericin B, respectively. The mating gene alleles of the clinical C. glabrata isolates were not associated with differences in the MICs of the tested antifungals, except for the MTL3 α-allele and echinocandins. The mating genotypes of the clinical C. glabrata isolates had no recognisable common effect on antifungal susceptibility.

The phosphorelay signal transduction system in Candida glabrata: an in silico analysis

Signaling systems allow microorganisms to sense and respond to different stimuli through the modification of gene expression. The phosphorelay signal transduction system in eukaryotes involves three proteins: a sensor protein, an intermediate protein and a response regulator, and requires the transfer of a phosphate group between two histidine-aspartic residues. The SLN1-YPD1-SSK1 system enables yeast to adapt to hyperosmotic stress through the activation of the HOG1-MAPK pathway. The genetic sequences available from Saccharomyces cerevisiae were used to identify orthologous sequences in Candida glabrata, and putative genes were identified and characterized by in silico assays. An interactome analysis was carried out with the complete genome of C. glabrata and the putative proteins of the phosphorelay signal transduction system. Next, we modeled the complex formed between the sensor protein CgSln1p and the intermediate CgYpd1p. Finally, phosphate transfer was examined by a molecular dynamic assay. Our in silico analysis showed that the putative proteins of the C. glabrata phosphorelay signal transduction system present the functional domains of histidine kinase, a downstream response regulator protein, and an intermediate histidine phosphotransfer protein. All the sequences are phylogenetically more related to S. cerevisiae than to C. albicans. The interactome suggests that the C. glabrata phosphorelay signal transduction system interacts with different proteins that regulate cell wall biosynthesis and responds to oxidative and osmotic stress the same way as similar systems in S. cerevisiae and C. albicans. Molecular dynamics simulations showed complex formation between the response regulator domain of histidine kinase CgSln1 and intermediate protein CgYpd1 in the presence of a phosphate group and interactions between the aspartic residue and the histidine residue. Overall, our research showed that C. glabrata harbors a functional SLN1-YPD1-SSK1 phosphorelay system.

Evolution of sexual reproduction: a view from the Fungal Kingdom supports an evolutionary epoch with sex before sexes

Sexual reproduction is conserved throughout each supergroup within the eukaryotic tree of life, and therefore thought to have evolved once and to have been present in the last eukaryotic common ancestor (LECA). Given the antiquity of sex, there are features of sexual reproduction that are ancient and ancestral, and thus shared in diverse extant organisms. On the other hand, the vast evolutionary distance that separates any given extant species from the LECA necessarily implies that other features of sex will be derived. While most types of sex we are familiar with involve two opposite sexes or mating types, recent studies in the fungal kingdom have revealed novel and unusual patterns of sexual reproduction, including unisexual reproduction. In this mode of reproduction a single mating type can on its own undergo self-fertile/homothallic reproduction, either with itself or with other members of the population of the same mating type. Unisexual reproduction has arisen independently as a derived feature in several different lineages. That a myriad of different types of sex determination and sex determinants abound in animals, plants, protists, and fungi suggests that sex specification itself may not be ancestral and instead may be a derived trait. If so, then the original form of sexual reproduction may have been unisexual, onto which sexes were superimposed as a later feature. In this model, unisexual reproduction is both an ancestral and a derived trait. In this review, we consider what is new and what is old about sexual reproduction from the unique vantage point of the fungal kingdom.

Phylogenetic and Transcripts Profiling of Glucose Sensing Related Genes in Candida glabrata

Background:

The sensing mechanism of glucose in Saccharomyces cerevisiae is well studied. However, such information is scarcely found in other yeast species such as Candida glabrata.

Objectives:

This study aimed to identify the glucose sensing pathway related genes of C. glabrata and to analyze the regulation pattern of these genes in response to different surrounding glucose concentrations through the quantitative real time polymerase chain reaction (qRT-PCR).

Materials and Methods:

Phylogenetic analysis was carried out on predicted amino acid sequences of C. glabrata and S. cerevisiae to compare their degree of similarity. In addition, the growth of C. glabrata in response to different amounts of glucose (0%, 0.01%, 0.1%, 1% and 2%) was evaluated via the spot dilution assay on prepared agar medium. Besides, the SNF3 and RGT2, which act as putative glucose sensors, and the RGT1 and MIG1, which act as putative transcriptional regulators and selected downstream hexose transporters (HXTs), were analysed through qRT-PCR analysis for the gene expression level under different glucose concentrations.

Results:

Comparative analysis of predicted amino acids in the phylogenetic tree showed high similarity between C. glabrata and S cerevisiae. Besides, C. glabrata demonstrated the capability to grow in glucose levels as low as 0.01% in the spot dilution assay. In qRT-PCR analysis, differential expressions were observed in selected genes when C. glabrata was subjected to different glucose concentrations.

Conclusions:

The constructed phylogenetic tree suggests the close evolutionary relationship between C. glabrata and S. cerevisiae. The capability of C. glabrata to grow in extremely low glucose environments and the differential expression of selected glucose-sensing related genes suggested the possible role of these genes in modulating the growth of C. glabrata in response to different glucose concentrations. This study helps deepen our understanding of the glucose sensing mechanism in C. glabrata and serves to provide fundamental data that may assist in unveiling this mechanism as a potential drug target.

Efficient Mating-Type Switching in Candida glabrata Induces Cell Death

Candida glabrata is an apparently asexual haploid yeast that is phylogenetically closer to Saccharomyces cerevisiae than to Candida albicans. Its genome contains three MAT-like cassettes, MAT, which encodes either MATa or MATalpha information in different strains, and the additional loci, HML and HMR. The genome also contains an HO gene homolog, but this yeast has never been shown to switch mating-types spontaneously, as S. cerevisiae does. We have recently sequenced the genomes of the five species that, together with C. glabrata, make up the Nakaseomyces clade. All contain MAT-like cassettes and an HO gene homolog. In this work, we express the HO gene of all Nakaseomyces and of S. cerevisiae in C. glabrata. All can induce mating-type switching, but, despite the larger phylogenetic distance, the most efficient endonuclease is the one from S. cerevisiae. Efficient mating-type switching in C. glabrata is accompanied by a high cell mortality, and sometimes results in conversion of the additional cassette HML. Mortality probably results from the cutting of the HO recognition sites that are present, in HML and possibly HMR, contrary to what happens naturally in S. cerevisiae. This has implications in the life-cycle of C. glabrata, as we show that efficient MAT switching is lethal for most cells, induces chromosomal rearrangements in survivors, and that the endogenous HO is probably rarely active indeed.

Utilising established SDL-screening methods as a tool for the functional genomic characterisation of model and non-model organisms

The trend for large-scale genetic and phenotypic screens has revealed a wealth of information on biological systems. A major challenge is understanding how genes function and putative roles in networks. The majority of current gene knowledge is garnered from studies utilising the model yeast Saccharomyces cerevisiae. We demonstrate that synthetic dosage lethal genetic array methodologies can be used to study genetic networks in other yeasts, namely the fungal pathogen Candida glabrata, which has limited forward genetic tools, due to the lack of 'natural' mating. We performed two SDL screens in S. cerevisiae, overexpressing the transcriptional regulator UME6 as bait in the first screen and its C. glabrata ortholog CAGL0F05357g in the second. Analysis revealed that SDL maps share 204 common interactors, with 10 genetic interactions unique to C. glabrata indicating a level of genetic rewiring, indicative of linking genotype to phenotype in fungal pathogens. This was further validated by incorporating our results into the global genetic landscape map of the cell from Costanzo et al. to identify common and novel gene attributes. This data demonstrated the utility large data sets and more robust analysis made possible by interrogating exogenous genes in the context of the eukaryotic global genetic landscape.

Evolution of regulatory networks in Candida glabrata: learning to live with the human host

The opportunistic human fungal pathogen Candida glabrata is second only to C. albicans as the cause of Candida infections and yet is more closely related to Saccharomyces cerevisiae. Recent advances in functional genomics technologies and computational approaches to decipher regulatory networks, and the comparison of these networks among these and other Ascomycete species, have revealed both unique and shared strategies in adaptation to a human commensal/opportunistic pathogen lifestyle and antifungal drug resistance in C. glabrata. Recently, several C. glabrata sister species in the Nakeseomyces clade representing both human associated (commensal) and environmental isolates have had their genomes sequenced and analyzed. This has paved the way for comparative functional genomics studies to characterize the regulatory networks in these species to identify informative patterns of conservation and divergence linked to phenotypic evolution in the Nakaseomyces lineage.

Comparative genome analysis of Pseudogymnoascus spp. reveals primarily clonal evolution with small genome fragments exchanged between lineages

Background

Pseudogymnoascus spp. is a wide group of fungi lineages in the family Pseudorotiaceae including an aggressive pathogen of bats P. destructans. Although several lineages of P. spp. were shown to produce ascospores in culture, the vast majority of P. spp. demonstrates no evidence of sexual reproduction. P. spp. can tolerate a wide range of different temperatures and salinities and can survive even in permafrost layer. Adaptability of P. spp. to different environments is accompanied by extremely variable morphology and physiology.

Results

We sequenced genotypes of 14 strains of P. spp., 5 of which were extracted from permafrost, 1 from a cryopeg, a layer of unfrozen ground in permafrost, and 8 from temperate surface environments. All sequenced genotypes are haploid. Nucleotide diversity among these genomes is very high, with a typical evolutionary distance at synonymous sites dS ≈ 0.5, suggesting that the last common ancestor of these strains lived >50Mya. The strains extracted from permafrost do not form a separate clade. Instead, each permafrost strain has close relatives from temperate environments.

We observed a strictly clonal population structure with no conflicting topologies for ~99% of genome sequences. However, there is a number of short (~100–10,000 nt) genomic segments with the total length of 67.6 Kb which possess phylogenetic patterns strikingly different from the rest of the genome. The most remarkable case is a MAT-locus, which has 2 distinct alleles interspersed along the whole-genome phylogenetic tree.

Conclusions

Predominantly clonal structure of genome sequences is consistent with the observations that sexual reproduction is rare in P. spp. Small number of regions with noncanonical phylogenies seem to arise due to some recombination events between derived lineages of P. spp., with MAT-locus being transferred on multiple occasions. All sequenced strains have heterothallic configuration of MAT-locus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1570-9) contains supplementary material, which is available to authorized users.

Parasexual Ploidy Reduction Drives Population Heterogeneity Through Random and Transient Aneuploidy in Candida albicans

The opportunistic pathogen Candida albicans has a large repertoire of mechanisms to generate genetic and phenotypic diversity despite the lack of meiosis in its life cycle. Its parasexual cycle enables shifts in ploidy, which in turn facilitate recombination, aneuploidy, and homozygosis of whole chromosomes to fuel rapid adaptation. Here we show that the tetraploid state potentiates ploidy variation and drives population heterogeneity. In tetraploids, the rate of losing a single heterozygous marker [loss of heterozygosity (LOH)] is elevated ∼30-fold higher than the rate in diploid cells. Furthermore, isolates recovered after selection for LOH of one, two, or three markers were highly aneuploid, with a broad range of karyotypes including strains with a combination of di-, tri-, and tetrasomic chromosomes. We followed the ploidy trajectories for these tetraploid- and aneuploid-derived isolates, using a combination of flow cytometry and double-digestion restriction-site-associated DNA analyzed with next-generation sequencing. Isolates derived from either tetraploid or aneuploid isolates predominately resolved to a stable euploid state. The majority of isolates reduced to the conventional diploid state; however, stable triploid and tetraploid states were observed in ∼30% of the isolates. Notably, aneuploid isolates were more transient than tetraploid isolates, resolving to a euploid state within a few passages. Furthermore, the likelihood that a particular isolate will resolve to the same ploidy state in replicate evolution experiments is only ∼50%, supporting the idea that the chromosome loss process of the parasexual cycle is random and does not follow trajectories involving specific combinations of chromosomes. Together, our results indicate that tetraploid progenitors can produce populations of progeny cells with a high degree of genomic diversity, from altered ploidy to homozygosis, providing an excellent source of genetic variation upon which selection can act.

Comparative phenotypic analysis of the major fungal pathogens Candida parapsilosis and Candida albicans

Candida parapsilosis and Candida albicans are human fungal pathogens that belong to the CTG clade in the Saccharomycotina. In contrast to C. albicans, relatively little is known about the virulence properties of C. parapsilosis, a pathogen particularly associated with infections of premature neonates. We describe here the construction of C. parapsilosis strains carrying double allele deletions of 100 transcription factors, protein kinases and species-specific genes. Two independent deletions were constructed for each target gene. Growth in >40 conditions was tested, including carbon source, temperature, and the presence of antifungal drugs. The phenotypes were compared to C. albicans strains with deletions of orthologous transcription factors. We found that many phenotypes are shared between the two species, such as the role of Upc2 as a regulator of azole resistance, and of CAP1 in the oxidative stress response. Others are unique to one species. For example, Cph2 plays a role in the hypoxic response in C. parapsilosis but not in C. albicans. We found extensive divergence between the biofilm regulators of the two species. We identified seven transcription factors and one protein kinase that are required for biofilm development in C. parapsilosis. Only three (Efg1, Bcr1 and Ace2) have similar effects on C. albicans biofilms, whereas Cph2, Czf1, Gzf3 and Ume6 have major roles in C. parapsilosis only. Two transcription factors (Brg1 and Tec1) with well-characterized roles in biofilm formation in C. albicans do not have the same function in C. parapsilosis. We also compared the transcription profile of C. parapsilosis and C. albicans biofilms. Our analysis suggests the processes shared between the two species are predominantly metabolic, and that Cph2 and Bcr1 are major biofilm regulators in C. parapsilosis.

Candida species are among the most common causes of fungal infection worldwide. Infections can be both community-based and hospital-acquired, and are particularly associated with immunocompromised individuals. Candida albicans is the most commonly isolated species and is the best studied. However, other species are becoming of increasing concern. Candida parapsilosis causes outbreaks of infection in neonatal wards, and is one of the few Candida species that is transferred from the hands of healthcare workers. C. parapsilosis, like C. albicans, grows as biofilms (cell communities) on the surfaces of indwelling medical devices like feeding tubes. We describe here the construction of a set of tools that allow us to characterize the virulence properties of C. parapsilosis, and in particular its ability to grow as biofilms. We find that some of the regulatory mechanisms are shared with C. albicans, but others are unique to each species. Our tools, based on selectively deleting regulatory genes, will provide a major resource to the fungal research community.

The Candida pathogenic species complex

Candida species are the most common causes of fungal infection. Approximately 90% of infections are caused by five species: Candida albicans, Candida glabrata, Candida tropicalis, Candida parapsilosis, and Candida krusei. Three (C. albicans, C. tropicalis, and C. parapsilosis) belong to the CTG clade, in which the CTG codon is translated as serine and not leucine. C. albicans remains the most commonly isolated but is decreasing relative to the other species. The increasing incidence of C. glabrata is related to its reduced susceptibility to azole drugs. Genome analysis suggests that virulence in the CTG clade is associated with expansion of gene families, particularly of cell wall genes. Similar independent processes took place in the C. glabrata species group. Gene loss and expansion in an ancestor of C. glabrata may have resulted in preadaptations that enabled pathogenicity.

Candida glabrata: a deadly companion?

The yeast Candida glabrata has become a major fungal opportunistic pathogen of humans since the 1980s. Contrary to what its name suggests, it is much closer, phylogenetically, to the model yeast Saccharomyces cerevisiae than to the most prevalent human fungal pathogen, Candida albicans. Its similarity to S. cerevisiae fortunately extends to their amenability to molecular genetics methods. C. glabrata is now described as part of the Nakaseomyces clade, which includes two new pathogens and other environmental species. C. glabrata is likely a commensal species of the human digestive tract, but systemic infections of immunocompromised patients are often fatal. In addition to being the subject of active medical research, other studies on C. glabrata focus on fundamental aspects of evolution of yeast genomes and adaptation. For example, the genome of C. glabrata has undergone major gene and intron loss compared to S. cerevisiae. It is also an apparently asexual species, a feature that inevitably leads to questions about the species' evolutionary past, present and future. On-going research with this yeast continues to address various aspects of adaptation to the human host and mechanisms of evolution in the Saccharomycetaceae, major model organisms for biology.

Sex or no sex: evolutionary adaptation occurs regardless

All species continuously evolve to adapt to changing environments. The genetic variation that fosters such adaptation is caused by a plethora of mechanisms, including meiotic recombination that generates novel allelic combinations in the progeny of two parental lineages. However, a considerable number of eukaryotic species, including many fungi, do not have an apparent sexual cycle and are consequently thought to be limited in their evolutionary potential. As such organisms are expected to have reduced capability to eliminate deleterious mutations, they are often considered as evolutionary dead ends. However, inspired by recent reports we argue that such organisms can be as persistent as organisms with conventional sexual cycles through the use of other mechanisms, such as genomic rearrangements, to foster adaptation.

The role of phenotypic switching in the basic biology and pathogenesis of Candida albicans

The "white-opaque" transition in Candida albicans was discovered in 1987. For the next fifteen years, a significant body of knowledge accumulated that included differences between the cell types in gene expression, cellular architecture and virulence in cutaneous and systemic mouse models. However, it was not until 2002 that we began to understand the role of switching in the life history of this pathogen, the role of the mating type locus and the molecular pathways that regulated it. Then in 2006, both the master switch locus WORI and the pheromone-induced white cell biofilm were discovered. Since that year, a number of new observations on the regulation and biology of switching have been made that have significantly increased the perceived complexity of this fascinating phenotypic transition.

ER stress response mechanisms in the pathogenic yeast Candida glabrata and their roles in virulence

The maintenance of endoplasmic reticulum (ER) homeostasis is critical for numerous aspects of cell physiology. Eukaryotic cells respond to the accumulation of misfolded proteins in the ER (ER stress) by activating the unfolded protein response (UPR), an intracellular signaling pathway that adjusts the folding capacity of the ER. Recent studies of several pathogenic fungi have revealed that the UPR is important for antifungal resistance and virulence; therefore, the pathway has attracted much attention as a potential therapeutic target. While the UPR is highly conserved among eukaryotes, our group recently discovered that the pathogenic yeast Candida glabrata lacks the typical fungal UPR, but possesses alternative mechanisms to cope with ER stress. This review summarizes how C. glabrata responds to ER stress and discusses the impacts of ER quality control systems on antifungal resistance and virulence.

Comparative genomics of emerging pathogens in the Candida glabrata clade

BACKGROUND

Candida glabrata follows C. albicans as the second or third most prevalent cause of candidemia worldwide. These two pathogenic yeasts are distantly related, C. glabrata being part of the Nakaseomyces, a group more closely related to Saccharomyces cerevisiae. Although C. glabrata was thought to be the only pathogenic Nakaseomyces, two new pathogens have recently been described within this group: C. nivariensis and C. bracarensis. To gain insight into the genomic changes underlying the emergence of virulence, we sequenced the genomes of these two, and three other non-pathogenic Nakaseomyces, and compared them to other sequenced yeasts.

RESULTS

Our results indicate that the two new pathogens are more closely related to the non-pathogenic N. delphensis than to C. glabrata. We uncover duplications and accelerated evolution that specifically affected genes in the lineage preceding the group containing N. delphensis and the three pathogens, which may provide clues to the higher propensity of this group to infect humans. Finally, the number of Epa-like adhesins is specifically enriched in the pathogens, particularly in C. glabrata.

CONCLUSIONS

Remarkably, some features thought to be the result of adaptation of C. glabrata to a pathogenic lifestyle, are present throughout the Nakaseomyces, indicating these are rather ancient adaptations to other environments. Phylogeny suggests that human pathogenesis evolved several times, independently within the clade. The expansion of the EPA gene family in pathogens establishes an evolutionary link between adhesion and virulence phenotypes. Our analyses thus shed light onto the relationships between virulence and the recent genomic changes that occurred within the Nakaseomyces.

SEQUENCE ACCESSION NUMBERS

Nakaseomyces delphensis: CAPT01000001 to CAPT01000179Candida bracarensis: CAPU01000001 to CAPU01000251Candida nivariensis: CAPV01000001 to CAPV01000123Candida castellii: CAPW01000001 to CAPW01000101Nakaseomyces bacillisporus: CAPX01000001 to CAPX01000186.

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

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

Stress response in Candida glabrata: pieces of a fragmented picture

Candida glabrata is closely related to yeast but obviously adapted to human commensalism. Communication with the environment is important to adjust allocation of resources between protection and proliferation in order to adapt to different situations in and outside of the host. Gene transcription regulated by environmental conditions is a major response strategy of simple fungal organisms. Differences to yeast include an extended repertoire of adhesive genes, and high drug, starvation and stress resistance. These properties largely do not originate from novel virulence genes but rather from adaptations of the transcriptional wiring. C. glabrata signaling pathways providing stress protection are adopted to meet conditions possibly encountered in a host-pathogen confrontation. The view on C. glabrata is getting clearer and points to a simple strategy combining resilience and a few adaptations.

Sex at the origin: an Asian population of the rice blast fungus Magnaporthe oryzae reproduces sexually

Sexual reproduction may be cryptic or facultative in fungi and therefore difficult to detect. Magnaporthe oryzae, which causes blast, the most damaging fungal disease of rice, is thought to originate from southeast Asia. It reproduces asexually in all rice-growing regions. Sexual reproduction has been suspected in limited areas of southeast Asia, but has never been demonstrated in contemporary populations. We characterized several M. oryzae populations worldwide both biologically and genetically, to identify candidate populations for sexual reproduction. The sexual cycle of M. oryzae requires two strains of opposite mating types, at least one of which is female-fertile, to come into contact. In one Chinese population, the two mating types were found to be present at similar frequencies and almost all strains were female-fertile. Compatible strains from this population completed the sexual cycle in vitro and produced viable progenies. Genotypic richness and linkage disequilibrium data also supported the existence of sexual reproduction in this population. We resampled this population the following year, and the data obtained confirmed the presence of all the biological and genetic characteristics of sexual reproduction. In particular, a considerable genetic reshuffling of alleles was observed between the 2 years. Computer simulations confirmed that the observed genetic characteristics were unlikely to have arisen in the absence of recombination. We therefore concluded that a contemporary population of M. oryzae, pathogenic on rice, reproduces sexually in natura in southeast Asia. Our findings provide evidence for the loss of sexual reproduction by a fungal plant pathogen outside its centre of origin.

Evolutionary erosion of yeast sex chromosomes by mating-type switching accidents

We investigate yeast sex chromosome evolution by comparing genome sequences from 16 species in the family Saccharomycetaceae, including data from genera Tetrapisispora, Kazachstania, Naumovozyma, and Torulaspora. We show that although most yeast species contain a mating-type (MAT) locus and silent HML and HMR loci structurally analogous to those of Saccharomyces cerevisiae, their detailed organization is highly variable and indicates that the MAT locus is a deletion hotspot. Over evolutionary time, chromosomal genes located immediately beside MAT have continually been deleted, truncated, or transposed to other places in the genome in a process that is gradually shortening the distance between MAT and HML. Each time a gene beside MAT is removed by deletion or transposition, the next gene on the chromosome is brought into proximity with MAT and is in turn put at risk for removal. This process has also continually replaced the triplicated sequence regions, called Z and X, that allow HML and HMR to be used as templates for DNA repair at MAT during mating-type switching. We propose that the deletion and transposition events are caused by evolutionary accidents during mating-type switching, combined with natural selection to keep MAT and HML on the same chromosome. The rate of deletion accelerated greatly after whole-genome duplication, probably because genes were redundant and could be deleted without requiring transposition. We suggest that, despite its mutational cost, switching confers an evolutionary benefit by providing a way for an isolated germinating spore to reform spores if the environment is too poor.

Conserved meiotic machinery in Glomus spp., a putatively ancient asexual fungal lineage

Arbuscular mycorrhizal fungi (AMF) represent an ecologically important and evolutionarily intriguing group of symbionts of land plants, currently thought to have propagated clonally for over 500 Myr. AMF produce multinucleate spores and may exchange nuclei through anastomosis, but meiosis has never been observed in this group. A provocative alternative for their successful and long asexual evolutionary history is that these organisms may have cryptic sex, allowing them to recombine alleles and compensate for deleterious mutations. This is partly supported by reports of recombination among some of their natural populations. We explored this hypothesis by searching for some of the primary tools for a sustainable sexual cycle—the genes whose products are required for proper completion of meiotic recombination in yeast—in the genomes of four AMF and compared them with homologs of representative ascomycete, basidiomycete, chytridiomycete, and zygomycete fungi. Our investigation used molecular and bioinformatic tools to identify homologs of 51 meiotic genes, including seven meiosis-specific genes and other “core meiotic genes” conserved in the genomes of the AMF Glomus diaphanum (MUCL 43196), Glomus irregulare (DAOM-197198), Glomus clarum (DAOM 234281), and Glomus cerebriforme (DAOM 227022). Homology of AMF meiosis-specific genes was verified by phylogenetic analyses with representative fungi, animals (Mus, Hydra), and a choanoflagellate (Monosiga). Together, these results indicate that these supposedly ancient asexual fungi may be capable of undergoing a conventional meiosis; a hypothesis that is consistent with previous reports of recombination within and across some of their populations.

Reinventing heterochromatin in budding yeasts: Sir2 and the origin recognition complex take center stage

The transcriptional silencing of the cryptic mating-type loci in Saccharomyces cerevisiae is one of the best-studied models of repressive heterochromatin. However, this type of heterochromatin, which is mediated by the Sir proteins, has a distinct molecular composition compared to the more ubiquitous type of heterochromatin found in Schizosaccharomyces pombe, other fungi, animals, and plants and characterized by the presence of HP1 (heterochromatin protein 1). This review discusses how the loss of important heterochromatin proteins, including HP1, in the budding yeast lineage presented an evolutionary opportunity for the development and diversification of alternative varieties of heterochromatin, in which the conserved deacetylase Sir2 and the replication protein Orc1 play key roles. In addition, we highlight how this diversification has been facilitated by gene duplications and has contributed to adaptations in lifestyle.

Haploinsufficiency and the sex chromosomes from yeasts to humans

BACKGROUND

Haploinsufficient (HI) genes are those for which a reduction in copy number in a diploid from two to one results in significantly reduced fitness. Haploinsufficiency is increasingly implicated in human disease, and so predicting this phenotype could provide insights into the genetic mechanisms behind many human diseases, including some cancers.

RESULTS

In the present work we show that orthologues of Saccharomyces cerevisiae HI genes are preferentially retained across the kingdom Fungi, and that the HI genes of S. cerevisiae can be used to predict haploinsufficiency in humans. Our HI gene predictions confirm known associations between haploinsufficiency and genetic disease, and predict several further disorders in which the phenotype may be relevant. Haploinsufficiency is also clearly relevant to the gene-dosage imbalances inherent in eukaryotic sex-determination systems. In S. cerevisiae, HI genes are over-represented on chromosome III, the chromosome that determines yeast's mating type. This may be a device to select against the loss of one copy of chromosome III from a diploid. We found that orthologues of S. cerevisiae HI genes are also over-represented on the mating-type chromosomes of other yeasts and filamentous fungi. In animals with heterogametic sex determination, accumulation of HI genes on the sex chromosomes would compromise fitness in both sexes, given X chromosome inactivation in females. We found that orthologues of S. cerevisiae HI genes are significantly under-represented on the X chromosomes of mammals and of Caenorhabditis elegans. There is no X inactivation in Drosophila melanogaster (increased expression of X in the male is used instead) and, in this species, we found no depletion of orthologues to yeast HI genes on the sex chromosomes.

CONCLUSION

A special relationship between HI genes and the sex/mating-type chromosome extends from S. cerevisiae to Homo sapiens, with the microbe being a useful model for species throughout the evolutionary range. Furthermore, haploinsufficiency in yeast can predict the phenotype in higher organisms.

From Saccharomyces cerevisiae to Candida glabratain a few easy steps: important adaptations for an opportunistic pathogen

The opportunistic human fungal pathogen Candida glabrata is closely related to Saccharomyces cerevisiae, yet it has evolved to survive within mammalian hosts. Which traits help C. glabrata to adapt to this different environment? Which specific responses are crucial for its survival in the host? The main differences seem to include an extended repertoire of adhesin genes, high drug resistance, an enhanced ability to sustain prolonged starvation and adaptations of the transcriptional wiring of key stress response genes. Here, we discuss the properties of C. glabrata with a focus on the differences to related fungi.

Fungal sex and pathogenesis

Human fungal pathogens are associated with diseases ranging from dandruff and skin colonization to invasive bloodstream infections. The major human pathogens belong to the Candida, Aspergillus, and Cryptococcus clades, and infections have high and increasing morbidity and mortality. Many human fungal pathogens were originally assumed to be asexual. However, recent advances in genome sequencing, which revealed that many species have retained the genes required for the sexual machinery, have dramatically influenced our understanding of the biology of these organisms. Predictions of a rare or cryptic sexual cycle have been supported experimentally for some species. Here, I examine the evidence that human pathogens reproduce sexually. The evolution of the mating-type locus in ascomycetes (including Candida and Aspergillus species) and basidiomycetes (Malassezia and Cryptococcus) is discussed. I provide an overview of how sex is suppressed in different species and discuss the potential associations with pathogenesis.

Organization and evolutionary trajectory of the mating type (MAT) locus in dermatophyte and dimorphic fungal pathogens

Sexual reproduction in fungi is governed by a specialized genomic region, the mating type (MAT) locus, whose gene identity, organization, and complexity are diverse. We identified the MAT locus of five dermatophyte fungal pathogens (Microsporum gypseum, Microsporum canis, Trichophyton equinum, Trichophyton rubrum, and Trichophyton tonsurans) and a dimorphic fungus, Paracoccidioides brasiliensis, and performed phylogenetic analyses. The identified MAT locus idiomorphs of M. gypseum control cell type identity in mating assays, and recombinant progeny were produced. Virulence tests in Galleria mellonella larvae suggest the two mating types of M. gypseum may have equivalent virulence. Synteny analysis revealed common features of the MAT locus shared among these five dermatophytes: namely, a small size ( approximately 3 kb) and a novel gene arrangement. The SLA2, COX13, and APN2 genes, which flank the MAT locus in other Ascomycota are instead linked on one side of the dermatophyte MAT locus. In addition, the transcriptional orientations of the APN2 and COX13 genes are reversed compared to the dimorphic fungi Histoplasma capsulatum, Coccidioides immitis, and Coccidioides posadasii. A putative transposable element, pogo, was found to have inserted in the MAT1-2 idiomorph of one P. brasiliensis strain but not others. In conclusion, the evolution of the MAT locus of the dermatophytes and dimorphic fungi from the last common ancestor has been punctuated by both gene acquisition and expansion, and asymmetric gene loss. These studies further support a foundation to develop molecular and genetic tools for dermatophyte and dimorphic human fungal pathogens.

Uneven distribution of mating types among genotypes of Candida glabrata isolates from clinical samples

In order to shed light on its basic biology, we initiated a population genetic analysis of Candida glabrata, an emerging pathogenic yeast with no sexual stage yet recognized. A worldwide collection of clinical strains was subjected to analysis using variable number of tandem repeats (VNTR) at nine loci. The clustering of strains obtained with this method was congruent with that obtained using sequence polymorphism of the NMT1 gene, a locus previously proposed for lineage assignment. Linkage disequilibrium supported the hypothesis of a mainly clonal reproduction. No heterozygous diploid genotype was found. Minimum-spanning tree analysis of VNTR data revealed clonal expansions and associated genotypic diversification. Mating type analysis revealed that 80% of the strains examined are MATa and 20% MATalpha and that the two alleles are not evenly distributed. The MATa genotype dominated within large clonal groups that contained only one or a few MATalpha types. In contrast, two groups were dominated by MATalpha strains. Our data are consistent with rare independent mating type switching events occurring preferentially from type a to alpha, although the alternative possibility of selection favoring type a isolates cannot be excluded.

Phylogenetic and functional assessment of orthologs inference projects and methods

Accurate genome-wide identification of orthologs is a central problem in comparative genomics, a fact reflected by the numerous orthology identification projects developed in recent years. However, only a few reports have compared their accuracy, and indeed, several recent efforts have not yet been systematically evaluated. Furthermore, orthology is typically only assessed in terms of function conservation, despite the phylogeny-based original definition of Fitch. We collected and mapped the results of nine leading orthology projects and methods (COG, KOG, Inparanoid, OrthoMCL, Ensembl Compara, Homologene, RoundUp, EggNOG, and OMA) and two standard methods (bidirectional best-hit and reciprocal smallest distance). We systematically compared their predictions with respect to both phylogeny and function, using six different tests. This required the mapping of millions of sequences, the handling of hundreds of millions of predicted pairs of orthologs, and the computation of tens of thousands of trees. In phylogenetic analysis or in functional analysis where high specificity is required, we find that OMA and Homologene perform best. At lower functional specificity but higher coverage level, OrthoMCL outperforms Ensembl Compara, and to a lesser extent Inparanoid. Lastly, the large coverage of the recent EggNOG can be of interest to build broad functional grouping, but the method is not specific enough for phylogenetic or detailed function analyses. In terms of general methodology, we observe that the more sophisticated tree reconstruction/reconciliation approach of Ensembl Compara was at times outperformed by pairwise comparison approaches, even in phylogenetic tests. Furthermore, we show that standard bidirectional best-hit often outperforms projects with more complex algorithms. First, the present study provides guidance for the broad community of orthology data users as to which database best suits their needs. Second, it introduces new methodology to verify orthology. And third, it sets performance standards for current and future approaches.

The asexual yeast Candida glabrata maintains distinct a and alpha haploid mating types

The genome of the type strain of Candida glabrata (CBS138, ATCC 2001) contains homologs of most of the genes involved in mating in Saccharomyces cerevisiae, starting with the mating pheromone and receptor genes. Only haploid cells are ever isolated, but C. glabrata strains of both mating types are commonly found, the type strain being MAT alpha and most other strains, such as BG2, being MATa. No sexual cycle has been documented for this species. In order to understand which steps of the mating pathway are defective, we have analyzed the expression of homologs of some of the key genes involved as well as the production of mating pheromones and the organism's sensitivity to artificial pheromones. We show that cells of opposite mating types express both pheromone receptor genes and are insensitive to pheromones. Nonetheless, cells maintain specificity through regulation of the alpha1 and alpha2 genes and, more surprisingly, through differential splicing of the a1 transcript.

Multiple-locus variable-number tandem-repeat analysis for rapid typing of Candida glabrata

A multiple-locus variable-number tandem-repeat analysis (MLVA) using six microsatellite markers was assessed in 127 Candida glabrata isolates. Thirty-seven different genotypes, stable both in vitro and in vivo, were observed. The highest discriminatory power (D = 0.902) was reached by using only four markers. MLVA seems to be relevant for C. glabrata typing.