Abstract
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.
1 Introduction
Mitochondria and plastids are semi-autonomous organelles of eukaryotic cells endowed with their own genome and molecular machineries for replication, transcription, and translation. While mitochondria originated from bacterial endosymbionts related to extant α-Proteobacteria, plastids share ancestry with Cyanobacteria. Across eukaryotes, the genomes of mitochondria (mtDNA) and plastids (ptDNA) vary considerably in size, architecture, and coding capacity. MtDNA encodes 5 to 100 genes, which play a role in oxidative phosphorylation, protein synthesis, protein transport and maturation, RNA processing, and in some rare instances, transcription (). PtDNAs encode the same types of genes plus those involved in photosynthesis, summing up to as many as 250 genes (; ).
Sequencing and complete assembly of eukaryotic organelle genomes has become routine and affordable. Yet, despite the relatively small coding capacity of organelle genomes compared to nuclear genomes, identifying genes and subsequently inferring their internal structure (e.g., exon-intron boundaries, herein referred to as ‘gene modelling’) can be challenging and time-consuming. Indeed, whereas organelle genome annotation typically involves automated gene prediction tools, manual analysis and curation by skilled experts are usually necessary to produce accurate results. In the case of mtDNA, the challenges stem from numerous Group I and Group II introns, twintrons (), difficult-to-recognize mini-exons, marginally conserved genes, such as rps3, rnpB, ssrA (; ; ; ), and structurally reduced rRNAs and tRNAs (). Furthermore, intron identification and classification is often only possible using elaborate and manually-refined computational models ().
Several tools have been developed to annotate organelle genomes, including DOGMA (Wyman et al., 2004), MOSAS (), MITOS2 (), Mitofy (), AGORA (), GeSeq (Tillich et al., 2017), and MFannot (). These tools have varying strengths and limitations and are specialized for different groups of organisms. DOGMA and MOSAS were the first to be developed for bilaterian animals. Yet, they produced incomplete gene models requiring substantial expert intervention for completion, and often failed to detect genes outside animals. The more recent tool MITOS2, also tailored to animals, has significantly improved prediction capabilities due to probabilistic inference methods (profile HMMs, CMs) for recognizing protein and ncRNA genes (; ). However, MITOS2 cannot model introns. Although rare in metazoans, introns are present in e.g., corals and sponges (; ; ; ; ). The tools Mitofy, AGORA, and GeSeq were initially optimized for plant organelle genomes which remains their principal strength.
Unfortunately, expert curation of results generated by the above-mentioned tools is not always performed. Consequently, a number of published mitogenomes, even records in the widely used NCBI RefSeq repository (), contain latent errors and deficiencies. Moreover, using such data for novel mitogenome annotations inherently propagates errors and deficiencies, particularly in the case of computational methods that use pairwise similarity searches. The obvious drawbacks of this situation are that researchers who download sequences for various comparative analyses and phylogenomics must curate datasets thoroughly.
A critical component of genome annotators is the algorithms employed for gene-model inference. All the tools mentioned above, except MITOS2 and MFannot, heavily use BLAST-like algorithms to search for sequence similarity with known genes, an approach that often has insufficient sensitivity and precision. A more suited approach involves profile HMMs (), i.e., Hmmsearch for proteins (; ), Cmsearch for ncRNA genes and introns [Infernal, ()] or as an alternative to Cmsearch, ERPIN (; ). Among current organelle annotators, only MITOS2 and MFannot use profile HMMs, Infernal or ERPIN (GeSeq applies HMMs only for prediction refinements).
MFannot, developed in our laboratory, is a comprehensive mitochondrial genome annotation pipeline available as a stand-alone software and a web service. It is optimized for annotating mitogenomes of eukaryotes other than bilaterian animals, but is also capable of annotating plastid genomes (although less accurately than GeSeq). For modelling of intron-containing protein-coding genes, MFannot employs Exonerate () and Hmmsearch (). In contrast, the tools used for identifying ncRNAs (including RNase P RNAs, 5S rRNAs and tmRNAs) are Infernal and ERPIN (; ; Valach et al., 2014). MFannot is unique compared to other annotators for employing probabilistic intron prediction, and its capability to detect mini-exons that are difficult to recognize with other tools.
2 Results and discussion
2.1 Two conceptually distinct approaches to organelle genome annotation
From an algorithmic point of view, organelle genome annotators come in two different flavours, one of which is a next-neighbour-guided annotation (e.g., DOGMA, MOSAS, Mitofy, AGORA, GeSeq), i.e., the identification of genes and genetic elements through comparison with and transfer from well-annotated genomes of very closely related species. For this, BLAST-type similarity-search algorithms (e.g., BLAST and Diamond (; )) are well suited. The advantage of next-neighbour-guided annotation is its very fast computational speed. However, it critically relies on a large and taxonomically broad collection of essentially error-free and complete genome annotations. Hence, the drawback is that this approach is prone to perpetuating occasional annotation errors and gene omissions. The procedure is also less effective for species in which well-annotated mitogenomes from close neighbours are currently unavailable. Therefore, high-quality expert curation of a large number and phylogenetically broad collection of model mitogenomes is a prerequisite this approach.
The alternative to next-neighbour guided annotation is the ab-initio inference of gene models using probabilistic methods (widely employed by MITOS2 and MFannot). For this, sensitive sequence search algorithms (profile HMMs, ERPIN, CMs), based on accurate and evolutionary broad multiple sequence alignments, are employed to identify and model even marginally conserved sequences without requiring annotated genomes of close relatives. This approach comes with longer computation times, notably 1-2 h for a small bilaterian mitogenome by MITOS2 (), but on the other hand, a reasonable 2-10 min for mitogenomes of ~20-200 kbp with the current version of MFannot. (Note that future versions will likely require more execution time, particularly for inferring complete gene structures of rRNA-encoding genes). Significant advantages of ab-initio approaches are high-quality genome annotations for species without close-neighbour information and a moderate requirement for expert curation.
2.2 RNA mapping evidence is of limited value in mitochondrial gene-model prediction
In contrast to nuclear genome annotation approaches, none of the (above-mentioned) organelle genome annotators use RNA data. Furthermore, RNA data are not reported in most organelle genome publications. The reason lies in the particularities of organelle transcript processing and intron splicing. Only a few species produce a high proportion of fully mature mitochondrial transcripts [e.g., the fungus Schizosaccharomyces pombe (; )], to be used to infer intron and gene boundaries from RNA-seq read coverage. In most other species, the transcript landscape is highly complex, particularly for intron-containing genes (e.g., Saccharomyces cerevisiae and most other fungal mtDNAs rich in introns). The complexity of the observed transcript population is due to the stability of intron RNAs that encode proteins [e.g., (; ; Turk et al., 2013)] or form thermodynamically stable RNA structures, as well as slow and partial intron splicing (Figure 1). Together, this leads to highly variable coverage of RNA-seq reads, often spanning exon-intron boundaries. As a result, mapping RNA-seq reads or splice-aligning assembled transcripts to the genome sequence often generates conflicting information that interferes with or misleads gene modelling. In the example shown in Figure 1 [data from ()], expert inference of the gene structure was only possible by sequence comparison with intron-less gene homologs in related species.
Figure 1
2.3 The MFannot annotation procedure
MFannot is written in Perl and was designed to annotate protein-coding and ncRNA genes in mitochondrial and chloroplast genomes. It uses the RNA/intron detection tools described below and is particularly helpful when organelle genomes contain numerous introns. Intron-exon boundaries of protein-coding genes are identified by sequence conservation of exons together with profiles of Group I and II intron-splice sites that, in most instances, can be precisely inferred without transcript data. The output of MFannot lists gene coordinates either in a format that can be directly loaded into NCBI sequence submission tools or in ‘masterfile’ format (a computer-parsable and simultaneously human-readable format developed in-house, with annotations embedded into the sequence as comment lines).
The current annotation procedure (Figure 2) starts with the conceptual translation of the mitogenome (step 1) into Open Reading Frames (ORFs) ≥40 amino acids long, using an in-house tool called Flip (). For this, the user supplies the genetic code. The genetic code in mitochondria varies substantially [e.g., (; )] and should be verified case-by-case, and may require the analysis of potential codon deviations with a dedicated tool ().
Figure 2
In step 2, all predicted ORFs are searched with BLAST against a broad collection of known mitochondrial and plastid proteins, to identify known organellar genes as well as typical intron-encoded ORFs (nucleases, intron maturases, reverse transcriptases) and genes (in particular those encoding RNA and DNA polymerases) derived from the insertion of mitochondrial plasmids into the mitogenome.
In step 3, MFannot uses full-length BLAST matches directly for annotating the respective genes. In addition, MFannot will identify and flag potential frameshifts indicative of either a pseudo-gene or sequencing errors. Partial matches indicate intron insertions, incomplete genes, or trans-splicing (i.e., exons encoded on different DNA strands or at positions that are too distant to account for intron insertions) and are dealt with in subsequent steps.
Notably, the information on the location of bona fide intron-containing genes will be analyzed with Exonerate in steps 4, and annotated in step 5 unless the gene models are incomplete due to the presence of mini-exons or other reasons such as genes in pieces or pseudo-genes (to be resolved in step 7). The specific intron Group (I or II) will be assigned with RNAweasel (), in step 6.
In step 7, splice sites in gene models derived from Exonerate are refined and mini-exons are identified. As Exonerate allows only a single intron splice-site matrix for splice junctions, combining Group I and Group II intron splice patterns results in an indistinctive, almost neutral matrix. Therefore, MFannot checks and potentially adjusts specific intron boundaries in step 7. Furthermore, Exonerate has difficulties with recognizing and modelling small exons shorter than ~4 codons in length (here referred to as ‘mini-exons’), which occur regularly in mitochondrial genes (for more specific information on Exonerate shortcomings, see ‘Future Developments’, below). As a result, small stretches of otherwise highly conserved amino acid positions can be missing in Exonerate alignments, which is also corrected in step 7 of MFannot (Figure 2). Detecting mini-exons is a complex procedure that involves the identification of missing conserved amino acid positions, and scanning of genomic regions that are predicted to contain mini-exons, for best-fitting sequences of the expected size. Technically, this is done by merging each candidate mini-exon with one of the flanking exons and identifying the best candidate from the profile HMM scores for the translated, merged sequences.
In step 8 of protein-gene modelling, MFannot adjusts translation start sites based on matches with profile HMMs, and assesses potential trans-spliced genes (‘genes in pieces’, e.g., (; Wissinger et al., 1991; ; )), frameshifts, and in-frame sequence insertions. Start codon identification is based on (i) the range of start positions in the curated multiple protein alignments (used to create the profile HMMs), and (ii) the presence of potential start codons that fall within or close to this range. If ATG codons are absent, known alternative start codons are considered in the order GTG, TTG and ATA, and if no match is found, MFannot leaves a respective comment in the masterfile record.
Step 9 of the procedure is dedicated to finding genes encoding tRNA and other ncRNA genes, such as rrn5, rnpB, ssrA, rnl and rns, using ERPIN models and CMs (see next section). Finally, step 10 uses profile HMMs to identify less well-conserved proteins among the predicted ORFs.
MFannot stands out from other annotators for listing E-values for intron and gene identification assigned by Infernal and HMM searches. The reporting of E-values is important because it allows the user to make an informed assessment of the results rather than relying on a yes-no answer.
2.4 Identification of well-conserved ncRNAs
MtDNA-encoded ncRNAs that are well conserved at the nucleotide level and easy to identify in most species (with the notable exception of Bilateria) are tRNAs and the small and large subunit rRNAs (rns and rnl, respectively, except for Bilateria and Euglenozoa). MFannot identifies the set of tRNAs with ERPIN models and predicts anticodons and the tRNA identity (amino-acid decoding specificity based on the genetic code that has to be supplied by the user) plus predicted anticodon-codon interactions based on super-wobble rules (for more details see ()). The ERPIN models are sufficiently flexible to recognize unusual structures such as the yeast tRNA(UAG) with an eight-nucleotide anticodon loop that decodes threonine instead of leucine (i.e., in this case, CUN codons have been reassigned from leucine to threonine ()). Note, however, that MFannot does not consider identity elements such as the G3:U70 base pair in the tRNA acceptor stem that is recognized by alanyl-tRNA synthetase (; ), which is the molecular basis for mitochondrial codon reassignment in certain yeast species, from Leu to Thr to Ala (; ). In other words, MFannot will err in rare cases of codon reassignment, both with respect to predicted tRNA identity and protein translation.
As to rRNAs, the small and large rRNAs are easily identified, even without consideration of 2D interactions. In contrast, the 5S rRNA varies substantially in primary sequence, and therefore requires the use of CMs as described below. In most instances, the small subunit rRNA can be precisely mapped with respect to 5′ and 3′ termini, using an ERPIN model for the corresponding regions. Yet, if introns are present, they are just detected without positioning them in a proper gene model (exon/intron structure). For the latter, manual expert work is required, either by sequence comparison with genes from closely related species (preferentially genes that contain no or few introns) or via RNA-seq data that allow positioning exons with reasonable confidence (but see caveats on the use of mitochondrial transcript mapping above). The same applies to identifying rnl, yet with one important exception. The 5′ and 3′ termini of the mitochondrial large subunit rRNA carry marginal conservation at the nucleotide level so that the gene’s ends can only be placed in a window of +/- 50 nt. A more precise prediction could be made by pinpointing a terminal helical structure that occurs in almost all instances (Figure 3). However, in the absence of significant sequence conservation in this helix, it is virtually impossible to build a eukaryote-wide CM or ERPIN model that uses just this base-pairing information (Figure 3). A promising approach to resolving this issue is the construction of clade-specific CMs, as a higher primary sequence conservation can be expected at a shorter evolutionary distance.
Figure 3
2.5 Prediction of less-well conserved ncRNAs
The three additional ncRNAs that are encoded sporadically in organelle genomes across eukaryotes are 5S rRNA (rrn5), tmRNA (ssrA), and RNase P RNAs (rnpB). Many of these genes remain unidentified in GenBank records, in particular, the genes for tmRNA (Figure 4) and RNase P RNA. Still unknown ncRNAs await detection by searching conserved orphan transcripts, followed by comparative phylogenetic modelling with either ERPIN or Infernal.
Among the three mentioned ncRNAs, MFannot identifies 5S rRNA most reliably, based on a CM we developed a few years ago (Valach et al., 2014). The model detects the highly derived structure of Acanthamoeba castellanii that was previously identified and characterized biochemically and proposed to represent a 5S rRNA based on expert structure modelling (). Identifying genes of tmRNA is similarly complex as for 5S rRNA but they can be effectively spotted when searching with a published mitochondrion-specific tmRNA covariance model (). As the tmRNA CM has so far not been integrated into MFannot, users who did not search separately for this gene could have easily missed it, as documented in Figure 4C. This shortcoming will be eliminated in the next version of MFannot. The major remaining challenge is the prediction of RNase P RNA, which according to our preliminary results will require several taxon-specific models (e.g., separate models for yeasts, as well as several other ascomycete and basidiomycete lineages, just to cover the fungal lineage). In fact, it is entirely possible that the reported sporadic presence of less-well conserved and structurally highly variable ncRNAs is in part due to missed identification rather than evolutionary loss. A point in case is fungal mitochondrial RNase P RNA that has an unprecedented structural variability (; ).
Figure 4
2.6 Current limitations and future developments of MFannot
At the time of writing, two major issues remain to be resolved: mini-exon predictions in protein-coding genes, and modelling of intron-containing genes encoding the small and large subunit rRNAs (rns and rnl).
For initial gene identification and gene structure modelling of protein-coding genes, the current implementation of MFannot heavily relies on BLAST and Exonerate. Yet, Exonerate has several shortcomings. Foremost, Exonerate employs pairwise sequence searches for gene-structure prediction rather than a profile built from several sequences that better represents the protein’s diversity. Furthermore, Exonerate uses a dynamic programming algorithm that maximizes a column-wise similarity score (), thus applying a generic substitution probability matrix that is unaware of the actual protein evolution (). Consequently, the tool cannot pinpoint highly conserved or invariable amino-acid positions that are missing in a gene prediction. MFannot compensates, to some degree, for Exonerate’s limited search algorithm and evolutionary model by initially identifying proteins from close neighbours with BLAST and handing them over to Exonerate, which improves the predictions. Still, we have noted that Exonerate sometimes does not find mini-exons (Figure 5A) and, in other instance, predicts an incorrect exon sequence while passing over the valid one (Figure 5B, exons 3 and 7). To resolve this issue, we have developed an MFannot-specific code (Figure 2, step 7) that identifies potential errors and omissions of mini-exons and makes corrections based on the profile HMM approach. While the current mini-exon detection works well in the majority of cases, for unknown reasons, it sometimes fails when gene structures are complex, calling for an investigation of conditions that cause such failures. Pitfalls include the failure to identify short N-terminal exons that are impossible to resolve with the current approach based on Exonerate and our custom mini-exon identification routine, justifying the development of a more robust algorithm. The presence of unidentified Group I introns exacerbates these issues in certain fungal lineages (e.g., Morchellaceae), which is the reason why we plan to improve our ERPIN intron models in the near future.
Figure 5
The second limitation of the current MFannot version is that ERPIN does not model intron-containing rRNA genes. We plan to develop an HMM-based procedure for identifying and modelling rRNA and protein-coding genes. Profile HMMs will be used to scan the genome for conserved regions (as shown in Figure 5). Subsequently, exon boundaries will be refined by identifying the best-fitting intron group, either I or II (as described in Figure 2, step 7). This requires prior intron identification (), which gives hints about the presence of either a mini-exon or, less frequently, twintrons when two introns are seemingly adjacent to each other but not separated by an exon ().
To summarize, in future versions of MFannot, the prediction of protein-coding genes will rely exclusively on HMMs, and that of ncRNA genes on HMMs, CMs or ERPIN.
3 Methods
3.1 Mapping of RNA-seq reads to the mitogenome of C. posadasii
To generate the read coverage map shown in Figure 1, paired-end Illumina RNA-Seq reads (average length of 170 nt per read) were mapped with HISAT2 (; ) to the indexed mitogenome assembly. The resulting uncompressed file was sorted and compressed to BAM format using Samtools (). The output from Samtools was then passed to StringTie () to assemble into contiguous transcripts.
The corresponding commands were:
hisat2-build Coccidioides-posadasii-TGC0611.fna hisat2-index
hisat2 –threads 30 –min-intronlen 20 –max-intronlen 30000 -x hisat2-index -1 RNAseq_mito_R1_001.fastq.gz -2 RNAseq_mito_R1_002.fastq.gz | nice -19 samtools sort -m 1G -@ 30 - | samtools view -bh > RNAseq_mito_R1.bam
stringtie RNAseq_mito_R1.bam -p 20 –rf -o test.gtf -f 0.1 -m 90 -j 3 -c 3 –conservative
The resulting files were loaded into IGV for visualization ().
3.2 Comparison of profile HMM and Exonerate for cox1 gene modelling
Cox1 exon sequences of the Morchella crassipes M10 and Allomyces macrogynus mitogenomes were identified with either Exonerate or profile HMM searches. For the profile HMM, the COX1 protein model currently used by MFannot was searched against the six conceptually translated reading frames using HMMER (v3.3.2) with and without the heuristic filters. For Exonerate, a procedure similar to that of MFannot was used. First, the best candidate protein was identified from the MFannot collection using BLAST+ (v2.13.0). This protein sequence was then searched against the mitogenome using Exonerate (v2.2.0), using parameters and splice models used in MFannot. For M. crassipes, the maximum intron length was set to 20,000 nucleotides to find the complete gene in a single hit.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.
Author contributions
BL, NB, PR, and GB developed the conceptual framework of MFannot, and NB wrote and then updated the Perl code over a 16 years period, and created a containerized version of MFannot and deployed it in CBRAIN (). MS provided informatics and coding support, GB tested the various MFannot versions and provided suggestions for improving the tool, and SP conducted the comparison between Exonerate and profile HMM exon predictions. All authors contributed to the article and approved the submitted version.
Funding
The authors acknowledge generous support by the Natural Sciences and Engineering Research Council of Canada (NSERC grant numbers RGPIN-2014-05286 and RGPIN- 2017-05411; plus an NSERC PGS D scholarship for SP) and by the Fond de Recherche Nature et Technologie, Quebec (FRQNT grant number 2018-PR-206806).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Availability of informatics tools
The code for MFannot, RNAweasel and Flip is available via separate GitHub repositories (; ; ). MFannot is also available within CBRAIN (), an open-source platform that allows using MFannot with various extended options, such as modification of E-value cutoffs. Note that the execution of all tools mentioned in this paper requires a Linux system, and knowledge of command-line program execution. The authors encourage people who wish to participate in the further development of MFannot to contact BFL via email.
References
1
AbboudT. G.ZubaerA.WaiA.HausnerG. (2018). The complete mitochondrial genome of the Dutch elm disease fungus Ophiostoma novo-ulmi subsp. novo-ulmi. Can. J. Microbiol.64 (5), 339–348. doi: 10.1139/cjm-2017-0605
2
AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool. J. Mol. Biol.215 (3), 403–410. doi: 10.1016/S0022-2836(05)80360-2
3
AlversonA. J.WeiX.RiceD. W.SternD. B.BarryK.PalmerJ. D. (2010). Insights into the evolution of mitochondrial genome size from complete sequences of citrullus lanatus and Cucurbita pepo (Cucurbitaceae). Mol. Biol. Evol.27 (6), 1436–1448. doi: 10.1093/molbev/msq029
4
AnzianoP. Q.HansonD. K.MahlerH. R.PerlmanP. S. (1982). Functional domains in introns: trans-acting and cis-acting regions of intron 4 of the cob gene. Cell30 (3), 925–932. doi: 10.1016/0092-8674(82)90297-5
5
BeckN.LangB. F. (2009) RNAweasel, a webserver for identification of mitochondrial, structured RNAs. Available at: https://github.com/BFL-lab/RNAweasel.
6
BeckN.LangB. F. (2010) MFannot, organelle genome annotation websever. Available at: https://github.com/BFL-lab/Mfannot.
7
BerntM.DonathA.JuhlingF.ExternbrinkF.FlorentzC.FritzschG.et al. (2013). MITOS: improved de novo metazoan mitochondrial genome annotation. Mol. Phylogenet. Evol.69 (2), 313–319. doi: 10.1016/j.ympev.2012.08.023
8
BondarenkoN.GlotovaA.NassonovaE.MasharskyA.PolevD.SmirnovA. (2019). The complete mitochondrial genome of Paravannella minima (Amoebozoa, discosea, vannellida). Eur. J. Protistol.68, 80–87. doi: 10.1016/j.ejop.2019.01.005
9
BonenL. (1993). Trans-splicing of pre-mRNA in plants, animals, and protists. FASEB J.7 (1), 40–46. doi: 10.1096/fasebj.7.1.8422973
10
BrossardN.LangB. F.BurgerG. (1996) FLIP, an ORF finder and translator. URL. Available at: https://github.com/BFL-lab/Mfannot.
11
BuchfinkB.XieC.HusonD. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods12 (1), 59–60. doi: 10.1038/nmeth.3176
12
BullerwellC. E.BurgerG.LangB. F. (2000). A novel motif for identifying rps3 homologs in fungal mitochondrial genomes. Trends Biochem. Sci.25 (8), 363–365. doi: 10.1016/S0968-0004(00)01612-1
13
BullerwellC. E.SchnareM. N.GrayM. W. (2003). Discovery and characterization of Acanthamoeba castellanii mitochondrial 5S rRNA. RNA9 (3), 287–292. doi: 10.1261/rna.2170803
14
BurgerG.GrayM. W.ForgetL.LangB. F. (2013). Strikingly bacteria-like and gene-rich mitochondrial genomes throughout jakobid protists. Genome Biol. Evol.5 (2), 418–438. doi: 10.1093/gbe/evt008
15
de Melo TeixeiraM.LangB. F.MatuteD. R.StajichJ. E.BarkerB. (2021). The mitochondrial genomes of the human pathogens Coccidioides immitis and C. posadasii. G3 (Bethesda)11 (7), jkab132. doi: 10.1093/g3journal/jkab132
16
DerelleR.LangB. F. (2011). Rooting the eukaryotic tree with mitochondrial and bacterial proteins. Mol. Biol. Evol.29 (4), 1277–1289. doi: 10.1093/molbev/msr295
17
DerelleR.López-GarcíaP.TimpanoH.MoreiraD. (2016). A phylogenomic framework to study the diversity and evolution of stramenopiles (=Heterokonts). Mol. Biol. Evol.33 (11), 2890–2898. doi: 10.1093/molbev/msw168
18
DerelleR.TorruellaG.KlimesV.BrinkmannH.KimE.VlcekC.et al. (2015). Bacterial proteins pinpoint a single eukaryotic root. Proc. Natl. Acad. Sci. USA112 (7), E693–E699. doi: 10.1073/pnas.1420657112
19
de VriesJ.ArchibaldJ. M. (2018). Plastid genomes. Curr. Biol.28 (8), R336–r337. doi: 10.1016/j.cub.2018.01.027
20
DonathA.JühlingF.Al-ArabM.BernhartS. H.ReinhardtF.StadlerP. F.et al. (2019). Improved annotation of protein-coding genes boundaries in metazoan mitochondrial genomes. Nucleic Acids Res.47 (20), 10543–10552. doi: 10.1093/nar/gkz833
21
EddyS. R. (1995). Multiple alignment using hidden Markov models. Proc. Int. Conf. Intell. Syst. Mol. Biol.3, 114–120.
22
EddyS. R. (2009). A new generation of homology search tools based on probabilistic inference. Genome Inform23 (1), 205–211. doi: 10.1142/9781848165632_0019
23
EddyS. R. (2011). Accelerated profile HMM searches. PloS Comput. Biol.7 (10), e1002195. doi: 10.1371/journal.pcbi.1002195
24
GalindoL. J.ProkinaK.TorruellaG.López-GarcíaP.MoreiraD. (2023). Maturases and group II introns in the mitochondrial genomes of the deepest jakobid branch. Genome Biol. Evol.15 (4), evad058. doi: 10.1093/gbe/evad058
25
GautheretD.LambertA. (2001). Direct RNA motif definition and identification from multiple sequence alignments using secondary structure profiles. J. Mol. Biol.313 (5), 1003–1011. doi: 10.1006/jmbi.2001.5102
26
GiegéR.SisslerM.FlorentzC. (1998). Universal rules and idiosyncratic features in tRNA identity. Nucleic Acids Res.26 (22), 5017–5035.
27
HafezM.BurgerG.SteinbergS. V.LangB. F. (2013a). A second eukaryotic group with mitochondrion-encoded tmRNA: in silico identification and experimental confirmation. RNA Biol.10 (7), 1117–1124. doi: 10.4161/rna.25376
28
HafezM.MajerA.SethuramanJ.RudskiS. M.MichelF.HausnerG. (2013b). The mtDNA rns gene landscape in the ophiostomatales and other fungal taxa: twintrons, introns, and intron-encoded proteins. Fungal Genet. Biol.53, 71–83. doi: 10.1016/j.fgb.2013.01.005
29
HansonD. K.LambM. R.MahlerH. R.PerlmanP. S. (1982). Evidence for translated intervening sequences in the mitochondrial genome of Saccharomyces cerevisiae. J. Biol. Chem.257 (6), 3218–3224. doi: 10.1016/S0021-9258(19)81098-0
30
HeissA. A.WarringS. D.LukacsK.FavateJ.YangA.GyaltshenY.et al. (2021). Description of Imasa heleensis, gen. nov., sp. nov. (Imasidae, fam. nov.), a deep-branching marine malawimonad and possible key taxon in understanding early eukaryotic evolution. J. Eukaryot Microbiol.68 (2), e12837. doi: 10.1111/jeu.12837
31
HenikoffS.HenikoffJ. G. (1992). Amino acid substitution matrices from protein blocks. Proc. Natl. Acad. Sci. U.S.A.89 (22), 10915–10919. doi: 10.1073/pnas.89.22.10915
32
JanouškovecJ.TikhonenkovD. V.BurkiF.HoweA. T.RohwerF. L.MylnikovA. P.et al. (2017). A new lineage of eukaryotes illuminates early mitochondrial genome reduction. Curr. Biol.27 (23), 3717–3724.e3715. doi: 10.1016/j.cub.2017.10.051
33
JungJ.KimJ. I.JeongY. S.YiG. (2018). AGORA: organellar genome annotation from the amino acid and nucleotide references. Bioinformatics34 (15), 2661–2663. doi: 10.1093/bioinformatics/bty196
34
KimD.LangmeadB.SalzbergS. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods12 (4), 357–360. doi: 10.1038/nmeth.3317
35
KimD.PaggiJ. M.ParkC.BennettC.SalzbergS. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol.37 (8), 907–915. doi: 10.1038/s41587-019-0201-4
36
LambertA.FontaineJ. F.LegendreM.LeclercF.PermalE.MajorF.et al. (2004). The ERPIN server: an interface to profile-based RNA motif identification. Nucleic Acids Res.32 (Web Server issue), W160–W165. doi: 10.1093/nar/gkh418
37
LangB. F.BurgerG. (2012). ''Mitochondrial and eukaryotic origins: a critical review' in Advances in Botanical Research63, 1-20 doi: 10.1016/B978-0-12-394279-1.00001-6
38
LangB. F.BurgerG.O'KellyC. J.CedergrenR.GoldingG. B.LemieuxC.et al. (1997). An ancestral mitochondrial DNA resembling a eubacterial genome in miniature. Nature387 (6632), 493–497. doi: 10.1038/387493a0
39
LangB. F.CedergrenR.GrayM. W. (1987). The mitochondrial genome of the fission yeast, Schizosaccharomyces pombe. sequence of the large-subunit ribosomal RNA gene, comparison of potential secondary structure in fungal mitochondrial large-subunit rRNAs and evolutionary considerations. Eur. J. Biochem.169 (3), 527–537.
40
LangB. F.LaforestM. J.BurgerG. (2007). Mitochondrial introns: a critical view. Trends Genet.23, 119–125. doi: 10.1016/j.tig.2007.01.006
41
LangB. F.LavrovD.BeckN.SteinbergV. (2011). “Mitochondrial tRNA structure, identity and evolution of the genetic code”. in Organelle genetics. Ed. BullerwellC. E. (Berlin, Heidelberg, 2012: Springer), 431–474.
42
LavrovD.LangB. F. (2013). ''Mitochondrial genomes in unicellular relatives of animals (ASBMB)''. in Molecular Life Sciences: An Encyclopedic Reference. Eds. WellsR. D.BondJ. S.KlinmanJ.MastersB. S. S.BellE. (New York : Springer) doi: 10.1007/978-1-4614-6436-5_178-2.
43
LavrovD. V.PettW. (2016). Animal mitochondrial DNA as we do not know it: mt-genome organization and evolution in nonbilaterian lineages. Genome Biol. Evol.8 (9), 2896–2913. doi: 10.1093/gbe/evw195
44
LiH.HandsakerB.WysokerA.FennellT.RuanJ.HomerN.et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics25 (16), 2078–2079. doi: 10.1093/bioinformatics/btp352
45
LingJ.DaoudR.LajoieM. J.ChurchG. M.SollD.LangB. F. (2014). Natural reassignment of CUU and CUA sense codons to alanine in Ashbya mitochondria. Nucleic Acids Res.42 (1), 499–508. doi: 10.1093/nar/gkt842
46
Muñoz-GómezS. A.Mejía-FrancoF. G.DurninK.ColpM.GrisdaleC. J.ArchibaldJ. M.et al. (2017). The new red algal subphylum proteorhodophytina comprises the largest and most divergent plastid genomes known. Curr. Biol.27 (11), 1677–1684.e1674. doi: 10.1016/j.cub.2017.04.054
47
Musier-ForsythK.UsmanN.ScaringeS.DoudnaJ.GreenR.SchimmelP. (1991). Specificity for aminoacylation of an RNA helix: an unpaired, exocyclic amino group in the minor groove. Science253 (5021), 784–786. doi: 10.1126/science.1876835
48
NadimiM.BeaudetD.ForgetL.HijriM.LangB. F. (2012). Group I intron-mediated trans-splicing in mitochondria of Gigaspora rosea and a robust phylogenetic affiliation of arbuscular mycorrhizal fungi with mortierellales. Mol. Biol. Evol.29 (9), 2199–2210. doi: 10.1093/molbev/mss088
49
NawrockiE. P.EddyS. R. (2013). Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics29 (22), 2933–2935. doi: 10.1093/bioinformatics/btt509
50
NoutahiE.CalderonV.BlanchetteM.LangF. B.El-MabroukN. (2017). CoreTracker: accurate codon reassignment prediction, applied to mitochondrial genomes. Bioinformatics33 (21), 3331–3339. doi: 10.1093/bioinformatics/btx421
51
OkimotoR.MacfarlaneJ. L.WolstenholmeD. R. (1994). The mitochondrial ribosomal RNA genes of the nematodes Caenorhabditis elegans and Ascaris suum: consensus secondary-structure models and conserved nucleotide sets for phylogenetic analysis. J. Mol. Evol.39 (6), 598–613. doi: 10.1007/BF00160405
52
Pereira de SouzaA.JubierM. F.DelcherE.LancelinD.LejeuneB. (1991). A trans-splicing model for the expression of the tripartite nad5 gene in wheat and maize mitochondria. Plant Cell3 (12), 1363–1378. doi: 10.2307/3869315
53
PerteaM.PerteaG. M.AntonescuC. M.ChangT. C.MendellJ. T.SalzbergS. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol.33 (3), 290–295. doi: 10.1038/nbt.3122
54
PrinceS.MunozC.Filion-BienvenueF.RiouxP.SarrasinM.LangB. F. (2022). Refining mitochondrial intron classification with ERPIN: identification based on conservation of sequence plus secondary structure motifs. Front. Microbiol.13. doi: 10.3389/fmicb.2022.866187
55
PruittK. D.TatusovaT.MaglottD. R. (2003). NCBI reference sequence project: update and current status. Nucleic Acids Res.31 (1), 34–37. doi: 10.1093/nar/gkg111
56
SchaferB.HansenM.LangB. F. (2005). Transcription and RNA-processing in fission yeast mitochondria. RNA11 (5), 785–795. doi: 10.1261/rna.7252205
57
SeifE. R.ForgetL.MartinN. C.LangB. F. (2003). Mitochondrial RNase p RNAs in ascomycete fungi: lineage-specific variations in RNA secondary structure. RNA9 (9), 1073–1083. doi: 10.1261/rna.5880403
58
SeifE.LeighJ.LiuY.RoewerI.ForgetL.LangB. F. (2005). Comparative mitochondrial genomics in zygomycetes: bacteria-like RNase p RNAs, mobile elements and a close source of the group I intron invasion in angiosperms. Nucleic Acids Res.33 (2), 734–744. doi: 10.1093/nar/gki199
59
ShangJ.YangY.WuL.ZouM.HuangY. (2018). The S. pombe mitochondrial transcriptome. RNA24 (9), 1241–1254. doi: 10.1261/rna.064477.117
60
SheffieldN. C.HiattK. D.ValentineM. C.SongH.WhitingM. F. (2010). Mitochondrial genomics in Orthoptera using MOSAS. Mitochondrial DNA21 (3-4), 87–104. doi: 10.3109/19401736.2010.500812
61
SherifT.RiouxP.RousseauM. E.KassisN.BeckN.AdalatR.et al. (2014). CBRAIN: a web-based, distributed computing platform for collaborative neuroimaging research. Front. Neuroinform8. doi: 10.3389/fninf.2014.00054
62
ShiratoriT.IshidaK. I. (2016). A new heterotrophic cryptomonad: Hemiarma marina n. g., n. sp. J. Eukaryot Microbiol.63 (6), 804–812. doi: 10.1111/jeu.12327
63
SlaterG. S. C.BirneyE. (2005). Automated generation of heuristics for biological sequence comparison. BMC Bioinf.6, 31. doi: 10.1186/1471-2105-6-31
64
StrassertJ. F.TikhonenkovD. V.PombertJ. F.KoliskoM.TaiV.MylnikovA. P.et al. (2016). Moramonas marocensis gen. nov., sp. nov.: a jakobid flagellate isolated from desert soil with a bacteria-like, but bloated mitochondrial genome. Open Biol.6 (2), 150239. doi: 10.1098/rsob.150239
65
SuD.LiebermanA.LangB. F.SimonovicM.SollD.LingJ. (2011). An unusual tRNAThr derived from tRNAHis reassigns in yeast mitochondria the CUN codons to threonine. Nucleic Acids Res.39 (11), 4866–4874. doi: 10.1093/nar/gkr073
66
TekleY. I.WangF.WoodF. C.AndersonO. R.SmirnovA. (2022). New insights on the evolutionary relationships between the major lineages of amoebozoa. Sci. Rep.12 (1), 11173. doi: 10.1038/s41598-022-15372-7
67
ThorvaldsdóttirH.RobinsonJ. T.MesirovJ. P. (2013). Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform.14 (2), 178–192. doi: 10.1093/bib/bbs017
68
TillichM.LehwarkP.PellizzerT.Ulbricht-JonesE. S.FischerA.BockR.et al. (2017). GeSeq - versatile and accurate annotation of organelle genomes. Nucleic Acids Res.45 (W1), W6–w11. doi: 10.1093/nar/gkx391
69
TurkE. M.DasV.SeibertR. D.AndrulisE. D. (2013). The mitochondrial RNA landscape of Saccharomyces cerevisiae. PloS One8 (10), e78105. doi: 10.1371/journal.pone.0078105
70
ValachM.BurgerG.GrayM. W.LangB. F. (2014). Widespread occurrence of organelle genome-encoded 5S rRNAs including permuted molecules. Nucleic Acids Res.42 (22), 13764–13777. doi: 10.1093/nar/gku1266
71
WissingerB.SchusterW.BrennickeA. (1991). Trans splicing in Oenothera mitochondria: nad1 mRNAs are edited in exon and trans-splicing group II intron sequences. Cell65 (3), 473–482. doi: 10.1016/0092-8674(91)90465-B
72
WymanS. K.JansenR. K.BooreJ. L. (2004). Automatic annotation of organellar genomes with DOGMA. Bioinformatics20 (17), 3252–3255. doi: 10.1093/bioinformatics/bth352
73
YabukiA.GyaltshenY.HeissA. A.FujikuraK.KimE. (2018). Ophirina amphinema n. gen., n. sp., a new deeply branching discobid with phylogenetic affinity to jakobids. Sci. Rep.8 (1), 16219. doi: 10.1038/s41598-018-34504-6
74
YazakiE.YabukiA.ImaizumiA.KumeK.HashimotoT.InagakiY. (2022). The closest lineage of archaeplastida is revealed by phylogenomics analyses that include Microheliella maris. Open Biol.12 (4), 210376. doi: 10.1098/rsob.210376
75
ZubaerA.WaiA.HausnerG. (2018). The mitochondrial genome of Endoconidiophora resinifera is intron rich. Sci. Rep.8 (1), 17591. doi: 10.1038/s41598-018-35926-y
Summary
Keywords
mitochondrial DNA, genome annotation, introns, profile HMMs, ERPIN, covariance models, RNA structure
Citation
Lang BF, Beck N, Prince S, Sarrasin M, Rioux P and Burger G (2023) Mitochondrial genome annotation with MFannot: a critical analysis of gene identification and gene model prediction. Front. Plant Sci. 14:1222186. doi: 10.3389/fpls.2023.1222186
Received
13 May 2023
Accepted
15 June 2023
Published
04 July 2023
Volume
14 - 2023
Edited by
Jeffrey P. Mower, University of Nebraska-Lincoln, United States
Reviewed by
Georg Hausner, University of Manitoba, Canada; Weilong Hao, Wayne State University, United States
Updates
Copyright
© 2023 Lang, Beck, Prince, Sarrasin, Rioux and Burger.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: B. Franz Lang, Franz.Lang@UMontreal.ca
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.