Categorizing 161 plant mitochondrial group II introns into 29 families of related paralogues finds only limited links between intron mobility and intron-borne maturases

doi:10.21203/rs.3.rs-2086450/v1

Download PDF

Research Article

Categorizing 161 plant mitochondrial group II introns into 29 families of related paralogues finds only limited links between intron mobility and intron-borne maturases

https://doi.org/10.21203/rs.3.rs-2086450/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 12 Mar, 2023

Read the published version in BMC Ecology and Evolution →

You are reading this latest preprint version

Group II introns are common in the two endosymbiotic organelle genomes of plants. Chloroplasts harbor 22 positionally conserved group II introns whereas their occurrence in land plant mitogenomes is highly variable and specific for the seven major embryophyte clades: liverworts, mosses, hornworts, lycophytes, ferns, gymnosperms and flowering plants. Each plant group features “signature selections” of ca. 20–30 paralogues from a superset of altogether 105 group II introns meantime identified in embryophyte mtDNAs, suggesting massive intron gains and losses along the backbone of plant phylogeny. We report on systematically categorizing plant mitochondrial group II introns into “families”, comprising evidently related paralogues at different insertion sites, which may even be more similar than their respective orthologues in phylogenetically distant taxa. Including streptophyte algae extends our sampling to 161 and we sort 104 streptophyte mitochondrial group II introns into 25 core families of related paralogues evidently arising from retrotransposition events. Adding to discoveries of only recently created intron paralogues, hypermobile introns and twintrons, our survey led to further discoveries including previously overlooked “fossil” introns in spacer regions or e.g., in the rps8 pseudogene of lycophytes. Initially excluding intron-borne maturase sequences for family categorization, we added an independent analysis of maturase phylogenies and find a surprising incongruence between intron mobility and the presence of intron-borne maturases. Intriguingly, however, we find that several examples of nuclear splicing factors meantime characterized simultaneously facilitate splicing of independent paralogues now placed into the same intron families. Altogether this suggests that plant group II intron mobility, in contrast to their bacterial counterparts, is not intimately linked to intron-encoded maturases.

Group II introns

Retrotransposition

Maturases

Group I introns

mtDNA

Mitogenomes

Horizontal Gene Transfer

At first glance, the two endosymbiont genomes in the plant cell seem to have followed very similar evolutionary trajectories with massive gene transfers into the nucleus, strongly reducing the ancestral genomes of an a-proteobacterium and a cyanobacterium that gave rise to mitochondria and chloroplasts, respectively. Even later emerging molecular mechanisms such as the site-specific C-to-U RNA editing characteristic for land plants affect the transcriptomes of both plant organelles equally [1]. In contrast, another characteristic feature of the two endosymbiont genomes in plants shows a striking discrepancy: the overwhelming stasis in the occurrence of group II introns in chloroplast DNAs for more than 500 million years, even extending into the streptophyte algal ancestors vs. the strikingly dynamic evolution of their counterparts in plant mitochondrial genomes. Examples for the latter are documented with the multiple splicing pathways of plant mitochondrial group II introns [2–4], their many transitions from cis- into trans-splicing arrangements [e.g. 5–12], the recently discovered functional or degenerated twintrons [13] or the evident degeneration of intron-borne maturase reading frames [5, 14] and the creation of nuclear-encoded maturase genes [15–20] in the multifarious pathways of evolution (Fig. 1A-E). Here, we focus on the occupation of new loci by mobile plant mitochondrial group II introns, including examples of introns “fossilized” in pseudogenes or intergenic regions (Fig. 1F-G).

The peculiar clade-specific patterns in the occurrence of plant mitochondrial group II introns [21–23] and the examples for evolutionary recent occupations of new intron insertion sites [24, 25] suggest that plant mitogenomes could offer a particularly rich and attractive data source to elucidate the mobility of group II introns on a broad evolutionary scale. Previous research on group II intron mobility has so far largely focused on functional studies with selected group II model introns in fungi or bacteria that have helped to elucidate the biochemical pathways of intron propagation and invasion of new loci [26–28]. Moreover, many features make group II introns particularly interesting ribozyme RNAs in general [29–38]: Firstly, the biochemical mechanism of their splicing, involving two transesterifications using a 5’-terminal guanosine and a looped out adenosine upstream of the 3’-splice acceptor site, is homologous to that of spliceosomal introns in the nuclear genomes of eukaryotes. In fact, all available evidence convincingly suggests group II introns to be the evolutionary ancestors of the eukaryotic spliceosome machinery [39–42]. Secondly, some group II introns have been demonstrated to have self-splicing activity in vitro, a notable case being the second intron in the mitochondrial rrnL gene of the brown alga Pylaiella littoralis [43] that splices efficiently even at low Mg²⁺ ion concentrations. Thirdly, in their likely evolutionary ancestral state, group II introns carry reading frames for maturases, intron-encoded proteins (IEPs) that are crucial to assist both in splicing but also in the mobility, allowing their host introns to occupy new insertion sites. In their most complete forms maturases consist of DNA-binding and endonucleolytic cleavage (D/En) and reverse transcriptase (RT) domains allowing the conversion of the RNA spliced into a novel DNA location, thereby creating a new intron paralogue.

However, none of the particularly abundant group II introns in plant organelle genomes have been experimentally demonstrated to be mobile and reports on self-splicing or in vitro splicing in the presence of chloroplast extracts, respectively, have only recently been published [44, 45]. The lack of plant mitochondrial group II intron self-splicing activity is likely associated with the diverse set of nuclear-encoded proteins that act as either intron-specific or more promiscuous splicing factors [46, 47] and which continue to be characterized in ever-increasing numbers by reverse genetic studies, as it is easily documented with only some most recently published examples [48–52].

While studies of group II intron mobility in plant mitogenomes are hampered by lacking experimental approaches to manipulate plant mtDNAs, the extraordinary diversity of introns in plant mitogenomes offers an intriguing evolutionary perspective. Sequence similarities between plant mitochondrial intron paralogues had become evident soon after the first complete plant mitogenome sequences became available [53]. For example, the overall sequence similarity between the first intron in the nad2 gene (nad2i156g2) and the second intron of the nad1 gene (nad1i477g2) in flowering plant mitogenomes has been recognized early [54]. Similarly, the case of nad5i392g2 is striking, an intron that is phylogenetically very restricted to the family Lycopodiaceae among the lycophytes but clearly related to downstream intron nad5i1242g2 in the nad5 gene, which has a much wider phylogenetic conservation among lycophytes and ferns [55]. More recently emerging examples for evidently related group II intron paralogue pairs are two cases in ferns: atp1i361g2 and rps3i249g2 [25] and rps1i25g2 and rpl2i846g2 [24]. The hitherto most striking example of extraordinary sequence similarity between two intron paralogues has only very recently been discovered in hornwort mtDNAs accompanying the discovery of twintrons [13]. Intron paralogue cox2i98g2 is exclusively present in the hornwort genus Anthoceros and shares more than 98% sequence similarity with the internal group II intron of a twintron in the atp1 gene (atp1i1050g2ii1536g2), which is widely distributed among the hornworts. This example immediately suggests a recent copying event from atp1 into cox2 late in the recent phylogeny of hornworts.

Here, we systematically scanned available mitogenomes of streptophytes (i.e., land plants and the phylogenetically most closely related “charophyte” green algae) and scored the relationships between their group II introns as an effort to elucidate how intron migrations and losses may have contributed to shape the mitogenome makeups in extant plants. We define 25 core “families” F01 - F25 of clearly sequence-related group II intron paralogues, which indicate ancient and recent intron copying events. With the extraordinary slow sequence evolution in mitogenomes of the plant lineage combined with an increasingly well-understood phylogeny of land plants and related streptophyte algae our approach will offer a new perspective on an evolutionary timescale covering more than 500 million years. Altogether, we were able to assign 104 from a total of 161 streptophyte mitochondrial group II introns to one of the now defined 25 families whereas the 63 remaining introns presently lack significant similarity to any other paralogues. We discuss possible intron migration scenarios, considering the likely important role of intron-encoded or free-standing maturases.

The streptophyte mitochondrial group II intron data set

We scanned plant mitochondrial genomes for the occurrence of group II introns, including the complete phylogenetic diversity representing the seven major land plant clades: Flowering plants (Angiosperms), Gymnosperms, Ferns (Monilophytes), Lycophytes, Hornworts (Anthocerotophyta), Mosses (Bryophyta) and Liverworts (Marchantiophyta). We additionally included the available mitogenomes of streptophyte algae (“Charophytes”) representing six classes: Zygnematophyceae (Desmidiales: Closterium bailyanum, Gonatozygon brebissonii; Zygnematales: Entransia fimbriata, Roya anglica, R. obtusa and Zygnema circumcarinatum) recently considered to be most closely related to the plant lineage as well as Charophyceae (Chara vulgaris, C. braunii, Nitella hyalina), Coleochaetophyceae (Chaetosphaeridium globosum, Coleochaete scutata), Chlorokybophyceae (Chlorokybus atmophyticus), Klebsormidiophyceae (Klebsormidium flaccidum) and the early-branching Mesostigmatophyceae (Mesostigma viride). For a clear designation of introns we use the previously suggested intron nomenclature based on the intron insertion site in a given gene behind the reference position in the respective homologue of the liverwort Marchantia polymorpha [56, 57].

The example of the cox2 gene (Fig. 2) is used to introduce into important issues of our analyses addressing the huge diversity of mitochondrial introns in the plant lineage. Altogether twelve different group II insertion sites are presently identified in the cox2 genes of land plants and streptophyte algae. Their phylogenetic distributions vary widely from introns present in most land plant lineages excluding liverworts (cox2i373g2 and cox2i691g2) to others presently only identified in the hornwort genus Anthoceros (cox2i98g2) or in the streptophyte alga Coleochaete scutata (cox2i550g2). The assignments of introns to core families, which comprise significantly similar intron paralogues as introduced in this work, are given below the intron names: cox2i381g2 in family F01, cox2i550g2 in F03, cox2i97g2 and cox2i98g2 in F04, cox2i564g2 in F16 and cox2i94g2 in F19 (see below). Including intron-based maturases allows for definition of maturase-based families and superfamilies, here cox2i373g2, in mF26 and cox2i127g2 in mF27. Remaining “solitary” introns lacking significant similarity to other paralogues are labeled “S”. Insertion sites must be considered very carefully and precisely, notably when introns occur differently in a small gene region such as introns cox2i94g2, cox2i97g2, cox2i98g2 and cox2i104g2, an issue occasionally overlooked and leading to mis-annotation in database entries. Introns in the same position were considered as orthologues also when lacking significant sequence similarities across large phylogenetic distances if they were not assigned to different families.

Our re-evaluation of database entries during intron sampling has made us re-consider several intron and splice site annotations and allowed us to suggest yet further, and very likely functionally splicing, introns that were previously unnoticed but would reconstitute important and conserved parts of their host genes (e.g. cox3i34g2 in Gonatozygon, nad2i81g2 in Zygnema or nad9i89g2 in Closterium and Gonatozygon) as well as dysfunctional, “fossil” introns (e.g. rps8i52g2f in Phlegmariurus), as we will discuss below. Our final total mitochondrial group II intron sampling in streptophyte mitogenomes comprised 161 group II intron paralogues defined by their unique insertion sites.

The occurrence of streptophyte mitochondrial group II intron paralogues in the major clades is displayed with Euler diagrams in Fig. 3. The striking discrepancy of mitochondrial vs. chloroplast group II introns is immediately apparent with 22 land plant chloroplast introns all of which have counterparts in streptophyte algae vs. altogether 161 mitochondrial group II introns, of which only 13 are shared between embryophytes and streptophyte algae (Fig. 3A). Differentiating among the embryophytes, larger intersections are found between bryophytes and tracheophytes (Fig. 3B) than with either group and the outgroup algae and among the latter between hornworts and mosses (Fig. 3C) and between hornworts, mosses and tracheophytes (Fig. 3D), respectively. Notably, of 101 group II introns identified in embryophyte mtDNAs, only one (atp9i87g2) is shared between all three bryophyte clades (Fig. 3C).

The group II intron family concept

Despite their highly conserved six-domain structures, group II intron paralogues a priori share no significant sequence similarities aside from mostly conserved sequence motifs at the 5’-end (mostly GUGCG), the 3’- termini (mostly AY) and their characteristic domain V comprising 34 nucleotides, which mostly folds into two base-paired regions of 9 and 5 base pairs, respectively, with a dinucleotide bulge and a terminal tetraloop (most often GNRA) in the overwhelming majority of group II introns. Only a consensus sequence representing domain V may in fact be used as an initial query to scan for group II intron candidates in sequence databases [58].

Our criteria for considering introns occupying different insertion sites for inclusion into families of related paralogues are detailed under methods. In brief, we used several rounds of identifying sequence-related intron paralogues that share significant sequence similarities (beyond short similarities of domain V and the immediate flanking regions), which can exceed those of evidently orthologous introns in the same positions in distant plant taxa. Naturally, some rare borderline cases are represented with introns occupying the same insertion site in phylogenetically very distant taxa. An independent insertion of a given position cannot be excluded as an alternative explanation to vertical transmission followed by sequence divergences obliterating recognizable similarities. Below, we will discuss such borderline cases in the context of our consideration of the now defined group II intron families. A cladogram of group II introns sorted into families is shown in Fig. 4. Subsequently, we will consider, where present, similarities of intron-borne maturases to independently verify family assignments based on the nucleotide sequence similarities alone. Furthermore, maturase similarities will occasionally allow for the inclusion of some previously “solitary” introns lacking significant nucleotide sequence similarities to other paralogues into extended families and for some fusions of primary intron families into “superfamilies” (Fig. 5).

Family F01

Among the here defined families of mitochondrial group II introns in the streptophyte lineage, family F01 is the largest one, comprising altogether 13 intron paralogues (Fig. 4, Fig. 6). It includes four isolated introns, each presently identified in only one streptophyte algae clade (atp9i209g2 and cox3i745g2 in Chlorokybus, atp9i214g2 in Charophyceae and cox1i1039g2 in Coleochaete, respectively). The “moss-specific” introns cox3i507g2 and nad9i283g2 of F01 are universally conserved in that plant clade. Phylogenetically more broadly distributed introns atp6i439g2, nad4i1399g2, nad5i230g2, nad7i209g2 and rps3i74g2 occur in up to five embryophyte clades. In two cases, the respective introns remain recognizable despite pseudogene degeneration of nad7 and rps3 in hornwort mitogenomes. Three introns in family F01 are present both in at least one embryophyte clade and in up to three classes of streptophyte algae: rps3i74g2, atp6i439g2 and cox2i381g2 (Fig. 4, Fig. 6). Notably though, family F01 does not include an intron paralogue present in liverworts. Only four of the intron paralogues in F01 carry maturases: atp9i209g2, cox1i1039g2, cox2i381g2 and nad3i211g2 (Fig. 4, Fig. 6).

Although outgroup rooting is necessarily difficult for the individual clades of intron families, the unrooted phylogeny suggests particular close relations of certain introns in F01, in this case of the “moss” intron cox3i507g2 with the “tracheophyte” intron nad4i1399g2 and of the phylogenetically wider distributed paralogues nad5i230g2 with rps3i74g2 (Fig. 6). Interestingly, the latter two and nad7i209g2 as a third intron paralogue of the F01 family have recently been demonstrated to be affected by the same nuclear-encoded splicing factor in the model angiosperm Arabidopsis thaliana, OZ2 [50]. Also quite notably, the maize NCS2 mutant (non-chromosomal stripe mutant) described long ago [59] is the result of a mitogenome recombination between two mitochondrial introns that we have now assigned to family F01: nad4i1399g2 and nad7i209g2.

The phylogeny of intron paralogues is in accordance with known species phylogenies for the well-supported nodes, while apparent phylogenetic “mis-placements” e.g., of the diverging lycophyte sequences, remain without support (Fig. 6). However, and notably so with respect to phenomena like shared splicing factors, the possibility of concerted and converging evolution must be considered. Whereas the case of atp6i439g2 is an example where even the distant algal orthologue in Zygnema is included with its counterparts in lycophytes and hornworts in one subclade, this is not the case for rps3i74g2 where the Charales orthologues are attracted to the isolated atp9i214g2 paralogue in Chara vulgaris. Our later additional considerations of maturase protein similarities that we will discuss below includes family F01 into a “superfamily” also including the small family F22 and two “solitary” introns (Fig. 5).

Family F02

As in family F01, some intron paralogues in family F02 likewise have a very restricted occurrence: cox3i625g2 is exclusively present in liverworts and nad9i246g2 is only identified in hornworts (Fig. 4). However, despite their phylogenetically disjunct occurrence, cox3i625g2 and nad9i246g2 are closely related paralogues in F02 (Fig. 7). Yet more phylogenetically restricted, nad5i392g2 is only present in the lycophyte order Lycopodiales, here represented by Phlegmariurus squarrosus. The extreme diversity of mitogenome evolution within the lycophyte clade with retained genes and introns and a conserved genetic synteny in the P. squarrosus mtDNA versus highly recombining mitogenomes with reduced gene sets in Isoetales and Selaginellales is also fully in line with intron rps10i235g2 being retained in Phlegmariurus and shared with the seed plants. The similarity of angiosperm rps10i235g2 to liverwort introns rrnLi833g2 and cox3i625g2 had already been noted right along with its initial discovery [60]. Other than rps10i235g2, P. squarrosus intron nad5i1242g2 is shared with ferns, similarly indicating an early vascular plant ancestry. Notably, we recently found intron rrnLi833g2 universally conserved among liverworts now to be also present in the mitogenome of the leptosporangiate fern Haplopteris ensiformis [61] and in a preliminary mtDNA assembly of the fern Dryopteris crassirhizoma (database accession MW732172). Intriguingly, however, the fern rrnLi833g2 introns cluster with their nad5i1242g2 paralogues, possibly indicating concerted evolution or a loss-and-regain scenario.

Intron family F02 neither includes an extant intron paralogue present in mosses nor one carrying a maturase ORF (Fig. 4). Notably though, we were able to identify several additional fossilized group II intron sequences “F02g2f” clearly tracing back to family F02 intron paralogues both in Phlegmariurus squarrosus and in hornwort mitogenomes (Suppl. Figure 1A & B). The series of intron copying events leading to the five F02 paralogues now recognized in the Phlegmariurus mitogenome (Fig. 7B) remains unclear except for nad5i392g2 present only in Lycopodiales and likely emerging late from the more ancestral nad5i1242g2 paralogue. The newly identified fossil intron rps8i51g2f receives the ‘f’ behind the g2 intron label because we could not confirm functional splicing despite significant similarity to its functional counterparts in nad5 and rps10 extending up to the very intron 5’- and 3’-ends as opposed to another intergenic F02 intron fossil in the tatC-cox2 spacer (Fig. 7C). Altogether five intergenic F02 intron fossils could be discovered in the hornwort mitogenomes, likely originating from nad9i246g2 or a common ancestor (Fig. 7D). In this case, the two intergenic intron fossils between rps4 and tatC and between atp6 and nad6 appear more closely related to each other than the remaining three. The extension of the intron sequence similarities ending with the respective 5’- or 3’-intron termini and their identification in the conservatively evolving mitogenomes of early land plants devoid of heavy recombination as in seed plant mitogenomes strongly argues for their origin by retrotransposition rather than mtDNA recombination.

Family F03

Intron family F03 contains two intron paralogues, nad4i976g2 and nad5i1477g2, shared between hornworts and all four tracheophyte clades, potentially supporting a hornwort-tracheophyte (HT) clade (Fig. 4). However, both nad4i976g2 and another intron paralogue in F03, nad3i140g2, also reveal counterparts in Charophyceae algae. No losses are identified for nad3i140g2 among available liverwort, hornwort or lycophyte mitogenomes. Similarly, two further intron paralogues are universally conserved in mosses and shared with either lycophytes (atp9i21g2) or hornworts (sdh3i100g2). Finally, three F03 introns have a restricted occurrence among streptophyte algae, of which only one is shared with a land plant clade, cox2i250g2 present in Closterium (Zygnematophyceae) and in the liverworts. This intron is an interesting exceptional case with a maturase carried in the liverwort orthologues, but not in the algal counterpart of cox2i250g2. Among the eight intron paralogues in F03 (Fig. 4), two pairs of introns show particularly high sequence similarities: nad3i140g2 with nad5i1477g2 and sdh3i100g2 with nad4i976g2. In this case, the intron paralogues occur with a phylogenetic overlap in the hornworts, a clade characterized by high intron mobility as also reflected below with the example of family F04. Again, we were able to identify fossilized group II intron sequences (“F03g2f”) in intergenic sequences of the lycophyte Phlegmariurus squarrosus clearly tracing back to family F03 (Suppl. Figure 1C).

Family F04

Group II intron family F04 is dominated by introns exclusively occurring in hornworts, indicating a pronounced intron mobility in that land plant clade: cox2i98g2, nad5i881g2, nad6i444g2 and nad9i502g2 (Fig. 4). This is further enhanced by the observation that the respective internal introns of two recently discovered twintrons [13], atp1i1050g2ii1536g2 and cox1i1116g1ii207g2 in the hornwort genus Anthoceros, likewise belong to the here defined F04 family. In fact, the extraordinary sequence similarity of 98% between atp1i1050g2ii1536g2 to cox2i98g2 is the most striking case of closely related paralogues in our entire intron sampling, indicating a very recent copying event. Considering that the atp1 twintron is present in all hornworts whereas cox2i98g2 exists only in the genus Anthoceros immediately suggests a copying from the former to the latter insertion site. Notably, however, internal splicing of zombie-twintron atp1i1050g2ii1536g2 could not be detected.

At the same time, family F04 is a prime example for the necessity of careful analyses to distinguish cox2i98g2 from its paralogue cox2i97g2 in liverworts that is inserted one nucleotide upstream in the cox2 gene (Fig. 2). Intron nad7i676g2 is conserved in all tracheophyte clades. Its former presence in hornworts, and accordingly in a possible HT stem lineage, is elusive owing to the pseudogene degeneration of the downstream part of nad7 in all hornwort mtDNAs. Notably, and despite the very clear evidence of recently active retrotransposition, none of the F04 members carries a maturase-ORF (Fig. 4).

Family F05

In contrast to F04 mainly comprising intron paralogues in hornworts, group II intron family F05 (Fig. 4) is dominated by intron paralogues exclusively occurring in liverworts: nad4i548g2, nad4Li100g2, nad7i336g2 and rpl2i28g2. Liverwort-specific intron nad4i548g2 had been introduced to elucidate the liverwort phylogeny [62] and is now found to be notably similar to rpl2i28g2. Other than the four liverwort paralogues, F05 also includes lycophyte introns recently discovered and characterized as the external intron of a “zombie” twintron upon closer reinspection of the Phlegmariurus squarrosus sdh3 gene [13].

Family F06

Apart from rps3i249g2 that is shared between lycophytes, ferns and gymnosperms, the group II intron paralogues in family F06 are phylogenetically very restricted in occurrence: atp1i361g2 exclusively in ferns, cobi783g2 only in liverworts and cox1i150g2 exclusively in hornworts (Fig. 4). A detailed study of atp1i361g2 concluded that this fern-specific intron has originated from the, likely more ancestral, paralogue rps3i249g2 and also found evidence for convergent evolution of specific intron structures, mainly group II intron domain III, in the two paralogues in later emerging fern lineages [25].

Interestingly, a truncated copy of cobi783g2 had been identified earlier in the spacer between nad5 and nad4 in liverworts in an early sampling of this intergenic region among bryophytes [63]. Our extended consideration of intron families now adds support to the idea that also this case of an intron fossil has arisen through copying, likely by a retrotransposition event. Similarity of the fossil sequence starts exactly from the intron 5’-end and extends for nearly 800 bp with 99% identity in the case of Treubia lacunosa representing an early liverwort branch (Suppl. Figure 1D) whereas a higher degree of degeneration is observed in derived taxa. Yet other intron fossils have now been discovered in the intergenic region between rrn5 and trnM-CAU of the hornwort Anthoceros agrestis and the moss Sphagnum palustre mitogenomes (Suppl. Figure 1E). The only evidence for a maturase among the four members of intron family F06 are traces of a former maturase-ORF in cox1i150g2.

Family F07

Intron family F07 contains three paralogues of phylogenetically wide distribution and only one member with isolated occurrence in the liverworts alone: cobi372g2 (Fig. 4). The similarity between introns nad1i477g2 and nad2i156g2, now found to be widely conserved among vascular plants (Fig. 4), has been recognized early after the complex structures of nad1 and nad2 in flowering plants had been elucidated [54]. Quite interestingly, the nuclear encoded splicing factor ODB1 has meantime been found to promote splicing of both of these two paralogues [64]. The fourth paralogue in family F07, intron atp9i87g2, is particularly noteworthy for (i) being an intron now identified to be shared between all three bryophyte clades, (ii) carrying an ancestral maturase that has independently degenerated in all four plant clades where it is present and (iii) carrying the internal intron atp9i87g2ii1114g2 of family F08 (see below) turning it into a twintron in the lycophyte Phlegmariurus squarrosus. A history of intron paralogues is immediately suggestive for the family F07 paralogues (Fig. 8). The ancestral maturase-carrying atp9i87g2 independently likely gave rise to cobi372g2 in liverworts and to nad2i156g2 in a possible non-liverwort lineage. After loss of nad2i156g2 in hornworts, it gave rise to nad1i477g2 only in the tracheophyte lineage.

Family F08

Two introns in family F08 have a phylogenetically striking distribution being conserved in liverworts but also occurring in lycophytes: nad7i1113g2 and rps14i114g2 (Fig. 4). Given that the host genes, nad7 and rps14, are frequently subject to EGT, nad7i1113g2 is at present only determined in Isoetes engelmannii and rps14i114g2 only in Phlegmariurus squarrosus among the lycophytes. However, we now found rps14i114g2 equally conserved in the mitogenome of the leptosporangiate fern Haplopteris ensiformis [61]. Yet more notably, the recently characterized inner intron of a twintron in the atp9 gene, atp9i87g2ii1114g2 [13] as a third member in family F08 is characteristically more similar to its rps14i114g2 counterpart in liverworts than in P. squarrosus.

Family F09

Intron ccmFCi829g2 in F09 has a phylogenetically wide distribution in bryophytes and tracheophytes (Fig. 4). Its former presence in lycophytes remains unclear, however, owing to the loss of the entire ccm gene suite for cytochrome c maturation in this clade. In contrast, pseudogene traces of ccmFC including ccmFCi829g2 are clearly detected in hornworts [65]. Intron rpl2i846g2 is evidently a gain in the tracheophyte stem lineage. This ancestral intron obviously gave rise to its closely related paralogue rps1i25g2 exclusively present in ferns [24]. As in family F08, no traces of maturase-ORFs are recognizable in any of the F09 intron paralogues. Highly interesting, however, the splicing of both F09 paralogues present in Arabidopsis thaliana, ccmFCi829g2 and rpl2i846g2, was found to be affected by the same nuclear-encoded splicing factor, WTF9 [66].

Family F10

Intron family F10 comprises group II introns with a particularly complex history (Fig. 4). It contains two intron paralogues with a phylogenetically disjunct distribution: cobi824g2 exclusively present in liverworts and nad1i669g2 previously assumed to be restricted to tracheophytes. Intron nad1i669g2 has received attention as being conserved in a trans-splicing arrangement in seed plants. A cis-arranged orthologue was initially identified in the fern Osmunda regalis, also noticing traces of a degenerated maturase [5]. We now found that intron nad1i669g2 has a clear orthologue also in the mtDNA of the alga Coleochaete scutata [67] with a maturase-ORF annotated in the corresponding database entry (MN613583). Intron family F10 also contains the “hypermobile invader” intron cox1i1149g2 in the lycophyte Phlegmariurus squarrosus, which gave rise to internal introns of two twintrons (in the newly determined sdh3i349g2 and in itself) and to seven further intron fossils in intergenic regions [13] (Suppl. Figure 1F).

Family F11

Intron family F11 contains introns atp6i80g2 and atp9i95g2 of phylogenetically wider distribution, shared between mosses, hornworts and lycophytes (Fig. 4). Like the above cases of the neighboring introns cox2i97g2 and cox2i98g2 in F04, also atp9i95g2 requires careful inspection given the closely neighboring intron atp9i87g2 in family F07 [65]. Despite the evident dynamics of intron insertions in atp9 of bryophytes and lycophytes, no evidence is ever found for an intron in the atp9 gene among ferns, gymnosperms or angiosperms. Intron cox1i1064g2 is exclusively present in mosses and cox1i653g2 is only present in hornworts. Again, the latter needs particular attention to distinguish it from an unrelated “solitary” intron cox1i652g2 inserted one nucleotide upstream in the Coleochaete mitogenome.

Family F12

The two intron paralogues in family F12 (Fig. 4) have a strikingly divergent distribution with cobi274g2 presently only identified in Charales algae and nad2i1282g2 present in hornworts and euphyllophytes (i.e., ferns and seed plants). No traces of former maturase reading frames can be detected in nad2i1282g2 or cobi274g2.

Family F13

Group II intron family F13 contains two intron paralogues conserved in liverworts: cox1i178g2 and cox3i171g2 (Fig. 4). The former has counterparts in the algae Klebsormidium flaccidum and Coleochaete scutata. However, only cox1i178g2 in Coleochaete shares significant sequence similarity with the liverwort homologues, clearly warranting family inclusion according to our criteria. The “positional homologue” in Klebsormidium neither shares similarity with the liverwort nor with the Coleochaete counterpart, leaving its status as a true orthologue vs. a possible “analogue” occupying the same insertion site open at present. The cox1i178g2 introns carry maturase reading frames both in algae and the liverworts, no intron-borne ORFs are present in the much smaller cox3i171g2 paralogues shared between liverworts and the lycophyte Phlegmariurus squarrosus.

Family F14

As in family F12, the intron members in family F14 are also very disjunct in occurrence with nad1i258g2 being restricted to the monilophyte (fern) clade and nad1i517g2 so far only identified in the alga Zygnema circumcarinatum (Fig. 4). The third paralogue, cox1i511g2, however, is present in liverworts, mosses, lycopyhtes and the algae Coleochaete and Zygnema. Like most group II introns in the algal mitogenomes also cox1i511g2 carries maturase reading frames. The counterparts in the land plant lineage, however, are frameshifted ORFs in liverworts and degenerated or unrecognizable in the cox1i511g2 orthologues of mosses or in the lycophyte Selaginella moellendorffii. We identified an F14-type intron fossil (F14g2f) in the cox1-rrnS spacer of the Coleochaete mitogenome (Suppl. Figure 1G).

Family F15

As in F14, the intron paralogues in family F15 are likewise phylogenetically disjunct (Fig. 4). Intron cox1i835g2 is presently only identified in the algal genus Chara (and absent in the Characeae genera Nitella and Nitellopsis). Intron nad1i287g2 had initially been identified serendipitously in a screen for ancestors of trans-splicing intron nad1i394g2 [5] and is meantime found to be universally conserved both in mosses and in hornworts [56, 65]. As in the above cases, a maturase reading frame is only present in the algal cox1i835g2 but not in the much smaller nad1i287g2 paralogues in the bryophytes of less than 800 bp. The subsequent additional analyses of intron-borne maturases (see below) add S-type intron atp9i145g2 in Charales to an extended family eF15, which is ultimately linked to the superfamily SF13-14 (Fig. 5).

Family F16

Family F16 contains two intron paralogues that occur exclusively in hornworts: cox2i564g2 and nad1i348g2 (Fig. 4). The latter intron is of exceptionally small size of less than 600 nt. and conserved in all hornwort genera. In contrast, cox2i564g2 is ca. five times larger and lost together with the upstream and downstream neighboring introns cox2i381g2 and cox2i691g2 in Nothoceros aenigmaticus. Despite their extended sizes of more than 2.5 kb no maturase traces are discernible in the hornwort cox2i564g2 copies.

Family F17

Intron family F17 contains two intron paralogues with a phylogenetic distribution that could have been taken as further support for an NLE (“Non-Liverwort Embryophyte”) clade: nad1i728g2 and nad4i461g2 (Fig. 4). Both introns are absent in liverworts but particularly well conserved in mosses, hornworts, lycophytes and vascular plants with only rare exceptions including the absence of nad4i461g2 in the hornwort Leiosporoceros dussii and of nad1i728g2 in the lycophyte Isoetes engelmannii. Intron nad4i461g2 counterparts in the algae Coleochaete scutata and Zygnema circumcarinatum are large introns of 3.1 and 5.7 kb with long maturase reading frames that are continuous with the upstream nad4 coding sequence. Matching the general observations, only small traces of former maturases remain in the land plant counterparts where the sizes of nad4i461g copies in the mosses are reduced to less than 800 bp.

Intron nad1i728g2 is a particularly interesting case, being the only mitochondrial intron carrying a maturase reading frame in flowering plants, widely known as “matR”, now systematically labelled mat-nad1i728g2. Moreover, intron nad1i728g2 is also known for multiple independent transitions from cis- to trans-splicing with intron disruptions either up- or downstream of the maturase ORF in flowering plants [6]. Intriguingly, a gene transfer of the matR/mat-nad1i728g2 reading frame into the nuclear genome has been identified in Pelargonium [68]. Most interestingly, three different nuclear-encoded splicing factors have already been identified, which affect the two closely related F17 angiosperm paralogues nad1i728g2 and nad4i461g2 simultaneously: EMP8 [69], DEK55 [70] and SMK3 [71]. The third paralogue in F17, cobi399g2, is so far only identified in the Coleochaete scutata mitogenome and, as in most cases in the algal mitogenomes, also carries a maturase reading frame. The intergenic region between trnM-CAU and trnA-UGC contains an F17-type intron fossil (F17g2f) in the Phlegmariurus squarrosus mitogenome (Suppl. Figure 1H).

Family F18

Family F18 contains two intron paralogues presently only identified in algae of the Zygnematophyceae (Fig. 4): nad5i537g2 is present exclusively in Zygnema circumcarinatum and nad7i777g2 is present in Z. circumcarinatum and in Closterium baillyanum, of which only the latter homologue is equipped with an intron-borne maturase ORF.

Family F19

Group II intron family F19 (Fig. 4) comprises two introns present in lycophytes where sequence similarities are blurred by the highly divergent mitogenomes in genera Isoetes and Selaginella and largely rely on the conserved mitogenome of Phlegmariurus squarrosus. While intron cox2i94g2 is exclusively present in all three orders of lycophytes, intron paralogue cox1i323g2 is also present in the mitogenome of Sphagnum, representing a very early branch in the phylogeny of mosses. In both, Phlegmariurus and Sphagnum cox1i323g2 contains maturase remnants. Notably, cox1i323g2 has extensive similarity with the extended intergenic region between nad9 and trnI-CAU in hornwort mtDNAs, representing yet another example for traces of an intron fossil (Suppl. Figure 1I).

Family F20

Intron family F20 comprises two intron paralogues so far exclusively identified in the mitogenome of the alga Closterium baillyanum of the Zygnematophyceae: cox3i641g2 and nad9i275g2 (Fig. 4). Neither of the two intron paralogues carries a maturase reading frame.

Family F21

Intron family F21 contains introns presently only found in Coleochaete scutata: atp1i66g2, nad5i1725g2, rrnSi435g2 and trnH-GUGi32g2 (Fig. 4). Only the nad5i1725g2 paralogue carries a maturase. The size increase of the Coleochaete scutata mitogenome in comparison to that of Chaetosphaeridium globosum, also in the Coleochaetales, had in part been ascribed to the presence of 57 vs. only 11 introns [67]. We now find the here defined family F21 particularly interesting, because additional intron homologies are also present in multiple intergenic regions of the Coleochaete mitogenome: trnMf-nad9, mttB-trnL and trnV-trnD. In all three cases, homologies sharply coincide with the intron 5’-end and extend for minimally 600 bp indicating retrotransposition rather than events of DNA recombination at their origins (Suppl. Figure 1J).

Family F22

As in the case of F21, intron family F22 also includes members presently only identified in the mitogenome of Coleochaete scutata: atp1i850g2 and nad5i362g2 (Fig. 4). Of these two, only atp1i850g2 carries a maturase reading frame.

Family F23

Intron family F23 comprises two intron paralogues that are hitherto likewise only identified in streptophyte algae: nad7i925g2 in Closterium baillyanum and rrnLi2032g2 identified in Coleochaete scutata and Nitella hyalina (Fig. 4). Only the latter carries a maturase reading frame.

Family F24

Intron family F24 comprises rrnSi1148g2 present in Coleochaete scutata, rrnLi1747g2 presently only identified in Entransia fimbriata, nad7i250g2 present in the Zygnematophyceae algae Closterium baillyanum and Gonatozygon brebissonii, rrnLi629g2 present in Entransia, Coleochaete and Nitella and trnS-GCUi43g2, the only F24 paralogue shared with embryophytes (Fig. 4). The latter is present in Chlorokybus and the Zygnematophyceae genera Closterium, Gonatozygon and Roya and highly conserved among liverworts. Moreover, trnS-GCUi43g2 is evidently present as a degenerated copy in the trnS-GCU pseudogene retained in the conserved intergenic space between trnA and trnD in the mitogenomes of mosses and the lycophyte Phlegmariurus squarrosus. An independent degeneration has occurred among hornworts with functional copies present in Leiosporoceros and Anthoceros but pseudogenes in Nothoceros and Phaeoceros. Hence three independent degenerations of trnS-GCU and its intron are evident in land plants in mosses, among hornworts and in the tracheophytes (Suppl. Figure 1K).

Family F25

Family F25 comprises two intron paralogues of phylogenetically disjunct distribution: cobi537g2 present in the Charales algae and nad5i1455g2, present in all embryophytes except the liverworts (Fig. 4).

Maturases in the plant mitochondrial lineage and extended maturase-based intron families

Our categorization of 100 group II introns into the 25 “core” families F01 - F25 outlined above is based on their primary nucleotide sequence similarities alone. To cross-check for independent confirmation and to explore further and deeper relationships, also outside of the streptophytes, we independently compiled the intron-borne maturase ORFs present in 43 streptophyte mitochondrial group II introns as seeds for identifying protein homologs. This seed query data set also included several maturases that remained hitherto unnoticed or not annotated in database entries (e.g., the spliced variants of mat-atp9i87g2 and mat-atp9i95g2 in Phlegmariurus). The search for homologs ultimately resulted in a large protein data set that also contained related maturases of distant chlorophyte algae as well several maturase proteins in red algae, stramenopiles, fungi, animals and bacteria. The independent phylogenetic analysis of the large maturase protein sequence collection (Suppl. Figure 2) fully confirmed the identified core families F01, F03, F10, F11, F17 and F25, all of which contain at least two intron paralogues carrying maturases (Fig. 5). Moreover, the maturase similarities identified four additional, “maturase-based” group II intron families mF26-mF29 and helped to define superfamilies (SF) of higher order that combine the core intron families and include additional, previously solitary, introns.

For space limitations, we here focus on the examples of the large superfamily SF01-22 comprising families F01, F03, F22 and the two previously solitary introns cox1i652g2 and cox1i769g2 (Fig. 9A) and on superfamily SF10-28 comprising families F10, F17 and mF28 (Fig. 9B). The independent protein analysis fully confirms monophyly of the maturases in F03 and a well-supported clade of maturases in F01, now extended to include to mat-cox1i769g2, mat-atp1i850g2, the free-standing maturases in the mitogenomes of liverworts and a nuclear maturase copy (mat-nuc1) in the moss Physcomitrium (Fig. 9A). Notably, the extended SF01-22 superfamily also includes homologs in fungi having identical insertion sites and clustering with mat-cox1i652g2 in Coleochaete with high support. Vice versa, the extended F03 maturase clade likewise includes fungal mitochondrial maturases and a rhodophyte plastid maturase and, maybe even more notable, a cluster of nuclear maturases in tracheophytes (Fig. 9A).

As in the above case, the independent maturase phylogeny perfectly confirms the intron assignments to families F10 and F17 and adds the maturase-based family mF28 for a joint inclusion in superfamily SF10-28 (Fig. 9B). Particularly intriguing further cases for introns giving rise to fossil paralogues are “liverwort” introns atp1i989g2 and atp1i1050g2 (also present in hornworts), which are now jointly placed in mF28 (Fig. 4, Fig. 9B). An extended and significant sequence similarity (with perfect intron domain V and VI ends) of atp1i989g2 is present in the mitogenome of the lycophyte Phlegmariurus squarrosus embedding the trnI-rps11 region with all nine genes in the same direction fitting the intron orientation (Suppl. Figure 1L). Hence, it appears that a huge block of genes was inserted into an intergenic intron fossil paralogue of atp1i989g2. Interestingly, intron atp1i989g2 is absent in the early branching liverwort genera Treubia and Haplomitrium, which could have indicated a gain within the liverworts only after split of the Haplomitriopsida, but this now seems unlikely given the unorthodox Phlegmariurus fossil paralogue. Along the same lines, intron atp1i1050g2 has fossil intron paralogues not only in liverwort mitogenomes but also in the Phlegmariurus mtDNA behind cox1 in the spacer towards the trnW gene running in the opposite direction (Suppl. Figure 1M).

A further, very notable insight emerges from the maturase phylogeny: Mitochondrial group II introns with intron-borne maturases at the same insertion sites appear distributed across very distantly related lineages of eukaryotes. The most striking example concerns group II introns inserted into position 1147 of the cox1 gene. Solitary type intron cox1i1147g2 inhabits mitogenomes of the streptophyte alga Coleochaete, but also in chlorophytes, rhodophytes, fungi and metazoa. The associated mat-cox1i1147g2 RT-domains, the X-domain and the D/En domain are highly conserved. Similarly, the peculiar case of cox1i748g2 in Equisetum arvense [72], but not E. diffusum, having no similarity to its Chlorokybus counterpart is also found in the brown alga Pylaiella littoralis and the red alga Pyropia fucicola. A third intriguing example along those lines is cox2i373g2 of Coleochaete (now placed in mF26) that has maturase-free orthologs in mosses, hornworts and tracheophytes but forms a well-supported maturase-based clade with mat-cox2i373g2 in ascomycetous fungi, e.g., the endophytic symbiont Epichloe.

The remaining “solitary” introns lacking significantly similar paralogs

Altogether 61 streptophyte group II introns lacked paralogues with significant nucleotide sequence similarity precluding their assignment into our 25 core families. Extending the analysis to characteristic maturase similarities placed 15 of them in the four additional families mF26-mF29 and included another four previously solitary introns into superfamilies (cox1i44g2, cox1i245g2, cox1i769g2 and cox1i652g2). This leaves 42 streptophyte mitochondrial group II intron solitary, lacking both primary nucleotide similarity to paralogues and an intron-borne maturase of significant similarity to protein homologues (Fig. 5).

Many of the core solitary group II introns are presently very restricted in occurrence, most notably many introns presently identified in only one algal genus. This category, for example, also includes the unique trans-splicing introns nad3i84g2 and nad3i301g2 in the nad3 gene in Mesostigma viride. Among land plants this includes solitary introns exclusively restricted to mosses (atp1i1127g2 and cox1i1200g2), to hornworts (atp1i805g2, atp1i1019g2, cobi838g2, cox1i1298g2, cox2i281g2 and cox3i109g2) and to lycophytes (cobi693g2, cox1i227g2, cox1i266g2 and cox1i995g2), respectively, whereas no solitary introns are identified that are restricted to either ferns, gymnosperms or angiosperms.

Other solitary introns, however, are shared between at least two major plant clades (Fig. 5) e.g., nad4Li283g2 in liverworts and mosses or nad3i52g2 in hornworts and lycophytes. Yet others are shared among all euphyllophytes (nad5i1872g2), all tracheophytes (nad1i394g2, nad2i542g2 and nad7i917g2) or tracheophytes and at least one bryophyte clade (cox2i691g2, nad2i709g2 and nad7i140g2). An interesting further case is trnN-GUUi38g2 which is present in three algal classes likely close to the land plant lineage, absent in liverworts but clearly recognizable in pseudogenized form in mosses, hornworts and Phlegmariurus among the lycophytes (Suppl. Figure 1N).

Group II intron mobility outside of land plants

It is reasonable to assume that an ancestral state of group II introns with fully equipped intron-borne maturases is at the origin of their mobility and diversity (Fig. 1A). After the first discoveries of group II introns in bacteria [58, 73–75], their collection has grown immensely [76]. Bacterial group II introns are frequently associated with mobile genetic elements and/or present on plasmids [77] whereas their presence in essential genes is quite rare in bacteria [78]. Accordingly, their routes of dispersal are somewhat hard to trace and additionally complicated by the generally abundant horizontal gene transfer (HGT) among prokaryotes. Notably, prokaryotic group II introns lacking maturases are rare with “only a handful of ORF-less introns in bacteria” [79] and the few examples suggest a (likely quick) degeneration from maturase-bearing counterparts (Fig. 1B).

It should be noted that a “family” terminology for introns has already been used earlier in scoring bacterial group II intron occurrences and their retroelement behavior [79–81]. Among the research on group II introns in the bacterial world [82], the studies of mobile group II introns in Wolbachia endosymbionts [83] or of the highly mobile RmInt1 intron shaping the genome of Sinorhizobium meliloti and related α-proteobacteria [84] are typical cases in point on low taxonomic levels.

Several reports on variability, and likely mobility, of group II introns have also been published for eukaryotic organelle genomes, e.g. for the chloroplast genome of Euglena species [85], for mitochondria of different isolates of the brown alga Pylaiella littoralis [86] or in chloroplasts and mitochondria of the red algal genus Porphyridium [87, 88]. On somewhat higher taxonomic levels, significant group II intron diversities have been reported for the plastid genomes of cryptophytes [89] or the mitogenomes of diatoms [90, 91]. Extraordinary similarity of cyanobacterial introns in Porphyra [92] or the chloroplast of Chlamydomonas [93] or Euglena [94] have been discussed as horizontal transfer events. Likewise, horizontal transfers were also invoked as the likely cause in other cases of striking intron similarities, e.g. between those of diatoms and a haptophyte [95], for diatoms as likely group II intron donors into a raphidophycean flagellate genus Chattonella [96] or to explain the discovery of a first group II intron in a bilaterian mitogenome [97].

The plant mitochondrial group II intron family concept

The extraordinary dynamics of plant mitochondrial introns is in stark contrast not only to the overall stability of the chloroplast intron complements [98, 99], but also to the stasis of nuclear introns since 500 million years and more [100]. Many review articles [21, e.g. 101] and reports on newly completed plant mitogenomes, respectively, often contain comparative summaries and updates on the striking diversity of plant mitogenome makeups with respect to their gene and intron complements [11, 12, 65, 102].

Here, we have added on the previous comparative compilations by providing the hitherto most extensive overview on mitochondrial group II introns by including mitogenome analyses of all embryophyte and streptophyte algae lineages, also taking care of pseudogenization events and the existence of fossil introns in intergenic spacer regions. More importantly, we provide a concept for categorizing group II introns into families of related paralogs as a foundation to explore intron “copying” retrotransposition events during plant evolution. Of altogether 161 streptophyte mitochondrial group II introns, we have assigned 100 into 25 “core” families of minimally two and up to 13 intron paralogues based on their nucleotide sequences alone (Fig. 4). The family assignments may certainly be subject to further changes with newly discovered mitochondrial group II introns in streptophytes added to the 25 defined families or with new families or superfamilies (Fig. 5) being created by newly identified introns linking existing families. While we do not expect many more introns to be discovered in the land plant (embryophyte) lineage, the highly diverse streptophyte algae are presently still underrepresented in the sampling for mitogenomes with only twelve genera representing the five classically distinguished classes and will likely reveal many more intron paralogues in the future.

Somewhat surprisingly, our intron inventory does not reveal a particular strong affiliation of land plants with Zygnematophyceae algae, currently considered to be the sister lineage of embryophytes despite this class being best sampled for mitogenomes with five genera representing four families. Similarly, mitochondrial group II introns have contributed to a model phylogeny of land plants assuming liverworts as sister to all other embryophytes and hornworts as sister to tracheophytes [8, 103, 104]. Alternative datasets, however, now favor the concept of monophyletic bryophytes with a sister group relationship of mosses and liverworts [105–108]. Both phylogenies must postulate massive gains and losses of mitochondrial group II introns in the early embryophyte lineages (Fig. 10). Testing the alternative phylogenetic concepts by Maximum Parsimony searches for the 101 mitochondrial group II introns present in embryophytes, we still find weak support for the NLE/HT topology requiring nine steps less than the monophyletic Bryophyte topology (Fig. 10).

Similarities of intron paralogues in the absence of intron-borne maturases

The dynamics of group II introns in plant mitochondrial vs. chloroplast genomes is puzzling at first sight. Key to an explanation for this observation may be the co-evolution of group II introns with concomitant splicing factors and the ancestrally intron-borne maturases certainly play an initial key role. In flowering plants, only one maturase each remains present in the organelle DNAs – “MatK” in the group II intron disrupting the trnK gene (mat-trnKi37g2) and “MatR” in the terminal group II intron of the nad1 gene of mitochondria (mat-nad1i728g2). Accordingly, the lack of evidently recent group II intron mobility in angiosperm organelles, and in fact even spermatophytes, is likely no coincidence. Notably, whereas MatK is also the only chloroplast maturase conserved in the other land plant clades, many more mitochondrial maturases exist in early land plant lineages (Fig. 4 and Fig. 5), which have likely contributed to intron retrotranspositions during establishment of the ancestral plant lineages. Additionally, four nuclear-encoded maturases (Fig. 1C) have been shown to functionally affect diverse mitochondrial, but not chloroplast, group II introns in Arabidopsis [15–20].

In some cases, it appears indeed suggestive that a maturase-containing intron paralogue is the likely source of the other paralogues in a family. Key examples are the large intron families F01 and F03 where five out of 13 and four out of eight paralogues, respectively, carry maturases at least in taxa representing early branching lineages (Fig. 4). Other examples are family F07 with its likely ancestral “mother intron” atp9i87g2, family F10 with cobi824g2 or nad1i669g2, family F13 with cox1i178g2 and family F14 with cox1i511g2 carrying a maturase in the alga Coleochaete scutata.

Vice versa, however, no traces of maturases are discernible in many other families (Fig. 1F) and even despite the existence of related paralogues in the same plant lineages such as the “liverwort-lycophyte” family F05, the “liverwort-lycophyte” family F08, the “moss-hornwort” family F11, the “hornwort” family F16 or the large and diverse family F02 (Fig. 4). Other striking examples include family F06 where it is highly likely that atp1i361g2 has arisen from rps3i249g2 exclusively in ferns, family F09 where the same holds true for rps1i25g2 likely originating from rpl2i846g2 and family F04 dominated by hornwort paralogues including recently emerged twintron arrangements (Fig. 4). The latter represents the most intriguing case of 98% sequence similarity between the maturase-free paralogues atp1i1050g2ii1536g2 and cox2i98g2, indicating a very recent origin of cox2i98g2 in the genus Anthoceros originating from the internal intron of the atp1 twintron present in all hornwort mitogenomes.

Similarly, the lack of intron-borne maturases holds true for most of the solitary, S-type introns lacking any evident similarity to other introns (Fig. 5). Only 14 of 61 S-type introns show up in more than one plant clade and any evidence for an external origin via HGT is missing. Like the above cases of maturase-less, but clearly related, paralogues their origin from maturase-bearing ancestors would necessitate a surprisingly quick degeneration of intron-borne reading frames that is incompatible with the overall slow sequence evolution in plant mitogenomes, which reveal the retention of pseudogene traces including former maturases in many other cases [14, 109, 110]. The numerous examples of group II intron fossils reported here are a further case in point (Fig. 1G). After an early discovery of a degenerate group II intron in Chlamydomonas [111] this issue may have received less attention than deserved in the exploration of newly determined organelle genomes.

The participation of external protein factors rather than only intron-borne maturases promoting retrotransposition, possibly including nuclear-encoded maturases in the early plant lineages [112], seems a more likely explanation for early group II intron propagations in plant mitochondria (Fig. 1E).

Nuclear splicing factors acting on paralogues in the same family and an outlook

Quite interestingly, we found that in several cases nuclear-encoded splicing factors have already been identified (Fig. 1E), which simultaneously act on related intron paralogues that we have now grouped into families. The nuclear-encoded splicing factor ODB1 (“Organelle DNA Binding 1”) containing a RAD52-like domain has been found to promote splicing of the two angiosperm intron paralogues in family F07, nad1i477g2 and nad2i156g2 [64]. Likewise, the two F09 introns present in angiosperms, ccmFci829g2 and rpl2i846g2, depend on WTF9, a splicing factor containing a PORR (Plant Organelle RNA Recognition) domain [66]. In family F01, even three different intron paralogues present in angiosperms (nad5i230g2, nad7i209g2 and rps3i74g2) rely on splicing factor OZ2, a RanBP2-type zinc finger protein [50]. Finally, the two F17 introns present in angiosperms (nad1i728g2 and nad4i461g2 have also been shown to both depend on three different splicing factors: the two pentatricopeptide repeat proteins DEK55 [70] and EMP8 [69] and the mTERF-type (Mitochondrial Transcription Termination Factor) protein SMK3 [71]. It will be highly interesting to see whether the splicing factor functionalities demonstrated in the model angiosperm species Arabidopsis thaliana or Zea mays are conserved deep in flowering plant evolution or even beyond angiosperms once adequate model taxa are established among gymnosperms, ferns and lycophytes.

Possibly even more interesting for the future developments will be whether these factors, aside from their evident function in facilitating proper RNA structures for forward-splicing, may also have contributed to the ancient retrotranspositions creating the related paralogues. The rapidly growing amount of data on introns in plant organelle genomes and on the functional characterization of nuclear splicing factors calls makes review compilations [46, 47] quickly outdated and rather calls for expert database systems in the future. Towards that end, we will soon establish a relational database system for plant group organelle introns (‘G12Base’) and their functionally characterized splicing factors (‘SpliciFacts’). In the long run, such an endeavor should clearly not remain restricted to group II introns in streptophyte mitogenomes but will certainly be extended in multiple ways. Firstly, it will certainly consider deeper phylogenetic branches in chlorophytes and beyond, already now justified for example by the tremendous variability of introns already identified in Ulvophyceae, including a likely mitochondrion-to-chloroplast intron transfer [113] or the examples of intron mobility in red algae [87]. Along the same lines, and in spite their striking conservation of insertion sites among land plants (and despite their lack of primary sequence conservation, which precludes any family assignments compared to their mitochondrial counterparts), chloroplast introns will not go unnoticed. Chloroplast intron dynamics on deeper taxonomic levels may be highly interesting. The gain of introns discovered early to define an early split among algae is but one interesting example [114], the proliferation of group II introns in the algal genus Oedocladium (Chlorophyceae) another [115].

Finally, and although group I introns rely on entirely different mechanisms of splicing and propagation and are rare in flowering plants, the example of the “rampant invader” group I intron cox1i726g1 in angiosperms of likely fungal origin [116–118] and the generally much more frequent occurrence of group I introns in early branching green (and non-green) lineages clearly warrants the future database inclusion of reports on their occurrence, invasion of new loci and likely mechanisms of splicing, too. In any case, we consider the diversity of organelle introns in the plant mitochondrial lineage a highly interesting evolutionary playground clearly worth of further investigation.

Primary nucleotide sequence data sampling

Streptophyte group II intron sequences were extracted from GenBank accessions available at the National Center for Biotechnology Information (see Table 1 for the primary selection). The sampling also included the recently determined mitogenome of the leptosporangiate fern Haplopteris ensiformis, which revealed introns that are not present in the two available mitogenomes of eusporangiate ferns [119]. Introns were labeled using the standard nomenclature proposal indicating the respective gene and the homologous nucleotide position in the Marchantia polymorpha mitogenome reference [56, 57]. Introns were manually re-checked for conserved sequence motifs at their terminal 5’- and 3’-splicing sites and corrections were applied where necessary (see examples under results). Mitogenome accessions were also checked for potentially overlooked group II introns by taking advantage of searching for their characteristic and mostly conserved domains V of 34 bp [58]. Individual introns were used as queries in sensitive BLASTN searches to identify significant similarities also in intergenic sequences (IGS) of mitogenomes to identify cases of “intron fossils” (see main text).

Table 1

List of Species-data used for different analyses List of taxa and mitogenomes for detailed initial scoring of mitochondrial group II intron presence. The respective numbers of group II introns are indicated. Taxon selection was aimed to maximize phylogenetic diversity with the respective clades. Care was taken to extract intron sequences with proper 5’- and 3’-termini. Asterisks indicate species not further considered for detailed group II intron sequence analyses owing to redundancy with other species within the respective group. A plus sign (+) in the IGS analysis column indicated careful inspection of intergenic sequences for the presence of intron fossils
Major lineage		Species	Acc. number	Mt- group II introns
ALG	Coleochaetophyceae	Chaetosphaeridium globosum	NC_004118.1	2 +
	Coleochaetophyceae	Coleochaete scutata	NC_045180.1	26 +
	Mesostigmatophyceae	Mesostigma viride	AF353999	3
	Chlorokybophyceae	Chlorokybus atmophyticus	EF463011.1	14 +
	Charophyceae	Nitella hyalina	NC_017598	8
	Charophyceae	Chara vulgaris	NC_005255.1	13 +
	Klebsormidiophyceae	Klebsormidium flaccidum	KP165386	4
	Zygnematophyceae	Closterium baillyanum	NC_022860.1	12 +
	Zygnematophyceae	Roya obtusa	NC_22863.1	2
	Zygnematophyceae	Entransia fimbriata	NC_22861.1	2
	Zygnematophyceae	Zygnema circumcarinatum	MT040698.1	10 +
	Zygnematophyceae	Gonatozygon brebissonii	NC_046951.1	6
LIV	Haplomitriopsida	Treubia lacunosa	NC_016122.1	21 +
	Marchantiopsida	Marchantia polymorpha	MK202951.1	25 +
	Jungermanniopsida	Calypogeia fissa	NC_035980	19 +
	Jungermanniopsida	Pleurozia purpurea	NC_013444	22
MOS	Sphagnopsida	Flatbergium novo-caledoniae*	KU725492	26 +
	Sphagnopsida	Sphagnum palustre	NC_024521.1	26 +
	Bryopsida	Physcomitrium patens	NC_007945.1	24 +
	Bryopsida	Anomodon rugelii	NC_016121	24
	Bryopsida	Ulota crispa	NC_031393	24
HOR	Leiosporocerotopsida	Leiosporoceros dussii	NC_039751.1	34 +
	Anthocerotopsida	Anthoceros agrestis	NC_049004.1	39 +
	Anthocerotopsida	Nothoceros aenigmaticus	EU660574	28 +
	Anthocerotopsida	Phaeoceros laevis	NC_013765	31
LYC	Lycopodiopsida	Phlegmariurus squarrosus	NC_017755.1	37 +
	Lycopodiopsida	Isoetes engelmannii	FJ390841.1, FJ176330.1, FJ010859.1, FJ628360.1	23
	Lycopodiopsida	Selaginella moellendorffii	JF338143.1 – JF338147.1	33
FER	Polypodiopsida	Ophioglossum californicum	NC_030900.1	20
	Polypodiopsida	Psilotum nudum	NC_030952 & KX171639.1	24 +
	Polypodiopsida	Haplopteris ensiformis	OM867545- OM867553	24 +
GYM	Ginkgoopsida	Ginkgo biloba	NC_027976.1	25 +
	Cycadopsida	Cycas taitungensis	NC_010303	25 +
	Gnetopsida	Welwitschia mirabilis	NC_029130	10
ANG	Magnoliopsida	Amborella trichopoda	KF754799- KF754803	24
	Magnoliopsida	Nelumbo nucifera	NC_030753	23 +
	Magnoliopsida	Zea mays	NC_007982.1	22 +
	Magnoliopsida	Arabidopsis thaliana	NC_037304.1	23 +

Maturase sequence sampling

Mitochondrial and nuclear encoded maturases were collected independently of the primary group II intron sequences, extending an earlier sampling [112] with sequences from streptophyte algae, also including evident maturases that have not been annotated in sequence entries (see Suppl. Figure 2). The conservation and disintegration of the respective maturases was evaluated with respect to subdomains RT0 - RT7, including signature peptide motifs [120] as well as the conservation of the X-domain and signature peptide motifs for DNA-endonucleases [121]. The ultimate set of streptophyte mitochondrial maturases was used for BLASTP searches for homologs with significant similarities also outside of the plant mitochondrial lineage in other eukaryotic lineages: Metazoa, Fungi, Stramenopiles, Oomycetes, Cryptophytes, Rhodophytes and Chlorophytes. The top ten significant hits for each query were retained after removal of duplicates for alignment and phylogenetic tree construction (see Suppl. Figure 3).

The “most-distant ortholog” concept to define core primary group II intron families

Starting from occasional observations that streptophyte mitochondrial group II intron paralogs in different insertion sites may share more sequence similarities than orthologs in phylogenetically distant clades, we developed a “most-distant-ortholog”-concept (MDO) to cluster the collection of 161 introns (781 representatives in the initial collection of organelle genomes, see Table 1) into families of related paralogs. The full collection of group II intron sequences was initially compared against itself using the stand-alone BLAST + tool [122]. Sensitive BLASTN search parameters were set as follows: word size of 7, match and mismatch values of 2 and − 3 and penalty values of 5 and 2 for gap opening and extension and a random expectancy cutoff value of 1e-5, respectively. Results were evaluated using the RStudio package “dplyr” (tidyverse library) [123] of the. Initial clustering following the MDO criterion used the BLASTN search bit score values for most distant orthologs as cut-off to consider paralog hits with higher bit scores for inclusion into a given family.

Refinement of the group II intron family concepts

Intron homologs present at evidently identical insertion sites that failed to be included in the initial orthologous clustering (“false negatives”) owing to large, “unbridged” phylogenetic distances owing to (multiple) losses (e.g., rps3i74g2 present in charophyceaen algae, hornworts and tracheophytes) were included among the respective family sampling if not showing higher similarity to members of another family. Taxonomically isolated introns receiving only low bit scores of less than 200 were manually re-inspected and only considered for family inclusion when independent sequence similarities were identified in multiple intron regions and/or clearly responsible for conserved RNA secondary structures (see next paragraph). Vice versa, we wished to exclude false positive family assignments for paralogs with arbitrarily high bit scores owing to taxon-specific sequence evolution e.g., due to independently arising repeat motifs.

Group II intron alignments, phylogenetic analyses and RNA structure modeling

The clustering of intron orthologs and related paralogs into families was manually re-evaluated and re-checked by careful inspection of the respective full alignments. Intron sequences were automatically aligned with MAFFT [124] followed by manual adjustment where necessary. Sequence conservation was visualized with GeneDoc 2.7 (https://genedoc.software.informer.com/). Alignment positions with 60% site coverage identified in MEGA 7.0 [125] were extracted for further phylogenetic analyses and modelling of RNA-secondary structures. The filtered group II intron family alignments were used for Maximum Likelihood phylogenetic analyses with IQ-TREE 2 [126]. Substitution model testing was done automatically, and node support was evaluated with 1,000 ultra-fast bootstrap and SH-aLRT repetitions, respectively. Phylogenetic Trees were visualized and edited with TreeGraph2 [127]. RNA secondary structures were modelled manually following established group II intron structure modeling routines [76] with predictions for substructure generated by the mxfold2 application [128]. Nucleotide positions from group II intron family alignments with 60% site coverage sites were preferably used for folding compared to non-conserved sites. RNA secondary structures were annotated in Dot-Bracket annotation and visualized with VARNA [129].

Evaluation of streptophyte mitochondrial group II intron presence and absence

Presence and absence of the full collection of streptophyte mitochondrial group II intron paralogues were coded as ‘1’ and ‘0’, respectively. The R4.1.2 packages “dplyr” and “eulerr” were used for evaluation and plotting of group II intron distributions in the major land plant lineages to create Venn/Euler diagrams [130, 131]. Intron distribution data was converted into the FASTA format to model gain and loss scenarios empirically (gain/loss costs of 4 and 1, respectively) for group II introns for different plant phylogeny scenarios with GLOOME [132].

Consent of participation

Not applicable

Consent of publication

Not applicable

Availability of data and materials

Nucleotide sequence data including accessions used in this study are listed in manuscript Table 1. Nucleotide sequence data is available at the National Center for Biotechnology Information (NCBI) cited 2022 September 26. Available from https://www.ncbi.nlm.nih.gov/nucleotide/.

Competing interests

The authors have no competing interests as defined by BMC, or other interests that might be perceived to influence the results and/or discussion reported in this paper. The results/data/figures in this manuscript have not been published elsewhere, nor are they under consideration (from you or one of your Contributing Authors) by another publisher.

Funding

No funding

Author contribution

SZ established bioinformatic pipelines and collected data. VK and SZ analyzed data and prepared figures. VK wrote the main manuscript text and all authors reviewed the manuscript.

Acknowledgements

We are grateful for inspiring discussions with members of the IZMB Bonn Molecular Evolution group and to Dietrich Dmitriev for comments on the manuscript.

Small ID, Schallenberg-Rüdinger M, Takenaka M, Mireau H, Ostersetzer-Biran O. Plant organellar RNA editing: what 30 years of research has revealed. Plant J. 2020;101:1040–56.
Li-Pook-Than J, Bonen L. Multiple physical forms of excised group II intron RNAs in wheat mitochondria. Nucleic Acids Res. 2006;34:2782–90.
Ngu M, Massel K, Bonen L. Group II introns in wheat mitochondria have degenerate structural features and varied splicing pathways. Int J Biochem Cell Biol. 2017;91:156–67.
Massel K, Silke JR, Bonen L. Multiple splicing pathways of group II trans-splicing introns in wheat mitochondria. Mitochondrion. 2016. https://doi.org/10.1016/j.mito.2016.03.002.
Malek O, Knoop V. Trans-splicing group II introns in plant mitochondria: the complete set of cis-arranged homologs in ferns, fern allies, and a hornwort. RNA. 1998;4:1599–609.
Qiu YL, Palmer JD. Many independent origins of trans splicing of a plant mitochondrial group II intron. J Mol Evol. 2004;59:80–9.
Guo W, Zhu A, Fan W, Adams RP, Mower JP. Extensive shifts from cis to trans splicing of gymnosperm mitochondrial introns. Mol Biol Evol. 2020. https://doi.org/10.1093/molbev/msaa029.
Groth-Malonek M, Pruchner D, Grewe F, Knoop V. Ancestors of trans-splicing mitochondrial introns support serial sister group relationships of hornworts and mosses with vascular plants. Mol Biol Evol. 2005;22:117–25.
Knoop V, Altwasser M, Brennicke A. A tripartite group II intron in mitochondria of an angiosperm plant. Mol Gen Genet MGG. 1997;255:269–76.
Kück U, Schmitt O. The Chloroplast Trans -Splicing RNA–Protein Supercomplex from the Green Alga Chlamydomonas reinhardtii. Cells (MDPI). 2021;10:290.
Yu R, Sun C, Liu Y, Zhou R. Shifts from cis- to trans-splicing of five mitochondrial introns in Tolypanthus maclurei. PeerJ. 2021;9:e12260.
Kan S-L, Shen T-T, Gong P, Ran J-H, Wang X-Q. The complete mitochondrial genome of Taxus cuspidata (Taxaceae): eight protein-coding genes have transferred to the nuclear genome. BMC Evol Biol. 2020;20:10.
Zumkeller SM, Gerke P, Knoop V. A functional twintron, “zombie” twintrons and a hypermobile group II intron invading itself in plant mitochondria. Nucleic Acids Res. 2020;48:2661–75.
Ahlert D, Piepenburg K, Kudla J, Bock R. Evolutionary origin of a plant mitochondrial group II intron from a reverse transcriptase/maturase-encoding ancestor. J Plant Res. 2006;119:363–71.
Keren I, Tal L, Des Francs-Small CC, Araújo WL, Shevtsov S, Shaya F, et al. NMAT1, a nuclear-encoded maturase involved in the trans-splicing of nad1 intron 1, is essential for mitochondrial complex i assembly and function. Plant J. 2012;71:413–26.
Keren I, Bezawork-Geleta A, Kolton M, Maayan I, Belausov E, Levy M, et al. AtnMat2, a nuclear-encoded maturase required for splicing of group II introns in Arabidopsis mitochondria. RNA. 2009;15:2299–311.
Shevtsov-Tal S, Best C, Matan R, Aldrin Chandran S, Brown GG, Ostersetzer‐Biran O. nMAT3 is an essential maturase splicing factor required for holo‐complex I biogenesis and embryo‐development in Arabidopsis thaliana plants. Plant J. 2021;:tpj.15225.
Cohen S, Zmudjak M, Colas Des Francs-Small C, Malik S, Shaya F, Keren I, et al. NMAT4, A maturase factor required for nad1 pre-mRNA processing and maturation, is essential for holocomplex i biogenesis in Arabidopsis mitochondria. Plant J. 2014;78:253–68.
Schmitz-Linneweber C, Lampe M-KK, Sultan LD, Ostersetzer-Biran O. Organellar maturases: A window into the evolution of the spliceosome. Biochim Biophys Acta. 2015;1847:798–808.
Mohr G, Lambowitz AM. Putative proteins related to group II intron reverse transcriptase/maturases are encoded by nuclear genes in higher plants. Nucleic Acids Res. 2003;31:647–52.
Mower JP. Variation in protein gene and intron content among land plant mitogenomes. Mitochondrion. 2020;53:203–13.
Pruchner D, Nassal B, Schindler M, Knoop V. Mosses share mitochondrial group II introns with flowering plants, not with liverworts. Mol Genet Genomics. 2001;266:608–13.
Pruchner D, Beckert S, Muhle H, Knoop V. Divergent intron conservation in the mitochondrial nad2 gene: signatures for the three bryophyte classes (mosses, liverworts, and hornworts) and the lycophytes. J Mol Evol. 2002;55:265–71.
Knie N, Grewe F, Knoop V. Monilophyte mitochondrial rps1 genes carry a unique group II intron that likely originated from an ancient paralog in rpl2. RNA. 2016;22:1338–48.
Zumkeller SM, Knoop V, Knie N. Convergent Evolution of Fern-Specific Mitochondrial Group II Intron atp1i361g2 and Its Ancient Source Paralogue rps3i249g2 and Independent Losses of Intron and RNA Editing among Pteridaceae. Genome Biol Evol. 2016;8:2505–19.
Mukhopadhyay J, Hausner G. Organellar introns in fungi, algae and plants. Cells (MDPI). 2021;10:1–35.
Lambowitz AMAM, Zimmerly S. Mobile group II introns. Annu Rev Genet. 2004;38:1–35.
Enyeart PJ, Mohr G, Ellington AD, Lambowitz AM. Biotechnological applications of mobile group II introns and their reverse transcriptases: Gene targeting, RNA-seq, and non-coding RNA analysis. Mob DNA. 2014;5:1–19.
Robart AR, Zimmerly S. Group II intron retroelements: function and diversity. Cytogenet Res. 2005;110:589–97.
Bonen L, Vogel J. The ins and outs of group II introns. Trends Genet. 2001;17:322–31.
McNeil BA, Semper C, Zimmerly S. Group II introns: versatile ribozymes and retroelements. Wiley Interdiscip Rev RNA. 2016;7:341–55.
Zimmerly S, Semper C. Evolution of group II introns. Mob DNA. 2015;6:7.
Lambowitz AM, Zimmerly S. Group II Introns: Mobile Ribozymes that Invade DNA. Cold Spring Harb Perspect Biol. 2011;3:a003616–6.
Lambowitz AM, Belfort M. Mobile Bacterial Group II Introns at the Crux of Eukaryotic Evolution. Microbiol Spectr. 2015;3.
Pyle AM. Group II Intron Self-Splicing. Annu Rev Biophys. 2016;45:183–205.
Pyle AM, Fedorova O, Waldsich C. Folding of group II introns: a model system for large, multidomain RNAs? Trends Biochem Sci. 2007;32:138–45.
Patel S, Fedorova O. Ribozyme Structural Elements: Group II Introns. In: Reference Module in Life Sciences. Elsevier; 2021.
Fedorova O, Zingler N. Group II introns: structure, folding and splicing mechanism. Biol Chem. 2007;388:665–78.
Peters JK, Toor N. Group II intron lariat: Structural insights into the spliceosome. RNA Biol. 2015;12:913–7.
Chalamcharla VR, Curcio MJ, Belfort M. Nuclear expression of a group II intron is consistent with spliceosomal intron ancestry. Genes Dev. 2010;24:827–36.
Smathers CM, Robart AR. The mechanism of splicing as told by group II introns: Ancestors of the spliceosome. Biochim Biophys Acta - Gene Regul Mech. 2019. https://doi.org/10.1016/j.bbagrm.2019.06.001.
Novikova O, Belfort M. Mobile Group II Introns as Ancestral Eukaryotic Elements. Trends Genet. 2017;33:773–83.
Costa M, Fontaine JM, Loiseaux-De Goër S, Michel F. A group II self-splicing intron from the brown alga Pylaiella littoralis is active at unusually low magnesium concentrations and forms populations of molecules with a uniform conformation. J Mol Biol. 1997;274:353–64.
Barthet MM, Pierpont CL, Tavernier EK. Unraveling the role of the enigmatic MatK maturase in chloroplast group IIA intron excision. Plant Direct. 2020;4:e00208.
Inaba-Hasegawa K, Ohmura A, Nomura M, Sugiura M. Development of an In Vitro Chloroplast Splicing System: Sequences Required for Correct pre-mRNA Splicing. Plant Cell Physiol. 2021;62:1311–20.
Brown GG, Colas des Francs-Small C, Ostersetzer-Biran O. Group II intron splicing factors in plant mitochondria. Front Plant Sci. 2014;5 February:35.
de Longevialle AF, Small ID, Lurin C. Nuclearly encoded splicing factors implicated in RNA splicing in higher plant organelles. MolPlant. 2010;3:691–705.
Qian J, Li M, Zheng M, Hsu Y-F. Arabidopsis SSB1, a Mitochondrial Single-Stranded DNA-Binding Protein, is Involved in ABA Response And Mitochondrial RNA Splicing. Plant Cell Physiol. 2021. https://doi.org/10.1093/PCP/PCAB097.
Cao S-K, Liu R, Sayyed A, Sun F, Song R, Wang X, et al. Regulator of Chromosome Condensation 1-Domain Protein DEK47 Functions on the Intron Splicing of Mitochondrial Nad2 and Seed Development in Maize. Front Plant Sci. 2021;0:1612.
Bentolila S, Gipson AB, Kehl AJ, Hamm LN, Hayes ML, Mulligan RM, et al. A RanBP2-type zinc finger protein functions in intron splicing in Arabidopsis mitochondria and is involved in the biogenesis of respiratory complex I St ephane. Nucleic Acids Res. 2021;:1–17.
Lee K, Leister D, Kleine T. Arabidopsis Mitochondrial Transcription Termination Factor mTERF2 Promotes Splicing of Group IIB Introns. Cells. 2021;10:315.
Yang Y-Z, Ding S, Wang Y, Wang H-C, Liu X-Y, Sun F, et al. PPR20 Is Required for the cis-Splicing of Mitochondrial nad2 Intron 3 and Seed Development in Maize. Plant Cell Physiol. 2020;61:370–80.
Ohyama K, Oda K, Ohta E, Takemura M. Gene organization and evolution of introns of a liverwort, Marchantia polymorpha, mitochondrial genome. In: Brennicke A, Kück U, editors. Weinheim: VCH Verlagsgesellschaft; 1993. pp. 115–29.
Lippok B, Brennicke A, Wissinger B. Differential RNA editing in closely related introns in Oenothera mitochondria. Mol Gen Genet. 1994;243:39–46.
Vangerow S, Teerkorn T, Knoop V. Phylogenetic information in the mitochondrial nad5 gene of pteridophytes: RNA editing and intron sequences. Plant Biol. 1999;1:235–43.
Dombrovska E, Qiu Y-L. Distribution of introns in the mitochondrial gene nad1 in land plants: phylogenetic and molecular evolutionary implications. Mol Phylogenet Evol. 2004;32:246–63.
Knoop V. The mitochondrial DNA of land plants: peculiarities in phylogenetic perspective. Curr Genet. 2004;46:123–39.
Knoop V, Kloska S, Brennicke A. On the identification of group II introns in nucleotide sequence data. J Mol Biol. 1994;242:389–96.
Marienfeld JR, Newton KJ. The Maize Ncs2 Abnormal Growth Mutant Has A Chimeric Nad4-Nad7 Mitochondrial Gene and Is Associated with Reduced Complex-I Function. Genetics. 1994;138:855–63.
Knoop V, Ehrhardt T, Lättig K, Brennicke A. The gene for ribosomal protein S10 is present in mitochondria of pea and potato but absent from those of Arabidopsis and Oenothera. Curr Genet. 1995;27:559–64.
Zumkeller SM, Polsakiewicz M, Knoop V. Rickettsial DNA invasions and a scrambled rRNA cluster with a trans-splicing group I intron: The highly unorthodox mitogenome of the fern Haplopteris ensiformis. Submitted. 2022.
Volkmar U, Groth-Malonek M, Heinrichs J, Muhle H, Polsakiewicz M, Knoop V. Exclusive conservation of mitochondrial group II intron nad4i548 among liverworts and its use for phylogenetic studies in this ancient plant clade. Plant Biol (Stuttg). 2012;14:382–91.
Groth-Malonek M, Rein T, Wilson R, Groth H, Heinrichs J, Knoop V. Different fates of two mitochondrial gene spacers in early land plant evolution. Int J Plant Sci. 2007;168:709–17.
Gualberto JM, Le Ret M, Beator B, Kühn K. The RAD52-like protein ODB1 is required for the efficient excision of two mitochondrial introns spliced via first-step hydrolysis. Nucleic Acids Res. 2015;43:6500–10.
Gerke P, Szövényi P, Neubauer A, Lenz H, Gutmann B, McDowell R, et al. Towards a plant model for enigmatic U-to-C RNA editing: the organelle genomes, transcriptomes, editomes and candidate RNA editing factors in the hornwort Anthoceros agrestis. New Phytol. 2020;225:1974–92.
Francs-Small CC des, Kroeger T, Zmudjak M, Ostersetzer-Biran O, Rahimi N, Small I, et al. A PORR domain protein required for rpl2 and ccmF(C) intron splicing and for the biogenesis of c-type cytochromes in Arabidopsis mitochondria. Plant J Cell Mol Biol. 2012;69:996–1005.
Turmel M, Otis C, Lemieux C. Complete mitogenome of the streptophyte green alga Coleochaete scutata (Coleochaetophyceae). Mitochondrial DNA Part B Resour. 2019;4:4209–10.
Grewe F, Zhu A, Mower JP. Loss of a Trans-Splicing nad1 Intron from Geraniaceae and Transfer of the Maturase Gene matR to the Nucleus in Pelargonium. Genome Biol Evol. 2016;8:3193–201.
Sun F, Zhang X, Shen Y, Wang H, Liu R, Wang X, et al. The pentatricopeptide repeat protein EMPTY PERICARP8 is required for the splicing of three mitochondrial introns and seed development in maize. Plant J. 2018. https://doi.org/10.1111/tpj.14030.
Ren RC, Yan XW, Zhao YJ, Wei YM, Lu X, Zang J, et al. The novel E-subgroup pentatricopeptide repeat protein DEK55 is responsible for RNA editing at multiple sites and for the splicing of nad1 and nad4 in maize. BMC Plant Biol. 2020;20:1–15.
Pan Z, Ren X, Zhao H, Liu L, Tan Z, Qiu F. A mitochondrial transcription termination factor, ZmSmk3, is required for nad1 Intron4 and nad4 Intron1 splicing and kernel development in maize. G3 Genes, Genomes. Genet. 2019;9:2677–86.
Bégu D, Araya A. The horsetail Equisetum arvense</i > mitochondria share two group I introns with the liverwort Marchantia, acquired a novel group II intron but lost intron-encoded ORFs. Curr Genet. 2009;55:69–79.
Ferat JL, Michel F. Group II self-splicing introns in bacteria. Nature. 1993;364:358–61.
Ferat JL, Legouar M, Michel F. Multiple Group-Ii Self-Splicing Introns in Mobile Dna from Escherichia-Coli. Comptes Rendus l Acad des Sci Ser Iii-Sciences la Vie-Life Sci. 1994;317:141–8.
Knoop V, Brennicke A. Evidence for a group II intron in Escherichia coli inserted into a highly conserved reading frame associated with mobile DNA sequences. Nucleic Acids Res. 1994;22:1167–71.
Candales MA, Duong A, Hood KS, Li T, Neufeld RAE, Sun R, et al. Database for bacterial group II introns. Nucleic Acids Res. 2012;40.
Dai L, Zimmerly S. The dispersal of five group II introns among natural populations of Escherichia coli. RNA. 2002;8:1294–307.
Ferat J-L, Le Gouar M, Michel F. A group II intron has invaded the genus Azotobacter and is inserted within the termination codon of the essential groEL gene. Mol Microbiol. 2003;49:1407–23.
Simon DM, Clarke NAC, McNeil BA, Johnson I, Pantuso D, Dai L, et al. Group II introns in eubacteria and archaea: ORF-less introns and new varieties. RNA. 2008;14:1704–13.
Toor N, Zimmerly S. Identification of a family of group II introns encoding LAGLIDADG ORFs typical of group I introns. RNA. 2002;8:1373–7.
Dai L, Zimmerly S. Compilation and analysis of group II intron insertions in bacterial genomes: evidence for retroelement behavior. Nucleic Acids Res. 2002;30:1091–102.
Martínez-Abarca F, Toro N. Group II introns in the bacterial world. Mol Microbiol. 2000;38:917–26.
Leclercq S, Giraud I, Cordaux R. Remarkable Abundance and Evolution of Mobile Group II Introns in Wolbachia Bacterial Endosymbionts. Mol Biol Evol. 2011;28:685–97.
Toro N, Martínez-Abarca F, Molina-Sánchez MD, García-Rodríguez FM, Nisa-Martínez R. Contribution of Mobile Group II Introns to Sinorhizobium meliloti Genome Evolution. Front Microbiol. 2018;9:627.
Thompson MD, Copertino DW, Thompson E, Favreau MR, Hallick RB. Evidence for the late origin of introns in chloroplast genes from an evolutionary analysis of the genus Euglena. Nucleic Acids Res. 1995;23:4745–52.
Ikuta K, Kawai H, Müller DG, Ohama T. Recurrent invasion of mitochondrial group II introns in specimens of Pylaiella littoralis (brown alga), collected worldwide. Curr Genet. 2008;53:207–16.
Perrineau M-M, Price DC, Mohr G, Bhattacharya D. Recent mobility of plastid encoded group II introns and twintrons in five strains of the unicellular red alga Porphyridium. PeerJ. 2015;3:e1017.
Kim D, Lee JM, Cho CH, Kim EJ, Bhattacharya D, Yoon HS. Group II intron and repeat-rich red algal mitochondrial genomes demonstrate the dynamic recent history of autocatalytic RNAs. BMC Biol. 2022;20:1–16.
Khan H, Archibald JM. Lateral transfer of introns in the cryptophyte plastid genome. Nucleic Acids Res. 2008;36:3043–53.
Pogoda CS, Keepers KG, Hamsher SE, Stepanek JG, Kane NC, Kociolek JP. Comparative analysis of the mitochondrial genomes of six newly sequenced diatoms reveals group II introns in the barcoding region of cox1. Mitochondrial DNA Part A, DNA mapping. Seq Anal. 2019;30:43–51.
Guillory WX, Onyshchenko A, Ruck EC, Parks M, Nakov T, Wickett NJ, et al. Recurrent Loss, Horizontal Transfer, and the Obscure Origins of Mitochondrial Introns in Diatoms (Bacillariophyta). Genome Biol Evol. 2018;10:1504–15.
Burger G, Saint-Louis D, Gray MW, Lang BF. Complete sequence of the mitochondrial DNA of the red alga Porphyra purpurea: Cyanobacterial introns and shared ancestry of red and green algae. Plant Cell. 1999;11:1675–94.
Odom OW, Shenkenberg DL, Garcia JA, Herrin DL. A horizontally acquired group II intron in the chloroplast psbA gene of a psychrophilic Chlamydomonas: In vitro self-splicing and genetic evidence for maturase activity. RNA. 2004;10:1097–107.
Sheveleva EV, Hallick RB. Recent horizontal intron transfer to a chloroplast genome. Nucleic Acids Res. 2004;32:803–10.
Ehara M, Watanabe KI, Ohama T. Distribution of cognates of group II introns detected in mitochondrial cox1 genes of a diatom and a haptophyte. Gene. 2000;256:157–67.
Kamikawa R, Masuda I, Demura M, Oyama K, Yoshimatsu S, Kawachi M, et al. Mitochondrial Group II Introns in the Raphidophycean Flagellate Chattonella spp. suggests a Diatom-to-Chattonella Lateral Group II Intron Transfer. Protist. 2009;160:364–75.
Vallès Y, Halanych KM, Boore JL. Group II introns break new boundaries: presence in a bilaterian’s genome. PLoS ONE. 2008;3:e1488.
Kwon E-C, Kim J-H, Kim N-S. Comprehensive genomic analyses with 115 plastomes from algae to seed plants: structure, gene contents, GC contents, and introns. Genes Genomics. 2020;42:553–70.
Lemieux C, Otis C, Turmel M. Comparative Chloroplast Genome Analyses of Streptophyte Green Algae Uncover Major Structural Alterations in the Klebsormidiophyceae, Coleochaetophyceae and Zygnematophyceae. Front Plant Sci. 2016;7:697.
Teich R, Grauvogel C, Petersen J. Intron distribution in Plantae: 500 million years of stasis during land plant evolution. Gene. 2007;394:96–104.
Bonen L. Cis- and trans-splicing of group II introns in plant mitochondria. Mitochondrion. 2008;8:26–34.
Orton LM, Fitzek E, Fitzek E, Feng X, Grayburn WS, Mower JP, et al. Zygnema circumcarinatum UTEX 1559 chloroplast and mitochondrial genomes provide insight into land plant evolution. J Exp Bot. 2020;71:3361–73.
Qiu YL, Cho Y, Cox JC, Palmer JD. The gain of three mitochondrial introns identifies liverworts as the earliest land plants. Nature. 1998;394:671–4.
Qiu Y-L, Li L, Wang B, Chen Z, Knoop V, Groth-Malonek M, et al. The deepest divergences in land plants inferred from phylogenomic evidence. Proc Natl Acad Sci U S A. 2006;103:15511–6.
Puttick MN, Morris JL, Williams TA, Cox CJ, Edwards D, Kenrick P, et al. The interrelationships of land plants and the nature of the ancestral embryophyte. Curr Biol. 2018;28:733–45.
Renzaglia KS, Villarreal Aguilar JC, Garbary DJ. Morphology supports the setaphyte hypothesis: mosses plus liverworts form a natural group. Bryophyt Divers Evol. 2018;40:11–7.
Morris JL, Puttick MN, Clark JW, Edwards D, Kenrick P, Pressel S, et al. The timescale of early land plant evolution. Proc Natl Acad Sci U S A. 2018;115:E2274–83.
Su D, Yang L, Shi X, Ma X, Zhou X, Hedges SB, et al. Large-Scale Phylogenomic Analyses Reveal the Monophyly of Bryophytes and Neoproterozoic Origin of Land Plants. Mol Biol Evol. 2021;38:3332–44.
Steinhauser S, Malek O, Schock I, Beckert S, Knoop V. Group II introns in plant mitochondria: Trans-splicing and transgenic approaches. In: Moller IM, Gardeström P, Glimelius K, Glaser E, editors. Leiden: Backhuys Publishers; 1998. pp. 201–6.
Groth-Malonek M, Wahrmund U, Polsakiewicz M, Knoop V. Evolution of a Pseudogene: Exclusive Survival of a Functional Mitochondrial nad7 Gene Supports Haplomitrium as the Earliest Liverwort Lineage and Proposes a Secondary Loss of RNA Editing in Marchantiidae. Mol Biol Evol. 2007;24:1068–74.
Nedelcu AM, Lee RW. A degenerate group II intron in the intronless mitochondrial genome of Chlamydomonas reinhardtii: Evolutionary implications. Mol Biol Evol. 1998;15:918–22.
Guo W, Mower JP. Evolution of plant mitochondrial intron-encoded maturases: frequent lineage-specific loss and recurrent intracellular transfer to the nucleus. J Mol Evol. 2013;77:43–54.
Turmel M, Otis C, Lemieux C. Mitochondrion-to-Chloroplast DNA. Transfers and Intragenomic Proliferation of Chloroplast Group II Introns in Gloeotilopsis Green Algae (Ulotrichales, Ulvophyceae). Genome Biol Evol. 2016;8:2789–805.
Manhart JR, Palmer JD. The gain of two chloroplast tRNA introns marks the green algal ancestors of land plants. Nature. 1990;345:268–70.
Brouard J-S, Turmel M, Otis C, Lemieux C. Proliferation of group II introns in the chloroplast genome of the green alga Oedocladium carolinianum (Chlorophyceae). PeerJ. 2016;4:e2627.
Cho Y, Palmer JD. Multiple acquisitions via horizontal transfer of a group I intron in the mitochondrial cox1 gene during evolution of the Araceae family. Mol Biol Evol. 1999;16:1155–65.
Cusimano N, Zhang LBL-B, Renner SS. Reevaluation of the cox1 group I intron in Araceae and angiosperms indicates a history dominated by loss rather than horizontal transfer. Mol Biol Evol. 2008;25:265–76.
Vaughn JC, Mason MT, Sper-Whitis GL, Kuhlman P, Palmer JD. Fungal origin by horizontal transfer of a plant mitochondrial group I intron in the chimeric CoxI gene of Peperomia. J Mol Evol. 1995;41:563–72.
Guo W, Zhu A, Fan W, Mower JP. Complete mitochondrial genomes from the ferns Ophioglossum californicum and Psilotum nudum are highly repetitive with the largest organellar introns. New Phytol. 2017;213:391–403.
Blocker FJH, Mohr G, Conlan LH, Qi L, Belfort M, Lambowitz AM. Domain structure and three-dimensional model of a group II intron-encoded reverse transcriptase. RNA. 2005;11:14–28.
San Filippo J, Lambowitz AM. Characterization of the C-Terminal DNA-binding/DNA Endonuclease Region of a Group II Intron-encoded Protein. J Mol Biol. 2002;324:933–51.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421.
Hadley Wickham R, François L, Henry, Kirill Müller. A Grammar of Data Manipulation [R package dplyr version 1.0.0]. Media. 2020.
Katoh K, Misawa K, Kuma KI, Miyata T. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30:3059–66.
Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016;33:1870–4.
Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol Biol Evol. 2020;37:1530–4.
Stöver BC, Müller KF. TreeGraph 2: combining and visualizing evidence from different phylogenetic analyses. BMC Bioinformatics. 2010;11 SUPPL. 6.
Sato K, Akiyama M, Sakakibara Y. RNA secondary structure prediction using deep learning with thermodynamic integration. Nat Commun. 2021;12:1–9.
Darty K, Denise A, Ponty Y. VARNA: Interactive drawing and editing of the RNA secondary structure. Bioinformatics. 2009;25:1974–5.
Larsson J, Gustafsson P, {eulerr}: Area-Proportional Euler and Venn Diagrams with Ellipses. Proceedings of {{International Workshop}} on {{Set Visualization}} and {{Reasoning}}. 2018;2116:84–91.
Wickman H, François R, Henry L, Muller K. dlpyr: A Grammar of Data Manipulation. CRAN Repository. 2021;:1–2.
Cohen O, Ashkenazy H, Belinky F, Huchon D, Pupko T. GLOOME: gain loss mapping engine. Bioinformatics. 2010;26:2914–5.
Marchetti F, Cainzos M, Shevtsov S, Córdoba JP, Sultan LD, Brennicke A, et al. Mitochondrial Pentatricopeptide Repeat Protein, EMB2794, Plays a Pivotal Role in NADH Dehydrogenase Subunit nad2 mRNA Maturation in Arabidopsis thaliana. Plant Cell Physiol. 2020;61:1080–94.

No competing interests reported.

Download PDF

Journal Publication

published 12 Mar, 2023

Read the published version in BMC Ecology and Evolution →

Editorial decision: Major revision
05 Dec, 2022
Reviews received at journal
02 Dec, 2022
Reviewers agreed at journal
16 Nov, 2022
Reviews received at journal
12 Oct, 2022
Reviewers agreed at journal
03 Oct, 2022
Reviewers invited by journal
29 Sep, 2022
Editor assigned by journal
29 Sep, 2022
Editor invited by journal
29 Sep, 2022
Submission checks completed at journal
29 Sep, 2022
First submitted to journal
20 Sep, 2022

You are reading this latest preprint version

Categorizing 161 plant mitochondrial group II introns into 29 families of related paralogues finds only limited links between intron mobility and intron-borne maturases

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Results

Family F01

Family F02

Family F03

Family F04

Family F05

Family F06

Family F07

Family F08

Family F09

Family F10

Family F11

Family F12

Family F13

Family F14

Family F15

Family F16

Family F17

Family F18

Family F19

Family F20

Family F21

Family F22

Family F23

Family F24

Family F25

Discussion

Material And Methods

Declarations

Consent of participation

Consent of publication

Availability of data and materials

Competing interests

Funding

Author contribution

Acknowledgements

References

Additional Declarations

Supplementary Files

Status:

Journal Publication

Version 1