# The Analysis of the Editing Defects in the dyw2 Mutant Provides New Clues for the Prediction of RNA Targets of Arabidopsis E+-Class PPR Proteins

Bastien Malbert, Matthias Burger, Mauricio Lopez Obando, Kévin Baudry, Alexandra Launay-Avon, Barbara Härtel, Daniil Verbitskiy, Anja Jörg, Richard Berthomé, Claire Lurin, et al.

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

## To cite this version:

Bastien Malbert, Matthias Burger, Mauricio Lopez Obando, Kévin Baudry, Alexandra Launay-Avon, et al.. The Analysis of the Editing Defects in the dyw2 Mutant Provides New Clues for the Prediction of RNA Targets of Arabidopsis E+-Class PPR Proteins. Plants, 2020, 9 (2), pp.280. 10.3390/plants9020280 . hal-02628277

HAL Id: hal-02628277

https://hal.inrae.fr/hal-02628277

Submitted on 26 May 2020

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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Article

# The Analysis of the Editing Defects in the *dyw2* Mutant Provides New Clues for the Prediction of RNA Targets of Arabidopsis E+-Class PPR Proteins

Bastien Malbert<sup>1,2,†</sup>, Matthias Burger<sup>3,†</sup>, Mauricio Lopez-Obando<sup>1,2,†,‡</sup>, Kevin Baudry<sup>1,2</sup>,
Alexandra Launay-Avon<sup>1,2</sup>, Barbara Härtel<sup>3</sup>, Daniil Verbitskiy<sup>3</sup>, Anja Jörg<sup>3</sup>,

Richard Berthomé<sup>1,2,§</sup>, Claire Lurin<sup>1,2,*</sup>①, Mizuki Takenaka<sup>3,4</sup> and Etienne Delannoy<sup>1,2</sup>①

Received: 31 January 2020; Accepted: 19 February 2020; Published: 21 February 2020

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

check for
updates

**Abstract:** C to U editing is one of the post-transcriptional steps which are required for the proper expression of chloroplast and mitochondrial genes in plants. It depends on several proteins acting together which include the PLS-class pentatricopeptide repeat proteins (PPR). DYW2 was recently shown to be required for the editing of many sites in both organelles. In particular almost all the sites associated with the E+ subfamily of PPR proteins are depending on DYW2, suggesting that DYW2 is required for the function of E+-type PPR proteins. Here we strengthened this link by identifying 16 major editing sites controlled by 3 PPR proteins: OTP90, a DYW-type PPR and PGN and MEF37, 2 E+-type PPR proteins. A re-analysis of the DYW2 editotype showed that the 49 sites known to be associated with the 18 characterized E+-type PPR proteins all depend on DYW2. Considering only the 288 DYW2-dependent editing sites as potential E+-type PPR sites, instead of the 795 known editing sites, improves the performances of binding predictions systems based on the PPR code for E+-type PPR proteins. However, it does not compensate for poor binding predictions.

**Keywords:** *Arabidopsis*; genome containing organelles; mutants; pentatricopeptide repeat (PPR) proteins; RNA editing

## 1. Introduction

The bacterial ancestry of plastids and mitochondria can be seen in their gene expression systems which rely on numerous post-transcriptional regulatory steps [1]. Among these, RNA editing converts selected cytidine (C) to uridine (U) at the RNA level [2,3]. In terrestrial plants, hundreds or thousands of C-to-U editing sites have been identified in plastids and mitochondria [4–6]. The biological significance of this mechanism and of these sites is still not entirely understood but the associated molecular

<footnote n="1"><sup>1</sup> Institute of Plant Sciences Paris-Saclay (IPS2), Université Paris-Saclay, CNRS, INRAE, Univ Evry, 91405 Orsay, France; bastien.malbert@ips2.universite-paris-saclay.fr (B.M.); mauricio.obando@vedascii.org (M.L.-O.); kevin.baudry@ips2.universite-paris-saclay.fr (K.B.); alexandra.launay-avon@inrae.fr (A.L.-A.); etienne.delannoy@inrae.fr (E.D.); richard.berthome@inrae.fr (R.B.)</footnote>

<footnote n="2"><sup>2</sup> Institute of Plant Sciences Paris-Saclay (IPS2), Université de Paris, CNRS, INRAE, 91405 Orsay, France</footnote>

<footnote n="3"><sup>3</sup> Molekulare Botanik, Universität Ulm, 89069 Ulm, Germany; matthias.burger@uni-ulm.de (M.B.); mail@barbara-haertel.de (B.H.); daniil.verbitskiy@me.com (D.V.); anja.joerg@gmx.net (A.J.); mizuki.takenaka@pmg bot.kyoto-u.ac.jp (M.T.)</footnote>

<footnote n="4"><sup>4</sup> Department of Botany, Graduate School of Science, Kyoto University, Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan</footnote>

<footnote>* Correspondence: claire.lurin@inrae.fr; Tel.: +33-(0)1-6915-7593</footnote>

<footnote>† These authors contributed equally to this work.</footnote>

<footnote>‡ Present address: Vedas Investigación e Innovación, Centro empresarial Puerto Seco, Medellín, Colombia.</footnote>

<footnote>§ Present address: Laboratoire des Interactions Plantes-Microorganismes (LIPM), INRA, CNRS, Université de Toulouse, 31326 Castanet-Tolosan, France.</footnote>
mechanisms have been widely studied in the last 30 years (for a review see [7]). The editing of each of these sites requires the recognition of its upstream cis-elements by a specific pentatricopeptide repeat (PPR) protein [8]. PPR proteins are characterized by a tandem array of degenerate 35 amino acid motifs [9] and are almost exclusively targeted to organelles [10]. The PPR proteins are sequence specific RNA binding proteins through the interaction of each PPR motif with one RNA base [11]. The specificity is driven by hydrogen bonds between two amino acids of the PPR motif and the RNA base [12]. For editing, the PPR proteins are binding upstream of the site they are associated to in a way that the last C-terminal PPR motif binds the base at position -4 [13]. Using the list of editing sites experimentally associated to characterized editing PPR proteins, a PPR RNA recognition code was discovered, refined and experimentally validated [14–16]. It can be used to predict the binding sites of uncharacterized PPR proteins [17]. Besides PPR proteins, several other proteins have been shown to be involved in organellar RNA editing suggesting the existence of protein complexes which were named editosomes [7]. In particular, members of a small family called Multiple Organellar RNA editing Factors (MORF/RIP) were found to be involved in RNA editing [18–20]. These proteins function as general editing factors, can interact with many PPR proteins [20,21] and were shown to increase PPR affinity for their target RNA [22].

The PPR proteins can be divided into two classes: the P-type PPR proteins and the PLS-type PPR proteins which include long (L) and short (S) variant of the canonical PPR motif [23,24]. P-type PPR proteins have been associated with a wide range of organellar RNA metabolisms [25] while PLS-type PPR proteins are almost only involved in C-to-U editing [26]. Most PLS-type PPR proteins also harbor conserved C-terminal domains: E1, E2, E+ and DYW which are used to classified PLS-type PPR proteins into corresponding subgroups [24]. While the function of the E1 and E2 domains remains unclear, the DYW domain shows strong similarities to cytidine deaminases [27]. Targeted mutagenesis experiments [28,29] as well as phylogenetic correlations [27] and the observation of C-to-U editing in *Escherichia coli* expressing a plant DYW-type PPR protein [30] strongly support that the DYW domain provides the cytidine deaminase activity required for the PPR editing factors carrying it. However, many PPR proteins involved in editing are truncated and lack a DYW domain. We and others demonstrated that they can recruit DYW domains in *trans* from members of a small PPR protein subgroup, the DYW1-like proteins [31,32]. These six proteins, which were named after DYW1 the first one to be characterized [31], carry very few PPR domains, no canonical E2 and E+ domains but end with a DYW domain [6]. In particular, DYW2 was shown to be associated with the editing of hundreds of sites in *Arabidopsis thaliana* [6]. The DYW2 dependent editing sites show a strong positive bias to editing sites associated to a particular class of PLS-type PPR proteins, the E+ PPR proteins. These PPR proteins harbor the E1, E2 and E+ domains but lack the DYW domain. The functional association of DYW2 and E+ PPR proteins as well as their physical interactions both *in vitro* and *in vivo* [6,32] suggest that the DYW domain is provided by DYW2 exclusively to E+ PPR proteins.

In this study, we strengthened the link between E+ PPR proteins and DYW2 in *A. thaliana* by the refined analysis the *dyw2* mutant and the identification of 16 editing sites associated to 3 PLS-type PPR proteins, PGN, MEF37 and OTP90. We used this functional association to improve the binding predictions of the *A. thaliana* E+ PPR proteins.

## 2. Results

### 2.1. *The Editing Defects in the dyw2 Mutant Define 2 Classes of RNA Editing Targets*

The RNA-seq analysis of the *dyw2-1* KO mutant identified 392 differentially edited sites compared to the wild type (WT) and 261 were considered to depend on DYW2 [6]. This analysis covered 95 sites associated with known PPR proteins and showed that none of the 47 analyzed sites associated with DYW-PPR proteins depend on DYW2 whereas all of the 35 analyzed sites associated with E+-PPR proteins depend on DYW2. It supported the requirement of DYW2 specifically for the function of E+ PPR proteins. However, 14 editing sites associated with PPR proteins were missing from this study.
Slight modifications of our bio-informatic pipeline (see Material and Methods) allowed the study of all 109 editing sites associated with PPR proteins using the Guillaumot et al. dataset [6]. This re-analysis identified 529 differentially edited sites in *dyw2-1* compared to the WT (Data S1). Among them, 487 were found in the mitochondrial genome and 42 in the chloroplastic genome. It should be noted that no genuine editing modification was found in nuclear encoded transcripts demonstrating the specificity of our improved pipeline (not shown). Among these 529 sites, all 39 sites associated with E+ PPR proteins showed a decrease of editing above 30% (ccmB_428, M18007 associated to CMW1 being the less affected with -31.4% of editing extent variation). Overall, 272 sites are differentially edited in *dyw2-1* with a decrease of editing extent above 25%. They don't include any of the 55 editing sites associated with DYW-PPR proteins but one site associated to MEF19 [33], an E-PPR (ccmB_566, M17869) and one site associated to PPME [34], a P-PPR (nad1_898, M234052) (Figure 1, Data S1).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** The editing defects in the *dyw2* mutant define two classes of RNA editing targets. The number of editing sites associated to a PPR of a given subfamily and depending on DYW2 (differentially edited in *dyw2-1* with a decrease of editing extent above 25%) is shown in black. The number of editing sites associated to a PPR of a given subfamily and independent of DYW2 is shown in white.</figcaption></figure>

To further study the specificity of DYW2 regarding PPR sub-families, we extended the repertoire of characterized editing sites by studying three editing PPR proteins: PGN [35], AT5G08305 and AT1G25360.

## 2.2. Characterization of AT1G56570 (PGN), AT5G08305 (MEF37) and AT1G25360 (OTP90)

PGN (AT1G56570) belongs to the E+ PPR protein subfamily with 13 predicted PPR motifs upstream of E1, E2 and E+ domains and is targeted to the mitochondria [35]. Similarly, AT5G08305 encodes an E+ PPR protein with 12 identified PPR motifs followed by E1, E2 and E+ domains. It is predicted to be targeted to mitochondria (Suba4 [36]) and was named Mitochondrial Editing Factor 37 (MEF37) (Figure 2A).

The PPR protein encoded by AT1G25360 was named ORGANELLAR TRANSCRIPT PROCESSING 90 (OTP90). Containing 15 identified PPR domains upstream of E1, E2 and DYW C-terminal domains (Figure 2A), OTP90 belongs to the PLS-DYW subgroup of the PPR family [23]. OTP90 was shown to be targeted to both organelles after expression of fusions of its N-terminal 100aa peptide as well as the full length protein to the fluorescent DsRed2 protein in transient expression assays [10]. To complete these data, we confirmed the subcellular localization of OTP90 in stable transgenic *Arabidopsis* lines (Figure S1). We also showed that OTP90 interacts with MORF proteins (Figures S2 and S3).

To decipher their molecular function, we identified mutants with T-DNA insertions in each of these three genes. The Salk_141937 line (*pgn*) was used by Laluk et al. [35] to characterize the function of PGN. We confirmed that *pgn* mutants did not show any major macroscopic phenotype in normal
growth conditions (Figure 2B) but are slightly more sensitive to *Botrytis cinerea* infection (Figure S4). An RNA-seq comparison of *pgn* leaves with Col-0 ones in normal growth conditions confirmed the absence of full length *PGN* transcripts in *pgn* and identified 98 differentially expressed genes (DEG; Data S2) but no significant gene ontology enrichment.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** The three PPR proteins and their mutants. (a) Predicted domains of the proteins according to [37]. The positions of the T-DNA insertions in the corresponding genes are shown. (b) Phenotype of 25 day-old mutants compared to Col-0.</figcaption></figure>

Two independent mutant lines were used to characterize the function of MEF37, Salk_058773 (*mef37-1*) and Salk_110420 (*mef37-2*). The analysis of the T-DNA flanking sequences identified an insertion 373 bp and 586bp from the initiation codon respectively. The *mef37-1* mutant did not show any major macroscopic phenotype in normal growth conditions (Figure 2B). Its RNA-seq analysis confirmed the absence of full length *MEF37* transcripts in leaves and identified 34 DEG (Data S2) showing a significant enrichment for genes with the GO term "Response to acid chemical".

The analysis of the T-DNA flanking sequences in the SM_3_29659 line (*otp90*) identified an insertion 806 bp after the initiation codon. The *otp90* mutant did not show any major macroscopic phenotype in normal growth conditions (Figure 2B) but we observed shorter primary roots in *in vitro*-germinated *otp90* mutants and smaller seeds were collected on adult plants (Figure S5). These phenotypes were not observed in OTP90 complemented *otp90* lines (Figure S5). The RNA-seq analysis of *otp90* mutant confirmed the absence of full length *OTP90* transcripts in *otp90* leaves and identified 12 DEG (Data S2) including two out of the three copper/zinc superoxide dismutases (*CSD1* and *CDS2*) which are downregulated in *otp90*.
## 2.3. *pgn*, *mef37* and *otp90* Mutants Are Impaired in Mitochondrial C to U Editing

PGN, MEF37 and OTP90 are proteins belonging to PPR subfamilies involved in RNA editing. We characterized the editotype of *pgn*, *mef37-1* and *otp90* by RNA-seq. Using the same analysis pipeline as for *dyw2*, we were able to quantify the editing rate at 204 to 453 editing sites in plastids and 1305 to 1441 editing sites in mitochondria. Among them, we identified several differentially edited sites in each mutant line (Table 1).

<table><caption>**Table 1.** Differentially edited sites identified by RNA-seq in the *mef37-1*, *pgn* and *otp90* mutants compared to edit extent at these sites in *dyw2* mutant.</caption><thead><tr><th>Mutant</th><th>Position <sup>1</sup></th><th>Site Name <sup>2</sup></th><th>WT <sup>3</sup></th><th>Mutant <sup>4</sup></th><th>ΔEE <sup>5</sup></th><th>Padj. <sup>6</sup></th><th>dyw2 <sup>7</sup></th><th>ΔEE dyw2 <sup>8</sup></th><th>Padj. dyw2 <sup>9</sup></th></tr></thead><tbody><tr><td rowspan="10">mef37-1</td><td>M17869</td><td>ccmB_566</td><td>90.5%</td><td>6.6%</td><td>-93%</td><td>0.000</td><td>13%</td><td>-86%</td><td>0.000</td></tr><tr><td>M17884</td><td>ccmB_551</td><td>96.4%</td><td>1.4%</td><td>-99%</td><td>0.000</td><td>4%</td><td>-96%</td><td>0.000</td></tr><tr><td>M23217</td><td>rps3_1470</td><td>72.4%</td><td>0.0%</td><td>-100%</td><td>0.000</td><td>1%</td><td>-99%</td><td>0.000</td></tr><tr><td>M49473</td><td>atp6_71</td><td>0.6%</td><td>0.0%</td><td>-100%</td><td>0.002</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>M189896</td><td>ccmFc_414</td><td>3.5%</td><td>0.0%</td><td>-100%</td><td>0.002</td><td>0%</td><td>-98%</td><td>0.000</td></tr><tr><td>M215126</td><td>nad4_437</td><td>98.1%</td><td>48.5%</td><td>-51%</td><td>0.000</td><td>71%</td><td>-28%</td><td>0.000</td></tr><tr><td>M219378</td><td>mttb_387</td><td>42.8%</td><td>5.7%</td><td>-87%</td><td>0.000</td><td>1%</td><td>-96%</td><td>0.000</td></tr><tr><td>M308476</td><td>ccmC_179</td><td>95.7%</td><td>0.9%</td><td>-99%</td><td>0.000</td><td>2%</td><td>-97%</td><td>0.000</td></tr><tr><td>M362007</td><td>nad4L_trailer_72</td><td>90.2%</td><td>1.8%</td><td>-98%</td><td>0.000</td><td>2%</td><td>-97%</td><td>0.000</td></tr><tr><td>M362343</td><td>atp4_138</td><td>92.1%</td><td>1.3%</td><td>-99%</td><td>0.000</td><td>5%</td><td>-95%</td><td>0.000</td></tr><tr><td>M362349</td><td>atp4_144</td><td>7.7%</td><td>27.3%</td><td>253%</td><td>0.000</td><td>46%</td><td>331%</td><td>0.000</td></tr><tr><td rowspan="2">pgn</td><td>M8348</td><td>cox2_742</td><td>99.6%</td><td>0.0%</td><td>-100%</td><td>0.000</td><td>11%</td><td>-89%</td><td>0.000</td></tr><tr><td>M165765</td><td>nad6_leader_-73</td><td>84.4%</td><td>0.2%</td><td>-100%</td><td>0.000</td><td>3%</td><td>-97%</td><td>0.000</td></tr><tr><td rowspan="14">otp90</td><td>M17839</td><td>ccmB_596</td><td>81.8%</td><td>23.0%</td><td>-72%</td><td>0.000</td><td>93%</td><td>11%</td><td>0.000</td></tr><tr><td>M18355</td><td>ccmB_80</td><td>83.9%</td><td>5.9%</td><td>-93%</td><td>0.000</td><td>89%</td><td>12%</td><td>0.000</td></tr><tr><td>M59321</td><td>nad1_500</td><td>73.9%</td><td>0.0%</td><td>-100%</td><td>0.000</td><td>87%</td><td>-3%</td><td>1.000</td></tr><tr><td>M191687</td><td>ccmFc_1246</td><td>67.6%</td><td>12.7%</td><td>-81%</td><td>0.000</td><td>77%</td><td>26%</td><td>0.000</td></tr><tr><td>M209816</td><td>nad4_194</td><td>9.6%</td><td>14.7%</td><td>54%</td><td>0.027</td><td>2%</td><td>-72%</td><td>0.000</td></tr><tr><td>M209909</td><td>nad4_111</td><td>68.9%</td><td>84.6%</td><td>23%</td><td>0.000</td><td>92%</td><td>5%</td><td>0.000</td></tr><tr><td>M219657</td><td>mttb_108</td><td>4.8%</td><td>0.2%</td><td>-95%</td><td>0.006</td><td>0%</td><td>-97%</td><td>0.000</td></tr><tr><td>M219668</td><td>mttb_97</td><td>81.6%</td><td>1.9%</td><td>-98%</td><td>0.000</td><td>77%</td><td>-4%</td><td>1.000</td></tr><tr><td>M233590</td><td>matR_1950</td><td>3.1%</td><td>7.8%</td><td>149%</td><td>0.049</td><td>0%</td><td>-94%</td><td>0.000</td></tr><tr><td>M308481</td><td>ccmC_184</td><td>83.3%</td><td>5.6%</td><td>-93%</td><td>0.000</td><td>87%</td><td>4%</td><td>1.000</td></tr><tr><td>M329728</td><td>cox3_trailer_248</td><td>0.6%</td><td>0.0%</td><td>-92%</td><td>0.000</td><td>2%</td><td>257%</td><td>0.000</td></tr><tr><td>M362007</td><td>nad4L_trailer_72</td><td>86.9%</td><td>92.6%</td><td>7%</td><td>0.000</td><td>2%</td><td>-97%</td><td>0.000</td></tr></tbody></table>

In *mef37-1*, 11 mitochondrial sites were found differentially edited. The edition of 10 sites was inhibited from 51 to 100%. Seven of them (*ccmB_566*, *ccmB_551*, *rps3_1470*, *mttb_387*, *ccmC_179*, *nad4L_trailer_72* and *atp4_138*) being edited with a rate of less than 7% in the mutant when they were 70-99% edited in wild type leaves. These sites include M17869 (*ccmB_566*) which is known to be edited by MEF19. In contrast, editing of *nad4_437* was decreased by only 50% suggesting that MEF37 is not required for this site and that another editing factor is also able to target this site. Interestingly, in addition to these sites with decreased editing in the *mef37-1* mutant, the edition of M362349 (*atp4_144*) was shown to be strongly increased from 7.7% to 27.3% in *mef37-1*. Except M49473 (*atp6_71*) which was not detected in the *dyw2* analysis because its editing is hardly detectable, all the editing defects of *mef37-1*, including the increase at *atp4_144*, were shown to be similarly impacted by the *dyw2* mutation (Table 1).

Only two mitochondrial editing sites were shown to be affected in *pgn* mutants. Editing of M8348 (*cox2_742*) and M165765 (*nad6_leader_-73*) was completely abolished with -100% of ΔEE. These 2 sites are similarly affected in the *dyw2* mutant with 11% and 3% editing, respectively (Table 1).

<footnote n="1">¹ M stands for "mitochondrion" followed by the position of the site in the mitochondrial genome of Col-0. ² the name of the transcript followed by the position of the site after the start codon except when "trailer" or "leader" is mentioned in which case it is the position after the stop codon or before the start codon respectively. ³ editing extent in Col-0. ⁴ editing extent in the mutant. ⁵ variation of editing extent between the mutant and Col-0. ⁶ P-value adjusted with the Bonferroni correction for the mutant/Col-0 comparison. Values below 0.05 are considered significant.</footnote>

<footnote n="7">⁷ editing extent in *dyw2*. ⁸ variation of editing extent between *dyw2* and Col-0. ⁹ Adjusted P-value for the *dyw2/Col-0* comparison. ND: not detected. The values for the detected editing sites are given in Data S1.</footnote>
Finally, the detailed editotype analysis of the *otp90* mutant by RNA-seq detected 12 mitochondrial sites differentially edited in *otp90*. The edition of seven of them was strongly inhibited (from 72 to 100%). These sites include six sites highly edited in Col-0: M59321 (*nad1_500*), M18355 (*ccmB_80*), M308481 (*ccmC_184*), M219668 (*mttB_97*), M17839 (*ccmB_596*) and M191687 (*ccmFc_1246*). On the other hand, similarly to what was observed for the *atp4_144* site in *mef37-1*, the edition of four sites was increased (from 7 to 149%). None of the six major editing defects of *otp90* mutant were observed in *dyw2*.

To confirm the molecular function of the three PPR proteins, the corresponding mutants were complemented by overexpression of ORFs under the 35S promotor (MEF37 and PGN) or expression of the complete genomic locus (OTP90). As described previously, the OTP90 complemented *otp90* lines did not show the root and seed phenotypes observed in the mutants (Figure S5). The editotypes of the three mutants and their corresponding complemented lines were analyzed by direct cDNA sequencing of the major editing sites identified during the RNAseq analysis (Table 2, Figures S6–S8).

<caption>**Table 2.** Mitochondrial C to U editing events analyzed in each mutant by direct cDNA sequencing. Editing levels are estimated for each gene in percentage related to the peak sizes (shown in Figures S6–S8). Arabidopsis mitochondrial nucleotide editing events analyzed are indicated by their position in the transcript sequence.</caption>

All major RNA editing defects observed in RNAseq experiments were confirmed by Sanger sequencing and peak height analysis (Table 2). All these defects (8 in *mef37-2*, 2 in *pgn* and 6 in *otp90*) were restored in the corresponding complemented lines (Table 2). In particular, 6, 2 and 4 editing sites were abolished in *mef37*, *pgn* and *otp90* mutants, respectively.

<table><thead><tr><th>Mutant</th><th>Position</th><th>Site Name</th><th>WT <sup>1</sup></th><th>Mutant <sup>2</sup></th><th>Compl <sup>3</sup></th></tr></thead><tbody><tr><td rowspan="9">mef37-2</td><td>M17869</td><td>ccmB_566</td><td>100%</td><td>&lt;5%</td><td>100%</td></tr><tr><td>M17884</td><td>ccmB_551</td><td>100%</td><td>&lt;5%</td><td>100%</td></tr><tr><td>M23217</td><td>rps3_1470</td><td>85%</td><td>&lt;5%</td><td>100%</td></tr><tr><td>M215126</td><td>nad4_437</td><td>100%</td><td>70%</td><td>100%</td></tr><tr><td>M219378</td><td>mttb_387</td><td>40%</td><td>&lt;5%</td><td>50%</td></tr><tr><td>M308476</td><td>ccmC_179</td><td>100%</td><td>&lt;5%</td><td>100%</td></tr><tr><td>M362343</td><td>atp4_138</td><td>100%</td><td>&lt;5%</td><td>100%</td></tr><tr><td>M362349</td><td>atp4_144</td><td>&lt;5%</td><td>35%</td><td>&lt;5%</td></tr><tr><td rowspan="2">pgn</td><td>M8348</td><td>cox2_742</td><td>85-90%</td><td>40-45%</td><td>85-90%</td></tr><tr><td>M165765</td><td>nad6_leader_-73</td><td>85-90%</td><td>15%</td><td>85-90%</td></tr><tr><td rowspan="7">otp90</td><td>M17839</td><td>ccmB_596</td><td>100%</td><td>30-35%</td><td>100%</td></tr><tr><td>M18355</td><td>ccmB_80</td><td>85-90%</td><td>&lt;5%</td><td>75%</td></tr><tr><td>M59321</td><td>nad1_500</td><td>100%</td><td>&lt;5%</td><td>100%</td></tr><tr><td>M191687</td><td>ccmFc_1246</td><td>55%</td><td>15%</td><td>90%</td></tr><tr><td>M219668</td><td>mttb_97</td><td>80%</td><td>&lt;5%</td><td>80-85%</td></tr><tr><td>M308481</td><td>ccmC_184</td><td>70%</td><td>&lt;5%</td><td>85%</td></tr></tbody></table>

## 2.4. Using the dyw2 Editing Defects to Refine the PPR Code for PLS PPR Proteins

With the molecular characterization of OTP90, MEF37 and PGN, 60 of the 205 PLS-type PPR proteins encoded in the genome of *Arabidopsis thaliana* have now been characterized and associated to 132 major editing sites both in plastids and mitochondria. In particular, 18 of the 59 E+ PPR proteins have now been associated to 49 major editing sites and all of them are DYW2-dependent reinforcing a model where the E+ PPR proteins require DYW2 for their function. This model can then be used to improve the predictions for uncharacterized E+ PPR proteins. The predictions are using a PPR code but also a list of potential binding sites. There are currently 795 identified editing sites corresponding to as many potential binding sites (Table S1). Given the previous results, for E+ PPR proteins, this list can be reduced to the 272 editing sites depending on DYW2. Using the binding site

<footnote n="1">¹ Editing level in Col-0; ² editing level in the mutant; ³ editing level in a corresponding complementation line. The symbol (-) means that no significant differences were found in mutant compared to Col-0.</footnote>
predictions of Kobayashi et al. [17] for the 18 characterized E+ PPR proteins showed that filtering for DYW2-dependent sites strongly improved the ranks of their associated editing sites (Table 3).

<table><caption>**Table 3.** Prediction of binding on their associated editing sites for the 18 characterized E+ PPR proteins using the system published in [17].</caption><thead><tr><th>PPR</th><th>Position <sup>1</sup></th><th>Site Name <sup>2</sup></th><th>Rank <sup>3</sup></th><th>Rank DYW2 <sup>4</sup></th></tr></thead><tbody><tr><td>AEF1</td><td>M26928</td><td>nad5_1580</td><td>1</td><td>1</td></tr><tr><td>AEF1</td><td>P12707</td><td>atpF_92</td><td>18</td><td>6</td></tr><tr><td>AHG11</td><td>M215187</td><td>nad4_376</td><td>1</td><td>1</td></tr><tr><td>CLB19</td><td>P69942</td><td>clpP_559</td><td>2</td><td>2</td></tr><tr><td>CLB19</td><td>P78691</td><td>rpoA_200</td><td>5</td><td>5</td></tr><tr><td>COD1</td><td>M6961</td><td>cox2_698</td><td>6</td><td>2</td></tr><tr><td>COD1</td><td>M6516</td><td>cox2_253</td><td>8</td><td>4</td></tr><tr><td>COD1</td><td>M209881</td><td>nad4_1129</td><td>34</td><td>14</td></tr><tr><td>CRR21</td><td>P116785</td><td>ndhD_383</td><td>1</td><td>1</td></tr><tr><td>CWM1</td><td>M235780</td><td>nad5_598</td><td>2</td><td>1</td></tr><tr><td>CWM1</td><td>M18007</td><td>ccmB_428</td><td>4</td><td>2</td></tr><tr><td>CWM1</td><td>M308760</td><td>ccmC_463</td><td>12</td><td>6</td></tr><tr><td>GRS1</td><td>M165940</td><td>nad6_103</td><td>1</td><td>1</td></tr><tr><td>GRS1</td><td>M361691</td><td>nad4L_55</td><td>2</td><td>2</td></tr><tr><td>GRS1</td><td>M160356</td><td>rps4_377</td><td>33</td><td>14</td></tr><tr><td>GRS1</td><td>M83057</td><td>nad1_265</td><td>109</td><td>41</td></tr><tr><td>MEF12</td><td>M235556</td><td>nad5_374</td><td>1</td><td>1</td></tr><tr><td>MEF13</td><td>M189532</td><td>ccmFc_50</td><td>1</td><td>1</td></tr><tr><td>MEF13</td><td>M189897</td><td>ccmFc_415</td><td>2</td><td>2</td></tr><tr><td>MEF13</td><td>M161857</td><td>nad2_59</td><td>3</td><td>3</td></tr><tr><td>MEF13</td><td>M215405</td><td>nad4_158</td><td>4</td><td>4</td></tr><tr><td>MEF13</td><td>M28242</td><td>nad5_1916</td><td>5</td><td>5</td></tr><tr><td>MEF13</td><td>M330460</td><td>cox3_314</td><td>6</td><td>6</td></tr><tr><td>MEF13</td><td>M27013</td><td>nad5_1665</td><td>47</td><td>21</td></tr><tr><td>MEF21</td><td>M330517</td><td>cox3_257</td><td>1</td><td>1</td></tr><tr><td>MEF25</td><td>M83014</td><td>nad1_308</td><td>12</td><td>7</td></tr><tr><td>MEF37</td><td>M215126</td><td>nad4_437</td><td>1</td><td>1</td></tr><tr><td>MEF37</td><td>M17884</td><td>ccmB_551</td><td>3</td><td>2</td></tr><tr><td>MEF37</td><td>M17869</td><td>ccmB_566</td><td>12</td><td>5</td></tr><tr><td>MEF37</td><td>M23217</td><td>rps3_1470</td><td>17</td><td>8</td></tr><tr><td>MEF37</td><td>M308476</td><td>ccmC_179</td><td>46</td><td>22</td></tr><tr><td>MEF37</td><td>M362343</td><td>atp4_138</td><td>57</td><td>24</td></tr><tr><td>MEF37</td><td>M219378</td><td>mttB_387</td><td>70</td><td>31</td></tr><tr><td>OTP72</td><td>M23724</td><td>rpl16_440</td><td>1</td><td>1</td></tr><tr><td>OTP80</td><td>P86055</td><td>rpl23_89</td><td>1</td><td>1</td></tr><tr><td>PGN</td><td>M8348</td><td>cox2_742</td><td>1</td><td>1</td></tr><tr><td>PGN</td><td>M165765</td><td>nad6 leader_73</td><td>2</td><td>2</td></tr><tr><td>SLG1</td><td>M288290</td><td>nad3_250</td><td>1</td><td>1</td></tr><tr><td>SLO1</td><td>M215114</td><td>nad4_449</td><td>1</td><td>1</td></tr><tr><td>SLO1</td><td>M24992</td><td>nad9_328</td><td>2</td><td>2</td></tr><tr><td>SLO2</td><td>M241512</td><td>nad7_739</td><td>2</td><td>1</td></tr><tr><td>SLO2</td><td>M219621</td><td>mttB_144</td><td>4</td><td>2</td></tr><tr><td>SLO2</td><td>M361746</td><td>nad4L_110</td><td>76</td><td>23</td></tr><tr><td>SLO2</td><td>M219099</td><td>mttB_666</td><td>139</td><td>47</td></tr><tr><td>SLO2</td><td>M219620</td><td>mttB_145</td><td>177</td><td>65</td></tr></tbody></table>

Table S2 provides the binding predictions extracted from Kobayashi et al. [17] for the 59 A. thaliana E+ PPR proteins to the DYW2 dependent sites. As at least two DYW2 dependent sites (M17869

<footnote n="1">¹ M stands for "mitochondrion" followed by the position of the site in the mitochondrial genome of Col-0. P stands for "plastid". ² the name of the transcript followed by the position of the site after the start codon except when "leader" is mentioned in which case it is the position before the start codon. ³ rank of the corresponding editing site using the full list of potential binding sites. ⁴ rank using only the DYW2-dependent sites.</footnote>
and M234052) are not associated to E+ PPR proteins (Figure 1), it is not possible to reduce the list of candidate PPR proteins binding to DYW2 dependent sites. For the same reason, it is not possible to reduce the list of potential editing sites for non E+ PPR proteins to DYW2 independent sites. However, as the characterized DYW2-independent sites are so far exclusively associated to non E+ PPR proteins (Figure 1), it is possible to reduce the list of candidate PPR proteins binding to DYW2-independent sites by removing the E+ PPR proteins (Table S3).

### 3. Discussion

A strong positive bias for sites associated to PPR-E+ proteins was recently identified in the RNAseq analysis of the *dyw2* mutant affected in a member of a small subfamily of PPR-DYW proteins [6]. Here, we strengthened this correlation by reanalyzing the *dyw2* dataset including 14 previously PPR-associated sites which were missed in the first analysis and by identifying the editing sites controlled by three new PPR proteins: OTP90, a DYW-type PPR and PGN and MEF37, and two E+-type PPR proteins. OTP90, PGN and MEF37 were shown to be mitochondrial editing PPR proteins involved in the editing of a total of 16 major editing sites.

#### 3.1. OTP90, PGN and MEF37 Are Major Mitochondrial Editing Factors

*pgn* mutants were previously shown to be more sensitive to *Botrytis cinerea* infection [35]. Here, we identify PGN molecular function: this E+-type PPR is required for the editing of two mitochondrial sites, M8348 (*cox2_742*) and M165765 (*nad6Leader*-73). The impacts of these editing sites on COX2 and NAD6 proteins as well as the activities of the mitochondrial complexes they belong, cytochrome c oxidase and complex I, respectively, were not characterized. Accordingly, it is not possible to speculate on the causality link between PGN molecular function and its involvement in pathogen sensitivity.

As PGN, MEF37 is a canonical E+-type PPR. While *mef37* mutants did not show any macroscopic phenotype, our RNAseq analysis shows that MEF37 is a major mitochondrial editing factor which is required for editing of 6 sites, M17869 (*ccmB_566*), M17884 (*ccmB_551*), M23217 (*rps3_1470*), M219378 (*mttb_387*), M308476 (*ccmC_179*) and M362343 (*atp4_138*), and which is also involved in editing of M362349 (*atp4_144*) and M215126 (*nad4_437*). Noticeably, two couples of these sites have very close positions (M362343/M362349 and M17869/M17884) questioning the binding site positions of MEF37 and the putative links between the corresponding editing defects (see below).

The mutation of OTP90, a DYW-type PPR, affects root growth at seedling stage and seed size. At the molecular scale, OTP90 is involved in the editing of 8 mitochondrial sites. Indeed, in *otp90*, editing is nearly abolished at 6 sites, M18355 (*ccmB_80*), M59321 (*nad1_500*), M219668 (*mttb_97*) and M17839 (*ccmB_596*), and strongly affected at M191687 (*ccmFc_1246*) and M308481 (*ccmC_184*). While the protein is dually targeted to mitochondria and plastids, we did not observe any defect in plastidial editing. We cannot exclude that OTP90 has another molecular function in these organelles. However, all DYW-type PPR identified so far are editing factors [26] with the exception of CRR2 which is involved in RNA cleavage in plastids [38] and our RNAseq analysis of *otp90* did not identify any defect in plastidial RNA. This absence of identified function of a PPR protein in a cell compartment where it is localized was already shown for OTP87 [39]. These silent editing factors could participate in the puzzling adaptability of the PPR family to organellar genome mutations especially at the RNA editing level.

#### 3.2. OTP90, MEF37 and PGN Share Editing Sites with Several Other Editing Factors

Among the 16 major editing sites affected in *otp90*, *mef37* or *pgn* mutants, several were previously identified in other mutants. This is the case of all the editing sites affected in *mef37* or *pgn* mutants which are also differentially edited in *dyw2* (Table 1). This is also the case of the six major targets of OTP90 which are affected in *morf1* and *morf8/rip1* mutants [18,20].

Sharing a common editing site may reflect the need for the interaction of several editing factors to achieve editing. This was shown for CRR4 and DYW1 proteins which physically interact in order to
edit plastidial ndhD_2 site [31] and, similarly, for DYW2 and E+-type PPR proteins such as CLB19 or SLO2 [6,32]. In this case, DYW1 and DYW2 proteins are proposed to bring a cytidine deaminase activity to the editosomes.

MORF proteins have been reported to bind PLS class PPR proteins to enhance their affinity to target RNAs. Direct binding of MORF1 to the OTP90 *in vivo* and in yeast is consistent with the editing deficient phenotypes. Similarly, MORF1 may be required for increasing affinity between OTP90 and its target RNA sequence. It is unclear how MORF8 is involved in the RNA editing at the target sites of OTP90, since it did not bind to OTP90 at least in yeast. MORF1 may intermediate or support interaction between OTP90 and MORF8 as shown that MORF1 enhances interaction between MEF13 and MORF3 proteins [40].

In contrast, sharing a common editing site may also reflect a competition for site binding between different RNA binding proteins as well as indirect effects of mutations on binding sites of other RNA binding proteins. M17869 (*ccmB_566*) is affected in *mef37* as well as in *mef19*, a previously studied PPR mutant [33]. The PPR protein MEF19, encoded by the gene At3g05240, has 13 PPR repeats and is a member of the E-PPR proteins. In the *mef19* mutant line, RNA editing at M17869 (*ccmB_566*) was completely absent indicating that MEF19 is the single responsible specificity factor for this site and no other PPR protein can take its place for *cis*-element recognition [33]. The nucleotide at M17884 (*ccmB_551*), which is one of the targets of MEF37, is bound by the second PPR motif of MEF19 according to the PPR binding model. The two key amino acid residues for the PPR code, positions 6 and 1' of the second PPR motif are an N and a D, respectively, which bind both U and C but preferentially to a U according to the PPR code. Therefore, transcripts of *ccmB* with an unedited C at position M17884 (*ccmB_551*) should be bound less efficiently by MEF19 than the edited version. These results strongly suggest that the RNA edition at M17884 (*ccmB_551*) is required for the binding of MEF19 and editing at M17869 (*ccmB_566*). Abolishment of an editing event leading to a lack of the downstream edition was also reported in *Physcomitrella patens*. While knock-out of PpPPR_71 gene completely lost RNA editing at *ccmFc-C122* [41], knock-down lines of PpPPR_65 reduced both *ccmFc-C103* and *ccmFc-C122* [42], suggesting editing at *ccmFc-C103* is required for the downstream site editing. Among the 795 putative *cis*-elements of editing sites in *Arabidopsis* organelles (pos -31 to -4 of known editing sites) 323 contain other editing sites which can influence the affinity of the PPR proteins as shown for the binding of MEF19 for the editing of M17869 (*ccmB_566*). Similarly, M362349 (*atp4_144*) is more edited in *mef37* while the edition of M362343 (*atp4_138*) which lies in the binding site of M362349 (*atp4_144*) is almost abolished. It suggests a higher affinity of the unknown PPR editing M362349 for the unedited M362343. Such secondary or even multiple sequential effects for downstream sites may be commonly present especially in mitochondria where many editing sites are closely located each other. This is somehow limiting reverse genetics based identification of PPR targets because in some cases it becomes impossible to separate the primary direct effect of a *ppr* mutant on its RNA target from secondary effects of sequence modification of the RNA target of another RNA binding protein.

This limitation can be exemplified by the study of PPME. PPME is so far one of the three non PLS-type PPR proteins involved in RNA editing. It is required for the editing of M234052 (*nad1_898*) and M234091 (*nad1_937*) but not M234082 (*nad1_928*) which stands between the 2 sites. In parallel, PPME was shown to bind M234052 (*nad1_898*) and M234082 (*nad1_928*) sites and not the downstream M234091 (*nad1_937*) site which is affected by its mutation [34]. Without further experiments, it remains unclear how this P-type PPR protein participates in the editing of mitochondrial transcripts.

### 3.3. Specificity of DYW2 for E+ PPR Proteins

All editing sites associated to E+ PPR so far are dependent on DYW2 for their edition but among the DYW2 dependent sites, M17869 (*ccmB_566*) is associated to MEF19, an E-PPR [33] and M234052 (*nad1_898*) is associated to PPME, a pure PPR [34].
As mentioned in the previous section, it is most likely that the impairment of M17869 (*ccmB_566*) edition in an E+ PPR mutant *mef37* and *dyw2* is a secondary effect of the lack of the direct target editing at M17884 (*ccmB_551*) which prevents the binding of MEF19.

As also previously discussed, if PPME is required for the editing of M234052 (*nad1_898*) and M234091 (*nad1_937*), its molecular function in editing is not yet understood. In particular, the PPR code is not well conserved in this P-type PPR protein and it is unclear whether its function is to replace a canonical PLS-type protein or whether, like NUWA, it performs an unknown function in editing [34]. Moreover, only one of its two sites (M234052) is DYW2-dependent suggesting that DYW2 is not fundamentally required for the editing function of PPME. It is possible that M234052 (*nad1_898*) is edited by an unknown E+ PPR which requires both PPME and DYW2 to work. The best candidate PPR proteins for targeting this site according to the PPR code are two E+ PPR, At3g28640 and At3g28660 [17].

As a consequence the DYW2 criterion cannot be used to reduce the number of potential binding sites for other PLS-type PPR because it is not possible to consider that editing sites associated with non E+ PPR proteins are exclusively non DYW2-dependent.

On the other hand, the DYW2 criterion is useful to limit the number of candidate targets of E+ PPR proteins. However, it needs to be noted that some E+ editing sites are not completely abolished in the *dyw2* mutant while they are in the corresponding E+ knock-out mutant. One explanation could be that the DYW2 protein was still present as traces in the *dyw2* mutant because of the use of the ABI3 promoter to bypass the embryo-lethality [6]. These traces could be sufficient to sustain the partial editing of a few sites. However, no DYW2 transcript was detected by RNA-seq in these plants. It does not rule out the presence of traces of DYW2 but makes it unlikely. There are 6 DYW1-like proteins (DYW-PPR with very few PPR domains, no PG box and lack of conserved residues in the E2 domain): DYW1, DYW2, MEF8, MEF8s, AT2G34370 (DYW3) and AT1G29710 (DYW4). DYW1 is associated to CRR4 [31] but MEF8, MEF8s, DYW3 and/or DYW4 could be associated to some E+ PPR proteins. The editotype of the *mef8* mutant was characterized [43]. About a third of the MEF8-dependent sites are also DYW2-dependent and in these cases the editing is only partially abolished in either mutants or both (Table S4).

### 3.4. Limitations of the Code for PLS-Type PPR Proteins

The current PPR codes are almost exclusively based on the known editing sites/PPR associations which allow deriving the binding preferences of each PPR motif [14,15,17,44]. This analysis relies on the assumptions that editing PPR proteins bind from -4 and upstream of the sites which edition is significantly reduced in the corresponding KO mutant, one PPR motif per base and that the RNA binding site is not edited.

While the binding of PLS type PPR proteins to the nucleotides located in position -4 and upstream of the editing site, one base per PPR motif seems demonstrated by several studies [13,45], care should be taken when identifying PPR binding sites through the reduction of editing in the corresponding KO mutant. First, it was shown that the proper binding of PLS PPR to its target RNA does not always result in editing [46–48] meaning that this approach can miss true binding sites. Second, as previously discussed the reduction of editing in KO mutants can be indirect and thus could not correspond to bona fide binding sites. These indirect effects can be suspected when a reduction in editing is associated to a close-by variation of editing or splicing. Furthermore, the advent of NGS allows not only the detection of completely abolished editing sites but also simple decrease in editing. We chose an arbitrary threshold of 25% of editing decrease to consider a site as dependent of a particular protein. Should these sites be considered as putative binding sites? And if yes, should they be given the same weight as abolished editing sites in the estimation of the binding preferences of each PPR motif? These remain open questions. They could be solved by dedicated RNA-PPR binding experiments such as RNA electromobility shift assay (REMSA) or analysis of PPR footprints on RNAs. However,
these methods remain challenging for editing PPR proteins which probably transiently interact with their targets.

Using the DYW2 criterion for the prediction E+-PPR binding implies that some potential sites are ignored even if they get very good prediction scores. A simple explanation is that current prediction systems are not accurate. However, factors other than PPR motifs can dictate the binding pattern of a PPR. For example, it was shown that RNA secondary structures can impede PPR binding *in vitro* [13,49,50] but at the same time, several PPR proteins are known to bind to sites predicted to form RNA hairpins [51,52] which suggests that other RNA binding proteins are necessary to make these sites available to PPR proteins. On the other hand, PPR proteins are probably competing with other RNA binding proteins (other PPR proteins, CRM proteins ...) so that PPR proteins are probably not binding *in vivo* to all their potential binding sites.

# 4. Materials and Methods

## 4.1. Plant Material, Phenotype Characterization and Complementation Assay

*Arabidopsis thaliana* ecotype Columbia-0 (Col-0) was used as wild type plant. The T-DNA insertional mutants were obtained from the ABRC stock center. T-DNA homozygous plants were selected by PCR genotyping using the primers described in Table S5. The localization of the insertions was verified by sequencing of the PCR products. Plants were grown in chambers at 16 h/8 h light/dark cycles, 22 °C at day and 20 °C at night. For *in vitro* phenotype characterization of *otp90*, seeds were surface sterilized, sown in Murashige and Skoog (MS) solid medium with 0.8% (W/V) of agar-agar type E (Sigma Aldrich, Saint-Louis, MO, USA) and grown in a growth chamber (16-h light/8-h dark cycle, 22 °C, 50% hygrometry) after cold treatment for 48 h at 4 °C. Primary root length elongation was measured on 12 day old plantlets after cold treatment of seeds using EZ-Rhizo software [53].

For the complementation of the *otp90* mutant, a genomic 3595 bp fragment containing the promoter, the Open Reading Frame (ORF) and the terminator region of the *OTP90* gene was amplified using the primers described in Table S5 and cloned into the pGWB1 binary vector [54] by Gateway™ (Invitrogen, Carlsbad, CA, USA). The *pgn* mutant was complemented with the ORF of *PGN* amplified with the primers described in Table S5 and cloned into the pGWB14 binary vector [55] by Gateway™ (Invitrogen, Carlsbad, CA, USA). For the complementation of *mef37-2*, the ORF of the MEF37 gene was amplified using the primers described in Table S5 and cloned into the pMpGWB102 binary vector [56] by the In-Fusion cloning system (Takara Bio Europe, Saint-Germain-en-Laye, France).

## 4.2. RNA Analysis

Total RNA was extracted from leaves using the NucleoZOL kit (Macherey-Nagel, Hœrdt, France) followed by a purification with the Agencourt RNA Cleanup XP beads (Beckmann-Coulter, Villepinte, France). The sequencing libraries were constructed using the TruSeq Stranded Total RNA with Ribo-Zero Plant (Illumina Inc., San Diego, CA, USA) and then sequenced with a Nextseq500 in single-end 75 bp. The RNA-seq data was analysed as described [57] except for the following modifications. The reads were mapped with STAR v2.7.0c [58] on the *A. thaliana* genome available from the release 39 of EnsemblGenomes but the mitochondrial genome was replaced by the Col-0 mitochondrial genome from the NCBI under the accession BK010421. The mapping parameters were -outSAMprimaryFlag AllBestScore -outFilterMultimapScoreRange 0-alignIntronMax 1.

## 4.3. Direct cDNA Sequence Analysis.

For direct cDNA sequencing, relevant cDNA fragments covering one or more editing sites were obtained by RT-PCR amplification according to published protocols [47]. At the RNA editing sites, cDNA sequences were evaluated for the respective C to T differences. RNA editing levels were estimated by the relative heights of the respective nucleotide peaks in the sequence analyses.
Ratios between heights were calculated with the DNA Dynamo program. Sequences were obtained commercially from 4base lab (Reutlingen, Germany) or from Macrogen (Seoul, South Korea).

### 4.4. Subcellular Protein Localization

The OTP90 full-length ORF without stop codon was cloned into p0229-DsRed2 destination vector using Gateway cloning and the primers described in Table S5. After transformation of A. thaliana Col-0 plants, Subcellular localization of the OTP90:DsRed2 fusion protein was observed in root and cotyledon cells of 10–15 days old primary transformants grown on MS solid medium using a spectral Leica SP2 AOBS confocal microscope (Leica Microsystems, Wetzlar, Germany) equipped with an argon laser and a HeNe laser. To confirm the mitochondrial localization, the seedlings were stained with 1 µM of Mito-tracker green marker (Invitrogen) for 15–30 min. Signals were detected using laser lines at 488 nm (Mito-tracker green; excitation/emission 488/510–530 nm) and 543 nm (DsRed2; excitation/emission 543/570–600 nm/chlorophyll autofluorescence; excitation/emission 543/600–700 nm). The images were coded green (Mitotraker green or chlorophyll autofluorescence) and red (DsRed2) giving yellow co-localization signals when green and red signals overlapped in merged images.

### 4.5. *Yeast Two Hybrid Assays*

Sequences coding for OTP90 and the various MORF were cloned with the In-Fusion HD cloning system (Takara Bio Europe, Saint-Germain-en-Laye, France). OTP90 inserts were integrated into the binding domain containing vector pGBKT7, and MORF into the activator domain containing vector pGADT7 of the GAL4 Two Hybrid System 3 (Takara Bio Europe, Saint-Germain-en-Laye, France). The vectors were co-transfected for expression into yeast cells (PJ69-4A) according to the protocol. Yeast cells with both bait and prey vectors were cultured in synthetic dropout medium without Leu and Trp. 5 µL of suspended cells with an OD600 of 0.3 were dropped onto the various selection media plates. To detect strong interactions, 2.5 mM of 3-Amino-1,2,4-triazole was added in the media plates.

### 4.6. *In Vivo Protein-Protein Interaction Assay with BiFC Analysis*

For in vivo protein-protein interaction assays, OTP90, MORF1, MORF2, MORF3, MORF4, MORF5 and MORF9 were fused with either YFP-N (1-155) or YFP-C (156-239) at their C terminus, respectively and cloned under 35S promoter. OTP90 was cloned in pGWB1 binary vector [54] by LR recombination reaction (Invitrogen) whereas MORF proteins were cloned in pMDC123 vector [59] containing a multiple cloning site from pET41 (Merck Millipore Novagen®, Darmstadt, Germany) with In-Fusion HD cloning system (Clontech). Combinations of agrobacteria containing OTP90-YFP-N and each MORF-YFP-C, on one hand and OTP90-YFP-C and each MORF-YFPN, on the other hand, were used at equal concentration to agroinfiltrate 3-weeks old tobacco leaves. After 48 h of incubation at 21 °C, transformed tobacco cells were analyzed using a confocal microscope Leica, TCS SP5 II.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2223-7747/9/2/280/s1, Figure S1: OTP90 is dually targeted to mitochondria and plastids in transgenic Arabidopsis thaliana plants, Figure S2: OTP90 interacts with MORF proteins in Yeast two hybrid, Figure S3: OTP90 interacts with MORF1 in BiFC, Figure S4: The pgn mutant is more sensitive to Botrytis cinerea, Figure S5: The otp90 mutant is impaired in root length and seed size, Figure S6: Complementation of the pgn mutant for RNA editing defects in two mitochondrial sites. Figure S7: Complementation of the mef37-2 mutant for RNA editing defects in eight mitochondrial sites. Figure S8: Complementation of the otp90 mutant for RNA editing defects in six mitochondrial sites, Data S1: Full editotypes of dyw2, pgn, mef37-1 and otp90 mutants by RNA-seq, Data S2: full expression profiles of dyw2, pgn, mef37-1 and otp90 mutants by RNA-seq, Table S1: Mitochondrial and plastidial editing sites identified by RNAseq analysis in Arabidopsis, Table S2: "E+PPR-to-DYW2 editing site" prediction data, Table S3: "non DYW2 Editing site-to-non E+ PPR" prediction data, Table S4: comparison of MEF8-dependent sites with dyw2 results, Table S5: primers used in this study.

**Author Contributions:** Conceptualization, E.D. and C.L.; formal analysis, E.D., B.M., M.L.-O., M.B.; investigation, B.M., K.B., M.L.-O., A.L.-A., M.B., A.J., B.H., D.V.; writing—original draft preparation, E.D., C.L.; writing—review and editing, E.D., B.M., M.L.-O., M.B., C.L., M.T.; supervision, C.L., M.T., E.D., R.B.; funding acquisition, C.L., M.T., E.D.. All authors have read and agreed to the published version of the manuscript.
**Funding:** M.L.-O. and K.B. research were supported by French Ph. D. fellowships from "Ministère de la Recherche et de l'Enseignement Supérieur". B.M. research was supported by a Labex Saclay Plant Sciences-SPS doctoral fellowship. The IPS2 benefits from the support of the Labex Saclay Plant Sciences-SPS (ANR-10-LABX-0040-SPS). M.T. research was supported by the Deutsche Forschungsgemeinschaft (TA621/10-1) and JSPS Grants-in-Aid for Scientific Research (18H02462).

**Conflicts of Interest:** The authors declare no conflict of interest.

**Data Accessibility:** The RNA-seq raw data were deposited at CATdb and GEO under the accession numbers NGS2018_04_Half_Edit, NGS2019_18_Eplus and GSE144560, GSE144674.

## References

1. Small, I.D.; Rackham, O.; Filipovska, A. Organelle transcriptomes: Products of a deconstructed genome. *Curr. Opin. Microbiol.* **2013**, *16*, 652–658. [CrossRef] [PubMed]

2. Gualberto, J.M.; Lamattina, L.; Bonnard, G.; Weil, J.H.; Grienenberger, J.M. RNA editing in wheat mitochondria results in the conservation of protein sequences. *Nature* **1989**, *341*, 660–662. [CrossRef]

3. Hoch, B.; Maier, R.M.; Appel, K.; Igloi, G.L.; Kössel, H. Editing of a chloroplast mRNA by creation of an initiation codon. *Nature* **1991**, *353*, 178–180. [CrossRef] [PubMed]

4. Edera, A.A.; Gandini, C.L.; Sanchez-Puerta, M.V. Towards a comprehensive picture of C-to-U RNA editing sites in angiosperm mitochondria. *Plant Mol. Biol.* **2018**, *97*, 215–231. [CrossRef] [PubMed]

5. Grewe, F.; Herres, S.; Viehöver, P.; Polsakiewicz, M.; Weisshaar, B.; Knoop, V. A unique transcriptome: 1782 positions of RNA editing alter 1406 codon identities in mitochondrial mRNAs of the lycophyte Isoetes engelmannii. *Nucleic Acids Res.* **2011**, *39*, 2890–2902. [CrossRef] [PubMed]

6. Guillaumot, D.; Lopez-Obando, M.; Baudry, K.; Avon, A.; Rigaill, G.; Falcon De Longevialle, A.; Broche, B.; Takenaka, M.; Berthomé, R.; De Jaeger, G.; et al. Two interacting PPR proteins are major Arabidopsis editing factors in plastid and mitochondria. *Proc. Natl. Acad. Sci. USA* **2017**, *114*, 8877–8882. [CrossRef]

7. Small, I.D.; Schallenberg-Rüdinger, M.; Takenaka, M.; Mireau, H.; Ostersetzer-Biran, O. Plant organellar RNA editing: What 30 years of research has revealed. *Plant J.* **2019**. [CrossRef]

8. Kotera, E.; Tasaka, M.; Shikanai, T. A pentatricopeptide repeat protein is essential for RNA editing in chloroplasts. *Nature* **2005**, *433*, 326–330. [CrossRef]

9. Small, I.D.; Peeters, N. The PPR motif—A TPR-related motif prevalent in plant organellar proteins. *Trends Biochem. Sci.* **2000**, *25*, 45–47. [CrossRef]

10. Colcombet, J.; Lopez-Obando, M.; Heurtevin, L.; Bernard, C.; Martin, K.; Berthomé, R.; Lurin, C. Systematic study of subcellular localization of Arabidopsis PPR proteins confirms a massive targeting to organelles. *RNA Biol.* **2013**, *10*, 1557–1575. [CrossRef]

11. Yin, P.; Li, Q.; Yan, C.; Liu, Y.; Liu, J.; Yu, F.; Wang, Z.; Long, J.; He, J.; Wang, H.-W.; et al. Structural basis for the modular recognition of single-stranded RNA by PPR proteins. *Nature* **2013**, *504*, 168–171. [CrossRef]

12. Shen, C.; Wang, X.; Liu, Y.; Li, Q.; Yang, Z.; Yan, N.; Zou, T.; Yin, P. Specific RNA recognition by designer pentatricopeptide repeat protein. *Mol. Plant* **2015**, *8*, 667–670. [CrossRef] [PubMed]

13. Kindgren, P.; Yap, A.; Bond, C.S.; Small, I. Predictable alteration of sequence recognition by RNA editing factors from Arabidopsis. *Plant Cell* **2015**, *27*, 403–416. [CrossRef] [PubMed]

14. Barkan, A.; Rojas, M.; Fujii, S.; Yap, A.; Chong, Y.S.; Bond, C.S.; Small, I. A Combinatorial Amino Acid Code for RNA Recognition by Pentatricopeptide Repeat Proteins. *PLoS Genet.* **2012**, *8*. [CrossRef] [PubMed]

15. Takenaka, M.; Zehrmann, A.; Brennicke, A.; Graichen, K. Improved computational target site prediction for pentatricopeptide repeat RNA editing factors. *PLoS ONE* **2013**, *8*. [CrossRef]

16. Yan, J.; Yao, Y.; Hong, S.; Yang, Y.; Shen, C.; Zhang, Q.; Zhang, D.; Zou, T.; Yin, P. Delineation of pentatricopeptide repeat codes for target RNA prediction. *Nucleic Acids Res.* **2019**, *47*, 3728–3738. [CrossRef]

17. Kobayashi, T.; Yagi, Y.; Nakamura, T. Comprehensive Prediction of Target RNA Editing Sites for PLS-Class PPR Proteins in Arabidopsis thaliana. *Plant Cell Physiol.* **2019**, *60*, 862–874. [CrossRef]

18. Bentolila, S.; Heller, W.P.; Sun, T.; Babina, A.M.; Friso, G.; Van Wijk, K.J.; Hanson, M.R. 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*. [CrossRef]

19. Bentolila, S.; Oh, J.; Hanson, M.R.; Bukowski, R. Comprehensive High-Resolution Analysis of the Role of an Arabidopsis Gene Family in RNA Editing. *PLoS Genet.* **2013**, *9*. [CrossRef]
20. Takenaka, M.; Zehrmann, A.; Verbitskiy, D.; Kugelmann, M.; Hartel, B.; Brennicke, A. 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–5109. [CrossRef]

21. Bayer-Császár, E.; Haag, S.; Jörg, A.; Glass, F.; Härtel, B.; Obata, T.; Meyer, E.H.; Brennicke, A.; Takenaka, M. The conserved domain in MORF proteins has distinct affinities to the PPR and E elements in PPR RNA editing factors. *Biochim. Biophys. Acta Gene Regul. Mech.* **2017**, *1860*, 813–828. [CrossRef]

22. Yan, J.; Zhang, Q.; Guan, Z.; Wang, Q.; Li, L.; Ruan, F.; Lin, R.; Zou, T.; Yin, P. MORF9 increases the RNA-binding activity of PLS-type pentatricopeptide repeat protein in plastid RNA editing. *Nat. Plants* **2017**, *3*, 17037. [CrossRef] [PubMed]

23. Lurin, C.; Andrés, C.; Aubourg, S.; Bellaoui, M.; Bitton, F.; Bruyère, C.; Caboche, M.; Debast, C.; Gualberto, J.; Hoffmann, B.; et al. Genome-wide analysis of Arabidopsis pentatricopeptide repeat proteins reveals their essential role in organelle biogenesis. *Plant Cell* **2004**, *16*, 2089–2103. [CrossRef] [PubMed]

24. Cheng, S.; Gutmann, B.; Zhong, X.; Ye, Y.; Fisher, M.F.; Bai, F.; Castleden, I.; Song, Y.; Song, B.; Huang, J.; 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–547. [CrossRef] [PubMed]

25. Waltz, F.; Nguyen, T.T.; Arrivé, M.; Bochler, A.; Chicher, J.; Hammann, P.; Kuhn, L.; Quadrado, M.; Mireau, H.; Hashem, Y.; et al. Small is big in Arabidopsis mitochondrial ribosome. *Nat. Plants* **2019**, *5*, 106–117. [CrossRef] [PubMed]

26. Barkan, A.; Small, I. Pentatricopeptide Repeat Proteins in Plants. *Annu. Rev. Plant Biol.* **2014**, *65*, 415–442. [CrossRef] [PubMed]

27. Salone, V.; Rüdinger, M.; Polsakiewicz, M.; Hoffmann, B.; Groth-Malonek, M.; Szurek, B.; Small, I.; Knoop, V.; Lurin, C. A hypothesis on the identification of the editing enzyme in plant organelles. *FEBS Lett.* **2007**, *581*, 4132–4138. [CrossRef]

28. Boussardon, C.; Avon, A.; Kindgren, P.; Bond, C.S.; Challenor, M.; Lurin, C.; Small, I. The cytidine deaminase signature HxE(x)nCxxC of DYW1 binds zinc and is necessary for RNA editing of ndhD-1. *New Phytol.* **2014**, *203*, 1090–1095. [CrossRef]

29. Wagoner, J.A.; Sun, T.; Lin, L.; Hanson, M.R. Cytidine deaminase motifs within the DYW domain of two pentatricopeptide repeat-containing proteins are required for site-specific chloroplast RNA editing. *J. Biol. Chem.* **2015**, *290*, 2957–2968. [CrossRef]

30. Oldenkott, B.; Yang, Y.; Lesch, E.; Knoop, V.; Schallenberg-Rüdinger, M. Plant-type pentatricopeptide repeat proteins with a DYW domain drive C-to-U RNA editing in Escherichia coli. *Commun. Biol.* **2019**, *2*. [CrossRef]

31. Boussardon, C.; Salone, V.; Avon, A.; Berthomé, R.; Hammani, K.; Okuda, K.; Shikanai, T.; Small, I.; Lurina, C. Two interacting proteins are necessary for the editing of the ndhD-1 site in Arabidopsis plastids. *Plant Cell* **2012**, *24*, 3684–3694. [CrossRef] [PubMed]

32. Andrés-Colás, N.; Zhu, Q.; Takenaka, M.; De Rybel, B.; Weijers, D.; Van Der Straeten, D. Multiple PPR protein interactions are involved in the RNA editing system in Arabidopsis mitochondria and plastids. *Proc. Natl. Acad. Sci. USA* **2017**, *114*, 8883–8888. [CrossRef] [PubMed]

33. Takenaka, M.; Verbitskiy, D.; Zehrmann, A.; Brennicke, A. Reverse genetic screening identifies five E-class PPR proteins involved in RNA editing in mitochondria of Arabidopsis thaliana. *J. Biol. Chem.* **2010**, *285*, 27122–27129. [CrossRef] [PubMed]

34. Leu, K.-C.; Hsieh, M.-H.; Wang, H.-J.; Hsieh, H.-L.; Jauh, G.-Y. Distinct role of Arabidopsis mitochondrial P-type pentatricopeptide repeat protein-modulating editing protein, PPME, in nad1 RNA editing. *RNA Biol.* **2016**, *13*, 593–604. [CrossRef]

35. Laluk, K.; Abuqamar, S.; Mengiste, T. The Arabidopsis mitochondria-localized pentatricopeptide repeat protein PGN functions in defense against necrotrophic fungi and abiotic stress tolerance. *Plant Physiol.* **2011**, *156*, 2053–2068. [CrossRef]

36. Hooper, C.M.; Castleden, I.R.; Tanz, S.K.; Aryamanesh, N.; Millar, A.H. SUBA4: The interactive data analysis centre for Arabidopsis subcellular protein locations. *Nucleic Acids Res.* **2017**, *45*, D1064–D1074. [CrossRef]

37. Gutmann, B.; Royan, S.; Schallenberg-Rüdinger, M.; Lenz, H.; Castleden, I.R.; McDowell, R.; Vacher, M.A.; Tonti-Filippini, J.; Bond, C.S.; Knoop, V.; et al. The expansion and diversification of pentatricopeptide repeat RNA editing factors in plants. *Mol. Plant* **2019**. [CrossRef]

38. Hashimoto, M.; Endo, T.; Peltier, G.; Tasaka, M.; Shikanai, T. A nucleus-encoded factor, CRR2, is essential for the expression of chloroplast ndhB in Arabidopsis. *Plant J.* **2003**, *36*, 541–549. [CrossRef]
39. Hammani, K.; Des Francs-Small, C.C.; Takenaka, M.; Tanz, S.K.; Okuda, K.; Shikanai, T.; Brennicke, A.; Small, I. The pentatricopeptide repeat protein OTP87 is essential for RNA editing of nad7 and atp1 transcripts in Arabidopsis mitochondria. *J. Biol. Chem.* **2011**, *286*, 21361–21371. [CrossRef]

40. Glass, F.; Härtel, B.; Zehrmann, A.; Verbitskiy, D.; Takenaka, M. MEF13 requires MORF3 and MORF8 for RNA editing at eight targets in mitochondrial mRNAs in Arabidopsis thaliana. *Mol. Plant* **2015**, *8*, 1466–1477. [CrossRef]

41. Tasaki, E.; Hattori, M.; Sugita, M. The moss pentatricopeptide repeat protein with a DYW domain is responsible for RNA editing of mitochondrial ccmFc transcript. *Plant J.* **2010**, *62*, 560–570. [CrossRef] [PubMed]

42. Ichinose, M.; Sugita, C.; Yagi, Y.; Nakamura, T.; Sugita, M. Two DYW subclass PPR proteins are involved in RNA editing of ccmFc and atp9 transcripts in the moss Physcomitrella patens: First complete set of PPR editing factors in plant mitochondria. *Plant Cell Physiol.* **2013**, *54*, 1907–1916. [CrossRef] [PubMed]

43. Diaz, M.F.; Bentolila, S.; Hayes, M.L.; Hanson, M.R.; Mulligan, R.M. A protein with an unusually short PPR domain, MEF8, affects editing at over 60 Arabidopsis mitochondrial C targets of RNA editing. *Plant J.* **2017**, *92*, 638–649. [CrossRef] [PubMed]

44. Yagi, Y.; Hayashi, S.; Kobayashi, K.; Hirayama, T.; Nakamura, T. Elucidation of the RNA recognition code for pentatricopeptide repeat proteins involved in organelle RNA editing in plants. *PLoS ONE* **2013**, *8*, e57286. [CrossRef]

45. Okuda, K.; Shoki, H.; Arai, M.; Shikanai, T.; Small, I.; Nakamura, T. Quantitative analysis of motifs contributing to the interaction between PLS-subfamily members and their target RNA sequences in plastid RNA editing. *Plant J.* **2014**, *80*, 870–882. [CrossRef]

46. Chateigner-Boutin, A.-L.; des Francs-Small, C.C.; Delannoy, E.; Kahlau, S.; Tanz, S.K.S.K.; de Longevialle, A.F.; Fujii, S.; Small, I. OTP70 is a pentatricopeptide repeat protein of the E subgroup involved in splicing of the plastid transcript rpoC1. *Plant J.* **2011**, *65*, 532–542. [CrossRef]

47. Hammani, K.; Takenaka, M.; Miranda, R.; Barkan, A. A PPR protein in the PLS subfamily stabilizes the 5'-end of processed rpl16 mRNAs in maize chloroplasts. *Nucleic Acids Res.* **2016**, *44*, 4278–4288. [CrossRef]

48. Ruwe, H.; Gutmann, B.; Schmitz-Linneweber, C.; Small, I.; Kindgren, P. The E domain of CRR2 participates in sequence-specific recognition of RNA in plastids. *New Phytol.* **2019**, *222*, 218–229. [CrossRef]

49. McDermott, J.J.; 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. [CrossRef]

50. Miranda, R.G.; Rojas, M.; Montgomery, M.P.; Gribbin, K.P.; Barkan, A. RNA-binding specificity landscape of the pentatricopeptide repeat protein PPR10. *RNA* **2017**, *23*, 586–599. [CrossRef]

51. 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–1941. [CrossRef] [PubMed]

52. Prikryl, J.; Rojas, M.; Schuster, G.; Barkan, A. Mechanism of RNA stabilization and translational activation by a pentatricopeptide repeat protein. *Proc. Natl. Acad. Sci. USA* **2011**, *108*, 415–420. [CrossRef]

53. Armengaud, P.; Zambaux, K.; Hills, A.; Sulpice, R.; Pattison, R.J.; Blatt, M.R.; Amtmann, A. EZ-Rhizo: Integrated software for the fast and accurate measurement of root system architecture. *Plant J.* **2009**, *57*, 945–956. [CrossRef] [PubMed]

54. Nakagawa, T.; Suzuki, T.; Murata, S.; Nakamura, S.; Hino, T.; Maeo, K.; Tabata, R.; Kawai, T.; Tanaka, K.; Niwa, Y.; et al. Improved gateway binary vectors: High-performance vectors for creation of fusion constructs in transgenic analysis of plants. *Biosci. Biotechnol. Biochem.* **2007**, *71*, 2095–2100. [CrossRef] [PubMed]

55. Nakagawa, T.; Kurose, T.; Hino, T.; Tanaka, K.; Kawamukai, M.; Niwa, Y.; Toyooka, K.; Matsuoka, K.; Jinbo, T.; Kimura, T. Development of series of gateway binary vectors, pGWBs, for realizing efficient construction of fusion genes for plant transformation. *J. Biosci. Bioeng.* **2007**, *104*, 34–41. [CrossRef] [PubMed]

56. Ishizaki, K.; Nishihama, R.; Ueda, M.; Inoue, K.; Ishida, S.; Nishimura, Y.; Shikanai, T.; Kohchi, T. Development of gateway binary vector series with four different selection markers for the liverwort marchantia polymorpha. *PLoS ONE* **2015**, *10*. [CrossRef] [PubMed]

57. Malbert, B.; Rigaill, G.; Brunaud, V.; Lurin, C.; Delannoy, E. Bioinformatic Analysis of Chloroplast Gene Expression and RNA Posttranscriptional Maturation Using RNA Sequencing. In *Plastids*; Springer: Berlin/Heidelberg, Germany, 2018; Volume 1829.
58. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. *Bioinformatics* **2013**, *29*, 15–21. [CrossRef]

59. Curtis, M.D.; Grossniklaus, U. A Gateway Cloning Vector Set for High-Throughput Functional Analysis of Genes in Planta. *Plant Physiol.* **2003**, *133*, 462–469. [CrossRef]

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

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).