A novel motif for identifying rps3 homologs in fungal mitochondrial genomes

Intraspecific and interspecific variations in the synonymous codon usage in mitochondrial genomes of 8 pleurotus strains

In this study, we investigated the codon bias of twelve mitochondrial core protein coding genes (PCGs) in eight Pleurotus strains, two of which are from the same species. The results revealed that the codons of all Pleurotus strains had a preference for ending in A/T. Furthermore, the correlation between codon base compositions and codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) indices was also detected, implying the influence of base composition on codon bias. The two P. ostreatus species were found to have differences in various base bias indicators. The average effective number of codons (ENC) of mitochondrial core PCGs of Pleurotus was found to be less than 35, indicating strong codon preference of mitochondrial core PCGs of Pleurotus. The neutrality plot analysis and PR2-Bias plot analysis further suggested that natural selection plays an important role in Pleurotus codon bias. Additionally, six to ten optimal codons (ΔRSCU > 0.08 and RSCU > 1) were identified in eight Pleurotus strains, with UGU and ACU being the most widely used optimal codons in Pleurotus. Finally, based on the combined mitochondrial sequence and RSCU value, the genetic relationship between different Pleurotus strains was deduced, showing large variations between them. This research has improved our understanding of synonymous codon usage characteristics and evolution of this important fungal group.

The second complete mitochondrial genome of Capillidium rhysosporum within the family Capillidiaceae, Entomophthorales

The complete mitochondrial genome of the entomophthoroid fungus Capillidium rhysosporum (strain no.: ATCC 12588) was sequenced using next-generation sequencing technology. The assembled circular genome has a length of 46,756 base pairs with a GC content of 27.06%. Gene prediction identified 15 core protein-coding genes (PCGs), two rRNA genes, and 27 tRNA genes. Phylogenetic analysis confirmed that C. rhysosporum belongs to the Zoopagomycota clade and is closely related to C. heterosporum. This study presents the second complete mitochondrial genome within the family Capillidiaceae, contributing to the mitochondrial DNA database of entomophthoroid fungi.

Mitochondrial genome diversity across the subphylum Saccharomycotina

Introduction

Eukaryotic life depends on the functional elements encoded by both the nuclear genome and organellar genomes, such as those contained within the mitochondria. The content, size, and structure of the mitochondrial genome varies across organisms with potentially large implications for phenotypic variance and resulting evolutionary trajectories. Among yeasts in the subphylum Saccharomycotina, extensive differences have been observed in various species relative to the model yeast Saccharomyces cerevisiae, but mitochondrial genome sampling across many groups has been scarce, even as hundreds of nuclear genomes have become available.

Methods

By extracting mitochondrial assemblies from existing short-read genome sequence datasets, we have greatly expanded both the number of available genomes and the coverage across sparsely sampled clades.

Results

Comparison of 353 yeast mitochondrial genomes revealed that, while size and GC content were fairly consistent across species, those in the genera Metschnikowia and Saccharomyces trended larger, while several species in the order Saccharomycetales, which includes S. cerevisiae, exhibited lower GC content. Extreme examples for both size and GC content were scattered throughout the subphylum. All mitochondrial genomes shared a core set of protein-coding genes for Complexes III, IV, and V, but they varied in the presence or absence of mitochondrially-encoded canonical Complex I genes. We traced the loss of Complex I genes to a major event in the ancestor of the orders Saccharomycetales and Saccharomycodales, but we also observed several independent losses in the orders Phaffomycetales, Pichiales, and Dipodascales. In contrast to prior hypotheses based on smaller-scale datasets, comparison of evolutionary rates in protein-coding genes showed no bias towards elevated rates among aerobically fermenting (Crabtree/Warburg-positive) yeasts. Mitochondrial introns were widely distributed, but they were highly enriched in some groups. The majority of mitochondrial introns were poorly conserved within groups, but several were shared within groups, between groups, and even across taxonomic orders, which is consistent with horizontal gene transfer, likely involving homing endonucleases acting as selfish elements.

Discussion

As the number of available fungal nuclear genomes continues to expand, the methods described here to retrieve mitochondrial genome sequences from these datasets will prove invaluable to ensuring that studies of fungal mitochondrial genomes keep pace with their nuclear counterparts.

The mitogenomes of Leptographium aureum, Leptographium sp., and Grosmannia fruticeta: expansion by introns

Introduction

Many members of the Ophiostomatales are of economic importance as they are bark-beetle associates and causative agents for blue stain on timber and in some instances contribute towards tree mortality. The taxonomy of these fungi has been challenging due to the convergent evolution of many traits associated with insect dispersal and a limited number of morphological characters that happen to be highly pleomorphic. This study examines the mitochondrial genomes for three members of Leptographium sensu lato [Leptographium aureum (also known as Grosmannia aurea), Grosmannia fruticeta (also known as Leptographium fruticetum), and Leptographium sp. WIN(M)1376)].

Methods

Illumina sequencing combined with gene and intron annotations and phylogenetic analysis were performed.

Results

Sequence analysis showed that gene content and gene synteny are conserved but mitochondrial genome sizes were variable: G. fruticeta at 63,821 bp, Leptographium sp. WIN(M)1376 at 81,823 bp and L. aureum at 104,547 bp. The variation in size is due to the number of introns and intron-associated open reading frames. Phylogenetic analysis of currently available mitochondrial genomes for members of the Ophiostomatales supports currently accepted generic arrangements within this order and specifically supports the separation of members with Leptographium-like conidiophores into two genera, with L. aureum grouping with Leptographium and G. fruticeta aligning with Grosmannia.

Discussion

Mitochondrial genomes are promising sequences for resolving evolutionary relationships within the Ophiostomatales.

Mitochondrial genome annotation with MFannot: a critical analysis of gene identification and gene model prediction

Compared to nuclear genomes, mitochondrial genomes (mitogenomes) are small and usually code for only a few dozen genes. Still, identifying genes and their structure can be challenging and time-consuming. Even automated tools for mitochondrial genome annotation often require manual analysis and curation by skilled experts. The most difficult steps are (i) the structural modelling of intron-containing genes; (ii) the identification and delineation of Group I and II introns; and (iii) the identification of moderately conserved, non-coding RNA (ncRNA) genes specifying 5S rRNAs, tmRNAs and RNase P RNAs. Additional challenges arise through genetic code evolution which can redefine the translational identity of both start and stop codons, thus obscuring protein-coding genes. Further, RNA editing can render gene identification difficult, if not impossible, without additional RNA sequence data. Current automated mito- and plastid-genome annotators are limited as they are typically tailored to specific eukaryotic groups. The MFannot annotator we developed is unique in its applicability to a broad taxonomic scope, its accuracy in gene model inference, and its capabilities in intron identification and classification. The pipeline leverages curated profile Hidden Markov Models (HMMs), covariance (CMs) and ERPIN models to better capture evolutionarily conserved signatures in the primary sequence (HMMs and CMs) as well as secondary structure (CMs and ERPIN). Here we formally describe MFannot, which has been available as a web-accessible service (https://megasun.bch.umontreal.ca/apps/mfannot/) to the research community for nearly 16 years. Further, we report its performance on particularly intron-rich mitogenomes and describe ongoing and future developments.

Functional annotation and comparative analysis of four Botrytis cinerea mitogenomes reported from Punjab, Pakistan

Botrytis cinerea is one of the top phytopathogenic fungus which ubiquitously cause grey mold on a variety of horticultural plants. The mechanism of respiration in the fungus occurs within the mitochondria. Mitogenomes serve as a key molecular marker for the investigation of fungal evolutionary patterns. This study aimed at the complete assembly, characterization, and comparative relationship of four mitogenomes of Botrytis cinerea strains including Kst5C, Kst14A, Kst32B, Kst33A, respectively. High throughput sequencing of four mitogenomes allowed the full assembly and annotation of these sequences. The total genome length of these 4 isolates Kst5C Kst14A, Kst32B, Kst33A was 69,986 bp, 77,303 bp, 76,204 bp and 55, 226 bp respectively. The distribution of features represented 2 ribosomal RNA genes,14 respiration encoding proteins, 1 mitochondrial ribosomal protein-encoding gene, along with varying numbers of transfer RNA genes, protein-coding genes, mobile intronic regions and homing endonuclease genes including LAGLIDADG and GIY-YIG domains were found in all four mitogenomes. The comparative analyses performed also decipher significant results for four mitogenomes among fungal isolates included in the study. This is the first report on the detailed annotation of mitogenomes as a proof for investigation of variation patterns present with in the B. cinerea causing grey mold on strawberries in Pakistan. This study will also contribute to the rapid evolutionary analysis and population patterns present among Botrytis cinerea.

Analysis of synonymous codon usage patterns in mitochondrial genomes of nine Amanita species

Introduction

Codon basis is a common and complex natural phenomenon observed in many kinds of organisms.

Methods

In the present study, we analyzed the base bias of 12 mitochondrial core protein-coding genes (PCGs) shared by nine Amanita species.

Results

The results showed that the codons of all Amanita species tended to end in A/T, demonstrating the preference of mitochondrial codons of Amanita species for a preference for this codon. In addition, we detected the correlation between codon base composition and the codon adaptation index (CAI), codon bias index (CBI), and frequency of optimal codons (FOP) indices, indicating the influence of base composition on codon bias. The average effective number of codons (ENC) of mitochondrial core PCGs of Amanita is 30.81, which is <35, demonstrating the strong codon preference of mitochondrial core PCGs of Amanita. The neutrality plot analysis and PR2-Bias plot analysis further demonstrated that natural selection plays an important role in Amanita codon bias. In addition, we obtained 5–10 optimal codons (ΔRSCU > 0.08 and RSCU > 1) in nine Amanita species, and GCA and AUU were the most widely used optimal codons. Based on the combined mitochondrial sequence and RSCU value, we deduced the genetic relationship between different Amanita species and found large variations between them.

Discussion

This study promoted the understanding of synonymous codon usage characteristics and evolution of this important fungal group.

Comprehensive analysis of codon bias in 13 Ganoderma mitochondrial genomes

Introduction

Codon usage bias is a prevalent phenomenon observed across various species and genes. However, the specific attributes of codon usage in the mitochondrial genome of Ganoderma species remain unknown.

Methods

In this study, we investigated the codon bias of 12 mitochondrial core protein-coding genes (PCGs) in 9 Ganoderma species, including 13 Ganoderma strains.

Results

The codons of all Ganoderma strains showed a preference for ending in A/T. Additionally, correlations between codon base composition and the codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) were identified, demonstrating the impact of base composition on codon bias. Various base bias indicators were found to vary between or within Ganoderma strains, including GC3s, the CAI, the CBI, and the FOP. The results also revealed that the mitochondrial core PCGs of Ganoderma have an average effective number of codons (ENC) lower than 35, indicating strong bias toward certain codons. Evidence from neutrality plot and PR2-bias plot analysis indicates that natural selection is a major factor affecting codon bias in Ganoderma. Additionally, 11 to 22 optimal codons (ΔRSCU>0.08 and RSCU>1) were identified in 13 Ganoderma strains, with GCA, AUC, and UUC being the most widely used optimal codons in Ganoderma. By analyzing the combined mitochondrial sequences and relative synonymous codon usage (RSCU) values, the genetic relationships between or within Ganoderma strains were determined, indicating variations between them. Nevertheless, RSCU-based analysis illustrated the intra- and interspecies relationships of certain Ganoderma species.

Discussion

This study deepens our insight into the synonymous codon usage characteristics, genetics, and evolution of this important fungal group.

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.

The Gene Rearrangement, Loss, Transfer, and Deep Intronic Variation in Mitochondrial Genomes of Conidiobolus

The genus Conidiobolus s.s. was newly delimited from Conidiobolus s.l. In order to gain insight into its mitochondrial genetic background, this study sequenced six mitochondrial genomes of the genus Conidiobolus s.s. These mitogenomes were all composed of circular DNA molecules, ranging from 29,253 to 48,417 bp in size and from 26.61 to 27.90% in GC content. The order and direction for 14 core protein-coding genes (PCGs) were identical, except for the atp8 gene lost in Conidiobolus chlamydosporus, Conidiobolus polyspermus, and Conidiobolus polytocus, and rearranged in the other Conidiobolus s.s. species. Besides, the atp8 gene split the cox1 gene in Conidiobolus taihushanensis. Phylogenomic analysis based on the 14 core PCGs confirmed that all Conidiobolus s.s. species formed a monophyly in the Entomophthoromycotina lineage. The number and length of introns were the main factors contributing to mitogenomic size, and deep variations and potential transfer were detected in introns. In addition, gene transfer occurred between the mitochondrial and nuclear genomes. This study promoted the understanding of the evolution and phylogeny of the Conidiobolus s.s. genus.

Pan-Mitogenomics Approach Discovers Diversity and Dynamism in the Prominent Brown Rot Fungal Pathogens

Monilinia fructicola and Monilinia laxa species are the most destructive and economically devastating fungal plant pathogens causing brown rot disease on stone and pome fruits worldwide. Mitochondrial genomes (mitogenomes) play critical roles influencing the mechanisms and directions of the evolution of fungal pathogens. The pan-mitogenomics approach predicts core and accessory regions of the mitochondrial genomes and explains the gain or loss of variation within and between species. The present study is a fungal pan-mitogenome of M. fructicola (N = 8) and M. laxa (N = 8) species. The completely sequenced and annotated mitogenomes showed high variability in size within and between the species. The mitogenomes of M. laxa were larger, ranging from 178,351 to 179,780bp, than the mitogenomes of M. fructicola, ranging from 158,607 to 167,838bp. However, size variation within the species showed that M. fructicola isolates were more variable in the size range than M. laxa isolates. All the mitogenomes included conserved mitochondrial genes, as well as variable regions including different mobile introns encoding homing endonucleases or maturase, non-coding introns, and repetitive elements. The linear model analysis supported the hypothesis that the mitogenome size expansion is due to presence of variable (accessory) regions. Gene synteny was mostly conserved among all samples, with the exception for order of the rps3 in the mitogenome of one isolate. The mitogenomes presented AT richness; however, A/T and G/C skew varied among the mitochondrial genes. The purifying selection was detected in almost all the protein-coding genes (PCGs) between the species. However, cytochrome b was the only gene showing a positive selection signal among the total samples. Combined datasets of amino acid sequences of 14 core mitochondrial PCGs and rps3 obtained from this study together with published mitochondrial genome sequences from some other species from Heliotales were used to infer a maximum likelihood (ML) phylogenetic tree. ML tree indicated that both Monilinia species highly diverged from each other as well as some other fungal species from the same order. Mitogenomes harbor much information about the evolution of fungal plant pathogens, which could be useful to predict pathogenic life strategies.

Comparative Mitogenomic Analysis Reveals Dynamics of Intron Within and Between Tricholoma Species and Phylogeny of Basidiomycota

The genus of Tricholoma is a group of important ectomycorrhizal fungi. The overlapping of morphological characteristics often leads to the confusion of Tricholoma species classification. In this study, the mitogenomes of five Tricholoma species were sequenced based on the next-generation sequencing technology, including T. matsutake SCYJ1, T. bakamatsutake, T. terreum, T. flavovirens, and T. saponaceum. These five mitogenomes were all composed of circular DNA molecules, with sizes ranging from 49,480 to 103,090 bp. Intergenic sequences were considered to be the main factor contributing to size variations of Tricholoma mitogenomes. Comparative mitogenomic analysis showed that the introns of the Agaricales mitogenome experienced frequent loss/gain events. In addition, potential gene transfer was detected between the mitochondrial and nuclear genomes of the five species of Tricholoma. Evolutionary analysis showed that the rps3 gene of the Tricholoma species was under positive selection or relaxed selection in the evolutionary process. In addition, large-scale gene rearrangements were detected between some Tricholoma species. Phylogenetic analysis using the Bayesian inference and maximum likelihood methods based on a combined mitochondrial gene set yielded identical and well-supported tree topologies. This study promoted the understanding of the genetics, evolution, and phylogeny of the Tricholoma genus and related species.

First characterization of the complete mitochondrial genome of fungal plant-pathogen Monilinia laxa which represents the mobile intron rich structure

Monilinia laxa is an important fungal plant pathogen causing brown rot on many stone and pome fruits worldwide. Mitochondrial genome (mitogenome) plays a critical role in evolutionary biology of the organisms. This study aimed to characterize the complete mitogenome of M. laxa by using next-generation sequencing and approaches of de novo assembly and annotation. The total length of the mitogenome of M. laxa was 178,357 bp, and its structure was circular. GC content of the mitogenome was 30.1%. Annotation of the mitogenome presented 2 ribosomal RNA (rRNA) genes, 32 transfer RNA genes (tRNA), 1 gene encoding mitochondrial ribosomal protein S3, 14 protein-coding genes and 15 open reading frame encoding hypothetical proteins. Moreover, the group I mobile introns encoding homing endonucleases including LAGLIDADG and GIY-YIG families were found both within coding regions (genic) and intergenic regions of the mitogenome, indicating an enlarged size and a dynamic structure of the mitogenome. Furthermore, a comparative mitogenomic analysis was performed between M. laxa and the three closely related fungal phytopathogen species (Botryotinia fuckeliana, Sclerotinia sclerotiorum and, S. borealis). Due to the number and distribution of introns, the large extent of structural rearrangements and diverse mitogenome sizes were detected among the species investigated. Monilinia laxa presented the highest number of homing endonucleases among the fungal species considered in the analyses. This study is the first to report a detailed annotation of the mitogenome of an isolate of M. laxa, providing a solid basis for further investigations of mitogenome variations for the other Monilinia pathogens causing brown rot disease.

Comparative Mitogenome Analysis Reveals Mitochondrial Genome Differentiation in Ectomycorrhizal and Asymbiotic Amanita Species

In this present study, we assembled and analyzed the mitogenomes of two asymbiotic and six ectomycorrhizal Amanita species based on next-generation sequencing data. The size of the eight Amanita mitogenomes ranged from 37,341 to 137,428 bp, and we considered introns to be one of the main factors contributing to the size variation of Amanita. The introns of the cox1 gene experienced frequent gain/loss events in Amanita; and the intron position class cox1P386 was lost in the six ectomycorrhizal Amanita species. In addition, ectomycorrhizal Amanita species had more repetitive sequences and fewer intergenic sequences than asymbiotic Amanita species in their mitogenomes. Large-scale gene rearrangements were detected in the Amanita species we tested, including gene displacements and inversions. On the basis of the combined mitochondrial gene set, we reconstructed the phylogenetic relationships of 66 Basidiomycetes. The six ectomycorrhizal Amanita species were of single origin, and the two saprophytic Amanita species formed two distinct clades. This study is the first to elucidate the functions of the mitogenome in the evolution and ecological adaptation of Amanita species.

Evolution of Fusarium tricinctum and Fusarium avenaceum mitochondrial genomes is driven by mobility of introns and of a new type of palindromic microsatellite repeats

Background

Increased contamination of European and Asian wheat and barley crops with “emerging” mycotoxins such as enniatins or beauvericin, produced by Fusarium avenaceum and Fusarium tricinctum, suggest that these phylogenetically close species could be involved in future food-safety crises.

Results

The mitochondrial genomes of F. tricinctum strain INRA104 and F. avenaceum strain FaLH27 have been annotated. A comparative analysis was carried out then extended to a set of 25 wild strains. Results show that they constitute two distinct species, easily distinguished by their mitochondrial sequences. The mitochondrial genetic variability is mainly located within the intergenic regions. Marks of variations show they have evolved (i) by Single Nucleotide Polymorphisms (SNPs), (ii) by length variations mediated by insertion/deletion sequences (Indels), and (iii) by length mutations generated by DNA sliding events occurring in mononucleotide (A)_n or (T)_n microsatellite type sequences arranged in a peculiar palindromic organization. The optionality of these palindromes between both species argues for their mobility. The presence of Indels and SNPs in palindrome neighbouring regions suggests their involvement in these observed variations. Moreover, the intraspecific and interspecific variations in the presence/absence of group I introns suggest a high mobility, resulting from several events of gain and loss during short evolution periods. Phylogenetic analyses of intron orthologous sequences suggest that most introns could have originated from lateral transfers from phylogenetically close or distant species belonging to various Ascomycota genera and even to the Basidiomycota fungal division.

Conclusions

Mitochondrial genome evolution between F. tricinctum and F. avenaceum is mostly driven by two types of mobile genetic elements, implicated in genome polymorphism. The first one is represented by group I introns. Indeed, both genomes harbour optional (inter- or intra-specifically) group I introns, all carrying putatively functional hegs, arguing for a high mobility of these introns during short evolution periods. The gain events were shown to involve, for most of them, lateral transfers between phylogenetically distant species. This study has also revealed a new type of mobile genetic element constituted by a palindromic arrangement of (A) n and (T) n microsatellite sequences whose presence was related to occurrence of SNPs and Indels in the neighbouring regions.

Complete Mitochondrial Genome of the Fungal Biocontrol Agent Trichoderma atroviride: Genomic Features, Comparative Analysis and Insight Into the Mitochondrial Evolution in Trichoderma

The improvement of biopesticides for use in the agriculture industry requires an understanding of the biological- and ecological principles underlying their behavior in natural environments. The nuclear genomes of members of the genus Trichoderma, which are representative fungal biocontrol agents, have been actively studied in relation to the unique characteristics of these species as effective producers of CAZymes/secondary metabolites and biopesticides, but their mitochondrial genomes have received much less attention. In this study, the mitochondrial genome of Trichoderma atroviride (Hypocreales, Sordariomycetes), which targets wood-decaying fungal pathogens and has the ability to degrade chemical fungicides, was assembled de novo. A 32,758 bp circular DNA molecule was revealed with specific features, such as a few more protein CDS and trn genes, two homing endonucleases (LAGLIDADG-/GIY-YIG-type), and even a putative overlapping tRNA gene, on a closer phylogenetic relationship with T. gamsii among hypocrealean fungi. Particularly, introns were observed with several footprints likely to be evolutionarily associated with the intron dynamics of the Trichoderma mitochondrial genomes. This study is the first to report the complete de novo mitochondrial genome of T. atroviride, while comparative analyses of Trichoderma mitochondrial genomes were also conducted from the perspective of mitochondrial evolution for the first time.

The complete mitochondrial genomes of two model ectomycorrhizal fungi (Laccaria): features, intron dynamics and phylogenetic implications

Laccaria amethystine and L. bicolor have served as model species for studying the life history and genetics of ectomycorrhizal fungi. However, the characterizations and variations of their mitogenomes are still unknown. In the present study, the mitogenomes of the two Laccaria species were assembled, annotated, and compared. The two mitogenomes of L. amethystine and L. bicolor comprised circular DNA molecules, with the sizes of 65,156 bp and 95,304 bp, respectively. Genome collinearity analysis revealed large-scale gene rearrangements between the two Laccaria species. Comparative mitogenome analysis indicated the introns of cox1 genes in Agaricales experienced frequent lost/gain eveants, which promoted the organization and size variations in Agaricales mitogenomes. Evolutionary analysis indicated the core protein-coding genes in the two mitogenomes were subject to strong pressure of purifying selection. Phylogenetic analysis using the Bayesian inference (BI) and Maximum likelihood (ML) methods based on a combined mitochondrial gene set resulted in identical and well-supported tree topologies, wherein the two Laccaria species were most closely related to Coprinopsis cinerea. This study severed as the first study on the mitogenomes of Laccaria species, which promoted a comprehensive understanding of the genetics and evolution of the model ectomycorrhizal fungi.

Mt-rps3 is an ancient gene which provides insight into the evolution of fungal mitochondrial genomes

The nuclear ribosomal protein S3 (Rps3) is implicated in the assembly of the ribosomal small subunit. Fungi and plants present a gene copy in their mitochondrial (mt) genomes. An analysis of 303 complete fungal mt genomes showed that, when rps3 is found, it is either a free-standing gene or an anchored gene within the omega intron of the rnl gene. Early divergent fungi, Basidiomycota and all yeasts but the CTG group belong to the first case, and Pezizomycotina to the second. Its position, size and genetic code employed are conserved within species of the same Order. Size variability is attributed to different number of repeats. These repeats consist of AT-rich sequences. MtRps3 proteins lack the KH domain, necessary for binding to rRNA, in their N-terminal region. Their C-terminal region is conserved in all Domains of life. Phylogenetic analysis showed that nuclear and mtRps3 proteins are descendants of archaeal and a-proteobacterial homologues, respectively. Thus, fungal mt-rps3 gene is an ancient gene which evolved within the endosymbiotic model and presents different evolutionary routes: (a) coming from a-proteobacteria, it was relocated to another region of the mt genome, (b) via its insertion to the omega intron, it was transferred to the nucleus and/or got lost, and (c) it was re-routed to the mt genome again. Today, Basidiomycota and Saccharomycetales seem to follow the first evolutionary route and almost all Pezizomycotina support the second scenario with their exceptions being the result of the third scenario, i.e., the gene's re-entry to the mt genome.

The structure of the yeast mitochondrial ribosome

Mitochondria have specialized ribosomes (mitoribosomes) dedicated to the expression of the genetic information encoded by their genomes. Here, using electron cryomicroscopy, we have determined the structure of the 75-component yeast mitoribosome to an overall resolution of 3.3 angstroms. The mitoribosomal small subunit has been built de novo and includes 15S ribosomal RNA (rRNA) and 34 proteins, including 14 without homologs in the evolutionarily related bacterial ribosome. Yeast-specific rRNA and protein elements, including the acquisition of a putatively active enzyme, give the mitoribosome a distinct architecture compared to the mammalian mitoribosome. At an expanded messenger RNA channel exit, there is a binding platform for translational activators that regulate translation in yeast but not mammalian mitochondria. The structure provides insights into the evolution and species-specific specialization of mitochondrial translation.

The mitochondrial genome of the plant-pathogenic fungus Stemphylium lycopersici uncovers a dynamic structure due to repetitive and mobile elements

Stemphylium lycopersici (Pleosporales) is a plant-pathogenic fungus that has been associated with a broad range of plant-hosts worldwide. It is one of the causative agents of gray leaf spot disease in tomato and pepper. The aim of this work was to characterize the mitochondrial genome of S. lycopersici CIDEFI-216, to use it to trace taxonomic relationships with other fungal taxa and to get insights into the evolutionary history of this phytopathogen. The complete mitochondrial genome was assembled into a circular double-stranded DNA molecule of 75,911 bp that harbors a set of 37 protein-coding genes, 2 rRNA genes (rns and rnl) and 28 tRNA genes, which are transcribed from both sense and antisense strands. Remarkably, its gene repertoire lacks both atp8 and atp9, contains a free-standing gene for the ribosomal protein S3 (rps3) and includes 13 genes with homing endonuclease domains that are mostly located within its 15 group I introns. Strikingly, subunits 1 and 2 of cytochrome oxidase are encoded by a single continuous open reading frame (ORF). A comparative mitogenomic analysis revealed the large extent of structural rearrangements among representatives of Pleosporales, showing the plasticity of their mitochondrial genomes. Finally, an exhaustive phylogenetic analysis of the subphylum Pezizomycotina based on mitochondrial data reconstructed their relationships in concordance with several studies based on nuclear data. This is the first report of a mitochondrial genome belonging to a representative of the family Pleosporaceae.

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.

High variability of mitochondrial gene order among fungi

From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

Massively convergent evolution for ribosomal protein gene content in plastid and mitochondrial genomes

Plastid and mitochondrial genomes have undergone parallel evolution to encode the same functional set of genes. These encode conserved protein components of the electron transport chain in their respective bioenergetic membranes and genes for the ribosomes that express them. This highly convergent aspect of organelle genome evolution is partly explained by the redox regulation hypothesis, which predicts a separate plastid or mitochondrial location for genes encoding bioenergetic membrane proteins of either photosynthesis or respiration. Here we show that convergence in organelle genome evolution is far stronger than previously recognized, because the same set of genes for ribosomal proteins is independently retained by both plastid and mitochondrial genomes. A hitherto unrecognized selective pressure retains genes for the same ribosomal proteins in both organelles. On the Escherichia coli ribosome assembly map, the retained proteins are implicated in 30S and 50S ribosomal subunit assembly and initial rRNA binding. We suggest that ribosomal assembly imposes functional constraints that govern the retention of ribosomal protein coding genes in organelles. These constraints are subordinate to redox regulation for electron transport chain components, which anchor the ribosome to the organelle genome in the first place. As organelle genomes undergo reduction, the rRNAs also become smaller. Below size thresholds of approximately 1,300 nucleotides (16S rRNA) and 2,100 nucleotides (26S rRNA), all ribosomal protein coding genes are lost from organelles, while electron transport chain components remain organelle encoded as long as the organelles use redox chemistry to generate a proton motive force.

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.

Complete mitochondrial genome of the aluminum-tolerant fungus Rhodotorula taiwanensis RS1 and comparative analysis of Basidiomycota mitochondrial genomes

The complete mitochondrial genome of Rhodotorula taiwanensis RS1, an aluminum-tolerant Basidiomycota fungus, was determined and compared with the known mitochondrial genomes of 12 Basidiomycota species. The mitochondrial genome of R. taiwanensis RS1 is a circular DNA molecule of 40,392 bp and encodes the typical 15 mitochondrial proteins, 23 tRNAs, and small and large rRNAs as well as 10 intronic open reading frames. These genes are apparently transcribed in two directions and do not show syntenies in gene order with other investigated Basidiomycota species. The average G+C content (41%) of the mitochondrial genome of R. taiwanensis RS1 is the highest among the Basidiomycota species. Two introns were detected in the sequence of the atp9 gene of R. taiwanensis RS1, but not in that of other Basidiomycota species. Rhodotorula taiwanensis is the first species of the genus Rhodotorula whose full mitochondrial genome has been sequenced; and the data presented here supply valuable information for understanding the evolution of fungal mitochondrial genomes and researching the mechanism of aluminum tolerance in microorganisms.

Mitochondrial genomes from the lichenized fungi Peltigera membranacea and Peltigera malacea: features and phylogeny

Mitochondrial genomes from the fungal partners of two terricolous foliose lichen symbioses, Peltigera membranacea and Peltigera malacea, have been determined using metagenomic approaches, including RNA-seq. The roughly 63 kb genomes show all the major features found in other Pezizomycotina, such as unidirectional transcription, 14 conserved protein genes, genes for the two subunit rRNAs and for a set of 26 tRNAs used in translating the 62 amino acid codons. In one of the tRNAs a CAU anticodon is proposed to be modified, via the action of the nuclear-encoded enzyme, tRNA Ile lysidine synthase, so that it recognizes the codon AUA (Ile) instead of AUG (Met). The overall arrangements and sequences of the two circular genomes are similar, the major difference being the inversion and deterioration of a gene encoding a type B DNA polymerase. Both genomes encode the RNA component of RNAse P, a feature seldom found in ascomycetes. The difference in genome size from the minimal ascomycete mitochondrial genomes is largely due to 17 and 20 group I introns, respectively, most associated with homing endonucleases and all found within protein-coding genes and the gene encoding the large subunit rRNA. One new intron insertion point was found, and an unusually small exon of seven nucleotides (nt) was identified and verified by RNA sequencing. Comparative analysis of mitochondrion-encoded proteins places the Peltigera spp., representatives of the class Lecanoromycetes, close to Leotiomycetes, Dothidiomycetes, and Sordariomycetes, in contrast to phylogenies found using nuclear genes.

Learning to live together: mutualism between self-splicing introns and their hosts

Group I and II introns can be considered as molecular parasites that interrupt protein-coding and structural RNA genes in all domains of life. They function as self-splicing ribozymes and thereby limit the phenotypic costs associated with disruption of a host gene while they act as mobile DNA elements to promote their spread within and between genomes. Once considered purely selfish DNA elements, they now seem, in the light of recent work on the molecular mechanisms regulating bacterial and phage group I and II intron dynamics, to show evidence of co-evolution with their hosts. These previously underappreciated relationships serve the co-evolving entities particularly well in times of environmental stress.

Analysis of the complete mitochondrial genome sequences of the soybean rust pathogens phakopsora pachyrhizi and p. meibomiae

The mitochondrial (mt) genomes of two soybean rust pathogens, Phakopsora pachyrhizi and P. meibomiae, have been sequenced. The mt genome of P. pachyrhizi is a circular 31 825-bp molecule with a mean GC content of 34.6%, while P. meibomiae possesses a 32 520-bp circular molecule with a mean GC content of 34.9%. Both mt genomes contain the genes encoding ATP synthase subunits 6, 8 and 9 (atp6, atp8 and atp9), cytochrome oxidase subunits I, II and III (cox1, cox2 and cox3), apocytochrome b (cob), reduced nicotinamide adenine dinucleotide ubiquinone oxidoreductase subunits (nad1, nad2, nad3, nad4, nad4L, nad5 and nad6), the large and small mt ribosomal RNA genes, one ORF coding for a ribosomal protein (rps3), and a set of 24 tRNA genes that recognize codons for all amino acids. The order of the protein-coding genes and tRNA is identical in the two Phakopsora species, and all genes are transcribed from the same DNA strand clockwise. Introns were identified in the cox1, cob and mnl genes of both species, with three of the introns having ORFs with motifs similar to the LAGLIDADG endonucleases of other fungi. Phylogenetic analysis of the 14 shared protein-coding genes agrees with commonly accepted fungal taxonomy.

RNA editing in six mitochondrial ribosomal protein genes of Didymium iridis

Similarity searches with Didymium iridis mitochondrial genomic DNA identified six possible ribosomal protein-coding regions, however, each region contained stop codons that would need to be removed by RNA editing to produce functional transcripts. RT-PCR was used to amplify these regions from total RNA for cloning and sequencing. Six functional transcripts were verified for the following ribosomal protein genes: rpS12, rpS7, rpL2, rpS19, rpS3, and rpL16. The editing events observed, such as single C and U nucleotide insertions and a dinucleotide insertion, were consistent with previously observed editing patterns seen in D. iridis. Additionally, a new form of insertional editing, a single A insertion, was observed in a conserved region of the rpL16 gene. While the majority of codons created by editing specify hydrophobic amino acids, a greater proportion of the codons created in these hydrophilic ribosomal proteins called for positively charged amino acids in comparison to the previously characterized hydrophobic respiratory protein genes. This first report of edited soluble mitochondrial ribosomal proteins in myxomycetes expands upon the RNA editing patterns previously seen; there was: a greater proportion of created codons specifying positively charged amino acids, a shift in the codon position edited, and the insertion of single A nucleotides.

The sporadic occurrence of a group I intron-like element in the mtDNA rnl gene of Ophiostoma novo-ulmi subsp. americana

The presence of group I intron-like elements within the U7 region of the mtDNA large ribosomal subunit RNA gene (rnl) was investigated in strains of Ophiostoma novo-ulmi subsp. americana from Canada, Europe and Eurasia, and in selected strains of O. ips, O. minus, O. piceae, O. ulmi, and O. himal-ulmi. This insertion is of interest as it has been linked previously to the generation of plasmid-like mtDNA elements in diseased strains of O. novo-ulmi. Among 197 O. novo-ulmi subsp. americana strains tested, 61 contained a 1.6kb insertion within the rnl-U7 region and DNA sequence analysis suggests the presence of a group I intron (IA1 type) that encodes a potential double motif LAGLIDADG homing endonuclease-like gene (HEG). Phylogenetic analysis of rnl-U7 intron encoded HEG-like elements supports the view that double motif HEGs originated from a duplication event of a single-motif HEG followed by a fusion event that combined the two copies into one open reading frame (ORF). The data also show that rnl-U7 intron encoded ORFs belong to a clade that includes ORFs inserted into different types of group I introns, e.g. IB, ID, IC3, IA1, present within a variety of different mtDNA genes, such as the small ribosomal subunit RNA gene (rns), apo-cytochrome b gene (cob), NADH dehydrogenase subunit 5 (nad5), cytochrome oxidase subunit 1 gene (coxI), and ATPase subunit 9 gene (atp9). We also compared the occurrence of the rnl-U7 intron in our collection of 227 strains with the presence of the rnl-U11 group I intron and concluded that the U7 intron appears to be an optional element and the U11 intron is probably essential among the strains tested.

Complete mitochondrial genome sequence of the wine yeast Candida zemplinina: intraspecies distribution of a novel group-IIB1 intron with eubacterial affiliations

The mtDNA of the ascomycetous wine yeast Candida zemplinina is a circularly mapping genome of 23,114 bp. It contains 35 genes coding for the seven basic subunits of oxidative phosporylation found in yeasts (the genes encoding for NADH oxidoreductase subunits are absent), the ribosomal protein Var1, two rRNAs and 25 tRNA genes. Although protein phylogenetic analysis showed a divergent mitochondrial genome, several traits appeared preserved. The conserved gene blocks between the mtDNAs of C. zemplinina and Candida glabrata were maintained and changes in gene order and putative promoters were due to restricted genome reshuffling. New heterogeneous hairpin elements were identified scattered throughout cox1 introns. The large subunit rRNA gene harboured the first group-IIB1 intron containing a putative active reverse transcriptase (RT) in mitochondrial genomes of fungi. Phylogenetic analysis of the RT protein confirmed its closer relationship to eubacterial intronic RTs, while being only distantly related to all other fungal mitochondrial group-II introns and RTs. The findings point towards an early migration event of a eubacterial group-II intron to the mitochondrial genome of C. zemplinina.

Reconstructing the evolution of the mitochondrial ribosomal proteome

For production of proteins that are encoded by the mitochondrial genome, mitochondria rely on their own mitochondrial translation system, with the mitoribosome as its central component. Using extensive homology searches, we have reconstructed the evolutionary history of the mitoribosomal proteome that is encoded by a diverse subset of eukaryotic genomes, revealing an ancestral ribosome of alpha-proteobacterial descent that more than doubled its protein content in most eukaryotic lineages. We observe large variations in the protein content of mitoribosomes between different eukaryotes, with mammalian mitoribosomes sharing only 74 and 43% of its proteins with yeast and Leishmania mitoribosomes, respectively. We detected many previously unidentified mitochondrial ribosomal proteins (MRPs) and found that several have increased in size compared to their bacterial ancestral counterparts by addition of functional domains. Several new MRPs have originated via duplication of existing MRPs as well as by recruitment from outside of the mitoribosomal proteome. Using sensitive profile-profile homology searches, we found hitherto undetected homology between bacterial and eukaryotic ribosomal proteins, as well as between fungal and mammalian ribosomal proteins, detecting two novel human MRPs. These newly detected MRPs constitute, along with evolutionary conserved MRPs, excellent new screening targets for human patients with unresolved mitochondrial oxidative phosphorylation disorders.

The complete mitochondrial genome of the vascular wilt fungus Verticillium dahliae: a novel gene order for Verticillium and a diagnostic tool for species identification

The complete sequence (27,184 bp) of the mitochondrial (mt) genome of the phytopathogenic fungus Verticillium dahliae has been determined. It contains 14 protein-coding genes related to oxidative phosphorylation, two rRNA genes and a set of 25 tRNA genes. A single intron, that harbors an intronic ORF coding for a putative ribosomal protein (rps), is located within the large rRNA gene (rnl). Gene order comparisons of V. dahliae mtDNA and complete mt genomes of Pezizomycotina revealed four units of synteny for Sordariomycetes, namely rnl-trn ((11-12))-nad2-nad3, nad4L-nad5-cob-cox1, nad1-nad4-atp8-atp6 and rns-trn ((1-5))-cox3-trn ((1-5))-nad6-trn ((2-5)). These four units, in different combinations, merged to single continuous unit in the orders of Hypocreales and Sordariales. V. dahliae (Phyllachorales) and all members of the genus showed a unique feature which is the translocation of the nad1-nad4-atp8-atp6-rns-cox3-nad6 region in between genes nad3 and atp9 of the Hypocreales mtDNA gene order. Analysis of mt intergenic sequences of Verticillium species permitted the design of a species-specific primer allowing the discrimination of V. longisporum against V. dahliae and V. albo-atrum. By considering the protein-coding gene sequences as one unit, a phylogenetic comparison with representatives of Ascomycota complete mtDNA was performed.

The complete mitochondrial genome of the entomopathogenic fungus Metarhizium anisopliae var. anisopliae: gene order and trn gene clusters reveal a common evolutionary course for all Sordariomycetes, while intergenic regions show variation

The mitochondrial genome (mtDNA) of the entomopathogenic fungus Metarhizium anisopliae var. anisopliae, with a total size of 24,673 bp, was one of the smallest known mtDNAs of Pezizomycotina. It contained the 14 typical genes coding for proteins related to oxidative phosphorylation, the two rRNA genes, a single intron that harbored an intronic ORF coding for a putative ribosomal protein (rps) within the large rRNA gene (rnl), and a set of 24 tRNA genes which recognized codons for all amino acids, except proline and valine. Gene order comparison with all known mtDNAs of Sordariomycetes illustrated a highly conserved genome organization for all the protein- and rRNA-coding genes, as well as three clusters of tRNA genes. By considering all mitochondrial essential protein-coding genes as one unit a phylogenetic study of these small genomes strongly supported the common evolutionary course of Sordariomycetes (100% bootstrap support) and highlighted the advantages of analyzing small genomes (mtDNA) over single genes. In addition, comparative analysis of three intergenic regions demonstrated sequence variability that can be exploited for intra- and inter-specific identification of Metarhizium.

Optional mitochondrial introns and evidence for a homing-endonuclease gene in the mtDNA rnl gene in Ophiostoma ulmi s. lat

Strains of Ophiostoma ulmi, O. novo-ulmi subsp. americana, O. novo-ulmi subsp. novo-ulmi and O. himal-ulmi were examined for optional introns/insertions within the following mitochondrial genes: small subunit RNA gene (rns), large ribosomal subunit gene (rnl) and the cytochrome oxidase subunit I gene (coxI). Insertions were noted in the rns and coxI genes in strains of O. ulmi, the less aggressive species, but absent in strains of the more aggressive O. novo-ulmi subsp. americana. Strains of all species examined had a group I intron present in the U11 region of the mitochondrial-rnl gene. In all but two strains of O. novo-ulmi subsp. americana, this rnl-U11 intron was about 1.5 kb in length whereas a 2.6 kb version of this element was present in all strains representing O. ulmi, O. novo-ulmi subsp. novo-ulmi, and Ophiostoma himal-ulmi. Irrespective of size, this intron based on RNA folds is a class IA1 group I intron and it encodes a putative ORF for the rps3 ribosomal protein. The size variation of the rnl-U11 intron was examined in detail for two strains of O. novo-ulmi subsp. americana and sequence data suggests the presence of a complex ORF within the 2.6 kb version of this intron; here a homing endonuclease-like gene has been inserted in frame and fused to the carboxyl-terminus of the putative rps3 coding region. The mitochondrial optional introns/insertions in combination with nuclear markers might be useful in distinguishing among the various species and subspecies of the O. ulmi s. lat. complex.

Fungal evolution: the case of the vanishing mitochondrion

Mitochondria, the energy-producing organelles of the eukaryotic cell, are derived from an ancient endosymbiotic alpha-Proteobacterium. These organelles contain their own genetic system, a remnant of the endosymbiont's genome, which encodes only a fraction of the mitochondrial proteome. The majority of mitochondrial proteins are translated from nuclear genes and are imported into mitochondria. Recent studies of phylogenetically diverse representatives of Fungi reveal that their mitochondrial DNAs are among the most highly derived, encoding only a limited set of genes. Much of the reduction in the coding content of the mitochondrial genome probably occurred early in fungal evolution. Nevertheless, genome reduction is an ongoing process. Fungi in the chytridiomycete order Neocallimastigales and in the pathogenic Microsporidia have taken mitochondrial reduction to the extreme and have permanently lost a mitochondrial genome. These organisms have organelles derived from mitochondria that retain traces of their mitochondrial ancestry.

The Rhodomonas salina mitochondrial genome: bacteria-like operons, compact gene arrangement and complex repeat region

To gain insight into the mitochondrial genome structure and gene content of a putatively ancestral group of eukaryotes, the cryptophytes, we sequenced the complete mitochondrial DNA of Rhodomonas salina. The 48 063 bp circular-mapping molecule codes for 2 rRNAs, 27 tRNAs and 40 proteins including 23 components of oxidative phosphorylation, 15 ribosomal proteins and two subunits of tat translocase. One potential protein (ORF161) is without assigned function. Only two introns occur in the genome; both are present within cox1 belong to group II and contain RT open reading frames. Primitive genome features include bacteria-like rRNAs and tRNAs, ribosomal protein genes organized in large clusters resembling bacterial operons and the presence of the otherwise rare genes such as rps1 and tatA. The highly compact gene organization contrasts with the presence of a 4.7 kb long, repeat-containing intergenic region. Repeat motifs ∼40–700 bp long occur up to 31 times, forming a complex repeat structure. Tandem repeats are the major arrangement but the region also includes a large, ∼3 kb, inverted repeat and several potentially stable ∼40–80 bp long hairpin structures. We provide evidence that the large repeat region is involved in replication and transcription initiation, predict a promoter motif that occurs in three locations and discuss two likely scenarios of how this highly structured repeat region might have evolved.

Comparative mitochondrial genomics in zygomycetes: bacteria-like RNase P RNAs, mobile elements and a close source of the group I intron invasion in angiosperms

To generate data for comparative analyses of zygomycete mitochondrial gene expression, we sequenced mtDNAs of three distantly related zygomycetes, Rhizopus oryzae, Mortierella verticillata and Smittium culisetae. They all contain the standard fungal mitochondrial gene set, plus rnpB, the gene encoding the RNA subunit of the mitochondrial RNase P (mtP-RNA) and rps3, encoding ribosomal protein S3 (the latter lacking in R.oryzae). The mtP-RNAs of R.oryzae and of additional zygomycete relatives have the most eubacteria-like RNA structures among fungi. Precise mapping of the 5′ and 3′ termini of the R.oryzae and M.verticillata mtP-RNAs confirms their expression and processing at the exact sites predicted by secondary structure modeling. The 3′ RNA processing of zygomycete mitochondrial mRNAs, SSU-rRNA and mtP-RNA occurs at the C-rich sequence motifs similar to those identified in fission yeast and basidiomycete mtDNAs. The C-rich motifs are included in the mature transcripts, and are likely generated by exonucleolytic trimming of RNA 3′ termini. Zygomycete mtDNAs feature a variety of insertion elements: (i) mtDNAs of R.oryzae and M.verticillata were subject to invasions by double hairpin elements; (ii) genes of all three species contain numerous mobile group I introns, including one that is closest to an intron that invaded angiosperm mtDNAs; and (iii) at least one additional case of a mobile element, characterized by a homing endonuclease insertion between partially duplicated genes [Paquin,B., Laforest,M.J., Forget,L., Roewer,I., Wang,Z., Longcore,J. and Lang,B.F. (1997) Curr. Genet., 31, 380–395]. The combined mtDNA-encoded proteins contain insufficient phylogenetic signal to demonstrate monophyly of zygomycetes.

Forces maintaining organellar genomes: is any as strong as genetic code disparity or hydrophobicity?

It remains controversial why mitochondria and chloroplasts retain the genes encoding a small subset of their constituent proteins, despite the transfer of so many other genes to the nucleus. Two candidate obstacles to gene transfer, suggested long ago, are that the genetic code of some mitochondrial genomes differs from the standard nuclear code, such that a transferred gene would encode an incorrect amino acid sequence, and that the proteins most frequently encoded in mitochondria are generally very hydrophobic, which may impede their import after synthesis in the cytosol. More recently it has been suggested that both these interpretations suffer from serious "false positives" and "false negatives": genes that they predict should be readily transferred but which have never (or seldom) been, and genes whose transfer has occurred often or early, even though this is predicted to be very difficult. Here I consider the full known range of ostensibly problematic such genes, with particular reference to the sequences of events that could have led to their present location. I show that this detailed analysis of these cases reveals that they are in fact wholly consistent with the hypothesis that code disparity and hydrophobicity are much more powerful barriers to functional gene transfer than any other. The popularity of the contrary view has led to the search for other barriers that might retain genes in organelles even more powerfully than code disparity or hydrophobicity; one proposal, concerning the role of proteins in redox processes, has received widespread support. I conclude that this abandonment of the original explanations for the retention of organellar genomes has been premature. Several other, relatively minor, obstacles to gene transfer certainly exist, contributing to the retention of relatively many organellar genes in most lineages compared to animal mtDNA, but there is no evidence for obstacles as severe as code disparity or hydrophobicity. One corollary of this conclusion is that there is currently no reason to suppose that engineering nuclear versions of the remaining mammalian mitochondrial genes, a feat that may have widespread biomedical relevance, should require anything other than sequence alterations obviating code disparity and causing modest reductions in hydrophobicity without loss of enzymatic function.

Mitochondria of Protists

Over the past several decades, our knowledge of the origin and evolution of mitochondria has been greatly advanced by determination of complete mitochondrial genome sequences. Among the most informative mitochondrial genomes have been those of protists (primarily unicellular eukaryotes), some of which harbor the most gene-rich and most eubacteria-like mitochondrial DNAs (mtDNAs) known. Comparison of mtDNA sequence data has provided insights into the radically diverse trends in mitochondrial genome evolution exhibited by different phylogenetically coherent groupings of eukaryotes, and has allowed us to pinpoint specific protist relatives of the multicellular eukaryotic lineages (animals, plants, and fungi). This comparative genomics approach has also revealed unique and fascinating aspects of mitochondrial gene expression, highlighting the mitochondrion as an evolutionary playground par excellence.