Reverse U-to-C editing exceeds C-to-U RNA editing in some ferns – a monilophyte-wide comparison of chloroplast and mitochondrial RNA editing suggests independent evolution of the two processes in both organelles
© The Author(s). 2016
Received: 23 April 2016
Accepted: 9 June 2016
Published: 21 June 2016
RNA editing by C-to-U conversions is nearly omnipresent in land plant chloroplasts and mitochondria, where it mainly serves to reconstitute conserved codon identities in the organelle mRNAs. Reverse U-to-C RNA editing in contrast appears to be restricted to hornworts, some lycophytes, and ferns (monilophytes). A well-resolved monilophyte phylogeny has recently emerged and now allows to trace the side-by-side evolution of both types of pyrimidine exchange editing in the two endosymbiotic organelles.
Our study of RNA editing in four selected mitochondrial genes show a wide spectrum of divergent RNA editing frequencies including a dominance of U-to-C over the canonical C-to-U editing in some taxa like the order Schizaeales. We find that silent RNA editing leaving encoded amino acids unchanged is highly biased with more than ten-fold amounts of silent C-to-U over U-to-C edits. In full contrast to flowering plants, RNA editing frequencies are low in early-branching monilophyte lineages but increase in later emerging clades. Moreover, while editing rates in the two organelles are usually correlated, we observe uncoupled evolution of editing frequencies in fern mitochondria and chloroplasts. Most mitochondrial RNA editing sites are shared between the recently emerging fern orders whereas chloroplast editing sites are mostly clade-specific. Finally, we observe that chloroplast RNA editing appears to be completely absent in horsetails (Equisetales), the sister clade of all other monilophytes.
C-to-U and U-to-C RNA editing in fern chloroplasts and mitochondria follow disinct evolutionary pathways that are surprisingly different from what has previously been found in flowering plants. The results call for careful differentiation of the two types of RNA editing in the two endosymbiotic organelles in comparative evolutionary studies.
RNA editing that converts specific cytidines into uridines in chloroplast and mitochondrial transcripts is nearly omnipresent among land plants [1–4]. As an apparently unique exception among plants, the marchantiid (complex-thalloid) liverworts have secondarily lost RNA editing [5, 6]. In most cases, RNA editing affects organelle mRNAs, where evolutionarily conserved codons are restored. However, numerous examples have also been found for RNA editing in tRNAs [7–10], UTRs , rRNAs , and introns [13–16]. Editing in the non-coding RNAs is likewise considered to be essential for the correct biological function of the respective molecules by re-establishing necessary base-pairing in secondary or tertiary RNA structures [15, 16].
Whereas the C-to-U type of RNA editing is nearly ubiquitous among land plants, the reverse process, U-to-C editing, appears to be more restricted in occurrence. U-to-C editing is unequivocally present in hornworts [17–19] and in monilophyte organelles [15, 20, 21]. U-to-C editing occurs in at least some lycophytes, the quillworts (Isoetales) and the club mosses (Lycopodiales) [10, 22–24]. It is surprisingly absent in Selaginellales, the third lycophyte order, despite having record numbers of C-to-U editing in both mitochondria and chloroplasts [12, 16].
The monilophytes are a morphologically heterogenous group that include the true eusporangiate ferns (sporangia with multicellular walls), such as Ophioglossales (moonworts) and Marattiales, the Equisetales (horsetails), the Psilotales (whisk ferns), and the species-rich group of leptosporangiate ferns (sporangia with unicellular walls). Early studies with rbcL, atpB, rps4, and nuclear 18S rDNA demonstrated the monophyly of the monilophytes and placed them as the sister group to the spermatophytes (seed plants) [25, 26]. However, the early dichotomies in the phylogeny of the monilophytes affecting the eusporangiate lineages and the horsetails relative to the leptosporangiate ferns remained unresolved. A recent study combining chloroplast loci atpA, atpB, rbcL, rps4, and matK with the mitochondrial loci nad2, nad5, atp1, and rpl2 confidently placed the horsetails as sister group to all other monilophytes and Ophioglossales/Psilotales as the sister group to a joint clade of Marattiales and leptosporangiate ferns .
Previous studies on mitochondrial and chloroplast RNA editing in ferns already indicated highly differing frequencies of RNA editing in different fern taxa. In the chloroplasts of Adiantum capillus-veneris  and Ophioglossum californicum  RNA editing is abundant (315/35 and 297/3 C-to-U/U-to-C editing sites), whereas in Psilotum nudum only 27 C-to-U and no U-to-C editing sites were detected . In the absence of a complete monilophyte mitochondrial genome, all mitochondrial RNA editing analyses in ferns are restricted to individual loci [15, 20, 21, 30].
The now available backbone phylogeny of monilophytes allows for phylogenetic insights into the evolution of C-to-U and U-to-C RNA editing among ferns, thus complementing similar studies restricted to C-to-U editing in flowering plants. To this end, we analyzed RNA editing in a phylogenetically wide sampling of monilophytes for four selected mitochondrial loci and for all available chloroplast genomic data. In contrast to an overall loss of RNA editing in angiosperms, we observed (i) an increase of RNA editing in monilophytes after the diversification of later emerging lineages, (ii) a mostly uncoupled evolution of editing frequencies in chloroplasts and mitochondria, and (iii) a largely uncoupled evolution of C-to-U and U-to-C editing. In some taxa, such as the order Schizaeales, U-to-C editing even exceeds the canonical C-to-U editing. Additionally, a complete chloroplast transcriptome analysis of the horsetail Equisetum hyemale confirmed our assumptions of total absence of RNA editing.
Mitochondrial RNA editing in monilophytes
Overall, the four different mitochondrial loci reflect comparable frequencies of RNA editing for each individual taxon (Fig. 1). However, the genus Gleichenia and Anemia phyllitidis are exceptions with unusual high numbers of RNA editing. Gleichenia has higher editing than all other examined taxa in two loci, the 1104 bp amplicon of the nad5 gene (92 C-to-U and 9 U-to-C edits) and the 411 bp amplicon of the rps1 gene (12 C-to-U edits). A. phyllitidis shows highest editing numbers in the other two loci, the rpl2 gene (16 C-to-U and 46 U-to-C edits) and the atp1 gene (53 C-to-U and 70 U-to-C edits).
Of the remaining 117 non-silent edit sites, 115 were initially correctly predicted by PREPACT with other monilophyte cDNAs as references. One predicted candidate editing event (atp1eU134TM) remained unconfirmed in the A. phyllitidis cDNA. This false positive candidate site was also predicted for Lygodium japonicum but likewise remained unconfirmed. Accordingly, the overall number of 116 predicted events was very close to the actual number of verified 117 non-silent edit sites. Therefore, also editing predictions for other monilophyte taxa could be considered as trustworthy (Fig. 1).
Contrary to the high levels of RNA editing in Anemia or Gleichenia, we observed very low numbers of editing in the early-branching eusporangiate fern lineages and the horsetails (Equisetales). In contrast to Anemia and Gleichenia that each feature more than 200 editing sites in the four sampled mitochondrial loci, only 27 editing sites were found in the same four genes of Angiopteris evecta.
Other than overall high RNA editing numbers, G. dicarpa also shows widely divergent ratios of C-to-U vs. U-to-C editing ranging from 12:0 in rps1 to 19:21 in rpl2. The latter example is the earliest emerging case of U-to-C editing exceeding C-to-U editing in our phylogenetic sampling (Fig. 1). The trend of U-to-C editing surpassing the canonical C-to-U editing becomes more pronounced in the later branching lineages, reaching a peak with an overall dominance of U-to-C over C-to-U editing at a ratio of 159:93 in A. phyllitidis (Schizaeales). While fewer U-to-C RNA editing sites were predicted (and confirmed) for the eusporangiate lineages Psilotales, Ophioglossales, and the Osmundales as the earliest-branching leptosporangiate lineage, we failed to find any evidence for U-to-C editing in Equisetales and Marattiales.
The ratio of C-to-U and U-to-C RNA editing becomes more balanced in the later emerging water ferns (Salviniales), tree ferns (Cyatheales) and in the large order of Polypodiales, which contains most of the extant fern diversity . Resolving the phylogeny of Cyatheales and Salviniales relative to the Polypodiales has been demanding, likely owing to the extremely divergent generation times resulting in very short branch lengths of the tree ferns  as opposed to the long branches of the water ferns. Noteworthy are the very similar frequencies of RNA editing of both types in these two morphologically and developmentally extremely divergent fern groups (Fig. 1).
Our extended data set of 1794 mitochondrial editing sites identified by cDNA analysis (Fig. 1) revealed a strong bias of silent editing when the two types of pyrimidine conversions are compared. Out of altogether 960 documented events of C-to-U editing, 114 are silent (11.9 %). In contrast, only 7 of 834 identified sites of U-to-C editing (0.8 %) are silent.
Comprehensive analysis of the Equisetum hyemale chloroplast transcriptome
In the course of our studies we observed very low numbers of predicted chloroplast RNA editing in the Equisetum hyemale chloroplast genome (accession number KC117177) . No in-frame stop codons are present in the protein coding genes, which are a good indicator for reverse editing. However, a conserved putative start codon in the accD gene is ACG, which may be subject to C-to-U editing. In order to investigate the issue, we performed an exhaustive cDNA analysis of the Equisetum hyemale chloroplast coding regions. All annotated introns are spliced out correctly and all negative controls did not show any traces of a potential contamination of genomic DNA. However, we did not find a single RNA editing site in any protein coding gene. Even the creation of the putative AUG start codon by C-to-U editing in accD could not verified. We therefore conclude that the Equisetum hyemale chloroplast is completely devoid of RNA editing.
Comparison of mitochondrial and chloroplastid RNA editing
For the comparison of mitochondrial and chloroplast RNA editing frequencies we investigated all available complete monilophyte chloroplast genomes. Extensive chloroplast editome analyses have previously been done for A. capillus-veneris , O. californicum and P. nudum  and are valuable additional references for predictions of RNA editing in the chloroplast genomes of Alsophila spinulosa , Diplopterygium glaucum, Osmundastrum cinnamomeum , Angiopteris evecta , Marsilea crenata and Lygodium japonicum .
We sampled the chloroplast genes ndhF, chlB and chlL to test for the low numbers of predicted RNA editing in Angiopteris evecta for comparison to the low levels of mitochondrial editing we identified (Fig. 1). Indeed, we identified only seven predicted C-to-U editing sites (chlBeU1268SL, chlLeU317SL, ndhFeU296SL, ndhFeU431PL, ndhFeU980PL, ndhFeU1001PL and ndhFeU1028SL) in the respective amplicons. Of these, only chlLeU317SL had previously been reported as a monilophyte editing site, in Adiantum capillus-veneris . The chlB edit in Angiopteris is exclusively shared with the lycophyte Selaginella uncinata  and edit ndhFeU1001PL is shared with the basal angiosperm Amborella  and the hornwort Anthoceros angustus . In contrast to the confirmed C-to-U editing sites, none of three predicted U-to-C edits (chlBeC923VA, ndhFeC826FL, ndhFeC847FL) were confirmed.
Conservation of editing sites in mitochondria and chloroplasts
This discrepancy in the evolution of mitochondrial vs. chloroplast editing patterns is still confirmed when taxa of the heavy-editing Schizaeales (Figs. 4c, d) and Gleicheniales (Figs. 4e, f) are included. Although Schizaeales and Gleicheniales feature large numbers of unique mitochondrial editing sites, a high number of edits (64 and 34, respectively) are shared between all four or five orders (Figs. 4c and e). In contrast, of more than 1,000 chloroplast editing sites only 8 are shared among the representatives of all five clades, respectively.
Phylogenetic studies of plant organelle RNA editing have so far focused on flowering plants for several obvious reasons [40–43]. First, angiosperms contain most plant model taxa like Arabidopsis, rice, pea, wheat or tobacco where RNA editing has been studied extensively. Second, much more organelle genome information is available for comparative studies among flowering plants than for other plant clades. Finally, a very a good phylogenetic framework has been available for angiosperms since a couple of years.
Monilophytes, the sister group to seed plants (spermatophytes), are particularly interesting for comparative studies of organelle RNA editing evolution for two main reasons. First, the evolution of editing spans a much larger time frame since the monilophyte clade is about three times as old as the angiosperms. Second, the plant-typical C-to-U-type of RNA editing is accompanied by U-to-C editing in the reverse direction in monilophytes. Moreover, a well-supported backbone phylogeny of ferns has been obtained very recently , now offering the opportunity to trace the evolution of chloroplast and mitochondrial RNA editing of both types. We here present such a side-by-side study on organelle RNA editing in monilophytes paying particular attention to consider chloroplast vs. mitochondrial editing on the one hand and C-to-U vs. U-to-C editing on the other hand separately.
Several interesting observations emerge. Our analysis shows a high variance of RNA editing in monilophytes. There is a clear trend towards more editing in the more diversified and species-rich orders of the leptosporangiate ferns compared to the basal eusporangiate fern lineages. This is in full contrast to the large-scale evolution of RNA editing in angiosperms, where RNA editing is frequent in early lineages but significantly lower in later emerging lineages, both in mitochondria [44–46] and in chloroplasts .
Despite extreme differences in morphology and development, the three extant lineages of crown leptosporangiates, the Polypodiales, the Cyatheales (tree ferns) and the Salviniales (water ferns) show comparable mitochondrial RNA editing frequencies of ca. 30 editing sites per kb both for C-to-U and U-to-C editing (Fig. 3). Peak mitochondrial editing frequencies, however, occur in earlier-branching leptosporangiate orders with C-to-U editing reaching a top value at 40 sites per kb in Gleichenia and U-to-C editing with more than 45 events per kb in Anemia.
Particularly noteworthy is the finding that in Anemia phyllitidis more than two thirds of all editing sites are of the U-to-C type. This strong dominance of reverse over the canonical C-to-U editing has to our knowledge never been observed before in a plant organelle. Intriguingly, we find that such a dominance of reverse U-to-C editing is apparently restricted to the mitochondrial lineage. Even in the case of the Schizaeales, the fraction of reverse editing in the chloroplast amounts to only 30 % of all candidate sites. Certainly, the different genera analyzed in our estimates of chloroplast vs. mitochondrial RNA editing are a remaining caveat.
Comparable numbers of co-existing U-to-C and C-to-U editing have previously been found (or predicted) in hornworts [17–19]. Hornworts (Anthocerotophyta) obviously are another highly interesting plant clade for comparative studies of the different types of editing in the future. However, the current scarcity of organelle sequence data (in particular plastome sequences) and the low numbers of extant hornwort taxa impede such analyses.
Other than the highly variable numbers of C-to-U and U-to-C editing in the two endosymbiotic organelles across monilophyte diversity, the entirely different patterns of conservation of individual editing sites are striking. Only 3.5 % of chloroplast RNA editing sites are shared between the representatives of Salviniales, Cyatheales and Polypodiales (Fig. 4b). In contrast a full 45 % of mitochondrial RNA editing sites are shared between the representatives of the three orders (Fig. 4a), indicating shared ancestry of those sites from their last common ancestor living approximately 218 million years ago. Similar pictures emerge when representatives of the Schizaeales and Gleicheniales are serially included for comparison of editing site conservations in the ancestors living approximately 266 or 282 million years ago, respectively. However, the numbers of shared sites also decline very quickly for the mitochondrial comparisons, very likely owing to dramatic changes in the organelle editomes with numerous gains and losses in that period of monilophyte evolution. Both the individual RNA editing sites as well as their total numbers become more conserved in the mitochondrial lineage with the LCA of Cyatheales, Polypodiales and Salviniales. No similar observation can be made for the chloroplast lineage, however, indicating a highly dynamic gain and loss scenario in the chloroplast editomes. Again, as a caveat, denser sampling of taxa in both organelles will be needed in the future to evaluate these insights.
Gains and losses of organelle RNA editing sites are expected to be accompanied by corresponding changes in the nuclear-encoded co-factors addressing individual RNA editing sites in chloroplast or mitochondria. Key editing factors are unique RNA-binding pentatricopeptide repeat (PPR) proteins of the so-called “PLS-type”, which are particularly abundant in the plant lineage [47–49]. The dynamic gain and loss phases of organelle editing sites in monilophytes and in particular in the leptosporangiate ferns are expected to be accompanied by a massive gain, loss, or re-assignments of the nuclear encoded PLS-type PPR proteins. Unfortunately, ferns are notorious for their large, polyploid and complex genomes [50–52]. Hitherto available monilophyte transcriptome data, e.g. in the OneKP project  are as yet of insufficient quality and preclude to give good estimates on the complexity and diversity of PPR gene families or to make RNA targeting prognoses for individual members based on the recently deduced PPR-RNA binding code [54, 55]. Accordingly, no monilophyte sample has as yet been included in a very recent novel approach to comprehensively identify and distinguish PPR proteins in available genomic data . This may hopefully change in the near future. A recent genome project sequencing the nuclear genome of Azolla filiculoides is ongoing and the nuclear genome of Ceratopteris richardii is proposed as a possible project .
Other than their tremendous diversity of editing patterns and the simultaneous existence of reverse U-to-C editing, ferns offer yet another feature distinguishing them from flowering plants with regard to organelle RNA editing. RNA editing in angiosperms is affected by additional protein factors called multiple organelle RNA editing factors (MORFs) or RNA editing factor interacting proteins (RIPs) [58–61], which assemble together with PPR proteins as complex editosomes on the target RNAs. No evidence for MORFs/RIPs is identified in the hitherto available fern genome or transcriptome data, very much like in the more ancient bryophyte or lycophyte lineages. Hence, the monilophyte clade may have retained a more ancient and simple RNA editing machinery devoid of helper proteins like MORFs/RIPs.
Most importantly, any forthcoming fern genome data will help to elucidate the yet enigmatic biochemical machinery of reverse U-to-C editing. The DYW domain with cytidine deaminase similarity present at the carboxyterminus of many PLS-type RNA editing factors is the bona fide candidate to perform the deamination of cytidine to uridine . No convincing protein candidate has a yet been proposed for the reverse reaction evidently requiring an amino group donor as a co-substrate for the amination of uridine. Any proteins proposed to be involved in U-to-C editing should be correspondingly diverse in taxa for which we here demonstrate a high proportion in this type of editing. Of particular interest in this respect is the much lower amount of silent editing for U-to-C conversions (0.8 %) in comparison to the phylogenetically more widely distributed C-to-U editing (11.9 %). A similar bias has previously emerged for silent editing in mitochondria of Isoetes engelmannii . Whether these findings may indicate a higher fidelity of sequence recognition specificity of the yet elusive U-to-C editing factors remains to be seen.
The recent editing analysis of the chloroplast genomes of Ophioglossum californicum and Psilotum nudum assumed that reverse editing is entirely lost in Psilotum and that it may also be lost in Equisetum, Angiopteris and Osmundastrum based on the absence of in-frame stop codons in chloroplastid genes . Here, we found that RNA editing in both organelles evolves independently and show experimental evidence for reverse editing in Psilotum nudum mitochondria (Fig. 1). The hypothesis of a complete loss of reverse editing in Angiopteris and Equisetum, and possibly in the orders Marattiales and Equisetales altogether, is supported by our study, however. Similarly, no evidence for reverse RNA editing had previously been found in the study of other mitochondrial genes in Equisetum .
Moreover, we could not find any evidence for RNA editing at all in our extensive transcriptome study of the chloroplast genome of Equisetum hyemale. To our knowledge this would be the first example of RNA editing documented in only one of the two endosymbiotic organelles. Extrapolating from the only 12 mitochondrial editing sites identified in the three amplicons investigated here (Fig. 1) one can assume that, depending on its total mitochondrial gene complement, Equisetum hyemale may have only around 100 mitochondrial editing sites. Even in the case of the model moss Physcomitrella patens with only 11 mitochondrial editing sites [63, 64], two editing sites exist in the chloroplast [65, 66]. Similarly reduced RNA editing is observed for sister taxa in other Funariaceae mosses, too .
The overall low amount of RNA editing in the early-branching eusporangiate lineages, the complete absence of chloroplast editing in Equisetum hyemale and the likely absence of reverse editing altogether in Marattiales and Equisetales may indicate that the last common ancestor of monilophytes had very low amounts of RNA editing in general and reverse U-to-C editing in particular. However, we are very reluctant to come to this conclusion given that members of the lycophytes, the sister lineages to euphyllophytes (comprising monilophytes and spermatophytes), feature highly frequent and highly diverse RNA editing in both organelles, exclusively of the C-to-U type in Selaginellales [12, 16] and in both direction of pyrimidine exchange in the Isoetales . Since fern mitochondrial RNA editing proved to be more deeply conserved we checked for its conservation among the lycophytes (Additional file 2: Figure S1). This comparison shows that the overwhelming majority of editing sites in the high-level editing lycophytes Isoetes and Selaginella are taxon-specific, notably also in comparison to the low-editing, early-branching lycophyte Phlegmariurus (Additional file 2: Figure S1B), making independent gains of most editing sites in Isoetes and Selaginella more likely.
One way or the other, the monilophytes have experienced dramatic changes in their editomes. Our results warrant for differentiated considerations of plant organelle RNA editing, both with respect to mitochondria vs. chloroplasts and with respect to C-to-U vs. U-to-C exchanges, for which we here show independent patterns of evolution. The most intriguing open question concerns the evolutionary forces driving the massive increases of C-to-U but notably also U-to-C RNA editing during the diversification of the early leptosporangiate lineages.
This study reports a broad-scale evolutionary comparative analysis of C-to-U and U-to-C RNA editing in the two endosymbiotic organelles of monilophytes. We find that mitochondrial RNA editing is highly diverse in monilophytes, including particularly low editing rates in early vs. late branching lineages, which is in full contrast to the evolutionary patterns previously observed among flowering plants. The reverse type of U-to-C RNA editing appears to be completely absent in Equisetales (horsetails) and in Marattiales. Within the leptosporangiate ferns, however, RNA editing of both types is highly abundant with record amounts of reverse editing in Schizaeales. Mitochondrial RNA editing sites in the leptosporangiate ferns are strikingly conserved in contrast to RNA editing in the chloroplasts. Hence, C-to-U and U-to-C RNA editing is evolving independently in the two organelles and the results call for careful differentiation between the different types of RNA editing in mitochondria and chloroplasts. Its great variability in the fern organelles promises an interesting field of co-evolution of C-to-U and U-to-C editing sites and their hitherto unknown nuclear-encoded specificity factors. Importantly, the biochemical mechanisms of “reverse” U-to-C editing, evolutionarily much more restricted than the more widespread C-to-U editing both within and outside of the plant kingdom , are completely unknown at present. Whereas cytidine deaminases are key to the latter editing type, a yet completely enigmatic transamination process involving an unknown amino-group donor must be postulated for the former. We believe that the here documented diversity of editing in both directions among monilophytes will ultimately help to identify the U-to-C editing factors and explain their unique features such as the here identified low amount of accompanying “superfluous” silent editing possibly indicating a higher degree of target specificity than the factors performing C-to-U editing.
Plant material and molecular work
Plant material was obtained from the Botanic Garden Bonn. Nucleic acids were isolated according to the CTAB method  followed by DNA digestion with DNAse I (Thermo Scientific/Fermentas). Total RNA was reverse transcribed into single-stranded cDNA with the RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific/Fermentas). Priming was performed with random hexamer primers (Roth). Subsequent RT-PCRs with Go-Taq (Promega) and/or Q5 Polymerase (New England Biolabs) were done with primers from [21, 27, 31]. For the transcriptome analysis of the Equisetum hyemale chloroplast primers binding in the UTRs of single or polycistronic genes were used. As a negative control samples without reverse transcriptase treatment were used to test for remnants of genomic DNA.
RT-PCR products were separated by agarose gel electrophoresis and recovered by the NucleoSpin Extract II Kit (Macherey Nagel). Gel-eluted products were either sequenced directly or after cloning into the pGEM-T Easy vector (Promega) and amplification in Escherichia coli. Sequencing was done by Macrogen Europe (Amsterdam, NL). New cDNA sequences were submitted to GenBank (see Additional file 1: Table S1). Sequences were analyzed with MEGA 5.05  and aligned using the implemented ClustalW algorithm.
RNA editing analysis
For each of the four investigated mitochondrial loci (atp1, nad5, rpl2, rps1), cDNA sequences were determined for at least 10 fern taxa (Fig. 1, Additional file 1: Table S1). The cDNA sequences were employed as references to predict mitochondrial RNA editing sites in the respective other monilophyte taxa for which we did not generate cDNA sequences using the alignment prediction tool of PREPACT . Sites were counted as candidate editing sites when predicted by at least 80 % of the references.
Chloroplast RNA editing sites were predicted with the BLASTX prediction tool of PREPACT. Again, a threshold for prediction of 80 % and ten selected reference sequences were used: Adiantum capillus-veneris, Amborella trichopoda, Arabidopsis thaliana, Chara vulgaris, Equisetum hyemale, Marchantia polymorpha, Ophioglossum californicum, Pellia endiviifolia, Physcomitrella patens and Psilotum nudum. All editing sites were labeled according to a previously published nomenclature proposal [64, 71] using the name of the gene, type of editing, the nucleotide position and the respective amino acid change.
The authors are very grateful to Dr. Wolfram Lobin, curator of the Bonn University Botanical Garden, and his team for providing plant material.
The study relied exclusively on basic funding through the University of Bonn only.
Availability of data and materials
All sequences generated in this study have been deposited in GenBank. For a detailed overview see Additional file 1: Table S1.
FG, NK and SF conducted the experimental work. NK and VK analyzed the data, designed the study and wrote the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Chateigner-Boutin A, Small I. Organellar RNA editing. Wiley Interdiscip Rev RNA. 2011;2:493–506. doi:10.1002/wrna.72.View ArticlePubMedGoogle Scholar
- Knoop V. When you can’t trust the DNA: RNA editing changes transcript sequences. Cell Mol Life Sci. 2011;68:567–86. doi:10.1007/s00018-010-0538-9.View ArticlePubMedGoogle Scholar
- Steinhauser S, Beckert S, Capesius I, Malek O, Knoop V. Plant mitochondrial RNA editing. J Mol Evol. 1999;48:303–12.View ArticlePubMedGoogle Scholar
- Freyer R, Kiefer-Meyer MC, Kössel H. Occurrence of plastid RNA editing in all major lineages of land plants. Proc Natl Acad Sci U S A. 1997;94:6285–90.View ArticlePubMedPubMed CentralGoogle Scholar
- 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. doi:10.1093/molbev/msm026.View ArticlePubMedGoogle Scholar
- Rüdinger M, Volkmar U, Lenz H, Groth-Malonek M, Knoop V. Nuclear DYW-Type PPR Gene Families Diversify with Increasing RNA Editing Frequencies in Liverwort and Moss Mitochondria. J Mol Evol. 2012;74:37–51. doi:10.1007/s00239-012-9486-3.View ArticlePubMedGoogle Scholar
- Binder S, Marchfelder A, Brennicke A. RNA editing of tRNA(Phe) and tRNA(Cys) in mitochondria of Oenothera berteriana is initiated in precursor molecules. Mol Gen Genet. 1994;244:67–74.View ArticlePubMedGoogle Scholar
- Marchfelder A, Brennicke A, Binder S. RNA editing is required for efficient excision of tRNA(Phe) from precursors in plant mitochondria. J Biol Chem. 1996;271:1898–903.View ArticlePubMedGoogle Scholar
- Grewe F, Viehoever P, Weisshaar B, Knoop V. A trans-splicing group I intron and tRNA-hyperediting in the mitochondrial genome of the lycophyte Isoetes engelmannii. Nucleic Acids Res. 2009;37:5093–104. doi:10.1093/nar/gkp532.View ArticlePubMedPubMed CentralGoogle Scholar
- Grewe F, Herres S, Viehöver P, Polsakiewicz M, Weisshaar B, Knoop V. A unique transcriptome: 1782 positions of RNA editing alter 1406 codon identities in mitochondrial mRNAs of the lycophyte Isoetes engelmannii. Nucleic Acids Res. 2011;39:2890–902. doi:10.1093/nar/gkq1227.View ArticlePubMedGoogle Scholar
- Schuster W, Hiesel R, Wissinger B, Brennicke A. RNA editing in the cytochrome b locus of the higher plant Oenothera berteriana includes a U-to-C transition. Mol Cell Biol. 1990;10:2428–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Hecht J, Grewe F, Knoop V. Extreme RNA Editing in Coding Islands and Abundant Microsatellites in Repeat Sequences of Selaginella moellendorffii Mitochondria: The Root of Frequent Plant mtDNA Recombination in Early Tracheophytes. Genome Biol Evol. 2011;3:344–58. doi:10.1093/gbe/evr027.View ArticlePubMedGoogle Scholar
- Castandet B, Choury D, Bégu D, Jordana X, Araya A. Intron RNA editing is essential for splicing in plant mitochondria. Nucleic Acids Res. 2010;38:7112–21. doi:10.1093/nar/gkq591.View ArticlePubMedPubMed CentralGoogle Scholar
- Carrillo C, Bonen L. RNA editing status of nad7 intron domains in wheat mitochondria. Nucleic Acids Res. 1997;25:403–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Bégu D, Castandet B, Araya A. RNA editing restores critical domains of a group I intron in fern mitochondria. Curr Genet. 2011;57:317–25. doi:10.1007/s00294-011-0349-z.View ArticlePubMedGoogle Scholar
- Oldenkott B, Yamaguchi K, Tsuji-Tsukinoki S, Knie N, Knoop V. Chloroplast RNA editing going extreme: more than 3400 events of C-to-U editing in the chloroplast transcriptome of the lycophyte Selaginella uncinata. RNA. 2014;20:1499–506. doi:10.1261/rna.045575.114.View ArticlePubMedPubMed CentralGoogle Scholar
- Xue J, Liu Y, Li L, Wang B, Qiu Y. The complete mitochondrial genome sequence of the hornwort Phaeoceros laevis: retention of many ancient pseudogenes and conservative evolution of mitochondrial genomes in hornworts. Curr Genet. 2010;56:53–61. doi:10.1007/s00294-009-0279-1.View ArticlePubMedGoogle Scholar
- Li L, Wang B, Liu Y, Qiu Y. The complete mitochondrial genome sequence of the hornwort Megaceros aenigmaticus shows a mixed mode of conservative yet dynamic evolution in early land plant mitochondrial genomes. J Mol Evol. 2009;68:665–78. doi:10.1007/s00239-009-9240-7.View ArticlePubMedGoogle Scholar
- Kugita M, Yamamoto Y, Fujikawa T, Matsumoto T, Yoshinaga K. RNA editing in hornwort chloroplasts makes more than half the genes functional. Nucleic Acids Res. 2003;31:2417–23. doi:10.1093/nar/gkg327.View ArticlePubMedPubMed CentralGoogle Scholar
- Bonavita S, Regina TMR. The evolutionary conservation of rps3 introns and rps19-rps3-rpl16 gene cluster in Adiantum capillus-veneris mitochondria. Curr Genet. 2015. doi:10.1007/s00294-015-0512-z.PubMedGoogle Scholar
- 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. doi:10.1111/j.1438-8677.1999.tb00249.x.View ArticleGoogle Scholar
- Wolf PG, Karol KG, Mandoli DF, Kuehl J, Arumuganathan K, Ellis MW, et al. The first complete chloroplast genome sequence of a lycophyte, Huperzia lucidula (Lycopodiaceae). Gene. 2005;350:117–28. doi:10.1016/j.gene.2005.01.018.View ArticlePubMedGoogle Scholar
- Karol KG, Arumuganathan K, Boore JL, Duffy AM, Everett KDE, Hall JD, et al. Complete plastome sequences of Equisetum arvense and Isoetes flaccida: implications for phylogeny and plastid genome evolution of early land plant lineages. BMC Evol Biol. 2010;10:321. doi:10.1186/1471-2148-10-321.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu Y, Wang B, Cui P, Li L, Xue J, Yu J, et al. The mitochondrial genome of the lycophyte Huperzia squarrosa: the most archaic form in vascular plants. PLoS ONE. 2012;7, e35168. doi:10.1371/journal.pone.0035168.View ArticlePubMedPubMed CentralGoogle Scholar
- Pryer KM, Schuettpelz E, Wolf PG, Schneider H, Smith AR, Cranfill R. Phylogeny and evolution of ferns (monilophytes) with a focus on the early leptosporangiate divergences. Am J Bot. 2004;91:1582–98. doi:10.3732/ajb.91.10.1582.View ArticlePubMedGoogle Scholar
- Pryer KM, Schneider H, Smith AR, Cranfill R, Wolf PG, Hunt JS, Sipes SD. Horsetails and ferns are a monophyletic group and the closest living relatives to seed plants. Nature. 2001;409:618–22. doi:10.1038/35054555.View ArticlePubMedGoogle Scholar
- Knie N, Fischer S, Grewe F, Polsakiewicz M, Knoop V. Horsetails are the sister group to all other monilophytes and Marattiales are sister to leptosporangiate ferns. Mol Phylogenet Evol. 2015;90:140–9. doi:10.1016/j.ympev.2015.05.008.View ArticlePubMedGoogle Scholar
- Wolf PG, Rowe CA, Hasebe M. High levels of RNA editing in a vascular plant chloroplast genome: analysis of transcripts from the fern Adiantum capillus-veneris. Gene. 2004;339:89–97. doi:10.1016/j.gene.2004.06.018.View ArticlePubMedGoogle Scholar
- Guo W, Grewe F, Mower JP. Variable frequency of plastid RNA editing among ferns and repeated loss of uridine-to-cytidine editing from vascular plants. PLoS ONE. 2015;10, e0117075. doi:10.1371/journal.pone.0117075.View ArticlePubMedPubMed CentralGoogle Scholar
- Bégu D, Araya A. The horsetail Equisetum arvense 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. doi:10.1007/s00294-008-0225-7.View ArticlePubMedGoogle Scholar
- Knie N, Grewe F, Knoop V. Monilophyte mitochondrial rps1 genes carry a unique group II intron that likely originated from an ancient paralogue in rpl2. 2016. RNA. In press.Google Scholar
- Wikström N, Pryer KM. Incongruence between primary sequence data and the distribution of a mitochondrial atp1 group II intron among ferns and horsetails. Mol Phylogenet Evol. 2005;36:484–93. doi:10.1016/j.ympev.2005.04.008.View ArticlePubMedGoogle Scholar
- Smith AR, Pryer KM, Schuettpelz E, Korall P, Schneider H, Wolf PG. A Classification for Extant Ferns. Taxon. 2006;55:705–31. doi:10.2307/25065646.View ArticleGoogle Scholar
- Zhong B, Fong R, Collins LJ, McLenachan PA, Penny D. Two new fern chloroplasts and decelerated evolution linked to the long generation time in tree ferns. Genome Biol Evol. 2014;6:1166–73. doi:10.1093/gbe/evu087.View ArticlePubMedPubMed CentralGoogle Scholar
- Grewe F, Guo W, Gubbels EA, Hansen AK, Mower JP. Complete plastid genomes from Ophioglossum californicum, Psilotum nudum, and Equisetum hyemale reveal an ancestral land plant genome structure and resolve the position of Equisetales among monilophytes. BMC Evol Biol. 2013;13:8. doi:10.1186/1471-2148-13-8.View ArticlePubMedPubMed CentralGoogle Scholar
- Gao L, Yi X, Yang Y, Su Y, Wang T. Complete chloroplast genome sequence of a tree fern Alsophila spinulosa: insights into evolutionary changes in fern chloroplast genomes. BMC Evol Biol. 2009;9:130. doi:10.1186/1471-2148-9-130.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim HT, Chung MG, Kim K. Chloroplast genome evolution in early diverged leptosporangiate ferns. Mol Cells. 2014;37:372–82. doi:10.14348/molcells.2014.2296.View ArticlePubMedPubMed CentralGoogle Scholar
- Roper JM, Kellon Hansen S, Wolf PG, Karol KG, Mandoli DF, Everett KDE, et al. The Complete Plastid Genome Sequence of Angiopteris evecta (G. Forst.) Hoffm. (Marattiaceae). Am Fern J. 2007;97:95–106. doi:10.1640/0002-8444(2007)97[95:TCPGSO]2.0.CO;2.View ArticleGoogle Scholar
- Gao L, Wang B, Wang Z, Zhou Y, Su Y, Wang T. Plastome sequences of Lygodium japonicum and Marsilea crenata reveal the genome organization transformation from basal ferns to core leptosporangiates. Genome Biol Evol. 2013;5:1403–7. doi:10.1093/gbe/evt099.View ArticlePubMedPubMed CentralGoogle Scholar
- Hein A, Polsakiewicz M, Knoop V. Frequent chloroplast RNA editing in early-branching flowering plants: pilot studies on angiosperm-wide coexistence of editing sites and their nuclear specificity factors. BMC Evol Biol. 2016;16:23. doi:10.1186/s12862-016-0589-0.View ArticlePubMedPubMed CentralGoogle Scholar
- Mower JP. Modeling sites of RNA editing as a fifth nucleotide state reveals progressive loss of edited sites from angiosperm mitochondria. Mol Biol Evol. 2008;25:52–61. doi:10.1093/molbev/msm226.View ArticlePubMedGoogle Scholar
- Tillich M, Funk HT, Schmitz-Linneweber C, Poltnigg P, Sabater B, Martin M, et al. Editing of plastid RNA in Arabidopsis thaliana ecotypes. Plant J. 2005;43:708–15. doi:10.1111/j.1365-313X.2005.02484.x.View ArticlePubMedGoogle Scholar
- Mulligan RM, Chang KLC, Chou CC. Computational analysis of RNA editing sites in plant mitochondrial genomes reveals similar information content and a sporadic distribution of editing sites. Mol Biol Evol. 2007;24:1971–81. doi:10.1093/molbev/msm125.View ArticlePubMedGoogle Scholar
- Rice DW, Alverson AJ, Richardson AO, Young GJ, Sanchez-Puerta MV, Munzinger J, et al. Horizontal transfer of entire genomes via mitochondrial fusion in the angiosperm Amborella. Science. 2013;342:1468–73. doi:10.1126/science.1246275.View ArticlePubMedGoogle Scholar
- Chaw S, Shih AC, Wang D, Wu Y, Liu S, Chou T. The mitochondrial genome of the gymnosperm Cycas taitungensis contains a novel family of short interspersed elements, Bpu sequences, and abundant RNA editing sites. Mol Biol Evol. 2008;25:603–15. doi:10.1093/molbev/msn009.View ArticlePubMedGoogle Scholar
- Giegé P, Brennicke A. RNA editing in Arabidopsis mitochondria effects 441 C to U changes in ORFs. Proc Natl Acad Sci U S A. 1999;96:15324–9.View ArticlePubMedPubMed CentralGoogle Scholar
- O’Toole N, Hattori M, Andres C, Iida K, Lurin C, Schmitz-Linneweber C, et al. On the expansion of the pentatricopeptide repeat gene family in plants. Mol Biol Evol. 2008;25:1120–8. doi:10.1093/molbev/msn057.View ArticlePubMedGoogle Scholar
- Fujii S, Small I. The evolution of RNA editing and pentatricopeptide repeat genes. New Phytol. 2011;191:37–47. doi:10.1111/j.1469-8137.2011.03746.x.View ArticlePubMedGoogle Scholar
- Schallenberg-Rüdinger M, Knoop V. Coevolution of Organelle RNA Editing and Nuclear Specificity Factors in Early Land Plants. In: Genomes and Evolution of Charophytes, Bryophytes, Lycophytes and Ferns: Elsevier; 2016. p. 37–93. doi:10.1016/bs.abr.2016.01.002.
- Wolf PG, Sessa EB, Marchant DB, Li F, Rothfels CJ, Sigel EM, et al. An Exploration into Fern Genome Space. Genome Biol Evol. 2015;7:2533–44. doi:10.1093/gbe/evv163.View ArticlePubMedPubMed CentralGoogle Scholar
- Dyer RJ, Pellicer J, Savolainen V, Leitch IJ, Schneider H. Genome size expansion and the relationship between nuclear DNA content and spore size in the Asplenium monanthes fern complex (Aspleniaceae). BMC Plant Biol. 2013;13:219. doi:10.1186/1471-2229-13-219.View ArticlePubMedPubMed CentralGoogle Scholar
- Obermayer R, Leitch IJ, Hanson L, Bennett MD. Nuclear DNA C-values in 30 species double the familial representation in pteridophytes. Ann Bot. 2002;90:209–17.View ArticlePubMedPubMed CentralGoogle Scholar
- Matasci N, Hung L, Yan Z, Carpenter EJ, Wickett NJ, Mirarab S, et al. Data access for the 1,000 Plants (1KP) project. Gigascience. 2014;3:17. doi:10.1186/2047-217X-3-17.View ArticlePubMedPubMed CentralGoogle Scholar
- Yagi Y, Hayashi S, Kobayashi K, Hirayama T, Nakamura T. Elucidation of the RNA recognition code for pentatricopeptide repeat proteins involved in organelle RNA editing in plants. PLoS ONE. 2013;8, e57286. doi:10.1371/journal.pone.0057286.View ArticlePubMedPubMed CentralGoogle Scholar
- Barkan A, Rojas M, Fujii S, Yap A, Chong YS, Bond CS, Small I. A combinatorial amino acid code for RNA recognition by pentatricopeptide repeat proteins. PLoS Genet. 2012;8, e1002910. doi:10.1371/journal.pgen.1002910.View ArticlePubMedPubMed CentralGoogle Scholar
- Cheng S, Gutmann B, Zhong X, Ye Y, Fisher MF, Bai F, et al. Redefining the structural motifs that determine RNA binding and RNA editing by pentatricopeptide repeat proteins in land plants. Plant J. 2016;85:532–47. doi:10.1111/tpj.13121.View ArticlePubMedGoogle Scholar
- Sessa EB, Banks JA, Barker MS, Der JP, Duffy AM, Graham SW, et al. Between two fern genomes. Gigascience. 2014;3:15. doi:10.1186/2047-217X-3-15.View ArticlePubMedPubMed CentralGoogle Scholar
- Shikanai T. RNA editing in plants: Machinery and flexibility of site recognition. Biochim Biophys Acta. 1847;2015:779–85. doi:10.1016/j.bbabio.2014.12.010.Google Scholar
- Bentolila S, Heller WP, Sun T, Babina AM, Friso G, van Wijk KJ, Hanson MR. RIP1, a member of an Arabidopsis protein family, interacts with the protein RARE1 and broadly affects RNA editing. Proc Natl Acad Sci U S A. 2012;109:E1453–61. doi:10.1073/pnas.1121465109.View ArticlePubMedPubMed CentralGoogle Scholar
- Takenaka M, Verbitskiy D, Zehrmann A, Härtel B, Bayer-Császár E, Glass F, Brennicke A. RNA editing in plant mitochondria -Connecting RNA target sequences and acting proteins. Mitochondrion. 2014. doi:10.1016/j.mito.2014.04.005.Google Scholar
- Zehrmann A, Härtel B, Glass F, Bayer-Császár E, Obata T, Meyer E, et al. Selective homo- and heteromer interactions between the multiple organellar RNA editing factor (MORF) proteins in Arabidopsis thaliana. J Biol Chem. 2015;290:6445–56. doi:10.1074/jbc.M114.602086.View ArticlePubMedPubMed CentralGoogle Scholar
- Salone V, Rüdinger M, Polsakiewicz M, Hoffmann B, Groth-Malonek M, Szurek B, et al. A hypothesis on the identification of the editing enzyme in plant organelles. FEBS Lett. 2007;581:4132–8. doi:10.1016/j.febslet.2007.07.075.View ArticlePubMedGoogle Scholar
- Tasaki E, Hattori M, Sugita M. The moss pentatricopeptide repeat protein with a DYW domain is responsible for RNA editing of mitochondrial ccmFc transcript. Plant J. 2010;62:560–70. doi:10.1111/j.1365-313X.2010.04175.x.View ArticlePubMedGoogle Scholar
- Rüdinger M, Funk HT, Rensing SA, Maier UG, Knoop V. RNA editing: only eleven sites are present in the Physcomitrella patens mitochondrial transcriptome and a universal nomenclature proposal. Mol Genet Genomics. 2009;281:473–81. doi:10.1007/s00438-009-0424-z.View ArticlePubMedGoogle Scholar
- Miyata Y, Sugiura C, Kobayashi Y, Hagiwara M, Sugita M. Chloroplast ribosomal S14 protein transcript is edited to create a translation initiation codon in the moss Physcomitrella patens. Biochim Biophys Acta. 2002;1576:346–9.View ArticlePubMedGoogle Scholar
- Miyata Y, Sugita M. Tissue- and stage-specific RNA editing of rps 14 transcripts in moss (Physcomitrella patens) chloroplasts. J Plant Physiol. 2004;161:113–5. doi:10.1078/0176-1617-01220.View ArticlePubMedGoogle Scholar
- Rüdinger M, Szövényi P, Rensing SA, Knoop V. Assigning DYW-type PPR proteins to RNA editing sites in the funariid mosses Physcomitrella patens and Funaria hygrometrica. Plant J. 2011;67:370–80. doi:10.1111/j.1365-313X.2011.04600.x.View ArticlePubMedGoogle Scholar
- Doyle LL, Doyle JL. Isolation of plant DNA from fresh tissue. Focus. 1990;12:13–5.Google Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9. doi:10.1093/molbev/msr121.View ArticlePubMedPubMed CentralGoogle Scholar
- Lenz H, Knoop V. PREPACT 2.0: Predicting C-to-U and U-to-C RNA Editing in Organelle Genome Sequences with Multiple References and Curated RNA Editing Annotation. Bioinform Biol Insights. 2013;7:1–19. doi:10.4137/BBI.S11059.View ArticlePubMedPubMed CentralGoogle Scholar
- Lenz H, Rüdinger M, Volkmar U, Fischer S, Herres S, Grewe F, Knoop V. Introducing the plant RNA editing prediction and analysis computer tool PREPACT and an update on RNA editing site nomenclature. Curr Genet. 2010;56:189–201. doi:10.1007/s00294-009-0283-5.View ArticlePubMedGoogle Scholar
- Schneider H, Schuettpelz E, Pryer KM, Cranfill R, Magallon S, Lupia R. Ferns diversified in the shadow of angiosperms. Nature. 2004;428:553–7. doi:10.1038/nature02361.View ArticlePubMedGoogle Scholar