# De novo RNA base editing in plant organelles with engineered synthetic P-type PPR editing factors

Sébastien Mathieu, Elena Lesch, Shahinez Garcia, Stéfanie Graindorge, Mareike Schallenberg-Rüdinger, Kamel Hammani

► To cite this version:

Sébastien Mathieu, Elena Lesch, Shahinez Garcia, Stéfanie Graindorge, Mareike Schallenberg-Rüdinger, et al. De novo RNA base editing in plant organelles with engineered synthetic P-type PPR editing factors. Nucleic Acids Research, 2025, 53 (18), pp.10297-10312. 10.1093/nar/gkaf279. hal-05191098

HAL Id: hal-05191098

https://hal.science/hal-05191098v1

Submitted on 29 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# De novo RNA base editing in plant organelles with engineered synthetic P-type PPR editing factors

Sébastien Mathieu<fnref n="1" />,†, Elena Lesch<fnref n="2" />,<fnref n="3" />,†, Shahinez Garcia<fnref n="1" />, Stéfanie Graindorge<fnref n="1" />, Mareike Schallenberg-Rüdinger ⑤, Kamel Hammani ⑥,<fnref n="1" />,*

## Abstract

In plant mitochondria and chloroplasts, cytidine-to-uridine RNA editing is necessary for the production of functional proteins. While natural PLS-type PPR proteins are specialized in this process, synthetic PPR proteins offer significant potential for targeted RNA editing. In this study, we engineered chimeric editing factors by fusing synthetic P-type PPR guides with the DYW cytidine deaminase domain of a moss mitochondrial editing factor, PPR56. These designer PPR editors (dPPRe) elicited efficient and precise *de novo* RNA editing in *Escherichia coli* as well as in the chloroplasts and mitochondria of *Nicotiana benthamiana*. Chloroplast transcriptome-wide analysis of the most efficient dPPRe revealed minimal off-target effects, with only three nontarget C sites edited due to sequence similarity with the intended target. This study introduces a novel and precise method for RNA base editing in plant organelles, paving the way for new approaches in gene regulation applicable to plants and potentially other organisms.

## Graphical abstract

<figure><img src="image_2.png" /></figure>

RNA editing is a post-transcriptional process that introduces specific nucleotide changes in RNA sequences, which is essential for the proper functioning of many eukaryotic cells [1, 2]. As a result, it can alter the genetic information encoded in the DNA, allowing for the synthesis of proteins with amino acid sequences that differ from those dictated by the genomic template [3]. The ability to reprogram RNA editing specificity offers significant potential for biotechnological applications

in basic research, medicine, and agriculture. It enables precise correction of genetic mutations, modulation of gene expression, regulation of protein function, and adaptation of cellular responses to various conditions.

In plants, RNA editing primarily involves the conversion of cytidines to uridines (C-to-U), and less frequently, uridines to cytidines (U-to-C) in organellar transcripts [4–9]. In angiosperms like *Arabidopsis thaliana*, over 30 cytidines in

<footnote n="1">¹Institut de Biologie Moléculaire des Plantes, Centre National de la Recherche Scientifique (CNRS), Université de Strasbourg, 12 rue du Général Zimmer, 67084 Strasbourg, France</footnote>

<footnote n="2">²Institut für Zelluläre und Molekulare Botanik, Abteilung Molekulare Evolution, Universität Bonn, 53115 Bonn, Germany</footnote>

<footnote n="3">³Present Address: Institut für Biologie und Biotechnologie der Pflanzen, Abteilung Grüne Biotechnologie, Universität Münster, 48143 Münster, Germany</footnote>

<footnote>*To whom correspondence should be addressed. Email: khammani@unistra.fr</footnote>

<footnote>†The first two authors should be regarded as Joint First Authors.</footnote>
chloroplast transcripts and >500 in mitochondrial transcripts are specifically edited [10, 11]. These modifications are crucial for the expression of functional proteins involved in essential metabolic and bioenergetic pathways [12–17]. The specificity and efficiency of RNA editing in plants is mediated by nuclear-encoded pentatricopeptide repeat (PPR) proteins. PPR proteins are a class of RNA-binding proteins composed of tandem helix-turn-helix repeats of ~35 amino acids [18–20]. They are specific to eukaryotes and besides RNA editing, they regulate various post-transcriptional steps in organelles, from RNA stabilization to translation [21]. In plants, PPR proteins can be subdivided into two subfamilies based on the type and arrangement of their PPR arrays [22, 23]. P-type PPR proteins contain long stretches of canonical 35 amino acid P-PPR motifs and are widespread in eukaryotes, while PLS-type PPR proteins consist of arrays of P, L (long, 35–36 aa), and S (short, 31–32 aa) motifs and are specific to land plants performing organellar RNA editing [23, 24]. Genetic studies have implicated P-type PPR proteins in various post-transcriptional processes such as RNA stabilization, RNA cleavage, intron splicing, and mRNA translation [21]. In contrast, PLS-type PPR proteins are the key players in RNA editing [8]. They bind precisely to an RNA sequence located in the cis-element 6 nucleotides upstream of the target C to guide the deamination reaction (Fig. 1A). Each PPR motif (P, L, and S) aligns with a ribonucleotide in the cis-element in a parallel orientation and the RNA base recognition is mainly determined by a combinatorial code involving the identity of two amino acids in each P- and S-type PPR motifs: amino acid in position five distinguishes purines (A, G) from pyrimidines (C, U) and the last amino acid defines preferences for amino (A, C) or keto ribonucleobases (G, U) [25–28]. The final P–L–S motif triplet (called P2–L2–S2) in PLS–PPR editing factors is often followed by E1 (for extension 1) and E2 motifs. The C-to-U conversion is catalyzed by a DYW deaminase domain of ~20 kDa, which is either naturally appended to the C-terminus of the PPR editing factor or recruited in trans [29–35]. Additionally, unusual SS-type PPR motifs were recently identified in several PLS-type PPR proteins of the lycophyte plant Selaginella and consists of tandem repeats of the S motif [22, 23] (Fig. 1A).

The modular architecture of PPR proteins and the elucidation of the PPR base recognition code have enabled the rational design of synthetic PPR proteins using consensus P, L, and S motifs with specific nucleotide binding capabilities [36–40]. *In vivo* applications of such synthetic P-PPR or PLS-PPR scaffolds for gene expression control in plants have been demonstrated in several contexts and include chloroplast RNA stabilization [41], translational activation [42], and RNA editing [39]. These studies have shown that synthetic PPR proteins can effectively mimic the functionality of their natural Arabidopsis counterparts in the respective mutant lines. Although these experiments provide proof-of-concept for synthetic PPR scaffolds being functional in chloroplasts, the potential to engineer these proteins for novel functions in gene expression by targeting new RNA sites in organelles has not yet been fully realized and requires further exploration. Additionally, the application of synthetic PPR scaffolds in mitochondria, which are associated with cytoplasmic male sterility in plants [43] and inherited myopathies in humans [44], remains largely unexplored.

In this study, we aimed to create *de novo* RNA editing in plant organelles by engineering chimeric editing factors. To this end, synthetic P-type PPR guides were fused to the DYW

cytidine deaminase domain of a moss mitochondrial editing factor, PPR56 [45]. We demonstrated that these designer PPR editors function independently of cofactors and successfully edit their RNA targets in *Escherichia coli*, as well as in the chloroplasts and mitochondria of *Nicotiana benthamiana* plants. Transcriptome-wide analysis of off-target effects of the most efficient dPPRe in *N. benthamiana* chloroplasts, which achieved an editing efficiency of up to ~70% at its designated target C site in *rpl2* gene, revealed minimal off-target activity, affecting only three C sites with conserved cis-elements similar to the target site. This study is the first to demonstrate precise *de novo* RNA editing in both plant organelles—chloroplasts and mitochondria—using engineered synthetic P-type PPR repeats. It also marks the first successful use of synthetic PPR in mitochondria, offering a novel strategy for regulating gene expression by targeting RNA effectors, with potential applications in both plant systems and other biological contexts.

## Materials and methods

Primers used in this study are listed in Supplementary Table S5 and DNA and protein sequences are provided in Supplementary Fig. S6.

### dPPRe transgene cloning

The DNA sequences encoding the Arabidopsis chloroplast RecA or yeast mitochondrial Cox4 transit peptide (68 and 29 amino acids, respectively) fused to a 3×HA tag and dPPRe were synthesized and cloned into pDONR™/Zeo by Thermo Fisher Scientific. For plant expression, the genes were subsequently cloned into binary vectors pMDC32 (2×35S promoter) [46] or pGWB2 (35S promoter) [47] using LR clonase II (Thermo Fisher Scientific, Waltham, MA, USA) following the manufacturer's protocol.

For expression in *E. coli*, the coding sequences of mature dPPRe without the organellar transit peptide and 3×HA tag were amplified from the pDONR™/Zeo templates. Target sequences comprising 100 bp up- and downstream of the cytidine to be edited, respectively, were amplified from *N. benthamiana* genomic DNA. dPPRe and RNA target encoding sequences were combined via Overlap Extension polymerase chain reaction (PCR) and inserted into petG41K_MCS using FastDigest restriction endonucleases KpnI and SgsI (Thermo Fisher Scientific).

### Bacterial RNA editing assays

*Escherichia coli* RNA editing assays were performed as outlined in Oldenkott *et al.* (2019). In short, constructs were transformed into *E. coli* Rosetta 2(DE3) competent cells (Novagen), which were grown in liquid culture until reaching an optical density (OD₆₀₀) of 0.4–0.6 before inducing recombinant dPPRe (rdPPRe) expression by the addition of 0.4 mM isopropyl β-d-1-thiogalactopyranoside (IPTG). Cells were harvested after 20 h of incubation at 16°C under continuous shaking. RNA editing was determined by reverse transcriptase-polymerase chain reaction (RT-PCR) and Sanger sequencing (Macrogen) and rdPPRe expression was verified by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS–PAGE) and Coomassie blue staining. The mean RNA editing rates (ERs) were calculated from three biological replicates (independent primary *E. coli* clones that underwent individual RNA editing assays). RNA editing frequencies in the
A

PLS-PPR editing factor

<figure><img src="image_3.png" /></figure>

SS-PPR editing factor

<figure><img src="image_4.png" /></figure>

B

dPPRe<sup>rp12</sup>

<figure><img src="image_5.png" /></figure>

dPPRe<sup>ndhB</sup>

<figure><img src="image_6.png" /></figure>

dPPRe<sup>nad7</sup>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

Shen et al., 2015

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

Royan et al., 2021

**Figure 1.** Design of precision RNA base PPR editors. (A) Motif composition and arrangement of natural PLS- and SS-type PPR editing factors in plant organelles. PPR motifs predicted to bind RNA bases are indicated by dotted black lines. Question marks indicate uncertain contribution to binding. The target cytidine is shown in red, and the designated PPR binding site ('cis-element') is highlighted in bold. The target cytidine is often preceded by a uridine in plant organellar genomes. TP: organellar transit peptide. (B) Engineering of chloroplastic and mitochondrial designer PPR editors (dPPRe) using arrays of synthetic P-type PPR motifs fused to the DYW cytidine deaminase domain from the mitochondrial editing factor PPR56. The specificity-determining amino acids at positions 5 and last of each PPR motif are indicated. The exact genomic position of the target cytidine in the *N. benthamiana* *rpl2* and *ndhB* chloroplast genes are shown beneath the RNA sequences, whereas the position within mitochondrial *nad7* gene was inferred from the mitochondrial genome sequence of *N. tabacum*. The codon containing the target cytidine is circled by a dotted box. CTP: chloroplast transit peptide, MTP: mitochondrial transit peptide. HA: hemagglutinin tag. NTD: N-terminal domain of PPR10 (amino acids 37–208).

complementary DNA (cDNA) were quantified through the MacVector application by measuring the heights of overlapping T and C peaks in Sanger chromatograms at the position of the target cytidine.

Transient expression in *N. benthamiana*

For transient protein expression, 6-week *N. benthamiana* plants were infiltrated with *Agrobacterium tumefaciens* GV3101 carrying the *dPPRe* transgenes. Agrobacterium cultures were grown overnight, resuspended in infiltration buffer (10 mM MgCl₂, 10 mM 2-(N-morpholino)ethanesulfonic acid, pH 5.6, and 100 µM acetosyringone) to an OD₆₀₀ of 0.7, and infiltrated into the abaxial side of fully expanded *N. benthamiana* leaves using a needleless syringe. A Tomato bushy stunt virus-encoded 35S:*P19-HA* silencing suppressor gene construct [48] (P19) was co-expressed with *dPPRe* transgenes to increase protein expression. Mixtures of P19/dPPRe agrobacteria were prepared in a 1:1 (v/v) ratio as required. P19 alone was used as the control. Infiltrated plants were maintained under controlled growth conditions and leaf discs were harvested at 2- and 3-day post-infiltration (dpi) for protein analysis by western blot using monoclonal anti-HA antibody (H9658 clone) from Sigma-Aldrich. Entire leaves were also collected at the same time points for RNA extraction using the TRIzol reagent, following the manufacturer's protocol, to assess RNA editing and gene expression levels by Sanger sequencing of RT-PCR products or RNA-seq. RNA editing frequencies were quantified as previously described or by GATK analysis as later described.

Plant RNA sequencing

Total RNA was extracted from *N. benthamiana* leaves infiltrated with dPPRe<sup>rp12</sup> and P19 constructs or P19 alone at 2-dpi or 3-dpi by Trizol extraction followed by phenol/chloroform extraction and ethanol precipitation. RNA was treated by DNase Turbo (Thermo Fisher Scientific) to remove contaminant DNA. Library preparation (ribosomal RNA-depleted) and 2×150-bp pair-ended Illumina sequencing were done by Novogene (https://www.novogene.com/), generating ~105–138 mio reads per library. Libraries were prepared and sequenced for three independent biological replicates of P19, dPPRe<sup>rp12</sup> (2×35S)/P19 and dPPRe<sup>rp12</sup> (35S)/P19 transient expression.

The quality of the reads was assessed with FastQC 0.12.1 and adapters were trimmed with Trim-galore 0.6.10. Sequences with a quality lower than 30, with >5 'N's or shorter than 20 bp after trimming were removed. The remaining reads were mapped against the *N. benthamiana* chloroplast genome SRR7540368 [49] (https://github.com/WangSoybean/Nicotiana_chloroplast_genomes) combined with the *dPPRe<sup>rp12</sup>* nucleotide sequence with HISAT2 2.2.1 using the -k 50 option and default parameters.
Mapped reads were counted with featureCounts 2.0.1 with
options -O -M -minOverlap 5, to consider multi-mapped
reads. Differential expression analysis was done using the
DESeq2 1.42.0 package.

Chloroplast RNA editing sites were identified from
the HISAT2 aligned reads with GATK 4.5.0.0 accord-
ing to the GATK best practice recommendations [50, 51]
(https://gencore.bio.nyu.edu/variant-calling-pipeline-gatk4/).
Initially, MarkDuplicatesSpark, HaplotypeCaller, and Select-
Variants were run with default parameters, and VariantFiltra-
tion with options "QD_filter" -filter "QD < 2.0", -filter-name
"FS_filter" -filter "FS > 60.0", -filter-name "MQ_filter" -filter
"MQ < 40.0", -filter-name "SOR_filter" -filter "SOR > 4.0",
-filter-name "MQRankSum_filter" -filter "MQRankSum <
-12.5", -filter-name "ReadPosRankSum_filter", -filter "Read-
PosRankSum < -8.0". This initial run served as input for
BaseRecalibrator, ApplyBQSR, and BaseRecalibrator steps.
These steps were then repeated with the same options to
call single nucleotide polymorphisms (SNPs) between the
chloroplast transcriptome and genome. RNA editing sites
were identified by examining C-to-T or G-to-A conversions.
*De novo* RNA editing sites induced by dPPRe<sup>rp12</sup>expression
were identified by detecting base conversions specifically
present in experimental libraries (dPPRe<sup>rp12</sup>/P19) but absent
in control libraries (P19). The editing sites introduced by the
expression of dPPRe<sup>rp12</sup> were considered significant if they
were detected in at least two independent libraries. ERs were
calculated as the ratio between the number of SNP reads
and read depth at the specific base position. The spearman
correlation test (cor.test in R 4.2.0) was used to assess the
correlation between the ER and expression level of dPPRe<sup>rp12</sup>.

The MA-plot and correlation test graphs were generated
using the ggplot2 3.5.1 R package. Secondary structures were
predicted with the mfold server (RNA folding form version
2.3) at http://www.unafold.org/ using default parameters and
a folding temperature of 24°C.

Expression of rdPPRe

The *E. coli* recombinant dPPRe (rdPPRe) expressing vec-
tors were cloned in one step by DNA assembly of overlap-
ping DNA fragments using the NEBuilder HiFi DNA As-
sembly Master Mix following the manufacturer's instruc-
tions (New England Biolabs). The nucleotide fragments cod-
ing for the predicted mature dPPRe (i.e. lacking the tran-
sit peptide) or the pMAL vector were amplified using Q5
High-Fidelity DNA polymerase (New England Biolabs) with
primers k1504/k1505 and k1502/k1503, respectively. The fi-
nal vectors express mature dPPRe<sup>rp12</sup> and dPPRe<sup>ndhB</sup> fused
in frame to the N-terminal maltose binding protein tag
(MBP) and were transformed into BL21 (DE3) pLysS cells.
The expression of rdPPRe was induced overnight at 16°C
by adding 0.5 mM IPTG, with constant agitation at 220
rpm. Bacterial cells were then lysed in a cold buffer con-
taining 30 mM Tris-HCl, pH 7.5, 400 mM NaCl, 10%
glycerol, 0.05% 3-[(3-cholamidopropyl)dimethylammonio]-
1-propanesulfonate (CHAPS), 5 mM β-mercaptoethanol, and
an ethylenediaminetetraacetic acid-free protease inhibitor
cocktail (Roche). Soluble rdPPRe was first affinity-purified us-
ing an amylose column (New England Biolabs) and subse-
quently resolved on a Superdex S200 column (GE Health-
care Life Sciences) in lysis buffer lacking protease inhibitors.
Peak fractions containing rdPPRe were pooled and dialyzed

overnight in elution buffer containing 50% glycerol. The pu-
rity of the tandem affinity-purified protein was assessed by
SDS-PAGE followed by Coomassie Brilliant Blue staining. Re-
combinant proteins were stored at -20°C and used within 10
days. Recombinant maltose-binding protein (rMBP) was pu-
rified following TEV digestion of rHCF107 as described [52].

Gel mobility shift assays

Synthetic RNAs (Integrated DNA Technologies) were puri-
fied on a denaturing polyacrylamide gel and 5'-end-labeled
with [γ-32P]-ATP and T4 polynucleotide kinase. Unincorpo-
rated radio-labelled nucleotides were removed by filtration on
illustra Microspin G-25 columns (GE Healthcare) followed
by phenol-chloroform extraction and ethanol precipitation.
Binding reactions contained 150 mM NaCl, 30 mM Tris-HCl,
pH 7.5, 4 mM Dithiothreitol, 0.04 mg/ml bovine serum al-
bumin, 0.5 mg/ml heparin, 10% glycerol, 0.01% CHAPS, 10
units RNaseOUT (Invitrogen), and 15 pM radiolabeled RNA.
Reactions were incubated for 60 min at 25°C and resolved on
8% native polyacrylamide gels. Results were visualized on a
Typhoon phosphorimager. Data quantification was performed
with ImageQuant TL software (Amersham).

Statistics and reproducibility

All experiments were repeated at least three times. The num-
ber of samples for each experience is provided in the fig-
ures. The statistical significance of the RNA ER differences
between experimental and control groups or between ex-
periments in Fig. 3C was assessed by *t*-test using GraphPad
Prism 8.4.0 (*t*-values and degrees of freedom can be found in
*Supplementary Table S2*).

Results

Engineering synthetic *dPPRe* for precise RNA base editing

To date, two synthetic editing factors with PLS-type or SS-type
PPR protein arrangements have been engineered to edit the
chloroplast *rpoA*-78,691 site *in vivo*, a site naturally targeted
by the endogenous CLB19 PPR editing factor in Arabidop-
sis [36, 39]. The design of these factors involved concatenat-
ing synthetic domains derived from consensus P, L, S, P2, L2,
S2, E1, E2, and DYW motifs to create protein scaffolds that
 mimic natural plant editing factors (Fig. 1A). Here, rather than
 replicating the natural architecture of plant editing factors, we
took a different approach to design PPR editors (dPPRe) for
 precise *de novo* RNA editing in plant organelles. We hypoth-
esized that any sufficiently long PPR tract with RNA-binding
 specificity for the sequence located six nucleotides upstream
 of a target cytidine (C) could effectively position the DYW
domain and enable cytidine deamination *in vivo*.

Unlike PLS-type PPR arrays, which contain L motifs whose
contribution to RNA binding is unpredictable, we reasoned
that P-type PPR arrays would be more suitable for reprogram-
ming RNA editing specificity *in vivo*. This is due to their clear
and predictable RNA-binding specificity, as all P motifs are in-
volved in sequence recognition [25, 26]. To test this, we con-
structed guides composed of a tandem array of 10 synthetic
P-type PPR motifs, termed designer PPR (dPPR) [37] (Fig.
1B), which we and others have previously demonstrated to
possess programmable RNA-binding specificity *in planta* [41,
42, 53]. A synthetic consensus P2–S2–L2 triplet was added
after the last P motif to replicate the structural arrangement found in natural PPR editing factors [39]. While the P2 motif is expected to confer RNA-binding specificity similar to P motifs, the roles of the L2 and S2 motifs in editing factors remain less clear [28, 54]. A guide composed of 11 P-type PPR repeats should be sufficient to specifically target an 11-nucleotide RNA sequence within plant chloroplast and mitochondrial genomes of few hundred kilobase pairs.

For the cytidine deamination reaction, we appended the C-terminal E1–E2–DYW domains from the *Physcomitrium patens* mitochondrial editing factor PPR56 to the synthetic [P]₁₀–P2–L2–S2 scaffold. This particular E-DYW domain was chosen due to PPR56’s demonstrated ability to catalyze cytidine deamination across various biological environments, including plant mitochondria [45] and cytoplasm [54], *E. coli* [28, 33], and human cell cytoplasm [55]. To stabilize the chimeric protein, an N-terminal domain from the maize chloroplast-localized PPR protein PPR10 [56] was added upstream of the first P motif, resulting in the engineered dPPRe. Finally, we applied the PPR code to the [P]₁₀-P2 motifs to target the RNA sequence located six nucleotides upstream of the selected cytidines, respectively.

For *in planta* editing assays, an N-terminal chloroplast transit peptide from Arabidopsis RecA (first 68 amino acids) or a mitochondrial presequence from yeast Cox4 (first 29 amino acids), both of which have been shown to efficiently direct fusion proteins to the respective organelles in plant cells [41, 57], were added to the dPPRe constructs. A 3×HA tag was also included for immunodetection.

Initially, we engineered two chloroplast dPPRe constructs, $dPPRe^{rpl2}$ and $dPPRe^{nddB}$, to target specific cytidines in the *rpl2* and *nddB* mRNAs of *N. benthamiana*, whose chloroplast genome has been recently sequenced [49] (Fig. 1B). These chloroplast genes encode a subunit of the large subunit of the chloroplast ribosome and the NADH dehydrogenase-like complex, respectively. Natively, these specific cytidines are not edited in the chloroplast transcriptome of *N. benthamiana* [49]. Editing the cytidine in the duplicated *rpl2* gene at positions 87,669 or 154,376 would affect the last base of the codon (CUC > CUU) without altering the amino acid sequence of the RPL2 protein (Leu > Leu) (Fig. 1B). In contrast, editing the cytidine at position 14,727 in *nddB* mRNA is expected to change a CAA codon to UAA, introducing a premature stop codon and truncating the NdhB protein.

Before expressing the constructs *in planta*, we verified their functionality using a well-established heterologous *E. coli* RNA editing assay that allows co-expression of rdPPRe without their transit peptide as N-terminal fusion proteins with a six-histidine and MBP tags, along with their RNA targets [33] (Fig. 2). Upon IPTG induction of $rdPPRe^{rpl2}$ and $rdPPRe^{nddB}$ expression in *E. coli*, a specific ~130 kDa protein band corresponding to the expected molecular weight of the recombinant proteins accumulated in the bacterial total protein extracts (Fig. 2A). Sanger sequencing of cDNA from these extracts confirmed that rdPPRe successfully edited their target cytidines in *E. coli* (Fig. 2B and Supplementary Table S1). Specifically, $rdPPRe^{rpl2}$ edited its target cytidine with an efficiency ranging from 49% to 62% (average of 57% ± 6.8%), while $rdPPRe^{nddB}$ exhibited a lower editing efficiency, averaging 33% ± 3.7%. These results demonstrate that our chimeric dPPRe are functional for RNA editing *in vivo* and that synthetic P-type PPR arrays can effectively guide RNA editing by the DYW domain, similar to PLS- and S-type scaffolds.

## Expression of dPPRe in plants induces *de novo* RNA editing in chloroplasts

Next, we evaluated the RNA editing activities of $dPPRe^{rpl2}$ and $dPPRe^{nddB}$ in chloroplasts through agroinfiltration-mediated transient gene expression in *N. benthamiana*. To enhance transient expression levels, dPPRe constructs were co-expressed with the tombusvirus suppressor of silencing, P19-HA [48]. Initially, we optimized assay conditions by testing *in vivo* protein accumulation of the most efficient dPPRe from *E. coli*, $dPPRe^{rpl2}$, when driven by a single or double 35S promoter (2×35S), a strong and constitutive promoter, at 2 and 3 days post-infiltration (dpi). Immunodetection with anti-HA antibodies revealed that the 2×35S promoter facilitated higher accumulation of $dPPRe^{rpl2}$ protein (expected molecular weight of ~97 kDa) in *N. benthamiana* leaves compared to the single 35S promoter, with protein levels peaking at 3 dpi for both promoters (Fig. 3A).

In parallel, we assessed RNA ERs at the *rpl2-87,669/154,376* sites in chloroplasts upon $dPPRe^{rpl2}$ expression by cDNA sequencing. The results demonstrated that $dPPRe^{rpl2}$ induced editing of the target cytidine in *rpl2* mRNA to an average of 38% at 2 dpi, increasing to 63% at 3 dpi while the expression of P19-HA alone did not cause editing (Fig. 3B and C, and Supplementary Table S2). A t-test confirmed the editing increase to be highly significant (*P*-value < 0.05). Notably, editing efficiency at the target site could reach up to 77% at 3 dpi (Fig. 3C and Supplementary Table S2). These findings suggest that $dPPRe^{rpl2}$ functions as an effective RNA editing factor in plant chloroplasts, with higher activity observed at 3 dpi, indicating a doseable effect of $dPPRe^{rpl2}$ on *rpl2* editing efficiency *in vivo*.

Subsequently, we expressed $dPPRe^{nddB}$ in *N. benthamiana* and examined RNA editing at *nddB*-14,727 site at 3 dpi. Like $dPPRe^{rpl2}$, $dPPRe^{nddB}$ showed higher protein accumulation with the 2×35S promoter (Fig. 3A). However, $dPPRe^{nddB}$ induced a lower ER *in planta*, achieving an average of 19% C-to-U conversion at its target site, reflecting the results obtained in *E. coli* (Fig. 3B and C, and Supplementary Table S2). This difference in editing efficiency between $dPPRe^{rpl2}$ and $dPPRe^{nddB}$ in *E. coli* and *in planta* was not due to differences in protein accumulation (Figs 2A and 3A). To determine whether this could be attributed to a difference in dPPRe affinity for the respective RNA molecules, electrophoretic mobility shift assays were performed using rdPPRe fused to the N-terminal MBP, and 19-nt RNAs spanning the -17/+1 cis-elements of the target cytidines, respectively, and with minimal expected secondary structures (Fig. 4A and B, and Supplementary Fig. S1). The results showed that $rdPPRe^{rpl2}$ and $rdPPRe^{nddB}$ bind their RNA targets with sequence specificity, whereas rMBP alone exhibited no RNA-binding activity (Fig. 4B). Quantification of the binding affinity revealed that both proteins bind their respective RNA targets with an approximate $K_D$ in the nanomolar range, with $rdPPRe^{nddB}$ displaying a higher affinity for its RNA sequence than $rdPPRe^{rpl2}$ ($K_D$ ~2.6 nM versus 14.6 nM) (Fig. 4C and Supplementary Fig. S2). These results indicate that RNA-binding affinity is not responsible for the lower *in vivo* editing activity of $dPPRe^{nddB}$ compared to $dPPRe^{rpl2}$.

Since PPR protein binding to RNA can be hindered by the formation of RNA secondary structures [15, 58, 59], we hypothesize that the stronger hairpin structure predicted at the $dPPRe^{nddB}$ binding site *in vivo* could interfere with RNA recognition (Supplementary Fig. S3). This structural
<figure><img src="image_11.png" /><figcaption>**Figure 2.** RNA editing activity of rdPPRe in *E. coli*. (A) Analysis of rdPPRe<sup>ndhB</sup>, rdPPRe<sup>rpl2</sup>, and rdPPRe<sup>nad7</sup> protein expression in *E. coli* by SDS–PAGE and Coomassie blue staining of total protein extracts from three independent bacterial cultures harvested before induction (bi) or 20 h after expression of constructs was induced. The arrowhead indicates protein bands whose apparent molecular weights correspond to the expected sizes for rdPPRe (131.5 kDa). (B) Sequencing electrophoregrams of cDNA from three independent *E. coli* cultures expressing rdPPRe<sup>ndhB</sup>, rdPPRe<sup>rpl2</sup>, or rdPPRe<sup>nad7</sup> for 20 h. Expression of rdPPRe in bacteria resulted in C-to-U editing at the co-expressed rpl2, ndhB, and nad7 target sites, as revealed by the appearance of a T peak in the cDNA sequence at these positions. The average ERs are shown below (n = 3).</figcaption></figure>

constraint may reduce the accessibility of the target sequence, thereby impairing $dPPRe^{ndhB}$ binding and leading to the lower RNA editing efficiency observed *in vivo*.

Overall, these findings demonstrate that chimeric dPPRe proteins, composed of synthetic P-PPR guides and a mitochondrial E1–E2–DYW domain from moss, are functional in chloroplasts of angiosperms for *de novo* RNA editing. The results also suggest that the local folding of RNA target sequences plays an important role in determining dPPRe editing efficiency *in vivo*, consistent with previous reports for natural PPR editing factors [28, 60].

### $dPPRe^{rpl2}$ exhibits minimal off-target activity in planta

To evaluate the off-target RNA editing activity of dPPRe, we focused on the most efficient variant, $dPPRe^{rpl2}$, and analyzed potential off-target sites in *N. benthamiana* chloroplasts. We performed RNA sequencing (RNA-seq) on samples from three independent agroinfiltration experiments. In these experiments, leaves were infiltrated either with the P19 protein alone (control) or with $dPPRe^{rpl2}$/P19 driven by either a 35S or 2×35S promoter, resulting in nine cDNA libraries. To assess the relationship between $dPPRe^{rpl2}$ expression levels and RNA ERs, RNA was extracted at two time points: 2 dpi for one cDNA library and 3 dpi for two cDNA libraries in each of the three experiments.

Differential expression analysis using DESeq2 revealed no significant changes in chloroplast gene expression across the samples, with the exception of the $dPPRe^{rpl2}$ nuclear transgene, which was significantly overexpressed in the experimental samples (Fig. 5A and Supplementary Table S3). This find-

ing indicates that $dPPRe^{rpl2}$ expression does not disrupt the overall regulation of chloroplast gene expression.

Next, we used the GATK4 toolkit to identify *de novo* C-to-U RNA editing sites triggered by $dPPRe^{rpl2}$ expression (Supplementary Table S4). As expected, the target cytidine in the *rpl2* gene was edited (Fig. 5B). However, we also detected three off-target editing sites in the *matK*, *rpoC1*, and *psbB* chloroplast genes. Editing at the *rpl2* target site and the *matK-2329* site was consistently observed across all $dPPRe^{rpl2}$-expressing libraries, while editing at the *rpoC1-22,069* and *psbB-75,097* sites was detected in only five of the six libraries (Fig. 5B).

We quantified the ER at these four sites in each experimental library in relation to the expression levels of $dPPRe^{rpl2}$ (Fig. 5B). Consistent with previous protein immunodetection results, the 2×35S promoter drove higher expression of $dPPRe^{rpl2}$ than the single 35S promoter, and $dPPRe^{rpl2}$ expression was higher at 3 dpi compared to 2 dpi for both promoters. The quantification revealed that the *rpl2* target site had the highest ER across all libraries, reaching up to 67.8%, followed by *matK-2329* (ER ≤ 46.3%), *rpoC1-22,069* (ER ≤ 31.4%), and *psbB-75,097* (ER ≤ 21.9%) (Fig. 5B). These results suggest that $dPPRe^{rpl2}$ preferentially guides RNA editing towards its intended cytidine target over off-target sites *in vivo*. Moreover, comparison of ERs across independent experiments indicated a positive correlation between $dPPRe^{rpl2}$ expression levels and RNA editing efficiency (Fig. 5B), particularly at the *rpl2* target site, as confirmed by a Spearman correlation test (*P*-value < 0.05) (Fig. 5C and Supplementary Fig. S4).

To further investigate the off-target activity of $dPPRe^{rpl2}$, we aligned the cis-element sequences spanning from nucleotide -16 to +1 relative to each off-target site with that of the *rpl2*
<figure><img src="image_12.png" /><figcaption>**Figure 3.** RNA editing activity of dPPRe *in planta*. (A) Immunodetection of chloroplast dPPRe<sup>rp12</sup>, dPPRe<sup>ndhB</sup>, and mitochondrial dPPRe<sup>nad7</sup> accumulation in *N. benthamiana* leaves at 2 and 3 days post-infiltration (dpi). Arrowheads indicate the protein bands corresponding to dPPRe (~97 kDa) and the suppressor of silencing, P19 protein (~19 kDa). (B) Sequencing electrophoregrams of cDNA from a replicate of each *in planta* experiment. Expression of dPPRe in *N. benthamiana* leaf led to *de novo* C-to-U RNA editing at the target sites in endogenous chloroplast *rpl2*, *ndhB*, and mitochondrial *nad7* transcripts as revealed by the appearance of a T peak in the cDNA sequence at these positions that is absent in the control cDNA (P19 alone). The average ERs with standard deviation (SD) and the number (*n*) of independently treated leaves are shown below (Supplementary Table S2). (C) ERs at the respective targeted dPPRe sites for individual replicates from each experimental set up. For P19/dPPRe<sup>rp12</sup>, the 2 dpi time point was assessed in three (*n* = 3) independently treated *N. benthamiana* leaves, and the 3 dpi time point was assessed in six (*n* = 6) leaves. The 3 dpi time point for P19/dPPRe<sup>ndhB</sup> and P19/dPPRe<sup>nad7</sup> was assessed in seven (*n* = 7) independently treated leaves. The P19 control for P19/dPPRe<sup>nad7</sup> at 3 dpi was assessed in six (*n* = 6) leaves. Data are presented as means ± SD and were compared using *t*-tests with *P*-values shown in the graph (Supplementary Table S2).</figcaption></figure>

target site. This analysis revealed a conserved nucleotide sequence among the off-target and *rpl2* cis-elements (Fig. 5D). Notably, the first 11 nucleotides, predicted to be bound by the synthetic [P]₁₀-P2 guide, exhibited high conservation, with all off-target sites sharing 8 nucleotides with the *rpl2* target site within this region. Additionally, a U was observed at position –3, and an A immediately downstream of the edited cytidine in all off-targets.

Further analysis of the nucleotide preferences of the 11 PPRs in dPPRe<sup>rp12</sup> showed that some motifs tolerate mismatches along the 11-nucleotide target sequence (Fig. 5D). The number and position of these mismatches within the PPR-RNA duplex appeared to influence *in vivo* editing efficiency, likely by affecting protein-RNA interactions. For example, the most efficiently edited off-target, *matK-2329*, contained two mismatches positioned one nucleotide apart in the central region of the dPPR binding site. In contrast, the less efficiently

edited sites, *rpoC1-22,069* and *psbB-75,097*, exhibited two to three mismatches distributed differently: *rpoC1-22,069* had three mismatches spread across the 5' and 3' ends of the RNA, while *psbB-75,097* had two mismatches located between the 5' end and the middle of the RNA. These findings suggest that the dPPR guide is more tolerant of mismatches in the central region of its binding site than at the 5' or 3' ends.

## dPPR editors are active in plant mitochondria

Finally, we evaluated the potential of dPPRe for *de novo* RNA editing in plant mitochondria. For this purpose, we engineered dPPRe<sup>nad7</sup> to edit a specific cytidine in the *nad7* mRNA within *N. benthamiana* mitochondria, which encodes a subunit of the respiratory complex I (Fig. 1B). Since the mitochondrial genome of *N. benthamiana* has not yet been sequenced, we selected this editing site based on the sequence of the mito-
<figure><img src="image_13.png" /><figcaption>**Figure 4.** RNA binding activity of rdPPRe<sup>rpl2</sup> and rdPPRe<sup>ndhb</sup>. (A) SDS-PAGE analysis of recombinant rdPPRe and rMBP proteins purified from *E. coli*. A total of 100 ng of each protein was loaded onto the gel. (B) Electrophoretic mobility shift assays showing the RNA binding specificity of rdPPRe. Increasing protein concentrations (0, 19, 38, 76 nM) are indicated by triangles. The RNA probe used is depicted at the bottom left of the gel and the target cytidine is underlined. B: bound RNA, U: unbound RNA. The sequences of the RNA probes are shown below. (C) Equilibrium dissociation constants (K<sub>D</sub>) of rdPPRe<sup>rpl2</sup> and rdPPRe<sup>ndhb</sup> for their respective rpl2 and ndhb RNA targets. Binding reactions were performed as in panel (A) (n = 2).</figcaption></figure>

chondrial *nad7* gene from *Nicotiana tabacum*, assuming that the sequence is conserved between the two species. Similar to the editing event at *ndhb*-14,727, RNA editing at the selected cytidine in the mitochondrial *nad7* gene is expected to create a premature stop codon in the mRNA (CAA > UAA) (Fig. 1B).

The editing activity of rdPPRe<sup>nad7</sup> was assessed in *E. coli* (Fig. 2). Upon induction, rdPPRe<sup>nad7</sup> accumulated to a much higher level in *E. coli* compared to rdPPRe<sup>rpl2</sup> and rdPPRe<sup>ndhb</sup>, leading to RNA editing at the *nad7* target site with an average efficiency of 74% (Fig. 2B). While the secondary structure of *nad7* RNA expressed in *E. coli* is not predicted to be stronger than that of rpl2 or ndhb RNAs (Supplementary Fig. S3), the higher accumulation of rdPPRe<sup>nad7</sup> compared to rdPPRe<sup>rpl2</sup> and rdPPRe<sup>ndhb</sup> may be responsible for its increased editing efficiency in bacteria (Fig. 2A). Finally, dPPRe<sup>nad7</sup> was expressed under the control of a 2×35S promoter in *N. benthamiana* leaves and successfully accumulated by 3 dpi, as confirmed by anti-HA immunodetection. cDNA sequencing confirmed *nad7* sequence and revealed that the expression of dPPRe<sup>nad7</sup> in plant cells induced RNA editing at the target cytidine within mitochondria, achieving an ER of up to 42% (Fig. 3B and C). This result demonstrates that synthetic dPPRe are functional ribocytidine editors in plant mitochondria.

## Discussion

The results of this study highlight the significant potential of synthetic designer PPR RNA editors, as a powerful tool for precise RNA base editing in plant organelles, with successful applications demonstrated in both chloroplasts and mitochondria. Our engineering of dPPRe constructs, which incorporate arrays of synthetic P-type PPR motifs [37] fused with the DYW cytidine deaminase domain from the moss PPR56 editing factor, represents a significant advancement over previous designs that relied on PLS-type or SS-type PPR arrangements.

Previously, the first synthetic PLS-PPR editing factor (referred to as dsnPLS<sup>CLB19</sup>) was capable to substitute the func-

tion of the natural CLB19 PLS-PPR editing factor at the rpoA-78,691 site in the respective Arabidopsis knockout line [39]. However, efficient editing was only possible in the presence of plant protein cofactors from the RIP/MORF family (Multiple Organellar RNA editing Factors). Structural and biochemical studies suggest that MORF interaction with PLS-PPR editing factors, particularly with the L motifs, enhances RNA-binding efficiency by inducing conformational changes [61]. Consequently, the application of synthetic PLS-PPR scaffolds may be limited to angiosperm plants that possess MORF proteins [62, 63].

To overcome the MORF-dependency, a subsequent SS-type PPR scaffold (dsnS<sup>CLB19</sup>) was developed using a similar approach. dsnS<sup>CLB19</sup> edited its RNA target in *E. coli* without plant MORF cofactors, achieving an ER below 50%, similar to dsnPLS<sup>CLB19</sup> [36]. However, the effectiveness of the dsnS scaffold *in planta* and the broader applicability of dsnPLS- and dsnS-type editors for reprogramming RNA editing specificity in plant organelles have not been tested yet. In contrast, our study demonstrates that synthetic P-type dPPR editors can induce *de novo*, site-specific RNA editing in both chloroplasts and mitochondria achieving high efficiency above 60%. This significantly expands the RNA editing toolkit for plant systems.

Why plants have evolved different types of PPR arrays to fulfill distinct functions remains an open question [64]. A direct comparison of the performance of designer P-, S-, and PLS-type PPR arrays fused to DYW deaminase domains, in terms of efficiency and specificity, would be a valuable experiment to conduct in the future.

However, the preference for PLS-type PPR arrays over classical P-type arrays in land plants for RNA editing site recognition may be linked to the function and positioning of L-type motifs. Unlike P- and S-type PPRs, L-type motifs do not follow a clear recognition code and do not seem to be directly involved in the recognition of individual nucleotides. This characteristic could have facilitated the conserved recognition of specific editing sites in messenger RNAs throughout
<figure><img src="image_14.png" /><figcaption>**Figure 5.** Off-target RNA editing activity of dPPRe<sup>rp12</sup> in *N. benthamiana* chloroplasts. (A) MA-plot displaying differentially expressed chloroplast genes and the dPPRe<sup>rp12</sup> transgene highlighted as a red dot. The plot displays log<sub>2</sub> fold-changes in expression for dPPRe<sup>rp12</sup>-expressing samples versus the normalized mean expression (log<sub>2</sub> scale) for each gene. The dPPRe<sup>rp12</sup> transgene is the only gene with statistically significant differential expression (adjusted P-value < 0.05). (B) Bar plots depicting DESeq2-normalized dPPRe<sup>rp12</sup> expression in each experimental library, along with the corresponding C-to-U editing efficiency at identified *de novo* RNA editing sites. dpi: days post-inoculation. (C) Spearman correlation analysis between dPPRe<sup>rp12</sup> expression level and editing efficiency at *rp12* target site. (D) Alignment between dPPRe<sup>rp12</sup> with its predicted binding sites upstream of the *de novo* chloroplast editing sites. A sequence logo was generated from the alignment of the cis-element sequences (–16/+1) at these sites. The PPR code determining amino acids five and last of each PPR motif are indicated, with nucleotide binding compatibility in the binding site sequences shaded from black to light gray, where light gray indicates no compatibility between the PPR code and the nucleotide in the binding site.</figcaption></figure>

plant evolution. Given that RNA editing sites frequently occur at the second position of codons and that PLS-type PPRs bind approximately six nucleotides upstream of the editing site, L-motifs are typically positioned opposite the third codon position—an area generally less conserved across species due to the degeneracy of the genetic code. By omitting these variable positions, PLS-type PPRs may have ensured greater target recognition stability over evolutionary time [64]. However, this evolutionary advantage does not necessarily translate into a benefit for the design of synthetic editing factors.

The DYW deaminase of PPR56 has been shown to function in various biological systems including human cells, plants, and bacteria without perturbating the target cell [33, 45, 55]. In this study, we demonstrated that when fused to a programmable P-type PPR array, this deaminase can rewrite the genetic information within a target RNA by enzymatically converting a specific cytidine in both plant chloroplasts and mitochondria. This makes dPPRe a versatile tool for site-directed RNA base editing. However, like other site-

directed RNA editing approaches, dPPRe face challenges related to editing efficiency and target scope. Although dPPRe can achieve a satisfactory RNA editing efficiency over 50% for biotechnological applications, certain plant PPR editing factors, such as the native PPR56, can nearly fully convert their cytidine target *in vivo* [33, 45]. The lower efficiency of the DYW domain of PPR56 in dPPRe may be due to the cooperative role of this domain with the PPR array for RNA recognition and editing [65, 66]. This cooperation might not be fully maintained in synthetic PPR scaffolds, potentially resulting in reduced *in vivo* activity.

In addition, similar to ADAR (Adenosine Deaminases Acting on RNA) and APOBEC (Apolipoprotein B Editing Catalytic polypeptide) base editors [67–69], the DYW deaminases of natural PPR editing factors and their associated E1–E2 domains have nucleotide context preferences, preferring cytidine that are embedded in specific RNA sequences [65, 66, 70]. For instance, the –3/+1 nucleotide sequences surrounding the dPPRe<sup>rp12</sup> off-target cytidines in chloroplasts are rich in U/C
and A (Fig. 5D), aligning with the described preferences of PPR56 [54, 65, 70]. This nucleotide context preference might also explain the lower editing efficiency of dPPRe<sup>ndhB</sup> *in vivo*: sequence analysis of dPPRe<sup>rpl2</sup> off-targets revealed a strict preference of PPR56's E1–E2–DYW domain for U at position –3, while the *ndhB* target C site has an A at this position (Fig. 1B). While this nucleotide preference is beneficial for the engineering of site-specific RNA editors by increasing editing precision, it limits the number of editable cytidines in a given transcriptome.

Furthermore, this study provides important insights into the precision and efficiency of RNA editing achieved by dPPRe *in planta*. While dPPRe<sup>rpl2</sup> exhibited robust ERs at its target site in the *rpl2* gene *in planta*, we also detected a few off-targets. These off-target sites, however, shared significant sequence similarities with the intended target, particularly in the regions predicted to interact with the synthetic PPR motifs. Interestingly, our results showed that when synthetic P-type PPR motifs are fused with the C-terminal P2–L2–S2–E1–E2–DYW domains to create dPPRe, these motifs primarily tolerate mismatches in the central region of their RNA binding site (Fig. 5D). This indicates that, in the context of dPPRe, the P-type PPR motifs may exhibit somewhat different binding behavior compared to purely synthetic P-type dPPR scaffold (composed entirely of P-type PPR motifs), which typically shows more permissiveness to mismatches at the 3' end of the RNA binding site [38]. In contrast, natural PLS-type PPR editing factors have been shown to tolerate mismatches more frequently at the 5' end of the binding site, with less tolerance towards the 3' end, close to the editing site [33, 71]. Given that dPPRe combines features of both P-type and PLS-type architectures, it is not surprising that mismatch tolerance is primarily constrained to the central region, while stricter complementarity is required at the 5' and 3' ends. This suggests that the hybrid nature of dPPRe imposes a unique pattern of mismatch tolerance, influenced by both the P-type and PLS-type elements, leading to a distinct RNA binding specificity. Although the synthetic S2 motif used in this study was previously reported to preferentially recognize *A in vivo* and contribute to the RNA site specificity in synthetic PLS editing factors [39], our findings suggest otherwise. In fact, all four RNA bases are represented at this position in the *in vivo* target sites of dPPRe<sup>rpl2</sup> (Fig. 5D).

These findings underscore the importance of carefully designing PPR motifs to minimize unintended interactions, as PPR-RNA interactions are not solely defined by the PPR RNA-binding code itself. For instance, when we applied the PPR code to predict binding sites for the [P]₁₀-P2 motifs of dPPRe<sup>rpl2</sup> and possible cytidine targets within the chloroplast genome of *N. benthamiana* using the Prepect Tool TargetScan [72], the designated cytidine target in *rpl2* was identified as the top hit (position 154,376; Supplementary Fig. S5). However, none of the other predicted targets were edited in plants expressing dPPRe<sup>rpl2</sup>, and the three off-targets observed *in planta* were not among the top 40 predicted matches. This discrepancy lies in the difficulty to evaluate the energetic costs of mismatches and their positional effects within a PPR-RNA duplex, competing RNA secondary structures or protein occupancy, and the influence of the C-terminal DYW domains on nucleotide specificity [54, 65, 70, 73]. While one possible strategy to enhance specificity could be to increase the length of synthetic PPR tracts, previous studies indicate that elongating the dPPR array beyond 11 motifs does not systematically

improve target discrimination. In fact, increasing the number of PPR motifs can, in some cases, lead to greater tolerance of mismatches rather than enhanced specificity [38, 53]. A deeper understanding of these factors will significantly enhance our ability to predict PPR-RNA interactions and improve the specificity of dPPRe design without compromising editing efficiency in the future.

Recently, the use of the *E. coli* editing system and a plant cytosolic editing assay, demonstrated that the DYW domains from plant editing factors besides PPR56 show some variations in their flexibility towards the nucleotide context around the cytidine target site [54, 70]. These DYW domains offer great promise for broadening the editing site-specific capability of dPPRe. Although synthetic E1, E2, and DYW domains engineered from consensus amino acid sequences are currently available [39], their preferred nucleotide context remains unknown and complicates their use for *de novo* RNA editing *in vivo*.

In conclusion, the successful application of dPPRe in *E. coli*, chloroplasts and mitochondria emphasizes their versatility for targeted RNA editing and provides the first mitochondrial RNA base editor tool kit. Mitochondrial dysfunction is associated with a range of phenotypic abnormalities in plants, including cytoplasmic male sterility, which is of considerable interest in hybrid seed production [74]. The ability to precisely edit mitochondrial genes could therefore have profound implications for the development of new plant varieties with improved traits. Beyond plant systems, dPPRe hold great promise for mitochondrial research and associated human diseases [44]. Given that PPR56 is active for RNA editing in human cells [55], dPPRe technology could pave the way for innovative therapeutic strategies targeting mitochondrial genes. As we continue to refine and expand the capabilities of dPPRe, this technology holds the potential to become a versatile and indispensable tool in organellar genetic research.

## Acknowledgements

The authors thank Volker Knoop for useful discussions and Sarah Schwirblat for technical assistance.

*Author contributions:* S.M., E.L., M.S.-R., and K.H. designed and conceived the experiments. E.L. performed the *E. coli* experiments. S.M., Sh.G., and K.H. performed the *in planta* work. K.H. performed EMSA assays. S.G. performed bioinformatic analysis. S.M., E.L., S.G., and K.H. analyzed the results. K.H. supervised the work, made the figures, and drafted the manuscript. S.M., E.L., and M.S.-R. edited the manuscript.

## Supplementary data

Supplementary data is available at NAR online.

## Conflict of interest

None declared.

## Funding

This work was supported by Agence National de la Recherche [ANR-18-CE20-0013 to K.H.] and Centre National de la Recherche Scientifique. Work of M.S-R. and E.L. was funded by the Deutsche Forschungsgemeinschaft [SCHA 1952/2-2].
Funding to pay the Open Access publication charges for this article was provided by Centre National de la Recherche Scientifique.

## Data availability

The read data were deposited to the NCBI Sequence Read Archive database with the Bioproject identifier GSE276199.

## References

1. Gray MW. Evolutionary origin of RNA editing. *Biochemistry* 2012;51:5235-42. https://doi.org/10.1021/bi300419r

2. Simpson L, Emeson RB. RNA editing. *Annu Rev Neurosci* 1996;19:27-52.
https://doi.org/10.1146/annurev.ne.19.030196.000331

3. Knoop V. When you can't trust the DNA: RNA editing changes transcript sequences. *Cell Mol Life Sci* 2011;68:567-86.
https://doi.org/10.1007/s00018-010-0538-9

4. Hoch B, Maier RM, Appel K et al. Editing of a chloroplast mRNA by creation of an initiation codon. *Nature* 1991;353:178-80.
https://doi.org/10.1038/353178a0

5. Hiesel R, Wissinger B, Schuster W et al. RNA editing in plant mitochondria. *Science* 1989;246:1632-4.
https://doi.org/10.1126/science.2480644

6. Gualberto JM, Lamattina L, Bonnard G et al. RNA editing in wheat mitochondria results in the conservation of protein sequences. *Nature* 1989;341:660-2.
https://doi.org/10.1038/341660a0

7. Knoop V. C-to-U and U-to-C: RNA editing in plant organelles and beyond. *J Exp Bot* 2023;74:2273-94.
https://doi.org/10.1093/jxb/erac488

8. Small ID, Schallenberg-Rudinger M, Takenaka M et al. Plant organellar RNA editing: what 30 years of research has revealed. *Plant J* 2020;101:1040-56. https://doi.org/10.1111/tpj.14578

9. Shikanai T. RNA editing in plant organelles: machinery, physiological function and evolution. *Cell Mol Life Sci* 2006;63:698-708. https://doi.org/10.1007/s00018-005-5449-9

10. Ruwe H, Castandet B, Schmitz-Linneweber C et al. Arabidopsis chloroplast quantitative editotype. *FEBS Lett* 2013;587:1429-33.
https://doi.org/10.1016/j.febslet.2013.03.022

11. Bentolila S, Elliott LE, Hanson MR. Genetic architecture of mitochondrial editing in *Arabidopsistaliana*. *Genetics* 2008;178:1693-708. https://doi.org/10.1534/genetics.107.073585

12. Schallenberg-Rudinger M, Kindgren P, Zehrmann A et al. A DYW-protein knockout in Physcomitrella affects two closely spaced mitochondrial editing sites and causes a severe developmental phenotype. *Plant J* 2013;76:420-32.
https://doi.org/10.1111/tpj.12304

13. Ichinose M, Sugita C, Yagi Y et al. Two DYW subclass PPR proteins are involved in RNA editing of *ccmFc* and *atp9* transcripts in the moss Physcomitrella patens: first complete set of PPR editing factors in plant mitochondria. *Plant Cell Physiol* 2013;54:1907-16. https://doi.org/10.1093/pcp/pct132

14. Chateigner-Boutin AL, Ramos-Vega M, Guevara-Garcia A et al. CLB19, a pentatricopeptide repeat protein required for editing of *rpoA* and *clpP* chloroplast transcripts. *Plant J* 2008;56:590-602.
https://doi.org/10.1111/j.1365-313X.2008.03634.x

15. Hammani K, des Francs-Small CC, Takenaka M et al. The pentatricopeptide repeat protein OTP87 is essential for RNA editing of *nad7* and *atp1* transcripts in Arabidopsis mitochondria. *J Biol Chem* 2011;286:21361-71.
https://doi.org/10.1074/jbc.M111.230516

16. Schmitz-Linneweber C, Kushnir S, Babiychuk E et al. Pigment deficiency in nightshade/tobacco cybrids is caused by the failure to edit the plastid ATPase alpha-subunit mRNA. *Plant Cell* 2005;17:1815-28. https://doi.org/10.1105/tpc.105.032474

17. Kotera E, Tasaka M, Shikanai T. A pentatricopeptide repeat protein is essential for RNA editing in chloroplasts. *Nature* 2005;433:326-30. https://doi.org/10.1038/nature03229

18. Peeters N, Small I. Dual targeting to mitochondria and chloroplasts. *Biochim Biophys Acta* 2001;1541:54-63.
https://doi.org/10.1016/S0167-4889(01)00146-X

19. Lurin C, Andres C, Aubourg S et al. Genome-wide analysis of Arabidopsis pentatricopeptide repeat proteins reveals their essential role in organelle biogenesis. *Plant Cell* 2004;16:2089-103. https://doi.org/10.1105/tpc.104.022236

20. Shen C, Zhang D, Guan Z et al. Structural basis for specific single-stranded RNA recognition by designer pentatricopeptide repeat proteins. *Nat Commun* 2016;7:11285.
https://doi.org/10.1038/ncomms11285

21. Barkan A, Small I. Pentatricopeptide repeat proteins in plants. *Annu Rev Plant Biol* 2014;65:415-42.
https://doi.org/10.1146/annurev-arplant-050213-040159

22. Cheng S, Gutmann B, Zhong X et al. Redefining the structural motifs that determine RNA binding and RNA editing by pentatricopeptide repeat proteins in land plants. *Plant J* 2016;85:532-47. https://doi.org/10.1111/tpj.13121

23. Gutmann B, Royan S, Schallenberg-Rudinger M et al. The expansion and diversification of pentatricopeptide repeat RNA-editing factors in plants. *Mol Plant* 2020;13:215-30.
https://doi.org/10.1016/j.molp.2019.11.002

24. O'Toole N, Hattori M, Andres C et al. On the expansion of the pentatricopeptide repeat gene family in plants. *Mol Biol Evol* 2008;25:1120-8. https://doi.org/10.1093/molbev/msn057

25. Barkan A, Rojas M, Fujii S et al. A combinatorial amino acid code for RNA recognition by pentatricopeptide repeat proteins. *PLoS Genet* 2012;8:e1002910.
https://doi.org/10.1371/journal.pgen.1002910

26. Kobayashi K, Kawabata M, Hisano K et al. Identification and characterization of the RNA binding surface of the pentatricopeptide repeat protein. *Nucleic Acids Res* 2012;40:2712-23. https://doi.org/10.1093/nar/gkr1084

27. Yan J, Yao Y, Hong S et al. Delineation of pentatricopeptide repeat codes for target RNA prediction. *Nucleic Acids Res* 2019;47:3728-38. https://doi.org/10.1093/nar/gkz075

28. Yang Y, Ritzenhofen K, Otrzonsek J et al. Beyond a PPR-RNA recognition code: many aspects matter for the multi-targeting properties of RNA editing factor PPR56. *PLoS Genet* 2023;19:e1010733. https://doi.org/10.1371/journal.pgen.1010733

29. Baudry K, Monachello D, Castandet B et al. Dissecting the molecular puzzle of the editosome core in Arabidopsis organelles. *Plant Sci* 2024;344:112101.
https://doi.org/10.1016/j.plantsci.2024.112101

30. Andres-Colas N, Zhu Q, Takenaka M et al. Multiple PPR protein interactions are involved in the RNA editing system in Arabidopsis mitochondria and plastids. *Proc Natl Acad Sci USA* 2017;114:8883-8. https://doi.org/10.1073/pnas.1705815114

31. Guillaumot D, Lopez-Obando M, Baudry K et al. Two interacting PPR proteins are major Arabidopsis editing factors in plastid and mitochondria. *Proc Natl Acad Sci USA* 2017;114:8877-82.
https://doi.org/10.1073/pnas.1705780114

32. Toma-Fukai S, Sawada Y, Maeda A et al. Structural insight into the activation of an Arabidopsis organellar C-to-U RNA editing enzyme by active site complementation. *Plant Cell* 2023;35:1888-900. https://doi.org/10.1093/plcell/koac318

33. Oldenkott B, Yang Y, Lesch E et al. Plant-type pentatricopeptide repeat proteins with a DYW domain drive C-to-U RNA editing in *Escherichia coli*. *Commun Biol* 2019;2:85.
https://doi.org/10.1038/s42003-019-0328-3

34. Takenaka M, Takenaka S, Barthel T et al. DYW domain structures imply an unusual regulation principle in plant organellar RNA editing catalysis. *Nat Catal* 2021;4:510-22.
https://doi.org/10.1038/s41929-021-00633-x

35. Hayes ML, Santibanez PI. A plant pentatricopeptide repeat protein with a DYW-deaminase domain is sufficient for catalyzing C-to-U
RNA editing *in vitro*. *J Biol Chem* 2020;295:3497-505.
https://doi.org/10.1074/jbc.RA119.011790

36. Bernath-Levin K, Schmidberger J, Honkanen S *et al.* Cofactor-independent RNA editing by a synthetic S-type PPR protein. *Synth Biol (Oxf)* 2022;7:ysab034.
https://doi.org/10.1093/synbio/ysab034

37. Shen C, Wang X, Liu Y *et al.* Specific RNA recognition by designer pentatricopeptide repeat protein. *Mol Plant* 2015;8:667-70.
https://doi.org/10.1016/j.molp.2015.01.001

38. Miranda RG, McDermott JJ, Barkan A. RNA-binding specificity landscapes of designer pentatricopeptide repeat proteins elucidate principles of PPR-RNA interactions. *Nucleic Acids Res* 2018;46:2613-23. https://doi.org/10.1093/nar/gkx1288

39. Royan S, Gutmann B, Colas des Francs-Small C *et al.* A synthetic RNA editing factor edits its target site in chloroplasts and bacteria. *Commun Biol* 2021;4:545.
https://doi.org/10.1038/s42003-021-02062-9

40. Ichinose M, Kawabata M, Akaiwa Y *et al.* U-to-C RNA editing by synthetic PPR-DYW proteins in bacteria and human culture cells. *Commun Biol* 2022;5:968.
https://doi.org/10.1038/s42003-022-03927-3

41. Manavski N, Mathieu S, Rojas M *et al.* *In vivo* stabilization of endogenous chloroplast RNAs by customized artificial pentatricopeptide repeat proteins. *Nucleic Acids Res* 2021;49:5985-97. https://doi.org/10.1093/nar/gkab390

42. Rojas M, Chotewutmontri P, Barkan A. Translational activation by a synthetic PPR protein elucidates control of psbA translation in Arabidopsis chloroplasts. *Plant Cell* 2024;36:4168-78.
https://doi.org/10.1093/plcell/koae112

43. Gautam R, Shukla P, Kirti PB. Male sterility in plants: an overview of advancements from natural CMS to genetically manipulated systems for hybrid seed production. *Theor Appl Genet* 2023;136:195. https://doi.org/10.1007/s00122-023-04444-5

44. Song M, Ye L, Yan Y *et al.* Mitochondrial diseases and mtDNA editing. *Genes Dis* 2024;11:101057.
https://doi.org/10.1016/j.gendis.2023.06.026

45. Ohtani S, Ichinose M, Tasaki E *et al.* Targeted gene disruption identifies three PPR-DYW proteins involved in RNA editing for five editing sites of the moss mitochondrial transcripts. *Plant Cell Physiol* 2010;51:1942-9. https://doi.org/10.1093/pcp/pcq142

46. Curtis MD, Grossniklaus U. A gateway cloning vector set for high-throughput functional analysis of genes in planta. *Plant Physiol* 2003;133:462-9. https://doi.org/10.1104/pp.103.027979

47. Nakamura S, Nakao A, Kawamukai M *et al.* Development of Gateway binary vectors, R4L1pGWBs, for promoter analysis in higher plants. *Biosci Biotechnol Biochem* 2009;73:2556-9.
https://doi.org/10.1271/bbb.90720

48. Incarbone M, Ritzenthaler C, Dunoyer P. Peroxisomal targeting as a sensitive tool to detect protein-small RNA interactions through *in vivo* piggybacking. *Front Plant Sci* 2018;9:135.
https://doi.org/10.3389/fpls.2018.00135

49. Wang S, Gao J, Chao H *et al.* Comparative chloroplast genomes of nicotiana species (Solanaceae): insights into the genetic variation, phylogenetic relationship, and polyploid speciation. *Front Plant Sci* 2022;13:899252. https://doi.org/10.3389/fpls.2022.899252

50. Auwera G, O'Connor BD. *Genomics in the Cloud : Using Docker, GATK, and WDL in Terra* (First edition). Sebastopol, CA: O'Reilly Media, 2020.

51. DePristo MA, Banks E, Poplin R *et al.* A framework for variation discovery and genotyping using next-generation DNA sequencing data. *Nat Genet* 2011;43:491-8. https://doi.org/10.1038/ng.806

52. Hammani K, Cook WB, Barkan A. RNA binding and RNA remodeling activities of the half-a-tetratricopeptide (HAT) protein HCF107 underlie its effects on gene expression. *Proc Natl Acad Sci USA* 2012;109:5651-6. https://doi.org/10.1073/pnas.1200318109

53. McDermott JJ, Watkins KP, Williams-Carrier R *et al.* Ribonucleoprotein capture by *in vivo* expression of a designer pentatricopeptide repeat protein in Arabidopsis. *Plant Cell* 2019;31:1723-33. https://doi.org/10.1105/tpc.19.00177

54. Thielen M, Gartner B, Knoop V *et al.* Conquering new grounds: plant organellar C-to-U RNA editing factors can be functional in the plant cytosol. *Plant J* 2024;119:895-915.
https://doi.org/10.1111/tpj.16804

55. Lesch E, Schilling MT, Brenner S *et al.* Plant mitochondrial RNA editing factors can perform targeted C-to-U editing of nuclear transcripts in human cells. *Nucleic Acids Res* 2022;50:9966-83.
https://doi.org/10.1093/nar/gkac752

56. Pfalz J, Bayraktar OA, Prikryl J *et al.* Site-specific binding of a PPR protein defines and stabilizes 5' and 3' mRNA termini in chloroplasts. *EMBO J* 2009;28:2042-52.
https://doi.org/10.1038/emboj.2009.121

57. Peeters NM, Chapron A, Giritch A *et al.* Duplication and quadruplication of *Arabidopsis thaliana* cysteinyl- and asparaginyl-tRNA synthetase genes of organellar origin. *J Mol Evol* 2000;50:413-23. https://doi.org/10.1007/s002390010044

58. McDermott JJ, Civic B, Barkan A. Effects of RNA structure and salt concentration on the affinity and kinetics of interactions between pentatricopeptide repeat proteins and their RNA ligands. *PLoS One* 2018;13:e0209713.
https://doi.org/10.1371/journal.pone.0209713

59. Williams-Carrier R, Kroeger T, Barkan A. Sequence-specific binding of a chloroplast pentatricopeptide repeat protein to its native group II intron ligand. *RNA* 2008;14:1930-41.
https://doi.org/10.1261/rna.1077708

60. Kindgren P, Yap A, Bond CS *et al.* Predictable alteration of sequence recognition by RNA editing factors from Arabidopsis. *Plant Cell* 2015;27:403-16.
https://doi.org/10.1105/tpc.114.134189

61. Yan J, Zhang Q, Guan Z *et al.* MORF9 increases the RNA-binding activity of PLS-type pentatricopeptide repeat protein in plastid RNA editing. *Nat Plants* 2017;3:17037.
https://doi.org/10.1038/nplants.2017.37

62. Takenaka M, Zehrmann A, Verbitskiy D *et al.* Multiple organellar RNA editing factor (MORF) family proteins are required for RNA editing in mitochondria and plastids of plants. *Proc Natl Acad Sci USA* 2012;109:5104-9. https://doi.org/10.1073/pnas.1202452109

63. Bentolila S, Heller WP, Sun T *et al.* RIP1, a member of an Arabidopsis protein family, interacts with the protein RARE1 and broadly affects RNA editing. *Proc Natl Acad Sci USA* 2012;109:E1453-1461. https://doi.org/10.1073/pnas.1121465109

64. Knoop V. Plant PPRs come in multiple flavors-but why? *Plant Cell Physiol* 2020;61:1685-6. https://doi.org/10.1093/pcp/pcaa096

65. Yang Y, Oldenkott B, Ramanathan S *et al.* DYW cytidine deaminase domains have a long-range impact on RNA recognition by the PPR array of chimeric plant C-to-U RNA editing factors and strongly affect target selection. *Plant J* 2023;116:840-54.
https://doi.org/10.1111/tpj.16412

66. Ichinose M, Sugita M. The DYW domains of pentatricopeptide repeat RNA editing factors contribute to discriminate target and non-target editing sites. *Plant Cell Physiol* 2018;59:1652-9.
https://doi.org/10.1093/pcp/pcy086

67. Polson AG, Bass BL. Preferential selection of adenosines for modification by double-stranded RNA adenosine deaminase. *EMBO J* 1994;13:5701-11.
https://doi.org/10.1002/j.1460-2075.1994.tb06908.x

68. Brown GW. The cytidine deaminase APOBEC3C has unique sequence and genome feature preferences. *Genetics* 2024;122:4.
https://doi.org/10.1093/genetics/iyae092

69. Hou S, Lee JM, Myint W *et al.* Structural basis of substrate specificity in human cytidine deaminase family APOBEC3s. *J Biol Chem* 2021;297:100909.
https://doi.org/10.1016/j.jbc.2021.100909

70. Maeda A, Takenaka S, Wang T *et al.* DYW deaminase domain has a distinct preference for neighboring nucleotides of the target RNA editing sites. *Plant J* 2022;111:756-67.
https://doi.org/10.1111/tpj.15850

71. Takenaka M, Zehrmann A, Brennicke A *et al.* Improved computational target site prediction for pentatricopeptide repeat
RNA editing factors. *PLoS One* 2013;8:e65343.
https://doi.org/10.1371/journal.pone.0065343

72. Lenz H, Hein A, Knoop V. Plant organelle RNA editing and its specificity factors: enhancements of analyses and new database features in PREPACT 3.0. *BMC Bioinformatics* 2018;19:255.
https://doi.org/10.1186/s12859-018-2244-9

73. Okuda K, Shoki H, Arai M *et al*. Quantitative analysis of motifs contributing to the interaction between PLS-subfamily members and their target RNA sequences in plastid RNA editing. *Plant J* 2014;80:870–82. https://doi.org/10.1111/tpj.12687

74. Chen L, Liu YG. Male sterility and fertility restoration in crops. *Annu Rev Plant Biol* 2014;65:579–606.
https://doi.org/10.1146/annurev-arplant-050213-040119