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

Characterization of the complete mitochondrial genome of Mucor indicus Lendn. 1930 (Mucorales: Mucoraceae), isolated from the wine fermentation system

Mucor indicus Lendn. 1930 has been widely used in food fermentation; however, its mitochondrial genome characteristics are not well understood. In this study, the complete mitochondrial genome of M. indicus was obtained, which was 61,400 bp in length with a GC content of 33%. The M. indicus mitochondrial genome was found to contain 14 core protein-coding genes, four free-standing ORFs, 18 intronic ORFs, 26 tRNAs, and two rRNA genes. Phylogenetic trees were generated for 25 early-differentiated fungi using the Bayesian inference (BI) method, which demonstrated that M. indicus is closely related to Mucor piriformis. This study provides useful information for the classification and evolution of Mucor species or other early-differentiated fungi.

Mitochondrial genomic characteristics and phylogenetic analysis of a brewing fungus, Rhizopus microsporus Tiegh. 1875 (Mucorales: Rhizopodaceae)

Rhizopus microsporus Tiegh. 1875 is widely used in a variety of industries, such as brewing, wine making, baking, and medicine production, as it has the capability to break down proteins and generate surface-active agents. To date, the mitochondrial genome features of early evolved fungi from the Rhizopus genus have not been extensively studied. Our research obtained a full mitochondrial genome of R. microsporus species, which was 43,837 bp in size and had a GC content of 24.93%. This genome contained 14 core protein-coding genes, 3 independent ORFs, 7 intronic ORFs, 24 tRNAs, and 2 rRNA genes. Through the use of the BI phylogenetic inference method, we were able to create phylogenetic trees for 25 early differentiation fungi which strongly supported the major clades; this indicated that R. microsporus is most closely related to Rhizopus oryzae.

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 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.

Comparative Genomics of Mortierellaceae Provides Insights into Lipid Metabolism: Two Novel Types of Fatty Acid Synthase

Fungal species in the family Mortierellaceae are important for their remarkable capability to synthesize large amounts of polyunsaturated fatty acids, especially arachidonic acid (ARA). Although many genomes have been published, the quality of these data is not satisfactory, resulting in an incomplete understanding of the lipid pathway in Mortierellaceae. We provide herein two novel and high-quality genomes with 55.32% of syntenic gene pairs for Mortierella alpina CGMCC 20262 and M. schmuckeri CGMCC 20261, spanning 28 scaffolds of 40.22 Mb and 25 scaffolds of 49.24 Mb, respectively. The relative smaller genome for the former is due to fewer protein-coding gene models (11,761 vs. 13,051). The former yields 45.57% of ARA in total fatty acids, while the latter 6.95%. The accumulation of ARA is speculated to be associated with delta-5 desaturase (Delta5) and elongation of very long chain fatty acids protein 3 (ELOVL3). A further genomic comparison of 19 strains in 10 species in three genera in the Mortierellaceae reveals three types of fatty acid synthase (FAS), two of which are new to science. The most common type I exists in 16 strains of eight species of three genera, and was discovered previously and consists of a single unit with eight active sites. The newly revealed type II exists only in M. antarctica KOD 1030 where the unit is separated into two subunits α and β comprised of three and five active sites, respectively. Another newly revealed type III exists in M. alpina AD071 and Dissophora globulifera REB-010B, similar to type II but different in having one more acyl carrier protein domain in the α subunit. This study provides novel insights into the enzymes related to the lipid metabolism, especially the ARA-related Delta5, ELOVL3, and FAS, laying a foundation for genetic engineering of Mortierellaceae to modulate yield in polyunsaturated fatty acids.

First two mitochondrial genomes for the order Filobasidiales reveal novel gene rearrangements and intron dynamics of Tremellomycetes

In the present study, two mitogenomes from the Filobasidium genus were assembled and compared with other Tremellomycetes mitogenomes. The mitogenomes of F. wieringae and F. globisporum both comprised circular DNA molecules, with sizes of 27,861 bp and 71,783 bp, respectively. Comparative mitogenomic analysis revealed that the genetic contents, tRNAs, and codon usages of the two Filobasidium species differed greatly. The sizes of the two Filobasidium mitogenomes varied greatly with the introns being the main factor contributing to mitogenome expansion in F. globisporum. Positive selection was observed in several protein-coding genes (PCGs) in the Agaricomycotina, Pucciniomycotina, and Ustilaginomycotina species, including cob, cox2, nad2, and rps3 genes. Frequent intron loss/gain events were detected to have occurred during the evolution of the Tremellomycetes mitogenomes, and the mitogenomes of 17 species from Agaricomycotina, Pucciniomycotina, and Ustilaginomycotina have undergone large-scale gene rearrangements. Phylogenetic analyses based on Bayesian inference and the maximum likelihood methods using a combined mitochondrial gene set generated identical and well-supported phylogenetic trees, wherein Filobasidium species had close relationships with Trichosporonales species. This study, which is the first report on mitogenomes from the order Filobasidiales, provides a basis for understanding the genomics, evolution, and taxonomy of this important fungal group.

Characterization of the mitogenome of Gongronella sp. w5 reveals substantial variation in Mucoromycota

Gongronella is a genus of fungi in Mucorales (Mucoromycota). Some of its members have important biotechnological applications, but until now, not a single mitogenome has been characterized in Gongronella. Here, we present the complete mitogenome assembly of Gongronella sp. w5, a soil isolate known to interact with plants and several fungi. Its 36,593-bp circular mitogenome encodes the large and small subunit rRNAs, 14 standard mitochondrial proteins, 24 tRNAs, three free-standing ORF proteins, and the RNA subunit of RNase P (rnpB). These genes arrange in an order novel to known fungal mitogenomes. Three group I introns are present in the cob, cox1, and nad5 genes, respectively, and they are probably acquired by horizontal gene transfer. Phylogenetic analysis based on mitochondrion-encoded proteins supports the grouping of Gongronella sp. w5 with Absidia glauca, forming the Cunninghamellaceae clade within Mucoromycota. Gongronella and most other Mucoromycota species are predicted to use the standard genetic code in mitochondrial translation, rather than code 4 assigned by GenBank. A comparison among seven publicly available mitogenomes in Mucoromycota reveals the presence of the same 14 typical protein-coding genes plus rnpB, yet substantial variation in mitogenome size, intron number, gene order, and orientation. In this comparison, the uniqueness of Gongronella is evident from similarly large differences to its closest phylogenetic neighbor, A. glauca. This study promotes our understanding of fungal evolution in Mucoromycota. KEY POINTS: • This study reports the first mitogenome in Gongronella, which presents a novel gene order. • Different Mucoromycota mitogenomes show substantial variation of gene organizations. • Most Mucoromycota species use the standard genetic code to translate mitochondrial genes.

The complete mitochondrial genome of Linnemannia amoeboidea (W. Gams) Vandepol & Bonito (Mortierellales: Mortierellaceae)

The complete mitochondrial genome of Linnemannia amoeboidea (W. Gams) Vandepol & Bonito 2020 (Strain no.: CBS 889.72) was sequenced under the next-generation sequencing platform. It was the second one in the family Mortierellaceae Luerss. 1877. The circular genome was 49,702 bp in size, with a GC content of 20.86%. Gene prediction revealed 15 PCGs, two rRNA genes, 26 tRNA genes, one rnpB gene and seven ORFs. Phylogenetic analyses showed that L. amoeboidea was closely related to Podila verticillate (Linnem.) Vandepol & Bonito 2020.

Mitochondrial Genomics of Six Cacao Pathogens From the Basidiomycete Family Marasmiaceae

Thread blight disease has recently been described as an emerging disease on cacao (Theobroma cacao) in Ghana. In Ghana, thread blight disease is caused by multiple species of the Marasmiaceae family: Marasmius tenuissimus, M. crinis-equi, M. palmivorus, and Marasmiellus scandens. Interestingly, two additional members of the Marasmiaceae; Moniliophthora roreri (frosty pod rot) and Moniliophthora perniciosa (witches’ broom disease), are major pathogens of cacao in the Western hemisphere. It is important to accurately characterize the genetic relationships among these economically important species in support of their disease management. We used data from Illumina NGS-based genome sequencing efforts to study the mitochondrial genomes (mitogenomes) of the four cacao thread blight associated pathogens from Ghana and compared them with published mitogenomes of Mon. roreri and Mon. perniciosa. There is a remarkable interspecies variation in mitogenome size within the six cacao-associated Marasmiaceae species, ranging from 43,121 to 109,103 bp. The differences in genome lengths are primarily due to the number and lengths of introns, differences in intergenic space, and differences in the size and numbers of unidentified ORFs (uORF). Among seven M. tenuissimus mitogenomes sequenced, there is variation in size and sequence pointing to divergent evolution patterns within the species. The intronic regions show a high degree of sequence variation compared to the conserved sequences of the 14 core genes. The intronic ORFs identified, regardless of species, encode GIY-YIG or LAGLIDADG domain-containing homing endonuclease genes. Phylogenetic relationships using the 14 core proteins largely mimic the phylogenetic relationships observed in gene order patterns, grouping M. tenuissimus with M. crinis-equi, and M. palmivorus with Mon. roreri and Mon. perniciosa, leaving Mar. scandens as an outlier. The results from this study provide evidence of independent expansion/contraction events and sequence diversification in each species and establish a foundation for further exploration of the evolutionary trajectory of the fungi in Marasmiaceae family.

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.

Evolutionary Insights Into Two Widespread Ectomycorrhizal Fungi (Pisolithus) From Comparative Analysis of Mitochondrial Genomes

The genus Pisolithus is a group of global ectomycorrhizal fungi. The characterizations of Pisolithus mitochondrial genomes have still been unknown. In the present study, the complete mitogenomes of two Pisolithus species, Pisolithus microcarpus, and Pisolithus tinctorius, were assembled and compared with other Boletales mitogenomes. Both Pisolithus mitogenomes comprised circular DNA molecules with sizes of 43,990 and 44,054 bp, respectively. Comparative mitogenomic analysis showed that the rps3 gene differentiated greatly between Boletales species, and this gene may be subjected to strong pressure of positive selection between some Boletales species. Several plasmid-derived genes and genes with unknown functions were detected in the two Pisolithus mitogenomes, which needs further analysis. The two Pisolithus species show a high degree of collinearity, which may represent the gene arrangement of the ancestors of ectomycorrhizal Boletales species. Frequent intron loss/gain events were detected in Boletales and basidiomycetes, and intron P717 was only detected in P. tinctorius out of the eight Boletales mitogenomes tested. We reconstructed phylogeny of 79 basidiomycetes based on combined mitochondrial gene dataset, and obtained well-supported phylogenetic topologies. This study served as the first report on the mitogenomes of the family Pisolithaceae, which will promote the understanding of the evolution of Pisolithus species.

Recent and Ongoing Horizontal Transfer of Mitochondrial Introns Between Two Fungal Tree Pathogens

Two recently introduced fungal plant pathogens (Ceratocystis lukuohia and Ceratocystis huliohia) are responsible for Rapid ‘ōhi‘a Death (ROD) in Hawai‘i. Despite being sexually incompatible, the two pathogens often co-occur in diseased ‘ōhi‘a sapwood, where genetic interaction is possible. We sequenced and annotated 33 mitochondrial genomes of the two pathogens and related species, and investigated 35 total Ceratocystis mitogenomes. Ten mtDNA regions [one group I intron, seven group II introns, and two autonomous homing endonuclease (HE) genes] were heterogeneously present in C. lukuohia mitogenomes, which were otherwise identical. Molecular surveys with specific primers showed that the 10 regions had uneven geographic distribution amongst populations of C. lukuohia. Conversely, identical orthologs of each region were present in every studied isolate of C. huliohia regardless of geographical origin. Close relatives of C. lukuohia lacked or, rarely, had few and dissimilar orthologs of the 10 regions, whereas most relatives of C. huliohia had identical or nearly identical orthologs. Each region included or worked in tandem with HE genes or reverse transcriptase/maturases that could facilitate interspecific horizontal transfers from intron-minus to intron-plus alleles. These results suggest that the 10 regions originated in C. huliohia and are actively moving to populations of C. lukuohia, perhaps through transient cytoplasmic contact of hyphal tips (anastomosis) in the wound surface of ‘ōhi‘a trees. Such contact would allow for the transfer of mitochondria followed by mitochondrial fusion or cytoplasmic exchange of intron intermediaries, which suggests that further genomic interaction may also exist between the two pathogens.

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.

Recent advances in understanding mitochondrial genome diversity

Ever since its discovery, the double-stranded DNA contained in the mitochondria of eukaryotes has fascinated researchers because of its bacterial endosymbiotic origin, crucial role in encoding subunits of the respiratory complexes, compact nature, and specific inheritance mechanisms. In the last few years, high-throughput sequencing techniques have accelerated the sequencing of mitochondrial genomes (mitogenomes) and uncovered the great diversity of organizations, gene contents, and modes of replication and transcription found in living eukaryotes. Some early divergent lineages of unicellular eukaryotes retain certain synteny and gene content resembling those observed in the genomes of alphaproteobacteria (the inferred closest living group of mitochondria), whereas others adapted to anaerobic environments have drastically reduced or even lost the mitogenome. In the three main multicellular lineages of eukaryotes, mitogenomes have pursued diverse evolutionary trajectories in which different types of molecules (circular versus linear and single versus multipartite), gene structures (with or without self-splicing introns), gene contents, gene orders, genetic codes, and transfer RNA editing mechanisms have been selected. Whereas animals have evolved a rather compact mitochondrial genome between 11 and 50 Kb in length with a highly conserved gene content in bilaterians, plants exhibit large mitochondrial genomes of 66 Kb to 11.3 Mb with large intergenic repetitions prone to recombination, and fungal mitogenomes have intermediate sizes of 12 to 236 Kb.

Conserved novel ORFs in the mitochondrial genome of the ctenophore Beroe forskalii

To date, five ctenophore species' mitochondrial genomes have been sequenced, and each contains open reading frames (ORFs) that if translated have no identifiable orthologs. ORFs with no identifiable orthologs are called unidentified reading frames (URFs). If truly protein-coding, ctenophore mitochondrial URFs represent a little understood path in early-diverging metazoan mitochondrial evolution and metabolism. We sequenced and annotated the mitochondrial genomes of three individuals of the beroid ctenophore Beroe forskalii and found that in addition to sharing the same canonical mitochondrial genes as other ctenophores, the B. forskalii mitochondrial genome contains two URFs. These URFs are conserved among the three individuals but not found in other sequenced species. We developed computational tools called pauvre and cuttlery to determine the likelihood that URFs are protein coding. There is evidence that the two URFs are under negative selection, and a novel Bayesian hypothesis test of trinucleotide frequency shows that the URFs are more similar to known coding genes than noncoding intergenic sequence. Protein structure and function prediction of all ctenophore URFs suggests that they all code for transmembrane transport proteins. These findings, along with the presence of URFs in other sequenced ctenophore mitochondrial genomes, suggest that ctenophores may have uncharacterized transmembrane proteins present in their mitochondria.

Evaluation of a Novel Mitochondrial Pan-Mucorales Marker for the Detection, Identification, Quantification, and Growth Stage Determination of Mucormycetes

Mucormycosis infections are infrequent yet aggressive and serious fungal infections. Early diagnosis of mucormycosis and its discrimination from other fungal infections is required for targeted treatment and more favorable patient outcomes. The majority of the molecular assays use 18 S rDNA. In the current study, we aimed to explore the potential of the mitochondrial rnl (encoding for large-subunit-ribosomal-RNA) gene as a novel molecular marker suitable for research and diagnostics. Rnl was evaluated as a marker for: (1) the Mucorales family, (2) species identification (Rhizopus arrhizus, R. microsporus, Mucor circinelloides, and Lichtheimia species complexes), (3) growth stage, and (4) quantification. Sensitivity, specificity, discriminatory power, the limit of detection (LoD), and cross-reactivity were evaluated. Assays were tested using pure cultures, spiked clinical samples, murine organs, and human paraffin-embedded-tissue (FFPE) samples. Mitochondrial markers were found to be superior to nuclear markers for degraded samples. Rnl outperformed the UMD universal® (Molyzm) marker in FFPE (71.5% positive samples versus 50%). Spiked blood samples highlighted the potential of rnl as a pan-Mucorales screening test. Fungal burden was reproducibly quantified in murine organs using standard curves. Identification of pure cultures gave a perfect (100%) correlation with the detected internal transcribed spacer (ITS) sequence. In conclusion, mitochondrial genes, such as rnl, provide an alternative to the nuclear 18 S rDNA genes and deserve further evaluation.

Complete mitogenome of the entomopathogenic fungus Sporothrix insectorum RCEF 264 and comparative mitogenomics in Ophiostomatales

The fungal order Ophiostomatales contains numerous species important in medical fields, agriculture, and forestry, and several species have had available mitogenome information. The nuclear genome of the entomopathogenic fungus Sporothrix insectorum has been reported, while its mitogenome remains unknown. Herein, we firstly described the mitogenome of S. insectorum RCEF 264 and then compared Ophiostomatales mitogenomes from both interspecific and intraspecific perspectives. The mitogenome of S. insectorum RCEF 264 was 31,454 bp in length, containing typical fungal mitochondrial genes plus rnpB. Four group I introns interrupted rnl and cox1. Phylogenetic analyses confirmed the placement of S. insectorum RCEF 264 in Ophiostomatales. Comparison of mitogenomes among seven Ophiostomatales species revealed conserved gene contents and a high synteny, although there were also some differences among them. Their mitogenomes showed more than two-fold variations (26.6-65.1 kb) in size, with a total of 37 intron insertional loci from 11 genes (1-25 introns per species). The sole intron shared by all species was an rps3-encoding intron in rnl (mL2450), and this intron-based phylogeny was highly consistent with those constructed using mitochondrial/nuclear genes, suggesting convergent evolution of this intron with Ophiostomatales species. The dendrogram based on presence/absence patterns at all intron loci was quite different from those based on mitochondrial/nuclear genes. Comparison of mitogenomes among two to three intraspecific individuals in Ophiostoma novo-ulmi subsp. novo-ulmi and Sporothrix schenckii revealed mitogenome size variations due to single-nucleotide polymorphisms (SNPs) and indels but without fluctuation of intron numbers for each species. This study greatly enhanced our understanding of mitogenome evolution in Ophiostomatales.

Mitovirus and Mitochondrial Coding Sequences from Basal Fungus Entomophthora muscae

Fungi constituting the Entomophthora muscae species complex (members of subphylum Entomophthoromycotina, phylum Zoopagamycota) commonly kill their insect hosts and manipulate host behaviors in the process. In this study, we made use of public transcriptome data to identify and characterize eight new species of mitoviruses associated with several different E. muscae isolates. Mitoviruses are simple RNA viruses that replicate in host mitochondria and are frequently found in more phylogenetically apical fungi (members of subphylum Glomeromyoctina, phylum Mucoromycota, phylum Basidiomycota and phylum Ascomycota) as well as in plants. E. muscae is the first fungus from phylum Zoopagomycota, and thereby the most phylogenetically basal fungus, found to harbor mitoviruses to date. Multiple UGA (Trp) codons are found not only in each of the new mitovirus sequences from E. muscae but also in mitochondrial core-gene coding sequences newly assembled from E. muscae transcriptome data, suggesting that UGA (Trp) is not a rarely used codon in the mitochondria of this fungus. The presence of mitoviruses in these basal fungi has possible implications for the evolution of these viruses.

Green potential of Pleurotus spp. in biotechnology

BACKGROUND

The genus Pleurotus is most exploitable xylotrophic fungi, with valuable biotechnological, medical, and nutritional properties. The relevant features of the representatives of this genus to provide attractive low-cost industrial tools have been reported in numerous studies to resolve the pressure of ecological issues. Additionally, a number of Pleurotus species are highly adaptive, do not require any special conditions for growth, and possess specific resistance to contaminating diseases and pests. The unique properties of Pleurotus species widely used in many environmental technologies, such as organic solid waste recycling, chemical pollutant degradation, and bioethanol production.

METHODOLOGY

The literature study encompasses peer-reviewed journals identified by systematic searches of electronic databases such as Google Scholar, NCBI, Springer, ResearchGate, ScienceDirect, and ISI Web of Knowledge. The search scheme was divided into several steps, as described below.

RESULTS

In this review, we describe studies examining the biotechnological feasibility of Pleurotus spp. to elucidate the importance of this genus for use in green technology. Here, we review areas of application of the genus Pleurotus as a prospective biotechnological tool.

CONCLUSION

The incomplete description of some fungal biochemical pathways emphasises the future research goals for this fungal culture.

Mystique of Phycomyces blakesleeanus is a peculiar mitochondrial genetic element that is highly variable in DNA sequence while subjected to strong negative selection

A DNA region in the mitochondrial genome of the fungus Phycomyces blakesleeanus (Mucorales, Mucoromycota) was characterized in a population of wild-type strains. The region encodes a predicted protein similar to the reverse transcriptases encoded by mitochondrial retroplasmids of Neurospora species and other Sordariomycetes (Ascomycota), but is uncommon in other fungi. DNA sequences of this element, named mystique, are highly variable between the strains, having greater than 2.5% divergence, yet most of the nucleotide differences fall in codon positions that do not change the amino acid sequence. The high proportion of polymorphisms coupled to the rarity of nonsynonymous changes suggests that mystique is subject to counteracting forces of hypermutation and purifying selection. However, while evidence for negative selection may infer that the element provides a fitness benefit, some strains of P. blakesleeanus do not have the element and grow equivalently well as those strains with it. A mechanism to explain the variability between the mystique alleles is proposed, of error-prone replication through an RNA intermediate, reverse transcription and reintegration of the element into the mitochondrial genome.

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.

Advances in Fungal Phylogenomics and Their Impact on Fungal Systematics

In the past decade, advances in next-generation sequencing technologies and bioinformatic pipelines for phylogenomic analysis have led to remarkable progress in fungal systematics and taxonomy. A number of long-standing questions have been addressed using comparative analysis of genome sequence data, resulting in robust multigene phylogenies. These have added to, and often surpassed traditional morphology or single-gene phylogenetic methods. In this chapter, we provide a brief history of fungal systematics and highlight some examples to demonstrate the impact of phylogenomics on this field. We conclude by discussing some of the challenges and promises in fungal biology posed by the ongoing genomics revolution.

The fate of mitochondria after infection of the Mucoralean fungus Absidia glauca by the fusion parasite Parasitella parasitica: comparison of mitochondrial genomes in zygomycetes

Absidia glauca and Parasitella parasitica constitute a versatile experimental system for studying horizontal gene transfer between a mucoralean host and its fusion parasite. The A. glauca chondriome has a length of approximately 63 kb and a GC content of 28%. The chondriome of P. parasitica is larger, 83 kb, and contains 31% GC base pairs. These mtDNAs contain the standard fungal mitochondrial gene set, small and large subunit rRNAs, plus ribonuclease P RNA. Comparing zygomycete chondriomes reveals an unusually high number of homing endonuclease genes in P. parasitica, substantiating the mobility of intron elements independent of host-parasite interactions.

The large (134.9 kb) mitochondrial genome of the glomeromycete Funneliformis mosseae

Funneliformis mosseae is among the most ecologically and economically important glomeromycete species and occurs both in natural and disturbed areas in a wide range of habitats and climates. In this study, we report the sequencing of the complete mitochondrial (mt) genome of F. mosseae isolate FL299 using 454 pyrosequencing and Illumina HiSeq technologies. This mt genome is a full-length circular chromosome of 134,925 bp, placing it among the largest mitochondrial DNAs (mtDNAs) in the fungal kingdom. A comparative analysis with publically available arbuscular mycorrhizal fungal mtDNAs revealed that the mtDNA of F. mosseae FL299 contained a very large number of insertions contributing to its expansion. The gene synteny was completely reshuffled compared to previously published glomeromycotan mtDNAs and several genes were oriented in an anti-sense direction. Furthermore, the presence of different types of introns and insertions in rnl (14 introns) made this gene very distinctive in Glomeromycota. The presence of alternative genetic codes in both initiation (GUG) and termination (UGA) codons was another new feature in this mtDNA compared to previously published glomeromycotan mt genomes. The phylogenetic analysis inferred from the analysis of 14 protein mt genes confirmed the position of the Glomeromycota clade as a sister group of Mortierellomycotina. This mt genome is the largest observed so far in Glomeromycota and the first mt genome within the Funneliformis clade, providing new opportunities to better understand their evolution and to develop molecular markers.

Mechanistic and Structural Studies of Protein-Only RNase P Compared to Ribonucleoproteins Reveal the Two Faces of the Same Enzymatic Activity

RNase P, the essential activity that performs the 5′ maturation of tRNA precursors, can be achieved either by ribonucleoproteins containing a ribozyme present in the three domains of life or by protein-only enzymes called protein-only RNase P (PRORP) that occur in eukaryote nuclei and organelles. A fast growing list of studies has investigated three-dimensional structures and mode of action of PRORP proteins. Results suggest that similar to ribozymes, PRORP proteins have two main domains. A clear functional analogy can be drawn between the specificity domain of the RNase P ribozyme and PRORP pentatricopeptide repeat domain, and between the ribozyme catalytic domain and PRORP N4BP1, YacP-like Nuclease domain. Moreover, both types of enzymes appear to dock with the acceptor arm of tRNA precursors and make specific contacts with the corner of pre-tRNAs. While some clear differences can still be delineated between PRORP and ribonucleoprotein (RNP) RNase P, the two types of enzymes seem to use, fundamentally, the same catalytic mechanism involving two metal ions. The occurrence of PRORP and RNP RNase P represents a remarkable example of convergent evolution. It might be the unique witness of an ongoing replacement of catalytic RNAs by proteins for enzymatic activities.

Complete Mitochondrial DNA Sequence of the Mucoralean Fungus Absidia glauca, a Model for Studying Host-Parasite Interactions

The mitochondrial DNA (mtDNA) ofAbsidia glaucahas been completely sequenced. It is 63,080 bp long, has a G+C content of 28%, and contains the standard fungal gene set.A. glaucais the recipient in a laboratory model for horizontal gene transfer withParasitella parasiticaas a donor of nuclei and mitochondria.

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

Background

Yeasts show remarkable variation in the organization of their mitochondrial genomes, yet there is little experimental data on organellar gene expression outside few model species. Candida albicans is interesting as a human pathogen, and as a representative of a clade that is distant from the model yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe. Unlike them, it encodes seven Complex I subunits in its mtDNA. No experimental data regarding organellar expression were available prior to this study.

Methods

We used high-throughput RNA sequencing and traditional RNA biology techniques to study the mitochondrial transcriptome of C. albicans strains BWP17 and SN148.

Results

The 14 protein-coding genes, two ribosomal RNA genes, and 24 tRNA genes are expressed as eight primary polycistronic transcription units. We also found transcriptional activity in the noncoding regions, and antisense transcripts that could be a part of a regulatory mechanism. The promoter sequence is a variant of the nonanucleotide identified in other yeast mtDNAs, but some of the active promoters show significant departures from the consensus. The primary transcripts are processed by a tRNA punctuation mechanism into the monocistronic and bicistronic mature RNAs. The steady state levels of various mature transcripts exhibit large differences that are a result of posttranscriptional regulation. Transcriptome analysis allowed to precisely annotate the positions of introns in the RNL (2), COB (2) and COX1 (4) genes, as well as to refine the annotation of tRNAs and rRNAs. Comparative study of the mitochondrial genome organization in various Candida species indicates that they undergo shuffling in blocks usually containing 2–3 genes, and that their arrangement in primary transcripts is not conserved. tRNA genes with their associated promoters, as well as GC-rich sequence elements play an important role in these evolutionary events.

Conclusions

The main evolutionary force shaping the mitochondrial genomes of yeasts is the frequent recombination, constantly breaking apart and joining genes into novel primary transcription units. The mitochondrial transcription units are constantly rearranged in evolution shaping the features of gene expression, such as the presence of secondary promoter sites that are inactive, or act as “booster” promoters, simplified transcriptional regulation and reliance on posttranscriptional mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2078-z) contains supplementary material, which is available to authorized users.

Distribution of Ribonucleoprotein and Protein-Only RNase P in Eukarya

RNase P is the endonuclease that removes 5' leader sequences from tRNA precursors. In Eukarya, separate RNase P activities exist in the nucleus and mitochondria/plastids. Although all RNase P enzymes catalyze the same reaction, the different architectures found in Eukarya range from ribonucleoprotein (RNP) enzymes with a catalytic RNA and up to 10 protein subunits to single-subunit protein-only RNase P (PRORP) enzymes. Here, analysis of the phylogenetic distribution of RNP and PRORP enzymes in Eukarya revealed 1) a wealth of novel P RNAs in previously unexplored phylogenetic branches and 2) that PRORP enzymes are more widespread than previously appreciated, found in four of the five eukaryal supergroups, in the nuclei and/or organelles. Intriguingly, the occurrence of RNP RNase P and PRORP seems mutually exclusive in genetic compartments of modern Eukarya. Our comparative analysis provides a global picture of the evolution and diversification of RNase P throughout Eukarya.

Extensive Horizontal Transfer and Homologous Recombination Generate Highly Chimeric Mitochondrial Genomes in Yeast

The frequency of horizontal gene transfer (HGT) in mitochondrial DNA varies substantially. In plants, HGT is relatively common, whereas in animals it appears to be quite rare. It is of considerable importance to understand mitochondrial HGT across the major groups of eukaryotes at a genome-wide level, but so far this has been well studied only in plants. In this study, we generated ten new mitochondrial genome sequences and analyzed 40 mitochondrial genomes from the Saccharomycetaceae to assess the magnitude and nature of mitochondrial HGT in yeasts. We provide evidence for extensive, homologous-recombination-mediated, mitochondrial-to-mitochondrial HGT occurring throughout yeast mitochondrial genomes, leading to genomes that are highly chimeric evolutionarily. This HGT has led to substantial intraspecific polymorphism in both sequence content and sequence divergence, which to our knowledge has not been previously documented in any mitochondrial genome. The unexpectedly high frequency of mitochondrial HGT in yeast may be driven by frequent mitochondrial fusion, relatively low mitochondrial substitution rates and pseudohyphal fusion to produce heterokaryons. These findings suggest that mitochondrial HGT may play an important role in genome evolution of a much broader spectrum of eukaryotes than previously appreciated and that there is a critical need to systematically study the frequency, extent, and importance of mitochondrial HGT across eukaryotes.

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.

Complete Mitochondrial DNA Sequence of the Mucoralean Fusion Parasite Parasitella parasitica

The complete mitochondrial DNA sequence of the Mucor-related fungus Parasitella parasitica has been sequenced. It has a G+C content of 30% and a total length of 83,361 bp. All protein-coding genes normally found in fungi are present in the sequence. A special feature is the remarkably high number of 27 homing endonucleases.

Complete Mitochondrial Genome Sequence of Lichtheimia ramosa (syn. Lichtheimia hongkongensis)

We report the complete mitochondrial genome sequence of Lichtheimia ramosa (syn. Lichtheimia hongkongensis), the first complete mitochondrial DNA sequence of the genus Lichtheimia. This 31.8-kb mitochondrial genome encodes 11 subunits of respiratory chain complexes, 3 ATP synthase subunits, 25 tRNAs, and small and large rRNAs, with the gene order atp9-cox2-atp6-cox3-cox1-nad2-nad3-cob-nad1-nad6-nad5-nad4l-nad4-atp8.

Strikingly bacteria-like and gene-rich mitochondrial genomes throughout jakobid protists

The most bacteria-like mitochondrial genome known is that of the jakobid flagellate Reclinomonas americana NZ. This genome also encodes the largest known gene set among mitochondrial DNAs (mtDNAs), including the RNA subunit of RNase P (transfer RNA processing), a reduced form of transfer-messenger RNA (translational control), and a four-subunit bacteria-like RNA polymerase, which in other eukaryotes is substituted by a nucleus-encoded, single-subunit, phage-like enzyme. Further, protein-coding genes are preceded by potential Shine-Dalgarno translation initiation motifs. Whether similarly ancestral mitochondrial characters also exist in relatives of R. americana NZ is unknown. Here, we report a comparative analysis of nine mtDNAs from five distant jakobid genera: Andalucia, Histiona, Jakoba, Reclinomonas, and Seculamonas. We find that Andalucia godoyi has an even larger mtDNA gene complement than R. americana NZ. The extra genes are rpl35 (a large subunit mitoribosomal protein) and cox15 (involved in cytochrome oxidase assembly), which are nucleus encoded throughout other eukaryotes. Andalucia cox15 is strikingly similar to its homolog in the free-living α-proteobacterium Tistrella mobilis. Similarly, a long, highly conserved gene cluster in jakobid mtDNAs, which is a clear vestige of prokaryotic operons, displays a gene order more closely resembling that in free-living α-proteobacteria than in Rickettsiales species. Although jakobid mtDNAs, overall, are characterized by bacteria-like features, they also display a few remarkably divergent characters, such as 3'-tRNA editing in Seculamonas ecuadoriensis and genome linearization in Jakoba libera. Phylogenetic analysis with mtDNA-encoded proteins strongly supports monophyly of jakobids with Andalucia as the deepest divergence. However, it remains unclear which α-proteobacterial group is the closest mitochondrial relative.

Rapid mitochondrial genome evolution through invasion of mobile elements in two closely related species of arbuscular mycorrhizal fungi

Arbuscular mycorrhizal fungi (AMF) are common and important plant symbionts. They have coenocytic hyphae and form multinucleated spores. The nuclear genome of AMF is polymorphic and its organization is not well understood, which makes the development of reliable molecular markers challenging. In stark contrast, their mitochondrial genome (mtDNA) is homogeneous. To assess the intra- and inter-specific mitochondrial variability in closely related Glomus species, we performed 454 sequencing on total genomic DNA of Glomus sp. isolate DAOM-229456 and we compared its mtDNA with two G. irregulare isolates. We found that the mtDNA of Glomus sp. is homogeneous, identical in gene order and, with respect to the sequences of coding regions, almost identical to G. irregulare. However, certain genomic regions vary substantially, due to insertions/deletions of elements such as introns, mitochondrial plasmid-like DNA polymerase genes and mobile open reading frames. We found no evidence of mitochondrial or cytoplasmic plasmids in Glomus species, and mobile ORFs in Glomus are responsible for the formation of four gene hybrids in atp6, atp9, cox2, and nad3, which are most probably the result of horizontal gene transfer and are expressed at the mRNA level. We found evidence for substantial sequence variation in defined regions of mtDNA, even among closely related isolates with otherwise identical coding gene sequences. This variation makes it possible to design reliable intra- and inter-specific markers.

RNase P branches out from RNP to protein: organelle-triggered diversification?

RNase P is the enzyme that removes 5' leader sequences from precursor tRNAs. Remarkably, in most organisms, RNase P is a ribonucleoprotein particle where the RNA component is responsible for catalysis. In this issue of Genes & Development, Gutmann and colleagues (pp. 1022-1027) report the first organism, Arabidopsis thaliana, to employ protein-only RNase P in both its nucleus and organelles. An intriguing possibility is that replacement of RNase P ribonucleoprotein particles (RNPs) by proteins may have been triggered by the acquisition of organelles.

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.

Group I intron-mediated trans-splicing in mitochondria of Gigaspora rosea and a robust phylogenetic affiliation of arbuscular mycorrhizal fungi with Mortierellales

Gigaspora rosea is a member of the arbuscular mycorrhizal fungi (AMF; Glomeromycota) and a distant relative of Glomus species that are beneficial to plant growth. To allow for a better understanding of Glomeromycota, we have sequenced the mitochondrial DNA of G. rosea. A comparison with Glomus mitochondrial genomes reveals that Glomeromycota undergo insertion and loss of mitochondrial plasmid-related sequences and exhibit considerable variation in introns. The gene order between the two species is almost completely reshuffled. Furthermore, Gigaspora has fragmented cox1 and rns genes, and an unorthodox initiator tRNA that is tailored to decoding frequent UUG initiation codons. For the fragmented cox1 gene, we provide evidence that its RNA is joined via group I-mediated trans-splicing, whereas rns RNA remains in pieces. According to our model, the two cox1 precursor RNA pieces are brought together by flanking cox1 exon sequences that form a group I intron structure, potentially in conjunction with the nad5 intron 3 sequence. Finally, we present analyses that address the controversial phylogenetic association of Glomeromycota within fungi. According to our results, Glomeromycota are not a separate group of paraphyletic zygomycetes but branch together with Mortierellales, potentially also Harpellales.

Genome characterization of the oleaginous fungus Mortierella alpina

Mortierella alpina is an oleaginous fungus which can produce lipids accounting for up to 50% of its dry weight in the form of triacylglycerols. It is used commercially for the production of arachidonic acid. Using a combination of high throughput sequencing and lipid profiling, we have assembled the M. alpina genome, mapped its lipogenesis pathway and determined its major lipid species. The 38.38 Mb M. alpina genome shows a high degree of gene duplications. Approximately 50% of its 12,796 gene models, and 60% of genes in the predicted lipogenesis pathway, belong to multigene families. Notably, M. alpina has 18 lipase genes, of which 11 contain the class 2 lipase domain and may share a similar function. M. alpina's fatty acid synthase is a single polypeptide containing all of the catalytic domains required for fatty acid synthesis from acetyl-CoA and malonyl-CoA, whereas in many fungi this enzyme is comprised of two polypeptides. Major lipids were profiled to confirm the products predicted in the lipogenesis pathway. M. alpina produces a complex mixture of glycerolipids, glycerophospholipids and sphingolipids. In contrast, only two major sterol lipids, desmosterol and 24(28)-methylene-cholesterol, were detected. Phylogenetic analysis based on genes involved in lipid metabolism suggests that oleaginous fungi may have acquired their lipogenic capacity during evolution after the divergence of Ascomycota, Basidiomycota, Chytridiomycota and Mucoromycota. Our study provides the first draft genome and comprehensive lipid profile for M. alpina, and lays the foundation for possible genetic engineering of M. alpina to produce higher levels and diverse contents of dietary lipids.

Of P and Z: mitochondrial tRNA processing enzymes

Mitochondrial tRNAs are generally synthesized as part of polycistronic transcripts. Release of tRNAs from these precursors is thus not only required to produce functional adaptors for translation, but also responsible for the maturation of other mitochondrial RNA species. Cleavage of mitochondrial tRNAs appears to be exclusively accomplished by endonucleases. 5'-end maturation in the mitochondria of different Eukarya is achieved by various kinds of RNase P, representing the full range of diversity found in this enzyme family. While ribonucleoprotein enzymes with RNA components of bacterial-like appearance are found in a few unrelated protists, algae, and fungi, highly degenerate RNAs of dramatic size variability are found in the mitochondria of many fungi. The majority of mitochondrial RNase P enzymes, however, appear to be pure protein enzymes. Human mitochondrial RNase P, the first to be identified and possibly the prototype of all animal mitochondrial RNases P, is composed of three proteins. Homologs of its nuclease subunit MRPP3/PRORP, are also found in plants, algae and several protists, where they are apparently responsible for RNase P activity in mitochondria (and beyond) without the help of extra subunits. The diversity of RNase P enzymes is contrasted by the uniformity of mitochondrial RNases Z, which are responsible for 3'-end processing. Only the long form of RNase Z, which is restricted to eukarya, is found in mitochondria, even when an additional short form is present in the same organism. Mitochondrial tRNA processing thus appears dominated by new, eukaryal inventions rather than bacterial heritage. This article is part of a Special Issue entitled: Mitochondrial Gene Expression.

Multiple recent horizontal transfers of the cox1 intron in Solanaceae and extended co-conversion of flanking exons

Background

The most frequent case of horizontal transfer in plants involves a group I intron in the mitochondrial gene cox1, which has been acquired via some 80 separate plant-to-plant transfer events among 833 diverse angiosperms examined. This homing intron encodes an endonuclease thought to promote the intron's promiscuous behavior. A promising experimental approach to study endonuclease activity and intron transmission involves somatic cell hybridization, which in plants leads to mitochondrial fusion and genome recombination. However, the cox1 intron has not yet been found in the ideal group for plant somatic genetics - the Solanaceae. We therefore undertook an extensive survey of this family to find members with the intron and to learn more about the evolutionary history of this exceptionally mobile genetic element.

Results

Although 409 of the 426 species of Solanaceae examined lack the cox1 intron, it is uniformly present in three phylogenetically disjunct clades. Despite strong overall incongruence of cox1 intron phylogeny with angiosperm phylogeny, two of these clades possess nearly identical intron sequences and are monophyletic in intron phylogeny. These two clades, and possibly the third also, contain a co-conversion tract (CCT) downstream of the intron that is extended relative to all previously recognized CCTs in angiosperm cox1. Re-examination of all published cox1 genes uncovered additional cases of extended co-conversion and identified a rare case of putative intron loss, accompanied by full retention of the CCT.

Conclusions

We infer that the cox1 intron was separately and recently acquired by at least three different lineages of Solanaceae. The striking identity of the intron and CCT from two of these lineages suggests that one of these three intron captures may have occurred by a within-family transfer event. This is consistent with previous evidence that horizontal transfer in plants is biased towards phylogenetically local events. The discovery of extended co-conversion suggests that other cox1 conversions may be longer than realized but obscured by the exceptional conservation of plant mitochondrial sequences. Our findings provide further support for the rampant-transfer model of cox1 intron evolution and recommend the Solanaceae as a model system for the experimental analysis of cox1 intron transfer in plants.