Assembly of a phased diploid Candida albicans genome facilitates allele-specific measurements and provides a simple model for repeat and indel structure

A Proteogenomic Pipeline for the Analysis of Protein Biosynthesis Errors in the Human Pathogen Candida albicans

Candida albicans is a diploid pathogen known for its ability to live as a commensal fungus in healthy individuals but causing both superficial infections and disseminated candidiasis in immunocompromised patients where it is associated with high morbidity and mortality. Its success in colonizing the human host is attributed to a wide range of virulence traits that modulate interactions between the host and the pathogen, such as optimal growth rate at 37 °C, the ability to switch between yeast and hyphal forms, and a remarkable genomic and phenotypic plasticity. A fascinating aspect of its biology is a prominent heterogeneous proteome that arises from frequent genomic rearrangements, high allelic variation, and high levels of amino acid misincorporations in proteins. This leads to increased morphological and physiological phenotypic diversity of high adaptive potential, but the scope of such protein mistranslation is poorly understood due to technical difficulties in detecting and quantifying amino acid misincorporation events in complex protein samples. We have developed and optimized mass spectrometry and bioinformatics pipelines capable of identifying rare amino acid misincorporation events at the proteome level. We have also analyzed the proteomic profile of an engineered C. albicans strain that exhibits high level of leucine misincorporation at protein CUG sites and employed an in vivo quantitative gain-of-function fluorescence reporter system to validate our LC-MS/MS data. C. albicans misincorporates amino acids above the background level at protein sites of diverse codons, particularly at CUG, confirming our previous data on the quantification of leucine incorporation at single CUG sites of recombinant reporter proteins, but increasing misincorporation of Leucine at these sites does not alter the translational fidelity of the other codons. These findings indicate that the C. albicans statistical proteome exceeds prior estimates, suggesting that its highly plastic phenome may also be modulated by environmental factors due to translational ambiguity.

Candida albicans natural diversity: a resource to dissect fungal commensalism and pathogenesis

Candida albicans is a ubiquitous fungus of humans. It is not only a component of the oral and intestinal microbiota of most healthy adults but also a major cause of mucosal disorders and life-threatening disseminated infections. Until recently, research on the biology and pathogenesis of the fungus was largely based on a single clinical isolate. We review investigations that have started to dissect a diverse set of C. albicans strains. Using different approaches to leverage the species' phenotypic and/or genetic diversity, these studies illuminate the wide range of interactions between fungus and host. While connecting genetic variants to phenotypes of interest remains challenging, research on C. albicans' natural diversity is central to understand fungal commensalism and pathogenesis.

Step-wise evolution of azole resistance through copy number variation followed by KSR1 loss of heterozygosity in Candida albicans

Antimicrobial drug resistance poses a global health threat, requiring a deeper understanding of the evolutionary processes that lead to its emergence in pathogens. Complex evolutionary dynamics involve multiple mutations that can result in cooperative or competitive (clonal interference) effects. Candida albicans, a major fungal pathogen, displays high rates of copy number variation (CNV) and loss of heterozygosity (LOH). CNV and LOH events involve large numbers of genes and could synergize during evolutionary adaptation. Understanding the contributions of CNV and LOH to antifungal drug adaptation is challenging, especially in the context of whole-population genome sequencing. Here, we document the sequential evolution of fluconazole tolerance and then resistance in a C. albicans isolate involving an initial CNV on chromosome 4, followed by an LOH on chromosome R that involves KSR1. Similar LOH events involving KSR1, which encodes a reductase in the sphingolipid biosynthesis pathway, were also detected in independently evolved fluconazole resistant isolates. We dissect the specific KSR1 codons that affect fluconazole resistance and tolerance. The combination of the chromosome 4 CNV and KSR1 LOH results in a >500-fold decrease in azole susceptibility relative to the progenitor, illustrating a compelling example of rapid, yet step-wise, interplay between CNV and LOH in drug resistance evolution.

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.

The anti-staphylococcal fusidic acid as an efflux pump inhibitor combined with fluconazole against vaginal candidiasis in mouse model

BACKGROUND

Candida albicans is the most common fungus that causes vaginal candidiasis in immunocompetent women and catastrophic infections in immunocompromised patients. The treatment of such infections is hindered due to the increasing emergence of resistance to azoles in C. albicans. New treatment approaches are needed to combat candidiasis especially in the dwindled supply of new effective and safe antifungals. The resistance to azoles is mainly attributed to export of azoles outside the cells by means of the efflux pump that confers cross resistance to all azoles including fluconazole (FLC).

OBJECTIVES

This study aimed to investigate the possible efflux pump inhibiting activity of fusidic acid (FA) in C. albicans resistant isolates and the potential use of Fusidic acid in combination with fluconazole to potentiate the antifungal activity of fluconazole to restore its activity in the resistant C. albicans isolates.

METHODS

The resistance of C. albicans isolates was assessed by determination of minimum inhibitory concentration. The effect of Fusidic acid at sub-inhibitory concentration on efflux activity was assayed by rhodamine 6G efflux assay and intracellular accumulation. Mice model studies were conducted to evaluate the anti-efflux activity of Fusidic acid and its synergistic effects in combination with fluconazole. Impact of Fusidic acid on ergosterol biosynthesis was quantified. The synergy of fluconazole when combined with Fusidic acid was investigated by determination of minimum inhibitory concentration. The cytotoxicity of Fusidic acid was tested against erythrocytes. The effect of Fusidic acid on efflux pumps was tested at the molecular level by real-time PCR and in silico study. In vivo vulvovaginitis mice model was used to confirm the activity of the combination in treating vulvovaginal candidiasis.

RESULTS

Fusidic acid showed efflux inhibiting activity as it increased the accumulation of rhodamine 6G, a substrate for ABC-efflux transporter, and decreased its efflux in C. albicans cells. The antifungal activity of fluconazole was synergized when combined with Fusidic acid. Fusidic acid exerted only minimal cytotoxicity on human erythrocytes indicating its safety. The FA efflux inhibitory activity could be owed to its ability to interfere with efflux protein transporters as revealed by docking studies and downregulation of the efflux-encoding genes of both ABC transporters and MFS superfamily. Moreover, in vivo mice model showed that using fluconazole-fusidic acid combination by vaginal route enhanced fluconazole antifungal activity as shown by lowered fungal burden and a negligible histopathological change in vaginal tissue.

CONCLUSION

The current findings highlight FA's potential as a potential adjuvant to FLC in the treatment of vulvovaginal candidiasis.

Use of a Candida albicans SC5314 PacBio HiFi reads dataset to close gaps in the reference genome assembly, reveal a subtelomeric gene family, and produce accurate phased allelic sequences

Candida albicans SC5314 is the most-often used strain for molecular manipulation of the species. The SC5314 reference genome sequence is the result of considerable effort from many scientists and has advanced research into fungal biology and pathogenesis. Although the resource is highly developed and presented in a phased diploid format, the sequence includes gaps and does not extend to the telomeres on its eight chromosome pairs. Accurate SC5314 genome assembly is complicated by the presence of extensive repeated sequences and considerable allelic length variation at some loci. Advances in genome sequencing technology provide the tools to obtain highly accurate long-read data that span even the most-difficult-to-assemble genome regions. Here, we describe derivation of a PacBio HiFi data set and creation of a collapsed haploid telomere-to-telomere assembly of the SC5314 genome (ASM3268872v1) that revealed previously unknown features of the strain. ASM3268872v1 subtelomeric distances were up to 19 kb larger than in the reference genome and revealed a family of highly conserved DNA helicase-encoding genes at 10 of the 16 chromosome ends. We also describe alignments of individual HiFi reads to deduce accurate diploid sequences for the most notoriously difficult-to-assemble C. albicans genes: the agglutinin-like sequence (ALS) gene family. We provide a tutorial that demonstrates how the HiFi reads can be visualized to explore any region of interest. Availability of the HiFi reads data set and the ASM3268872v1 comparative guide assembly will streamline research efforts because accurate diploid sequences can be derived using simple in silico methods rather than time-consuming laboratory-bench approaches.

Recurrent copy number variations in the human fungal pathogen Candida parapsilosis

Candida parapsilosis is an opportunistic fungal pathogen with increasing incidence in hospital settings worldwide; however, we lack a comprehensive understanding of the mechanisms promoting its virulence and drug resistance. Bergin et al. systematically quantify the frequency and effect of copy number variation (CNV) across 170 diverse clinical and environmental isolates of C. parapsilosis (Bergin SA, Zhao F, Ryan AP, Müller CA, Nieduszynski CA, Zhai B, Rolling T, Hohl TM, Morio F, Scully J, Wolfe KH, Butler G, 2022, mBio, https://doi.org/10.1128/mbio.01777-22). Using a combination of both short- and long-read whole genome sequencing techniques, they determine the structure and copy number of two CNVs that arose recurrently throughout the evolution of these isolates. Each CNV predominantly amplifies one coding sequence (ARR3 or RTA3); however, the amplitude and recombination breakpoints are variable across the isolates. Amplification of RTA3 correlates with drug resistance and deletion causes drug susceptibility. This study highlights the need for further research into the mechanisms and dynamics of CNV formation and the impact of these CNVs on virulence and drug resistance across diverse fungal pathogens.

Inter-species association mapping links splice site evolution to METTL16 and SNRNP27K

Eukaryotic genes are interrupted by introns that are removed from transcribed RNAs by splicing. Patterns of splicing complexity differ between species, but it is unclear how these differences arise. We used inter-species association mapping with Saccharomycotina species to correlate splicing signal phenotypes with the presence or absence of splicing factors. Here, we show that variation in 5' splice site sequence preferences correlate with the presence of the U6 snRNA N6-methyladenosine methyltransferase METTL16 and the splicing factor SNRNP27K. The greatest variation in 5' splice site sequence occurred at the +4 position and involved a preference switch between adenosine and uridine. Loss of METTL16 and SNRNP27K orthologs, or a single SNRNP27K methionine residue, was associated with a preference for +4 U. These findings are consistent with splicing analyses of mutants defective in either METTL16 or SNRNP27K orthologs and models derived from spliceosome structures, demonstrating that inter-species association mapping is a powerful orthogonal approach to molecular studies. We identified variation between species in the occurrence of two major classes of 5' splice sites, defined by distinct interaction potentials with U5 and U6 snRNAs, that correlates with intron number. We conclude that variation in concerted processes of 5' splice site selection by U6 snRNA is associated with evolutionary changes in splicing signal phenotypes.

Genome-Wide DNA Changes Acquired by Candida albicans Caspofungin-Adapted Mutants

Drugs from the echinocandin (ECN) class are now recommended 'front-line' treatments of infections caused by a prevailing fungal pathogen, C. albicans. However, the increased use of ECNs is associated with a rising resistance to ECNs. As the acquisition of ECN resistance in C. albicans is viewed as a multistep evolution, determining factors that are associated with the decreased ECN susceptibility is of importance. We have recently identified two cohorts of genes that are either up- or downregulated in concert in order to control remodeling of cell wall, an organelle targeted by ECNs, in laboratory mutants with decreased ECN susceptibility. Here, we profiled the global DNA sequence of four of these adapted mutants in search of DNA changes that are associated with decreased ECN susceptibility. We find a limited number of 112 unique mutations representing two alternative mutational pathways. Approximately half of the mutations occurred as hotspots. Approximately half of mutations and hotspots were shared by ECN-adapted mutants despite the mutants arising as independent events and differing in some of their phenotypes, as well as in condition of chromosome 5. A total of 88 mutations are associated with 43 open reading frames (ORFs) and occurred inside of an ORF or within 1 kb of an ORF, predominantly as single-nucleotide substitution. Mutations occurred more often in the 5'-UTR than in the 3'-UTR by a 1.67:1 ratio. A total of 16 mutations mapped to eight genomic features that were not ORFs: Tca4-4 retrotransposon; Tca2-7 retrotransposon; lambda-4a long terminal repeat; mu-Ra long terminal repeat; MRS-7b Major Repeat Sequence; MRS-R Major Repeat Sequence; RB2-5a repeat sequence; and tL (CAA) leucine tRNA. Finally, eight mutations are not associated with any ORF or other genomic feature. Repeated occurrence of single-nucleotide substitutions in non-related drug-adapted mutants strongly indicates that these DNA changes are accompanying drug adaptation and could possibly influence ECN susceptibility, thus serving as factors facilitating evolution of ECN drug resistance due to classical mutations in FKS1.

High-throughput functional profiling of the human fungal pathogen Candida albicans genome

Candida albicans is a major fungal pathogen of humans. Although its genome has been sequenced more than two decades ago, there are still over 4300 uncharacterized C. albicans genes. We previously generated an ORFeome as well as a collection of destination vectors to facilitate overexpression of C. albicans ORFs. Here, we report the construction of ∼2500 overexpression mutants and their evaluation by in vitro spotting on rich medium and in a liquid pool experiment in rich medium, allowing the identification of genes whose overexpression has a fitness cost. The candidates were further validated at the individual strain level. This new resource allows large-scale screens in different growth conditions to be performed routinely. Altogether, based on the concept of identifying functionally related genes by cluster analysis, the availability of this overexpression mutant collection will facilitate the characterization of gene functions in C. albicans.

Arthrocolins Synergizing with Fluconazole Inhibit Fluconazole-Resistant Candida albicans by Increasing Riboflavin Metabolism and Causing Mitochondrial Dysfunction and Autophagy

Our previous study reported that seminaturally occurring arthrocolins A to C with unprecedented carbon skeletons could restore the antifungal activity of fluconazole against fluconazole-resistant Candida albicans. Here, we showed that arthrocolins synergized with fluconazole, reducing the fluconazole minimum and dramatically augmenting the survivals of 293T human cells and nematode Caenorhabditis elegans infected with fluconazole-resistant C. albicans. Mechanistically, fluconazole can induce fungal membrane permeability to arthrocolins, leading to the intracellular arthrocolins that were critical to the antifungal activity of the combination therapy by inducing abnormal cell membranes and mitochondrial dysfunctions in the fungus. Transcriptomics and reverse transcription-quantitative PCR (qRT-PCR) analysis indicated that the intracellular arthrocolins induced the strongest upregulated genes that were involved in membrane transports while the downregulated genes were responsible for fungal pathogenesis. Moreover, riboflavin metabolism and proteasomes were the most upregulated pathways, which were accompanied by inhibition of protein biosynthesis and increased levels of reactive oxygen species (ROS), lipids, and autophagy. Our results suggested that arthrocolins should be a novel class of synergistic antifungal compounds by inducing mitochondrial dysfunctions in combination with fluconazole and provided a new perspective for the design of new bioactive antifungal compounds with potential pharmacological properties. IMPORTANCE The prevalence of antifungal-resistant Candida albicans, which is a common human fungal pathogen causing life-threatening systemic infections, has become a challenge in the treatment of fungal infections. Arthrocolins are a new type of xanthene obtained from Escherichia coli fed with a key fungal precursor toluquinol. Different from those artificially synthesized xanthenes used as important medications, arthrocolins can synergize with fluconazole against fluconazole-resistant Candida albicans. Fluconazole can induce the fungal permeability of arthrocolins into fungal cells, and then the intracellular arthrocolins exerted detrimental effects on the fungus by inducing fungal mitochondrial dysfunctions, leading to dramatically reduced fungal pathogenicity. Importantly, the combination of arthrocolins and fluconazole are effective against C. albicans in two models, including human cell line 293T and nematode Caenorhabditis elegans. Arthrocolins should be a novel class of antifungal compounds with potential pharmacological properties.

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.

Small Insertions and Deletions Drive Genomic Plasticity during Adaptive Evolution of Yersinia pestis

The life cycle of Yersinia pestis has changed a lot to adapt to flea-borne transmission since it evolved from an enteric pathogen, Yersinia pseudotuberculosis. Small insertions and deletions (indels), especially frameshift mutations, can have major effects on phenotypes and contribute to virulence and host adaptation through gene disruption and inactivation. Here, we analyzed 365 Y. pestis genomes and identified 2,092 genome-wide indels on the core genome. As recently reported in Mycobacterium tuberculosis, we also detected "indel pockets" in Y. pestis, with average complexity scores declining around indel positions, which we speculate might also exist in other prokaryotes. Phylogenic analysis showed that indel-based phylogenic tree could basically reflect the phylogenetic relationships of major phylogroups in Y. pestis, except some inconsistency around the Big Bang polytomy. We observed 83 indels arising in the trunk of the phylogeny, which played a role in accumulation of pseudogenes related to key metabolism and putatively pathogenicity. We also discovered 32 homoplasies at the level of phylogroups and 7 frameshift scars (i.e., disrupted reading frame being rescued by a second frameshift). Additionally, our analysis showed evidence of parallel evolution at the level of genes, with sspA, rpoS, rnd, and YPO0624, having enriched mutations in Brazilian isolates, which might be advantageous for Y. pestis to cope with fluctuating environments. The diversified selection signals observed here demonstrates that indels are important contributors to the adaptive evolution of Y. pestis. Meanwhile, we provide potential targets for further exploration, as some genes/pseudogenes with indels we focus on remain uncharacterized. IMPORTANCE Yersinia pestis, the causative agent of plague, is a highly pathogenic clone of Yersinia pseudotuberculosis. Previous genome-wide SNP analysis provided few adaptive signatures during its evolution. Here by investigating 365 public genomes of Y. pestis, we give a comprehensive overview of general features of genome-wide indels on the core genome and their roles in Y. pestis evolution. Detection of "indel pockets," with average complexity scores declining around indel positions, in both Mycobacterium tuberculosis and Y. pestis, gives us a clue that this phenomenon might appear in other bacterial genomes. Importantly, the identification of four different forms of selection signals in indels would improve our understanding on adaptive evolution of Y. pestis, and provide targets for further physiological mechanism researches of this pathogen. As evolutionary research based on genome-wide indels is still rare in bacteria, our study would be a helpful reference in deciphering the role of indels in other species.

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

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

Molecular cloning and functional characterization of UBC13 and MMS2 from Candida albicans

To maintain genome stability, eukaryotes have evolved a powerful DNA damage response system called DNA-damage tolerance (DDT) to deal with replication-blocking lesions. In the budding yeast Saccharomyces cerevisiae, K63-linked polyubiquitination of proliferating cell nuclear antigen (PCNA) is mediated by a Ubc13-Mms2 heterodimer, leading to error-free DDT. Candida albicans is one of the most studied fungal pathogens and to date no data regarding K63-linked ubiquitination or error-free DDT has been available. Here we report the identification and functional characterization of UBC13 and MMS2 genes from C. albicans. Both genes are highly conserved between S. cerevisiae and C. albicans. However, CaUbc13 differs from all other eukaryotes in that it contains a 21-amino acid tail that appears to attenuate its interaction with CaMms2, suggesting a possible regulatory mechanism in C. albicans. Both CaUBC13 and CaMMS2 genes can functionally rescue the corresponding budding yeast mutants from increased spontaneous mutagenesis and killing by DNA-damaging agents, indicating an error-free DDT pathway in C. albicans. Indeed Caubc13Δ/Δ and Camms2Δ/Δ null mutants were constructed and displayed characteristic sensitivity to DNA-damaging agents.

Proteogenomics analysis of CUG codon translation in the human pathogen Candida albicans

Background

Yeasts of the CTG-clade lineage, which includes the human-infecting Candida albicans, Candida parapsilosis and Candida tropicalis species, are characterized by an altered genetic code. Instead of translating CUG codons as leucine, as happens in most eukaryotes, these yeasts, whose ancestors are thought to have lost the relevant leucine-tRNA gene, translate CUG codons as serine using a serine-tRNA with a mutated anticodon, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathrm{tRNA}}_{\mathrm{CAG}}^{\mathrm{Ser}} $$\end{document}tRNACAGSer. Previously reported experiments have suggested that 3–5% of the CTG-clade CUG codons are mistranslated as leucine due to mischarging of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathrm{tRNA}}_{\mathrm{CAG}}^{\mathrm{Ser}} $$\end{document}tRNACAGSer. The mistranslation was suggested to result in variable surface proteins explaining fast host adaptation and pathogenicity.

Results

In this study, we reassess this potential mistranslation by high-resolution mass spectrometry-based proteogenomics of multiple CTG-clade yeasts, including various C. albicans strains, isolated from colonized and from infected human body sites, and C. albicans grown in yeast and hyphal forms. Our data do not support a bias towards CUG codon mistranslation as leucine. Instead, our data suggest that (i) CUG codons are mistranslated at a frequency corresponding to the normal extent of ribosomal mistranslation with no preference for specific amino acids, (ii) CUG codons are as unambiguous (or ambiguous) as the related CUU leucine and UCC serine codons, (iii) tRNA anticodon loop variation across the CTG-clade yeasts does not result in any difference of the mistranslation level, and (iv) CUG codon unambiguity is independent of C. albicans’ strain pathogenicity or growth form.

Conclusions

Our findings imply that C. albicans does not decode CUG ambiguously. This suggests that the proposed misleucylation of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathrm{tRNA}}_{\mathrm{CAG}}^{\mathrm{Ser}} $$\end{document}tRNACAGSer might be as prevalent as every other misacylation or mistranslation event and, if at all, be just one of many reasons causing phenotypic diversity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01197-9.

Interactions of the Intracellular Bacterium Cardinium with Its Host, the House Dust Mite Dermatophagoides farinae, Based on Gene Expression Data

Dermatophagoides farinae is inhabited by an intracellular bacterium, Cardinium. Using correlations between host and symbiont gene expression profiles, we identified several important molecular pathways that potentially regulate/facilitate their interactions. The expression of Cardinium genes collectively explained 95% of the variation in the expression of mite genes assigned to pathways for phagocytosis, apoptosis, the MAPK signaling cascade, endocytosis, the tumor necrosis factor (TNF) pathway, the transforming growth factor beta (TGF-β) pathway, lysozyme, and the Toll/Imd pathway. In addition, expression of mite genes explained 76% of the variability in Cardinium gene expression. In particular, the expression of the Cardinium genes encoding the signaling molecules BamD, LepA, SymE, and VirD4 was either positively or negatively correlated with the expression levels of mite genes involved in endocytosis, phagocytosis, and apoptosis. We also found that Cardinium possesses a complete biosynthetic pathway for lipoic acid and may provide lipoate, but not biotin, to mites. Cardinium gene expression collectively explained 84% of the variation in expression related to several core mite metabolic pathways, and, most notably, a negative correlation was observed between bacterial gene expression and expression of mite genes assigned to the glycolysis and citric acid cycle pathways. Furthermore, we showed that Cardinium gene expression is correlated with expression levels of genes associated with terpenoid backbone biosynthesis. This pathway is important for the synthesis of pheromones, thus providing an opportunity for Cardinium to influence mite reproductive behavior to facilitate transmission of the bacterium. Overall, our study provided correlational gene expression data that can be useful for future research on mite-Cardinium interactions. IMPORTANCE The molecular mechanisms of mite-symbiont interactions and their impacts on human health are largely unknown. Astigmatid mites, such as house dust and stored-product mites, are among the most significant allergen sources worldwide. Although mites themselves are the main allergen sources, recent studies have indicated that mite-associated microbiomes may have implications for allergen production and human health. The major medically important house dust mite, D. farinae, is known to harbor a highly abundant intracellular bacterium belonging to the genus Cardinium. Expression analysis of the mite and symbiont genes can identify key mite molecular pathways that facilitate interactions with this endosymbiont and possibly shed light on how this bacterium affects mite allergen production and physiology in general.

Allele-specific assembly of a eukaryotic genome corrects apparent frameshifts and reveals a lack of nonsense-mediated mRNA decay

To date, most reference genomes represent a mosaic consensus sequence in which the homologous chromosomes are collapsed into one sequence. This approach produces sequence artefacts and impedes analyses of allele-specific mechanisms. Here, we report an allele-specific genome assembly of the diploid parasite Trypanosoma brucei and reveal allelic variants affecting gene expression. Using long-read sequencing and chromosome conformation capture data, we could assign 99.5% of all heterozygote variants to a specific homologous chromosome and build a 66 Mb long allele-specific genome assembly. The phasing of haplotypes allowed us to resolve hundreds of artefacts present in the previous mosaic consensus assembly. In addition, it revealed allelic recombination events, visible as regions of low allelic heterozygosity, enabling the lineage tracing of T. brucei isolates. Interestingly, analyses of transcriptome and translatome data of genes with allele-specific premature termination codons point to the absence of a nonsense-mediated decay mechanism in trypanosomes. Taken together, this study delivers a reference quality allele-specific genome assembly of T. brucei and demonstrates the importance of such assemblies for the study of gene expression control. We expect the new genome assembly will increase the awareness of allele-specific phenomena and provide a platform to investigate them.

Extreme diversification driven by parallel events of massive loss of heterozygosity in the hybrid lineage of Candida albicans

Candida albicans is the most commonly reported species causing candidiasis. The taxonomic classification of C. albicans and related lineages is controversial, with Candida africana (syn. C. albicans var. africana) and Candida stellatoidea (syn. C. albicans var. stellatoidea) being considered different species or C. albicans varieties depending on the authors. Moreover, recent genomic analyses have suggested a shared hybrid origin of C. albicans and C. africana, but the potential parental lineages remain unidentified. Although the genomes of C. albicans and C. africana have been extensively studied, the genome of C. stellatoidea has not been sequenced so far. In order to get a better understanding of the evolution of the C. albicans clade, and to assess whether C. stellatoidea could represent one of the unknown C. albicans parental lineages, we sequenced C. stellatoidea type strain (CBS 1905). This genome was compared to that of C. albicans and of the closely related lineage C. africana. Our results show that, similarly to C. africana, C. stellatoidea descends from the same hybrid ancestor as other C. albicans strains and that it has undergone a parallel massive loss of heterozygosity.

Reversible gene silencing through frameshift indels and frameshift scars provide adaptive plasticity for Mycobacterium tuberculosis

Mycobacterium tuberculosis can adapt to changing environments by non-heritable mechanisms. Frame-shifting insertions and deletions (indels) may also participate in adaptation through gene disruption, which could be reversed by secondary introduction of a frame-restoring indel. We present ScarTrek, a program that scans genomic data for indels, including those that together disrupt and restore a gene's reading frame, producing "frame-shift scars" suggestive of reversible gene inactivation. We use ScarTrek to analyze 5977 clinical M. tuberculosis isolates. We show that indel frequency inversely correlates with genomic linguistic complexity and varies with gene-position and gene-essentiality. Using ScarTrek, we detect 74 unique frame-shift scars in 48 genes, with a 3.74% population-level incidence of unique scar events. We find multiple scars in the ESX-1 gene cluster. Six scars show evidence of convergent evolution while the rest shared a common ancestor. Our results suggest that sequential indels are a mechanism for reversible gene silencing and adaptation in M. tuberculosis.

CRISPR mediated genome editing, a tool to dissect RNA modification processes

Though over 100 distinct RNA modifications have been identified, the roles for many of these modifications in vivo remain unknown. Genome editing is one tool investigators are using to better understand the roles these modifications play and the consequences of their absence. In this chapter, we describe how CRISPR mediated genome editing can be used to interrogate the process of RNA modification in C. albicans. Furthermore, we discuss how the protocols described can be altered to meet experimental demands. The underlying theory on which these protocols are based are applicable to a variety of model systems. The protocols described utilize the widely used S. pyogenes Cas9, but the field of genome editing is quickly evolving. We discuss the recent developments of more flexible CRISPR systems that can target a greater number of sites in the genome. These and other advancements make CRISPR mediated genome editing a practical methodology to investigate RNA modification.

The fitness costs and benefits of trisomy of each Candida albicans chromosome

Candida albicans is a prevalent human fungal pathogen. Rapid genomic change, due to aneuploidy, is a common mechanism that facilitates survival from multiple types of stresses including the few classes of available antifungal drugs. The stress survival of aneuploids occurs despite the fitness costs attributed to most aneuploids growing under idealized lab conditions. Systematic study of the aneuploid state in C. albicans has been hindered by the lack of a comprehensive collection of aneuploid strains. Here, we describe a collection of diploid C. albicans aneuploid strains, each carrying one extra copy of each chromosome, all from the same genetic background. We tested the fitness of this collection under several physiological conditions including shifts in pH, low glucose, oxidative stress, temperature, high osmolarity, membrane stress, and cell wall stress. We found that most aneuploids, under most conditions, were less fit than their euploid parent, yet there were specific conditions under which specific aneuploid isolates provided a fitness benefit relative to the euploid parent strain. Importantly, this fitness benefit was attributable to the change in the copy number of specific chromosomes. Thus, C. albicans can tolerate aneuploidy of each chromosome and some aneuploids confer improved growth under conditions that the yeast encounters in its host niches.

DeORFanizing Candida albicans Genes using Coexpression

Candida albicans is a common and deadly fungal pathogen of humans, yet the genome of this organism contains many genes of unknown function. By determining gene function, we can help identify essential genes, new virulence factors, or new regulators of drug resistance, and thereby give new targets for antifungal development.

Functional characterization of open reading frames in nonmodel organisms, such as the common opportunistic fungal pathogen Candida albicans, can be labor-intensive. To meet this challenge, we built a comprehensive and unbiased coexpression network for C. albicans, which we call CalCEN, from data collected from 853 RNA sequencing runs from 18 large-scale studies deposited in the NCBI Sequence Read Archive. Retrospectively, CalCEN is highly predictive of known gene function annotations and can be synergistically combined with sequence similarity and interaction networks in Saccharomyces cerevisiae through orthology for additional accuracy in gene function prediction. To prospectively demonstrate the utility of the coexpression network in C. albicans, we predicted the function of underannotated open reading frames (ORFs) and identified CCJ1 as a novel cell cycle regulator in C. albicans. This study provides a tool for future systems biology analyses of gene function in C. albicans. We provide a computational pipeline for building and analyzing the coexpression network and CalCEN itself at http://github.com/momeara/CalCEN.

IMPORTANCECandida albicans is a common and deadly fungal pathogen of humans, yet the genome of this organism contains many genes of unknown function. By determining gene function, we can help identify essential genes, new virulence factors, or new regulators of drug resistance, and thereby give new targets for antifungal development. Here, we use information from large-scale RNA sequencing (RNAseq) studies and generate a C. albicans coexpression network (CalCEN) that is robust and able to predict gene function. We demonstrate the utility of this network in both retrospective and prospective testing and use CalCEN to predict a role for C4_06590W/CCJ1 in cell cycle. This tool will allow for a better characterization of underannotated genes in pathogenic yeasts.

Comparative Genome-Wide Characterization of Microsatellites in Candida albicans and Candida dubliniensis Leading to the Development of Species-Specific Marker

BACKGROUND

Microsatellites or simple sequence repeats (SSR) are related to genomic structure, function, and certain diseases of taxonomically different organisms.

OBJECTIVE

To characterize microsatellites in two closely related Candida species by searching and comparing 1-6 bp nucleotide motifs and utilizing them to develop species-specific markers.

METHODS

Whole-genome sequence was downloaded from the public domain, microsatellites were mined and analyzed, and primers were synthesized.

RESULTS

A total of 15,821 and 7,868 microsatellites, with mono-nucleotides (8,679) and trinucleotides (3,156) as most frequent microsatellites, were mined in Candida dubliniensis and Candida albicans, respectively. Chromosome size was found positively correlated with microsatellite number in both the species, whereas it was negatively correlated with the relative abundance and density of microsatellites. A number of unique motifs were also found in both the species. Overall, microsatellite frequencies of each chromosome in C. dubliniensis were higher than in C. albicans.

CONCLUSION

The features of microsatellite distribution in the two species' genomes revealed that it is probably not conserved in the genus Candida. Data generated in this article could be used for comparative genome mapping and understanding the distribution of microsatellites and genome structure between these closely related and phenotypically misidentified species and may provide a foundation for the development of a new set of species-specific microsatellite markers. Here, we also report a novel microsatellite-based marker for C. dubliniensis-specific identification.

HaploTypo: a variant-calling pipeline for phased genomes

Summary

An increasing number of phased (i.e. with resolved haplotypes) reference genomes are available. However, the most genetic variant calling tools do not explicitly account for haplotype structure. Here, we present HaploTypo, a pipeline tailored to resolve haplotypes in genetic variation analyses. HaploTypo infers the haplotype correspondence for each heterozygous variant called on a phased reference genome.

Availability and implementation

HaploTypo is implemented in Python 2.7 and Python 3.5, and is freely available at https://github.com/gabaldonlab/haplotypo, and as a Docker image.

Supplementary information

Supplementary data are available at Bioinformatics online.

Antifungal Activity of Gepinacin Scaffold Glycosylphosphatidylinositol Anchor Biosynthesis Inhibitors with Improved Metabolic Stability

The glycosylphosphatidylinositol anchor biosynthesis inhibitor gepinacin demonstrates broad-spectrum antifungal activity and negligible mammalian toxicity in culture but is metabolically labile. The stability and bioactivity of 39 analogs were tested in vitro to identify LCUT-8, a stabilized lead with increased potency and promising single-dose pharmacokinetics. Unfortunately, no antifungal activity was seen at the maximum dosing achievable in a neutropenic rabbit model. Nevertheless, structure-activity relationships identified here suggest strategies to further improve compound potency, solubility, and stability.

Expanded role of the Cu-sensing transcription factor Mac1p in Candida albicans

As part of the innate immune response, the host withholds metal micronutrients such as Cu from invading pathogens, and microbes respond through metal starvation stress responses. With the opportunistic fungal pathogen Candida albicans, the Cu-sensing transcription factor Mac1p governs the cellular response to Cu starvation by controlling Cu import. Mac1p additionally controls reactive oxygen species (ROS) homeostasis by repressing a Cu-containing superoxide dismutase (SOD1) and inducing Mn-containing SOD3 as a non-Cu alternative. We show here that C. albicans Mac1p is essential for virulence in a mouse model for disseminated candidiasis and that the cellular functions of Mac1p extend beyond Cu uptake and ROS homeostasis. Specifically, mac1∆/∆ mutants are profoundly deficient in mitochondrial respiration and Fe accumulation, both Cu-dependent processes. Surprisingly, these deficiencies are not simply the product of impaired Cu uptake; rather mac1∆/∆ mutants appear defective in Cu allocation. The respiratory defect of mac1∆/∆ mutants was greatly improved by a sod1∆/∆ mutation, demonstrating a role for SOD1 repression by Mac1p in preserving respiration. Mac1p downregulates the major Cu consumer SOD1 to spare Cu for respiration that is essential for virulence of this fungal pathogen. The implications for such Cu homeostasis control in other pathogenic fungi are discussed.

Virulence phenotypes result from interactions between pathogen ploidy and genetic background

Studying fungal virulence is often challenging and frequently depends on many contexts, including host immune status and pathogen genetic background. However, the role of ploidy has often been overlooked when studying virulence in eukaryotic pathogens. Since fungal pathogens, including the human opportunistic pathogen Candida albicans, can display extensive ploidy variation, assessing how ploidy impacts virulence has important clinical relevance. As an opportunistic pathogen, C. albicans causes nonlethal, superficial infections in healthy individuals, but life-threatening bloodstream infections in individuals with compromised immune function. Here, we determined how both ploidy and genetic background of C. albicans impacts virulence phenotypes in healthy and immunocompromised nematode hosts by characterizing virulence phenotypes in four near-isogenic diploid and tetraploid pairs of strains, which included both laboratory and clinical genetic backgrounds. We found that C. albicans infections decreased host survival and negatively impacted host reproduction, and we leveraged these two measures to survey both lethal and nonlethal virulence phenotypes across the multiple C. albicans strains. In this study, we found that regardless of pathogen ploidy or genetic background, immunocompromised hosts were susceptible to fungal infection compared to healthy hosts. Furthermore, for each host context, we found a significant interaction between C. albicans genetic background and ploidy on virulence phenotypes, but no global differences between diploid and tetraploid pathogens were observed.

The STRIPAK signaling complex regulates dephosphorylation of GUL1, an RNA-binding protein that shuttles on endosomes

The striatin-interacting phosphatase and kinase (STRIPAK) multi-subunit signaling complex is highly conserved within eukaryotes. In fungi, STRIPAK controls multicellular development, morphogenesis, pathogenicity, and cell-cell recognition, while in humans, certain diseases are related to this signaling complex. To date, phosphorylation and dephosphorylation targets of STRIPAK are still widely unknown in microbial as well as animal systems. Here, we provide an extended global proteome and phosphoproteome study using the wild type as well as STRIPAK single and double deletion mutants (Δpro11, Δpro11Δpro22, Δpp2Ac1Δpro22) from the filamentous fungus Sordaria macrospora. Notably, in the deletion mutants, we identified the differential phosphorylation of 129 proteins, of which 70 phosphorylation sites were previously unknown. Included in the list of STRIPAK targets are eight proteins with RNA recognition motifs (RRMs) including GUL1. Knockout mutants and complemented transformants clearly show that GUL1 affects hyphal growth and sexual development. To assess the role of GUL1 phosphorylation on fungal development, we constructed phospho-mimetic and -deficient mutants of GUL1 residues. While S180 was dephosphorylated in a STRIPAK-dependent manner, S216, and S1343 served as non-regulated phosphorylation sites. While the S1343 mutants were indistinguishable from wild type, phospho-deficiency of S180 and S216 resulted in a drastic reduction in hyphal growth, and phospho-deficiency of S216 also affects sexual fertility. These results thus suggest that differential phosphorylation of GUL1 regulates developmental processes such as fruiting body maturation and hyphal morphogenesis. Moreover, genetic interaction studies provide strong evidence that GUL1 is not an integral subunit of STRIPAK. Finally, fluorescence microscopy revealed that GUL1 co-localizes with endosomal marker proteins and shuttles on endosomes. Here, we provide a new mechanistic model that explains how STRIPAK-dependent and -independent phosphorylation of GUL1 regulates sexual development and asexual growth.

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

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

Cloning, functional characterization and expression analysis of LoTPS5 from Lilium 'Siberia'

Lilium 'Siberia' is a perennial herbaceous plant that is commercially significant because of its snowy white floral color and appealing scent which is mainly due to the presence of monoterpenes and benzoids compounds in floral volatile profile. In the current study, LoTPS5 was cloned and functionally characterized. Results revealed that LoTPS5 specifically generates squalene from FPP, whereas no product was produced when it was incubated with GPP or GGPP. The subcellular localization experiment showed that LoTPS5 was located in plastids. Furthermore, LoTPS5 showed its high expression in the leaf followed by petals and sepals of the flower. Moreover, the expression of LoTPS5 gradually increased from the bud stage and peak at the full-bloom stage. Besides, LoTPS5 showed a diurnal circadian rhythmic pattern with a peak in the afternoon (16:00) followed by deep night (24:00) and morning (8:00), respectively. LoTPS5 is highly responsive to mechanical wounding by rapidly elevating its mRNA transcript level. The current study will provide significant information for future studies of terpenoid and squalene biosynthesis in Lilium 'Siberia'.

Genomic evidence for a hybrid origin of the yeast opportunistic pathogen Candida albicans

Background

Opportunistic yeast pathogens of the genus Candida are an important medical problem. Candida albicans, the most prevalent Candida species, is a natural commensal of humans that can adopt a pathogenic behavior. This species is highly heterozygous and cannot undergo meiosis, adopting instead a parasexual cycle that increases genetic variability and potentially leads to advantages under stress conditions. However, the origin of C. albicans heterozygosity is unknown, and we hypothesize that it could result from ancestral hybridization. We tested this idea by analyzing available genomes of C. albicans isolates and comparing them to those of hybrid and non-hybrid strains of other Candida species.

Results

Our results show compelling evidence that C. albicans is an evolved hybrid. The genomic patterns observed in C. albicans are similar to those of other hybrids such as Candida orthopsilosis MCO456 and Candida inconspicua, suggesting that it also descends from a hybrid of two divergent lineages. Our analysis indicates that most of the divergence between haplotypes in C. albicans heterozygous blocks was already present in a putative heterozygous ancestor, with an estimated 2.8% divergence between homeologous chromosomes. The levels and patterns of ancestral heterozygosity found cannot be fully explained under the paradigm of vertical evolution and are not consistent with continuous gene flux arising from lineage-specific events of admixture.

Conclusions

Although the inferred level of sequence divergence between the putative parental lineages (2.8%) is not clearly beyond current species boundaries in Saccharomycotina, we show here that all analyzed C. albicans strains derive from a single hybrid ancestor and diverged by extensive loss of heterozygosity. This finding has important implications for our understanding of C. albicans evolution, including the loss of the sexual cycle, the origin of the association with humans, and the evolution of virulence traits.

Enhanced Efflux Pump Expression in Candida Mutants Results in Decreased Manogepix Susceptibility

Manogepix is a broad-spectrum antifungal agent that inhibits glycosylphosphatidylinositol (GPI) anchor biosynthesis. Using whole-genome sequencing, we characterized two efflux-mediated mechanisms in the fungal pathogens Candida albicans and Candida parapsilosis that resulted in decreased manogepix susceptibility. In C. albicans, a gain-of-function mutation in the transcription factor gene ZCF29 activated expression of ATP-binding cassette transporter genes CDR11 and SNQ2.

Manogepix is a broad-spectrum antifungal agent that inhibits glycosylphosphatidylinositol (GPI) anchor biosynthesis. Using whole-genome sequencing, we characterized two efflux-mediated mechanisms in the fungal pathogens Candida albicans and Candida parapsilosis that resulted in decreased manogepix susceptibility. In C. albicans, a gain-of-function mutation in the transcription factor gene ZCF29 activated expression of ATP-binding cassette transporter genes CDR11 and SNQ2. In C. parapsilosis, a mitochondrial deletion activated expression of the major facilitator superfamily transporter gene MDR1.

The Impact of Gene Dosage and Heterozygosity on The Diploid Pathobiont Candida albicans

Candida albicans is a fungal species that can colonize multiple niches in the human host where it can grow either as a commensal or as an opportunistic pathogen. The genome of C. albicans has long been of considerable interest, given that it is highly plastic and can undergo a wide variety of alterations. These changes play a fundamental role in determining C. albicans traits and have been shown to enable adaptation both to the host and to antifungal drugs. C. albicans isolates contain a heterozygous diploid genome that displays variation from the level of single nucleotides to largescale rearrangements and aneuploidy. The heterozygous nature of the genome is now increasingly recognized as being central to C. albicans biology, as the relative fitness of isolates has been shown to correlate with higher levels of overall heterozygosity. Moreover, loss of heterozygosity (LOH) events can arise frequently, either at single polymorphisms or at a chromosomal level, and both can alter the behavior of C. albicans cells during infection or can modulate drug resistance. In this review, we examine genome plasticity in this pathobiont focusing on how gene dosage variation and loss of heterozygosity events can arise and how these modulate C. albicans behavior.

The Magnitude of Candida albicans Stress-Induced Genome Instability Results from an Interaction Between Ploidy and Antifungal Drugs

Organismal ploidy and environmental stress impact the rates and types of mutational events. The opportunistic fungal pathogen Candida albicans, serves as a clinically relevant model for studying the interaction between eukaryotic ploidy and drug-induced mutagenesis. In this study, we compared the rates and types of genome perturbations in diploid and tetraploid C. albicans following exposure to two different classes of antifungal drugs; azoles and echinocandins. We measured mutations at three different scales: point mutation, loss-of-heterozygosity (LOH), and total DNA content for cells exposed to fluconazole and caspofungin. We found that caspofungin induced higher mutation rates than fluconazole, although this is likely an indirect consequence of stress-associated cell wall perturbations, rather than an inherent genotoxicity. Surprisingly, we found that antifungal drugs disproportionately elevated genome and ploidy instability in tetraploid C. albicans compared to diploids. Taken together, our results suggest that the magnitude of stress-induced mutagenesis results from an interaction between ploidy and antifungal drugs. These findings have both clinical and evolutionary implications for how fungal pathogens generate mutations in response to antifungal drug stress and how these mutations may facilitate the emergence of drug resistance.

Phased Diploid Genome Assemblies for Three Strains of Candida albicans from Oak Trees

Although normally a harmless commensal, Candida albicans, it is also one of the most common causes of bloodstream infections in the U.S. Candida albicans has long been considered an obligate commensal, however, recent studies suggest it can live outside animal hosts. Here, we have generated PacBio sequences and phased genome assemblies for three C. albicans strains from oak trees (NCYC 4144, NCYC 4145, and NCYC 4146). PacBio datasets are high depth (over 400 fold coverage) and more than half of the sequencing data are contained in reads longer than 15 kb. Primary assemblies showed high contiguity with several chromosomes for each strain recovered as single contigs, and greater than half of the alternative haplotype sequence was assembled in haplotigs at least 174 kb long. Using these assemblies we were able to identify structural polymorphisms, including a polymorphic inversion over 100 kb in length. These results show that phased de novo diploid assemblies for C. albicans can enable the study of genomic variation within and among strains of an important fungal pathogen.

To Repeat or Not to Repeat: Repetitive Sequences Regulate Genome Stability in Candida albicans

Genome instability often leads to cell death but can also give rise to innovative genotypic and phenotypic variation through mutation and structural rearrangements. Repetitive sequences and chromatin architecture in particular are critical modulators of recombination and mutability. In Candida albicans, four major classes of repeats exist in the genome: telomeres, subtelomeres, the major repeat sequence (MRS), and the ribosomal DNA (rDNA) locus. Characterization of these loci has revealed how their structure contributes to recombination and either promotes or restricts sequence evolution. The mechanisms of recombination that give rise to genome instability are known for some of these regions, whereas others are generally unexplored. More recent work has revealed additional repetitive elements, including expanded gene families and centromeric repeats that facilitate recombination and genetic innovation. Together, the repeats facilitate C. albicans evolution through construction of novel genotypes that underlie C. albicans adaptive potential and promote persistence across its human host.

Protein-Protein Interactions in Candida albicans

Despite being one of the most important human fungal pathogens, Candida albicans has not been studied extensively at the level of protein-protein interactions (PPIs) and data on PPIs are not readily available in online databases. In January 2018, the database called "Biological General Repository for Interaction Datasets (BioGRID)" that contains the most PPIs for C. albicans, only documented 188 physical or direct PPIs (release 3.4.156) while several more can be found in the literature. Other databases such as the String database, the Molecular INTeraction Database (MINT), and the Database for Interacting Proteins (DIP) database contain even fewer interactions or do not even include C. albicans as a searchable term. Because of the non-canonical codon usage of C. albicans where CUG is translated as serine rather than leucine, it is often problematic to use the yeast two-hybrid system in Saccharomyces cerevisiae to study C. albicans PPIs. However, studying PPIs is crucial to gain a thorough understanding of the function of proteins, biological processes and pathways. PPIs can also be potential drug targets. To aid in creating PPI networks and updating the BioGRID, we performed an exhaustive literature search in order to provide, in an accessible format, a more extensive list of known PPIs in C. albicans.

Chromatin Profiling of the Repetitive and Nonrepetitive Genomes of the Human Fungal Pathogen Candida albicans

The fungus Candida albicans is an opportunistic pathogen that normally lives on the human body without causing any harm. However, C. albicans is also a dangerous pathogen responsible for millions of infections annually. C. albicans is such a successful pathogen because it can adapt to and thrive in different environments. Chemical modifications of chromatin, the structure that packages DNA into cells, can allow environmental adaptation by regulating gene expression and genome organization. Surprisingly, the contribution of chromatin modification to C. albicans biology is still largely unknown. For the first time, we analyzed C. albicans chromatin modifications on a genome-wide basis. We demonstrate that specific chromatin states are associated with distinct regions of the C. albicans genome and identify the roles of the chromatin modifiers Sir2 and Set1 in shaping C. albicans chromatin and gene expression.

Eukaryotic genomes are packaged into chromatin structures that play pivotal roles in regulating all DNA-associated processes. Histone posttranslational modifications modulate chromatin structure and function, leading to rapid regulation of gene expression and genome stability, key steps in environmental adaptation. Candida albicans, a prevalent fungal pathogen in humans, can rapidly adapt and thrive in diverse host niches. The contribution of chromatin to C. albicans biology is largely unexplored. Here, we generated the first comprehensive chromatin profile of histone modifications (histone H3 trimethylated on lysine 4 [H3K4me³], histone H3 acetylated on lysine 9 [H3K9Ac], acetylated lysine 16 on histone H4 [H4K16Ac], and γH2A) across the C. albicans genome and investigated its relationship to gene expression by harnessing genome-wide sequencing approaches. We demonstrated that gene-rich nonrepetitive regions are packaged into canonical euchromatin in association with histone modifications that mirror their transcriptional activity. In contrast, repetitive regions are assembled into distinct chromatin states; subtelomeric regions and the ribosomal DNA (rDNA) locus are assembled into heterochromatin, while major repeat sequences and transposons are packaged in chromatin that bears features of euchromatin and heterochromatin. Genome-wide mapping of γH2A, a marker of genome instability, identified potential recombination-prone genomic loci. Finally, we present the first quantitative chromatin profiling in C. albicans to delineate the role of the chromatin modifiers Sir2 and Set1 in controlling chromatin structure and gene expression. This report presents the first genome-wide chromatin profiling of histone modifications associated with the C. albicans genome. These epigenomic maps provide an invaluable resource to understand the contribution of chromatin to C. albicans biology and identify aspects of C. albicans chromatin organization that differ from that of other yeasts.

Sef1-Regulated Iron Regulon Responds to Mitochondria-Dependent Iron-Sulfur Cluster Biosynthesis in Candida albicans

Iron homeostasis mechanisms allow the prime commensal-pathogen Candida albicans to cope with the profound shift in iron levels in the mammalian host. The regulators, Sef1 and Sfu1 influence activation and repression of genes required for iron uptake and acquisition by inducing the expression of iron regulon genes in iron-deplete conditions and inactivating them in iron-replete condition. Our study for the first time shows that C. albicans coordinates the activation of the iron regulon with the mitochondrial use of iron for Fe–S cluster biosynthesis, a cellular process that is connected to cellular iron metabolism. We took advantage of a mutant defective in mitochondrial biogenesis (fzo1Δ/Δ) to assess the aforesaid link as this mutant exhibited sustained expression of the Sef1 iron regulon, signifying an iron-starved state in the mutant. Our analysis demonstrates that mitochondrion is pivotal for regulation of Fe–S cluster synthesis such that the disruption of this cellular process in fzo1Δ/Δ cells lead to excessive mitochondrial iron accumulation and reduced activity of the Fe–S cluster-containing enzyme aconitase. Sef1 responds to defective Fe–S cluster synthesis by regulated changes in its subcellular localization; it was retained in the nucleus resulting in the induced expression of the iron regulon. We predict that the mitochondrial Fe–S assembly generates a molecule that is critical for ensuring iron-responsive transcriptional activation of the Sef1 regulon. All told, our data marks Fe–S biogenesis as a mechanism that meshes cellular iron procurement with mitochondrial iron metabolism resulting in regulating the Sef1 regulon in C. albicans.

ParaDB: A manually curated database containing genomic annotation for the human pathogenic fungi Paracoccidioides spp

BACKGROUND

The genus Paracoccidioides consists of thermodymorphic fungi responsible for Paracoccidioidomycosis (PCM), a systemic mycosis that has been registered to affect ~10 million people in Latin America. Biogeographical data subdivided the genus Paracoccidioides in five divergent subgroups, which have been recently classified as different species. Genomic sequencing of five Paracoccidioides isolates, representing each of these subgroups/species provided an important framework for the development of post-genomic studies with these fungi. However, functional annotations of these genomes have not been submitted to manual curation and, as a result, ~60-90% of the Paracoccidioides protein-coding genes (depending on isolate/annotation) are currently described as responsible for hypothetical proteins, without any further functional/structural description.

PRINCIPAL FINDINGS

The present work reviews the functional assignment of Paracoccidioides genes, reducing the number of hypothetical proteins to ~25-28%. These results were compiled in a relational database called ParaDB, dedicated to the main representatives of Paracoccidioides spp. ParaDB can be accessed through a friendly graphical interface, which offers search tools based on keywords or protein/DNA sequences. All data contained in ParaDB can be partially or completely downloaded through spreadsheet, multi-fasta and GFF3-formatted files, which can be subsequently used in a variety of downstream functional analyses. Moreover, the entire ParaDB environment has been configured in a Docker service, which has been submitted to the GitHub repository, ensuring long-term data availability to researchers. This service can be downloaded and used to perform fully functional local installations of the database in alternative computing ecosystems, allowing users to conduct their data mining and analyses in a personal and stable working environment.

CONCLUSIONS

These new annotations greatly reduce the number of genes identified solely as hypothetical proteins and are integrated into a dedicated database, providing resources to assist researchers in this field to conduct post-genomic studies with this group of human pathogenic fungi.

Mechanisms of genome evolution in Candida albicans

The fungus Candida albicans exists as a prevalent commensal and an important opportunistic pathogen that can infect multiple niches of its human host. Recent studies have examined the diploid genome of C. albicans by performing both short-term microevolution studies and comparative genomics on collections of clinical isolates. Common mechanisms driving genome dynamics include accumulation of point mutations, loss of heterozygosity (LOH) events, large-scale chromosomal rearrangements, and even ploidy change, with important consequences for both drug resistance and host adaptation. Evidence for recombination between C. albicans lineages also highlights a role for (para)sex in shaping the species population structure. Ongoing work will continue to define the contributions of genome evolution to phenotypic variation and the role of host pressures in driving adaptive processes.

High throughput gene expression profiling of yeast colonies with microgel-culture Drop-seq

Yeast can be engineered into "living foundries" for non-natural chemical production by reprogramming them via a "design-build-test" cycle. While methods for "design" and "build" are relatively scalable and efficient, "test" remains a bottleneck, limiting the effectiveness of the procedure. Here we describe isogenic colony sequencing (ICO-seq), a massively-parallel strategy to assess the gene expression, and thus engineered pathway efficacy, of large numbers of genetically distinct yeast colonies. We use the approach to characterize opaque-white switching in 658 C. albicans colonies. By profiling the transcriptomes of 1642 engineered S. cerevisiae strains, we assess gene expression heterogeneity in a protein mutagenesis library. Our approach will accelerate synthetic biology by allowing facile and cost-effective transcriptional profiling of large numbers of genetically distinct yeast strains.

Candida albicans: An Emerging Yeast Model to Study Eukaryotic Genome Plasticity

Saccharomyces cerevisiae and Schizosaccharomyces pombe have served as uncontested unicellular model organisms, as major discoveries made in the field of genome biology using yeast genetics have proved to be relevant from yeast to humans. The yeast Candida albicans has attracted much attention because of its ability to switch between a harmless commensal and a dreaded human pathogen. C. albicans bears unique features regarding its life cycle, genome structure, and dynamics, and their links to cell biology and adaptation to environmental challenges. Examples include a unique reproduction cycle with haploid, diploid, and tetraploid forms; a distinctive organisation of chromosome hallmarks; a highly dynamic genome, with extensive karyotypic variations, including aneuploidies, isochromosome formation, and loss-of-heterozygosity; and distinctive links between the response to DNA alterations and cell morphology. These features have made C. albicans emerge as a new and attractive unicellular model to study genome biology and dynamics in eukaryotes.

Diversity of cytosine methylation across the fungal tree of life

The generation of thousands of fungal genomes is leading to a better understanding of genes and genomic organization within the kingdom. However, the epigenome, which includes DNA and chromatin modifications, remains poorly investigated in fungi. Large comparative studies in animals and plants have deepened our understanding of epigenomic variation, particularly of the modified base 5-methylcytosine (5mC), but taxonomic sampling of disparate groups is needed to develop unifying explanations for 5mC variation. Here we utilize the largest phylogenetic resolution of 5mC methyltransferases (5mC MTases) and genome evolution to better understand levels and patterns of 5mC across fungi. We show that extant 5mC MTase genotypes are descendent from ancestral maintenance and de novo genotypes, whereas the 5mC MTases DIM-2 and RID are more recently derived, and that 5mC levels are correlated with 5mC MTase genotype and transposon content. Our survey also revealed that fungi lack canonical gene body methylation, which distinguishes fungal epigenomes from certain insect and plant species. However, some fungal species possess independently derived clusters of contiguous 5mC encompassing many genes. In some cases, DNA repair pathways and the N6-methyladenine (6mA) DNA modification negatively coevolved with 5mC pathways, which additionally contributed to interspecific epigenomic variation across fungi.

Identification of Recessive Lethal Alleles in the Diploid Genome of a Candida albicans Laboratory Strain Unveils a Potential Role of Repetitive Sequences in Buffering Their Deleterious Impact

The heterozygous diploid genome of Candida albicans is highly plastic, with frequent loss of heterozygosity (LOH) events. In the SC5314 laboratory strain, while LOH events are ubiquitous, a chromosome homozygosis bias is observed for certain chromosomes, whereby only one of the two homologs can occur in the homozygous state. This suggests the occurrence of recessive lethal allele(s) (RLA) preventing large-scale LOH events on these chromosomes from being stably maintained. To verify the presence of an RLA on chromosome 7 (Chr7), we utilized a system that allows (i) DNA double-strand break (DSB) induction on Chr7 by the I-SceI endonuclease and (ii) detection of the resulting long-range homozygosis. I-SceI successfully induced a DNA DSB on both Chr7 homologs, generally repaired by gene conversion. Notably, cells homozygous for the right arm of Chr7B were not recovered, confirming the presence of RLA(s) in this region. Genome data mining for RLA candidates identified a premature nonsense-generating single nucleotide polymorphism (SNP) within the HapB allele of C7_03400c whose Saccharomyces cerevisiae ortholog encodes the essential Mtr4 RNA helicase. Complementation with a wild-type copy of MTR4 rescued cells homozygous for the right arm of Chr7B, demonstrating that the mtr4^K880* RLA is responsible for the Chr7 homozygosis bias in strain SC5314. Furthermore, we observed that the major repeat sequences (MRS) on Chr7 acted as hot spots for interhomolog recombination. Such recombination events provide C. albicans with increased opportunities to survive DNA DSBs whose repair can lead to homozygosis of recessive lethal or deleterious alleles. This might explain the maintenance of MRS in this species.IMPORTANCE Candida albicans is a major fungal pathogen, whose mode of reproduction is mainly clonal. Its genome is highly tolerant to rearrangements, in particular loss of heterozygosity events, known to unmask recessive lethal and deleterious alleles in heterozygous diploid organisms such as C. albicans By combining a site-specific DSB-inducing system and mining genome sequencing data of 182 C. albicans isolates, we were able to ascribe the chromosome 7 homozygosis bias of the C. albicans laboratory strain SC5314 to an heterozygous SNP introducing a premature STOP codon in the MTR4 gene. We have also proposed genome-wide candidates for new recessive lethal alleles. We additionally observed that the major repeat sequences (MRS) on chromosome 7 acted as hot spots for interhomolog recombination. Maintaining MRS in C. albicans could favor haplotype exchange, of vital importance to LOH events, leading to homozygosis of recessive lethal or deleterious alleles that inevitably accumulate upon clonality.

Diverse Lineages of Candida albicans Live on Old Oaks

The human pathogen Candida albicans is considered an obligate commensal of animals, yet it is occasionally isolated from trees, shrubs, and grass. We generated genome sequence data for three strains of C. albicans that we isolated from oak trees in an ancient wood pasture, and compared these to the genomes of over 200 clinical strains. C. albicans strains from oak are similar to clinical C. albicans in that they are predominantly diploid and can become homozygous at the mating locus through whole-chromosome loss of heterozygosity. Oak strains differed from clinical strains in showing slightly higher levels of heterozygosity genome-wide. Using phylogenomic analyses and in silico chromosome painting, we show that each oak strain is more closely related to strains from humans and other animals than to strains from other oaks. The high genetic diversity of C. albicans from old oaks shows that they can live in this environment for extended periods of time.

Squalene synthase cloning and functional identification in wintersweet plant (Chimonanthus zhejiangensis)

BACKGROUND

Three species of wintersweets: Chimonanthus salicifolius S. Y. Hu, Chimonanthus zhejiangensis M. C. Liu and Chimonanthus grammalus M. C. Liu are widely distributed in China. The three wintersweets belonging to the genus of Chimonanthus that can synthesize abundant terpenoids that are beneficial to human health. Their buds and leaves are traditional Chinese herb applied by the 'She' ethnic minority in southeast of China. Squalene is a multi-functional and ubiquitous triterpene in plants, which is biosynthesized by squalene synthase (SQS) using farnesyl diphosphate (FPP) as the substrate. The synthesis of squalene in wintersweet was not clearly. This work would provide us much help to further understand the terpene metabolism in wintersweet and its health function to people at phytochemistry and molecular levels.

RESULTS

In this study, we identified squalene component in the extractions of leaves of three wintersweets and isolated SQS genes from leaf transcriptomes. The three SQSs were highly conservative, so CzSQS from C. zhejiangensis was just determined the enzymatic activity. The in vitro expressed CzSQS that deleted two transmembrane domains could catalyze FPP to generate squalene with the presence of NADPH and Mg²⁺.

CONCLUSIONS

The squalene was one of wintersweet leaves phytochemicals. The squalene synthases of three wintersweet plants were highly conserved. The CzSQS was capable to catalyze two FPP molecules to squalene.