CandidaDB: a multi-genome database for Candida species and related Saccharomycotina

CandidaDB (http://genodb.pasteur.fr/CandidaDB) was established in 2002 to provide the first genomic database for the human fungal pathogen Candida albicans. The availability of an increasing number of fully or partially completed genome sequences of related fungal species has opened the path for comparative genomics and prompted us to migrate CandidaDB into a multi-genome database. The new version of CandidaDB houses the latest versions of the genomes of C. albicans strains SC5314 and WO-1 along with six genome sequences from species closely related to C. albicans that all belong to the CTG clade of Saccharomycotina—Candida tropicalis, Candida (Clavispora) lusitaniae, Candida (Pichia) guillermondii, Lodderomyces elongisporus, Debaryomyces hansenii, Pichia stipitis—and the reference Saccharomyces cerevisiae genome. CandidaDB includes sequences coding for 54 170 proteins with annotations collected from other databases, enriched with illustrations of structural features and functional domains and data of comparative analyses. In order to take advantage of the integration of multiple genomes in a unique database, new tools using pre-calculated or user-defined comparisons have been implemented that allow rapid access to comparative analysis at the genomic scale

Identification and Characterization of Mediators of Fluconazole Tolerance in Candida albicans.

Candida albicans is an important human pathogen and a major concern in intensive care units around the world. C. albicans infections are associated with a high mortality despite the use of antifungal treatments. One of the causes of therapeutic failures is the acquisition of antifungal resistance by mutations in the C. albicans genome. Fluconazole (FLC) is one of the most widely used antifungal and mechanisms of FLC resistance occurring by mutations have been extensively investigated. However, some clinical isolates are known to be able to survive at high FLC concentrations without acquiring resistance mutations, a phenotype known as tolerance. Mechanisms behind FLC tolerance are not well studied, mainly due to the lack of a proper way to identify and quantify tolerance in clinical isolates. We proposed here culture conditions to investigate FLC tolerance as well as an easy and efficient method to identity and quantify tolerance to FLC. The screening of C. albicans strain collections revealed that FLC tolerance is pH- and strain-dependent, suggesting the involvement of multiple mechanisms. Here, we addressed the identification of FLC tolerance mediators in C. albicans by an overexpression strategy focusing on 572 C. albicans genes. This strategy led to the identification of two transcription factors, CRZ1 and GZF3. CRZ1 is a C2H2-type transcription factor that is part of the calcineurin-dependent pathway in C. albicans, while GZF3 is a GATA-type transcription factor of unknown function in C. albicans. Overexpression of each gene resulted in an increase of FLC tolerance, however, only the deletion of CRZ1 in clinical FLC-tolerant strains consistently decreased their FLC tolerance. Transcription profiling of clinical isolates with variable levels of FLC tolerance confirmed a calcineurin-dependent signature in these isolates when exposed to FLC

A Novel Regulator Couples Sporogenesis and Trehalose Biogenesis in Aspergillus nidulans

Trehalose is a compatible osmolyte produced by bacteria, fungi, insects and plants to protect the integrity of cells against various environmental stresses. Spores, the reproductive, survival and infection bodies of fungi require high amounts of trehalose for long-term survival. Here, via a gain-of-function genetic screen, we identify the novel regulator VosA that couples the formation of spores and focal trehalose biogenesis in the model fungus Aspergillus nidulans. The vosA gene is expressed specifically during the formation of both sexual and asexual spores (conidia). Levels of vosA mRNA and protein are high in both types of spore. The deletion of vosA results in the lack of trehalose in spores, a rapid loss of the cytoplasm, organelles and viability of spores, and a dramatic reduction in tolerance of conidia to heat and oxidative stress. Moreover, the absence of vosA causes uncontrolled activation of asexual development, whereas the enhanced expression of vosA blocks sporulation, suggesting that VosA also functions in negative-feedback regulation of sporogenesis. VosA localizes in the nucleus of mature conidia and its C-terminal region contains a potential transcription activation domain, indicating that it may function as a transcription factor primarily controlling the late process of sporulation including trehalose biogenesis. VosA is conserved in most fungi and may define a new fungus-specific transcription factor family

CandidaDB: a genome database for Candida albicans pathogenomics

CandidaDB is a database dedicated to the genome of the most prevalent systemic fungal pathogen of humans, Candida albicans. CandidaDB is based on an annotation of the Stanford Genome Technology Center C.albicans genome sequence data by the European Galar Fungail Consortium. CandidaDB Release 2.0 (June 2004) contains information pertaining to Assembly 19 of the genome of C.albicans strain SC5314. The current release contains 6244 annotated entries corresponding to 130 tRNA genes and 5917 protein-coding genes. For these, it provides tentative functional assignments along with numerous pre-run analyses that can assist the researcher in the evaluation of gene function for the purpose of specific or large-scale analysis. CandidaDB is based on GenoList, a generic relational data schema and a World Wide Web interface that has been adapted to the handling of eukaryotic genomes. The interface allows users to browse easily through genome data and retrieve information. CandidaDB also provides more elaborate tools, such as pattern searching, that are tightly connected to the overall browsing system. As the C.albicans genome is diploid and still incompletely assembled, CandidaDB provides tools to browse the genome by individual supercontigs and to examine information about allelic sequences obtained from complementary contigs. CandidaDB is accessible at http://genolist.pasteur.fr/CandidaDB

Hsp90 governs dispersion and drug resistance of fungal biofilms

Fungal biofilms are a major cause of human mortality and are recalcitrant to most treatments due to intrinsic drug resistance. These complex communities of multiple cell types form on indwelling medical devices and their eradication often requires surgical removal of infected devices. Here we implicate the molecular chaperone Hsp90 as a key regulator of biofilm dispersion and drug resistance. We previously established that in the leading human fungal pathogen, Candida albicans, Hsp90 enables the emergence and maintenance of drug resistance in planktonic conditions by stabilizing the protein phosphatase calcineurin and MAPK Mkc1. Hsp90 also regulates temperature-dependent C. albicans morphogenesis through repression of cAMP-PKA signalling. Here we demonstrate that genetic depletion of Hsp90 reduced C. albicans biofilm growth and maturation in vitro and impaired dispersal of biofilm cells. Further, compromising Hsp90 function in vitro abrogated resistance of C. albicans biofilms to the most widely deployed class of antifungal drugs, the azoles. Depletion of Hsp90 led to reduction of calcineurin and Mkc1 in planktonic but not biofilm conditions, suggesting that Hsp90 regulates drug resistance through different mechanisms in these distinct cellular states. Reduction of Hsp90 levels led to a marked decrease in matrix glucan levels, providing a compelling mechanism through which Hsp90 might regulate biofilm azole resistance. Impairment of Hsp90 function genetically or pharmacologically transformed fluconazole from ineffectual to highly effective in eradicating biofilms in a rat venous catheter infection model. Finally, inhibition of Hsp90 reduced resistance of biofilms of the most lethal mould, Aspergillus fumigatus, to the newest class of antifungals to reach the clinic, the echinocandins. Thus, we establish a novel mechanism regulating biofilm drug resistance and dispersion and that targeting Hsp90 provides a much-needed strategy for improving clinical outcome in the treatment of biofilm infections

Comparing genomic variant identification protocols for Candida auris.

Genomic analyses are widely applied to epidemiological, population genetic and experimental studies of pathogenic fungi. A wide range of methods are employed to carry out these analyses, typically without including controls that gauge the accuracy of variant prediction. The importance of tracking outbreaks at a global scale has raised the urgency of establishing high-accuracy pipelines that generate consistent results between research groups. To evaluate currently employed methods for whole-genome variant detection and elaborate best practices for fungal pathogens, we compared how 14 independent variant calling pipelines performed across 35 Candida auris isolates from 4 distinct clades and evaluated the performance of variant calling, single-nucleotide polymorphism (SNP) counts and phylogenetic inference results. Although these pipelines used different variant callers and filtering criteria, we found high overall agreement of SNPs from each pipeline. This concordance correlated with site quality, as SNPs discovered by a few pipelines tended to show lower mapping quality scores and depth of coverage than those recovered by all pipelines. We observed that the major differences between pipelines were due to variation in read trimming strategies, SNP calling methods and parameters, and downstream filtration criteria. We calculated specificity and sensitivity for each pipeline by aligning three isolates with chromosomal level assemblies and found that the GATK-based pipelines were well balanced between these metrics. Selection of trimming methods had a greater impact on SAMtools-based pipelines than those using GATK. Phylogenetic trees inferred by each pipeline showed high consistency at the clade level, but there was more variability between isolates from a single outbreak, with pipelines that used more stringent cutoffs having lower resolution. This project generated two truth datasets useful for routine benchmarking of C. auris variant calling, a consensus VCF of genotypes discovered by 10 or more pipelines across these 35 diverse isolates and variants for 2 samples identified from whole-genome alignments. This study provides a foundation for evaluating SNP calling pipelines and developing best practices for future fungal genomic studies

A Human-Curated Annotation of the Candida albicans Genome

Recent sequencing and assembly of the genome for the fungal pathogen Candida albicans used simple automated procedures for the identification of putative genes. We have reviewed the entire assembly, both by hand and with additional bioinformatic resources, to accurately map and describe 6,354 genes and to identify 246 genes whose original database entries contained sequencing errors (or possibly mutations) that affect their reading frame. Comparison with other fungal genomes permitted the identification of numerous fungus-specific genes that might be targeted for antifungal therapy. We also observed that, compared to other fungi, the protein-coding sequences in the C. albicans genome are especially rich in short sequence repeats. Finally, our improved annotation permitted a detailed analysis of several multigene families, and comparative genomic studies showed that C. albicans has a far greater catabolic range, encoding respiratory Complex 1, several novel oxidoreductases and ketone body degrading enzymes, malonyl-CoA and enoyl-CoA carriers, several novel amino acid degrading enzymes, a variety of secreted catabolic lipases and proteases, and numerous transporters to assimilate the resulting nutrients. The results of these efforts will ensure that the Candida research community has uniform and comprehensive genomic information for medical research as well as for future diagnostic and therapeutic applications

Pleiotropic effects of the vacuolar ABC transporter MLT1 of Candida albicans on cell function and virulence.

Among the several mechanisms that contribute to MDR (multidrug resistance), the overexpression of drug-efflux pumps belonging to the ABC (ATP-binding cassette) superfamily is the most frequent cause of resistance to antifungal agents. The multidrug transporter proteins Cdr1p and Cdr2p of the ABCG subfamily are major players in the development of MDR in Candida albicans Because several genes coding for ABC proteins exist in the genome of C. albicans, but only Cdr1p and Cdr2p have established roles in MDR, it is implicit that the other members of the ABC family also have alternative physiological roles. The present study focuses on an ABC transporter of C. albicans, Mlt1p, which is localized in the vacuolar membrane and specifically transports PC (phosphatidylcholine) into the vacuolar lumen. Transcriptional profiling of the mlt1∆/∆ mutant revealed a down-regulation of the genes involved in endocytosis, oxidoreductase activity, virulence and hyphal development. High-throughput MS-based lipidome analysis revealed that the Mlt1p levels affect lipid homoeostasis and thus lead to a plethora of physiological perturbations. These include a delay in endocytosis, inefficient sequestering of reactive oxygen species (ROS), defects in hyphal development and attenuated virulence. The present study is an emerging example where new and unconventional roles of an ABC transporter are being identified

Fusion and Fission of Genes Define a Metric between Fungal Genomes

Gene fusion and fission events are key mechanisms in the evolution of gene architecture, whose effects are visible in protein architecture when they occur in coding sequences. Until now, the detection of fusion and fission events has been performed at the level of protein sequences with a post facto removal of supernumerary links due to paralogy, and often did not include looking for events defined only in single genomes. We propose a method for the detection of these events, defined on groups of paralogs to compensate for the gene redundancy of eukaryotic genomes, and apply it to the proteomes of 12 fungal species. We collected an inventory of 1,680 elementary fusion and fission events. In half the cases, both composite and element genes are found in the same species. Per-species counts of events correlate with the species genome size, suggesting a random mechanism of occurrence. Some biological functions of the genes involved in fusion and fission events are slightly over- or under-represented. As already noted in previous studies, the genes involved in an event tend to belong to the same functional category. We inferred the position of each event in the evolution tree of the 12 fungal species. The event localization counts for all the segments of the tree provide a metric that depicts the “recombinational” phylogeny among fungi. A possible interpretation of this metric as distance in adaptation space is proposed

The gene structure and expression of human ABHD1: overlapping polyadenylation signal sequence with Sec12

BACKGROUND: Overlapping sense/antisense genes orientated in a tail-to-tail manner, often involving only the 3'UTRs, form the majority of gene pairs in mammalian genomes and can lead to the formation of double-stranded RNA that triggers the destruction of homologous mRNAs. Overlapping polyadenylation signal sequences have not been described previously. RESULTS: An instance of gene overlap has been found involving a shared single functional polyadenylation site. The genes involved are the human alpha/beta hydrolase domain containing gene 1 (ABHD1) and Sec12 genes. The nine exon human ABHD1 gene is located on chromosome 2p23.3 and encodes a 405-residue protein containing a catalytic triad analogous to that present in serine proteases. The Sec12 protein promotes efficient guanine nucleotide exchange on the Sar1 GTPase in the ER. Their sequences overlap for 42 bp in the 3'UTR in an antisense manner. Analysis by 3' RACE identified a single functional polyadenylation site, ATTAAA, within the 3'UTR of ABHD1 and a single polyadenylation signal, AATAAA, within the 3'UTR of Sec12. These polyadenylation signals overlap, sharing three bp. They are also conserved in mouse and rat. ABHD1 was expressed in all tissues and cells examined, but levels of ABHD1 varied greatly, being high in skeletal muscle and testis and low in spleen and fibroblasts. CONCLUSIONS: Mammalian ABHD1 and Sec12 genes contain a conserved 42 bp overlap in their 3'UTR, and share a conserved TTTATTAAA/TTTAATAAA sequence that serves as a polyadenylation signal for both genes. No inverse correlation between the respective levels of ABHD1 and Sec12 RNA was found to indicate that any RNA interference occurred