Analysis of gene evolution and metabolic pathways using the Candida Gene Order Browser

Biofilm-associated metabolism via ERG251 in Candida albicans

Biofilm formation by the fungal pathogen Candida albicans is the basis for its ability to infect medical devices. The metabolic gene ERG251 has been identified as a target of biofilm transcriptional regulator Efg1, and here we report that ERG251 is required for biofilm formation but not conventional free-living planktonic growth. An erg251Δ/Δ mutation impairs biofilm formation in vitro and in an in vivo catheter infection model. In both in vitro and in vivo biofilm contexts, cell number is reduced and hyphal length is limited. To determine whether the mutant defect is in growth or some other aspect of biofilm development, we examined planktonic cell features in a biofilm-like environment, which was approximated with sealed unshaken cultures. Under those conditions, the erg251Δ/Δ mutation causes defects in growth and hyphal extension. Overexpression in the erg251Δ/Δ mutant of the paralog ERG25, which is normally expressed more weakly than ERG251, partially improves biofilm formation and biofilm hyphal content, as well as growth and hyphal extension in a biofilm-like environment. GC-MS analysis shows that the erg251Δ/Δ mutation causes a defect in ergosterol accumulation when cells are cultivated under biofilm-like conditions, but not under conventional planktonic conditions. Overexpression of ERG25 in the erg251Δ/Δ mutant causes some increase in ergosterol levels. Finally, the hypersensitivity of efg1Δ/Δ mutants to the ergosterol inhibitor fluconazole is reversed by ERG251 overexpression, arguing that reduced ERG251 expression contributes to this efg1Δ/Δ phenotype. Our results indicate that ERG251 is required for biofilm formation because its high expression levels are necessary for ergosterol synthesis in a biofilm-like environment.

A gain-of-function mutation in zinc cluster transcription factor Rob1 drives Candida albicans adaptive growth in the cystic fibrosis lung environment

Candida albicans chronically colonizes the respiratory tract of patients with Cystic Fibrosis (CF). It competes with CF-associated pathogens (e.g. Pseudomonas aeruginosa) and contributes to disease severity. We hypothesize that C. albicans undergoes specific adaptation mechanisms that explain its persistence in the CF lung environment. To identify the underlying genetic and phenotypic determinants, we serially recovered 146 C. albicans clinical isolates over a period of 30 months from the sputum of 25 antifungal-naive CF patients. Multilocus sequence typing analyses revealed that most patients were individually colonized with genetically close strains, facilitating comparative analyses between serial isolates. We strikingly observed differential ability to filament and form monospecies and dual-species biofilms with P. aeruginosa among 18 serial isolates sharing the same diploid sequence type, recovered within one year from a pediatric patient. Whole genome sequencing revealed that their genomes were highly heterozygous and similar to each other, displaying a highly clonal subpopulation structure. Data mining identified 34 non-synonymous heterozygous SNPs in 19 open reading frames differentiating the hyperfilamentous and strong biofilm-former strains from the remaining isolates. Among these, we detected a glycine-to-glutamate substitution at position 299 (G299E) in the deduced amino acid sequence of the zinc cluster transcription factor ROB1 (ROB1G299E), encoding a major regulator of filamentous growth and biofilm formation. Introduction of the G299E heterozygous mutation in a co-isolated weak biofilm-former CF strain was sufficient to confer hyperfilamentous growth, increased expression of hyphal-specific genes, increased monospecies biofilm formation and increased survival in dual-species biofilms formed with P. aeruginosa, indicating that ROB1G299E is a gain-of-function mutation. Disruption of ROB1 in a hyperfilamentous isolate carrying the ROB1G299E allele abolished hyperfilamentation and biofilm formation. Our study links a single heterozygous mutation to the ability of C. albicans to better survive during the interaction with other CF-associated microbes and illuminates how adaptive traits emerge in microbial pathogens to persistently colonize and/or infect the CF-patient airways.

Candida auris-macrophage cellular interactions and transcriptional response

The pathogenic yeast Candida auris represents a global threat of the utmost clinical relevance. This emerging fungal species is remarkable in its resistance to commonly used antifungal agents and its persistence in the nosocomial settings. The innate immune system is one the first lines of defense preventing the dissemination of pathogens in the host. C. auris is susceptible to circulating phagocytes, and understanding the molecular details of these interactions may suggest routes to improved therapies. In this work, we examined the interactions of this yeast with macrophages. We found that macrophages avidly phagocytose C. auris; however, intracellular replication is not inhibited, indicating that C. auris resists the killing mechanisms imposed by the phagocyte. Unlike Candida albicans, phagocytosis of C. auris does not induce macrophage lysis. The transcriptional response of C. auris to macrophage phagocytosis is very similar to other members of the CUG clade (C. albicans, C. tropicalis, C. parapsilosis, C. lusitaniae), i.e., downregulation of transcription/translation and upregulation of alternative carbon metabolism pathways, transporters, and induction of oxidative stress response and proteolysis. Gene family expansions are common in this yeast, and we found that many of these genes are induced in response to macrophage co-incubation. Among these, amino acid and oligopeptide transporters, as well as lipases and proteases, are upregulated. Thus, C. auris shares key transcriptional signatures shared with other fungal pathogens and capitalizes on the expansion of gene families coding for potential virulence attributes that allow its survival, persistence, and evasion of the innate immune system.

The root pathogen Aphanomyces euteiches secretes modular proteases in pea apoplast during host infection

To successfully colonize the host, phytopathogens have developed a large repertoire of components to both combat the host plant defense mechanisms and to survive in adverse environmental conditions. Microbial proteases are predicted to be crucial components of these systems. In the present work, we aimed to identify active secreted proteases from the oomycete Aphanomyces euteiches, which causes root rot diseases on legumes. Genome mining and expression analysis highlighted an overrepresentation of microbial tandemly repeated proteases, which are upregulated during host infection. Activity Based Protein Profiling and mass spectrometry (ABPP-MS) on apoplastic fluids isolated from pea roots infected by the pathogen led to the identification of 35 active extracellular microbial proteases, which represents around 30% of the genes expressed encoding serine and cysteine proteases during infection. Notably, eight of the detected active secreted proteases carry an additional C-terminal domain. This study reveals novel active modular extracellular eukaryotic proteases as potential pathogenicity factors in Aphanomyces genus.

Recognition of galactose by a scaffold protein recruits a transcriptional activator for the GAL regulon induction in Candida albicans

The GAL pathway of yeasts has long served as a model system for understanding of how regulatory mode of eukaryotic metabolic pathways evolves. While Gal4 mode has been well-characterized in Saccharomycetaceae clade, little is known about the regulation of the GAL pathway in other yeasts. Here, we find that Rep1, a Ndt80-like family transcription factor, serves as a galactose sensor in the commensal-pathogenic fungus Candida albicans. It is presented at the GAL gene promoters independent of the presence of galactose. Rep1 recognizes galactose via a direct physical interaction. The net result of this interaction is the recruitment of a transcriptional activator Cga1 (Candida galactose gene activator, orf19.4959) and transcription of the GAL genes proceeds. Rep1 and Cga1 are conserved across the CTG species. Rep1 itself does not possess transcriptional activity. Instead, it provides a scaffold to recruit different factors for transcriptional regulation. Rep1-Cga1 mode of regulation represents a new example of network rewiring in fungi, which provides insight into how C. albicans evolves transcriptional programs to colonize diverse host niches.

An orthologous gene coevolution network provides insight into eukaryotic cellular and genomic structure and function

The evolutionary rates of functionally related genes often covary. We present a gene coevolution network inferred from examining nearly 3 million orthologous gene pairs from 332 budding yeast species spanning ~400 million years of evolution. Network modules provide insight into cellular and genomic structure and function. Examination of the phenotypic impact of network perturbation using deletion mutant data from the baker's yeast Saccharomyces cerevisiae, which were obtained from previously published studies, suggests that fitness in diverse environments is affected by orthologous gene neighborhood and connectivity. Mapping the network onto the chromosomes of S. cerevisiae and Candida albicans revealed that coevolving orthologous genes are not physically clustered in either species; rather, they are often located on different chromosomes or far apart on the same chromosome. The coevolution network captures the hierarchy of cellular structure and function, provides a roadmap for genotype-to-phenotype discovery, and portrays the genome as a linked ensemble of genes.

How to Use the Candida Genome Database

The Candida Genome Database provides access to biological information about genes and proteins of several medically important Candida species. The website is organized into easily navigable pages that enable data retrieval and analysis. This chapter shows how to explore the CGD Home page and Locus Summary pages, which are the main access points to the database. It also provides a description of how to use the GO analysis tools, GO Term Finder, and GO Slim Mapper and how to browse large-scale datasets using the JBrowse genome browser. Finally, it shows how to search and retrieve data for user-defined sets of genes using the Advanced Search and Batch Download tools.

Non-conventional yeasts for food and additives production in a circular economy perspective

Yeast species have been spontaneously participating in food production for millennia, but the scope of applications was greatly expanded since their key role in beer and wine fermentations was clearly acknowledged. The workhorse for industry and scientific research has always been Saccharomyces cerevisiae. It occupies the largest share of the dynamic yeast market, that could further increase thanks to the better exploitation of other yeast species. Food-related 'non-conventional' yeasts (NCY) represent a treasure trove for bioprospecting, with their huge untapped potential related to a great diversity of metabolic capabilities linked to niche adaptations. They are at the crossroad of bioprocesses and biorefineries, characterized by low biosafety risk and produce food and additives, being also able to contribute to production of building blocks and energy recovered from the generated waste and by-products. Considering that the usual pattern for bioprocess development focuses on single strains or species, in this review we suggest that bioprospecting at the genus level could be very promising. Candida, Starmerella, Kluyveromyces and Lachancea were briefly reviewed as case studies, showing that a taxonomy- and genome-based rationale could open multiple possibilities to unlock the biotechnological potential of NCY bioresources.

Lineage-specific selection and the evolution of virulence in the Candida clade

Of the many microbial species on earth, only a small number are able to thrive in humans and cause disease. Comparison of closely related pathogenic and nonpathogenic species can therefore be useful in identifying key features that contribute to virulence. We created interspecies hybrids between Candida albicans, a prevalent fungal pathogen of humans, and Candida dubliniensis, a close, but much less pathogenic, relative. By comparing genome-wide expression differences between the two genomes in the same cell, we surmised that since the two species diverged from a common ancestor, natural selection has acted upon the expression level of an ancient metabolic pathway, illustrating that pathogenicity traits can arise over evolutionary timescales through small expression changes in deeply conserved proteins.

Candida albicans is the most common cause of systemic fungal infections in humans and is considerably more virulent than its closest known relative, Candida dubliniensis. To investigate this difference, we constructed interspecies hybrids and quantified mRNA levels produced from each genome in the hybrid. This approach systematically identified expression differences in orthologous genes arising from cis-regulatory sequence changes that accumulated since the two species last shared a common ancestor, some 10 million y ago. We documented many orthologous gene-expression differences between the two species, and we pursued one striking observation: All 15 genes coding for the enzymes of glycolysis showed higher expression from the C. albicans genome than the C. dubliniensis genome in the interspecies hybrid. This pattern requires evolutionary changes to have occurred at each gene; the fact that they all act in the same direction strongly indicates lineage-specific natural selection as the underlying cause. To test whether these expression differences contribute to virulence, we created a C. dubliniensis strain in which all 15 glycolysis genes were produced at modestly elevated levels and found that this strain had significantly increased virulence in the standard mouse model of systemic infection. These results indicate that small expression differences across a deeply conserved set of metabolism enzymes can play a significant role in the evolution of virulence in fungal pathogens.

Population genomics of the pathogenic yeast Candida tropicalis identifies hybrid isolates in environmental samples

Candida tropicalis is a human pathogen that primarily infects the immunocompromised. Whereas the genome of one isolate, C. tropicalis MYA-3404, was originally sequenced in 2009, there have been no large-scale, multi-isolate studies of the genetic and phenotypic diversity of this species. Here, we used whole genome sequencing and phenotyping to characterize 77 isolates of C. tropicalis from clinical and environmental sources from a variety of locations. We show that most C. tropicalis isolates are diploids with approximately 2-6 heterozygous variants per kilobase. The genomes are relatively stable, with few aneuploidies. However, we identified one highly homozygous isolate and six isolates of C. tropicalis with much higher heterozygosity levels ranging from 36-49 heterozygous variants per kilobase. Our analyses show that the heterozygous isolates represent two different hybrid lineages, where the hybrids share one parent (A) with most other C. tropicalis isolates, but the second parent (B or C) differs by at least 4% at the genome level. Four of the sequenced isolates descend from an AB hybridization, and two from an AC hybridization. The hybrids are MTLa/α heterozygotes. Hybridization, or mating, between different parents is therefore common in the evolutionary history of C. tropicalis. The new hybrids were predominantly found in environmental niches, including from soil. Hybridization is therefore unlikely to be associated with virulence. In addition, we used genotype-phenotype correlation and CRISPR-Cas9 editing to identify a genome variant that results in the inability of one isolate to utilize certain branched-chain amino acids as a sole nitrogen source.

Transcription Profiles Associated with Inducible Adhesion in Candida parapsilosis

Candida parapsilosis has emerged as a frequent cause of invasive candidiasis with increasing evidence of unique biological features relative to C. albicans As it adapts to conditions within a mammalian host, rapid changes in gene expression are necessary to facilitate colonization and persistence in this environment. Adhesion of the organism to biological surfaces is a key first step in this process and is the focus of this study. Building on previous observations showing the importance of a member of the ALS gene family in C. parapsilosis adhesion, three clinical isolates were cultured under two conditions that mimic the mammalian host and promote adhesion, incubation at 37°C in tissue culture medium 199 or in human plasma. Transcriptional profiles using RNA-seq were obtained in these adhesion-inducing conditions and compared to profiles following growth in yeast media that suppress adhesion to identify gene expression profiles associated with adhesion. Overall gene expression profiles among the three strains were similar in both adhesion-inducing conditions and distinct from adhesion-suppressing conditions. Pairwise analysis among the three growth conditions identified 133 genes that were differentially expressed at a cutoff of ±4-fold, with the most upregulated genes significantly enriched in iron acquisition and transmembrane transport, while the most downregulated genes were enriched in oxidation-reduction processes. Gene family enrichment analysis identified gene families with diverse functions that may have an important role in this important step for colonization and disease.IMPORTANCE Invasive Candida infections are frequent complications of the immunocompromised and are associated with substantive morbidity and mortality. Although C. albicans is the best-studied species, emerging infections by non-albicans Candida species have led to increased efforts to understand aspects of their pathogenesis that are unique from C. albicans C. parapsilosis is a frequent cause of invasive infections, particularly among premature infants. Recent efforts have identified important virulence mechanisms that have features distinct from C. albicans C. parapsilosis can exist outside a host environment and therefore requires rapid modifications when it encounters a mammalian host to prevent its clearance. An important first step in the process is adhesion to host surfaces. This work takes a global, nonbiased approach to investigate broad changes in gene expression that accompany efficient adhesion. As such, biological pathways and individual protein targets are identified that may be amenable to manipulation to reduce colonization and disease from this organism.

Systematic Analysis of Functionally Related Gene Clusters in the Opportunistic Pathogen, Candida albicans

The proper balance of gene expression is essential for cellular health, organismal development, and maintaining homeostasis. In response to complex internal and external signals, the cell needs to modulate gene expression to maintain proteostasis and establish cellular identity within its niche. On a genome level, single-celled prokaryotic microbes display clustering of co-expressed genes that are regulated as a polycistronic RNA. This phenomenon is largely absent from eukaryotic microbes, although there is extensive clustering of co-expressed genes as functional pairs spread throughout the genome in Saccharomyces cerevisiae. While initial analysis demonstrated conservation of clustering in divergent fungal lineages, a comprehensive analysis has yet to be performed. Here we report on the prevalence, conservation, and significance of the functional clustering of co-regulated genes within the opportunistic human pathogen, Candida albicans. Our analysis reveals that there is extensive clustering within this organism-although the identity of the gene pairs is unique compared with those found in S. cerevisiae-indicating that this genomic arrangement evolved after these microbes diverged evolutionarily, rather than being the result of an ancestral arrangement. We report a clustered arrangement in gene families that participate in diverse molecular functions and are not the result of a divergent orientation with a shared promoter. This arrangement coordinates the transcription of the clustered genes to their neighboring genes, with the clusters congregating to genomic loci that are conducive to transcriptional regulation at a distance.

Pursuing Advances in DNA Sequencing Technology to Solve a Complex Genomic Jigsaw Puzzle: The Agglutinin-Like Sequence (ALS) Genes of Candida tropicalis

The agglutinin-like sequence (ALS) gene family encodes cell-surface adhesins that interact with host and abiotic surfaces, promoting colonization by opportunistic fungal pathogens such as Candida tropicalis. Studies of Als protein contribution to C. tropicalis adhesion would benefit from an accurate catalog of ALS gene sequences as well as insight into relative gene expression levels. Even in the genomics era, this information has been elusive: genome assemblies are often broken within ALS genes because of their extensive regions of highly conserved, repeated DNA sequences and because there are many similar ALS genes at different chromosomal locations. Here, we describe the benefit of long-read DNA sequencing technology to facilitate characterization of C. tropicalis ALS loci. Thirteen ALS loci in C. tropicalis strain MYA-3404 were deduced from a genome assembly constructed from Illumina MiSeq and Oxford Nanopore MinION data. Although the MinION data were valuable, PCR amplification and Sanger sequencing of ALS loci were still required to complete and verify the gene sequences. Each predicted Als protein featured an N-terminal binding domain, a central domain of tandemly repeated sequences, and a C-terminal domain rich in Ser and Thr. The presence of a secretory signal peptide and consensus sequence for addition of a glycosylphosphatidylinositol (GPI) anchor was consistent with predicted protein localization to the cell surface. TaqMan assays were designed to recognize each ALS gene, as well as both alleles at the divergent CtrALS3882 locus. C. tropicalis cells grown in five different in vitro conditions showed differential expression of various ALS genes. To place the C. tropicalis data into a larger context, TaqMan assays were also designed and validated for analysis of ALS gene expression in Candida albicans and Candida dubliniensis. These comparisons identified the subset of highly expressed C. tropicalis ALS genes that were predicted to encode proteins with the most abundant cell-surface presence, prioritizing them for subsequent functional analysis. Data presented here provide a solid foundation for future experimentation to deduce ALS family contributions to C. tropicalis adhesion and pathogenesis.

Polymorphic centromere locations in the pathogenic yeast Candida parapsilosis

Centromeres pose an evolutionary paradox: strongly conserved in function but rapidly changing in sequence and structure. However, in the absence of damage, centromere locations are usually conserved within a species. We report here that isolates of the pathogenic yeast species Candida parapsilosis show within-species polymorphism for the location of centromeres on two of its eight chromosomes. Its old centromeres have an inverted-repeat (IR) structure, whereas its new centromeres have no obvious structural features but are located within 30 kb of the old site. Centromeres can therefore move naturally from one chromosomal site to another, apparently spontaneously and in the absence of any significant changes in DNA sequence. Our observations are consistent with a model in which all centromeres are genetically determined, such as by the presence of short or long IRs or by the ability to form cruciforms. We also find that centromeres have been hotspots for genomic rearrangements in the C. parapsilosis clade.

Fungal biotin homeostasis is essential for immune evasion after macrophage phagocytosis and virulence

Biotin is an important cofactor for multiple enzymes in central metabolic processes. While many bacteria and most fungi are able to synthesise biotin de novo, Candida spp. are auxotrophic for this vitamin and thus require efficient uptake systems to facilitate biotin acquisition during infection. Here we show that Candida glabrata and Candida albicans use a largely conserved system for biotin uptake and regulation, consisting of the high-affinity biotin transporter Vht1 and the transcription factor Vhr1. Both species induce expression of biotin-metabolic genes upon in vitro biotin depletion and following phagocytosis by macrophages, indicating low biotin levels in the Candida-containing phagosome. In line with this, we observed reduced intracellular proliferation of both Candida cells pre-starved of biotin and deletion mutants lacking VHR1 or VHT1 genes. VHT1 was essential for the full virulence of C. albicans during systemic mouse infections, and the lack of VHT1 led to reduced fungal burden in C. glabrata-infected brains and C. albicans-infected brains and kidneys. Together, our data suggest a critical role of Vht1-mediated biotin acquisition for C. glabrata and C. albicans during intracellular growth in macrophages and systemic infections.

N-Acetylglucosamine Regulates Morphogenesis and Virulence Pathways in Fungi

N-acetylglucosamine (GlcNAc) is being increasingly recognized for its ability to stimulate cell signaling. This amino sugar is best known as a component of cell wall peptidoglycan in bacteria, cell wall chitin in fungi and parasites, exoskeletons of arthropods, and the extracellular matrix of animal cells. In addition to these structural roles, GlcNAc is now known to stimulate morphological and stress responses in a wide range of organisms. In fungi, the model organisms Saccharomyces cerevisiae and Schizosaccharomyces pombe lack the ability to respond to GlcNAc or catabolize it, so studies with the human pathogen Candida albicans have been providing new insights into the ability of GlcNAc to stimulate cellular responses. GlcNAc potently induces C. albicans to transition from budding to filamentous hyphal growth. It also promotes an epigenetic switch from White to Opaque cells, which differ in morphology, metabolism, and virulence properties. These studies have led to new discoveries, such as the identification of the first eukaryotic GlcNAc transporter. Other results have shown that GlcNAc can induce signaling in C. albicans in two ways. One is to act as a signaling molecule independent of its catabolism, and the other is that its catabolism can cause the alkalinization of the extracellular environment, which provides an additional stimulus to form hyphae. GlcNAc also induces the expression of virulence genes in the C. albicans, indicating it can influence pathogenesis. Therefore, this review will describe the recent advances in understanding the role of GlcNAc signaling pathways in regulating C. albicans morphogenesis and virulence.

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.

Correlating Genotype and Phenotype in the Asexual Yeast Candida orthopsilosis Implicates ZCF29 in Sensitivity to Caffeine

Candida orthopsilosis is diploid asexual yeast that causes human disease. Most C. orthopsilosis isolates arose from at least four separate hybridizations between related, but not identical, parents. Here, we used population genomics data to correlate genotypic and phenotypic variation in 28 C. orthopsilosis isolates. We used cosine similarity scores to identify 65 variants with potential high-impact (deleterious effects) that correlated with specific phenotypes. Of these, 19 were Single Nucleotide Polymorphisms (SNPs) that changed stop or start codons, or splice sites. One variant resulted in a premature stop codon in both alleles of the gene ZCF29 in C. orthopsilosis isolate 185, which correlated with sensitivity to nystatin and caffeine. We used CRISPR-Cas9 editing to introduce this polymorphism into two resistant C. orthopsilosis isolates. Introducing the stop codon resulted in sensitivity to caffeine and to ketoconazole, but not to nystatin. Our analysis shows that it is possible to associate genomic variants with phenotype in asexual Candida species, but that only a small amount of genomic variation can be easily explored.

Recent trends in molecular diagnostics of yeast infections: from PCR to NGS

The incidence of opportunistic yeast infections in humans has been increasing over recent years. These infections are difficult to treat and diagnose, in part due to the large number and broad diversity of species that can underlie the infection. In addition, resistance to one or several antifungal drugs in infecting strains is increasingly being reported, severely limiting therapeutic options and showcasing the need for rapid detection of the infecting agent and its drug susceptibility profile. Current methods for species and resistance identification lack satisfactory sensitivity and specificity, and often require prior culturing of the infecting agent, which delays diagnosis. Recently developed high-throughput technologies such as next generation sequencing or proteomics are opening completely new avenues for more sensitive, accurate and fast diagnosis of yeast pathogens. These approaches are the focus of intensive research, but translation into the clinics requires overcoming important challenges. In this review, we provide an overview of existing and recently emerged approaches that can be used in the identification of yeast pathogens and their drug resistance profiles. Throughout the text we highlight the advantages and disadvantages of each methodology and discuss the most promising developments in their path from bench to bedside.

Conservation of location of several specific inhibitory codon pairs in the Saccharomyces sensu stricto yeasts reveals translational selection

Synonymous codons provide redundancy in the genetic code that influences translation rates in many organisms, in which overall codon use is driven by selection for optimal codons. It is unresolved if or to what extent translational selection drives use of suboptimal codons or codon pairs. In Saccharomyces cerevisiae, 17 specific inhibitory codon pairs, each comprised of adjacent suboptimal codons, inhibit translation efficiency in a manner distinct from their constituent codons, and many are translated slowly in native genes. We show here that selection operates within Saccharomyces sensu stricto yeasts to conserve nine of these codon pairs at defined positions in genes. Conservation of these inhibitory codon pairs is significantly greater than expected, relative to conservation of their constituent codons, with seven pairs more highly conserved than any other synonymous pair. Conservation is strongly correlated with slow translation of the pairs. Conservation of suboptimal codon pairs extends to two related Candida species, fungi that diverged from Saccharomyces ∼270 million years ago, with an enrichment for codons decoded by I•A and U•G wobble in both Candida and Saccharomyces. Thus, conservation of inhibitory codon pairs strongly implies selection for slow translation at particular gene locations, executed by suboptimal codon pairs.

Unique, Diverged, and Conserved Mitochondrial Functions Influencing Candida albicans Respiration

Candida albicans is an opportunistic fungal pathogen of major clinical concern. The virulence of this pathogen is intimately intertwined with its metabolism. Mitochondria, which have a central metabolic role, have undergone many lineage-specific adaptations in association with their eukaryotic host. A screen for lineage-specific genes identified seven such genes specific to the CTG clade of fungi, of which C. albicans is a member. Each is required for respiratory growth and is integral to expression of complex I, III, or IV of the electron transport chain. Two genes, NUO3 and NUO4, encode supernumerary subunits of complex I, whereas NUE1 and NUE2 have nonstructural roles in expression of complex I. Similarly, the other three genes have nonstructural roles in expression of complex III (QCE1) or complex IV (COE1 and COE2). In addition to these novel additions, an alternative functional assignment was found for the mitochondrial protein encoded by MNE1 MNE1 was required for complex I expression in C. albicans, whereas the distantly related Saccharomyces cerevisiae ortholog participates in expression of complex III. Phenotypic analysis of deletion mutants showed that fermentative metabolism is unable to support optimal growth rates or yields of C. albicans However, yeast-hypha morphogenesis, an important virulence attribute, did not require respiratory metabolism under hypoxic conditions. The inability to respire also resulted in hypersensitivity to the antifungal fluconazole and in attenuated virulence in a Galleria mellonella infection model. The results show that lineage-specific adaptations have occurred in C. albicans mitochondria and highlight the significance of respiratory metabolism in the pathobiology of C. albicans IMPORTANCE Candida albicans is an opportunistic fungal pathogen of major clinical concern. The virulence of this pathogen is intimately intertwined with its metabolic behavior, and mitochondria have a central role in that metabolism. Mitochondria have undergone many evolutionary changes, which include lineage-specific adaptations in association with their eukaryotic host. Seven lineage-specific genes required for electron transport chain function were identified in the CTG clade of fungi, of which C. albicans is a member. Additionally, examination of several highly diverged orthologs encoding mitochondrial proteins demonstrated functional reassignment for one of these. Deficits imparted by deletion of these genes revealed the critical role of respiration in virulence attributes of the fungus and highlight important evolutionary adaptations in C. albicans metabolism.

Characterization of the Candida orthopsilosis agglutinin-like sequence (ALS) genes

Agglutinin like sequence (Als) cell-wall proteins play a key role in adhesion and virulence of Candida species. Compared to the well-characterized Candida albicans ALS genes, little is known about ALS genes in the Candida parapsilosis species complex. Three incomplete ALS genes were identified in the genome sequence for Candida orthopsilosis strain 90–125 (GenBank assembly ASM31587v1): CORT0C04210 (named CoALS4210), CORT0C04220 (CoALS4220) and CORT0B00800 (CoALS800). To complete the gene sequences, new data were derived from strain 90–125 using Illumina (short-read) and Oxford Nanopore (long-read) methods. Long-read sequencing analysis confirmed the presence of 3 ALS genes in C. orthopsilosis 90–125 and resolved the gaps located in repetitive regions of CoALS800 and CoALS4220. In the new genome assembly (GenBank PQBP00000000), the CoALS4210 sequence was slightly longer than in the original assembly. C. orthopsilosis Als proteins encoded features well-known in C. albicans Als proteins such as a secretory signal peptide, N-terminal domain with a peptide-binding cavity, amyloid-forming region, repeated sequences, and a C-terminal site for glycosylphosphatidylinositol anchor addition that, in yeast, suggest localization of the proteins in the cell wall. CoAls4210 and CoAls800 lacked the classic C. albicans Als tandem repeats, instead featuring short, imperfect repeats with consensus motifs such as SSSEPP and GSGN. Quantitative RT-PCR showed differential regulation of CoALS genes by growth stage in six genetically diverse C. orthopsilosis clinical isolates, which also exhibited length variation in the ALS alleles, and strain-specific gene expression patterns. Overall, long-read DNA sequencing methodology was instrumental in generating an accurate assembly of CoALS genes, thus revealing their unconventional features and first insights into their allelic variability within C. orthopsilosis clinical isolates.

Identification of an Exceptionally Long Intron in the HAC1 Gene of Candida parapsilosis

The unfolded protein response (UPR) responds to the build-up of misfolded proteins in the endoplasmic reticulum. The UPR has wide-ranging functions from fungal pathogenesis to applications in biotechnology. The UPR is regulated through the splicing of an unconventional intron in the HAC1 gene. This intron has been described in many fungal species and is of variable length. Until now it was believed that some members of the CTG-Ser1 clade such as C. parapsilosis did not contain an intron in HAC1, suggesting that the UPR was regulated in a different manner. Here we demonstrate that HAC1 plays an important role in regulating the UPR in C. parapsilosis. We also identified an unusually long intron (626 bp) in C. parapsilosisHAC1. Further analysis showed that HAC1 orthologs in several species in the CTG-Ser1 clade contain long introns.

The unfolded protein response (UPR) in the endoplasmic reticulum (ER) is well conserved in eukaryotes from metazoa to yeast. The transcription factor HAC1 is a major regulator of the UPR in many eukaryotes. Deleting HAC1 in the yeast Candida parapsilosis rendered cells more sensitive to DTT, a known inducer of the UPR. The deletion strain was also sensitive to Congo red, calcofluor white, and the antifungal drug ketoconazole, indicating that HAC1 has a role in cell wall maintenance. Transcriptomic analysis revealed that treatment of the wild type with DTT resulted in the increased expression of 368 genes. Comparison with mutant cells treated with DTT reveals that expression of 137 of these genes requires HAC1. Enriched GO term analysis includes response to ER stress, cell wall biogenesis and glycosylation. Orthologs of many of these are associated with UPR in Saccharomyces cerevisiae and Candida albicans. Unconventional splicing of an intron from HAC1 mRNA is required to produce a functional transcription factor. The spliced intron varies in length from 19 bases in C. albicans to 379 bases in Candida glabrata, but has not been previously identified in Candida parapsilosis and related species. We used RNA-seq data and in silico analysis to identify the HAC1 intron in 12 species in the CTG-Ser1 clade. We show that the intron has undergone major contractions and expansions in this clade, reaching up to 848 bases. Exposure to DTT induced splicing of the long intron in C. parapsilosisHAC1, inducing the UPR.

IMPORTANCE The unfolded protein response (UPR) responds to the build-up of misfolded proteins in the endoplasmic reticulum. The UPR has wide-ranging functions from fungal pathogenesis to applications in biotechnology. The UPR is regulated through the splicing of an unconventional intron in the HAC1 gene. This intron has been described in many fungal species and is of variable length. Until now it was believed that some members of the CTG-Ser1 clade such as C. parapsilosis did not contain an intron in HAC1, suggesting that the UPR was regulated in a different manner. Here we demonstrate that HAC1 plays an important role in regulating the UPR in C. parapsilosis. We also identified an unusually long intron (626 bp) in C. parapsilosisHAC1. Further analysis showed that HAC1 orthologs in several species in the CTG-Ser1 clade contain long introns.

Role of the Inducible Adhesin CpAls7 in Binding of Candida parapsilosis to the Extracellular Matrix under Fluid Shear

The yeast Candida parapsilosis is an increasingly common cause of systemic fungal infections among immunocompromised individuals, including premature infants. Adhesion to host surfaces is an important step in pathogenesis, but this process has not been extensively studied in this organism. A microfluidics assay was developed to test the ability of C. parapsilosis to adhere to immobilized host extracellular matrix proteins under physiological fluid shear conditions. Growth in mammalian tissue culture medium at 37°C for 3 to 6 h led to the induction of an adhesive phenotype at shear forces of 1 to 5 dynes/cm² in some isolates of C. parapsilosis Glutamic acid, proline, and calcium appeared to be the minimally necessary requirements for increased adhesion in these assays. To determine whether genes homologous to the ALS gene family of C. albicans were important for the adhesive phenotype, the expression levels of 5 homologous C. parapsilosis genes were quantified by using quantitative PCR (qPCR) under conditions leading to increased adhesion. CPAR2_404800 (CpALS7) and CPAR2_404780 showed increased expression levels compared to those in control yeast. The extent of adhesion was variable among different isolates, and linear regression identified the expression of CpALS7 but not CPAR2_404780 as having a strong positive correlation with adhesion. A homozygous CpALS7 deletion strain was deficient in adhesion, whereas the expression of CpALS7 in Saccharomyces cerevisiae resulted in increased adhesion. Together, these data provide strong evidence that CpAls7 aids in the adherence of C. parapsilosis to the extracellular matrix under shear forces and support its previously reported role in virulence.

Dal81 Regulates Expression of Arginine Metabolism Genes in Candida parapsilosis

Fungi can use a wide variety of nitrogen sources. In the absence of preferred sources such as ammonium, glutamate, and glutamine, secondary sources, including most other amino acids, are used. Expression of the nitrogen utilization pathways is very strongly controlled at the transcriptional level. Here, we investigated the regulation of nitrogen utilization in the pathogenic yeast Candida parapsilosis. We found that the functions of many regulators are conserved with respect to Saccharomyces cerevisiae and other fungi. For example, the core GATA activators GAT1 and GLN3 have a conserved role in nitrogen catabolite repression (NCR). There is one ortholog of GZF3 and DAL80, which represses expression of genes in preferred nitrogen sources. The regulators PUT3 and UGA3 are required for metabolism of proline and γ-aminobutyric acid (GABA), respectively. However, the role of the Dal81 transcription factor is distinctly different. In S. cerevisiae, Dal81 is a positive regulator of acquisition of nitrogen from GABA, allantoin, urea, and leucine, and it is required for maximal induction of expression of the relevant pathway genes. In C. parapsilosis, induction of GABA genes is independent of Dal81, and deleting DAL81 has no effect on acquisition of nitrogen from GABA or allantoin. Instead, Dal81 represses arginine synthesis during growth under preferred nitrogen conditions. IMPORTANCE Utilization of nitrogen by fungi is controlled by nitrogen catabolite repression (NCR). Expression of many genes is switched off during growth on nonpreferred nitrogen sources. Gene expression is regulated through a combination of activation and repression. Nitrogen regulation has been studied best in the model yeast Saccharomyces cerevisiae. We found that although many nitrogen regulators have a conserved function in Saccharomyces species, some do not. The Dal81 transcriptional regulator has distinctly different functions in S. cerevisiae and C. parapsilosis. In the former, it regulates utilization of nitrogen from GABA and allantoin, whereas in the latter, it regulates expression of arginine synthesis genes. Our findings make an important contribution to our understanding of nitrogen regulation in a human-pathogenic fungus.

Using the Candida Genome Database

Studying Candida biology requires access to genomic sequence data in conjunction with experimental information that together provide functional context to genes and proteins, and aid in interpreting newly generated experimental data. The Candida Genome Database (CGD) curates the Candida literature, and integrates functional information about Candida genes and their products with a set of analysis tools that facilitate searching for sets of genes and exploring their biological roles. This chapter describes how the various types of information available at CGD can be searched, retrieved, and analyzed. Starting with the guided tour of the CGD Home page and Locus Summary page, this unit shows how to navigate the various assemblies of the C. albicans genome, how to use Gene Ontology tools to make sense of large-scale data, and how to access the microarray data archived at CGD, as well as visualize high-throughput sequencing data through the use of JBrowse.

Multiple Approaches to Phylogenomic Reconstruction of the Fungal Kingdom

Fungi are possibly the most diverse eukaryotic kingdom, with over a million member species and an evolutionary history dating back a billion years. Fungi have been at the forefront of eukaryotic genomics, and owing to initiatives like the 1000 Fungal Genomes Project the amount of fungal genomic data has increased considerably over the last 5 years, enabling large-scale comparative genomics of species across the kingdom. In this chapter, we first review fungal evolution and the history of fungal genomics. We then review in detail seven phylogenomic methods and reconstruct the phylogeny of 84 fungal species from 8 phyla using each method. Six methods have seen extensive use in previous fungal studies, while a Bayesian supertree method is novel to fungal phylogenomics. We find that both established and novel phylogenomic methods can accurately reconstruct the fungal kingdom. Finally, we discuss the accuracy and suitability of each phylogenomic method utilized.

An Update on Candida tropicalis Based on Basic and Clinical Approaches

Candida tropicalis has emerged as one of the most important Candida species. It has been widely considered the second most virulent Candida species, only preceded by C. albicans. Besides, this species has been recognized as a very strong biofilm producer, surpassing C. albicans in most of the studies. In addition, it produces a wide range of other virulence factors, including: adhesion to buccal epithelial and endothelial cells; the secretion of lytic enzymes, such as proteinases, phospholipases, and hemolysins, bud-to-hyphae transition (also called morphogenesis) and the phenomenon called phenotypic switching. This is a species very closely related to C. albicans and has been easily identified with both phenotypic and molecular methods. In addition, no cryptic sibling species were yet described in the literature, what is contradictory to some other medically important Candida species. C. tropicalis is a clinically relevant species and may be the second or third etiological agent of candidemia, specifically in Latin American countries and Asia. Antifungal resistance to the azoles, polyenes, and echinocandins has already been described. Apart from all these characteristics, C. tropicalis has been considered an osmotolerant microorganism and this ability to survive to high salt concentration may be important for fungal persistence in saline environments. This physiological characteristic makes this species suitable for use in biotechnology processes. Here we describe an update of C. tropicalis, focusing on all these previously mentioned subjects.

Eukaryotic transporters for hydroxyderivatives of benzoic acid

Several yeast species catabolize hydroxyderivatives of benzoic acid. However, the nature of carriers responsible for transport of these compounds across the plasma membrane is currently unknown. In this study, we analyzed a family of genes coding for permeases belonging to the major facilitator superfamily (MFS) in the pathogenic yeast Candida parapsilosis. Our results revealed that these transporters are functionally equivalent to bacterial aromatic acid: H⁺ symporters (AAHS) such as GenK, MhbT and PcaK. We demonstrate that the genes HBT1 and HBT2 encoding putative transporters are highly upregulated in C. parapsilosis cells assimilating hydroxybenzoate substrates and the corresponding proteins reside in the plasma membrane. Phenotypic analyses of knockout mutants and hydroxybenzoate uptake assays provide compelling evidence that the permeases Hbt1 and Hbt2 transport the substrates that are metabolized via the gentisate (3-hydroxybenzoate, gentisate) and 3-oxoadipate pathway (4-hydroxybenzoate, 2,4-dihydroxybenzoate and protocatechuate), respectively. Our data support the hypothesis that the carriers belong to the AAHS family of MFS transporters. Phylogenetic analyses revealed that the orthologs of Hbt permeases are widespread in the subphylum Pezizomycotina, but have a sparse distribution among Saccharomycotina lineages. Moreover, these analyses shed additional light on the evolution of biochemical pathways involved in the catabolic degradation of hydroxyaromatic compounds.

Functional specialization of two paralogous TAF12 variants by their selective association with SAGA and TFIID transcriptional regulatory complexes

TATA box-binding protein (TBP)-associated factors (TAFs), evolutionarily conserved from yeast to humans, play a central role during transcription initiation. A subset of TAF proteins is shared in transcription factor II D (TFIID) and SAGA transcription regulatory complexes. Although higher eukaryotes contain multiple TAF variants that specify tissue- and developmental stage-specific organization of TFIID or SAGA complexes, in unicellular genomes, however, each TAF is encoded by a single gene. Surprisingly, we found that the genome of Candida albicans, the predominant human fungal pathogen, contains two paralogous TAF12 genes, CaTAF12L and CaTAF12, encoding H2B-like histone-fold domain-containing variants. Of the available fungal genome sequences, only seven other closely related diploid pathogenic Candida genomes encode the two TAF12 paralogs. Using affinity purifications from C. albicans cell extracts, we demonstrate that CaTAF12L uniquely associates with the SAGA complex and CaTAF12 associates with the TFIID complex. We further show that CaTAF12, but not CaTAF12L, is essential for C. albicans growth. Conditional depletion of the two TAF12 variant proteins caused distinct cellular and colony phenotypes. Together our results define a specialized organization of the TAF12 variants and non-redundant roles for the two TAF12 variants in the unicellular C. albicans genome.

Specific Human and Candida Cellular Interactions Lead to Controlled or Persistent Infection Outcomes during Granuloma-Like Formation

A delayed type of multicellular process could be crucial during chronic candidiasis in determining the course of infection. This reaction, consisting of organized immune cells surrounding the pathogen, initiates an inflammatory response to avoid fungal dissemination. The goal of the present study was to examine, at an in vitro cellular scale, Candida and human immune cell interaction dynamics during a long-term period. By challenging human peripheral blood immune cells from 10 healthy donors with 32 Candida albicans and non-albicans (C. glabrata, C. tropicalis, C. parapsilosis, C. dubliniensis, C. lusitaniae, C. krusei, and C. kefyr) clinical isolates, we showed that Candida spp. induced the formation of granuloma-like structures within 6 days after challenge, but their sizes and the respective fungal burdens differed according to the Candida species. These two parameters are positively correlated. Phenotypic characteristics, such as hypha formation and higher axenic growth rate, seem to contribute to yeast persistence within granuloma-like structures. We showed an interindividual variability of the human response against Candida spp. Higher proportions of neutrophils and elevated CD4⁺/CD8⁺ T cell ratios during the first days after challenge were correlated with early production of gamma interferon (IFN-γ) and associated with controlled infection. In contrast, the persistence of Candida could result from upregulation of proinflammatory cytokines such as interleukin-6 (IL-6), IFN-γ, and tumor necrosis factor alpha (TNF-α) and a poor anti-inflammatory negative feedback (IL-10). Importantly, regulatory subsets of NK cells and CD4^lo CD8^hi doubly positive (DP) lymphocytes at late stage infiltrate granuloma-like structures and could correlate with the IL-10 and TNF-α production. These data offer a base frame to explain cellular events that guide infection control or fungal persistence.

Mitochondrial Carriers Link the Catabolism of Hydroxyaromatic Compounds to the Central Metabolism in Candida parapsilosis

The pathogenic yeast Candida parapsilosis metabolizes hydroxyderivatives of benzene and benzoic acid to compounds channeled into central metabolism, including the mitochondrially localized tricarboxylic acid cycle, via the 3-oxoadipate and gentisate pathways. The orchestration of both catabolic pathways with mitochondrial metabolism as well as their evolutionary origin is not fully understood. Our results show that the enzymes involved in these two pathways operate in the cytoplasm with the exception of the mitochondrially targeted 3-oxoadipate CoA-transferase (Osc1p) and 3-oxoadipyl-CoA thiolase (Oct1p) catalyzing the last two reactions of the 3-oxoadipate pathway. The cellular localization of the enzymes indicates that degradation of hydroxyaromatic compounds requires a shuttling of intermediates, cofactors, and products of the corresponding biochemical reactions between cytosol and mitochondria. Indeed, we found that yeast cells assimilating hydroxybenzoates increase the expression of genes SFC1, LEU5, YHM2, and MPC1 coding for succinate/fumarate carrier, coenzyme A carrier, oxoglutarate/citrate carrier, and the subunit of pyruvate carrier, respectively. A phylogenetic analysis uncovered distinct evolutionary trajectories for sparsely distributed gene clusters coding for enzymes of both pathways. Whereas the 3-oxoadipate pathway appears to have evolved by vertical descent combined with multiple losses, the gentisate pathway shows a striking pattern suggestive of horizontal gene transfer to the evolutionarily distant Mucorales.

N-acetylglucosamine sensing by a GCN5-related N-acetyltransferase induces transcription via chromatin histone acetylation in fungi

N-acetylglucosamine (GlcNAc) exists ubiquitously as a component of the surface on a wide range of cells, from bacteria to humans. Many fungi are able to utilize environmental GlcNAc to support growth and induce cellular development, a property important for their survival in various host niches. However, how the GlcNAc signal is sensed and subsequently transduced is largely unknown. Here, we identify a gene that is essential for GlcNAc signalling (NGS1) in Candida albicans, a commensal and pathogenic yeast of humans. Ngs1 can bind GlcNAc through the N-terminal β-N-acetylglucosaminidase homology domain. This binding activates N-acetyltransferase activity in the C-terminal GCN5-related N-acetyltransferase domain, which is required for GlcNAc-induced promoter histone acetylation and transcription. Ngs1 is targeted to the promoters of GlcNAc-inducible genes constitutively by the transcription factor Rep1. Ngs1 is conserved in diverse fungi that have GlcNAc catabolic genes. Thus, fungi use Ngs1 as a GlcNAc-sensor and transducer for GlcNAc-induced transcription.

Many fungi are able to metabolise environmental N-acetylglucosamine, however the mechanism by which this molecule is sensed is unclear. Su and Lu et al. show that Candida albicans NGS1 fulfils this function by mediating N-acetylglucosamine-dependent histone acetylation at target genes.

Metabolic profiles of planktonic and biofilm cells of Candida orthopsilosis

AIM

This study aims to understand which Candida orthopsilosis protein aids fungus adaptation upon its switching from planktonic growth to biofilm.

MATERIALS & METHODS

Ion mobility separation within mass spectrometry analysis combination were used.

RESULTS

Proteins mapped for different biosynthetic pathways showed that selective ribosome autophagy might occur in biofilms. Glucose, used as a carbon source in the glycolytic flux, changed to glycogen and trehalose.

CONCLUSION

Candida orthopsilosis expresses proteins that combine a variety of mechanisms to provide yeasts with the means to adjust the catalytic properties of enzymes. Adjustment of the enzymes helps modulate the biosynthesis/degradation rates of the available nutrients, in order to control and coordinate the metabolic pathways that enable cells to express an adequate response to nutrient availability.

Identification of Non-Coding RNAs in the Candida parapsilosis Species Group

The Candida CTG clade is a monophyletic group of fungal species that translates CTG as serine, and includes the pathogens Candida albicans and Candida parapsilosis. Research has typically focused on identifying protein-coding genes in these species. Here, we use bioinformatic and experimental approaches to annotate known classes of non-coding RNAs in three CTG-clade species, Candida parapsilosis, Candida orthopsilosis and Lodderomyces elongisporus. We also update the annotation of ncRNAs in the C. albicans genome. The majority of ncRNAs identified were snoRNAs. Approximately 50% of snoRNAs (including most of the C/D box class) are encoded in introns. Most are within mono- and polycistronic transcripts with no protein coding potential. Five polycistronic clusters of snoRNAs are highly conserved in fungi. In polycistronic regions, splicing occurs via the classical pathway, as well as by nested and recursive splicing. We identified spliceosomal small nuclear RNAs, the telomerase RNA component, signal recognition particle, RNase P RNA component and the related RNase MRP RNA component in all three genomes. Stem loop IV of the U2 spliceosomal RNA and the associated binding proteins were lost from the ancestor of C. parapsilosis and C. orthopsilosis, following the divergence from L. elongisporus. The RNA component of the MRP is longer in C. parapsilosis, C. orthopsilosis and L. elongisporus than in S. cerevisiae, but is substantially shorter than in C. albicans.

Biotin Auxotrophy and Biotin Enhanced Germ Tube Formation in Candida albicans

Due to the increased number of immunocompromised patients, infections with the pathogen Candida albicans have significantly increased in recent years. C. albicans transition from yeast to germ tubes is one of the essential factors for virulence. In this study we noted that Lee's medium, commonly used to induce filamentation, contained 500-fold more biotin than needed for growth and 40-fold more biotin than is typically added to growth media. Thus, we investigated the effects of excess biotin on growth rate and filamentation by C. albicans in different media. At 37 °C, excess biotin (4 µM) enhanced germ tube formation (GTF) ca. 10-fold in both Lee's medium and a defined glucose-proline medium, and ca. 4-fold in 1% serum. Two biotin precursors, desthiobiotin and 7-keto-8-aminopelargonic acid (KAPA), also stimulated GTF. During these studies we also noted an inverse correlation between the number of times the inoculum had been washed and the concentration of serum needed to stimulate GTF. C. albicans cells that had been washed eight times achieved 80% GTF with only 0.1% sheep serum. The mechanism by which 1-4 µM biotin enhances GTF is still unknown except to note that equivalent levels of biotin are needed to create an internal supply of stored biotin and biotinylated histones. Biotin did not restore filamentation for any of the four known filamentation defective mutants tested. C. albicans is auxotrophic for biotin and this biotin auxotrophy was fulfilled by biotin, desthiobiotin, or KAPA. However, biotin auxotrophy is not temperature dependent or influenced by the presence of 5% CO₂. Biotin starvation upregulated the biotin biosynthetic genes BIO2, BIO3, and BIO4 by 11-, 1500-, and 150-fold, respectively, and BIO2p is predicted to be mitochondrion-localized. Based on our findings, we suggest that biotin has two roles in the physiology of C. albicans, one as an enzymatic cofactor and another as a morphological regulator. Finally, we found no evidence supporting prior claims that C. albicans only forms hyphae at very low biotin (0.1 nM) growth conditions.

Transcriptional rewiring over evolutionary timescales changes quantitative and qualitative properties of gene expression

Evolutionary changes in transcription networks are an important source of diversity across species, yet the quantitative consequences of network evolution have rarely been studied. Here we consider the transcriptional 'rewiring' of the three GAL genes that encode the enzymes needed for cells to convert galactose to glucose. In Saccharomyces cerevisiae, the transcriptional regulator Gal4 binds and activates these genes. In the human pathogen Candida albicans (which last shared a common ancestor with S. cerevisiae some 300 million years ago), we show that different regulators, Rtg1 and Rtg3, activate the three GAL genes. Using single-cell dynamics and RNA-sequencing, we demonstrate that although the overall logic of regulation is the same in both species-the GAL genes are induced by galactose-there are major differences in both the quantitative response of these genes to galactose and in the position of these genes in the overall transcription network structure of the two species.

Candida albicans Agglutinin-Like Sequence (Als) Family Vignettes: A Review of Als Protein Structure and Function

Approximately two decades have passed since the description of the first gene in the Candida albicans ALS (agglutinin-like sequence) family. Since that time, much has been learned about the composition of the family and the function of its encoded cell-surface glycoproteins. Solution of the structure of the Als adhesive domain provides the opportunity to evaluate the molecular basis for protein function. This review article is formatted as a series of fundamental questions and explores the diversity of the Als proteins, as well as their role in ligand binding, aggregative effects, and attachment to abiotic surfaces. Interaction of Als proteins with each other, their functional equivalence, and the effects of protein abundance on phenotypic conclusions are also examined. Structural features of Als proteins that may facilitate invasive function are considered. Conclusions that are firmly supported by the literature are presented while highlighting areas that require additional investigation to reveal basic features of the Als proteins, their relatedness to each other, and their roles in C. albicans biology.

Using RNA-seq for Analysis of Differential Gene Expression in Fungal Species

The ability to extract, identify and annotate large amounts of biological data is a key feature of the "omics" era, and has led to an explosion in the amount of data available. One pivotal advance is the use of Next-Generation Sequencing (NGS) techniques such as RNA-sequencing (RNA-seq). RNA-seq uses data from millions of small mRNA transcripts or "reads" which are aligned to a reference genome. Comparative transcriptomics analyses using RNA-seq can provide the researcher with a comprehensive view of the cells' response to a given environment or stimulus.Here, we describe the NGS techniques (based on Illumina technology) that are routinely used for comparative transcriptome analysis of fungal species. We describe the entire process from isolation of RNA to computational identification of differentially expressed genes. We provide instructions to allow the beginner to implement packages in R such as Bioconductor. The methods described are not limited to yeast, and can also be applied to other eukaryotic organisms.

How to Use the Candida Genome Database

Studying Candida biology requires access to genomic sequence data in conjunction with experimental information that provides functional context to genes and proteins. The Candida Genome Database (CGD) integrates functional information about Candida genes and their products with a set of analysis tools that facilitate searching for sets of genes and exploring their biological roles. This chapter describes how the various types of information available at CGD can be searched, retrieved, and analyzed. Starting with the guided tour of the CGD Home page and Locus Summary page, this unit shows how to navigate the various assemblies of the C. albicans genome, how to use Gene Ontology tools to make sense of large-scale data, and how to access the microarray data archived at CGD.

Models for the evolution of GC content in asexual fungi Candida albicans and C. dubliniensis

Although guanine–cytosine (GC)-biased gene conversion (gBGC) following meiotic recombination seems the most probable mechanism accounting for large-scale variations in GC content for many eukaryotes, it cannot explain such variations for organisms belonging to ancient asexual lineages, such as the pathogenic fungi Candida albicans and C. dubliniensis. Analysis of the substitution patterns for these two species reveals a strong anticorrelation between the synonymous transition rates at third codon positions. I propose two models that can account for this observation. According to the first model, the evolution of GC content is driven by gBGC linked to mitotic recombination, either associated with parasexuality or with damage repair. Variations in the GC content thus reflect variations in the strength of gBGC, presumably variations in the mitotic recombination rate. According to the second model, the evolution of GC content is driven by misincorporation errors during the process of DNA replication in S phase. This model proposes that variations in GC content are due to variations in the proportions of dCTPs and dGTPs at the time when sequences are replicated. Experimental data regarding mitotic recombination rates or the variations of dCTPs and dGTPs during S phase are required to validate definitively one of the two models, but in any case, the fit of the models to the data suggests that C. albicans and C. dubliniensis constitute so far unique examples of GC content evolution driven either by mitotic recombination or replicative errors.

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.

The complex evolutionary dynamics of Hsp70s: a genomic and functional perspective

Hsp70 molecular chaperones are ubiquitous. By preventing aggregation, promoting folding, and regulating degradation, Hsp70s are major factors in the ability of cells to maintain proteostasis. Despite a wealth of functional information, little is understood about the evolutionary dynamics of Hsp70s. We undertook an analysis of Hsp70s in the fungal clade Ascomycota. Using the well-characterized 14 Hsp70s of Saccharomyces cerevisiae, we identified 491 orthologs from 53 genomes. Saccharomyces cerevisiae Hsp70s fall into seven subfamilies: four canonical-type Hsp70 chaperones (SSA, SSB, KAR, and SSC) and three atypical Hsp70s (SSE, SSZ, and LHS) that play regulatory roles, modulating the activity of canonical Hsp70 partners. Each of the 53 surveyed genomes harbored at least one member of each subfamily, and thus establishing these seven Hsp70s as units of function and evolution. Genomes of some species contained only one member of each subfamily that is only seven Hsp70s. Overall, members of each subfamily formed a monophyletic group, suggesting that each diversified from their corresponding ancestral gene present in the common ancestor of all surveyed species. However, the pattern of evolution varied across subfamilies. At one extreme, members of the SSB subfamily evolved under concerted evolution. At the other extreme, SSA and SSC subfamilies exhibited a high degree of copy number dynamics, consistent with a birth–death mode of evolution. KAR, SSE, SSZ, and LHS subfamilies evolved in a simple divergent mode with little copy number dynamics. Together, our data revealed that the evolutionary history of this highly conserved and ubiquitous protein family was surprising complex and dynamic.

Filling annotation gaps in yeast genomes using genome-wide contact maps

MOTIVATIONS

De novo sequencing of genomes is followed by annotation analyses aiming at identifying functional genomic features such as genes, non-coding RNAs or regulatory sequences, taking advantage of diverse datasets. These steps sometimes fail at detecting non-coding functional sequences: for example, origins of replication, centromeres and rDNA positions have proven difficult to annotate with high confidence. Here, we demonstrate an unconventional application of Chromosome Conformation Capture (3C) technique, which typically aims at deciphering the average 3D organization of genomes, by showing how functional information about the sequence can be extracted solely from the chromosome contact map.

RESULTS

Specifically, we describe a combined experimental and bioinformatic procedure that determines the genomic positions of centromeres and ribosomal DNA clusters in yeasts, including species where classical computational approaches fail. For instance, we determined the centromere positions in Naumovozyma castellii, where these coordinates could not be obtained previously. Although computed centromere positions were characterized by conserved synteny with neighboring species, no consensus sequences could be found, suggesting that centromeric binding proteins or mechanisms have significantly diverged. We also used our approach to refine centromere positions in Kuraishia capsulata and to identify rDNA positions in Debaryomyces hansenii. Our study demonstrates how 3C data can be used to complete the functional annotation of eukaryotic genomes.

AVAILABILITY AND IMPLEMENTATION

The source code is provided in the Supplementary Material. This includes a zipped file with the Python code and a contact matrix of Saccharomyces cerevisiae.

CONTACT

romain.koszul@pasteur.fr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A family of glutathione peroxidases contributes to oxidative stress resistance in Candida albicans

Candida albicans is a well-adapted human commensal but is also a facultative pathogen that can cause superficial and systemic infections. Its remarkable capacity to thrive within the human host relies on its ability to adapt and respond to the local environment of different niches. C. albicans is able to cope with oxidative stress in a coordinated fashion via upregulation of different protective mechanisms. Here, we unravel the role of a family of glutathione peroxidase (GPx), designated Gpx31, Gpx32, and Gpx33, in oxidative stress resistance. We show that GPx activity in C. albicans is induced upon exposure to peroxides and that this enzymatic activity is required for full resistance to oxidative stress. The GPx activity relies on the presence of GPX31, with no apparent contribution from GPX32 and GPX33 during in vitro short-term (3 h) exposure to peroxides. However, a triple gpx31-33Δ/Δ mutant exhibited a more pronounced sensitivity than a single gpx31Δ/Δ mutant on solid media in the presence of oxidants, suggesting that GPX32 and GPX33 may be involved in long-term adaptation to oxidative stress. Interestingly, reintegration of a single allele of GPX31 was sufficient to restore the wild-type phenotype in both the single and triple mutants. We found that mutants lacking GPX31-33 were more susceptible to killing by phagocytic cells, suggesting that GPxs are required for full resistance to innate immune effector cells. Despite the sensitivity to oxidative stress and phagocytes, these mutants were not affected in their virulence in the chicken embryo model of candidiasis.

A pluralistic account of homology: adapting the models to the data

Defining homologous genes is important in many evolutionary studies but raises obvious issues. Some of these issues are conceptual and stem from our assumptions of how a gene evolves, others are practical, and depend on the algorithmic decisions implemented in existing software. Therefore, to make progress in the study of homology, both ontological and epistemological questions must be considered. In particular, defining homologous genes cannot be solely addressed under the classic assumptions of strong tree thinking, according to which genes evolve in a strictly tree-like fashion of vertical descent and divergence and the problems of homology detection are primarily methodological. Gene homology could also be considered under a different perspective where genes evolve as “public goods,” subjected to various introgressive processes. In this latter case, defining homologous genes becomes a matter of designing models suited to the actual complexity of the data and how such complexity arises, rather than trying to fit genetic data to some a priori tree-like evolutionary model, a practice that inevitably results in the loss of much information. Here we show how important aspects of the problems raised by homology detection methods can be overcome when even more fundamental roots of these problems are addressed by analyzing public goods thinking evolutionary processes through which genes have frequently originated. This kind of thinking acknowledges distinct types of homologs, characterized by distinct patterns, in phylogenetic and nonphylogenetic unrooted or multirooted networks. In addition, we define “family resemblances” to include genes that are related through intermediate relatives, thereby placing notions of homology in the broader context of evolutionary relationships. We conclude by presenting some payoffs of adopting such a pluralistic account of homology and family relationship, which expands the scope of evolutionary analyses beyond the traditional, yet relatively narrow focus allowed by a strong tree-thinking view on gene evolution.