The transcriptome of Candida albicans mitochondria and the evolution of organellar transcription units in yeasts

The Identification of the Mitochondrial DNA Polymerase γ (Mip1) of the Entomopathogenic Fungus Metarhizium brunneum

Replication of the mitochondrial (mt) genome in filamentous fungi is under-studied, and knowledge is based mainly on data from yeasts and higher eukaryotes. In this study, the mitochondrial DNA polymerase γ (Mip1) of the entomopathogenic fungus Metarhizium brunneum is characterized and analyzed with disruption experiments and its in silico interactions with key proteins implicated in mt gene transcription, i.e., mt RNA polymerase Rpo41 and mt transcription factor Mtf1. Disruption of mip1 gene and its partial expression influences cell growth, morphology, germination and stress tolerance. A putative in silico model of Mip1-Rpo41-Mtf1, which is known to be needed for the initiation of replication, was proposed and helped to identify potential amino acid residues of Mip1 that interact with the Rpo41-Mtf1 complex. Moreover, the reduced expression of mip1 indicates that Mip1 is not required for efficient transcription but only for replication. Functional differences between the M. brunneum Mip1 and its counterparts from Saccharomyces cerevisiae and higher eukaryotes are discussed.

RNA processing and degradation mechanisms shaping the mitochondrial transcriptome of budding yeasts

Yeast mitochondrial genes are expressed as polycistronic transcription units that contain RNAs from different classes and show great evolutionary variability. The promoters are simple, and transcriptional control is rudimentary. Posttranscriptional mechanisms involving RNA maturation, stability, and degradation are thus the main force shaping the transcriptome and determining the expression levels of individual genes. Primary transcripts are fragmented by tRNA excision by RNase P and tRNase Z, additional processing events occur at the dodecamer site at the 3' end of protein-coding sequences. groups I and II introns are excised in a self-splicing reaction that is supported by protein splicing factors encoded by the nuclear genes, or by the introns themselves. The 3'-to-5' exoribonucleolytic complex called mtEXO is the main RNA degradation activity involved in RNA turnover and processing, supported by an auxiliary 5'-to-3' exoribonuclease Pet127p. tRNAs and, to a lesser extent, rRNAs undergo several different base modifications. This complex gene expression system relies on the coordinated action of mitochondrial and nuclear genes and undergoes rapid evolution, contributing to speciation events. Moving beyond the classical model yeast Saccharomyces cerevisiae to other budding yeasts should provide important insights into the coevolution of both genomes that constitute the eukaryotic genetic system.

Respiration supports intraphagosomal filamentation and escape of Candida albicans from macrophages

IMPORTANCE

Candida albicans is a leading human fungal pathogen that often causes life-threatening infections in immunocompromised individuals. The ability of C. albicans to transition between yeast and filamentous forms is key to its virulence, and this occurs in response to many host-relevant cues, including engulfment by host macrophages. While previous efforts identified C. albicans genes required for filamentation in other conditions, the genes important for this morphological transition upon internalization by macrophages remained largely enigmatic. Here, we employed a functional genomic approach to identify genes that enable C. albicans filamentation within macrophages and uncovered a role for the mitochondrial ribosome, respiration, and the SNF1 AMP-activated kinase complex. Additionally, we showed that glucose uptake and glycolysis by macrophages support C. albicans filamentation. This work provides insights into the metabolic dueling that occurs during the interaction of C. albicans with macrophages and identifies vulnerabilities in C. albicans that could serve as promising therapeutic targets.

Mitochondrial genome of Cordyceps blackwelliae: organization, transcription, and evolutionary insights into Cordyceps

Cordyceps is a diverse genus of insect pathogenic fungi, with about 180 accepted species, including some well-known ones used as ethnic medicine and/or functional food. Nevertheless, mitogenomes are only available for four members of the genus. The current study reports the mitogenome of Cordyceps blackwelliae, a newly described entomopathogenic fungus. The 42,257-bp mitogenome of the fungus encoded genes typically found in fungal mitogenomes, and a total of 14 introns inserted into seven genes, including cob (1 intron), cox1 (4), cox3 (3), nad1 (1), nad4 (1), nad5 (1), and rnl (3). RNA-Seq analysis revealed differential expression of mitochondrial genes and supported annotations resulting from in silico analysis. There was clear evidence for polycistronic transcription and alternative splicing of mitochondrial genes. Comparison among mitogenomes of five different Cordyceps species (i.e., C. blackwelliae, C. chanhua, C. militaris, C. pruinosa, and C. tenuipes) revealed a high synteny, with mitogenome size expansion correlating with intron insertions. Different mitochondrial protein-coding genes showed variable degrees of genetic differentiation among these species, but they were all under purifying selection. Mitochondrial phylogeny based on either nucleotide or amino acid sequences confirmed the taxonomic position of C. blackwelliae in Cordycipitaceae, clustering together with C. chanhua. This study promotes our understanding of fungal evolution in Cordyceps.

Horizontal transfer of tRNA genes to mitochondrial plasmids facilitates gene loss from fungal mitochondrial DNA

Fungal and plant mitochondria are known to exchange DNA with retroviral plasmids. Transfer of plasmid DNA to the organellar genome is best known and occurs through wholesale insertion of the plasmid. Less well known is the transfer of organellar DNA to plasmids, in particular tRNA genes. Presently, it is unknown whether fungal plasmids can adopt mitochondrial functions such as tRNA production through horizontal gene transfer. In this paper, we studied the exchange of DNA between fungal linear plasmids and fungal mtDNA, mainly focusing on the basidiomycete family Lyophyllaceae. We report at least six independent transfers of complete tRNA genes to fungal plasmids. Furthermore, we discovered two independent cases of loss of a tRNA gene from a fungal mitochondrial genome following transfer of such a gene to a linear mitochondrial plasmid. We propose that loss of a tRNA gene from mtDNA following its transfer to a plasmid creates a mutualistic dependency of the host mtDNA on the plasmid. We also find that tRNA genes transferred to plasmids encode codons that occur at the lowest frequency in the host mitochondrial genomes, possibly due to a higher number of unused transcripts. We discuss the potential consequences of mtDNA transfer to plasmids for both the host mtDNA and the plasmid.

The Pet127 protein is a mitochondrial 5'-to-3' exoribonuclease from the PD-(D/E)XK superfamily involved in RNA maturation and intron degradation in yeasts

Pet127 is a mitochondrial protein found in multiple eukaryotic lineages, but absent from several taxa, including plants and animals. Distant homology suggests that it belongs to the divergent PD-(D/E)XK superfamily which includes various nucleases and related proteins. Earlier yeast genetics experiments suggest that it plays a nonessential role in RNA degradation and 5' end processing. Our phylogenetic analysis suggests that it is a primordial eukaryotic invention that was retained in diverse groups, and independently lost several times in the evolution of other organisms. We demonstrate for the first time that the fungal Pet127 protein in vitro is a processive 5'-to-3' exoribonuclease capable of digesting various substrates in a sequence nonspecific manner. Mutations in conserved residues essential in the PD-(D/E)XK superfamily active site abolish the activity of Pet127. Deletion of the PET127 gene in the pathogenic yeast Candida albicans results in a moderate increase in the steady-state levels of several transcripts and in accumulation of unspliced precursors and intronic sequences of three introns. Mutations in the active site residues result in a phenotype identical to that of the deletant, confirming that the exoribonuclease activity is related to the physiological role of the Pet127 protein. Pet127 activity is, however, not essential for maintaining the mitochondrial respiratory activity in C. albicans.

Mitochondrial Transcription of Entomopathogenic Fungi Reveals Evolutionary Aspects of Mitogenomes

Entomopathogenic fungi and more specifically genera Beauveria and Metarhizium have been exploited for the biological control of pests. Genome analyses are important to understand better their mode of action and thus, improve their efficacy against their hosts. Until now, the sequences of their mitochondrial genomes were studied, but not at the level of transcription. Except of yeasts and Neurospora crassa, whose mt gene transcription is well described, in all other Ascomycota, i.e., Pezizomycotina, related information is extremely scarce. In this work, mt transcription and key enzymes of this function were studied. RT-PCR experiments and Northern hybridizations reveal the transcriptional map of the mt genomes of B. bassiana and M. brunneum species. The mt genes are transcribed in six main transcripts and undergo post-transcriptional modifications to create single gene transcripts. Promoters were determined in both mt genomes with a comparative in silico analysis, including all known information from other fungal mt genomes. The promoter consensus sequence is 5'-ATAGTTATTAT-3' which is in accordance with the definition of the polycistronic transcripts determined with the experiments described above. Moreover, 5'-RACE experiments in the case of premature polycistronic transcript nad1-nad4-atp8-atp6 revealed the 5' end of the RNA transcript immediately after the in silico determined promoter, as also found in other fungal species. Since several conserved elements were retrieved from these analyses compared to the already known data from yeasts and N. crassa, the phylogenetic analyses of mt RNA polymerase (Rpo41) and its transcriptional factor (Mtf1) were performed in order to define their evolution. As expected, it was found that fungal Rpo41 originate from the respective polymerase of T7/T3 phages, while the ancestor of Mtf1 is of alpha-proteobacterial origin. Therefore, this study presents insights about the fidelity of the mt single-subunit phage-like RNA polymerase during transcription, since the correct identification of mt promoters from Rpo41 requires an ortholog to bacterial sigma factor, i.e., Mtf1. Thus, a previously proposed hypothesis of a phage infected alpha-proteobacterium as the endosymbiotic progenitor of mitochondrion is confirmed in this study and further upgraded by the co-evolution of the bacterial (Mtf1) and viral (Rpo41) originated components in one functional unit.

Global view of dynamic expression and precise mapping of mitochondrial tRNAs-derived fragments during stressed conditions in S. pombe

In this study, we provide a global view of population and processing of mitochondrial tRNAs-derived fragments (mt-tRFs) in fission yeast Schizosaccharomyces pombe. Here, mt-tRFs of 15-30 nucleotides were retrieved from S. pombe small RNA libraries obtained from unstressed, stress, and during stationary phase conditions. We demonstrate that production of these fragments increase during heat stress and stationary phase conditions in S. pombe, especially (most notably) in stationary phase. Analysis of data also reveals depending on the tRNA, either 5'-mt-tRF or 3'-mt-tRF was found and major mt-tRNA processing sites have been precisely identified. Furthermore, RNA-seq reveals that inactivation of trz2 encoding S. pombe mitochondrial tRNase Z^L globally impairs mt-tRF processing. Finally, our result showed mt-tRFs were predicted to target mitochondrial genome mapping mtDNA-encoded protein gene. These observations suggest that mitochondrial tRFs may play an important regulatory role in response to stress and development.

Pervasive transcription of the mitochondrial genome in Candida albicans is revealed in mutants lacking the mtEXO RNase complex

The mitochondrial genome of the pathogenic yeast Candida albicans displays a typical organization of several (eight) primary transcription units separated by noncoding regions. Presence of genes encoding Complex I subunits and the stability of its mtDNA sequence make it an attractive model to study organellar genome expression using transcriptomic approaches. The main activity responsible for RNA degradation in mitochondria is a two-component complex (mtEXO) consisting of a 3ʹ-5ʹ exoribonuclease, in yeasts encoded by the DSS1 gene, and a conserved Suv3p helicase. In C. albicans, deletion of either DSS1 or SUV3 gene results in multiple defects in mitochondrial genome expression leading to the loss of respiratory competence. Transcriptomic analysis reveals pervasive transcription in mutants lacking the mtEXO activity, with evidence of the entire genome being transcribed, whereas in wild-type strains no RNAs corresponding to a significant fraction of the noncoding genome can be detected. Antisense (‘mirror’) transcripts, absent from normal mitochondria are also prominent in the mutants. The expression of multiple mature transcripts, particularly those translated from bicistronic mRNAs, as well as those that contain introns is affected in the mutants, resulting in a decreased level of proteins and reduced respiratory complex activity. The phenotype is most severe in the case of Complex IV, where a decrease of mature COX1 mRNA level to ~5% results in a complete loss of activity. These results show that RNA degradation by mtEXO is essential for shaping the mitochondrial transcriptome and is required to maintain the functional demarcation between transcription units and non-coding genome segments.

The mitochondrial genome of the grape powdery mildew pathogen Erysiphe necator is intron rich and exhibits a distinct gene organization

Powdery mildews are notorious fungal plant pathogens but only limited information exists on their genomes. Here we present the mitochondrial genome of the grape powdery mildew fungus Erysiphe necator and a high-quality mitochondrial gene annotation generated through cloning and Sanger sequencing of full-length cDNA clones. The E. necator mitochondrial genome consists of a circular DNA sequence of 188,577 bp that harbors a core set of 14 protein-coding genes that are typically present in fungal mitochondrial genomes, along with genes encoding the small and large ribosomal subunits, a ribosomal protein S3, and 25 mitochondrial-encoded transfer RNAs (mt-tRNAs). Interestingly, it also exhibits a distinct gene organization with atypical bicistronic-like expression of the nad4L/nad5 and atp6/nad3 gene pairs, and contains a large number of 70 introns, making it one of the richest in introns mitochondrial genomes among fungi. Sixty-four intronic ORFs were also found, most of which encoded homing endonucleases of the LAGLIDADG or GIY-YIG families. Further comparative analysis of five E. necator isolates revealed 203 polymorphic sites, but only five were located within exons of the core mitochondrial genes. These results provide insights into the organization of mitochondrial genomes of powdery mildews and represent valuable resources for population genetic and evolutionary studies.

Clade II Candida auris possess genomic structural variations related to an ancestral strain

Candida auris is an invasive and multidrug-resistant ascomycetous yeast that is under global surveillance. All clinical cases of C. auris infection diagnosed from 1997 to 2019 in Japan were non-invasive and sporadic otitis media cases. In the present study, we performed whole-genome sequencing of seven C. auris strains isolated from patients with otitis media in Japan, all of which belonged to clade II. Comparative genome analysis using the high-quality draft genome sequences JCM 15448T revealed that single nucleotide variations (SNVs), clade-specific accessory genes, and copy number variations (CNVs) were identified in each C. auris clade. A total of 61 genes involved in cell wall and stress response-related functions was absent in clade II, and the pattern of conserved CNVs in each clade was more stable in clade II than in other clades. Our data suggest that the genomic structural diversity is stable in C. auris isolated from each biogeographic location, and Japanese strains isolated from patients with otitis media might belong to an ancestral type of C. auris. One Japanese strain, TWCC 58362, with reduced susceptibility to fluconazole, exhibited no mutation in ergosterol biosynthesis-related genes (ERG). However, TWCC 58362-specific variations, including SNVs, indels, and CNVs were detected, suggesting that gene duplication events in C. auris might contribute to antifungal drug resistance. Taken together, we demonstrated that genomic structural variations in C. auris could correlate to geographical dissemination, epidemiology, lesions in the host, and antifungal resistance.

Differential actinodin1 regulation in embryonic development and adult fin regeneration in Danio rerio

Actinotrichia are the first exoskeletal elements formed during zebrafish fin development. These rigid fibrils serve as skeletal support for the fin fold and as substrates for mesenchymal cell migration. In the adult intact fins, actinotrichia are restricted to the distal domain of the fin. Following fin amputation, actinotrichia also reform during regeneration. The actinodin gene family codes for structural proteins of actinotrichia. We have previously identified cis-acting regulatory elements in a 2kb genomic region upstream of the first exon of actinodin1, termed 2P, required for tissue-specific expression in the fin fold ectoderm and mesenchyme during embryonic development. Indeed, 2P contains an ectodermal enhancer in a 150bp region named epi. Deletion of epi from 2P results in loss of ectodermal-specific activity. In the present study, we sought to further characterize the activity of these regulatory sequences throughout fin development and during adult fin regeneration. Using a reporter transgenic approach, we show that a site within the epi region, termed epi3, contains an early mesenchymal-specific repressor. We also show that the larval fin fold ectodermal enhancer within epi3 remains functional in the basal epithelial layer during fin regeneration. We show that the first non-coding exon and first intron of actinodin1 contains a transcriptional enhancer and an alternative promoter that are necessary for the persistence of reporter expression reminiscent of actinodin1 expression during adulthood. Altogether, we have identified cis-acting regulatory elements that are required for tissue-specific expression as well as full recapitulation of actinodin1 expression during adulthood. Furthermore, the characterization of these elements provides us with useful molecular tools for the enhancement of transgene expression in adulthood.

The S. pombe mitochondrial transcriptome

Mitochondrial gene expression is largely controlled through post-transcriptional processes including mitochondrial RNA (mt-RNA) processing, modification, decay, and quality control. Defective mitochondrial gene expression results in mitochondrial oxidative phosphorylation (OXPHOS) deficiency and has been implicated in human disease. To fully understand mitochondrial transcription and RNA processing, we performed RNA-seq analyses of mt-RNAs from the fission yeast Schizosaccharomyces pombe RNA-seq analyses show that the abundance of mt-RNAs vary greatly. Analysis of data also reveals mt-RNA processing sites including an unusual RNA cleavage event by mitochondrial tRNA (mt-tRNA) 5'-end processing enzyme RNase P. Additionally, this analysis reveals previously unknown mitochondrial transcripts including the rnpB-derived fragment, mitochondrial small RNAs (mitosRNAs) such as mt-tRNA-derived fragments (mt-tRFs) and mt-tRNA halves, and mt-tRNAs marked with 3'-CCACCA/CCACC in S. pombe Finally, RNA-seq reveals that inactivation of trz2 encoding S. pombe mitochondrial tRNA 3'-end processing enzyme globally impairs mt-tRNA 3'-end processing, inhibits mt-mRNA 5'-end processing, and causes accumulation of unprocessed transcripts, demonstrating the feasibility of using RNA-seq to examine the protein known or predicted to be involved in mt-RNA processing in S. pombe Our work uncovers the complexity of a fungal mitochondrial transcriptome and provides a framework for future studies of mitochondrial gene expression using S. pombe as a model system.

Evidence for Mitochondrial Genome Methylation in the Yeast Candida albicans: A Potential Novel Epigenetic Mechanism Affecting Adaptation and Pathogenicity?

The commensal yeast Candida albicans is an opportunistic pathogen. In order to successfully colonize or infect the human body, the fungus must adapt to the host’s environmental conditions, such as low oxygen tension (hypoxia), temperature (37°C), and the different carbon sources available. Previous studies demonstrated the adaptive importance of C. albicans genetic variability for its pathogenicity, although the contributions of epigenetic and the influence of environmental factors are not fully understood. Mitochondria play important roles in fungal energetic metabolism, regulation of nuclear epigenetic mechanisms and pathogenicity. However, the specific impact of inter-strain mitochondrial genome variability and mitochondrial epigenetics in pathogenicity is unclear. Here, we draw attention to this relevant organelle and its potential role in C. albicans pathogenicity and provide preliminary evidence, for the first time, for methylation of the yeast mitochondrial genome. Our results indicate that environmental conditions, such as continuous exposure for 12 weeks to hypoxia and 37°C, decrease the mitochondrial genome methylation in strains SC5314 and L757. However, the methylation decrease is quantitatively different in specific genome positions when strains SC5314 and L757 are compared. We hypothesize that this phenomenon can be promising for future research to understand how physical factors of the host affect the C. albicans mitochondrial genome and its possible impact on adaptation and pathogenicity.

Getting into mitochondria

The human mitochondrial glutamate dehydrogenase isoenzymes (hGDH1 and hGDH2) are abundant matrix-localized proteins encoded by nuclear genes. The proteins are synthesized in the cytoplasm, with an atypically long N-terminal mitochondrial targeting sequence (MTS). The results of secondary structure predictions suggest the presence of two α-helices within the N-terminal region of the MTS. Results from deletion analyses indicate that individual helices have limited ability to direct protein import and matrix localization, but that there is a synergistic interaction when both helices are present [Biochem. J. (2016) 473: , 2813-2829]. Mutagenesis of the MTS cleavage sites blocked post-import removal of the presequences, but did not impede import. The authors propose that the high matrix levels of hGDH can be attributed to the unusual length and secondary structure of the MTS.

Traditional and novel tools to probe the mitochondrial metabolism in health and disease

Mitochondrial metabolism links energy production to other essential cellular processes such as signaling, cellular differentiation, and apoptosis. In addition to producing adenosine triphosphate (ATP) as an energy source, mitochondria are responsible for the synthesis of a myriad of important metabolites and cofactors such as tetrahydrofolate, α-ketoacids, steroids, aminolevulinic acid, biotin, lipoic acid, acetyl-CoA, iron-sulfur clusters, heme, and ubiquinone. Furthermore, mitochondria and their metabolism have been implicated in aging and several human diseases, including inherited mitochondrial disorders, cardiac dysfunction, heart failure, neurodegenerative diseases, diabetes, and cancer. Therefore, there is great interest in understanding mitochondrial metabolism and the complex relationship it has with other cellular processes. A large number of studies on mitochondrial metabolism have been conducted in the last 50 years, taking a broad range of approaches. In this review, we summarize and discuss the most commonly used tools that have been used to study different aspects of the metabolism of mitochondria: ranging from dyes that monitor changes in the mitochondrial membrane potential and pharmacological tools to study respiration or ATP synthesis, to more modern tools such as genetically encoded biosensors and trans-omic approaches enabled by recent advances in mass spectrometry, computation, and other technologies. These tools have allowed the large number of studies that have shaped our current understanding of mitochondrial metabolism. WIREs Syst Biol Med 2017, 9:e1373. doi: 10.1002/wsbm.1373 For further resources related to this article, please visit the WIREs website.

Pentatricopeptide Motifs in the N-Terminal Extension Domain of Yeast Mitochondrial RNA Polymerase Rpo41p Are Not Essential for Its Function

The core mitochondrial RNA polymerase is a single-subunit enzyme that in yeast Saccharomyces cerevisiae is encoded by the nuclear RPO41 gene. It is an evolutionary descendant of the bacteriophage RNA polymerases, but it includes an additional unconserved N-terminal extension (NTE) domain that is unique to the organellar enzymes. This domain mediates interactions between the polymerase and accessory regulatory factors, such as yeast Sls1p and Nam1p. Previous studies demonstrated that deletion of the entire NTE domain results only in a temperature-dependent respiratory deficiency. Several sequences related to the pentatricopeptide (PPR) motifs were identified in silico in Rpo41p, three of which are located in the NTE domain. PPR repeat proteins are a large family of organellar RNA-binding factors, mostly involved in posttranscriptional gene expression mechanisms. To study their function, we analyzed the phenotype of strains bearing Rpo41p variants where each of these motifs was deleted. We found that deletion of any of the three PPR motifs in the NTE domain does not affect respiratory growth at normal temperature, and it results in a moderate decrease in mtDNA stability. Steady-state levels of COX1 and COX2 mRNAs are also moderately affected. Only the deletion of the second motif results in a partial respiratory deficiency, manifested only at elevated temperature. Our results thus indicate that the PPR motifs do not play an essential role in the function of the NTE domain of the mitochondrial RNA polymerase.

The (in)complete organelle genome: exploring the use and nonuse of available technologies for characterizing mitochondrial and plastid chromosomes

Not long ago, scientists paid dearly in time, money and skill for every nucleotide that they sequenced. Today, DNA sequencing technologies epitomize the slogan 'faster, easier, cheaper and more', and in many ways, sequencing an entire genome has become routine, even for the smallest laboratory groups. This is especially true for mitochondrial and plastid genomes. Given their relatively small sizes and high copy numbers per cell, organelle DNAs are currently among the most highly sequenced kind of chromosome. But accurately characterizing an organelle genome and the information it encodes can require much more than DNA sequencing and bioinformatics analyses. Organelle genomes can be surprisingly complex and can exhibit convoluted and unconventional modes of gene expression. Unravelling this complexity can demand a wide assortment of experiments, from pulsed-field gel electrophoresis to Southern and Northern blots to RNA analyses. Here, we show that it is exactly these types of 'complementary' analyses that are often lacking from contemporary organelle genome papers, particularly short 'genome announcement' articles. Consequently, crucial and interesting features of organelle chromosomes are going undescribed, which could ultimately lead to a poor understanding and even a misrepresentation of these genomes and the genes they express. High-throughput sequencing and bioinformatics have made it easy to sequence and assemble entire chromosomes, but they should not be used as a substitute for or at the expense of other types of genomic characterization methods.