Transcriptome and proteome profiling reveals complex adaptations of Candida parapsilosis cells assimilating hydroxyaromatic carbon sources.

Many fungal species utilize hydroxyderivatives of benzene and benzoic acid as carbon sources. The yeast Candida parapsilosis metabolizes these compounds via the 3-oxoadipate and gentisate pathways, whose components are encoded by two metabolic gene clusters. In this study, we determine the chromosome level assembly of the C. parapsilosis strain CLIB214 and use it for transcriptomic and proteomic investigation of cells cultivated on hydroxyaromatic substrates. We demonstrate that the genes coding for enzymes and plasma membrane transporters involved in the 3-oxoadipate and gentisate pathways are highly upregulated and their expression is controlled in a substrate-specific manner. However, regulatory proteins involved in this process are not known. Using the knockout mutants, we show that putative transcriptional factors encoded by the genes OTF1 and GTF1 located within these gene clusters function as transcriptional activators of the 3-oxoadipate and gentisate pathway, respectively. We also show that the activation of both pathways is accompanied by upregulation of genes for the enzymes involved in ß-oxidation of fatty acids, glyoxylate cycle, amino acid metabolism, and peroxisome biogenesis. Transcriptome and proteome profiles of the cells grown on 4-hydroxybenzoate and 3-hydroxybenzoate, which are metabolized via the 3-oxoadipate and gentisate pathway, respectively, reflect their different connection to central metabolism. Yet we find that the expression profiles differ also in the cells assimilating 4-hydroxybenzoate and hydroquinone, which are both metabolized in the same pathway. This finding is consistent with the phenotype of the Otf1p-lacking mutant, which exhibits impaired growth on hydroxybenzoates, but still utilizes hydroxybenzenes, thus indicating that additional, yet unidentified transcription factor could be involved in the 3-oxoadipate pathway regulation. Moreover, we propose that bicarbonate ions resulting from decarboxylation of hydroxybenzoates also contribute to differences in the cell responses to hydroxybenzoates and hydroxybenzenes. Finally, our phylogenetic analysis highlights evolutionary paths leading to metabolic adaptations of yeast cells assimilating hydroxyaromatic substrates.

The yeast mitochondrial succinylome: Implications for regulation of mitochondrial nucleoids.

Acylation modifications, such as the succinylation of lysine, are post-translational modifications and a powerful means of regulating protein activity. Some acylations occur nonenzymatically, driven by an increase in the concentration of acyl group donors. Lysine succinylation has a profound effect on the corresponding site within the protein, as it dramatically changes the charge of the residue. In eukaryotes, it predominantly affects mitochondrial proteins because the donor of succinate, succinyl-CoA, is primarily generated in the tricarboxylic acid cycle. Although numerous succinylated mitochondrial proteins have been identified in Saccharomyces cerevisiae, a more detailed characterization of the yeast mitochondrial succinylome is still lacking. Here, we performed a proteomic MS analysis of purified yeast mitochondria and detected 314 succinylated mitochondrial proteins with 1763 novel succinylation sites. The mitochondrial nucleoid, a complex of mitochondrial DNA and mitochondrial proteins, is one of the structures whose protein components are affected by succinylation. We found that Abf2p, the principal component of mitochondrial nucleoids responsible for compacting mitochondrial DNA in S. cerevisiae, can be succinylated in vivo on at least thirteen lysine residues. Abf2p succinylation in vitro inhibits its DNA-binding activity and reduces its sensitivity to digestion by the ATP-dependent ScLon protease. We conclude that changes in the metabolic state of a cell resulting in an increase in the concentration of tricarboxylic acid intermediates may affect mitochondrial functions.

Role of folding kinetics of secondary structures in telomeric G-overhangs in the regulation of telomere maintenance in Saccharomyces cerevisiae.

The ends of eukaryotic chromosomes typically contain a 3' ssDNA G-rich protrusion (G-overhang). This overhang must be protected against detrimental activities of nucleases and of the DNA damage response machinery and participates in the regulation of telomerase, a ribonucleoprotein complex that maintains telomere integrity. These functions are mediated by DNA-binding proteins, such as Cdc13 in Saccharomyces cerevisiae, and the propensity of G-rich sequences to form various non-B DNA structures. Using CD and NMR spectroscopies, we show here that G-overhangs of S. cerevisiae form distinct Hoogsteen pairing-based secondary structures, depending on their length. Whereas short telomeric oligonucleotides form a G-hairpin, their longer counterparts form parallel and/or antiparallel G-quadruplexes (G4s). Regardless of their topologies, non-B DNA structures exhibited impaired binding to Cdc13 in vitro as demonstrated by electrophoretic mobility shift assays. Importantly, whereas G4 structures formed relatively quickly, G-hairpins folded extremely slowly, indicating that short G-overhangs, which are typical for most of the cell cycle, are present predominantly as single-stranded oligonucleotides and are suitable substrates for Cdc13. Using ChIP, we show that the occurrence of G4 structures peaks at the late S phase, thus correlating with the accumulation of long G-overhangs. We present a model of how time- and length-dependent formation of non-B DNA structures at chromosomal termini participates in telomere maintenance.

OCT1 - a yeast mitochondrial thiolase involved in the 3-oxoadipate pathway.

The 3-oxoacyl-CoA thiolases catalyze the last step of the fatty acid ß-oxidation pathway. In yeasts and plants, this pathway takes place exclusively in peroxisomes, whereas in animals it occurs in both peroxisomes and mitochondria. In contrast to baker's yeast Saccharomyces cerevisiae, yeast species from the Debaryomycetaceae family also encode a thiolase with predicted mitochondrial localization. These yeasts are able to utilize a range of hydroxyaromatic compounds via the 3-oxoadipate pathway the last step of which is catalyzed by 3-oxoadipyl-CoA thiolase and presumably occurs in mitochondria. In this work, we studied Oct1p, an ortholog of this enzyme from Candida parapsilosis. We found that the cells grown on a 3-oxoadipate pathway substrate exhibit increased levels of the OCT1 mRNA. Deletion of both OCT1 alleles impairs the growth of C. parapsilosis cells on 3-oxoadipate pathway substrates and this defect can be rescued by expression of the OCT1 gene from a plasmid vector. Subcellular localization experiments and LC-MS/MS analysis of enriched organellar fraction-proteins confirmed the presence of Oct1p in mitochondria. Phylogenetic profiling of Oct1p revealed an intricate evolutionary pattern indicating multiple horizontal gene transfers among different fungal groups.

A SARS-CoV-2 mutant from B.1.258 lineage with ∆H69/∆V70 deletion in the Spike protein circulating in Central Europe in the fall 2020.

SARS-CoV-2 mutants carrying the ∆H69/∆V70 deletion in the amino-terminal domain of the Spike protein emerged independently in at least six lineages of the virus (namely, B.1.1.7, B.1.1.298, B.1.160, B.1.177, B.1.258, B.1.375). We analyzed SARS-CoV-2 samples collected from various regions of Slovakia between November and December 2020 that were presumed to contain B.1.1.7 variant due to drop-out of the Spike gene target in an RT-qPCR test caused by this deletion. Sequencing of these samples revealed that although in some cases the samples were indeed confirmed as B.1.1.7, a substantial fraction of samples contained another ∆H69/∆V70 carrying mutant belonging to the lineage B.1.258, which has been circulating in Central Europe since August 2020, long before the import of B.1.1.7. Phylogenetic analysis shows that the early sublineage of B.1.258 acquired the N439K substitution in the receptor-binding domain (RBD) of the Spike protein and, later on, also the deletion ∆H69/∆V70 in the Spike N-terminal domain (NTD). This variant was particularly common in several European countries including the Czech Republic and Slovakia but has been quickly replaced by B.1.1.7 early in 2021.

Candida parapsilosis: from Genes to the Bedside.

Patients with suppressed immunity are at the highest risk for hospital-acquired infections. Among these, invasive candidiasis is the most prevalent systemic fungal nosocomial infection. Over recent decades, the combined prevalence of non-albicans Candida species outranked Candida albicans infections in several geographical regions worldwide, highlighting the need to understand their pathobiology in order to develop effective treatment and to prevent future outbreaks. Candida parapsilosis is the second or third most frequently isolated Candida species from patients. Besides being highly prevalent, its biology differs markedly from that of C. albicans, which may be associated with C. parapsilosis' increased incidence. Differences in virulence, regulatory and antifungal drug resistance mechanisms, and the patient groups at risk indicate that conclusions drawn from C. albicans pathobiology cannot be simply extrapolated to C. parapsilosis Such species-specific characteristics may also influence their recognition and elimination by the host and the efficacy of antifungal drugs. Due to the availability of high-throughput, state-of-the-art experimental tools and molecular genetic methods adapted to C. parapsilosis, genome and transcriptome studies are now available that greatly contribute to our understanding of what makes this species a threat. In this review, we summarize 10 years of findings on C. parapsilosis pathogenesis, including the species' genetic properties, transcriptome studies, host responses, and molecular mechanisms of virulence. Antifungal susceptibility studies and clinician perspectives are discussed. We also present regional incidence reports in order to provide an updated worldwide epidemiology summary.

Co-evolution in the Jungle: From Leafcutter Ant Colonies to Chromosomal Ends.

Biological entities are multicomponent systems where each part is directly or indirectly dependent on the others. In effect, a change in a single component might have a consequence on the functioning of its partners, thus affecting the fitness of the entire system. In this article, we provide a few examples of such complex biological systems, ranging from ant colonies to a population of amino acids within a single-polypeptide chain. Based on these examples, we discuss one of the central and still challenging questions in biology: how do such multicomponent consortia co-evolve? More specifically, we ask how telomeres, nucleo-protein complexes protecting the integrity of linear DNA chromosomes, originated from the ancestral organisms having circular genomes and thus not dealing with end-replication and end-protection problems. Using the examples of rapidly evolving topologies of mitochondrial genomes in eukaryotic microorganisms, we show what means of co-evolution were employed to accommodate various types of telomere-maintenance mechanisms in mitochondria. We also describe an unprecedented runaway evolution of telomeric repeats in nuclei of ascomycetous yeasts accompanied by co-evolution of telomere-associated proteins. We propose several scenarios derived from research on telomeres and supported by other studies from various fields of biology, while emphasizing that the relevant answers are still not in sight. It is this uncertainty and a lack of a detailed roadmap that makes the journey through the jungle of biological systems still exciting and worth undertaking.

Genome sequence of the opportunistic human pathogen Magnusiomyces capitatus.

The yeast Magnusiomyces capitatus is an opportunistic human pathogen causing rare yet severe infections, especially in patients with hematological malignancies. Here, we report the 20.2 megabase genome sequence of an environmental strain of this species as well as the genome sequences of eight additional isolates from human and animal sources providing an insight into intraspecies variation. The distribution of single-nucleotide variants is indicative of genetic recombination events, supporting evidence for sexual reproduction in this heterothallic yeast. Using RNAseq-aided annotation, we identified genes for 6518 proteins including several expanded families such as kexin proteases and Hsp70 molecular chaperones. Several of these families are potentially associated with the ability of M. capitatus to infect and colonize humans. For the purpose of comparative analysis, we also determined the genome sequence of a closely related yeast, Magnusiomyces ingens. The genome sequences of M. capitatus and M. ingens exhibit many distinct features and represent a basis for further comparative and functional studies.

Correction: The Genomic Aftermath of Hybridization in the Opportunistic Pathogen Candida metapsilosis.

[This corrects the article DOI: 10.1371/journal.pgen.1005626.].

Programmed translational bypassing elements in mitochondria: structure, mobility, and evolutionary origin.

Programmed translational bypassing enables ribosomes to 'ignore' a precise mRNA interval of several dozen nucleotides. Well-characterized bypassed sequences include hop and byp elements, present in bacteriophage T4 and mitochondria of the yeast Magnusiomyces capitatus, respectively. The bypassing mechanism of byps is probably similar to that of hop, yet the former appears more effective and less constrained as to sequence context. Furthermore, both elements are mobile but hop moves as part of a cassette including a homing endonuclease, whereas byps seem to spread like miniature DNA transposable elements known as GC clusters. Here, we argue that hop and byps arose independently by convergent evolution, and that byps evolved in magnusiomycete mitochondria due to (as yet unknown) alterations of the mitochondrial translation machinery.

The structure and DNA-binding properties of Mgm101 from a yeast with a linear mitochondrial genome.

To study the mechanisms involved in the maintenance of a linear mitochondrial genome we investigated the biochemical properties of the recombination protein Mgm101 from Candida parapsilosis. We show that CpMgm101 complements defects associated with the Saccharomyces cerevisiae mgm101-1(ts) mutation and that it is present in both the nucleus and mitochondrial nucleoids of C. parapsilosis. Unlike its S. cerevisiae counterpart, CpMgm101 is associated with the entire nucleoid population and is able to bind to a broad range of DNA substrates in a non-sequence specific manner. CpMgm101 is also able to catalyze strand annealing and D-loop formation. CpMgm101 forms a roughly C-shaped trimer in solution according to SAXS. Electron microscopy of a complex of CpMgm101 with a model mitochondrial telomere revealed homogeneous, ring-shaped structures at the telomeric single-stranded overhangs. The DNA-binding properties of CpMgm101, together with its DNA recombination properties, suggest that it can play a number of possible roles in the replication of the mitochondrial genome and the maintenance of its telomeres.

The Genomic Aftermath of Hybridization in the Opportunistic Pathogen Candida metapsilosis.

Candida metapsilosis is a rarely-isolated, opportunistic pathogen that belongs to a clade of pathogenic yeasts known as the C. parapsilosis sensu lato species complex. To gain insight into the recent evolution of C. metapsilosis and the genetic basis of its virulence, we sequenced the genome of 11 clinical isolates from various locations, which we compared to each other and to the available genomes of the two remaining members of the complex: C. orthopsilosis and C. parapsilosis. Unexpectedly, we found compelling genomic evidence that C. metapsilosis is a highly heterozygous hybrid species, with all sequenced clinical strains resulting from the same past hybridization event involving two parental lineages that were approximately 4.5% divergent in sequence. This result indicates that the parental species are non-pathogenic, but that hybridization between them formed a new opportunistic pathogen, C. metapsilosis, that has achieved a worldwide distribution. We show that these hybrids are diploid and we identified strains carrying loci for both alternative mating types, which supports mating as the initial mechanism for hybrid formation. We trace the aftermath of this hybridization at the genomic level, and reconstruct the evolutionary relationships among the different strains. Recombination and introgression -resulting in loss of heterozygosis- between the two subgenomes have been rampant, and includes the partial overwriting of the MTLa mating locus in all strains. Collectively, our results shed light on the recent genomic evolution within the C. parapsilosis sensu lato complex, and argue for a re-definition of species within this clade, with at least five distinct homozygous lineages, some of which having the ability to form hybrids.

Massive programmed translational jumping in mitochondria.

Programmed translational bypassing is a process whereby ribosomes "ignore" a substantial interval of mRNA sequence. Although discovered 25 y ago, the only experimentally confirmed example of this puzzling phenomenon is expression of the bacteriophage T4 gene 60. Bypassing requires translational blockage at a "takeoff codon" immediately upstream of a stop codon followed by a hairpin, which causes peptidyl-tRNA dissociation and reassociation with a matching "landing triplet" 50 nt downstream, where translation resumes. Here, we report 81 translational bypassing elements (byps) in mitochondria of the yeast Magnusiomyces capitatus and demonstrate in three cases, by transcript analysis and proteomics, that byps are retained in mitochondrial mRNAs but not translated. Although mitochondrial byps resemble the bypass sequence in the T4 gene 60, they utilize unused codons instead of stops for translational blockage and have relaxed matching rules for takeoff/landing sites. We detected byp-like sequences also in mtDNAs of several Saccharomycetales, indicating that byps are mobile genetic elements. These byp-like sequences lack bypassing activity and are tolerated when inserted in-frame in variable protein regions. We hypothesize that byp-like elements have the potential to contribute to evolutionary diversification of proteins by adding new domains that allow exploration of new structures and functions.

Replication intermediates of the linear mitochondrial DNA of Candida parapsilosis suggest a common recombination based mechanism for yeast mitochondria.

Variation in the topology of mitochondrial DNA (mtDNA) in eukaryotes evokes the question if differently structured DNAs are replicated by a common mechanism. RNA-primed DNA synthesis has been established as a mechanism for replicating the circular animal/mammalian mtDNA. In yeasts, circular mtDNA molecules were assumed to be templates for rolling circle DNA-replication. We recently showed that in Candida albicans, which has circular mapping mtDNA, recombination driven replication is a major mechanism for replicating a complex branched mtDNA network. Careful analyses of C. albicans-mtDNA did not reveal detectable amounts of circular DNA molecules. In the present study we addressed the question of how the unit sized linear mtDNA of Candida parapsilosis terminating at both ends with arrays of tandem repeats (mitochondrial telomeres) is replicated. Originally, we expected to find replication intermediates diagnostic of canonical bi-directional replication initiation at the centrally located bi-directional promoter region. However, we found that the linear mtDNA of Candida parapsilosis also employs recombination for replication initiation. The most striking findings were that the mitochondrial telomeres appear to be hot spots for recombination driven replication, and that stable RNA:DNA hybrids, with a potential role in mtDNA replication, are also present in the mtDNA preparations.

Saccharomyces cerevisiae as a model for the study of extranuclear functions of mammalian telomerase.

The experimental evidence from the last decade made telomerase a prominent member of a family of moonlighting proteins performing different functions at various cellular loci. However, the study of extratelomeric functions of the catalytic subunit of mammalian telomerase (TERT) is often complicated by the fact that it is sometimes difficult to distinguish them from its role(s) at the chromosomal ends. Here, we present an experimental model for studying the extranuclear function(s) of mammalian telomerase in the yeast Saccharomyces cerevisiae. We demonstrate that the catalytic subunit of mammalian telomerase protects the yeast cells against oxidative stress and affects the stability of the mitochondrial genome. The advantage of using S. cerevisiae to study of mammalian telomerase is that (1) mammalian TERT does not interfere with its yeast counterpart in the maintenance of telomeres, (2) yeast telomerase is not localized in mitochondria and (3) it does not seem to be involved in the protection of cells against oxidative stress and stabilization of mtDNA. Thus, yeast cells can be used as a 'test tube' for reconstitution of mammalian TERT extranuclear function(s).

Mdm31 protein mediates sensitivity to potassium ionophores but does not regulate mitochondrial morphology or phospholipid trafficking in Schizosaccharomyces pombe.

Mdm31p is an inner mitochondrial membrane (IMM) protein with unknown function in Saccharomyces cerevisiae. Mutants lacking Mdm31p contain only a few giant spherical mitochondria with disorganized internal structure, altered phospholipid composition and disturbed ion homeostasis, accompanied by increased resistance to the electroneutral K+ /H+ ionophore nigericin. These phenotypes are interpreted as resulting from diverse roles of Mdm31p, presumably in linking mitochondrial DNA (mtDNA) to the machinery involved in segregation of mitochondria, in mediating cation transport across IMM and in phospholipid shuttling between mitochondrial membranes. To investigate which of the roles of Mdm31p are conserved in ascomycetous yeasts, we analysed the Mdm31p orthologue in Schizosaccharomyces pombe. Our results demonstrate that, similarly to its S. cerevisiae counterpart, SpMdm31 is a mitochondrial protein and its absence results in increased resistance to nigericin. However, in contrast to S. cerevisiae, Sz. pombe cells lacking SpMdm31 are also less sensitive to the electrogenic K+ ionophore valinomycin. Moreover, mitochondria of the fission yeast mdm31Δ mutant display no changes in morphology or phospholipid composition. Therefore, in terms of function, the two orthologous proteins appear to have considerably diverged between these two evolutionarily distant yeast species, possibly sharing only their participation in ion homeostasis.

Metabolic gene clusters encoding the enzymes of two branches of the 3-oxoadipate pathway in the pathogenic yeast Candida albicans.

The pathogenic yeast Candida albicans utilizes hydroxyderivatives of benzene via the catechol and hydroxyhydroquinone branches of the 3-oxoadipate pathway. The genetic basis and evolutionary origin of this catabolic pathway in yeasts are unknown. In this study, we identified C. albicans genes encoding the enzymes involved in the degradation of hydroxybenzenes. We found that the genes coding for core components of the 3-oxoadipate pathway are arranged into two metabolic gene clusters. Our results demonstrate that C. albicans cells cultivated in media containing hydroxybenzene substrates highly induce the transcription of these genes as well as the corresponding enzymatic activities. We also found that C. albicans cells assimilating hydroxybenzenes cope with the oxidative stress by upregulation of cellular antioxidant systems such as alternative oxidase and catalase. Moreover, we investigated the evolution of the enzymes encoded by these clusters and found that most of them share a particularly sparse phylogenetic distribution among Saccharomycotina, which is likely to have been caused by extensive gene loss. We exploited this fact to find co-evolving proteins that are suitable candidates for the missing enzymes of the pathway.

The strictly aerobic yeast Yarrowia lipolytica tolerates loss of a mitochondrial DNA-packaging protein.

Mitochondrial DNA (mtDNA) is highly compacted into DNA-protein structures termed mitochondrial nucleoids (mt-nucleoids). The key mt-nucleoid components responsible for mtDNA condensation are HMG box-containing proteins such as mammalian mitochondrial transcription factor A (TFAM) and Abf2p of the yeast Saccharomyces cerevisiae. To gain insight into the function and organization of mt-nucleoids in strictly aerobic organisms, we initiated studies of these DNA-protein structures in Yarrowia lipolytica. We identified a principal component of mt-nucleoids in this yeast and termed it YlMhb1p (Y. lipolytica mitochondrial HMG box-containing protein 1). YlMhb1p contains two putative HMG boxes contributing both to DNA binding and to its ability to compact mtDNA in vitro. Phenotypic analysis of a Δmhb1 strain lacking YlMhb1p resulted in three interesting findings. First, although the mutant exhibits clear differences in mt-nucleoids accompanied by a large decrease in the mtDNA copy number and the number of mtDNA-derived transcripts, its respiratory characteristics and growth under most of the conditions tested are indistinguishable from those of the wild-type strain. Second, our results indicate that a potential imbalance between subunits of the respiratory chain encoded separately by nuclear DNA and mtDNA is prevented at a (post)translational level. Third, we found that mtDNA in the Δmhb1 strain is more prone to mutations, indicating that mtHMG box-containing proteins protect the mitochondrial genome against mutagenic events.

Mitochondrial genome of the basidiomycetous yeast Jaminaea angkorensis.

Jaminaea angkorensis is an anamorphic basidiomycetous yeast species originally isolated from decaying leaves in Cambodia. Taxonomically, J. angkorensis is affiliated with Microstromatales (Exobasidiomycetes, Ustilaginomycotina, Basidiomycota) and represents a basal phylogenetic lineage of this fungal order. To perform a comparative analysis of J. angkorensis with other basidiomycetes, we determined and analyzed its complete mitochondrial DNA sequence. The mitochondrial genome is represented by 29,999 base pairs long, circular DNA containing 32 % guanine and cytosine residues. Its genetic organization is relatively compact and comprises typical genes for 15 conserved proteins involved in oxidative phosphorylation (atp6, 8, and 9; cob; cox1, 2, and 3; and nad1, 2, 3, 4, 4L, 5, and 6) and translation (rps3), two ribosomal RNAs (rnl and rns) and twenty-two transfer RNAs (trnA-Y). Although the gene content is similar to other basidiomycetes, the gene orders in the examined species exhibit only a limited synteny, reflecting their phylogenetic distances and extensive genome rearrangements. In addition, a comparative analysis of basidiomycete mitochondrial genomes indicates that stop-to-tryptophan reassignment of the UGA codon was accompanied by structural alterations of tRNA-Trp(CCA). These results provide an insight into the evolution of the genetic code in fungal mitochondria.