- Research Article
- Open Access
Plant organelle RNA editing and its specificity factors: enhancements of analyses and new database features in PREPACT 3.0
BMC Bioinformatics volume 19, Article number: 255 (2018)
Gene expression in plant chloroplasts and mitochondria is affected by RNA editing. Numerous C-to-U conversions, accompanied by reverse U-to-C exchanges in some plant clades, alter the genetic information encoded in the organelle genomes. Predicting and analyzing RNA editing, which ranges from only few sites in some species to thousands in other taxa, is bioinformatically demanding.
Here, we present major enhancements and extensions of PREPACT, a WWW-based service for analysing, predicting and cataloguing plant-type RNA editing. New features in PREPACT’s core include direct GenBank accession query input and options to restrict searches to candidate U-to-C editing or to sites where editing has been documented previously in the references. The reference database has been extended by 20 new organelle editomes. PREPACT 3.0 features new modules “EdiFacts” and “TargetScan”. EdiFacts integrates information on pentatricopeptide repeat (PPR) proteins characterized as site-specific RNA editing factors. PREPACT’s editome references connect into EdiFacts, linking editing events to specific co-factors where known. TargetScan allows position-weighted querying for sequence motifs in the organelle references, optionally restricted to coding regions or sequences around editing sites, or in queries uploaded by the user. TargetScan is mainly intended to evaluate and further refine the proposed PPR-RNA recognition code but may be handy for other tasks as well. We present an analysis for the immediate sequence environment of more than 15,000 documented editing sites finding strong and different bias in the editome data sets.
We exemplarily present the novel features of PREPACT 3.0 aimed to enhance the analyses of plant-type RNA editing, including its new modules EdiFacts integrating information on characterized editing factors and TargetScan aimed to analyse RNA editing site recognition specificities.
Nearly 30 years after the discovery of C-to-U RNA editing in plant mitochondria [1,2,3] and quickly thereafter also in chloroplasts , the field has recently expanded tremendously in several directions of research [5,6,7]. After the initial characterization of a first chloroplast  and a first mitochondrial RNA editing factor  numerous such proteins continue to be characterized, quickly outdating published compilations [5, 10,11,12] by ever more new additions [13,14,15]. The key factors in RNA editing site recognition are pentatricopeptide repeat (PPR) proteins, which are encoded by tremendously enlarged gene families with hundreds of members in plants [16,17,18,19].
The arrays of PPRs are key to specifically recognizing the RNA sequences upstream of cytidines targeted for conversion into uridines via deamination. PPR proteins serving as editing factors have a unique makeup of alternating P-, L- and S-type PPRs with distinct amino acid conservation profiles. Moreover, PPR proteins characterized as editing factors carry carboxyterminal protein domain additions, minimally “E” (extension) domains, frequently followed by the so-called “DYW” domain. The latter in particular is of fundamental interest owing to its significant similarity to cytidine deaminases, which likely provides the biochemical activity for C-to-U conversion [20,21,22,23].
Intriguingly, DYW-type PPR proteins that were previously believed to be plant-specific, have recently been identified in very distant evolutionary lineages of eukaryotes where their presence likewise seems to be connected to mitochondrial RNA editing of the C-to-U type [24,25,26,27].
A PPR-RNA recognition code has been proposed , which is currently subject to further amendments and experimental testing in vivo and in vitro [29,30,31,32]. Linking RNA editing events or other transcript targets to specific PPR array sequences and vice versa is becoming an exciting field for bioinformatic approaches and for potential future applications using artificially designed PPR arrays [31, 33, 34]. The former issue becomes obvious, for example, when numbers of editing events both in mitochondria and in chloroplasts literally run into thousands, such as in the lycophytes [35,36,37].
The PREPACT WWW service developed in our group [38, 39] aimed for (i) standardizing RNA editing annotation and nomenclature, (ii) making the vast and ever-increasing amount of editing information available with manually curated reference editomes (i.e. the sets of editing sites determined with extensive cDNA analysis for organelle genomes), and (iii) helping to analyze and predict RNA editing in organelle sequence data. We here demonstrate an update of PREPACT in version 3.0 with respect to its “classic” features, but now also aiming to address the interplay between RNA editing sites and their cognate PPR-type specificity factors. Information on the latter are now included in a novel database module “EdiFacts” and the possibility to experimentally scan for potential RNA targets is realized with the new “TargetScan” module. We present the new features of PREPACT’s core functionalities, 20 new editome reference addendums, demonstrate the functionalities of EdiFacts and TargetScan and discuss future issues and developments of RNA editing analysis, especially those related to PPR-RNA recognition.
PREPACT editome reference extensions
One core component of PREPACT’s functionality is a set of mitochondrial and chloroplast genomes with curated and standardized RNA editing site information . A user-defined selection of these reference editomes can be used to simultaneously identify organelle protein-coding genes and candidate RNA editing sites in an unannotated organelle nucleotide sequence query using PREPACT’s BLASTX mode (Fig. 1).
We have added several new organelle editome references with reliably determined editing site identifications (Table 1). The 14 chloroplast editomes of the flowering plants Amborella trichopoda , Aegilops tauschii , coconut Cocos nucifera , cucumber Cucumis sativus [40, 43], cotton Gossypium hirsutum , the orchid Phalaenopsis aphrodite , the duckweed Spirodela polyrhiza  and the mung bean Vigna radiata , the gymnosperm Ginkgo biloba , the liverwort Apopellia endiviifolia , the lycophyte Selaginella uncinata , the horsetail Equisetum hyemale  and the ferns Ophioglossum californicum and Psilotum nudum  have been added to the plastome references previously included in PREPACT 2.0 . The gymnosperm Ginkgo, the early-branching angiosperm Amborella and the jungermanniid liverwort Apopellia fill important taxonomic gaps. Similarly, the horsetail Equisetum and the club moss Selaginella are important addendums as taxa representing the full range of a taxon lacking chloroplast editing altogether  and the most heavily edited organelle transcriptome known so far with more than 3400 sites of C-to-U editing .
Among the mitochondrial references, Ophioglossum californicum and Psilotum nudum are particularly valuable additions as the first fern mitochondrial editomes . The editomes of Liriodendron tulipifera , Amborella trichopoda  and Cocos nucifera  are interesting additions representing early diverging angiosperm and monocot lineages. Moreover, we have added the mitochondrial DNA of the protist Acrasis kona where two events of plant-type C-to-U editing have recently been identified  as a further mitochondrial editome reference.
In some cases we refrained from adding further reference data owing to an evident lack of documented editing sites in the editomes at this stage like in the case of the rubber tree Hevea brasiliensis mitochondria , the soy bean Glycine max  and the Utricularia reniformis chloroplast  as well as the Saccharum officinarum  and Cycas organelle editomes [60, 61]. Likewise, we avoided editome data, which obviously seem to be affected by artefacts including non-canonical types of editing, which we could not reproduce in independent cDNA analyses, like in the chloroplast transcriptome studies of Ipomoea batata , Deschampsia antarctica  or Elaeis guineensis . Altogether, the editome references now available in the updated PREPACT 3.0 database comprise 27 chloroplast and 25 mitochondrial entries.
PREPACT input enhancements
The enhanced query interface of PREPACT 3.0 has several new options (Fig. 1). Searches for reverse U-to-C editing, previously only implemented as an optional addition, are now offered as an individual option allowing to restrict searches to U-to-C editing sites exclusively. Moreover, users may choose to restrict searches for candidate editing sites to positions where RNA editing at an orthologous position has previously been identified in at least one of the chosen references. A further option allows to always include such sites in the commons output even when below the overall threshold settings. The sequence input has been redesigned for dynamic handling of queries, now also allowing to simply enter database accession numbers, which are directly retrieved from GenBank/NCBI. Certainly, uploading or copy-pasting of FASTA-formatted data remains possible, too. Sequences are now checked on-the-fly to report formatting errors and to allow for immediate re-upload of corrected data. For the “Alignment” and “cDNA” analysis modes of PREPACT, multiple uploaded sequences can be sorted, deleted and rearranged between query and reference side using drag-and-drop, as we will detail below (see Additional file 1).
PREPACT output enhancements
Examples for the PREPACT “commons” table output summarizing the RNA editing events predicted from comparisons to the selected references are shown in Fig. 2. PREPACT applies our previously suggested RNA editing site nomenclature, which is composed of the affected gene followed by an ‘e’ for editing, the nucleotide introduced by editing (C or U), the nucleotide position in the coding sequence and the resulting codon identity (if applicable) before and after editing to label editing sites . The BLASTX output comprises the editing predictions from the selected references as individual tabs and a summary prediction as the final commons tab (Fig. 2a). The commons tab output now also displays amino acid identities in references that do not contribute to editing site prognoses either because of retention of the unedited state or due to an inconvertible codon at the corresponding position. This new feature helps significantly in the interpretation of the output because it immediately shows the variability of amino acids present for a candidate site predicted by only some of the references. For example, in the case of predicting editing event petLeU5PL converting a proline into a leucine (Fig. 2a), a chemically similar isoleucine (I) is present in the alga Chaetosphaeridium and in the liverwort Apopellia, which can be taken to further corroborate the likelihood of editing in the query (here Wollemia nobilis KP259800). Similarly, in the case of a leucine codon in rps18 that is widely conserved in plants and algae and which requires editing from a serine codon in several references, a phenylalanine codon (F) is present in the Poales Oryza and Zea (Fig. 2b). Removal of the polar amino acid serine may be more important than the presence of either an aliphatic leucine or an aromatic phenylalanine at this position in the protein. A hyphen in the commons output is now restricted to cases where homology is lacking in a reference, e. g. in the case of the ndh genes lost altogether in the orchid Phalaenopsis aphrodite plastome (Fig. 2c).
Amending editome data
With an enlarged data set of editome references several cases became evident where editing sites may have been missed in previous analyses or where unexpected “orphan” editing events reported previously are restricted to individual taxa. The enhanced output now displaying non-edited or non-editable codons combined with the red highlighting of known editing events in the references facilitates interpretation of the results in the commons tab. For example, the presence of editing site ndhHeU505HY in phylogenetically distant taxa including Amborella, Anthoceros, Cocos and Hevea (red) and conservation of a genomically encoded tyrosine in other taxa (black; Fig. 2c) strongly suggested that the editing event was missed in the early Cucumis transcriptome study . We recently checked upon such cases extensively in Cucumis confirming this and several other candidate sites to extend its chloroplast editome . Here, we took the opportunity to selectively also investigate other cases, such as rps2eU134TI in Atropa belladonna (Fig. 2d) by independent cDNA analyses and indeed found that many such sites have apparently been overlooked in the previous transcriptome studies. Altogether we already confirmed 56 additional events of RNA editing in 10 species by our independent cDNA analyses (Additional file 2). The newly confirmed events of RNA editing now identified were incorporated into the updated PREPACT 3.0 references.
Less conserved RNA editing sites and shifts in amino acid conservation
In some cases, it becomes apparent that a shift in amino acid conservation has obviously affected RNA editing sites during plant evolution. The rps2 gene is a case in point, exemplarily shown for the Atropa belladonna rps2 query (Fig. 2d). Editing of a serine codon in position 248 is fundamental in several dicot angiosperms to convert it into a leucine codon conserved in all taxa. In contrast, editing rps2eU107SF in the hornwort Anthoceros appears to reflect an ancestral state to reconstitute a conserved phenylalanine (F) codon in algae, liverworts, mosses and ferns (here represented by Chara, Apopellia, Physcomitrella and Adiantum) but not in the angiosperms, which lack editing and retain the genomically encoded serine codon. Editing event rps2eU134TI converting a threonine into an isoleucine codon is among those now to be added to the Atropa chloroplast editome (and the ones of Oryza and Zea) upon our re-investigation of cDNAs (Additional file 2). At the same time, this site is a further example for an editing event where another, but chemically similar, amino acid – in this case leucine (L) – is present at the corresponding position in early-branching taxa (Fig. 2d).
A similar case is the hornwort RNA editing atpIeU158SL serving to reconstitute a leucine codon conserved as a plesiomorphy in the early-branching plant lineages whereas a serine codon remains unaltered in the flowering plants (Fig. 2e). These are interesting cases of differential conservation of RNA editing sites and amino acid signatures, possibly indicating functional protein adaptations during evolution.
In few cases, codon changes introduced by RNA editing in the chloroplast seem erratic. The case of the “orphan” editing psbZeU11AV introducing a valine codon exclusively in the Selaginella psbZ transcript (Fig. 2f) is an example where a genomically encoded alanine is present and retained in all other references. Cases like these, in which chemically similar amino acids are exchanged, may reflect tolerable mis-firings of the editing machinery in taxa like the lycopyhtes where editing is particularly abundant. More complex are the divergent editing patterns even among angiosperms alone, like in the case of ndhG (Fig. 2g). Whether such cases reflect less conservation of the individual protein subunits, and accordingly of editing, or may rather indicate adaptations to interaction partners in the protein complexes, here possibly associated with a loss of ndhG editing in the Solanaceae (Atropa and Nicotiana), remains to be seen. While edit ndhGeU332SL in Anthoceros (Fig. 2g) could be an orphan edit like psbZeU11AV in Selaginella (Fig. 2f), this edit is in fact shared with Selaginella and may rather reflect conservation of a leucine in the early-branching land plants (not shown). Examples like these emphasize the importance of taxonomically diverse editome data sets.
New features for multiple sequence comparisons
Aside from its BLASTX mode to identify coding regions and candidate RNA editing sites de novo in uncharacterized organelle sequence queries, PREPACT offers analyses of multiple sequence alignments for comparative analysis of RNA editing. The “cDNA” and the “alignment prediction” mode are intended for comparative analyses and graphic display for a set of homologous sequences including one or multiple references. We here demonstrate the new functionalities using aligned sequences of the small mitochondrial atp9 gene for a phylogenetically wide sampling as an example for the alignment prediction mode (Additional file 1). Uploading a multiple-sequence FASTA file now displays the names of all individual sequences, which can be re-sorted in order, dragged-and-dropped between the collection of references and entries for prediction or can be individually deleted (Additional file 1). When multiple references are used, the output is organized into separate tabs for the individual references plus the comparative commons tab (if selected), analogous to the BLASTX mode (Additional file 3). We have now introduced the “pie chart” mode for more informative graphic display with three sections of a circle distinguishing editings in the three different codon positions (Fig. 3). Silent editing events in 3rd codon positions (and 1st position leucine YUR codons) is of relevance for the cDNA analysis mode only (not shown). The atp9 example reflects several cases where simultaneous editing in codon positions 1 and 2 is needed to reconstitute conserved codons (codons 28, 53, 62 and 68) in the heavily editing lycophytes and gymnosperms. RNA editing frequencies vary significantly in each plant clade. The Physcomitrella patens RNA editing event atp9eU92SL affecting atp9 codon 31, which we will also discuss below in the context of the new TargetScan feature, is shared by all other mosses, the hornworts, lycopyhtes, three gymnosperms and the two angiosperms in our example (Fig. 3).
The EdiFacts module
As of writing of this manuscript, more than 70 nuclear-encoded RNA editing factors targeting specific chloroplast or mitochondrial RNA editing sites in plants have already been characterized. These site-specific RNA-editing factors are “PLS-type” PPR proteins with carboxyterminal E1, E2 and DYW domains. We have summarized the available information on the hitherto known editing factors as individual entries in “EdiFacts”, a database and query extension to PREPACT (Fig. 4). EdiFacts includes information about species, target genes and editing sites in the two organelles and links to the respective editing factor protein sequences and the corresponding literature reports. All data can be queried with Boolean AND/OR logic with options to choose multiple entries in fields where appropriate. Optional query restrictions can be made for the number of repeats in the PPR arrays, carboxyterminal protein domains or authors of the corresponding publications (Fig. 4a). This is exemplarily demonstrated with a query for “Physcomitrella patens” AND “ccmF” (Fig. 4a). This search retrieves the EdiFacts entries with ID 44 and 45 (Fig. 4b), corresponding to the characterization of Physcomitrella patens PPR proteins PPR_71 and PPR_65 [66, 67].
If a characterized factor is known for a given editing site, the PREPACT output of that editing event is now highlighted with italic font and underlining, and dynamically cross-linked to the respective EdiFacts entry as exemplarily demonstrated for the two Physcomitrella ccmFC editing sites (Fig. 4c). This allows users not only to identify candidate sites of editing where orthologous editing events have been seen in other taxa (red font), but immediately also reveals the information on co-factors when already known (italics and underlining). The two mitochondrial RNA editing sites, ccmFCeU103PS and ccmFCeU122SF, are widely conserved in the plant kingdom, but pre-edited as serine and phenylalanine codons in the (non-editing) alga Chara and the liverwort Marchantia in our sampling. Our ccmFC example (Fig. 4c) also illustrates further examples of phylogenetically restricted (ccmFCeU104/107TI in Amborella and Liriodendron) or “orphan” RNA editings ccb452eU925PS in Arabidopsis and ccmFCeU706PS and ccmFCeU707SF in Nicotiana.
The TargetScan module
A new module TargetScan has been added to allow position-weighted querying of the PREPACT organelle reference data sets or of user-uploaded sequences for sequence motif matches (Fig. 5). The TargetScan interface allows users to define an oligonucleotide sequence motif with individual weighting of base preferences for A, C, G and T(U) using integers that automatically add up to 100 (%). Accordingly, a weighting of 25–25–25-25 reflects no nucleotide selectivity (equals “N” ambiguity) whereas 0–40–0-60 would, for example, reflect a strong selectivity for pyrimidines with a slight preference of T over C. Single input weights can be locked by clicking onto the respective nucleotide (switches from green to red background), distributing the remaining percentage evenly across the non-locked variants (Fig. 5b). Matrix input may be saved by download and re-uploaded. Users may select any combination of PREPACT references and/or other uploaded data for querying, hence allowing to scan for sequence targets across arbitrarily selected organelle references. Additional options allow to restrict the search for sequence targets within coding sequences or to regions around known editing sites (only in annotated references). While TargetScan may be helpful for diverse other issues, the latter options are mainly intended to identify and rank candidate targets of PPR-type editing factors.
We here exemplarily demonstrate the use of TargetScan for the Physcomitrella patens editing factor PPR_98 (Fig. 5a), which has been characterized as the specificity factor binding to the target sequence upstream of editing site atp9eU92SL (see Fig. 3) in vivo and in vitro [67, 68]. PPR_98 is a canonical (PLS)7-type PPR protein with a terminal DYW domain.
At the core of the PPR-RNA binding code  are position 5 selecting purines vs. pyrimidines with amino acid residues T (or S) vs. N and position L (‘Last’) selecting the keto nucleotides G or U vs. the amino nucleotides A or C with amino acids D vs. N in the PPRs (Fig. 5a). The first selection mechanism for purines vs. pyrimidines appears to be stronger and so is the distinction between the two purines as compared to the two pyrimidines. Moreover, the suggested code only fits the P- and S-type but not to the L-type PPRs, the functions of which remain to be explored. Accordingly, we arbitrarily weighted 90 vs. 10 for purine and 70 vs. 30 for pyrimidine selection in the canonical T/S + D, T/S + N and N + N/S or N + D-carrying P- and S-type motifs with pyrimidine recognition weighted as 100% and purine recognitions weighted as 200%. Additional weight was given to the position immediately upstream (− 1) of the cytidine editing target, here set arbitrarily to 15, 35, 5 and 45% for A, C, G and T (Fig. 5b) based on empirical observations that purines, and especially guanosines, occur only rarely upstream of an editing site (see also the new investigations outlined in the following chapter). All other positions were weighted 0% (Fig. 5b). In the case of PPR_98 and its target, nine PPRs fit perfectly to the above concept, whereas the binding code would suggest other nucleotide preferences for one P- and two S-type PPRs (Fig. 5a).
With the simple arbitrary weights outlined above, the correct RNA editing target of PPR_98 is identified with a top score of 1385 (calculated from the sum of all positions with their individual scores multiplied by their weight) set apart from the second-best hit (the ccmFCeU103PS editing site scoring 945) when scanning upstream of all Physcomitrella chloroplast and mitochondrial RNA editing sites (Fig. 6a). Such a ranking of matches among known editing sites is certainly helpful to identify best candidate targets of uncharacterized PPR proteins when editing sites are known. However, any amendments to the PPR-RNA recognition mechanisms must always explain why any similar transcribed sequences (at least in the same organelle) are not targeted. Extending the above search to all Physcomitrella mitochondrial coding sequences (CDS) with the above settings (plus an added arbitrary 200% for C at the candidate editing position to place potential cytidine targets top of the list) still identifies the atp9eU92SL editing site as the top-ranking target (Fig. 6b). However, the next-best candidate cytidines for editing in orf533, rps4, rpl10 and rps4 receive nearly as good scorings for upstream binding using the simple scoring scheme outlined above. Such cases warrant for (i) re-inspection for RNA editing at such candidate alternative target sites, (ii) testing for different binding weights of the different PPR motifs and target positions and (iii) subsequent testing with recombinant proteins in model systems like Physcomitrella and Arabidopsis. Evidently, any TargetScan query motif has a different a priori probability to identify top-scoring hits in different organelle genomes depending on nucleotide composition (mainly GC content). However, any results with identical scores identified in different genomes should indicate equally good matches, e.g. as likely targets for a given PPR protein like in the example discussed above.
Using TargetScan to explore the immediate editing site environments
The growing amount of complete and reliable organelle editome data now also allows to identify potential nucleotide bias in the immediate sequence environment of editing sites. These RNA positions are currently not assumed to be targets for recognition by PPRs, but nucleotide preferences could result from close interactions with the downstream E1, E2 and DYW domains of editing factors. We here investigated the immediate sequence vicinity of editing sites in positions − 4 to + 3 for nucleotide bias in six different organelle editome data sets now included in PREPACT: angiosperm plastomes, angiosperm chondromes, the Selaginella uncinata plastome, the Selaginella moellendorffii chondrome and the available chloroplast editomes of bryophytes and ferns, here considering C-to-U and U-to-C editing events separately (Fig. 7).
As an alternative to more demanding script programming approaches, users may employ simple weight matrices in TargetScan for such and similar purposes. Using arbitrary weightings for the four different nucleotides in each position (e.g. A:40, C:30, G:20, T:10), one can quickly obtain sorted results for each nucleotide identity as exemplarily shown for position − 1 in the angiosperm plastome example (Additional file 4). Compiling the data for positions − 4 to + 3 reveals striking differences for the six different editome sets concerning nucleotide bias in the positions surrounding the editing sites (Fig. 7).
It had been noted early that guanosines occur only very rarely upstream of edits. We now find that the bias against G is most pronounced in angiosperm chloroplast editomes (Fig. 7a) with only 0.3% G (and only 7.8% A) vs. a strong preference for U (69.9%) in position − 1. A similar picture strongly avoiding purines directly upstream of edits also emerges for the angiosperm chondromes (Fig. 7b) with only 2.3% G and 4.6% A. Likewise, a similar bias is also found in the bryophyte and fern chloroplast editomes, but significantly less so for the U-to-C editing events co-existing with C-to-U editing in these taxa (Fig. 7e and f). Surprisingly, for the reverse U-to-C editing sites in the hornwort and fern plastomes (Fig. 7f), another strong nucleotide bias in position + 2 behind the editing site (53.5% A vs. 9.6% C) even outnumbers the less pronounced bias in position − 1. While still biased against G, cytidines even dominate over uridines in position − 1 in the GC-rich organelle genomes of Selaginella (Fig. 7c-d).
Notably, our survey identifies yet further strong bias in other positions. Most prominent is the strong bias against C (7.1%) in position − 4 and the dominance of A (67.5%) vs. C (6.4%) in position + 1 for the angiosperm chloroplast editomes (Fig. 7a). These preferences are less pronounced, but recognizable, in the bryophyte and fern plastomes (Fig. 7e and f), but not in the angiosperm chondromes (Fig. 7b), possibly indicating slight differences between the editing machineries in the two endosymbiotic organelles.
Moreover, the U vs. G bias in position − 1 is now also becoming evident for position − 2 both in the chloroplast and in the mitochondrial editomes of angiosperms (Fig. 7a and b) and likewise also for C-to-U editing in bryophyte plastomes (Fig. 7e). In position − 3, however, only bryophytes and ferns, but not angiosperms, show a clear uridine over guanosine bias. Again, the latter findings cannot be generalized for the exceptionally GC-rich organelle genomes of Selaginella (Fig. 7c and d).
The functional extensions of PREPACT presented here serve different purposes for the analyses of “plant-type” C-to-U and U-to-C RNA editing in organelles. The upgrade of PREPACT’s repository of organelle editomes extends its functionality in the prognosis and comparative analysis of RNA editing. The success of RNA editing prognoses can be expected to increase with a wider sampling of references, especially of taxa more closely related to the query sequence taxon. Denser taxon sampling will help to identify both overlooked and orphan cases of RNA editing in individual taxa, as we have here shown exemplarily. Given the likely ever-growing set of editome information, a future update of PREPACT aims to include an option to select sets of references based on higher taxonomic ranks. In the future, we also hope to include additional complementary information from the literature for individual editing sites such as e.g. conflicting reports, strain variabilities or variable editing frequencies depending on development. Especially regarding a better understanding of the PPR-RNA recognition code discussed below, it is particularly important that editome analyses are correct and complete, i.e. without false positives and without actual editing events being missed. We found that some previous studies of organelle editomes based on early RNA-seq data require re-analyses. Meantime, suggestions for adequate design and improved analysis of such studies have been made [69, 70].
The implemented editome references will make PREPACT an expanding database of reported C-to-U and U-to-C RNA editing sites. Moreover, with the increasing number of discoveries of C-to-U RNA editing also outside of plants [24, 26], PREPACT will likely extend its scope in the future to also include yet more non-plant taxa. We hope to be able to provide timely updates of PREPACT and further reference editomes in the future.
Equally important, we have now integrated data on the increasing number of PPR-type editing factors that have been functionally characterized in the new database module EdiFacts. Accordingly, the user directly obtains information that an editing factor has been characterized in a reference taxon for a predicted editing position. This aspect is of considerable interest given the co-evolution of organelle editing sites and their nuclear-encoded specificity factors. A PPR-RNA binding code has been proposed, which is currently being refined and experimentally tested [28, 30, 31, 71, 72]. The proposed PPR-RNA recognition code offers an exciting field for bioinformatic and subsequent reverse genetic testing of proposed interactions between PPR arrays and RNA targets. However, other protein features outside of the immediate PPR-RNA interaction surface such as the immediate nucleotide environment of the editing sites  need to be taken into account, too. Future amendments, refinements and additional assumptions on how the RNA recognition factors recognize their target(s) need to be tested against the actual in vivo situation, i.e. all alternative transcriptome targets in a given organelle. We hope that the here described TargetScan module will be of help towards this issue. With respect to accompanying molecular studies, the moss Physcomitrella patens now occupies a unique position after completing the assignment of all its organelle RNA editing sites to specific cofactors [67, 68, 74]. Given its features as a genetic model organism it may be particularly attractive for transformation with mutated versions of RNA editing factors. For example, it has been shown that the terminal (and likely cytidine deaminase) DYW domains of different editing factors are exchangeable in some , but not in other cases , indicating at least a partial preference for their native targets, likely depending on the immediate sequence environment of the cytidine to be edited. Much further understanding is needed to adequately tune the (here arbitrarily selected) weights for those positions currently understood to contribute to sequence recognition (the P- and S-type PPRs) and ascribing proper weights to those elements (like the L-type PPRs, the E1, E2 and the DYW domains), for which a contribution to confer sequence specificity has yet to be elucidated.
The manual identification of potentially relevant amino acid positions from the loosely conserved P-, L- and S-type repeats of PPR proteins is tedious and cumbersome. The TPRpred tool served as a publicly available bioinformatic service for de novo identification of tetratricopeptide (TPR) and pentatricopeptide (PPR) repeats . A new WWW-service (initially under www.plantppr.com, now available at ppr.plantenergy.uwa.edu.au) allows to distinguish the P-, L- and S-type PPRs of plant PLS-type editing factors specifically after carefully reconsidering the domains of plant-type PPR proteins . An automatic extraction of the key residues like PPR positions 5 and L and their direct translation, e.g. via appropriate lookup-tables, for direct input into the new TargetScan of PREPACT is a future goal. Independent of this approach, the mutual assignment of editing sites and editing factors as implemented with EdiFacts will hopefully already now help to further explore the yet enigmatic co-evolution of organelle RNA editing sites and their nuclear co-factors.
Over the recent years, research on plant-type RNA editing has extended to the characterization of the specific, RNA-binding pentatricopeptide repeat (PPR) protein factors addressing individual editing sites in the endosymbiotic organelles. We have accordingly extended our WWW service PREPACT to include information on PPR-type editing factors in an additional database module EdiFacts. As a further extension of PREPACT, the new TargetScan module allows to search for position-weighted motifs in the now extended reference editome set of PREPACT or in user-defined references. The novel feature now implemented in version 3 of PREPACT should be of use to integrate information of RNA editing sites and corresponding specificity factors and help to further elucidate how PPR-type editing factors recognize their respective RNA targets.
The core functional implementation of PREPACT using PHP and MySQL has been described earlier [38, 39]. Basic functions have been revised to yield higher performance and to cope with growing query complexity. This included consistent translation of different sequence/feature numbering schemes on a global and local scale to be able to match information in partial hits and globally numbered features. The internal GenBank engine has been extended to also handle remote locations (in other accessions) and partial CDS features with annotated editing sites locally as well as in the remote part. This was necessary to also deal with complex genomes split across multiple accessions in parallel with trans-splicing as e.g. in the Amborella trichopoda mitochondrial DNA. The reference tabs of the BLASTX output (see Fig. 2) now offer an option for download of the individual references in a GenBank-style flat file format including the standardized annotation of RNA editing sites with the additional “RNA_editing” feature we had introduced previously .
EdiFacts is an addition to the relational database with data collected manually from publications. New items are continually identified by routine literature searches, journal publication alerts and journal scanning services such as “PubCrawler”  using appropriate key words. Literature references are downloaded, parsed and stored locally for search purposes and linked to respective external NCBI PubMed and protein source entries. Editing sites affected by listed factors are referenced in the “RNA_editing” feature introduced in PREPACT2  using a “db_xref” qualifier. This internal crosslink is used for highlighting editing sites with known editing factors in the “commons” output. The EdiFacts input form is the graphical representation of the internal query builder which translates various combinations of selected filters/options into efficient MySQL queries combining all available data.
The TargetScan module is comparing the user-defined weight matrix in a sliding window approach to the selected sequences or sequence parts extracted from the internal GenBank database. As such, TargetScan is a custom-made and easy-to use alternative to more sophisticated motif identification algorithms such as FIMO  or PWMscan . Scores for each sub-sequence are calculated by multiplying the base value (percent) with the position weight and summing up. Results are ranked by descending score down to a certain number of results or optionally to all results with the same score after this number of results to avoid arbitrary cut-offs of equally good matching sub-sequences. In the output individual base stretches are listed with their position/features according to the selected mode and single base scores are colour coded from green (maximum score at this position = perfectly matching) to red (minimum score at this position), with mixed colours in between. Positions with no weight are excluded from colour coding to have less clutter. Editing sites are highlighted in the sequence in blue (C-to-U) or red (U-to-C) respectively. To be in line with other sequence features, the selection of sub-sequences for searching in different modes (“Genome”, “CDS”, “Around editing sites”) is internally implemented as an extension to the GenBank format defining “Search_range” and “Search_result” as GenBank features.
For detection of previously overlooked RNA editing sites, individual chloroplast references were run against all other available reference editomes. Strongly predicted editing sites (i.e. with a ‘commons’ score of at least 80% or at least one edited reference species) previously not reported not to be edited were rechecked in selected cases (Additional file 2). To that end, plant material was obtained from the Bonn University Botanic Garden Bonn and RNA was prepared by the CTAB method, the TRI Reagent Protocol (Sigma Aldrich) or with the NucleoSpin® Plant RNA II Kit (Macherey-Nagel). Subsequently, cDNA synthesis was performed with Revert Aid First Strand cDNA Synthesis Kit (Thermo Fisher) using random hexamer primers. The relevant regions were amplified by RT-PCR with gene-specific primers and products recovered from agarose gel with NucleoSpin® Extract II Kit (Macherey-Nagel). PCR products were sequenced directly after gel elution or after cloning into pGEM-T Easy (Promega).
Cetyl-trimethyl ammonium bromide
Plant RNA editing prediction and analysis computer tool
Hiesel R, Wissinger B, Schuster W, Brennicke A. RNA editing in plant mitochondria. Science. 1989;246:1632–4.
Gualberto JM, Lamattina L, Bonnard G, Weil JH, Grienenberger JM. RNA editing in wheat mitochondria results in the conservation of protein sequences. Nature. 1989;341:660–2.
Covello PS, Gray MW. RNA editing in plant mitochondria. Nature. 1989;341:662–6.
Hoch B, Maier RM, Appel K, Igloi GL, Kössel H. Editing of a chloroplast mRNA by creation of an initiation codon. Nature. 1991;353:178–80.
Shikanai T. RNA editing in plants: machinery and flexibility of the site recognition. Biochim Biophys Acta Bioenerg. 2015;1874:779–85.
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;19(Pt B):191–7.
Takenaka M, Zehrmann A, Verbitskiy D, Härtel B, Brennicke A. RNA editing in plants and its evolution. Annu Rev Genet. 2013;47:335–52.
Kotera E, Tasaka M, Shikanai T. A pentatricopeptide repeat protein is essential for RNA editing in chloroplasts. Nature. 2005;433:326–30.
Zehrmann A, Verbitskiy D, van der Merwe JA, Brennicke A, Takenaka M. A DYW domain-containing pentatricopeptide repeat protein is required for RNA editing at multiple sites in mitochondria of Arabidopsis thaliana. Plant Cell. 2009;21:558–67.
Finster S, Legen J, Qu Y, Schmitz-Linneweber C. Land plant RNA editing or: Don’t be fooled by plant Organellar DNA sequences. In: Bock R, Knoop V, Dordrecht NL, editors. Genomics of chloroplasts and mitochondria. Netherlands: Springer; 2012. p. 293–321.
Knoop V. When you can’t trust the DNA: RNA editing changes transcript sequences. Cell Mol Life Sci. 2011;68:567–86.
Barkan A, Small I. Pentatricopeptide repeat proteins in plants. Annu Rev Plant Biol. 2014;65:415–42.
Yang Y-Z, Ding S, Wang H-C, Sun F, Huang W-L, Song S, Xu C, Tan B-C. The pentatricopeptide repeat protein EMP9 is required for mitochondrial ccmB and rps4 transcript editing, mitochondrial complex biogenesis and seed development in maize. New Phytol. 2017;214:782–95.
Qi W, Tian Z, Lu L, Chen X, Chen X, Zhang W, Song R. Editing of mitochondrial transcripts nad3 and cox2 by Dek10 is essential for mitochondrial function and maize plant development. Genetics. 2017;205:1489–501.
Wang G, Zhong M, Shuai B, Song J, Zhang J, Han L, Ling H, Tang Y, Wang G, Song R. E+ subgroup PPR protein defective kernel 36 is required for multiple mitochondrial transcripts editing and seed development in maize and Arabidopsis. New Phytol. 2017;214:1563–78.
Small I, Peeters N. The PPR motif - a TPR-related motif prevalent in plant organellar proteins. Trends Biochem Sci. 2000;25:46–7.
Andrés C, Lurin C, Small ID. The multifarious roles of PPR proteins in plant mitochondrial gene expression. Physiol Plant. 2007;129:14–22.
O’Toole N, Hattori M, Andres C, Iida K, Lurin C, Schmitz-Linneweber C, Sugita M, Small I. On the expansion of the Pentatricopeptide repeat gene family in plants. Mol Biol Evol. 2008;25:1120–8.
Lurin C, Andrés C, Aubourg S, Bellaoui M, Bitton F, Bruyère C, Caboche M, Debast C, Gualberto J, Hoffmann B, Lecharny A, Le Ret M, Martin-Magniette M-L, Mireau H, Peeters N, Renou J-P, Szurek B, Taconnat L, Small I. Genome-wide analysis of Arabidopsis pentatricopeptide repeat proteins reveals their essential role in organelle biogenesis. Plant Cell. 2004;16:2089–103.
Hayes ML, Giang K, Berhane B, Mulligan RM. Identification of two pentatricopeptide repeat genes required for RNA editing and zinc binding by C-terminal cytidine deaminase-like domains. J Biol Chem. 2013;288:36519–29.
Boussardon C, Avon A, Kindgren P, Bond CS, Challenor M, Lurin C, Small I. The cytidine deaminase signature HxE(x)nCxxC of DYW1 binds zinc and is necessary for RNA editing of ndhD-1. New Phytol. 2014;203:1090–5.
Iyer LM, Zhang D, Rogozin IB, Aravind L. Evolution of the deaminase fold and multiple origins of eukaryotic editing and mutagenic nucleic acid deaminases from bacterial toxin systems. Nucleic Acids Res. 2011;39:9473–97.
Salone V, Rüdinger M, Polsakiewicz M, Hoffmann B, Groth-Malonek M, Szurek B, Small I, Knoop V, Lurin C. A hypothesis on the identification of the editing enzyme in plant organelles. FEBS Lett. 2007;581:4132–8.
Schallenberg-Rüdinger M, Lenz H, Polsakiewicz M, Gott JM, Knoop V. A survey of PPR proteins identifies DYW domains like those of land plant RNA editing factors in diverse eukaryotes. RNA Biol. 2013;10:1549–56.
Rüdinger M, Fritz-Laylin L, Polsakiewicz M, Knoop V. Plant-type mitochondrial RNA editing in the protist Naegleria gruberi. RNA. 2011;17:2058–62.
Fu C-J, Sheikh S, Miao W, Andersson SGE, Baldauf SL. Missing genes, multiple ORFs and C-to-U type RNA editing in Acrasis kona (Heterolobosea, Excavata) mitochondrial DNA. Genome Biol Evol. 2014;6:2240–57.
Knoop V, Rüdinger M. DYW-type PPR proteins in a heterolobosean protist: plant RNA editing factors involved in an ancient horizontal gene transfer? FEBS Lett. 2010;584:4287–91.
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.
Nakamura T, Yagi Y, Kobayashi K. Mechanistic insight into pentatricopeptide repeat proteins as sequence-specific RNA-binding proteins for organellar RNAs in plants. Plant Cell Physiol. 2012;53:1171–9.
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.
Kindgren P, Yap A, Bond CS, Small I. Predictable alteration of sequence recognition by RNA editing factors from Arabidopsis. Plant Cell. 2015;27:403–16.
Okuda K, Shoki H, Arai M, Shikanai T, Small I, Nakamura T. Quantitative analysis of motifs contributing to the interaction between PLS-subfamily members and their target RNA sequences in plastid RNA editing. Plant J. 2014;80:870–82.
Yagi Y, Nakamura T, Small I. The potential for manipulating RNA with pentatricopeptide repeat proteins. Plant J. 2014;78:772–82.
Coquille S, Filipovska A, Chia T, Rajappa L, Lingford JP, Razif MFM, Thore S, Rackham O. An artificial PPR scaffold for programmable RNA recognition. Nat Commun. 2014;5:5729.
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.
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.
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.
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.
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.
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.
Wang M, Liu H, Ge L, Xing G, Wang M, Weining S, Nie X. Identification and analysis of RNA editing sites in the chloroplast transcripts of Aegilops tauschii L. Genes (Basel). 2017;8:13.
Huang Y-Y, Matzke AJM, Matzke M. Complete sequence and comparative analysis of the chloroplast genome of coconut palm (Cocos nucifera). PLoS One. 2013;8:e74736.
Guzowska-Nowowiejska M, Fiedorowicz E, Plader W. Cucumber, melon, pumpkin, and squash: are rules of editing in flowering plants chloroplast genes so well known indeed? Gene. 2009;434:1–8.
Jiang Y, Fan SL, Song MZ, Yu JN, Yu SX. Identification of RNA editing sites in cotton (Gossypium hirsutum) chloroplasts and editing events that affect secondary and three-dimensional protein structures. Genet Mol Res. 2012;11:987–1001.
Zeng WH, Liao SC, Chang CC. Identification of RNA editing sites in chloroplast transcripts of Phalaenopsis aphrodite and comparative analysis with those of other seed plants. Plant Cell Physiol. 2007;48:362–8.
Wang W, Zhang W, Wu Y, Maliga P, Messing J. RNA editing in chloroplasts of Spirodela polyrhiza, an aquatic Monocotelydonous species. PLoS One. 2015;10:e0140285.
Lin C-P, Ko C-Y, Kuo C-I, Liu M-S, Schafleitner R, Chen L-FO. Transcriptional slippage and RNA editing increase the diversity of transcripts in chloroplasts: insight from deep sequencing of Vigna radiata genome and transcriptome. PLoS One. 2015;10:e0129396.
He P, Huang S, Xiao G, Zhang Y, Yu J. Abundant RNA editing sites of chloroplast protein-coding genes in Ginkgo biloba and an evolutionary pattern analysis. BMC Plant Biol. 2016;16:257.
Grosche C, Funk HT, Maier UG, Zauner S. The chloroplast genome of Pellia endiviifolia: gene content, RNA-editing pattern, and the origin of chloroplast editing. Genome Biol Evol. 2012;4:1349–57.
Knie N, Grewe F, Fischer S, Knoop V. 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. BMC Evol Biol. 2016;16:134.
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.
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.
Richardson AO, Rice DW, Young GJ, Alverson AJ, Palmer JD. The “fossilized” mitochondrial genome of Liriodendron tulipifera: ancestral gene content and order, ancestral editing sites, and extraordinarily low mutation rate. BMC Biol. 2013;11:29.
Rice DW, Alverson AJ, Richardson AO, Young GJ, Sanchez-Puerta MV, Munzinger J, Barry K, Boore JL, Zhang Y, DePamphilis CW, Knox EB, Palmer JD. Horizontal transfer of entire genomes via mitochondrial fusion in the angiosperm Amborella. Science. 2013;342:1468–73.
Aljohi HA, Liu W, Lin Q, Zhao Y, Zeng J, Alamer A, Alanazi IO, Alawad AO, Al-Sadi AM, Hu S, Yu J. Complete sequence and analysis of coconut palm (Cocos nucifera) mitochondrial genome. PLoS One. 2016;11:e0163990.
Shearman JR, Sangsrakru D, Ruang-areerate P, Sonthirod C, Uthaipaisanwong P, Yoocha T, Poopear S, Theerawattanasuk K, Tragoonrung S, Tangphatsornruang S. Assembly and analysis of a male sterile rubber tree mitochondrial genome reveals DNA rearrangement events and a novel transcript. BMC Plant Biol. 2014;14:45.
Rodrigues NF, da Fonseca GC, Kulcheski FR, Margis R, Rodrigues NF, da Fonseca GC, Kulcheski FR, Margis R. Salt stress affects mRNA editing in soybean chloroplasts. Genet Mol Biol. 2017;40:200–8.
Silva SR, Diaz YCA, Penha HA, Pinheiro DG, Fernandes CC, Miranda VFO, Michael TP, Varani AM. The chloroplast genome of Utricularia reniformis sheds light on the evolution of the ndh gene complex of terrestrial carnivorous plants from the Lentibulariaceae Family. PLoS One. 2016;11:e0165176.
Calsa Júnior T, Carraro DM, Benatti MR, Barbosa AC, Kitajima JP, Carrer H. Structural features and transcript-editing analysis of sugarcane (Saccharum officinarum L.) chloroplast genome. Curr Genet. 2004;46:366–73.
Chen H, Deng L, Jiang Y, Lu P, Yu J. RNA editing sites exist in protein-coding genes in the chloroplast genome of Cycas taitungensis. J Integr Plant Biol. 2011;53:961–70.
Salmans ML, Chaw S-M, Lin C-P, Shih AC-C, Wu Y-W, Mulligan RM. Editing site analysis in a gymnosperm mitochondrial genome reveals similarities with angiosperm mitochondrial genomes. Curr Genet. 2010;56:439–46.
Yan L, Lai X, Li X, Wei C, Tan X, Zhang Y. Analyses of the complete genome and gene expression of chloroplast of sweet potato [Ipomoea batata]. PLoS One. 2015;10:e0124083.
Lee J, Kang Y, Shin SC, Park H, Lee H. Combined analysis of the chloroplast genome and transcriptome of the antarctic vascular plant Deschampsia antarctica Desv. PLoS One. 2014;9:e92501.
Uthaipaisanwong P, Chanprasert J, Shearman JR, Sangsrakru D, Yoocha T, Jomchai N, Jantasuriyarat C, Tragoonrung S, Tangphatsornruang S. Characterization of the chloroplast genome sequence of oil palm (Elaeis guineensis Jacq.). Gene. 2012;500:172–80.
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 Gen Genomics. 2009;281:473–81.
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 Cell Mol Biol. 2010;62:560–70.
Schallenberg-Rüdinger M, Kindgren P, Zehrmann A, Small I, Knoop V. A DYW-protein knockout in Physcomitrella affects two closely spaced mitochondrial editing sites and causes a severe developmental phenotype. Plant J. 2013;76:420–32.
Ichinose M, Sugita C, Yagi Y, Nakamura T, Sugita M. Two DYW subclass PPR proteins are involved in RNA editing of ccmFc and atp9 transcripts in the moss Physcomitrella patens: first complete set of PPR editing factors in plant mitochondria. Plant Cell Physiol. 2013;54:1907–16.
Stone JD, Storchova H. The application of RNA-seq to the comprehensive analysis of plant mitochondrial transcriptomes. Mol Gen Genomics. 2015;290:1–9.
Picardi E, D’Erchia AM, Montalvo A, Pesole G. Using REDItools to detect RNA editing events in NGS datasets. Curr Protoc Bioinforma. 2015;49:12.1–12.15.
Takenaka M, Zehrmann A, Brennicke A, Graichen K. Improved computational target site prediction for pentatricopeptide repeat RNA editing factors. PLoS One. 2013;8:e65343.
Filipovska A, Rackham O. Pentatricopeptide repeats: modular blocks for building RNA-binding proteins. RNA Biol. 2013;10:1426–32.
Castandet B, Araya A. The RNA editing pattern of cox2 mRNA is affected by point mutations in plant mitochondria. PLoS One. 2011;6:e20867.
Ichinose M, Uchida M, Sugita M. Identification of a pentatricopeptide repeat RNA editing factor in Physcomitrella patens chloroplasts. FEBS Lett. 2014;588:4060–4.
Schallenberg-Rüdinger M, Oldenkott B, Hiss M, Le TP, Knoop V, Rensing SA. A single-target mitochondrial RNA editing factor of Funaria hygrometrica can fully reconstitute RNA editing at two sites in Physcomitrella patens. Plant Cell Physiol. 2017;58:496–507.
Ichinose M, Sugita M. The DYW domains of pentatricopeptide repeat RNA editing factors contribute to discriminate target and non-target editing sites. Plant Cell Physiol. 2018;
Karpenahalli MR, Lupas AN, Soding J. TPRpred: a tool for prediction of TPR-, PPR- and SEL1-like repeats from protein sequences. BMC Bioinformatics. 2007;8:2.
Cheng S, Gutmann B, Zhong X, Ye Y, Fisher MF, Bai F, Castleden I, Song Y, Song B, Huang J, Liu X, Xu X, Lim BL, Bond CS, Yiu S-M, Small I. Redefining the structural motifs that determine RNA binding and RNA editing by pentatricopeptide repeat proteins in land plants. Plant J. 2016;85:532–47.
Hokamp K, Wolfe KH. PubCrawler: keeping up comfortably with PubMed and GenBank. Nucleic Acids Res. 2004;32(Web Server):W16–9.
Grant CE, Bailey TL, Noble WS. FIMO: scanning for occurrences of a given motif. Bioinformatics. 2011;27:1017–8.
Ambrosini G, Groux R, Bucher P. PWMScan: a fast tool for scanning entire genomes with a position-specific weight matrix. Bioinformatics. 2018;
The authors wish to acknowledge the contributions of the bioinformatics community to share useful tools which have helped to improve PREPACT as described under methods. V.K. is grateful to his university for providing the basic financial and infrastructural support sustaining work like the ongoing further development of PREPACT.
Availability of data and materials
All data generated or analysed during this study are included in this published article or integrated into the editome references of PREPACT as described, respectively.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1:
Multiple sequence input. An example for multiple query and reference sequence input in PREPACT’s alignment modes as discussed in the text. (DOCX 198 kb)
Additional file 2:
Table of re-checked edits. Verification of additional RNA editing events previously overlooked in editome references. (DOCX 45 kb)
Additional file 3:
Alignment prediction output. An example for the output of a multiple-query alignment for different references. (DOCX 79 kb)
Additional file 4:
TargetScan of Editing environment. An example illustrating the use of TargetScan to identify nucleotide bias in the immediate environment of editing sites in position -4 to +3 as discussed in the text. (DOCX 87 kb)
Rights and permissions
Open Access This 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.
About this article
Cite this article
Lenz, H., Hein, A. & Knoop, V. Plant organelle RNA editing and its specificity factors: enhancements of analyses and new database features in PREPACT 3.0. BMC Bioinformatics 19, 255 (2018). https://doi.org/10.1186/s12859-018-2244-9
- Pentatricopeptide repeat (PPR) proteins
- Pyrimidine exchange RNA editing
- RNA-binding proteins