Frequent chloroplast RNA editing in early-branching flowering plants: pilot studies on angiosperm-wide coexistence of editing sites and their nuclear specificity factors

Conservation of the moss RNA editing factor PPR78 despite the loss of its known cytidine-to-uridine editing sites is explained by a hidden extra target

Cytidine (C)-to-uridine (U) RNA editing in plant organelles relies on specific RNA-binding pentatricopeptide repeat (PPR) proteins. In the moss Physcomitrium patens, all such RNA editing factors feature a C-terminal DYW domain that acts as the cytidine deaminase for C-to-U conversion. PPR78 of Physcomitrium targets 2 mitochondrial editing sites, cox1eU755SL and rps14eU137SL. Remarkably, the latter is edited to highly variable degrees in different mosses. Here, we aimed to unravel the coevolution of PPR78 and its 2 target sites in mosses. Heterologous complementation in a Physcomitrium knockout line revealed that the variable editing of rps14eU137SL depends on the PPR arrays of different PPR78 orthologues but not their C-terminal domains. Intriguingly, PPR78 has remained conserved despite the simultaneous loss of editing at both known targets among Hypnales (feather mosses), suggesting it serves an additional function. Using a recently established RNA editing assay in Escherichia coli, we confirmed site-specific RNA editing by PPR78 in the bacterium and identified 4 additional off-targets in the bacterial transcriptome. Based on conservation profiles, we predicted ccmFNeU1465RC as a candidate editing target of PPR78 in moss mitochondrial transcriptomes. We confirmed editing at this site in several mosses and verified that PPR78 targets ccmFNeU1465RC in the bacterial editing system, explaining the conservation and functional adaptation of PPR78 during moss evolution.

Emergence of Novel RNA-Editing Sites by Changes in the Binding Affinity of a Conserved PPR Protein

RNA editing converts cytidines to uridines in plant organellar transcripts. Editing typically restores codons for conserved amino acids. During evolution, specific C-to-U editing sites can be lost from some plant lineages by genomic C-to-T mutations. By contrast, the emergence of novel editing sites is less well documented. Editing sites are recognized by pentatricopeptide repeat (PPR) proteins with high specificity. RNA recognition by PPR proteins is partially predictable, but prediction is often inadequate for PPRs involved in RNA editing. Here we have characterized evolution and recognition of a recently gained editing site. We demonstrate that changes in the RNA recognition motifs that are not explainable with the current PPR code allow an ancient PPR protein, QED1, to uniquely target the ndhB-291 site in Brassicaceae. When expressed in tobacco, the Arabidopsis QED1 edits 33 high-confident off-target sites in chloroplasts and mitochondria causing a spectrum of mutant phenotypes. By manipulating the relative expression levels of QED1 and ndhB-291, we show that the target specificity of the PPR protein depends on the RNA:protein ratio. Finally, our data suggest that the low expression levels of PPR proteins are necessary to ensure the specificity of editing site selection and prevent deleterious off-target editing.

A Comprehensive Evolutionary Study of Chloroplast RNA Editing in Gymnosperms: A Novel Type of G-to-A RNA Editing Is Common in Gymnosperms

Although more than 9100 plant plastomes have been sequenced, RNA editing sites of the whole plastome have been experimentally verified in only approximately 21 species, which seriously hampers the comprehensive evolutionary study of chloroplast RNA editing. We investigated the evolutionary pattern of chloroplast RNA editing sites in 19 species from all 13 families of gymnosperms based on a combination of genomic and transcriptomic data. We found that the chloroplast C-to-U RNA editing sites of gymnosperms shared many common characteristics with those of other land plants, but also exhibited many unique characteristics. In contrast to that noted in angiosperms, the density of RNA editing sites in ndh genes was not the highest in the sampled gymnosperms, and both loss and gain events at editing sites occurred frequently during the evolution of gymnosperms. In addition, GC content and plastomic size were positively correlated with the number of chloroplast RNA editing sites in gymnosperms, suggesting that the increase in GC content could provide more materials for RNA editing and facilitate the evolution of RNA editing in land plants or vice versa. Interestingly, novel G-to-A RNA editing events were commonly found in all sampled gymnosperm species, and G-to-A RNA editing exhibits many different characteristics from C-to-U RNA editing in gymnosperms. This study revealed a comprehensive evolutionary scenario for chloroplast RNA editing sites in gymnosperms, and reported that a novel type of G-to-A RNA editing is prevalent in gymnosperms.

Plastid RNA editing reduction accompanied with genetic variations in Cymbidium, a genus with diverse lifestyle modes

Recent sequencing efforts have broadly uncovered the evolutionary trajectory of plastid genomes (plastomes) of flowering plants in diverse habitats, yet our knowledge of the evolution of plastid posttranscriptional modifications is limited. In this study, we generated 11 complete plastomes and performed ultra-deep transcriptome sequencing to investigate the co-evolution of plastid RNA editing and genetic variation in Cymbidium, a genus with diverse trophic lifestyles. Genome size and gene content is reduced in terrestrial and green mycoheterotrophic orchids relative to their epiphytic relatives. This could be partly due to extensive losses and pseudogenization of ndh genes for the plastid NADH dehydrogenase-like complex, but independent pseudogenization of ndh genes has also occurred in the epiphyte C. mannii, which was reported to use strong crassulacean acid metabolism photosynthesis. RNA editing sites are abundant but variable in number among Cymbidium plastomes. The nearly twofold variation in editing abundance is mainly due to extensive reduction of ancestral editing sites in ndh transcripts of terrestrial, mycoheterotrophic, and C. mannii plastomes. The co-occurrence of editing reduction and pseudogenization in ndh genes suggests functional constraints on editing machinery may be relaxed, leading to nonrandom loss of ancestral edited sites via reduced editing efficiency. This study represents the first systematic examination of RNA editing evolution linked to plastid genome variation in a single genus. We also propose an explanation for how genomic and posttranscriptional variations might be affected by lifestyle-associated ecological adaptation strategies in Cymbidium.

The complete plastid genome sequence of the enigmatic moss, Takakia lepidozioides (Takakiopsida, Bryophyta): evolutionary perspectives on the largest collection of genes in mosses and the intensive RNA editing

Complete chloroplast genome sequence of a moss, Takakia lepidozioides (Takakiopsida) is reported. The largest collection of genes in mosses and the intensive RNA editing were discussed from evolutionary perspectives. We assembled the entire plastid genome sequence of Takakia lepidozioides (Takakiopsida), emerging from the first phylogenetic split among extant mosses. The genome sequences were assembled into a circular molecule 149,016 bp in length, with a quadripartite structure comprising a large and a small single-copy region separated by inverted repeats. It contained 88 genes coding for proteins, 32 for tRNA, four for rRNA, two open reading frames, and at least one pseudogene (tufA). This is the largest number of genes of all sequenced plastid genomes in mosses and Takakia is the only moss that retains the seven coding genes ccsA, cysA, cysT, petN rpoA, rps16 and trnP^GGG. Parsimonious interpretation of gene loss suggests that the last common ancestor of bryophytes had all seven genes and that mosses lost at least three of them during their diversification. Analyses of the plastid transcriptome identified the extraordinary frequency of RNA editing with more than 1100 sites. We indicated a close correlation between the monoplastidy of vegetative tissue and the intensive RNA editing sites in the plastid genome in land plant lineages. Here, we proposed a hypothesis that the small population size of plastids in each vegetative cell of some early diverging land plants, including Takakia, might cause the frequent fixation of mutations in plastid genome through the intracellular genetic drift and that deleterious mutations might be continuously compensated by RNA editing during or following transcription.

Moss PPR-SMR protein PpPPR_64 influences the expression of a psaA-psaB-rps14 gene cluster and processing of the 23S-4.5S rRNA precursor in chloroplasts

Moss PPR-SMR protein PpPPR_64 is a pTAC2 homolog but is functionally distinct from pTAC2. PpPPR_64 is required for psaA gene expression and its function may have evolved in mosses. The pentatricopeptide repeat (PPR) proteins are key regulatory factors responsible for the control of plant organellar gene expression. A small subset of PPR proteins possess a C-terminal small MutS-related (SMR) domain and have diverse roles in plant organellar biogenesis. However, the function of PPR-SMR proteins is not fully understood. Here, we report the function of PPR-SMR protein PpPPR_64 in the moss Physcomitrium patens. Phylogenetic analysis indicated that PpPPR_64 belongs to the same clade as the Arabidopsis PPR-SMR protein pTAC2. PpPPR_64 knockout (KO) mutants grew autotrophically but with reduced protonemata growth and the poor formation of photosystems' antenna complexes. Quantitative reverse transcription-polymerase chain reaction and RNA gel blot hybridization analyses revealed a significant reduction in transcript levels of the psaA-psaB-rps14 gene cluster but no alteration to transcript levels of most photosynthesis- and non-photosynthesis-related genes. In addition, RNA processing of 23S-4.5S rRNA precursor was impaired in the PpPPR_64 KO mutants. This suggests that PpPPR_64 is specifically involved in the expression level of the psaA-psaB-rps14 gene and in processing of the 23S-4.5S rRNA precursor. Our results indicate that PpPPR_64 is functionally distinct from pTAC2 and is a novel PPR-SMR protein required for proper chloroplast biogenesis in P. patens.

Advancing organelle genome transformation and editing for crop improvement

Plant cells contain three organelles that harbor DNA: the nucleus, plastids, and mitochondria. Plastid transformation has emerged as an attractive platform for the generation of transgenic plants, also referred to as transplastomic plants. Plastid genomes have been genetically engineered to improve crop yield, nutritional quality, and resistance to abiotic and biotic stresses, as well as for recombinant protein production. Despite many promising proof-of-concept applications, transplastomic plants have not been commercialized to date. Sequence-specific nuclease technologies are widely used to precisely modify nuclear genomes, but these tools have not been applied to edit organelle genomes because the efficient homologous recombination system in plastids facilitates plastid genome editing. Unlike plastid transformation, successful genetic transformation of higher plant mitochondrial genome transformation was tested in several research group, but not successful to date. However, stepwise progress has been made in modifying mitochondrial genes and their transcripts, thus enabling the study of their functions. Here, we provide an overview of advances in organelle transformation and genome editing for crop improvement, and we discuss the bottlenecks and future development of these technologies.

RNA editing may stabilize membrane-embedded proteins by increasing phydrophobicity: A study of Zanthoxylum piperitum and Z. schinifolium chloroplast NdhG

Most chloroplast genes in Zanthoxylum schinifolium (Korean pepper) and Z. piperitum (Japanese pepper) are subject to neutral or purifying selection (Ka/Ks values < 1); however, NAD(P)H dehydrogenase subunit G (ndhG) has a Ka/Ks value of 1.43, which may indicate positive selection. Here, we modeled the ZsNdhG and ZpNdhG structures by comparing them with the NuoJ subunit of respiratory complex I in Escherichia coli, revealing the locations of four amino acid differences between ZsNdhG and ZpNdhG. As these polymorphisms were located at the end of a membrane-spanning α-helix or in peptide loops external to the membrane, they are not expected to have major effects on the membrane-embedding properties of these proteins. However, we found that C-to-U RNA editing occurred at the ndhG-50 sites of ndhG (uCa to uUa, Ser to Leu) in both species, resulting in changes to an amino acid located in the middle of a membrane-spanning α-helix, which may maintain its hydrophobicity. RNA editing at the ndhG-50 site was conserved in many plant species, and the modeled structures of Anthoceros formosae NdhG and Spirodela polyrhiza NdhB provided further evidence that RNA editing increases the hydrophobicity of membrane-embedded proteins. We also speculate that the polar residues inside membrane-spanning α-helices serve to support the protein structure. This report represents the first RNA-editing site identified in Zanthoxylum and points to the importance of considering RNA editing when identifying positively selected genes based on Ka/Ks values.

Plant organellar RNA editing: what 30 years of research has revealed

The central dogma in biology defines the flow of genetic information from DNA to RNA to protein. Accordingly, RNA molecules generally accurately follow the sequences of the genes from which they are transcribed. This rule is transgressed by RNA editing, which creates RNA products that differ from their DNA templates. Analyses of the RNA landscapes of terrestrial plants have indicated that RNA editing (in the form of C-U base transitions) is highly prevalent within organelles (that is, mitochondria and chloroplasts). Numerous C→U conversions (and in some plants also U→C) alter the coding sequences of many of the organellar transcripts and can also produce translatable mRNAs by creating AUG start sites or eliminating premature stop codons, or affect the RNA structure, influence splicing and alter the stability of RNAs. RNA-binding proteins are at the heart of post-transcriptional RNA expression. The C-to-U RNA editing process in plant mitochondria involves numerous nuclear-encoded factors, many of which have been identified as pentatricopeptide repeat (PPR) proteins that target editing sites in a sequence-specific manner. In this review we report on major discoveries on RNA editing in plant organelles, since it was first documented 30 years ago.

The Expansion and Diversification of Pentatricopeptide Repeat RNA-Editing Factors in Plants

The RNA-binding pentatricopeptide repeat (PPR) family comprises hundreds to thousands of genes in most plants, but only a few dozen in algae, indicating massive gene expansions during land plant evolution. The nature and timing of these expansions has not been well defined due to the sparse sequence data available from early-diverging land plant lineages. In this study, we exploit the comprehensive OneKP datasets of over 1000 transcriptomes from diverse plants and algae toward establishing a clear picture of the evolution of this massive gene family, focusing on the proteins typically associated with RNA editing, which show the most spectacular variation in numbers and domain composition across the plant kingdom. We characterize over 2 250 000 PPR motifs in over 400 000 proteins. In lycophytes, polypod ferns, and hornworts, nearly 10% of expressed protein-coding genes encode putative PPR editing factors, whereas they are absent from algae and complex-thalloid liverworts. We show that rather than a single expansion, most land plant lineages with high numbers of editing factors have continued to generate novel sequence diversity. We identify sequence variations that imply functional differences between PPR proteins in seed plants versus non-seed plants and variations we propose to be linked to seed-plant-specific editing co-factors. Finally, using the sequence variations across the datasets, we develop a structural model of the catalytic DYW domain associated with C-to-U editing and identify a clade of unique DYW variants that are strong candidates as U-to-C RNA-editing factors, given their phylogenetic distribution and sequence characteristics.

The complete chloroplast genome of greater duckweed (Spirodela polyrhiza 7498) using PacBio long reads: insights into the chloroplast evolution and transcription regulation

Background

Duckweeds (Lemnaceae) are aquatic plants distributed all over the world. The chloroplast genome, as an efficient solar-powered reactor, is an invaluable resource to study biodiversity and to carry foreign genes. The chloroplast genome sequencing has become routine and less expensive with the delivery of high-throughput sequencing technologies, allowing us to deeply investigate genomics and transcriptomics of duckweed organelles.

Results

Here, the complete chloroplast genome of Spirodela polyrhiza 7498 (SpV2) is assembled by PacBio sequencing. The length of 168,956 bp circular genome is composed of a pair of inverted repeats of 31,844 bp, a large single copy of 91,210 bp and a small single copy of 14,058 bp. Compared to the previous version (SpV1) assembled from short reads, the integrity and quality of SpV2 are improved, especially with the retrieval of two repeated fragments in ycf2 gene. There are a number of 107 unique genes, including 78 protein-coding genes, 25 tRNA genes and 4 rRNA genes. With the evidence of full-length cDNAs generated from PacBio isoform sequencing, seven genes (ycf3, clpP, atpF, rpoC1, rpl2, rps12 and ndhA) are detected to contain type-II introns. The ndhA intron has 50% more sequence divergence than the species-barcoding marker of atpF-atpH, showing the potential power to discriminate close species. A number of 37 RNA editing sites are recognized to have cytosine (C) to uracil (U) substitutions, eight of which are newly defined including six from the intergenic regions and two from the coding sequences of rpoC2 and ndhA genes. In addition, nine operon classes are identified using transcriptomic data. It is found that the operons contain multiple subunit genes encoding the same functional complexes comprising of ATP synthase, photosynthesis system, ribosomal proteins, et.al., which could be simultaneously transcribed and coordinately translated in response to the cell stimuli.

Conclusions

The understanding of the chloroplast genomics and the transcriptomics of S.polyrhiza would greatly facilitate the study of phylogenetic evolution and the application of genetically engineering duckweeds.

The dual-targeted RNA editing factor AEF1 is universally conserved among angiosperms and reveals only minor adaptations upon loss of its chloroplast or its mitochondrial target

Upon loss of either its chloroplast or mitochondrial target, a uniquely dual-targeted factor for C-to-U RNA editing in angiosperms reveals low evidence for improved molecular adaptation to its remaining target. RNA-binding pentatricopeptide repeat (PPR) proteins specifically recognize target sites for C-to-U RNA editing in the transcriptomes of plant chloroplasts and mitochondria. Among more than 80 PPR-type editing factors that have meantime been characterized, AEF1 (or MPR25) is a special case given its dual targeting to both organelles and addressing an essential mitochondrial (nad5eU1580SL) and an essential chloroplast (atpFeU92SL) RNA editing site in parallel in Arabidopsis. Here, we explored the angiosperm-wide conservation of AEF1 and its two organelle targets. Despite numerous independent losses of the chloroplast editing site by C-to-T conversion and at least four such conversions at the mitochondrial target site in other taxa, AEF1 remains consistently conserved in more than 120 sampled angiosperm genomes. Not a single case of simultaneous loss of the chloroplast and mitochondrial editing target or of AEF1 disintegration or loss could be identified, contrasting previous findings for editing factors targeted to only one organelle. Like in most RNA editing factors, the PPR array of AEF1 reveals potential for conceptually "improved fits" to its targets according to the current PPR-RNA binding code. Surprisingly, we observe only minor evidence for adaptation to the mitochondrial target also after deep losses of the chloroplast target among Asterales, Caryophyllales and Poales or, vice versa, for the remaining chloroplast target after a deep loss of the mitochondrial target among Malvales. The evolutionary observations support the notion that PPR-RNA mismatches may be essential for proper function of editing factors.

Variation in plastid genomes in the gynodioecious species Silene vulgaris

BACKGROUND

Gynodioecious species exist in two sexes - male-sterile females and hermaphrodites. Male sterility in higher plants often results from mitonuclear interaction between the CMS (cytoplasmic male sterility) gene(s) encoded by mitochondrial genome and by nuclear-encoded restorer genes. Mitochondrial and nuclear-encoded transcriptomes in females and hermaphrodites are intensively studied, but little is known about sex-specific gene expression in plastids. We have compared plastid transcriptomes between females and hermaphrodites in two haplotypes of a gynodioecious species Silene vulgaris with known CMS candidate genes.

RESULTS

We generated complete plastid genome sequences from five haplotypes S. vulgaris including the haplotypes KRA and KOV, for which complete mitochondrial genome sequences were already published. We constructed a phylogenetic tree based on plastid sequences of S. vulgaris. Whereas lowland S. vulgaris haplotypes including KRA and KOV clustered together, the accessions from high European mountains diverged early in the phylogram. S. vulgaris belongs among Silene species with slowly evolving plastid genomes, but we still detected 212 substitutions and 112 indels between two accessions of this species. We estimated elevated Ka/Ks in the ndhF gene, which may reflect the adaptation of S. vulgaris to high altitudes, or relaxed selection. We compared depth of coverage and editing rates between female and hermaphrodite plastid transcriptomes and found no significant differences between the two sexes. We identified 51 unique C to U editing sites in the plastid genomes of S. vulgaris, 38 of them in protein coding regions, 2 in introns, and 11 in intergenic regions. The editing site in the psbZ gene was edited only in one of two plastid genomes under study.

CONCLUSIONS

We revealed no significant differences between the sexes in plastid transcriptomes of two haplotypes of S. vulgaris. It suggests that gene expression of plastid genes is not affected by CMS in flower buds of S. vulgaris, although both sexes may still differ in plastid gene expression in specific tissues. We revealed the difference between the plastid transcriptomes of two S. vulgaris haplotypes in editing rate and in the coverage of several antisense transcripts. Our results document the variation in plastid genomes and transcriptomes in S. vulgaris.

Establishment of a Heterologous RNA Editing Event in Chloroplasts

In chloroplasts and plant mitochondria, specific cytidines in mRNAs are posttranscriptionally converted to uridines by RNA editing. Editing sites are recognized by nucleus-encoded RNA-binding proteins of the pentatricopeptide repeat (PPR) family, which bind upstream of the editing site in a sequence-specific manner and direct the editing activity to the target position. Editing sites have been lost many times during evolution by C-to-T mutations. Loss of an editing site is thought to be accompanied by loss or degeneration of its cognate PPR protein. Consequently, foreign editing sites are usually not recognized when introduced into species lacking the site. Previously, the spinach (Spinacia oleracea) psbF-26 editing site was introduced into the tobacco (Nicotiana tabacum) plastid genome. Tobacco lacks the psbF-26 site and cannot edit it. Expression of the "unedited" PsbF protein resulted in impaired PSII function. In Arabidopsis (Arabidopsis thaliana), the PPR protein LPA66 is required for editing at psbF-26. Here, we show that introduction of the Arabidopsis LPA66 reconstitutes editing of the spinach psbF-26 site in tobacco and restores a wild-type-like phenotype. Our findings define the minimum requirements for establishing new RNA editing sites and suggest that the evolutionary dynamics of editing patterns is largely explained by coevolution of editing sites and PPR proteins.

Evolutionary Model of Plastidial RNA Editing in Angiosperms Presumed from Genome-Wide Analysis of Amborella trichopoda

Amborella trichopoda is placed close to the base of the angiosperm lineage (basal angiosperm). By genome-wide RNA sequencing, we identified 184C-to-U RNA editing sites in the plastid genome of Amborella. This number is much higher than that observed in other angiosperms including maize (44 sites), rice (39 sites) and grape (115 sites). Despite the high frequency of RNA editing, the biased distribution of RNA editing sites in the genome, target codon preference and nucleotide preference adjacent to the edited cytidine are similar to that in other angiosperms, suggesting a common editing machinery. Consistent with this idea, the Amborella nuclear genome encodes 2-3 times more of the E- and DYW-subclass members of pentatricopeptide repeat proteins responsible for RNA editing site recognition in plant organelles. Among 165 editing sites in plastid protein coding sequences in Amborella, 100 sites were conserved at least in one out of 38 species selected to represent key branching points of the angiosperm phylogenetic tree. We assume these 100 sites represent at least a subset of the sites in the plastid editotype of ancestral angiosperms. We then mapped the loss and gain of editing sites on the phylogenetic tree of angiosperms. Our results support the idea that the evolution of angiosperms has led to the loss of RNA editing sites in plastids.

Identification and analysis of RNA editing sites in chloroplast transcripts of kenaf (Hibiscus cannabinus L.)

RNA editing is one of the post-transcriptional modification processes and can lead to changes in sequencing and functioning of corresponding proteins and genetic information. To reveal the composition and characteristic of RNA editing of kenaf chloroplast genome, the RNA editing sites in kenaf chloroplast were predicted and identified using bioinformatics and RT-PCR analysis. The prediction results showed a total of 48 editing sites distributed in 22 genes, all of them were C to U conversion leading to amino acid changes. Further analysis of the position of RNA editing sites revealed that except 11 editing sites located at the first codon base, the other editing sites were found at the second codon base. Then four genes were randomly selected to validate the editing sites. Results showed that it was accurate to study the chloroplast RNA editing sites by bioinformatics method accompanied with cloning sequencing. Furthermore, the protein secondary structure and transmembrane domain of ndhD and atpA that had undergone gene editing also changed after editing. This implied that proteins with structural changes may have an impact on kenaf growth. Meanwhile, the differential editing site was found in chloroplast transcripts in kenaf CMS line and its maintainer line, indicating that chloroplast RNA editing could be associated with kenaf CMS. Therefore, the present study laid a foundation to further reveal the biological functioning of chloroplast RNA editing in CMS and its maintainer lines in kenaf.

The complete chloroplast genome sequence of watercress (Nasturtium officinale R. Br.): Genome organization, adaptive evolution and phylogenetic relationships in Cardamineae

Watercress (Nasturtium officinale R. Br.), an aquatic leafy vegetable of the Brassicaceae family, is known as a nutritional powerhouse. Here, we de novo sequenced and assembled the complete chloroplast (cp) genome of watercress based on combined PacBio and Illumina data. The cp genome is 155,106 bp in length, exhibiting a typical quadripartite structure including a pair of inverted repeats (IRA and IRB) of 26,505 bp separated by a large single copy (LSC) region of 84,265 bp and a small single copy (SSC) region of 17,831 bp. The genome contained 113 unique genes, including 79 protein-coding genes, 30 tRNAs and 4 rRNAs, with 20 duplicate in the IRs. Compared with the prior cp genome of watercress deposited in GenBank, 21 single nucleotide polymorphisms (SNPs) and 27 indels were identified, mainly located in noncoding sequences. A total of 49 repeat structures and 71 simple sequence repeats (SSRs) were detected. Codon usage showed a bias for A/T-ending codons in the cp genome of watercress. Moreover, 45 RNA editing sites were predicted in 16 genes, all for C-to-U transitions. A comparative plastome study with Cardamineae species revealed a conserved gene order and high similarity of protein-coding sequences. Analysis of the Ka/Ks ratios of Cardamineae suggested positive selection exerted on the ycf2 gene in watercress, which might reflect specific adaptations of watercress to its particular living environment. Phylogenetic analyses based on complete cp genomes and common protein-coding genes from 56 species showed that the genus Nasturtium was a sister to Cardamine in the Cardamineae tribe. Our study provides valuable resources for future evolution, population genetics and molecular biology studies of watercress.

Multifarious Evolutionary Pathways of a Nuclear RNA Editing Factor: Disjunctions in Coevolution of DOT4 and Its Chloroplast Target rpoC1eU488SL

Nuclear-encoded pentatricopeptide repeat (PPR) proteins are site-specific factors for C-to-U RNA editing in plant organelles coevolving with their targets. Losing an editing target by C-to-T conversion allows for eventual loss of its editing factor, as recently confirmed for editing factors CLB19, CRR28, and RARE1 targeting ancient chloroplast editing sites in flowering plants. Here, we report on alternative evolutionary pathways for DOT4 addressing rpoC1eU488SL, a chloroplast editing site in the RNA polymerase β' subunit mRNA. Upon loss of rpoC1eU488SL by C-to-T conversion, DOT4 got lost multiple times independently in angiosperm evolution with intermediate states of DOT4 orthologs in various stages of degeneration. Surprisingly, we now also observe degeneration and loss of DOT4 despite retention of a C in the editing position (in Carica, Coffea, Vicia, and Spirodela). We find that the cytidine remains unedited, proving that DOT4 was not replaced by another editing factor. Yet another pathway of DOT4 evolution is observed among the Poaceae. Although the rpoC1eU488SL edit has been lost through C-to-T conversion, DOT4 orthologs not only remain conserved but also have their array of PPRs extended by six additional repeats. Here, the loss of the ancient target has likely allowed DOT4 to adapt for a new function. We suggest rps3 antisense transcripts as previously demonstrated in barley (Hordeum vulgare) arising from promotor sequences newly emerging in the rpl16 intron of Poaceae as a new candidate target for the extended PPR stretch of DOT4. Altogether, DOT4 and its target show more flexible pathways for evolution than the previously explored editing factors CLB19, CRR28, and RARE1. Certain plant clades (e.g., Amaranthus, Vaccinium, Carica, the Poaceae, Fabales, and Caryophyllales) show pronounced dynamics in the evolution of editing sites and corresponding factors.

Divergence of RNA editing among Arabidopsis species

RNA editing altered the RNA sequence by replacing the C nucleotide to U in the organellar genomes of plants. RNA editing status sometimes differed among distant species. The pattern of conservation and variation of RNA editing status made it possible to evaluate evolutionary mechanisms impacting functional aspects of RNA editing. In this study, divergence of RNA editing in the chloroplast genome among Arabidopsis species was analyzed to determine 9 losses and 1 gain in RNA editing. All changes in A. thaliana lineage resulted from changes to the chloroplast genome sequence, whereas changes in the A. lyrata / halleri lineage were possibly due to exclusive changes in the nuclear editing factors. One loss of RNA editing in A. lyrata was caused by a deficiency in the PPR gene OTP80. The changes in RNA editing occurred approximately every two million years and were not observed at functionally important sites. These results highlight the conserved nature of RNA editing status suggesting the importance of RNA editing during evolution.

ReFernment: An R package for annotating RNA editing in plastid genomes

PREMISE OF THE STUDY

In the absence of cDNA, the annotation of RNA editing in plastomes must be done manually, representing a significant time cost to those studying the organellar genomes of ferns and hornworts.

METHODS AND RESULTS

We developed an R package to automatically annotate apparent nonsense mutations in plastid genomes. The software successfully annotates such sites and results in no false positives for data with no sequencing or assembly errors.

CONCLUSIONS

Compared to manual annotation, ReFernment offers greater speed and accuracy for annotating RNA editing sites. This software should be especially useful for researchers generating large numbers of plastome sequences for taxa with high levels of RNA editing.

Complete Chloroplast Genome Sequence of Broomcorn Millet (Panicum miliaceum L.) and Comparative Analysis with Other Panicoideae Species

Broomcorn millet (Panicum miliaceum L.) is one of the earliest domesticated cereals worldwide, holding significant agricultural, historical, and evolutionary importance. However, our genomic knowledge of it is rather limited at present, hampering further genetic and evolutionary studies. Here, we sequenced and assembled the chloroplast genome (cp) of broomcorn millet and compared it with five other Panicoideae species. Results showed that the cp genome of broomcorn millet was 139,826 bp in size, with a typical quadripartite structure. In total, 108 genes were annotated and 18 genes were duplicated in the IR (inverted region) region, which was similar to other Panicoideae species. Comparative analysis showed a rather conserved genome structure between them, with three common regions. Furthermore, RNA editing, codon usage, and expansion of the IR, as well as simple sequence repeat (SSR) elements, were systematically investigated and 13 potential DNA markers were developed for Panicoideae species identification. Finally, phylogenetic analysis implied that broomcorn millet was a sister species to Panicum virgatum within the tribe Paniceae, and supported a monophyly of the Panicoideae. This study has reported for the first time the genome organization, gene content, and structural features of the chloroplast genome of broomcorn millet, which provides valuable information for genetic and evolutionary studies in the genus Panicum and beyond.

Genome-wide organellar analyses from the hornwort Leiosporoceros dussii show low frequency of RNA editing

Because hornworts occupy a pivotal position in early land colonization as sister to other bryophytes, sister to tracheophytes, or sister to all other land plants, a renewed interest has arisen in their phylogenetic diversity, morphology, and genomes. To date, only five organellar genome sequences are available for hornworts. We sequenced the plastome (155,956 bp) and mitogenome (212,153 bp) of the hornwort Leiosporoceros dussii, the sister taxon to all hornworts. The Leiosporoceros organellar genomes show conserved gene structure and order with respect to the other hornworts and other bryophytes. Additionally, using RNA-seq data we quantified the frequency of RNA-editing events (the canonical C-to-U and the reverse editing U-to-C) in both organellar genomes. In total, 109 sites were found in the plastome and 108 in the mitogenome, respectively. The proportion of edited sites corresponds to 0.06% of the plastome and 0.05% of the mitogenome (in reference to the total genome size), in contrast to 0.58% of edited sites in the plastome of Anthoceros angustus (161,162 bp). All edited sites in the plastome and 88 of 108 sites in the mitogenome are C-to-U conversions. Twenty reverse edited sites (U-to-C conversions) were found in the mitogenome (17.8%) and none in the plastome. The low frequency of RNA editing in Leiosporoceros, which is nearly 88% less than in the plastome of Anthoceros and the mitogenome of Nothoceros, indicates that the frequency of RNA editing has fluctuated during hornwort diversification. Hornworts are a pivotal land plant group to unravel the genomic implications of RNA editing and its maintenance despite the evident evolutionary disadvantages.

Plant organelle RNA editing and its specificity factors: enhancements of analyses and new database features in PREPACT 3.0

Background

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.

Results

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.

Conclusions

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.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2244-9) contains supplementary material, which is available to authorized users.

Expected and unexpected evolution of plant RNA editing factors CLB19, CRR28 and RARE1: retention of CLB19 despite a phylogenetically deep loss of its two known editing targets in Poaceae

BACKGROUND

C-to-U RNA editing in mitochondria and chloroplasts and the nuclear-encoded, RNA-binding PPR proteins acting as editing factors present a wide field of co-evolution between the different genetic systems in a plant cell. Recent studies on chloroplast editing factors RARE1 and CRR28 addressing one or two chloroplast editing sites, respectively, found them strictly conserved among 65 flowering plants as long as one of their RNA editing targets remained present.

RESULTS

Extending the earlier sampling to 117 angiosperms with high-quality genome or transcriptome data, we find more evidence confirming previous conclusions but now also identify cases for expected evolutionary transition states such as retention of RARE1 despite loss of its editing target or the degeneration of CRR28 truncating its carboxyterminal DYW domain. The extended angiosperm set was now used to explore CLB19, an "E+"-type PPR editing factor targeting two chloroplast editing sites, rpoAeU200SF and clpPeU559HY, in Arabidopsis thaliana. We found CLB19 consistently conserved if one of the two targets was retained and three independent losses of CLB19 after elimination of both targets. The Ericales show independent regains of the ancestrally lost clpPeU559HY editing, further explaining why multiple-target editing factors are lost much more rarely than single target factors like RARE1. The retention of CLB19 despite loss of both editing targets in some Ericaceae, Apocynaceae and in Camptotheca (Nyssaceae) likely represents evolutionary transitions. However, the retention of CLB19 after a phylogenetic deep loss in the Poaceae rather suggests a yet unrecognized further editing target, for which we suggest editing event ndhAeU473SL.

CONCLUSION

Extending the scope of studies on plant organelle RNA editing to further taxa and additional nuclear cofactors reveals expected evolutionary transitions, strikingly different evolutionary dynamics for multiple-target editing factors like CLB19 and CRR28 and suggests additional functions for editing factor CLB19 among the Poaceae.

Comparative assessment of chloroplast transcriptional responses highlights conserved and unique patterns across Triticeae members under salt stress

Chloroplast functional genomics, in particular understanding the chloroplast transcriptional response is of immense importance mainly due to its role in oxygenic photosynthesis. As a photosynthetic unit, its efficiency and transcriptional activity is directly regulated by reactive oxygen species during abiotic and biotic stress and subsequently affects carbon assimilation, and plant biomass. In crops, understanding photosynthesis is crucial for crop domestication by identifying the traits that could be exploited for crop improvement. Transcriptionally and translationally active chloroplast plays a key role by regulating the PSI and PSII photo-reaction centres, which ubiquitously affects the light harvesting. Using a comparative transcriptomics mapping approach, we identified differential regulation of key chloroplast genes during salt stress across Triticeae members with potential genes involved in photosynthesis and electron transport system such as CytB6f. Apart from differentially regulated genes involved in PSI and PSII, we found widespread evidence of intron splicing events, specifically uniquely spliced petB and petD in Triticum aestivum and high proportion of RNA editing in ndh genes across the Triticeae members during salt stress. We also highlight the role and differential regulation of ATP synthase as member of CF₀CF₁ and also revealed the effect of salt stress on the water-splitting complex under salt stress. It is worthwhile to mention that the observed conserved down-regulation of psbJ across the Triticeae is limiting the assembly of water-splitting complexes and thus making the BEP clade Triticeae members more vulnerable to high light during the salt stress. Comparative understanding of the chloroplast transcriptional dynamics and photosynthetic regulation will improve the approaches for improved crop domestication.

The complete plastome of macaw palm [Acrocomia aculeata (Jacq.) Lodd. ex Mart.] and extensive molecular analyses of the evolution of plastid genes in Arecaceae

The plastome of macaw palm was sequenced allowing analyses of evolution and molecular markers. Additionally, we demonstrated that more than half of plastid protein-coding genes in Arecaceae underwent positive selection. Macaw palm is a native species from tropical and subtropical Americas. It shows high production of oil per hectare reaching up to 70% of oil content in fruits and an interesting plasticity to grow in different ecosystems. Its domestication and breeding are still in the beginning, which makes the development of molecular markers essential to assess natural populations and germplasm collections. Therefore, we sequenced and characterized in detail the plastome of macaw palm. A total of 221 SSR loci were identified in the plastome of macaw palm. Additionally, eight polymorphism hotspots were characterized at level of subfamily and tribe. Moreover, several events of gain and loss of RNA editing sites were found within the subfamily Arecoideae. Aiming to uncover evolutionary events in Arecaceae, we also analyzed extensively the evolution of plastid genes. The analyses show that highly divergent genes seem to evolve in a species-specific manner, suggesting that gene degeneration events may be occurring within Arecaceae at the level of genus or species. Unexpectedly, we found that more than half of plastid protein-coding genes are under positive selection, including genes for photosynthesis, gene expression machinery and other essential plastid functions. Furthermore, we performed a phylogenomic analysis using whole plastomes of 40 taxa, representing all subfamilies of Arecaceae, which placed the macaw palm within the tribe Cocoseae. Finally, the data showed here are important for genetic studies in macaw palm and provide new insights into the evolution of plastid genes and environmental adaptation in Arecaceae.

Editing of Chloroplast rps14 by PPR Editing Factor EMB2261 Is Essential for Arabidopsis Development

RNA editing in plastids is known to be required for embryogenesis, but no single editing event had been shown to be essential. We show that the emb2261-2 mutation is lethal through a failure to express an editing factor that specifically recognizes the rps14-2 site. EMB2261 was predicted to bind the cis-element upstream of the rps14-2 site and genetic complementation with promoters of different strength followed by RNA-seq analysis was conducted to test the correlation between rps14-2 editing and EMB2261 expression. Rps14-2 is the only editing event in Arabidopsis chloroplasts that correlates with EMB2261 expression. Sequence divergence between the cis-element and the EMB2261 protein sequence in plants where rps14-2 editing is not required adds support to the association between them. We conclude that EMB2261 is the specificity factor for rps14-2 editing. This editing event converts P51 in Rps14 to L51, which is conserved among species lacking RNA editing, implying the importance of the editing event to Rps14 function. Rps14 is an essential ribosomal subunit for plastid translation, which, in turn, is essential for Arabidopsis embryogenesis.

Whole plastid transcriptomes reveal abundant RNA editing sites and differential editing status in Phalaenopsis aphrodite subsp. formosana

BACKGROUND

RNA editing is a process of post-transcriptional level of gene regulation by nucleotide modification. Previously, the chloroplast DNA of Taiwan endemic moth orchid, P. aphrodite subsp. formosana was determined, and 44 RNA editing sites were identified from 24 plastid protein-coding transcripts of leaf tissue via RT-PCR and then conventional Sanger sequencing. However, the RNA editing status of whole-plastid transcripts in leaf and other distinct tissue types in moth orchids has not been addressed. To sensitively and extensively examine the plastid RNA editing status of moth orchid, RNA-Seq was used to investigate the editing status of whole-plastid transcripts from leaf and floral tissues by mapping the sequence reads to the corresponding cpDNA template. With the threshold of at least 5% C-to-U or U-to-C conversion events observed in sequence reads considered as RNA editing sites.

RESULTS

In total, 137 edits with 126 C-to-U and 11 U-to-C conversions, including 93 newly discovered edits, were identified in plastid transcripts, representing an average of 0.09% of the nucleotides examined in moth orchid. Overall, 110 and 106 edits were present in leaf and floral tissues, respectively, with 79 edits in common. As well, 79 edits were involved in protein-coding transcripts, and the 58 nucleotide conversions caused the non-synonymous substitution. At least 32 edits showed significant (≧20%) differential editing between leaf and floral tissues. Finally, RNA editing in trnM is required for the formation of a standard clover-leaf structure.

CONCLUSIONS

We identified 137 edits in plastid transcripts of moth orchid, the highest number reported so far in monocots. The consequence of RNA editing in protein-coding transcripts mainly cause the amino acid change and tend to increase the hydrophobicity as well as conservation among plant phylogeny. RNA editing occurred in non-protein-coding transcripts such as tRNA, introns and untranslated regulatory regions could affect the formation and stability of secondary structure, which might play an important role in the regulation of gene expression. Furthermore, some unidentified tissue-specific factors might be required for regulating RNA editing in moth orchid.

A Single-Target Mitochondrial RNA Editing Factor of Funaria hygrometrica Can Fully Reconstitute RNA Editing at Two Sites in Physcomitrella patens

Nuclear-encoded pentatricopeptide repeat (PPR) proteins are key factors for site-specific RNA editing, converting cytidines into uridines in plant mitochondria and chloroplasts. All editing factors in the model moss Physcomitrella patens have a C-terminal DYW domain with similarity to cytidine deaminase. However, numerous editing factors in flowering plants lack such a terminal DYW domain, questioning its immediate role in the pyrimidine base conversion process. Here we further investigate the Physcomitrella DYW-type PPR protein PPR_78, responsible for mitochondrial editing sites cox1eU755SL and rps14eU137SL. Complementation assays with truncated proteins demonstrate that the DYW domain is essential for full PPR_78 editing functionality. The DYW domain can be replaced, however, with its counterpart from another editing factor, PPR_79. The PPR_78 ortholog of the related moss Funaria hygrometrica fully complements the Physcomitrella mutant for editing at both sites, although the editing site in rps14 is lacking in Funaria. Editing factor orthologs in different taxa may thus retain editing capacity for multiple sites despite the absence of editing requirement.

Chloroplast Genome Analysis of Resurrection Tertiary Relict Haberlea rhodopensis Highlights Genes Important for Desiccation Stress Response

Haberlea rhodopensis is a paleolithic tertiary relict species, best known as a resurrection plant with remarkable tolerance to desiccation. When exposed to severe drought stress, H. rhodopensis shows an ability to maintain the structural integrity of its photosynthetic apparatus, which re-activates easily upon rehydration. We present here the results from the assembly and annotation of the chloroplast (cp) genome of H. rhodopensis, which was further subjected to comparative analysis with the cp genomes of closely related species. H. rhodopensis showed a cp genome size of 153,099 bp, harboring a pair of inverted repeats (IR) of 25,415 bp separated by small and large copy regions (SSC and LSC) of 17,826 and 84,443 bp. The genome structure, gene order, GC content and codon usage are similar to those of the typical angiosperm cp genomes. The genome hosts 137 genes representing 70.66% of the plastome, which includes 86 protein-coding genes, 36 tRNAs, and 4 rRNAs. A comparative plastome analysis with other closely related Lamiales members revealed conserved gene order in the IR and LSC/SSC regions. A phylogenetic analysis based on protein-coding genes from 33 species defines this species as belonging to the Gesneriaceae family. From an evolutionary point of view, a site-specific selection analysis detected positively selected sites in 17 genes, most of which are involved in photosynthesis (e.g., rbcL, ndhF, accD, atpE, etc.). The observed codon substitutions may be interpreted as being a consequence of molecular adaptation to drought stress, which ensures an evolutionary advantage to H. rhodopensis.

RNA Editing and Its Molecular Mechanism in Plant Organelles

RNA editing by cytidine (C) to uridine (U) conversions is widespread in plant mitochondria and chloroplasts. In some plant taxa, “reverse” U-to-C editing also occurs. However, to date, no instance of RNA editing has yet been reported in green algae and the complex thalloid liverworts. RNA editing may have evolved in early land plants 450 million years ago. However, in some plant species, including the liverwort, Marchantia polymorpha, editing may have been lost during evolution. Most RNA editing events can restore the evolutionarily conserved amino acid residues in mRNAs or create translation start and stop codons. Therefore, RNA editing is an essential process to maintain genetic information at the RNA level. Individual RNA editing sites are recognized by plant-specific pentatricopeptide repeat (PPR) proteins that are encoded in the nuclear genome. These PPR proteins are characterized by repeat elements that bind specifically to RNA sequences upstream of target editing sites. In flowering plants, non-PPR proteins also participate in multiple RNA editing events as auxiliary factors. C-to-U editing can be explained by cytidine deamination. The proteins discovered to date are important factors for RNA editing but a bona fide RNA editing enzyme has yet to be identified.

The Unexpected Diversity of Plant Organelle RNA Editosomes

Flowering plants convert many hundreds of organelle cytidines (Cs) to uridines (Us) during post-transcriptional RNA editing. Pentatricopeptide repeat (PPR) proteins dictate specificity by recognizing RNA sequences near C targets. However, the complete mechanism of the editing machinery is not yet understood. Recently, non-PPR editing factors [RNA editing factor interacting proteins (RIPs)/multiple organellar RNA editing factors (MORFs), organelle RNA recognition motif (ORRM) proteins, organelle zinc-finger (OZ) proteins, and protoporphyrinogen oxidase 1 (PPO1)] have been identified as components of the plant RNA editosome, which is a small RNA-protein complex. Surprisingly, plant editosomes are highly diverse not only with regard to the PPR proteins they contain but also in the non-PPR components that are present. Here we review the most recent progress in the field and discuss the implications of the diversity of plant editosomes for the evolution of RNA editing and for possible future applications.

Genome-wide analysis of pentatricopeptide-repeat proteins of an aquatic plant

A large proportion of genes in plant genomes are organized as gene families. Whereas most gene families in the aquative plant Spirodela are reduced in their copy number, the PPR gene family is expanded, which match the RNA editing sites in organelles, providing us with new insights in the evolution of flowering plants. Pentatricopeptide-repeat proteins (PPRs) are nuclear-encoded proteins that are targeted to mitochondria and plastids to stabilize and edit mRNA transcribed from organellar genomes. They have been described for many terrestrial plant species from a diverse spectrum of sequenced genomes. To further increase our understanding of the evolution of this gene family across angiosperms, we analyzed the PPR genes in the aquatic species Spirodela polyrhiza in the order of the Alismatales (monocotyledonous plants). Because we had generated next generation sequencing data from transcripts and had sequenced the genome of Spirodela polyrhiza, we were able to identify its PPR genes and determine the level of their expression. In total, we could identify 556 PPR proteins, of which 238 members belong to the P (P motif) subfamily that is mainly involved in RNA stabilization and 318 ones to the PLS (P, Longer P, shorter P motif) subfamily responsible for RNA editing. Compared to other angiosperms, this is a large increase in the copy number of the PLS-PPRs subfamily and the expansion correlates with the increase of the number of RNA editing sites of organellar transcripts. Expression of PPR was generally stable even during growing and dormant stages, indicating that their function was critical throughout development. However, PPRs, especially those of the PLS subfamily, were expressed at relatively low levels, suggesting a delicate fine-tuning of its trans-acting function in the post-transcriptional regulation of gene expression. Thus, understanding PPR evolution and expression will help decipher the PPR code for their binding sites, which could genetically engineer RNA-binding proteins toward desired sequence.

Convergent Evolution of Fern-Specific Mitochondrial Group II Intron atp1i361g2 and Its Ancient Source Paralogue rps3i249g2 and Independent Losses of Intron and RNA Editing among Pteridaceae

Mitochondrial intron patterns are highly divergent between the major land plant clades. An intron in the atp1 gene, atp1i361g2, is an example for a group II intron specific to monilophytes (ferns). Here, we report that atp1i361g2 is lost independently at least 4 times in the fern family Pteridaceae. Such plant organelle intron losses have previously been found to be accompanied by loss of RNA editing sites in the flanking exon regions as a consequence of genomic recombination of mature cDNA. Instead, we now observe that RNA editing events in both directions of pyrimidine exchange (C-to-U and U-to-C) are retained in atp1 exons after loss of the intron in Pteris argyraea/biaurita and in Actiniopteris and Onychium We find that atp1i361g2 has significant similarity with intron rps3i249g2 present in lycophytes and gymnosperms, which we now also find highly conserved in ferns. We conclude that atp1i361g2 may have originated from the more ancestral rps3i249g2 paralogue by a reverse splicing copy event early in the evolution of monilophytes. Secondary structure elements of the two introns, most characteristically their domains III, show strikingly convergent evolution in the monilophytes. Moreover, the intron paralogue rps3i249g2 reveals relaxed evolution in taxa where the atp1i361g2 paralogue is lost. Our findings may reflect convergent evolution of the two related mitochondrial introns exerted by co-evolution with an intron-binding protein simultaneously acting on the two paralogues.

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

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

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.