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.

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 (Lavrov and Lang, 2013). PtDNAs encode the same types of genes plus those involved in photosynthesis, summing up to as many as 250 genes (Muñoz-Gómez et al., 2017; de Vries and Archibald, 2018).

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 (Hafez et al., 2013b), difficult-to-recognize mini-exons, marginally conserved genes, such as rps3, rnpB, ssrA (Bullerwell et al., 2000; Seif et al., 2003; Hafez et al., 2013a; Donath et al., 2019), and structurally reduced rRNAs and tRNAs (Okimoto et al., 1994). Furthermore, intron identification and classification is often only possible using elaborate and manually-refined computational models (Prince et al., 2022).

Several tools have been developed to annotate organelle genomes, including DOGMA (Wyman et al., 2004), MOSAS (Sheffield et al., 2010), MITOS2 (Donath et al., 2019), Mitofy (Alverson et al., 2010), AGORA (Jung et al., 2018), GeSeq (Tillich et al., 2017), and MFannot (Beck and Lang, 2010). 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 (Bernt et al., 2013; Donath et al., 2019). However, MITOS2 cannot model introns. Although rare in metazoans, introns are present in e.g., corals and sponges (Lang and Burger, 2012; Bernt et al., 2013; Lavrov and Pett, 2016; Donath et al., 2019; Prince et al., 2022). 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 (Pruitt et al., 2003), 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 (Eddy, 1995), i.e., Hmmsearch for proteins (Eddy, 2009; Eddy, 2011), Cmsearch for ncRNA genes and introns [Infernal, (Nawrocki and Eddy, 2013)] or as an alternative to Cmsearch, ERPIN (Lambert et al., 2004; Prince et al., 2022). 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 (Slater and Birney, 2005) and Hmmsearch (Eddy, 2009). In contrast, the tools used for identifying ncRNAs (including RNase P RNAs, 5S rRNAs and tmRNAs) are Infernal and ERPIN (Lang et al., 2007; Hafez et al., 2013a; 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 (Altschul et al., 1990; Buchfink et al., 2015)) 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 (Bernt et al., 2013), 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 (Schafer et al., 2005; Shang et al., 2018)], 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., (Anziano et al., 1982; Hanson et al., 1982; 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 (de Melo Teixeira et al., 2021)], expert inference of the gene structure was only possible by sequence comparison with intron-less gene homologs in related species.

FIGURE 1

Figure 1 Mapping of RNA-seq data to the mitochondrial genome of Coccidioides posadasii. The figure depicts a 2,579 nt window over the coordinates of exons 4 (6172-6736) and 5 (8458-8504) of the mitochondrial gene encoding the large subunit rRNA (rnl) of C. posadasii (see (de Melo Teixeira et al., 2021) for details on genome sequencing) as determined by rRNA sequence conservation and structural modelling of exons 4 and 5 plus adjacent introns. The read-mapping coverage distribution (grey) is shown below the annotated exons. Red and blue arches indicate forward and reverse split reads, respectively. The bottom part of the upper panel depicts the splice junctions inferred from the RNA-Seq read to genome mapping. The coverage profile reveals a substantial number of reads that map within intron 4 due to stable transcripts encoding intron ORFs and an elevated level of un-spliced RNA precursors, which is common in mitochondria. The two distinct splice junctions inferred at the 5’ and 3’ end of intron 4 are supported in roughly equal proportions by mapped reads. The bottom track depicts exons inferred from three distinct transcripts (marked 1, 2 and 3) that were reconstructed from the mapped reads using StringTie. The two lower panels show read mapping in higher resolution. The vertical dotted line indicates the predicted precise splice junction.

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 (Brossard et al., 1996). For this, the user supplies the genetic code. The genetic code in mitochondria varies substantially [e.g., (Su et al., 2011; Ling et al., 2014)] and should be verified case-by-case, and may require the analysis of potential codon deviations with a dedicated tool (Noutahi et al., 2017).

FIGURE 2

Figure 2 Flowchart of the MFannot annotation procedure. The figure summarizes the analytical steps and the external tools employed (indicated to the right). The external tools are Flip (conceptual translation), BLAST (fast sequence similarity search), Exonerate (annotation of genes with exons/intron structure), HMMER (high sensitivity sequence matching), and RNAweasel (introns Group I and II, ncRNAs). When no specific tool is mentioned, the corresponding step is executed by MFannot-specific code. The RNAweasel tool (step 6 and 9) uses ERPIN (and more recently, Infernal) as a search engine. ERPIN’s search algorithm is based on RNA secondary structure profiles computed from RNA sequence alignments plus user-defined secondary-structure information as an input. Much of ERPIN’s efficiency stems from the definition of precisely delimited structural elements that can be searched individually or in combination using a particular search order (‘search strategy’). In general, ERPIN (Gautheret and Lambert, 2001) is the second-most sensitive search algorithm for structured RNAs [following the covariance-based Infernal program (Nawrocki and Eddy, 2013)]. Still, it distinguishes more reliably certain mitochondrial Group I and II introns.

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 (Beck and Lang, 2009), 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., (Pereira de Souza et al., 1991; Wissinger et al., 1991; Bonen, 1993; Nadimi et al., 2012)), 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 (Lang et al., 2011)). 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 (Su et al., 2011)). 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 (Musier-Forsyth et al., 1991; Giegé et al., 1998), which is the molecular basis for mitochondrial codon reassignment in certain yeast species, from Leu to Thr to Ala (Su et al., 2011; Ling et al., 2014). 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

Figure 3 Helical RNA structures joining the 5′ and 3′-terminal regions of mitochondrial large subunit rRNAs. In the three given examples, the RNA termini have been determined experimentally at single-nucleotide resolution (Lang et al., 1987; Lang et al., 1997). The RNA termini are either adjacent to, or within a few nucleotide positions outside the helical region. Note the lack of significantly conserved sequence motifs and the variability of helix shape and length, which renders in-silico identification of the termini difficult.

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 (Bullerwell et al., 2003). 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 (Hafez et al., 2013a). 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 (Seif et al., 2003; Seif et al., 2005).

FIGURE 4

Figure 4 Identification of mitochondrial tmRNA (ssrA) genes across eukaryotes. (A) Examples of mitochondrial tmRNA 2D structures. The helices P1, P2 and P3 are indicated in the A. godoyi 2D structure, in accordance with the sequence alignment in (C). The broken line in the 2D structure marks the sites of endonucleolytic processing that give rise to mature tRNA-like 5′ and 3′ ends. The 3′ end undergoes further modification by adding a CCA [51], as indicated in the figure. The invariant G-U pair in helix 3 (boxed) is required to recognize the tmRNA by tRNA synthase alanine. The GAC motif that is invariant across tmRNAs is marked in red. (B) Schematic eukaryotic tree (consensus derived from topologies in (Derelle and Lang, 2011; Derelle et al., 2015; Derelle et al., 2016; Janouškovec et al., 2017; Heiss et al., 2021; Tekle et al., 2022; Yazaki et al., 2022)). Names of species containing a mt ssrA gene are marked in mauve, and the three species with 2D structures shown in (A) are boxed. Mt ssrA genes were identified by an exhaustive search of published mitogenomes with a previously published covariance model based on genes from jakobids plus stramenopiles (Hafez et al., 2013a). The search identifies homologs in four eukaryotic lineages formerly thought to lack mt ssrA (cryptophytes, haptophytes, malawimonads, and amoebozoans), as well as in numerous additional Oomycota (represented in the tree by a single species, Phytophthora infestans). (C) Structured alignment of tmRNA sequences. Brown text colour marks the six species for which an ssrA annotation in GenBank records is lacking (Shiratori and Ishida, 2016; Strassert et al., 2016; Janouškovec et al., 2017; Yabuki et al., 2018; Bondarenko et al., 2019; Heiss et al., 2021). For further details on previous results and the structural annotation, see (Hafez et al., 2013a). Note that a new version of MFannot (available in June 2023) includes tmRNA-CM searches. During the preparation of the manuscript, we noted a recent publication on two additional jakobid mitogenomes that does not mention the presence of ssrA genes (Galindo et al., 2023). As expected due to its almost ubiquitous presence across jakobids (the only current exception is J. bahamiensis), ssrA genes are present in both Agogonia voluta and Ophirina chinija mtDNAs (between rps2 and trnS(gcu); E-values of 6.7e-11 and 4.4e-13, respectively).

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 (Slater and Birney, 2005), thus applying a generic substitution probability matrix that is unaware of the actual protein evolution (Henikoff and Henikoff, 1992). 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

Figure 5 Comparison of exon identification with either HMM or Exonerate. Exons in cox1 genes (encoding cytochrome c oxidoreductase subunit 1) from (A) Allomyces macrogynus (GenBank U41288.1) and (B) Morchella crassipes (GenBank MN542893.2) were identified with either Exonerate (using a single Cox1 protein for the comparison), or HMMER (search with a Cox1 profile HMM in conceptually translated proteins of mitochondrial genes). Lines labelled ‘Sequence’ display the expected, expert-curated protein sequence. Intron positions are marked by a ‘^’. ‘PP’ indicates the posterior probability values of HMM searches (‘8’=0.8; ‘9’=0.9; ‘*’=1.0). When codons are split by an intron, the sequence of the split nucleotide is shown above the ‘Sequence’ line. Exons identified with either the default Exonerate settings or with an Hmmsearch at a reporting threshold of <= 1e-5 are shown in green. Additional exons (mauve sequence) were found with low-stringency parameters: Exonerate with the option –proteinwordlen 3 detects the A. macrogynus exon 5 and the M. crassipes exon 4. An HMM search using a cut-off e-value between 1e-5 and 0.1 detects the A. macrogynus exon 5. Note that cox1 exons 3 and 7 of M. crassipes are missed by both approaches. Exonerate at these low-stringency settings erroneously adds a single cysteine to exon 3 and an incorrect sequence for exon 7 which is in phase with exon 6. The incorrect Exonerate predictions are shown in red.

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 (Prince et al., 2022), 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 (Hafez et al., 2013b).

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 (Kim et al., 2015; Kim et al., 2019) to the indexed mitogenome assembly. The resulting uncompressed file was sorted and compressed to BAM format using Samtools (Li et al., 2009). The output from Samtools was then passed to StringTie (Pertea et al., 2015) to assemble into contiguous transcripts.

The corresponding commands were:

1. hisat2-build Coccidioides-posadasii-TGC0611.fna hisat2-index

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

3. 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 (Thorvaldsdóttir et al., 2013).

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.

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 (Sherif et al., 2014). 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.

References

Abboud, T. G., Zubaer, A., Wai, A., Hausner, G. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Altschul, S. F., Gish, W., Miller, W., Myers, E. W., Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215 (3), 403–410. doi: 10.1016/S0022-2836(05)80360-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Alverson, A. J., Wei, X., Rice, D. W., Stern, D. B., Barry, K., Palmer, J. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Anziano, P. Q., Hanson, D. K., Mahler, H. R., Perlman, P. S. (1982). Functional domains in introns: trans-acting and cis-acting regions of intron 4 of the cob gene. Cell 30 (3), 925–932. doi: 10.1016/0092-8674(82)90297-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Beck, N., Lang, B. F. (2009) RNAweasel, a webserver for identification of mitochondrial, structured RNAs. Available at: https://github.com/BFL-lab/RNAweasel.

Google Scholar

Beck, N., Lang, B. F. (2010) MFannot, organelle genome annotation websever. Available at: https://github.com/BFL-lab/Mfannot.

Google Scholar

Bernt, M., Donath, A., Juhling, F., Externbrink, F., Florentz, C., Fritzsch, G., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Bondarenko, N., Glotova, A., Nassonova, E., Masharsky, A., Polev, D., Smirnov, A. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonen, L. (1993). Trans-splicing of pre-mRNA in plants, animals, and protists. FASEB J. 7 (1), 40–46. doi: 10.1096/fasebj.7.1.8422973

PubMed Abstract | CrossRef Full Text | Google Scholar

Brossard, N., Lang, B. F., Burger, G. (1996) FLIP, an ORF finder and translator. URL. Available at: https://github.com/BFL-lab/Mfannot.

Google Scholar

Buchfink, B., Xie, C., Huson, D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12 (1), 59–60. doi: 10.1038/nmeth.3176

PubMed Abstract | CrossRef Full Text | Google Scholar

Bullerwell, C. E., Burger, G., Lang, B. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Bullerwell, C. E., Schnare, M. N., Gray, M. W. (2003). Discovery and characterization of Acanthamoeba castellanii mitochondrial 5S rRNA. RNA 9 (3), 287–292. doi: 10.1261/rna.2170803

PubMed Abstract | CrossRef Full Text | Google Scholar

Burger, G., Gray, M. W., Forget, L., Lang, B. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

de Melo Teixeira, M., Lang, B. F., Matute, D. R., Stajich, J. E., Barker, B. (2021). The mitochondrial genomes of the human pathogens Coccidioides immitis and C. posadasii. G3 (Bethesda) 11 (7), jkab132. doi: 10.1093/g3journal/jkab132

PubMed Abstract | CrossRef Full Text | Google Scholar

Derelle, R., Lang, B. F. (2011). Rooting the eukaryotic tree with mitochondrial and bacterial proteins. Mol. Biol. Evol. 29 (4), 1277–1289. doi: 10.1093/molbev/msr295

PubMed Abstract | CrossRef Full Text | Google Scholar

Derelle, R., López-García, P., Timpano, H., Moreira, D. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Derelle, R., Torruella, G., Klimes, V., Brinkmann, H., Kim, E., Vlcek, C., et al. (2015). Bacterial proteins pinpoint a single eukaryotic root. Proc. Natl. Acad. Sci. USA 112 (7), E693–E699. doi: 10.1073/pnas.1420657112

CrossRef Full Text | Google Scholar

de Vries, J., Archibald, J. M. (2018). Plastid genomes. Curr. Biol. 28 (8), R336–r337. doi: 10.1016/j.cub.2018.01.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Donath, A., Jühling, F., Al-Arab, M., Bernhart, S. H., Reinhardt, F., Stadler, P. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Eddy, S. R. (1995). Multiple alignment using hidden Markov models. Proc. Int. Conf. Intell. Syst. Mol. Biol. 3, 114–120.

PubMed Abstract | Google Scholar

Eddy, S. R. (2009). A new generation of homology search tools based on probabilistic inference. Genome Inform 23 (1), 205–211. doi: 10.1142/9781848165632_0019

PubMed Abstract | CrossRef Full Text | Google Scholar

Eddy, S. R. (2011). Accelerated profile HMM searches. PloS Comput. Biol. 7 (10), e1002195. doi: 10.1371/journal.pcbi.1002195

PubMed Abstract | CrossRef Full Text | Google Scholar

Galindo, L. J., Prokina, K., Torruella, G., López-García, P., Moreira, D. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Gautheret, D., Lambert, A. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Giegé, R., Sissler, M., Florentz, C. (1998). Universal rules and idiosyncratic features in tRNA identity. Nucleic Acids Res. 26 (22), 5017–5035.

PubMed Abstract | Google Scholar

Hafez, M., Burger, G., Steinberg, S. V., Lang, B. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Hafez, M., Majer, A., Sethuraman, J., Rudski, S. M., Michel, F., Hausner, G. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanson, D. K., Lamb, M. R., Mahler, H. R., Perlman, P. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Heiss, A. A., Warring, S. D., Lukacs, K., Favate, J., Yang, A., Gyaltshen, Y., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Henikoff, S., Henikoff, J. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Janouškovec, J., Tikhonenkov, D. V., Burki, F., Howe, A. T., Rohwer, F. L., Mylnikov, A. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Jung, J., Kim, J. I., Jeong, Y. S., Yi, G. (2018). AGORA: organellar genome annotation from the amino acid and nucleotide references. Bioinformatics 34 (15), 2661–2663. doi: 10.1093/bioinformatics/bty196

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D., Langmead, B., Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12 (4), 357–360. doi: 10.1038/nmeth.3317

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D., Paggi, J. M., Park, C., Bennett, C., Salzberg, S. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Lambert, A., Fontaine, J. F., Legendre, M., Leclerc, F., Permal, E., Major, F., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, B. F., Burger, G. (2012). ''Mitochondrial and eukaryotic origins: a critical review' in Advances in Botanical Research 63, 1-20 doi: 10.1016/B978-0-12-394279-1.00001-6

CrossRef Full Text | Google Scholar

Lang, B. F., Burger, G., O'Kelly, C. J., Cedergren, R., Golding, G. B., Lemieux, C., et al. (1997). An ancestral mitochondrial DNA resembling a eubacterial genome in miniature. Nature 387 (6632), 493–497. doi: 10.1038/387493a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, B. F., Cedergren, R., Gray, M. 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.

PubMed Abstract | Google Scholar

Lang, B. F., Laforest, M. J., Burger, G. (2007). Mitochondrial introns: a critical view. Trends Genet. 23, 119–125. doi: 10.1016/j.tig.2007.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, B. F., Lavrov, D., Beck, N., Steinberg, V. (2011). “Mitochondrial tRNA structure, identity and evolution of the genetic code”. in Organelle genetics. Ed. Bullerwell, C. E. (Berlin, Heidelberg, 2012: Springer), 431–474.

Google Scholar

Lavrov, D., Lang, B. F. (2013). ''Mitochondrial genomes in unicellular relatives of animals (ASBMB)''. in Molecular Life Sciences: An Encyclopedic Reference. Eds. Wells, R. D., Bond, J. S., Klinman, J., Masters, B. S. S., Bell, E. (New York : Springer) doi: 10.1007/978-1-4614-6436-5_178-2.

CrossRef Full Text | Google Scholar

Lavrov, D. V., Pett, W. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics 25 (16), 2078–2079. doi: 10.1093/bioinformatics/btp352

PubMed Abstract | CrossRef Full Text | Google Scholar

Ling, J., Daoud, R., Lajoie, M. J., Church, G. M., Soll, D., Lang, B. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Muñoz-Gómez, S. A., Mejía-Franco, F. G., Durnin, K., Colp, M., Grisdale, C. J., Archibald, J. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Musier-Forsyth, K., Usman, N., Scaringe, S., Doudna, J., Green, R., Schimmel, P. (1991). Specificity for aminoacylation of an RNA helix: an unpaired, exocyclic amino group in the minor groove. Science 253 (5021), 784–786. doi: 10.1126/science.1876835

PubMed Abstract | CrossRef Full Text | Google Scholar

Nadimi, M., Beaudet, D., Forget, L., Hijri, M., Lang, B. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Nawrocki, E. P., Eddy, S. R. (2013). Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics 29 (22), 2933–2935. doi: 10.1093/bioinformatics/btt509

PubMed Abstract | CrossRef Full Text | Google Scholar

Noutahi, E., Calderon, V., Blanchette, M., Lang, F. B., El-Mabrouk, N. (2017). CoreTracker: accurate codon reassignment prediction, applied to mitochondrial genomes. Bioinformatics 33 (21), 3331–3339. doi: 10.1093/bioinformatics/btx421

PubMed Abstract | CrossRef Full Text | Google Scholar

Okimoto, R., Macfarlane, J. L., Wolstenholme, D. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Pereira de Souza, A., Jubier, M. F., Delcher, E., Lancelin, D., Lejeune, B. (1991). A trans-splicing model for the expression of the tripartite nad5 gene in wheat and maize mitochondria. Plant Cell 3 (12), 1363–1378. doi: 10.2307/3869315

PubMed Abstract | CrossRef Full Text | Google Scholar

Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T. C., Mendell, J. T., Salzberg, S. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33 (3), 290–295. doi: 10.1038/nbt.3122

PubMed Abstract | CrossRef Full Text | Google Scholar

Prince, S., Munoz, C., Filion-Bienvenue, F., Rioux, P., Sarrasin, M., Lang, B. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Pruitt, K. D., Tatusova, T., Maglott, D. R. (2003). NCBI reference sequence project: update and current status. Nucleic Acids Res. 31 (1), 34–37. doi: 10.1093/nar/gkg111

PubMed Abstract | CrossRef Full Text | Google Scholar

Schafer, B., Hansen, M., Lang, B. F. (2005). Transcription and RNA-processing in fission yeast mitochondria. RNA 11 (5), 785–795. doi: 10.1261/rna.7252205

PubMed Abstract | CrossRef Full Text | Google Scholar

Seif, E. R., Forget, L., Martin, N. C., Lang, B. F. (2003). Mitochondrial RNase p RNAs in ascomycete fungi: lineage-specific variations in RNA secondary structure. RNA 9 (9), 1073–1083. doi: 10.1261/rna.5880403

PubMed Abstract | CrossRef Full Text | Google Scholar

Seif, E., Leigh, J., Liu, Y., Roewer, I., Forget, L., Lang, B. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Shang, J., Yang, Y., Wu, L., Zou, M., Huang, Y. (2018). The S. pombe mitochondrial transcriptome. RNA 24 (9), 1241–1254. doi: 10.1261/rna.064477.117

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheffield, N. C., Hiatt, K. D., Valentine, M. C., Song, H., Whiting, M. F. (2010). Mitochondrial genomics in Orthoptera using MOSAS. Mitochondrial DNA 21 (3-4), 87–104. doi: 10.3109/19401736.2010.500812

PubMed Abstract | CrossRef Full Text | Google Scholar

Sherif, T., Rioux, P., Rousseau, M. E., Kassis, N., Beck, N., Adalat, R., et al. (2014). CBRAIN: a web-based, distributed computing platform for collaborative neuroimaging research. Front. Neuroinform 8. doi: 10.3389/fninf.2014.00054

CrossRef Full Text | Google Scholar

Shiratori, T., Ishida, K. I. (2016). A new heterotrophic cryptomonad: Hemiarma marina n. g., n. sp. J. Eukaryot Microbiol. 63 (6), 804–812. doi: 10.1111/jeu.12327

PubMed Abstract | CrossRef Full Text | Google Scholar

Slater, G. S. C., Birney, E. (2005). Automated generation of heuristics for biological sequence comparison. BMC Bioinf. 6, 31. doi: 10.1186/1471-2105-6-31

CrossRef Full Text | Google Scholar

Strassert, J. F., Tikhonenkov, D. V., Pombert, J. F., Kolisko, M., Tai, V., Mylnikov, A. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, D., Lieberman, A., Lang, B. F., Simonovic, M., Soll, D., Ling, J. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Tekle, Y. I., Wang, F., Wood, F. C., Anderson, O. R., Smirnov, A. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Thorvaldsdóttir, H., Robinson, J. T., Mesirov, J. P. (2013). Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 14 (2), 178–192. doi: 10.1093/bib/bbs017

PubMed Abstract | CrossRef Full Text | Google Scholar

Tillich, M., Lehwark, P., Pellizzer, T., Ulbricht-Jones, E. S., Fischer, A., Bock, R., et al. (2017). GeSeq - versatile and accurate annotation of organelle genomes. Nucleic Acids Res. 45 (W1), W6–w11. doi: 10.1093/nar/gkx391

PubMed Abstract | CrossRef Full Text | Google Scholar

Turk, E. M., Das, V., Seibert, R. D., Andrulis, E. D. (2013). The mitochondrial RNA landscape of Saccharomyces cerevisiae. PloS One 8 (10), e78105. doi: 10.1371/journal.pone.0078105

PubMed Abstract | CrossRef Full Text | Google Scholar

Valach, M., Burger, G., Gray, M. W., Lang, B. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Wissinger, B., Schuster, W., Brennicke, A. (1991). Trans splicing in Oenothera mitochondria: nad1 mRNAs are edited in exon and trans-splicing group II intron sequences. Cell 65 (3), 473–482. doi: 10.1016/0092-8674(91)90465-B

PubMed Abstract | CrossRef Full Text | Google Scholar

Wyman, S. K., Jansen, R. K., Boore, J. L. (2004). Automatic annotation of organellar genomes with DOGMA. Bioinformatics 20 (17), 3252–3255. doi: 10.1093/bioinformatics/bth352

PubMed Abstract | CrossRef Full Text | Google Scholar

Yabuki, A., Gyaltshen, Y., Heiss, A. A., Fujikura, K., Kim, E. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Yazaki, E., Yabuki, A., Imaizumi, A., Kume, K., Hashimoto, T., Inagaki, Y. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zubaer, A., Wai, A., Hausner, G. (2018). The mitochondrial genome of Endoconidiophora resinifera is intron rich. Sci. Rep. 8 (1), 17591. doi: 10.1038/s41598-018-35926-y

PubMed Abstract | CrossRef Full Text | Google Scholar