Towards a plant model for enigmatic U-to-C RNA editing: the organelle genomes, transcriptomes, editomes and candidate RNA editing factors in the hornwort Anthoceros agrestis

Crosslinking of base modified RNAs by synthetic DYW-KP base editors implicates an enzymatic lysine as the nitrogen donor for U-to-C RNA editing

Base editing mechanisms are being investigated as potential therapeutic tools to alleviate genetic diseases. Sequence specific C-to-U and A-to-I base editing tools are capable of altering RNA and DNA sequences and utilize a hydrolytic deamination mechanism requiring an active site zinc ion and a glutamate residue. In plant organelles, DYW-PG domain containing enzymes catalyze C-to-U edits and likely use the canonical deamination mechanism. Proteins developed from consensus sequences from the DYW-KP domain family catalyze what initially appeared to be U-to-C edits leading to this investigation into mechanistic insights into U-to-C editing. The synthetic DYW-KP enzyme KP6 was found sufficient for C-to-U editing activity stimulated by the addition of carboxylic acids in vitro. Despite addition of 14 different putative amine/amide donors in vitro U-to-C editing could not be observed. C-to-U editing was found not to be concomitant with U-to-C editing thus discounting a pyrimidine transaminase mechanism. RNAs containing base modifications were covalently crosslinked to KP6, KP2, and KP3 proteins. Mass spectrometry of purified KP2 and KP6 proteins revealed an additional mass of 319 Da. A U-to-C crosslinking mechanism was projected to explain the link between crosslinking, RNA base changes, and the ∼319 Da mass. In this model an enzymatic lysine attacks C4 of uridine to form a Schiff base RNA-protein conjugate. Sequenced RT-PCR products from the fern Ceratopteris richardii indicate U-to-C base edits do not preserve proteinaceous crosslinks in planta. Hydrolysis of a protonated Schiff base conjugate releasing cytidine is hypothesized to explain the completed pathway in plants.

Is RNA editing truly absent in the complex thalloid liverworts (Marchantiopsida)? Evidence of extensive RNA editing from Cyathodium cavernarum

RNA editing is a crucial modification in plants' organellar transcripts that converts cytidine to uridine (C-to-U; and sometimes uridine to cytidine) in RNA molecules. This post-transcriptional process is controlled by the PLS-class protein with a DYW domain, which belongs to the pentatricopeptide repeat (PPR) protein family. RNA editing is widespread in land plants; however, complex thalloid liverworts (Marchantiopsida) are the only group reported to lack both RNA editing and DYW-PPR protein. The liverwort Cyathodium cavernarum (Marchantiopsida, Cyathodiaceae), typically found in cave habitats, was newly found to have 129 C-to-U RNA editing sites in its chloroplast and 172 sites in its mitochondria. The Cyathodium genus, specifically C. cavernarum, has a large number of PPR editing factor genes, including 251 DYW-type PPR proteins. These DYW-type PPR proteins may be responsible for C-to-U RNA editing in C. cavernarum. Cyathodium cavernarum possesses both PPR DYW proteins and RNA editing. Our analysis suggests that the remarkable RNA editing capability of C. cavernarum may have been acquired alongside the emergence of DYW-type PPR editing factors. These findings provide insight into the evolutionary pattern of RNA editing in land plants.

Applications of Synthetic Pentatricopeptide Repeat Proteins

RNA-binding proteins play integral roles in the regulation of essential processes in cells and as such are attractive targets for engineering to manipulate gene expression at the RNA level. Expression of transcripts in chloroplasts and mitochondria is heavily regulated by pentatricopeptide repeat (PPR) proteins. The diverse roles of PPR proteins and their naturally modular architecture make them ideal candidates for engineering. Synthetic PPR proteins are showing great potential to become valuable tools for controlling the expression of plastid and mitochondrial transcripts. In this review, by 'synthetic', we mean both rationally modified natural PPR proteins and completely novel proteins designed using the principles learned from their natural counterparts. We focus on the many different applications of synthetic PPR proteins, covering both their use in basic research to learn more about protein-RNA interactions and their use to achieve specific outcomes in RNA processing and the control of gene expression. We describe the challenges associated with the design, construction and deployment of synthetic PPR proteins and provide perspectives on how they might be assembled and used in future biotechnology applications.

Assembly and comparative analysis of the complete mitochondrial genome of white towel gourd (Luffa cylindrica)

Mitochondria play an important role in the energy production of plant cells through independent genetic systems. This study has aimed to assemble and annotate the functions of the mitochondrial (mt) genome of Luffa cylindrica. The mt genome of L. cylindrica contained two chromosomes with lengths of 380,879 bp and 67,982 bp, respectively. Seventy-seven genes including 39 protein-coding genes, 34 tRNA genes, 3 rRNA genes, and 1 pseudogene, were identified. About 90.63% of the codons ended with A or U bases, and 98.63% of monomers contained A/T, which contributed to the high A/T content (55.91%) of the complete mt genome. Six genes (ATP8, CCMFC, NAD4, RPL10, RPL5 and RPS4) showed positive selection. Phylogenetic analysis indicates that L. cylindrica is closely related to L. acutangula. The present results provide the mt genome of L. cylindrica, which may facilitate possible genetic variation, evolutionary, and molecular breeding studies of L. cylindrica.

Insights into U-to-C RNA editing from the lycophyte Phylloglossum drummondii

The lycophyte Phylloglossum drummondii is the sole inhabitant of its genus in the Huperzioideae group and one of a small minority of plants which perform uridine to cytidine RNA editing. We assembled the P. drummondii chloroplast and mitochondrial genomes and used RNA sequence data to build a comprehensive profile of organellar RNA editing events. In addition to many C-to-U editing events in both organelles, we found just four U-to-C editing events in the mitochondrial transcripts cob, nad1, nad5 and rpl2. These events are conserved in related lycophytes in the genera Huperzia and Phlegmariurus. De novo transcriptomes for three of these lycophytes were assembled to search for putative U-to-C RNA editing enzymes. Four putative U-to-C editing factors could be matched to the four mitochondrial U-to-C editing sites. Due to the unusually few numbers of U-to-C RNA editing sites, P. drummondii and related lycophytes are useful models for studying this poorly understood mechanism.

Profiling of RNA editing events in plant organellar transcriptomes with high-throughput sequencing

RNA editing is a crucial post-transcriptional modification process in plant organellar RNA metabolism. rRNA removal-based total RNA-seq is one of the most common methods to study this event. However, the lack of commercial kits to remove rRNAs limits the usage of this method, especially for non-model plant species. DSN-seq is a transcriptome sequencing method utilizing duplex-specific nuclease (DSN) to degrade highly abundant cDNA species especially those from rRNAs while keeping the robustness of transcript levels of the majority of other mRNAs, and has not been applied to study RNA editing in plants before. In this study, we evaluated the capability of DSN-seq to reduce rRNA content and profile organellar RNA editing events in plants, as well we used commercial Ribo-off-seq and standard mRNA-seq as comparisons. Our results demonstrated that DSN-seq efficiently reduced rRNA content and enriched organellar transcriptomes in rice. With high sensitivity to RNA editing events, DSN-seq and Ribo-off-seq provided a more complete and accurate RNA editing profile of rice, which was further validated by Sanger sequencing. Furthermore, DSN-seq also demonstrated efficient organellar transcriptome enrichment and high sensitivity for profiling RNA editing events in Arabidopsis thaliana. Our study highlights the capability of rRNA removal-based total RNA-seq for profiling RNA editing events in plant organellar transcriptomes and also suggests DSN-seq as a widely accessible RNA editing profiling method for various plant species.

The newly assembled chloroplast genome of Aeluropus littoralis: molecular feature characterization and phylogenetic analysis with related species

Aeluropus littoralis, a halophyte grass, is widely distributed from the Mediterranean to the Indian subcontinent through the Mongolian Gobi. This model halophyte has garnered increasing attention owing to its use as forage and its high tolerance to environmental stressors. The chloroplast genomes of many plants have been extensively examined for molecular, phylogenetic and transplastomic applications. However, no published research on the A. littoralis chloroplast (cp) genome was discovered. Here, the entire chloroplast genome of A. littoralis was assembled implementing accurate long-read sequences. The entire chloroplast genome, with an estimated length of 135,532 bp (GC content: 38.2%), has a quadripartite architecture and includes a pair of inverted repeat (IR) regions, IRa and IRb (21,012 bp each), separated by a large and a small single-copy regions (80,823 and 12,685 bp, respectively). The features of A. littoralis consist of 133 genes that synthesize 87 peptides, 38 transfer RNAs, and 8 ribosomal RNAs. Of these genes, 86 were unique, whereas 19 were duplicated in IR regions. Additionally, a total of forty-six simple sequence repeats, categorized into 32-mono, four-di, two-tri, and eight-tetranucleotides, were discovered. Furthermore, ten sets of repeats greater than 20 bp were located primarily in the LSC region. Evolutionary analysis based on chloroplast sequence data revealed that A. littoralis with A. lagopoides and A. sinensis belong to the Aeluropodinae subtribe, which is a sister to the Eleusininae in the tribe Cynodonteae and the subfamily Chloridoideae. This subfamily belongs to the PACMAD clade, which contains the majority of the C4 photosynthetic plants in the Poaceae. The newly constructed A. littoralis cp genome offers valuable knowledge for DNA barcoding, phylogenetic, transplastomic research, and other biological studies.

The complete mitochondrial genome of Aglaia odorata, insights into its genomic structure and RNA editing sites

Aglaia odorata, native to Guangdong, Guangxi, and Hainan provinces in China, has long been utilized as an herbal remedy in ancient China. In this study, we assembled and annotated the complete mitochondrial genome (mitogenome) of A. odorata, which spans a total length of 537,321 bp. Conformation of the A. odorata recombination was verified through PCR experiments and Sanger sequencing. We identified and annotated 35 protein-coding genes (PCGs), 22 tRNA genes, and 3 rRNA genes within the mitogenome. Analysis of repeated elements revealed the presence of 192 SSRs, 29 pairs of tandem repeats, and 333 pairs of dispersed repeats in the A. odorata mitogenome. Additionally, we analyzed codon usage and mitochondrial plastid DNAs (MTPTs). Twelve MTPTs between the plastome and mitogenome of A. odorata were identified, with a combined length of 2,501 bp, accounting for 0.47% of the mitogenome. Furthermore, 359 high-confidence C to U RNA editing sites were predicted on PCGs, and four selected RNA editing sites were specially examined to verify the creation of start and/or stop codons. Extensive genomic rearrangement was observed between A. odorata and related mitogenomes. Phylogenetic analysis based on mitochondrial PCGs were conducted to elucidate the evolutionary relationships between A. odorata and other angiosperms.

The first two whole mitochondrial genomes for the genus Dactylis species: assembly and comparative genomics analysis

BACKGROUND

Orchardgrass (Dactylis glomerata L.), a perennial forage, has the advantages of rich leaves, high yield, and good quality and is one of the most significant forage for grassland animal husbandry and ecological management in southwest China. Mitochondrial (mt) genome is one of the major genetic systems in plants. Studying the mt genome of the genus Dactylis could provide more genetic information in addition to the nuclear genome project of the genus.

RESULTS

In this study, we sequenced and assembled two mitochondrial genomes of Dactylis species of D. glomerata (597, 281 bp) and D. aschersoniana (613, 769 bp), based on a combination of PacBio and Illumina. The gene content in the mitochondrial genome of D. aschersoniana is almost identical to the mitochondrial genome of D. glomerata, which contains 22-23 protein-coding genes (PCGs), 8 ribosomal RNAs (rRNAs) and 30 transfer RNAs (tRNAs), while D. glomerata lacks the gene encoding the Ribosomal protein (rps1) and D. aschersoniana contains one pseudo gene (atp8). Twenty-three introns were found among eight of the 30 protein-coding genes, and introns of three genes (nad 1, nad2, and nad5) were trans-spliced in Dactylis aschersoniana. Further, our mitochondrial genome characteristics investigation of the genus Dactylis included codon usage, sequences repeats, RNA editing and selective pressure. The results showed that a large number of short repetitive sequences existed in the mitochondrial genome of D. aschersoniana, the size variation of two mitochondrial genomes is due largely to the presence of a large number of short repetitive sequences. We also identified 52-53 large fragments that were transferred from the chloroplast genome to the mitochondrial genome, and found that the similarity was more than 70%. ML and BI methods used in phylogenetic analysis revealed that the evolutionary status of the genus Dactylis.

CONCLUSIONS

Thus, this study reveals the significant rearrangements in the mt genomes of Pooideae species. The sequenced Dactylis mt genome can provide more genetic information and improve our evolutionary understanding of the mt genomes of gramineous plants.

The first mitochondrial genome of Calophyllum soulattri Burm.f

Calophyllum soulattri Burm.f. is traditionally used to treat skin infections and reduce rheumatic pain, yet genetic and genomic studies are still limited. Here, we present the first complete mitochondrial genome of C. soulattri. It is 378,262 bp long with 43.97% GC content, containing 55 genes (30 protein-coding, 5 rRNA, and 20 tRNA). Repeat analysis of the mitochondrial genome revealed 194 SSRs, mostly mononucleotides, and 266 pairs of dispersed repeats ( ≥ 30 bp) that were predominantly palindromic. There were 23 homologous fragments found between the mitochondrial and plastome genomes. We also predicted 345 C-to-U RNA editing sites from 30 protein-coding genes (PCGs) of the C. soulatrii mitochondrial genome. These RNA editing events created the start codon of nad1 and the stop codon of ccmFc. Most PCGs of the C. soulattri mitochondrial genome underwent negative selection, but atp4 and ccmB experienced positive selection. Phylogenetic analyses showed C. soulattri is a sister taxon of Garcinia mangostana. This study has shed light on C. soulattri's evolution and Malpighiales' phylogeny. As the first complete mitochondrial genome in Calophyllaceae, it can be used as a reference genome for other medicinal plant species within the family for future genetic studies.

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.

Multichromosomal mitochondrial genome of Punica granatum: comparative evolutionary analysis and gene transformation from chloroplast genomes

BACKGROUND

Punica granatum is a fundamentally important fruit tree that has important economic, medicinal and ornamental properties. At present, there are few reports on the mitochondrial genome of pomegranate. Hence, in this study the P. granatum mitogenome was sequenced and assembled to further understanding of organization, variation, and evolution of mitogenomes of this tree species.

RESULTS

The genome structure was multi-chromosomes with seven circular contigs, measuring 382,774 bp in length with a 45.91% GC content. It contained 74 genes, including 46 protein-coding genes, 25 tRNA genes, and three rRNA genes. There were 188 pairs of dispersed repeats with lengths of 30 or greater, primarily consisting of reverse complementary repeats. The mitogenome analysis identified 114SSRs and 466 RNA editing sites. Analyses of codon usage, nucleotide diversity and gene migration from chloroplast to mitochondrial were also conducted. The collinear and comparative analysis of mitochondrial structures between P. granatum and its proximal species indicated that P. granatum 'Taishanhong' was closely related to P. granatum 'Qingpitian' and Lagerstroemia indica. Phylogenetic examination based on the mitogenome also confirmed the evolutionary relationship.

CONCLUSION

The results offered crucial information on the evolutionary biology of pomegranate and highlighted ways to promote the utilization of the species' germplasm.

Assembly and comparative analysis of the complete mitochondrial genome of Viburnum chinshanense

BACKGROUND

Viburnum chinshanense is an endemic species found exclusively in the North-Central and South-Central regions of China. This species is a lush garden ornamental tree and is extensively utilized for vegetation restoration in rocky desertification areas.

RESULTS

In this study, we obtained 13.96 Gb of Oxford Nanopore data for the whole genome, and subsequently, by combining Illumina short-reads, we successfully assembled the complete mitochondrial genome (mitogenome) of the V. chinshanense using a hybrid assembly strategy. The assembled genome can be described as a circular genome. The total length of the V. chinshanense mitogenome measures 643,971 bp, with a GC content of 46.18%. Our annotation efforts have revealed a total of 39 protein-coding genes (PCGs), 28 tRNA genes, and 3 rRNA genes within the V. chinshanense mitogenome. The analysis of repeated elements has identified 212 SSRs, 19 long tandem repeat elements, and 325 pairs of dispersed repeats in the V. chinshanense mitogenome. Additionally, we have investigated mitochondrial plastid DNAs (MTPTs) and identified 21 MTPTs within the mitogenome and plastidial genome. These MTPTs collectively span a length of 9,902 bp, accounting for 1.54% of the mitogenome. Moreover, employing Deepred-mt, we have confidently predicted 623 C to U RNA editing sites across the 39 protein-coding genes. Furthermore, extensive genomic rearrangements have been observed between V. chinshanense and the mitogenomes of related species. Interestingly, we have also identified a bacterial-derived tRNA gene (trnC-GCA) in the V. chinshanense mitogenome. Lastly, we have inferred the phylogenetic relationships of V. chinshanense with other angiosperms based on mitochondrial PCGs.

CONCLUSIONS

This study marks the first report of a mitogenome from the Viburnum genus, offering a valuable genomic resource for exploring the evolution of mitogenomes within the Dipsacales order.

Chloroplast gene expression: Recent advances and perspectives

Chloroplasts evolved from an ancient cyanobacterial endosymbiont more than 1.5 billion years ago. During subsequent coevolution with the nuclear genome, the chloroplast genome has remained independent, albeit strongly reduced, with its own transcriptional machinery and distinct features, such as chloroplast-specific innovations in gene expression and complicated post-transcriptional processing. Light activates the expression of chloroplast genes via mechanisms that optimize photosynthesis, minimize photodamage, and prioritize energy investments. Over the past few years, studies have moved from describing phases of chloroplast gene expression to exploring the underlying mechanisms. In this review, we focus on recent advances and emerging principles that govern chloroplast gene expression in land plants. We discuss engineering of pentatricopeptide repeat proteins and its biotechnological effects on chloroplast RNA research; new techniques for characterizing the molecular mechanisms of chloroplast gene expression; and important aspects of chloroplast gene expression for improving crop yield and stress tolerance. We also discuss biological and mechanistic questions that remain to be answered in the future.

Beyond a PPR-RNA recognition code: Many aspects matter for the multi-targeting properties of RNA editing factor PPR56

The mitochondrial C-to-U RNA editing factor PPR56 of the moss Physcomitrium patens is an RNA-binding pentatricopeptide repeat protein equipped with a terminal DYW-type cytidine deaminase domain. Transferred into Escherichia coli, PPR56 works faithfully on its two native RNA editing targets, nad3eU230SL and nad4eU272SL, and also converts cytidines into uridines at over 100 off-targets in the bacterial transcriptome. Accordingly, PPR56 is attractive for detailed mechanistic studies in the heterologous bacterial setup, allowing for scoring differential RNA editing activities of many target and protein variants in reasonable time. Here, we report (i) on the effects of numerous individual and combined PPR56 protein and target modifications, (ii) on the spectrum of off-target C-to-U editing in the bacterial background transcriptome for PPR56 and two variants engineered for target re-direction and (iii) on combinations of targets in tandem or separately at the 5'- and 3'-ends of large mRNAs. The latter experimentation finds enhancement of RNA editing at weak targets in many cases, including cox3eU290SF as a new candidate mitogenome target. We conclude that C-to-U RNA editing can be much enhanced by transcript features also outside the region ultimately targeted by PPRs of a plant editing factor, possibly facilitated by its enrichment or scanning along transcripts.

Comparative analysis of the chloroplast genomes of Rosa species and RNA editing analysis

BACKGROUND

The genus Rosa (Rosaceae) contains approximately 200 species, most of which have high ecological and economic values. Chloroplast genome sequences are important for studying species differentiation, phylogeny, and RNA editing.

RESULTS

In this study, the chloroplast genomes of three Rosa species, Rosa hybrida, Rosa acicularis, and Rosa rubiginosa, were assembled and compared with other reported Rosa chloroplast genomes. To investigate the RNA editing sites in R. hybrida (commercial rose cultivar), we mapped RNA-sequencing data to the chloroplast genome and analyzed their post-transcriptional features. Rosa chloroplast genomes presented a quadripartite structure and had highly conserved gene order and gene content. We identified four mutation hotspots (ycf3-trnS, trnT-trnL, psbE-petL, and ycf1) as candidate molecular markers for differentiation in the Rosa species. Additionally, 22 chloroplast genomic fragments with a total length of 6,192 bp and > 90% sequence similarity with their counterparts were identified in the mitochondrial genome, representing 3.96% of the chloroplast genome. Phylogenetic analysis including all sections and all subgenera revealed that the earliest divergence in the chloroplast phylogeny roughly distinguished species of sections Pimpinellifoliae and Rosa and subgenera Hulthemia. Moreover, DNA- and RNA-sequencing data revealed 19 RNA editing sites, including three synonymous and 16 nonsynonymous, in the chloroplast genome of R. hybrida that were distributed among 13 genes.

CONCLUSIONS

The genome structure and gene content of Rosa chloroplast genomes are similar across various species. Phylogenetic analysis based on the Rosa chloroplast genomes has high resolution. Additionally, a total of 19 RNA editing sites were validated by RNA-Seq mapping in R. hybrida. The results provide valuable information for RNA editing and evolutionary studies of Rosa and a basis for further studies on genomic breeding of Rosa species.

Plant organellar RNA maturation

Plant organellar RNA metabolism is run by a multitude of nucleus-encoded RNA-binding proteins (RBPs) that control RNA stability, processing, and degradation. In chloroplasts and mitochondria, these post-transcriptional processes are vital for the production of a small number of essential components of the photosynthetic and respiratory machinery-and consequently for organellar biogenesis and plant survival. Many organellar RBPs have been functionally assigned to individual steps in RNA maturation, often specific to selected transcripts. While the catalog of factors identified is ever-growing, our knowledge of how they achieve their functions mechanistically is far from complete. This review summarizes the current knowledge of plant organellar RNA metabolism taking an RBP-centric approach and focusing on mechanistic aspects of RBP functions and the kinetics of the processes they are involved in.

Assembly of the Complete Mitochondrial Genome of Pereskia aculeata Revealed That Two Pairs of Repetitive Elements Mediated the Recombination of the Genome

Pereskia aculeata is a potential new crop species that has both food and medicinal (antinociceptive activity) properties. However, comprehensive genomic research on P. aculeata is still lacking, particularly concerning its organelle genome. In this study, P. aculeata was studied to sequence the mitochondrial genome (mitogenome) and to ascertain the assembly, informational content, and developmental expression of the mitogenome. The findings revealed that the mitogenome of P. aculeata is circular and measures 515,187 bp in length with a GC content of 44.05%. It contains 52 unique genes, including 33 protein-coding genes, 19 tRNA genes, and three rRNA genes. Additionally, the mitogenome analysis identified 165 SSRs, primarily consisting of tetra-nucleotides, and 421 pairs of dispersed repeats with lengths greater than or equal to 30, which were mainly forward repeats. Based on long reads and PCR experiments, we confirmed that two pairs of long-fragment repetitive elements were highly involved with the mitogenome recombination process. Furthermore, there were 38 homologous fragments detected between the mitogenome and chloroplast genome, and the longest fragment was 3962 bp. This is the first report on the mitogenome in the family Cactaceae. The decoding of the mitogenome of P. aculeata will provide important genetic materials for phylogenetic studies of Cactaceae and promote the utilization of species germplasm resources.

Rickettsial DNA and a trans-splicing rRNA group I intron in the unorthodox mitogenome of the fern Haplopteris ensiformis

Plant mitochondrial genomes can be complex owing to highly recombinant structures, lack of gene syntenies, heavy RNA editing and invasion of chloroplast, nuclear or even foreign DNA by horizontal gene transfer (HGT). Leptosporangiate ferns remained the last major plant clade without an assembled mitogenome, likely owing to a demanding combination of the above. We here present both organelle genomes now for Haplopteris ensiformis. More than 1,400 events of C-to-U RNA editing and over 500 events of reverse U-to-C edits affect its organelle transcriptomes. The Haplopteris mtDNA is gene-rich, lacking only the ccm gene suite present in ancestral land plant mitogenomes, but is highly unorthodox, indicating extraordinary recombinogenic activity. Although eleven group II introns known in disrupted trans-splicing states in seed plants exist in conventional cis-arrangements, a particularly complex structure is found for the mitochondrial rrnL gene, which is split into two parts needing reassembly on RNA level by a trans-splicing group I intron. Aside from ca. 80 chloroplast DNA inserts that complicated the mitogenome assembly, the Haplopteris mtDNA features as an idiosyncrasy 30 variably degenerated protein coding regions from Rickettiales bacteria indicative of heavy bacterial HGT on top of tRNA genes of chlamydial origin.

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

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

Mitochondrial Genome Sequence of Salvia officinalis (Lamiales: Lamiaceae) Suggests Diverse Genome Structures in Cogeneric Species and Finds the Stop Gain of Genes through RNA Editing Events

Our previous study was the first to confirm that the predominant conformation of mitochondrial genome (mitogenome) sequence of Salvia species contains two circular chromosomes. To further understand the organization, variation, and evolution of Salvia mitogenomes, we characterized the mitogenome of Salvia officinalis. The mitogenome of S. officinalis was sequenced using Illumina short reads and Nanopore long reads and assembled using a hybrid assembly strategy. We found that the predominant conformation of the S. officinalis mitogenome also had two circular chromosomes that were 268,341 bp (MC1) and 39,827 bp (MC2) in length. The S. officinalis mitogenome encoded an angiosperm-typical set of 24 core genes, 9 variable genes, 3 rRNA genes, and 16 tRNA genes. We found many rearrangements of the Salvia mitogenome through inter- and intra-specific comparisons. A phylogenetic analysis of the coding sequences (CDs) of 26 common protein-coding genes (PCGs) of 11 Lamiales species and 2 outgroup taxa strongly indicated that the S. officinalis was a sister taxon to S. miltiorrhiza, consistent with the results obtained using concatenated CDs of common plastid genes. The mapping of RNA-seq data to the CDs of PCGs led to the identification of 451 C-to-U RNA editing sites from 31 PCGs of the S. officinalis mitogenome. Using PCR amplification and Sanger sequencing methods, we successfully validated 113 of the 126 RNA editing sites from 11 PCGs. The results of this study suggest that the predominant conformation of the S. officinalis mitogenome are two circular chromosomes, and the stop gain of rpl5 was found through RNA editing events of the Salvia mitogenome.

What can hornworts teach us?

The hornworts are a small group of land plants, consisting of only 11 families and approximately 220 species. Despite their small size as a group, their phylogenetic position and unique biology are of great importance. Hornworts, together with mosses and liverworts, form the monophyletic group of bryophytes that is sister to all other land plants (Tracheophytes). It is only recently that hornworts became amenable to experimental investigation with the establishment of Anthoceros agrestis as a model system. In this perspective, we summarize the recent advances in the development of A. agrestis as an experimental system and compare it with other plant model systems. We also discuss how A. agrestis can help to further research in comparative developmental studies across land plants and to solve key questions of plant biology associated with the colonization of the terrestrial environment. Finally, we explore the significance of A. agrestis in crop improvement and synthetic biology applications in general.

Lactobacillus for ribosome peptide editing cancer

This study reviews newly discovered insect peptide point mutations as new possible cancer research targets. To interpret newly discovered peptide point mutations in insects as new possible cancer research targets, we focused on the numerous peptide changes found in the 'CSP' family on the sex pheromone gland of the female silkworm moth Bombyx mori. We predict that the Bombyx peptide modifications will have a significant effect on cancer CUP (cancers of unknown primary) therapy and that bacterial peptide editing techniques, specifically Lactobacillus combined to CRISPR, will be used to regulate ribosomes and treat cancer in humans.

Assembly and comparative analysis of the complete mitochondrial genome of Salix wilsonii using PacBio HiFi sequencing

Salix L. (willows) is one of the most taxonomically complex genera of flowering plants, including shrubs, tall trees, bushes, and prostrate plants. Despite the high species diversity, only five mitochondrial genomes (mitogenomes) have been released in this genus. Salix wilsonii is an important ornamental and economic willow tree in section Wilsonia of the genus Salix. In this study, the S. wilsonii mitogenome was assembled into a typical circular structure with a size of 711,456 bp using PacBio HiFi sequencing. A total of 58 genes were annotated in the S. wilsonii mitogenome, including 33 protein-coding genes (PCGs), 22 tRNAs, and 3 rRNAs. In the S. wilsonii mitogenome, four genes (mttB, nad3, nad4, and sdh4) were found to play important roles in its evolution through selection pressure analysis. Collinearity analysis of six Salix mitogenomes revealed high structural variability. To determine the evolutionary position of S. wilsonii, we conducted a phylogenetic analysis of the mitogenomes of S. wilsonii and 12 other species in the order Malpighiales. Results strongly supported the segregation of S. wilsonii and other five Salix species with 100% bootstrap support. The comparative analysis of the S. wilsonii mitogenome not only sheds light on the functional and structural features of S. wilsonii but also provides essential information for genetic studies of the genus Salix.

A Comparative Genomics Approach for Analysis of Complete Mitogenomes of Five Actinidiaceae Plants

Actinidiaceae, an economically important plant family, includes the Actinidia, Clematoclethra and Saurauia genus. Kiwifruit, with remarkably high vitamin C content, is an endemic species widely distributed in China with high economic value. Although many Actinidiaceae chloroplast genomes have been reported, few complete mitogenomes of Actinidiaceae have been studied. Here, complete circular mitogenomes of the four kiwifruit species and Saurauia tristyla were assembled. Codon usage, sequence repeats, RNA editing, gene transfers, selective pressure, and phylogenetic relationships in the four kiwifruit species and S. tristyla were comparatively analyzed. This research will contribute to the study of phylogenetic relationships within Actiniaceae and molecular barcoding in kiwifruit.

U-to-C RNA editing by synthetic PPR-DYW proteins in bacteria and human culture cells

Programmable RNA editing offers significant therapeutic potential for a wide range of genetic diseases. Currently, several deaminase enzymes, including ADAR and APOBEC, can perform programmable adenosine-to-inosine or cytidine-to-uridine RNA correction. However, enzymes to perform guanosine-to-adenosine and uridine-to-cytidine (U-to-C) editing are still lacking to complete the set of transition reactions. It is believed that the DYW:KP proteins, specific to seedless plants, catalyze the U-to-C reactions in mitochondria and chloroplasts. In this study, we designed seven DYW:KP domains based on consensus sequences and fused them to a designer RNA-binding pentatricopeptide repeat (PPR) domain. We show that three of these PPR-DYW:KP proteins edit targeted uridine to cytidine in bacteria and human cells. In addition, we show that these proteins have a 5′ but not apparent 3′ preference for neighboring nucleotides. Our results establish the DYW:KP aminase domain as a potential candidate for the development of a U-to-C editing tool in human cells.

DYW:KP domains, designed on proteins found in the mitochondria and chloroplasts of seedless plants, are fused to a designer RNA-binding pentatricopeptide repeat (PPR) domain to edit targeted uridine to cytidine in bacteria and human cells.

Plant mitochondrial RNA editing factors can perform targeted C-to-U editing of nuclear transcripts in human cells

RNA editing processes are strikingly different in animals and plants. Up to thousands of specific cytidines are converted into uridines in plant chloroplasts and mitochondria whereas up to millions of adenosines are converted into inosines in animal nucleo-cytosolic RNAs. It is unknown whether these two different RNA editing machineries are mutually incompatible. RNA-binding pentatricopeptide repeat (PPR) proteins are the key factors of plant organelle cytidine-to-uridine RNA editing. The complete absence of PPR mediated editing of cytosolic RNAs might be due to a yet unknown barrier that prevents its activity in the cytosol. Here, we transferred two plant mitochondrial PPR-type editing factors into human cell lines to explore whether they could operate in the nucleo-cytosolic environment. PPR56 and PPR65 not only faithfully edited their native, co-transcribed targets but also different sets of off-targets in the human background transcriptome. More than 900 of such off-targets with editing efficiencies up to 91%, largely explained by known PPR-RNA binding properties, were identified for PPR56. Engineering two crucial amino acid positions in its PPR array led to predictable shifts in target recognition. We conclude that plant PPR editing factors can operate in the entirely different genetic environment of the human nucleo-cytosol and can be intentionally re-engineered towards new targets.

Comparative chloroplast genome and transcriptome analysis on the ancient genus Isoetes from China

Isoetes is a famous living fossil that plays a significant role in the evolutionary studies of the plant kingdom. To explore the adaptive evolution of the ancient genus Isoetes from China, we focused on Isoetes yunguiensis (Q.F. Wang and W.C. Taylor), I. shangrilaensis (X. Li, Y.Q. Huang, X.K. Dai & X. Liu), I. taiwanensis (DeVol), I. sinensis (T.C. Palmer), I. hypsophila_GHC (Handel-Mazzetti), and I. hypsophila_HZS in this study. We sequenced, assembled, and annotated six individuals' chloroplast genomes and transcriptomes, and performed a series of analyses to investigate their chloroplast genome structures, RNA editing events, and adaptive evolution. The six chloroplast genomes of Isoetes exhibited a typical quadripartite structure with conserved genome sequence and structure. Comparative analyses of Isoetes species demonstrated that the gene organization, genome size, and GC contents of the chloroplast genome are highly conserved across the genus. Besides, our positive selection analyses suggested that one positively selected gene was statistically supported in Isoetes chloroplast genomes using the likelihood ratio test (LRT) based on branch-site models. Moreover, we detected positive selection signals using transcriptome data, suggesting that nuclear-encoded genes involved in the adaption of Isoetes species to the extreme environment of the Qinghai-Tibetan Plateau (QTP). In addition, we identified 291-579 RNA editing sites in the chloroplast genomes of six Isoetes based on transcriptome data, well above the average of angiosperms. RNA editing in protein-coding transcripts results from amino acid changes to increase their hydrophobicity and conservation in Isoetes, which may help proteins form functional three-dimensional structure. Overall, the results of this study provide comprehensive transcriptome and chloroplast genome resources and contribute to a better understanding of adaptive evolutionary and molecular biology in Isoetes.

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.

DYW domain structures imply an unusual regulation principle in plant organellar RNA editing catalysis

RNA editosomes selectively deaminate cytidines to uridines in plant organellar transcripts–mostly to restore protein functionality and consequently facilitate mitochondrial and chloroplast function. The RNA editosomal pentatricopeptide repeat proteins serve target RNA recognition, whereas the intensively studied DYW domain elicits catalysis. Here we present structures and functional data of a DYW domain in an inactive ground state and activated. DYW domains harbour a cytidine deaminase fold and a C-terminal DYW motif, with catalytic and structural zinc atoms, respectively. A conserved gating domain within the deaminase fold regulates the active site sterically and mechanistically in a process that we termed gated zinc shutter. Based on the structures, an autoinhibited ground state and its activation are cross-validated by RNA editing assays and differential scanning fluorimetry. We anticipate that, in vivo, the framework of an active plant RNA editosome triggers the release of DYW autoinhibition to ensure a controlled and coordinated cytidine deamination playing a key role in mitochondrial and chloroplast homeostasis.

Organelle Genomes and Transcriptomes of Nymphaea Reveal the Interplay between Intron Splicing and RNA Editing

Posttranscriptional modifications, including intron splicing and RNA editing, are common processes during regulation of gene expression in plant organelle genomes. However, the intermediate products of intron-splicing, and the interplay between intron-splicing and RNA-editing were not well studied. Most organelle transcriptome analyses were based on the Illumina short reads which were unable to capture the full spectrum of transcript intermediates within an organelle. To fully investigate the intermediates during intron splicing and the underlying relationships with RNA editing, we used PacBio DNA-seq and Iso-seq, together with Illumina short reads genome and transcriptome sequencing data to assemble the chloroplast and mitochondrial genomes of Nymphaea 'Joey Tomocik' and analyze their posttranscriptional features. With the direct evidence from Iso-seq, multiple intermediates partially or fully intron-spliced were observed, and we also found that both cis- and trans-splicing introns were spliced randomly. Moreover, by using rRNA-depleted and non-Oligo(dT)-enrichment strand-specific RNA-seq data and combining direct SNP-calling and transcript-mapping methods, we identified 98 and 865 RNA-editing sites in the plastome and mitogenome of N. 'Joey Tomocik', respectively. The target codon preference, the tendency of increasing protein hydrophobicity, and the bias distribution of editing sites are similar in both organelles, suggesting their common evolutionary origin and shared editing machinery. The distribution of RNA editing sites also implies that the RNA editing sites in the intron and exon regions may splice synchronously, except those exonic sites adjacent to intron which could only be edited after being intron-spliced. Our study provides solid evidence for the multiple intermediates co-existing during intron-splicing and their interplay with RNA editing in organelle genomes of a basal angiosperm.

Deepred-Mt: Deep representation learning for predicting C-to-U RNA editing in plant mitochondria

In land plant mitochondria, C-to-U RNA editing converts cytidines into uridines at highly specific RNA positions called editing sites. This editing step is essential for the correct functioning of mitochondrial proteins. When using sequence homology information, edited positions can be computationally predicted with high precision. However, predictions based on the sequence contexts of such edited positions often result in lower precision, which is limiting further advances on novel genetic engineering techniques for RNA regulation. Here, a deep convolutional neural network called Deepred-Mt is proposed. It predicts C-to-U editing events based on the 40 nucleotides flanking a given cytidine. Unlike existing methods, Deepred-Mt was optimized by using editing extent information, novel strategies of data augmentation, and a large-scale training dataset, constructed with deep RNA sequencing data of 21 plant mitochondrial genomes. In comparison to predictive methods based on sequence homology, Deepred-Mt attains significantly better predictive performance, in terms of average precision as well as F1 score. In addition, our approach is able to recognize well-known sequence motifs linked to RNA editing, and shows that the local RNA structure surrounding editing sites may be a relevant factor regulating their editing. These results demonstrate that Deepred-Mt is an effective tool for predicting C-to-U RNA editing in plant mitochondria. Source code, datasets, and detailed use cases are freely available at https://github.com/aedera/deepredmt.

A synthetic RNA editing factor edits its target site in chloroplasts and bacteria

Members of the pentatricopeptide repeat (PPR) protein family act as specificity factors in C-to-U RNA editing. The expansion of the PPR superfamily in plants provides the sequence variation required for design of consensus-based RNA-binding proteins. We used this approach to design a synthetic RNA editing factor to target one of the sites in the Arabidopsis chloroplast transcriptome recognised by the natural editing factor CHLOROPLAST BIOGENESIS 19 (CLB19). We show that our synthetic editing factor specifically recognises the target sequence in in vitro binding assays. The designed factor is equally specific for the target rpoA site when expressed in chloroplasts and in the bacterium E. coli. This study serves as a successful pilot into the design and application of programmable RNA editing factors based on plant PPR proteins.

Genome-Wide Identification of U-To-C RNA Editing Events for Nuclear Genes in Arabidopsis thaliana

Cytosine-to-Uridine (C-to-U) RNA editing involves the deamination phenomenon, which is observed in animal nucleus and plant organelles; however, it has been considered the U-to-C is confined to the organelles of limited non-angiosperm plant species. Although previous RNA-seq-based analysis implied U-to-C RNA editing events in plant nuclear genes, it has not been broadly accepted due to inadequate confirmatory analyses. Here we examined the U-to-C RNA editing in Arabidopsis tissues at different developmental stages of growth. In this study, the high-throughput RNA sequencing (RNA-seq) of 12-day-old and 20-day-old Arabidopsis seedlings was performed, which enabled transcriptome-wide identification of RNA editing sites to analyze differentially expressed genes (DEGs) and nucleotide base conversions. The results showed that DEGs were expressed to higher levels in 12-day-old seedlings than in 20-day-old seedlings. Additionally, pentatricopeptide repeat (PPR) genes were also expressed at higher levels, as indicated by the log₂FC values. RNA-seq analysis of 12-day- and 20-day-old Arabidopsis seedlings revealed candidates of U-to-C RNA editing events. Sanger sequencing of both DNA and cDNA for all candidate nucleotide conversions confirmed the seven U-to-C RNA editing sites. This work clearly demonstrated presence of U-to-C RNA editing for nuclear genes in Arabidopsis, which provides the basis to study the mechanism as well as the functions of the unique post-transcriptional modification.

The hornworts: morphology, evolution and development

Extant land plants consist of two deeply divergent groups, tracheophytes and bryophytes, which shared a common ancestor some 500 million years ago. While information about vascular plants and the two of the three lineages of bryophytes, the mosses and liverworts, is steadily accumulating, the biology of hornworts remains poorly explored. Yet, as the sister group to liverworts and mosses, hornworts are critical in understanding the evolution of key land plant traits. Until recently, there was no hornwort model species amenable to systematic experimental investigation, which hampered detailed insight into the molecular biology and genetics of this unique group of land plants. The emerging hornwort model species, Anthoceros agrestis, is instrumental in our efforts to better understand not only hornwort biology but also fundamental questions of land plant evolution. To this end, here we provide an overview of hornwort biology and current research on the model plant A. agrestis to highlight its potential in answering key questions of land plant biology and evolution.

Deciphering the organelle genomes and transcriptomes of a common ornamental plant Ligustrum quihoui reveals multiple fragments of transposable elements in the mitogenome

Ligustrum quihoui (L. quihoui) is an important hedge material for landscaping and also possesses medicinal value. To generate genomic resources for better understanding the evolutionary history of this important plant, the organelle genomes of L. quihoui are de novo assembled and functionally annotated. Compared with other Oleaceae species, the 163,069 bp chloroplast genome of L. quihoui exhibits a typical quadripartite structure with highly conserved gene content and gene order, while the 848,451 bp mitochondrial genome of L. quihoui exhibits highly divergent genome size and gene content. Codon usage analyses show that genes related with photosynthesis and mitochondrial respiratory chain show inconsistent codon biases. A total of 48,760 bp transposable elements (TEs) fragments and 41,887 bp chloroplast-like sequences are found in the L. quihoui mitochondrial genome. A striking discrepancy of RNA editing between the two organelle genomes is found in L. quihoui, in which 146 mitochondrial editing sites coexist with only 43 such sites in chloroplast. Based on DNA and RNA-Seq data, we propose that GTG may act as the start codon of mitochondrial rpl16 in Oleaceae species. Phylogenetic analysis based on chloroplast genome shows that L. quihoui and L. japonicum form a sister clade within the genus Ligustrum.

Two Novel PLS-Class Pentatricopeptide Repeat Proteins Are Involved in the Group II Intron Splicing of Mitochondrial Transcripts in the Moss Physcomitrella patens

Pentatricopeptide repeat (PPR) proteins are RNA-binding proteins that function in posttranscriptional regulation as gene-specific regulators of RNA metabolism in plant organelles. Plant PPR proteins are divided into four classes: P, PLS, E and DYW. The E- and DYW-class proteins are mainly implicated in RNA editing, whereas most of the P-class proteins predominantly participate in RNA cleavage, splicing and stabilization. In contrast, the functions of PLS-class proteins still remain obscure. Here, we report the function of PLS-class PpPPR_31 and PpPPR_9 in Physcomitrella patens. The knockout (KO) mutants of PpPPR_31 and PpPPR_9 exhibited slower protonema growth compared to the wild type. The PpPPR_31 KO mutants showed a considerable reduction in the splicing of nad5 intron 3 and atp9 intron 1. The PpPPR_9 KO mutants displayed severely reduced splicing of cox1 intron 3. An RNA electrophoresis mobility shift assay showed that the recombinant PpPPR_31 protein bound to the 5' region of nad5 exon 4 and the bulged A region in domain VI of atp9 group II intron 1 while the recombinant PpPPR_9 bound to the translated region of ORF622 in cox1 intron 3. These results suggest that a certain set of PLS-class PPR proteins may influence the splicing efficiency of mitochondrial group II introns.

Extreme plastid RNA editing may confound phylogenetic reconstruction: A case study of Selaginella (lycophytes)

Cytidine-to-uridine (C-to-U) RNA editing is common in coding regions of organellar genomes throughout land plants. In most cases RNA editing alters translated amino acids or creates new start codons, potentially confounds phylogenetic reconstructions. In this study, we used the spike moss genus Selaginella (lycophytes), which has the highest frequency of RNA editing, as a model to test the effects of extreme RNA editing on phylogenetic reconstruction. We predicted the C-to-U RNA editing sites in coding regions of 18 Selaginella plastomes, and reconstructed the phylogenetic relationships within Selaginella based on three data set pairs consisted of plastome or RNA-edited coding sequences, first and second codon positions, and translated amino acid sequences, respectively. We predicted between 400 and 3100 RNA editing sites of 18 Selaginella plastomes. The numbers of RNA editing sites in plastomes were highly correlated with the GC content of first and second codon positions, but not correlated with the GC content of plastomes as a whole. Contrast phylogenetic analyses showed that there were substantial differences (e.g., the placement of clade B in Selaginella) between the phylogenies generated by the plastome and RNA-edited data sets. This empirical study provides evidence that extreme C-to-U RNA editing in the coding regions of organellar genomes alters the sequences used for phylogenetic reconstruction, and might even confound phylogenetic reconstruction. Therefore, RNA editing sites should be corrected when plastid or mitochondrial genes are used for phylogenetic studies, particularly in those lineages with abundant organellar RNA editing sites, such as hornworts, quillworts, spike mosses, and some seed plants.

Variation in protein gene and intron content among land plant mitogenomes

The functional content of the mitochondrial genome (mitogenome) is highly diverse across eukaryotes. Among land plants, our understanding of the variation in mitochondrial gene and intron content is improving from concerted efforts to densely sample mitogenomes from diverse land plants. Here I review the current state of knowledge regarding the diversity in content of protein genes and introns in the mitogenomes of all major land plant lineages. Mitochondrial protein gene content is largely conserved among mosses and liverworts, but it varies substantially among and within other land plant lineages due to convergent losses of genes encoding ribosomal proteins and, to a lesser extent, genes for proteins involved in cytochrome c maturation and oxidative phosphorylation. Mitochondrial intron content is fairly stable within each major land plant lineage, but highly variable among lineages, resulting from occasional gains and many convergent losses over time. Trans-splicing has evolved dozens of times in various vascular plant lineages, particularly those with relatively higher rates of mitogenomic rearrangement. Across eukaryotes, mitochondrial protein gene and intron content has been shaped massive convergent evolution.

A functional twintron, 'zombie' twintrons and a hypermobile group II intron invading itself in plant mitochondria

The occurrence of group II introns in plant mitochondrial genomes is strikingly different between the six major land plant clades, contrasting their highly conserved counterparts in chloroplast DNA. Their present distribution likely reflects numerous ancient intron gains and losses during early plant evolution before the emergence of seed plants. As a novelty for plant organelles, we here report on five cases of twintrons, introns-within-introns, in the mitogenomes of lycophytes and hornworts. An internal group II intron interrupts an intron-borne maturase of an atp9 intron in Lycopodiaceae, whose splicing precedes splicing of the external intron. An invasive, hypermobile group II intron in cox1, has conquered nine further locations including a previously overlooked sdh3 intron and, most surprisingly, also itself. In those cases, splicing of the external introns does not depend on splicing of the internal introns. Similar cases are identified in the mtDNAs of hornworts. Although disrupting a group I intron-encoded protein in one case, we could not detect splicing of the internal group II intron in this ‘mixed’ group I/II twintron. We suggest the name ‘zombie’ twintrons (half-dead, half-alive) for such cases where splicing of external introns does not depend any more on prior splicing of fossilized internal introns.

A plant pentatricopeptide repeat protein with a DYW-deaminase domain is sufficient for catalyzing C-to-U RNA editing in vitro

Pentatricopeptide repeat (PPR) proteins with C-terminal DYW domains are present in organisms that undergo C-to-U editing of organelle RNA transcripts. PPR domains act as specificity factors through electrostatic interactions between a pair of polar residues and the nitrogenous bases of an RNA target. DYW-deaminase domains act as the editing enzyme. Two moss (Physcomitrella patens) PPR proteins containing DYW-deaminase domains, PPR65 and PPR56, can convert Cs to Us in cognate, exogenous RNA targets co-expressed in Escherichia coli We show here that purified, recombinant PPR65 exhibits robust editase activity on synthetic RNAs containing cognate, mitochondrial PpccmFC sequences in vitro, indicating that a PPR protein with a DYW domain is solely sufficient for catalyzing C-to-U RNA editing in vitro Monomeric fractions possessed the highest conversion efficiency, and oligomeric fractions had reduced activity. Inductively coupled plasma (ICP)-MS analysis indicated a stoichiometry of two zinc ions per highly active PPR65 monomer. Editing activity was sensitive to addition of zinc acetate or the zinc chelators 1,10-o-phenanthroline and EDTA. Addition of ATP or nonhydrolyzable nucleotide analogs stimulated PPR65-catalyzed RNA-editing activity on PpccmFC substrates, indicating potential allosteric regulation of PPR65 by ATP. Unlike for bacterial cytidine deaminase, addition of two putative transition-state analogs, zebularine and tetrahydrouridine, failed to disrupt RNA-editing activity. RNA oligonucleotides with a single incorporated zebularine also did not disrupt editing in vitro, suggesting that PPR65 cannot bind modified bases due to differences in the structure of the active site compared with other zinc-dependent nucleotide deaminases.