A conserved glutamate residue in the C-terminal deaminase domain of pentatricopeptide repeat proteins is required for RNA editing activity

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.

Epigenetic changes induced by chronic and acute chromium stress treatments in Arabidopsis thaliana identified by the MSAP-Seq

Chromium (Cr) is an highly toxic metal to plants and causes severe damage to their growth, development, and reproduction. Plant exposure to chronic and acute Cr stress treatments results in significant changes at short time in the gene expression profile and at long time in the genomic DNA methylation profile at a transgenerational level and, consequently, in gene expression. These epigenetic modifications and their implications imposed by the Cr stress are not yet completely known in plants. Herein, were identified the epigenetic changes induced by chronic and acute Cr stress treatments in Arabidopsis thaliana plants using Methylation Sensitive Amplification Polymorphism coupled with next-generation sequencing (MSAP-Seq). First-generation Arabidopsis plants (termed F0 plants) kept under hoagland solution were subjected to Cr stress treatments. For chronic Cr stress, plants were treated through hoagland solution with 2.5 μM Cr during the entire cultivation period until seed harvest. Meanwhile, for acute Cr stress, plants were treated with 5 μM Cr during the first three weeks and returned to unstressful control condition until seed harvest. Seeds from F0 plants were sown and F1 plants were re-submitted to the same Cr stress treatments. The seed germination rate was evaluated from F-2 seeds harvested of F1 plants kept under different Cr stress treatments (0, 10, 20, and 40 μM) compared to the unstressful control condition. These data showed significant changes in the germination rate of F-2 seeds originating from stressed F1 plants compared to F-2 seeds harvested from unstressful control plants. Given this data, F1 plants kept under these chronic and acute Cr stress treatments and unstressful control condition were evaluated for the transgenerational epigenetic modifications using MSAP-Seq. The MSAP-Seq data showed that several genes were modified in their methylation status as a consequence of chronic and acute Cr stress treatment to maintain plant defenses activated. In particular, RNA processing, protein translation, photorespiration, energy production, transmembrane transport, DNA transcription, plant development, and plant resilience were the major biological processes modulated by epigenetic mechanisms identified in F1 plants kept under chronic and acute Cr stress. Therefore, collective data suggested that Arabidopsis plants kept under Cr stress regulate their epigenetic status over generations based on DNA methylation to modulate defense and resilience mechanisms.

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.

The Gating Domain of MEF28 Is Essential for Editing Two Contiguous Cytidines in nad2 mRNA in Arabidopsis thaliana

In plant organelles, each C-to-U RNA-editing site is specifically recognized by pentatricopeptide repeat (PPR) proteins with E1-E2, E1-E2-E+ or E1-E2-DYW domain extensions at the C-terminus. The distance between the PPR domain-binding site and the RNA-editing site is usually fixed at four bases, increasing the specificity of target-site recognition in this system. We here report, in contrast to the general case, on MEF28, which edits two adjacent mitochondrial sites, nad2-89 and nad2-90. When the sDYW domain of MEF28 was replaced with one derived from MEF11 or CRR22, the ability to edit downstream sites was lost, suggesting that the DYW domain of MEF28 provides unique target flexibility for two continuous cytidines. By contrast, substitutions of the entire E1-E2-DYW domains by MEF19E1-E2, SLO2E1-E2-E+ or CRR22E1-E2-E+ target both nad2 sites. In these cases, access to the contiguous sites in the chimeric PPR proteins is likely to be provided by the trans-associated DYW1-like proteins via the replaced E1-E2 or E1-E2-E+ domains. Furthermore, we demonstrated that the gating domain of MEF28 plays an important role in specific target-site recognition of the DYW domain. This finding suggests that the DYW domain and its internal gating domain fine-tune the specificity of the target site, which is valuable information for designing specific synthetic RNA-editing tools based on plant RNA-editing factors.

OsPGL3A encodes a DYW-type pentatricopeptide repeat protein involved in chloroplast RNA processing and regulated chloroplast development

The chloroplast serves as the primary site of photosynthesis, and its development plays a crucial role in regulating plant growth and morphogenesis. The Pentatricopeptide Repeat Sequence (PPR) proteins constitute a vast protein family that function in the post-transcriptional modification of RNA within plant organelles. In this study, we characterized mutant of rice with pale green leaves (pgl3a). The chlorophyll content of pgl3a at the seedling stage was significantly reduced compared to the wild type (WT). Transmission electron microscopy (TEM) and quantitative PCR analysis revealed that pgl3a exhibited aberrant chloroplast development compared to the wild type (WT), accompanied by significant alterations in gene expression levels associated with chloroplast development and photosynthesis. The Mutmap analysis revealed that a single base deletionin the coding region of Os03g0136700 in pgl3a. By employing CRISPR/Cas9 mediated gene editing, two homozygous cr-pgl3a mutants were generated and exhibited a similar phenotype to pgl3a, thereby confirming that Os03g0136700 was responsible for pgl3a. Consequently, it was designated as OsPGL3A. OsPGL3A belongs to the DYW-type PPR protein family and is localized in chloroplasts. Furthermore, we demonstrated that the RNA editing efficiency of rps8-182 and rpoC2-4106, and the splicing efficiency of ycf3-1 were significantly decreased in pgl3a mutants compared to WT. Collectively, these results indicate that OsPGL3A plays a crucial role in chloroplast development by regulating the editing and splicing of chloroplast genes in rice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01468-7.

The dicot homolog of maize PPR103 carries a C-terminal DYW domain and may have a role in C-to-U editing of some chloroplast RNA transcripts

In plants, cytidine-to-uridine (C-to-U) editing is a crucial step in processing mitochondria- and chloroplast-encoded transcripts. This editing requires nuclear-encoded proteins including members of the pentatricopeptide (PPR) family, especially PLS-type proteins carrying the DYW domain. IPI1/emb175/PPR103 is a nuclear gene encoding a PLS-type PPR protein essential for survival in Arabidopsis thaliana and maize. Arabidopsis IPI1 was identified as likely interacting with ISE2, a chloroplast-localized RNA helicase associated with C-to-U RNA editing in Arabidopsis and maize. Notably, while the Arabidopsis and Nicotiana IPI1 orthologs possess complete DYW motifs at their C-termini, the maize homolog, ZmPPR103, lacks this triplet of residues which are essential for editing. In this study we examined the function of IPI1 in chloroplast RNA processing in N. benthamiana to gain insight into the importance of the DYW domain to the function of the EMB175/PPR103/ IPI1 proteins. Structural predictions suggest that evolutionary loss of residues identified as critical for catalyzing C-to-U editing in other members of this class of proteins, were likely to lead to reduced or absent editing activity in the Nicotiana and Arabidopsis IPI1 orthologs. Virus-induced gene silencing of NbIPI1 led to defects in chloroplast ribosomal RNA processing and changes to stability of rpl16 transcripts, revealing conserved function with its maize ortholog. NbIPI1-silenced plants also had defective C-to-U RNA editing in several chloroplast transcripts, a contrast from the finding that maize PPR103 had no role in editing. The results indicate that in addition to its role in transcript stability, NbIPI1 may contribute to C-to-U editing in N. benthamiana chloroplasts.

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.

A ribonuclease activity linked to DYW1 in vitro is inhibited by RIP/MORF proteins

Organellar C-to-U RNA editing in plants occurs in complexes composed of various classes of nuclear-encoded proteins. DYW-deaminases are zinc metalloenzymes that catalyze hydrolytic deamination required for C-to-U modification editing. Solved crystal structures for DYW-deaminase domains display all structural features consistent with a canonical cytidine deamination mechanism. However, some recombinant DYW-deaminases from plants have been associated with ribonuclease activity in vitro. Direct ribonuclease activity by an editing factor is confounding since it is not required for deamination of cytosine, theoretically would be inimical for mRNA editing, and does not have a clear physiological function in vivo. His-tagged recombinant DYW1 from Arabidopsis thaliana (rAtDYW1) was expressed and purified using immobilized metal affinity chromatography (IMAC). Fluorescently labeled RNA oligonucleotides were incubated with recombinant AtDYW1 under different conditions. Percent relative cleavage of RNA probes was recorded at multiple time points from triplicate reactions. The effects of treatment with zinc chelators EDTA and 1, 10-phenanthroline were examined for rAtDYW1. Recombinant His-tagged RNA editing factors AtRIP2, ZmRIP9, AtRIP9, AtOZ1, AtCRR4, and AtORRM1 were expressed in E. coli and purified. Ribonuclease activity was assayed for rAtDYW1 in the presence of different editing factors. Lastly, the effects on nuclease activity in the presence of nucleotides and modified nucleosides were investigated. RNA cleavage observed in this study was linked to the recombinant editing factor rAtDYW1 in vitro. The cleavage reaction is sensitive to high concentrations of zinc chelators, indicating a role for zinc ions for activity. The addition of equal molar concentrations of recombinant RIP/MORF proteins reduced cleavage activity associated with rAtDYW1. However, addition of equal molar concentrations of purified recombinant editing complex proteins AtCRR4, AtORRM1, and AtOZ1 did not strongly inhibit ribonuclease activity on RNAs lacking an AtCRR4 cis-element. Though AtCRR4 inhibited AtDYW1 activity for oligonucleotides with a cognate cis-element. The observation that editing factors limit ribonuclease activity of rAtDYW1 in vitro, suggests that nuclease activities are limited to RNAs in absence of native editing complex partners. Purified rAtDYW1 was associated with the hydrolysis of RNA in vitro, and activity was specifically inhibited by RNA editing factors.

Structural insight into the activation of an Arabidopsis organellar C-to-U RNA editing enzyme by active site complementation

RNA-binding pentatricopeptide repeat (PPR) proteins catalyze hundreds of cytidine to uridine RNA editing events in plant organelles; these editing events are essential for proper gene expression. More than half of the PPR-type RNA editing factors, however, lack the DYW cytidine deaminase domain. Genetic analyses have suggested that their cytidine deaminase activity arises by association with a family of DYW1-like proteins that contain an N-terminally truncated DYW domain, but their molecular mechanism has been unclear. Here, we report the crystal structure of the Arabidopsis thaliana DYW1 deaminase domain at 1.8 Å resolution. DYW1 has a cytidine deaminase fold lacking the PG box. The internal insertion within the deaminase fold shows an α-helical fold instead of the β-finger reported for the gating domain of the A. thaliana ORGANELLE TRANSCRIPT PROCESSING 86. The substrate-binding pocket is incompletely formed and appears to be complemented in the complex by the E2 domain and the PG box of the interacting PPR protein. In vivo RNA editing assays corroborate the activation model for DYW1 deaminase. Our study demonstrates the common activation mechanism of the DYW1-like proteins by molecular complementation of the DYW domain and reconstitution of the substrate-binding pocket.

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.

Characterization of the Liriodendron chinense Pentatricopeptide Repeat (PPR) Gene Family and Its Role in Osmotic Stress Response

The Pentatricopeptide repeat (PPR) superfamily is a large gene family in plants that regulates organelle RNA metabolism, which is important for plant growth and development. However, a genome-wide analysis of the PPR gene family and its response to abiotic stress has not been reported for the relict woody plant Liriodendron chinense. In this paper, we identified 650 PPR genes from the L. chinense genome. A phylogenetic analysis showed that the LcPPR genes could roughly be divided into the P and PLS subfamilies. We found that 598 LcPPR genes were widely distributed across 19 chromosomes. An intraspecies synteny analysis indicated that duplicated genes from segmental duplication contributed to the expansion of the LcPPR gene family in the L. chinense genome. In addition, we verified the relative expression of Lchi03277, Lchi06624, Lchi18566, and Lchi23489 in the roots, stems, and leaves and found that all four genes had the highest expression in the leaves. By simulating a drought treatment and quantitative reverse transcription PCR (qRT-PCR) analysis, we confirmed the drought-responsive transcriptional changes in four LcPPR genes, two of which responded to drought stress independent of endogenous ABA biosynthesis. Thus, our study provides a comprehensive analysis of the L. chinense PPR gene family. It contributes to research into their roles in this valuable tree species' growth, development, and stress resistance.

The dicot homolog of maize PPR103 carries a C-terminal DYW domain and is required for C-to-U editing of chloroplast RNA transcripts

In plants, cytidine-to-uridine (C-to-U) editing is a crucial step in processing mitochondria and chloroplast-encoded transcripts. This editing requires nuclear-encoded proteins including members of the pentatricopeptide (PPR) family, especially PLS-type proteins carrying the DYW domain. IPI1/emb175/PPR103 is a nuclear gene encoding a PLS-type PPR protein essential for survival in Arabidopsis thaliana and maize. Arabidopsis IPI1 was identified as likely interacting with ISE2, a chloroplast-localized RNA helicase associated with C-to-U RNA editing in Arabidopsis and maize. Notably, while the Arabidopsis and Nicotiana IPI1 homologs possess complete DYW motifs at their C-termini, the maize homolog, ZmPPR103, lacks this triplet of residues which are essential for editing. We examined the function of ISE2 and IPI1 in chloroplast RNA processing in N. benthamiana. A combination of deep sequencing and Sanger sequencing revealed C-to-U editing at 41 sites in 18 transcripts, with 34 sites conserved in the closely related N. tabacum. Virus induced gene silencing of NbISE2 or NbIPI1 led to defective C-to-U revealed that they have overlapping roles at editing a site in the rpoB transcript but have distinct roles in editing other transcripts. This finding contrasts with maize ppr103 mutants that showed no defects in editing. The results indicate that NbISE2 and NbIPI1 are important for C-to-U editing in N. benthamiana chloroplasts, and they may function in a complex to edit specific sites while having antagonistic effects on editing others. That NbIPI1, carrying a DYW domain, is involved in organelle C-to-U RNA editing supports previous work showing that this domain catalyzes RNA editing.

SNP Detection in Pinus pinaster Transcriptome and Association with Resistance to Pinewood Nematode

Pinewood nematode (PWN, Bursaphelenchus xylophilus) is the causal agent of pine wilt disease (PWD), which severely affects Pinus pinaster stands in southwestern Europe. Despite the high susceptibility of P. pinaster, individuals of selected half-sib families have shown genetic variability in survival after PWN inoculation, indicating that breeding for resistance can be a valuable strategy to control PWD. In this work, RNA-seq data from susceptible and resistant plants inoculated with PWN were used for SNP discovery and analysis. A total of 186,506 SNPs were identified, of which 31 were highly differentiated between resistant and susceptible plants, including SNPs in genes involved in cell wall lignification, a process previously linked to PWN resistance. Fifteen of these SNPs were selected for validation through Sanger sequencing and 14 were validated. To evaluate SNP-phenotype associations, 40 half-sib plants were genotyped for six validated SNPs. Associations with phenotype after PWN inoculation were found for two SNPs in two different genes (MEE12 and PCMP-E91), as well as two haplotypes of HIPP41, although significance was not maintained following Bonferroni correction. SNPs here detected may be useful for the development of molecular markers for PWD resistance and should be further investigated in future association studies.

EMP80 mediates the C-to-U editing of nad7 and atp4 and interacts with ZmDYW2 in maize mitochondria

RNA C-to-U editing is important to the expression and function of organellar genes in plants. Although several families of proteins have been identified to participate in this process, the underlying mechanism is not fully understood. Here we report the function of EMP80 in the C-to-U editing at the nad7-769 and atp4-118 sites, and the potential recruitment of ZmDYW2 as a trans deaminase in maize (Zea mays) mitochondria. Loss of EMP80 function arrests embryogenesis and endosperm development in maize. EMP80 is a PPR-E+ protein localised to mitochondria. An absence of EMP80 abolishes the C-to-U RNA editing at nad7-769 and atp4-118 sites, resulting in a cysteine-to-arginine (Cys→Arg) change in Nad7 and Atp4 in the emp80 mutant. The amino acid change consequently reduces the assembly of complexes I and V, leading to an accumulation of the F₁ subcomplex of complex V. EMP80 was found to interact with atypical DYW-type PPR protein ZmDYW2, which interacts with ZmNUWA. Co-expression of ZmNUWA enhances the interaction between EMP80 and ZmDYW2, suggesting that EMP80 potentially recruits ZmDYW2 as a trans deaminase through protein-protein interaction, and ZmNUWA may function as an enhancer of this interaction.

OTP970 Is Required for RNA Editing of Chloroplast ndhB Transcripts in Arabidopsis thaliana

RNA editing is essential for compensating for defects or mutations in haploid organelle genomes and is regulated by numerous trans-factors. Pentatricopeptide repeat (PPR) proteins are the prime factors that are involved in RNA editing; however, many have not yet been identified. Here, we screened the plastid-targeted PLS-DYW subfamily of PPR proteins belonging to Arabidopsis thaliana and identified ORGANELLE TRANSCRIPT PROCESSING 970 (OTP970) as a key player in RNA editing in plastids. A loss-of-function otp970 mutant was impaired in RNA editing of ndhB transcripts at site 149 (ndhB-C149). RNA-immunoprecipitation analysis indicated that OTP970 was associated with the ndhB-C149 site. The complementation of the otp970 mutant with OTP970 lacking the DYW domain (OTP970^∆DYW) failed to restore the RNA editing of ndhB-C149. ndhB gene encodes the B subunit of the NADH dehydrogenase-like (NDH) complex; however, neither NDH activity and stability nor NDH-PSI supercomplex formation were affected in otp970 mutant compared to the wild type, indicating that alteration in amino acid sequence is not necessary for NdhB function. Together, these results suggest that OTP970 is involved in the RNA editing of ndhB-C149 and that the DYW domain is essential for its function.

Research Progress of PPR Proteins in RNA Editing, Stress Response, Plant Growth and Development

RNA editing is a posttranscriptional phenomenon that includes gene processing and modification at specific nucleotide sites. RNA editing mainly occurs in the genomes of mitochondria and chloroplasts in higher plants. In recent years, pentatricopeptide repeat (PPR) proteins, which may act as trans-acting factors of RNA editing have been identified, and the study of PPR proteins has become a research focus in molecular biology. The molecular functions of these proteins and their physiological roles throughout plant growth and development are widely studied. In this minireview, we summarize the current knowledge of the PPR family, hoping to provide some theoretical reference for future research and applications.

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.

Functions of PPR Proteins in Plant Growth and Development

Pentatricopeptide repeat (PPR) proteins form a large protein family in land plants, with hundreds of different members in angiosperms. In the last decade, a number of studies have shown that PPR proteins are sequence-specific RNA-binding proteins involved in multiple aspects of plant organellar RNA processing, and perform numerous functions in plants throughout their life cycle. Recently, computational and structural studies have provided new insights into the working mechanisms of PPR proteins in RNA recognition and cytidine deamination. In this review, we summarized the research progress on the functions of PPR proteins in plant growth and development, with a particular focus on their effects on cytoplasmic male sterility, stress responses, and seed development. We also documented the molecular mechanisms of PPR proteins in mediating RNA processing in plant mitochondria and chloroplasts.

GmPGL2, Encoding a Pentatricopeptide Repeat Protein, Is Essential for Chloroplast RNA Editing and Biogenesis in Soybean

Chloroplast biogenesis and development are highly complex processes requiring interactions between plastids and nuclear genomic products. Pentatricopeptide repeat (PPR) proteins play an essential role in the development of chloroplasts; however, it remains unclear how RNA editing factors influence soybean development. In this study, a Glycine max pale green leaf 2 mutant (Gmpgl2) was identified with decreased chlorophyll contents. Genetic mapping revealed that a single-nucleotide deletion at position 1949 bp in the Glyma.05g132700 gene in the Gmpgl2 mutant, resulting in a truncated GmPGL2 protein. The nuclear-encoded GmPGL2 is a PLS-type PPR protein that localizes to the chloroplasts. The C-to-U editing efficiencies of rps16, rps18, ndhB, ndhD, ndhE, and ndhF were reduced in the Gmpgl2 mutant. RNA electrophoresis mobility shift assay (REMSA) analysis further revealed that GmPGL2 binds to the immediate upstream sequences at RNA editing sites of rps16 and ndhB in vitro, respectively. In addition, GmPGL2 was found to interact with GmMORF8, GmMORF9, and GmORRM6. These results suggest that GmPGL2 participates in C-to-U RNA editing via the formation of a complex RNA editosome in soybean chloroplasts.

PPR-DYW Protein EMP17 Is Required for Mitochondrial RNA Editing, Complex III Biogenesis, and Seed Development in Maize

The conversion of cytidines to uridines (C-to-U) at specific sites in mitochondrial and plastid transcripts is a post-transcriptional processing event that is important to the expression of organellar genes. Pentatricopeptide repeat (PPR) proteins are involved in this process. In this study, we report the function of a previously uncharacterized PPR-DYW protein, Empty pericarp17 (EMP17), in the C-to-U editing and kernel development in maize. EMP17 is targeted to mitochondria. The loss-function of EMP17 arrests maize kernel development, abolishes the editing at ccmF _C -799 and nad2-677 sites, and reduces the editing at ccmF _C -906 and -966 sites. The absence of editing causes amino acid residue changes in CcmF_C-267 (Ser to Pro) and Nad2-226 (Phe to Ser), respectively. As CcmF_C functions in cytochrome c (Cytc) maturation, the amount of Cytc and Cytc ₁ protein is drastically reduced in emp17, suggesting that the CcmF_C-267 (Ser to Pro) change impairs the CcmF_C function. As a result, the assembly of complex III is strikingly decreased in emp17. In contrast, the assembly of complex I appears less affected, suggesting that the Nad2-226 (Phe to Ser) change may have less impact on Nad2 function. Together, these results indicate that EMP17 is required for the C-to-U editing at several sites in mitochondrial transcripts, complex III biogenesis, and seed development in maize.

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.

Accumulation of the RNA polymerase subunit RpoB depends on RNA editing by OsPPR16 and affects chloroplast development during early leaf development in rice

Plastid-encoded genes are coordinately transcribed by the nucleus-encoded RNA polymerase (NEP) and the plastid-encoded RNA polymerase (PEP). Resulting primary transcripts are frequently subject to RNA editing by cytidine-to-uridine conversions at specific sites. The physiological role of many editing events is largely unknown. Here, we have used the CRISPR/Cas9 technique in rice to knock out a member of the PLS-DYW subfamily of pentatricopeptide repeat (PPR) proteins. We found that OsPPR16 is responsible for a single editing event at position 545 in the chloroplast rpoB messenger RNA (mRNA), resulting in an amino acid change from serine to leucine in the β-subunit of the PEP. In striking contrast to loss-of-function mutations of the putative orthologue in Arabidopsis, which were reported to have no visible phenotype, knockout of OsPPR16 leads to impaired accumulation of RpoB, reduced expression of PEP-dependent genes, and a pale phenotype during early plant development. Thus, by editing the rpoB mRNA, OsPPR16 is required for faithful plastid transcription, which in turn is required for Chl synthesis and efficient chloroplast development. Our results provide new insights into the interconnection of the finely tuned regulatory mechanisms that operate at the transcriptional and post-transcriptional levels of plastid gene expression.

A Genome-Wide Analysis of the Pentatricopeptide Repeat (PPR) Gene Family and PPR-Derived Markers for Flesh Color in Watermelon (Citrullus lanatus)

Watermelon (Citrullus lanatus) is an economically important fruit crop grown for consumption of its large edible fruit flesh. Pentatricopeptide-repeat (PPR) encoding genes, one of the large gene families in plants, are important RNA-binding proteins involved in the regulation of plant growth and development by influencing the expression of organellar mRNA transcripts. However, systematic information regarding the PPR gene family in watermelon remains largely unknown. In this comprehensive study, we identified and characterized a total of 422 C. lanatus PPR (ClaPPR) genes in the watermelon genome. Most ClaPPRs were intronless and were mapped across 12 chromosomes. Phylogenetic analysis showed that ClaPPR proteins could be divided into P and PLS subfamilies. Gene duplication analysis suggested that 11 pairs of segmentally duplicated genes existed. In-silico expression pattern analysis demonstrated that ClaPPRs may participate in the regulation of fruit development and ripening processes. Genotyping of 70 lines using 4 single nucleotide polymorphisms (SNPs) from 4 ClaPPRs resulted in match rates of over 0.87 for each validated SNPs in correlation with the unique phenotypes of flesh color, and could be used in differentiating red, yellow, or orange watermelons in breeding programs. Our results provide significant insights for a comprehensive understanding of PPR genes and recommend further studies on their roles in watermelon fruit growth and ripening, which could be utilized for cultivar development of watermelon.

Identification of soybean trans-factors associated with plastid RNA editing sites

RNA editing is a posttranscriptional process that changes nucleotide sequences, among which cytosine-to-uracil by a deamination reaction can revert non-neutral codon mutations. Pentatricopeptide repeat (PPR) proteins comprise a family of RNA-binding proteins, with members acting as editing trans-factors that recognize specific RNA cis-elements and perform the deamination reaction. PPR proteins are classified into P and PLS subfamilies. In this work, we have designed RNA biotinylated probes based in soybean plastid RNA editing sites to perform trans-factor specific protein isolation. Soybean cis-elements from these three different RNA probes show differences in respect to other species. Pulldown samples were submitted to mass spectrometry for protein identification. Among detected proteins, five corresponded to PPR proteins. More than one PPR protein, with distinct functional domains, was pulled down with each one of the RNA probes. Comparison of the soybean PPR proteins to Arabidopsis allowed identification of the closest homologous. Differential gene expression analysis demonstrated that the PPR locus Glyma.02G174500 doubled its expression under salt stress, which correlates with the increase of its potential rps14 editing. The present study represents the first identification of RNA editing trans-factors in soybean. Data also indicated that potential multiple trans-factors should interact with RNA cis-elements to perform the RNA editing.

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.

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

Hornworts are crucial to understand the phylogeny of early land plants. The emergence of 'reverse' U-to-C RNA editing accompanying the widespread C-to-U RNA editing in plant chloroplasts and mitochondria may be a molecular synapomorphy of a hornwort-tracheophyte clade. C-to-U RNA editing is well understood after identification of many editing factors in models like Arabidopsis thaliana and Physcomitrella patens, but there is no plant model yet to investigate U-to-C RNA editing. The hornwort Anthoceros agrestis is now emerging as such a model system. We report on the assembly and analyses of the A. agrestis chloroplast and mitochondrial genomes, their transcriptomes and editomes, and a large nuclear gene family encoding pentatricopeptide repeat (PPR) proteins likely acting as RNA editing factors. Both organelles in A. agrestis feature high amounts of RNA editing, with altogether > 1100 sites of C-to-U and 1300 sites of U-to-C editing. The nuclear genome reveals > 1400 genes for PPR proteins with variable carboxyterminal DYW domains. We observe significant variants of the 'classic' DYW domain, in the meantime confirmed as the cytidine deaminase for C-to-U editing, and discuss the first attractive candidates for reverse editing factors given their excellent matches to U-to-C editing targets according to the PPR-RNA binding code.

Seed comparative genomics in three coffee species identify desiccation tolerance mechanisms in intermediate seeds

In contrast to desiccation-tolerant 'orthodox' seeds, so-called 'intermediate' seeds cannot survive complete drying and are short-lived. All species of the genus Coffea produce intermediate seeds, but they show a considerable variability in seed desiccation tolerance (DT), which may help to decipher the molecular basis of seed DT in plants. We performed a comparative transcriptome analysis of developing seeds in three coffee species with contrasting desiccation tolerance. Seeds of all species shared a major transcriptional switch during late maturation that governs a general slow-down of metabolism. However, numerous key stress-related genes, including those coding for the late embryogenesis abundant protein EM6 and the osmosensitive calcium channel ERD4, were up-regulated during DT acquisition in the two species with high seed DT, C. arabica and C. eugenioides. By contrast, we detected up-regulation of numerous genes involved in the metabolism, transport, and perception of auxin in C. canephora seeds with low DT. Moreover, species with high DT showed a stronger down-regulation of the mitochondrial machinery dedicated to the tricarboxylic acid cycle and oxidative phosphorylation. Accordingly, respiration measurements during seed dehydration demonstrated that intermediate seeds with the highest DT are better prepared to cease respiration and avoid oxidative stresses.

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.

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.

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.

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.

Stable native RIP9 complexes associate with C-to-U RNA editing activity, PPRs, RIPs, OZ1, ORRM1 and ISE2

The mitochondrial and chloroplast mRNAs of the majority of land plants are modified through cytidine to uridine (C-to-U) RNA editing. Previously, forward and reverse genetic screens demonstrated a requirement for pentatricopeptide repeat (PPR) proteins for RNA editing. Moreover, chloroplast editing factors OZ1, RIP2, RIP9 and ORRM1 were identified in co-immunoprecipitation (co-IP) experiments, albeit the minimal complex sufficient for editing activity was never deduced. The current study focuses on isolated, intact complexes that are capable of editing distinct sites. Peak editing activity for four sites was discovered in size-exclusion chromatography (SEC) fractions ≥ 670 kDa, while fractions estimated to be approximately 413 kDa exhibited the greatest ability to convert a substrate containing the editing site rps14 C80. RNA content peaked in the ≥ 670 kDa fraction. Treatment of active chloroplast extracts with RNase A abolished the relationship of editing activity with high-MW fractions, suggesting a structural RNA component in native complexes. By immunoblotting, RIP9, OTP86, OZ1 and ORRM1 were shown to be present in active gel filtration fractions, though OZ1 and ORRM1 were mainly found in low-MW inactive fractions. Active editing factor complexes were affinity-purified using anti-RIP9 antibodies, and orthologs to putative Arabidopsis thaliana RNA editing factor PPR proteins, RIP2, RIP9, RIP1, OZ1, ORRM1 and ISE2 were identified via mass spectrometry. Western blots from co-IP studies revealed the mutual association of OTP86 and OZ1 with native RIP9 complexes. Thus, RIP9 complexes were discovered to be highly associated with C-to-U RNA editing activity and other editing factors indicative of their critical role in vascular plant editosomes.

Empty Pericarp21 encodes a novel PPR-DYW protein that is required for mitochondrial RNA editing at multiple sites, complexes I and V biogenesis, and seed development in maize

C-to-U editing is an important event in post-transcriptional RNA processing, which converts a specific cytidine (C)-to-uridine (U) in transcripts of mitochondria and plastids. Typically, the pentatricopeptide repeat (PPR) protein, which specifies the target C residue by binding to its upstream sequence, is involved in the editing of one or a few sites. Here we report a novel PPR-DYW protein EMP21 that is associated with editing of 81 sites in maize. EMP21 is localized in mitochondria and loss of the EMP21 function severely inhibits the embryogenesis and endosperm development in maize. From a scan of 35 mitochondrial transcripts produced by the Emp21 loss-of-function mutant, the C-to-U editing was found to be abolished at five sites (nad7-77, atp1-1292, atp8-437, nad3-275 and rps4-870), while reduced at 76 sites in 21 transcripts. In most cases, the failure to editing resulted in the translation of an incorrect residue. In consequence, the mutant became deficient with respect to the assembly and activity of mitochondrial complexes I and V. As six of the decreased editing sites in emp21 overlap with the affected editing sites in emp5-1, and the editing efficiency at rpl16-458 showed a substantial reduction in the emp21-1 emp5-4 double mutant compared with the emp21-1 and emp5-4 single mutants, we explored their interaction. A yeast two hybrid assay suggested that EMP21 does not interact with EMP5, but both EMP21 and EMP5 interact with ZmMORF8. Together, these results indicate that EMP21 is a novel PPR-DYW protein required for the editing of ~17% of mitochondrial target Cs, and the editing process may involve an interaction between EMP21 and ZmMORF8 (and probably other proteins).

The E domain of CRR2 participates in sequence-specific recognition of RNA in plastids

Pentatricopeptide repeat (PPR) proteins are modular RNA-binding proteins involved in different aspects of RNA metabolism in organelles. PPR proteins of the PLS subclass often contain C-terminal domains that are important for their function, but the role of one of these domains, the E domain, is far from resolved. Here, we elucidate the role of the E domain in CRR2 in plastids. We identified a surprisingly large number of small RNAs that represent in vivo footprints of the Arabidopsis PLS-class PPR protein CRR2. An unexpectedly strong base conservation was found in the nucleotides aligned to the E domain. We used both in vitro and in vivo experiments to reveal the role of the E domain of CRR2. The E domain of CRR2 can be predictably altered to prefer different nucleotides in its RNA ligand, and position 5 of the E1-motif is biologically important for the PPR-RNA interaction. The 'code' of the E domain PPR motifs is different from that of P- and S-motifs. The findings presented here show that the E domain of CRR2 is involved in sequence-specific interaction with its RNA ligand and have implications for our ability to predict RNA targets for PLS-PPRs and their use as biotechnological tools to manipulate specific RNAs in vivo.

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.

Plant-type pentatricopeptide repeat proteins with a DYW domain drive C-to-U RNA editing in Escherichia coli

RNA editing converting cytidines into uridines is a hallmark of gene expression in land plant chloroplasts and mitochondria. Pentatricopeptide repeat (PPR) proteins have a key role in target recognition, but the functional editosome in the plant organelles has remained elusive. Here we show that individual Physcomitrella patens DYW-type PPR proteins alone can perform efficient C-to-U editing in Escherichia coli reproducing the moss mitochondrial editing. Single amino acid exchanges in the DYW domain abolish RNA editing, confirming it as the functional cytidine deaminase. The modification of RNA targets and the identification of numerous off-targets in the E. coli transcriptome reveal nucleotide identities critical for RNA recognition and cytidine conversion. The straightforward amenability of the new E. coli setup will accelerate future studies on RNA target recognition through PPRs, on the C-to-U editing deamination machinery and towards future establishment of transcript editing in other genetic systems.

Bastian Oldenkott et al. show that single moss pentatricopeptide repeat proteins with a DYW domain are sufficient to drive efficient C-to-U RNA editing in Escherichia coli. They demonstrate that the E.coli system is an easy to manipulate platform for future studies on RNA target recognition and C-to-U RNA editing.

Molecular and Functional Diversity of RNA Editing in Plant Mitochondria

RNA editing is a fundamental biochemical process relating to the modification of nucleotides in messenger RNAs of functional genes in cells. RNA editing leads to re-establishment of conserved amino acid residues for functional proteins in nuclei, chloroplasts, and mitochondria. Identification of RNA editing factors that contributes to target site recognition increases our understanding of RNA editing mechanisms. Significant progress has been made in recent years in RNA editing studies for both animal and plant cells. RNA editing in nuclei and mitochondria of animal cells and in chloroplast of plant cells has been extensively documented and reviewed. RNA editing has been also extensively documented on plant mitochondria. However, functional diversity of RNA editing factors in plant mitochondria is not overviewed. Here, we review the biological significance of RNA editing, recent progress on the molecular mechanisms of RNA editing process, and function diversity of editing factors in plant mitochondrial research. We will focus on: (1) pentatricopeptide repeat proteins in Arabidopsis and in crop plants; (2) the progress of RNA editing process in plant mitochondria; (3) RNA editing-related RNA splicing; (4) RNA editing associated flower development; (5) RNA editing modulated male sterile; (6) RNA editing-regulated cell signaling; and (7) RNA editing involving abiotic stress. Advances described in this review will be valuable in expanding our understanding in RNA editing. The diverse functions of RNA editing in plant mitochondria will shed light on the investigation of molecular mechanisms that underlies plant development and abiotic stress tolerance.

Genome-wide analysis of the rice PPR gene family and their expression profiles under different stress treatments

Background

Pentatricopeptide-repeat proteins (PPRs) are characterized by tandem arrays of a degenerate 35-amino-acid (PPR motifs), which can bind RNA strands and participate in post-transcription. PPR proteins family is one of the largest families in land plants and play important roles in organelle RNA metabolism and plant development. However, the functions of PPR genes involved in biotic and abiotic stresses of rice (Oryza sativa L.) remain largely unknown.

Results

In the present study, a comprehensive genome-wide analysis of PPR genes was performed. A total of 491 PPR genes were found in the rice genome, of which 246 PPR genes belong to the P subfamily, and 245 genes belong to the PLS subfamily. Gene structure analysis showed that most PPR genes lack intron. Chromosomal location analysis indicated that PPR genes were widely distributed in all 12 rice chromosomes. Phylogenetic relationship analysis revealed the distinct difference between the P and PLS subfamilies. Many PPR proteins are predicted to target chloroplasts or mitochondria, and a PPR protein (LOC_Os10g34310) was verified to localize in mitochondria. Furthermore, three PPR genes (LOC_Os03g17634,LOC_Os07g40820,LOC_Os04g51350) were verified as corresponding miRNA targets. The expression pattern analysis showed that many PPR genes could be induced under biotic and abiotic stresses. Finally, seven PPR genes were confirmed with their expression patterns under salinity or drought stress.

Conclusions

We found 491 PPR genes in the rice genome, and our genes structure analysis and syntenic analysis indicated that PPR genes might be derived from amplification by retro-transposition. The expression pattern present here suggested that PPR proteins have crucial roles in response to different abiotic stresses in rice. Taken together, our study provides a comprehensive analysis of the PPR gene family and will facilitate further studies on their roles in rice growth and development.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5088-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.

Novel DYW-type pentatricopeptide repeat (PPR) protein BLX controls mitochondrial RNA editing and splicing essential for early seed development of Arabidopsis

In plants, RNA editing is a post-transcriptional process that changes specific cytidine to uridine in both mitochondria and plastids. Most pentatricopeptide repeat (PPR) proteins are involved in organelle RNA editing by recognizing specific RNA sequences. We here report the functional characterization of a PPR protein from the DYW subclass, Baili Xi (BLX), which contains five PPR motifs and a DYW domain. BLX is essential for early seed development, as plants lacking the BLX gene was embryo lethal and the endosperm failed to initiate cellularization. BLX was highly expressed in the embryo and endosperm, and the BLX protein was specifically localized in mitochondria, which is essential for BLX function. We found that BLX was required for the efficient editing of 36 editing sites in mitochondria. Moreover, BLX was involved in the splicing regulation of the fourth intron of nad1 and the first intron of nad2. The loss of BLX function impaired the mitochondrial function and increased the reactive oxygen species (ROS) level. Genetic complementation with truncated variants of BLX revealed that, in addition to the DYW domain, only the fifth PPR motif was essential for BLX function. The upstream sequences of the BLX-targeted editing sites are not conserved, suggesting that BLX serves as a novel and major mitochondrial editing factor (MEF) via a new non-RNA-interacting manner. This finding provides new insights into how a DYW-type PPR protein with fewer PPR motifs regulates RNA editing in plants.

Genome-wide investigation of pentatricopeptide repeat gene family in poplar and their expression analysis in response to biotic and abiotic stresses

Pentatricopeptide repeat (PPR) proteins, which are characterized by tandem 30-40 amino acid sequence motifs, constitute of a large gene family in plants. Some PPR proteins have been identified to play important roles in organellar RNA metabolism and organ development in Arabidopsis and rice. However, functions of PPR genes in woody species remain largely unknown. Here, we identified and characterized a total of 626 PPR genes containing PPR motifs in the Populus trichocarpa genome. A comprehensive genome-wide analysis of the poplar PPR gene family was performed, including chromosomal location, phylogenetic relationships and gene duplication. Genome-wide transcriptomic analysis showed that 154 of the PtrPPR genes were induced by biotic and abiotic treatments, including Marssonina brunnea, salicylic acid (SA), methyl jasmonate (MeJA), mechanical wounding, cold and salinity stress. Quantitative RT-PCR analysis further investigated the expression profiles of 11 PtrPPR genes under different stresses. Our results contribute to a comprehensive understanding the roles of PPR proteins and provided an insight for improving the stress tolerance in poplar.

The Pentatricopeptide Repeat Protein MEF31 is Required for Editing at Site 581 of the Mitochondrial tatC Transcript and Indirectly Influences Editing at Site 586 of the Same Transcript

Pentatricopeptide repeat (PPR) proteins constitute the largest family of proteins in angiosperms, and most members are predicted to play roles in the maturation of organellar RNAs. Here we describe the novel mitochondrial editing factor 31 (MEF31), an E-PPR protein involved in editing at two close sites in the same transcript encoding subunit C of the twin-arginine translocation (tat) pathway. MEF31 is essential for editing at site tatC-581 and application of the recently proposed amino acid code for RNA recognition by PPR proteins supports the view that MEF31 directly targets this site by recognizing its cis sequence. In contrast, editing at site tatC-586 five nucleotides downstream is only partially affected in plants lacking MEF31, being restored to wild-type levels in complemented plants. Application of the amino acid code and analysis of individual RNA molecules for editing at sites 581 and 586 suggest that MEF31 does not directly target site tatC-586, and only indirectly influences editing at this site. It is likely that editing at site tatC-581 improves recognition of the site tatC-586 cis sequence by a second unknown PPR protein.

A protein with an unusually short PPR domain, MEF8, affects editing at over 60 Arabidopsis mitochondrial C targets of RNA editing

An RNA-seq approach was used to investigate the role of a PLS-subfamily pentatricopeptide repeat protein, Mitochondrial Editing Factor 8 (MEF8), on editing in Arabidopsis mitochondria and plastids. MEF8 has an intact DYW domain, but possesses an unusually short PLS repeat region of only five repeats. The MEF8 T-DNA insertion (mef8) line exhibited reduced editing at 38 mitochondrial editing sites and increased editing at 24 sites; therefore the absence of MEF8 affects 11% of the mitochondrial editome. Notably, 60% of the matR transcripts' sites showed a decrease of editing extent in the mef8 mutant. An E549A substitution in the MEF8 protein replaced the putatively catalytic glutamate of the HXE motif in the DYW domain. Complementation with MEF8-E549A failed to restore editing at the main target sites but was able to restore editing at the matR transcript; it also decreased the editing extent of most of the C targets exhibiting an increase of editing extent in the mef8 mutant plant. Thus, MEF8 has two antagonistic effects on mitochondrial editing: stimulatory, which requires a catalytic glutamate for most of the targets except for the matR transcript, and inhibitory, for which glutamate is dispensable.

Whole-transcriptome RNA-seq, gene set enrichment pathway analysis, and exon coverage analysis of two plastid RNA editing mutants

In land plants, plastid and mitochondrial RNAs are subject to post-transcriptional C-to-U RNA editing. T-DNA insertions in the ORGANELLE RNA RECOGNITION MOTIF PROTEIN6 gene resulted in reduced photosystem II (PSII) activity and smaller plant and leaf sizes. Exon coverage analysis of the ORRM6 gene showed that orrm6-1 and orrm6-2 are loss-of-function mutants. Compared to other ORRM proteins, ORRM6 affects a relative small number of RNA editing sites. Sanger sequencing of reverse transcription-PCR products of plastid transcripts revealed 2 plastid RNA editing sites that are substantially affected in the orrm6 mutants: psbF-C77 and accD-C794. The psbF gene encodes the β subunit of cytochrome b₅₅₉, an essential component of PSII. The accD gene encodes the β subunit of acetyl-CoA carboxylase, a protein required in plastid fatty acid biosynthesis. Whole-transcriptome RNA-seq demonstrated that editing at psbF-C77 is nearly absent and the editing extent at accD-C794 was significantly reduced. Gene set enrichment pathway analysis showed that expression of multiple gene sets involved in photosynthesis, especially photosynthetic electron transport, is significantly upregulated in both orrm6 mutants. The upregulation could be a mechanism to compensate for the reduced PSII electron transport rate in the orrm6 mutants. These results further demonstrated that Organelle RNA Recognition Motif protein ORRM6 is required in editing of specific RNAs in the Arabidopsis (Arabidopsis thaliana) plastid.

RNA editing of cytochrome c maturation transcripts is responsive to the energy status of leaf cells in Arabidopsis thaliana

Overexpression of AtPAP2, a phosphatase located on the outer membranes of chloroplasts and mitochondria, leads to higher energy outputs from these organelles. AtPAP2 interacts with seven MORF proteins of the editosome complex. RNA-sequencing analysis showed that the editing degrees of most sites did not differ significantly between OE and WT, except some sites on the transcripts of several cytochrome c maturation (Ccm) genes. Western blotting of 2D BN-PAGE showed that the patterns of CcmF_N1 polypeptides were different between the lines. We proposed that AtPAP2 may influence cytochrome c biogenesis by modulating RNA editing through its interaction with MORF proteins.

In Arabidopsis thaliana distinct alleles encoding mitochondrial RNA PROCESSING FACTOR 4 support the generation of additional 5' termini of ccmB transcripts

In plant mitochondria, the 5' ends of many transcripts are generated post-transcriptionally. We show that the pentatricopeptide repeat (PPR) protein RNA PROCESSING FACTOR 4 (RPF4) supports the generation of extra 5' ends of ccmB transcripts in Landsberg erecta (Ler) and a number of other Arabidopsis thaliana ecotypes. RPF4 was identified in Ler applying a forward genetic approach supported by complementation studies of ecotype Columbia (Col), which generates the Ler-type extra ccmB 5' termini only after the introduction of the RPF4 allele from Ler. Studies with chimeric RPF4 proteins composed of various parts of the RPF4 proteins from Ler and Col identified differences in the N-terminal and central PPR motifs that explain ecotype-specific variations in ccmB processing. These results fit well with binding site predictions in ccmB transcripts based on the known determinants of nucleotide base recognition by PPR motifs.

The pentatricopeptide repeat protein EMP9 is required for mitochondrial ccmB and rps4 transcript editing, mitochondrial complex biogenesis and seed development in maize

Pentatricopeptide repeat (PPR) proteins comprise a large family of sequence-specific RNA binding proteins in land plants. Because of its large family size and frequent embryo lethality in the mutants, molecular functions and physiological roles of many PPR proteins are unknown. Through characterization of an empty pericarp9 (emp9) mutant in maize (Zea mays), we defined the functions of EMP9 in mitochondrial RNA editing, respiratory complex formation and seed development. Mu insertions in different regions of Emp9 facilitated dissection of the domain functions of the EMP9. Through genetic and functional analyses of multiple alleles, we showed that deletions of two N-terminal PPR motifs and partial E+ domain do not eliminate the editing function of EMP9. Emp9 encodes an E+ subclass PPR protein that is localized in mitochondria. Loss of EMP9 function abolishes the C-to-U editing of ccmB-43 and rps4-335 sites in mitochondria. The loss of editing at ccmB-43 and rps4-335 affects the maturation of cytochrome c and impairs the biogenesis of mitochondrial respiratory complexes, particularly complex III. This work extends our understanding of PPR-E+ protein in editing function and seed development, and provides insights into the molecular function of mitochondrial CcmB protein in higher plants.

A genome-wide identification and analysis of the DYW-deaminase genes in the pentatricopeptide repeat gene family in cotton (Gossypium spp.)

The RNA editing occurring in plant organellar genomes mainly involves the change of cytidine to uridine. This process involves a deamination reaction, with cytidine deaminase as the catalyst. Pentatricopeptide repeat (PPR) proteins with a C-terminal DYW domain are reportedly associated with cytidine deamination, similar to members of the deaminase superfamily. PPR genes are involved in many cellular functions and biological processes including fertility restoration to cytoplasmic male sterility (CMS) in plants. In this study, we identified 227 and 211 DYW deaminase-coding PPR genes for the cultivated tetraploid cotton species G. hirsutum and G. barbadense (2n = 4x = 52), respectively, as well as 126 and 97 DYW deaminase-coding PPR genes in the ancestral diploid species G. raimondii and G. arboreum (2n = 26), respectively. The 227 G. hirsutum PPR genes were predicted to encode 52-2016 amino acids, 203 of which were mapped onto 26 chromosomes. Most DYW deaminase genes lacked introns, and their proteins were predicted to target the mitochondria or chloroplasts. Additionally, the DYW domain differed from the complete DYW deaminase domain, which contained part of the E domain and the entire E+ domain. The types and number of DYW tripeptides may have been influenced by evolutionary processes, with some tripeptides being lost. Furthermore, a gene ontology analysis revealed that DYW deaminase functions were mainly related to binding as well as hydrolase and transferase activities. The G. hirsutum DYW deaminase expression profiles varied among different cotton tissues and developmental stages, and no differentially expressed DYW deaminase-coding PPRs were directly associated with the male sterility and restoration in the CMS-D2 system. Our current study provides an important piece of information regarding the structural and evolutionary characteristics of Gossypium DYW-containing PPR genes coding for deaminases and will be useful for characterizing the DYW deaminase gene family in cotton biology and breeding.

Salt stress affects mRNA editing in soybean chloroplasts

Soybean, a crop known by its economic and nutritional importance, has been the subject of several studies that assess the impact and the effective plant responses to abiotic stresses. Salt stress is one of the main environmental stresses and negatively impacts crop growth and yield. In this work, the RNA editing process in the chloroplast of soybean plants was evaluated in response to a salt stress. Bioinformatics approach using sRNA and mRNA libraries were employed to detect specific sites showing differences in editing efficiency. RT-qPCR was used to measure editing efficiency at selected sites. We observed that transcripts of NDHA, NDHB, RPS14 and RPS16 genes presented differences in coverage and editing rates between control and salt-treated libraries. RT-qPCR assays demonstrated an increase in editing efficiency of selected genes. The salt stress enhanced the RNA editing process in transcripts, indicating responses to components of the electron transfer chain, photosystem and translation complexes. These increases can be a response to keep the homeostasis of chloroplast protein functions in response to salt stress.

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.