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

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

Although many factors have been identified to function in the C-to-U editing, the exact mechanism remains to be deciphered Small et al., 2020). Pentatricopeptide repeat (PPR) proteins have been identified to play an essential role in editing site recognition. PPR proteins belong to one of the largest protein families in plants, with over 400 members in terrestrial plants (Lurin et al., 2004;O'Toole et al., 2008;Fujii and Small, 2011;Wei and Han, 2016). PPRs are classified into two classes: P-class and PLS-class (Small and Peeters, 2000;Lurin et al., 2004;Cheng et al., 2016). The P-class PPR proteins consist of bona fide P-motifs with 35 amino acids, while the PLS-class proteins harbor P-, L-(35-36 amino acids), and S-motifs (31 amino acids), and often carry an E, E+, and/or DYW domain at the C-terminus. As ubiquitous RNA binding factors, PPR proteins were found to take part in almost all of the post-transcriptional processing in mitochondria and plastids (Barkan and Small, 2014). The PLSclass PPR proteins mediate the C-to-U editing by specific binding to the nucleotide sequence upstream of the editing site through the PPR motifs (Tasaki et al., 2010;Okuda and Shikanai, 2012;Jiang et al., 2018). Furthermore, C-to-U editing also involves deamination of the cytidine (Blanc et al., 1995;Yu and Schuster, 1995), and the DYW domain possesses the cytidine deaminase (CDA) activity as demonstrated in PpPPR65 and PpPPR56 (Oldenkott et al., 2019;Hayes and Santibanez, 2020). The DYW domain also contains the motif of zinc-binding signature residues [HxE(x)nCxxC] commonly found in deaminases (Bhattacharya et al., 1994). Most DYW domains in higher plant PPR-DYW proteins contain the CDAs-like signature motif, which has been found to be indispensable for RNA editing (Boussardon et al., 2014;Hayes et al., 2015;Wagoner et al., 2015).
Accumulating evidence suggests that C-to-U editing is carried out by a protein complex, termed "editosome" Small et al., 2020). Recent studies have shown that proteins of distinct families are involved in RNA editing in addition to the PPRs. These include multiple organelle RNA editing factors (MORFs)/RNA-editing factor interacting proteins (RIPs), organelle RRM proteins (ORRMs), organelle zinc-finger 1 (OZ1), and other proteins as reviewed by Small et al. (2020). MORFs/RIPs harboring a conserved MORF/RIP motif have been shown to be required for the editing at a large number of sites in mitochondria and nearly all sites in plastids in Arabidopsis (Bentolila et al., 2013). MORFs also selectively interact with PPR proteins and form homo-or heterodimers (Bentolila et al., 2012;Takenaka et al., 2012;Glass et al., 2015;Zehrmann et al., 2015;Haag et al., 2017). ORRM proteins, on the other hand, are required for the editing at many sites in mitochondria (ORRM2, ORRM3, ORRM4, and ORRM5) (Shi et al., , 2016, and more than half of the sites in plastids in Arabidopsis (ORRM1 and ORRM6) (Sun et al., 2013;Hackett et al., 2017). ORRM proteins can interact with themselves, other ORRMs, and MORFs (Sun et al., 2013;Shi et al., 2015Shi et al., , 2016. ORRM1, specifically, is copurified with OZ1, which is required for the editing at 81% plastid target Cs . A recent report showed that an active editing complex contains PPRs, RIPs, ORRM1, OZ1, and ISE2 in maize chloroplasts (Sandoval et al., 2019).
Plant mitochondria harbor two mono-hemic c-type cytochromes: Cytochrome c (Cytc) and Cytochrome c1 (Cytc 1 ). They are the essential factors of the oxidative phosphorylation (OXPHOS) chain. Cytc 1 , anchored in the inner membrane, is a core subunit of complex III. On the other hand, Cytc, peripherally associated with the inner membrane, shuttles electrons from complex III to complex IV . After translation, apo-cytochrome c peptides are no longer functional. They have to undergo a maturation process in which the heme prosthetic groups are covalently attached to apo-Cytc 1 and apo-Cytc via thioether bonds. In gram-negative bacteria, 8-9-ccm genes (ccmA to ccmI) are involved in the cytochrome c maturation (CCM), referred to as the CCM pathway or system I (Thony- Meyer et al., 1995). Plant mitochondria inherit the CCM pathway from the prokaryote ancestor , and in maize specifically, the CCM pathway involves seven proteins: three (CCMA, CCME, and CCMH) encoded by the nuclear genes and four (CcmB, CcmC, CcmF N , and CcmF C ) by the mitochondrial genes (Clifton et al., 2004;Meyer et al., 2005;Rayapuram et al., 2008). Furthermore, loss-of-function of CCM factors often causes embryo lethality in plants. For example, a deficiency of CcmF N protein in the emp7 and ppr27 mutants results in a deficiency of Cytc and Cytc 1 , leading to impaired assembly of complex III and arrested seed development in maize (Sun F. et al., 2015;Liu et al., 2020).
In this study, we report the function of an uncharacterized mitochondrion-targeted PPR-DYW protein, erythropoietinmimetic peptide 17 (EMP17), in maize. The results demonstrate that EMP17 is required for the editing at four sites of ccmF C and nad2 transcripts in mitochondria. Deficiency in the editing at ccmF C -799 and nad2-677 causes amino acid changes in the encoded protein. We provide further evidence that the Ser-to-Pro change at CcmF C -267 impairs the CcmF C activity, blocks the maturation of Cytc and Cytc 1 , and disrupts the assembly of complex III, which attributes to the arrest of seed development in the emp17 maize mutant.

Plant Materials and Growth Conditions
The emp17 mutant was isolated from the UniformMu mutagenic population in nearly isogenic W22 background (McCarty et al., 2005). The plants were grown in the experimental field of Shandong University in Qingdao, Shandong province under natural conditions.

Construction of Emp17 Overexpression Plants
The protein-coding region of Emp17 complementary DNA (cDNA) was placed downstream of the maize ubiquitin 1 (Ubi-1) promoter in the pUNTF vector. This pUNTF-Emp17 construct was transformed into maize inbred KN5585 via callus transformation. Positive transgenic plants of Emp17 were screened by using the Bar gene, and further verified by PCR using the Ubi-F primer anchored to the ubiquitin sequence in the vector and the Emp17 specific primer Emp17-R1.

Subcellular Localization
The 420 bp, 5 -sequence of Emp17 encoding the 140 amino acids of the N-terminus region was cloned into pENTR/D-TOPO vector (Invitrogen in ThermoFisher Scientific, http:// www.thermofisher.com), and then transferred to the pBI221 vector to create the Emp17 N 140 -GFP fusion by the Gateway site-specific recombination. The resulting construct, pBI221-Emp17 N 140 -GFP, was introduced into Arabidopsis protoplasts. MitoTracker Red was used as a marker for mitochondria. The fluorescence signals were detected under a ZEISS LSM 880 confocal microscope.

Light Microscopy of Cytological Sections
Wild type and emp17 kernels were harvested at 10 and 14 days after pollination (DAP) from the selfed ears of emp17 (+/−) heterozygous plants. The kernels were fixed in 4% paraformaldehyde at 4 • C for 24 h. The fixed kernels were dehydrated, cleared, infiltrated, embedded, sectioned, stained, and observed as described by Liu et al. (2013).

RNA Extraction, RT-PCR, and qRT-PCR
Total RNA was extracted from developing kernels and other tissues of maize with the RNeasy Plant Mini Kit (Qiagen, http: //www.qiagen.com) according to the manufacturer's instruction. The potential gDNA contamination was removed through Dnase I (New England Biolabs, www.neb.sg) treatment. The single-stranded cDNA was obtained by reverse transcription reaction using the Transcriptor First Strand cDNA Synthesis kit (Thermo Fisher Scientific, http://www.thermofisher.com). Quantitative real-time PCR (qRT-PCR) with SYBR green (Bio-Rad, http://www.bio-rad.com) was performed in a Roche Light Cycler 96. The relative transcript level was calculated as described in previous studies . ZmActin (GRMZM2G126010) was used as a control in RT-PCR and qRT-PCR. Detailed information of the primers was listed in Supplementary Table 1.

Direct Sequencing of RT-PCR Amplicons
RNA sample preparation of the wild type and emp17 kernels and reverse transcription were performed as described above. The resulting cDNA was used as a template to amplify the 35 mitochondrial predicted-protein coding genes by RT-PCR. RNA editing was analyzed by directly sequencing the RT-PCR amplicons. The primers were listed in Supplementary Table 1. Three biological repetitions were analyzed.

Blue Native-PAGE and Mitochondrial Complexes Assembly and Activity Assay
Crude mitochondria were extracted from the wild type and emp17 kernels at 14 DAP as described previously (Li et al., 2014). Mitochondrial complexes were separated by the blue native polyacrylamide gel electrophoresis (BN-PAGE) assay. Ingel staining of complex I and supercomplex I + III 2 activity was performed as described previously (Meyer et al., 2009). The assembly of complex III and complex V was detected by transferring the protein complexes to nitrocellulose membranes and hybridized with anti-Cytc 1 and anti-Atp1 antibodies. Ingel staining of complex IV was detected as described by Wang et al. (2020). Crude mitochondria were extracted and protein concentration was determined by the Bradford assay kit (Bio-Rad, http://www.bio-rad.com). For Western detection of complex proteins, 8 µg protein from emp17, and 8 µg, 4 µg (1/2), and 2 µg (1/4) protein from wild type were subjected to sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). Proteins were transferred to nitrocellulose membranes and detected by anti-Nad9, anti-Cytc 1 , anti-Cytc, anti-Cox2, and anti-Atp1 as described previously (Sun F. et al., 2015).

Respiration Rate Assay
The respiration rates were determined according to the protocol described by Wang et al. (2015). O 2 consumption of the wild type and emp17 kernels was measured at dark by using a Chlorolab II liquid-phase oxygen electrode (Hansatech, http://www.hansatech-instruments.com/) in a reaction buffer containing 10 mM HEPES, 10 mM MES, and 2 mM CaCl 2 (where pH was adjusted to 6.8 with KOH). Total respiration rate (V t ) was defined by the O 2 consumption before adding any respiration inhibitor. The cytochrome pathway capacity (V cyt ) was determined by V t minus the O 2 consumption in the presence of 2 mM salicylhydroxamic acid (SHAM) and the alternative pathway capacity (V alt ) was defined by V t minus the O 2 consumption in the presence of 2 mM potassium cyanide (KCN). All three respiration rates were expressed as nmol O 2 min −1 g −1 fresh weight. The protein structure of EMP17. aa, amino acid; TP, targeted peptide; L, PPR-like long motif; S, PPR-like short motif. (B) Evolutionary relationships of putative EMP17 orthologs. The evolutionary history was inferred using the Neighbor-Joining method. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1,000 replicates) are shown next to the branches. (C) The subcellular localization of EMP17. The 140 amino acids fragment from the N-terminus of EMP17 was fused with green fluorescence protein (EMP17 N 140 :GFP) and transiently expressed in Arabidopsis protoplasts. Fluorescence signals were observed by confocal microscope ZEISS LSM 880. The mitochondria were stained by MitoTracker Red. DIC, differential interference contrast; Bar = 10 µm.

EMP17 Is a Mitochondrion-Targeted PPR-DYW Protein
PPR proteins form one of the largest protein families in maize with over 520 members annotated in the B73 genome.
EMP17 was predicted to localize in mitochondria by the TargetP 1 and Predotar algorithms. 2 To experimentally localize EMP17, the N-terminal region containing 140 aa of EMP17 was fused with the green fluorescent protein (GFP) and then transiently expressed in Arabidopsis protoplasts. GFP signals were detected in punctated dots that merged with the mitochondria stained by the MitoTracker Red ( Figure 1C). No GFP signals were found in other cellular compartments ( Figure 1C), indicating that EMP17 is specifically localized in mitochondria.

Loss of the EMP17 Function Severely Arrests Embryogenesis and Endosperm Development in Maize
To study the function of EMP17, a Mutator (Mu) insertion mutant (emp17) was isolated from the UniformMu mutagenic population (McCarty et al., 2005). A Mu3 element was confirmed to be inserted at +660 bp from the translation start codon of Emp17 (Figure 2A). The selfed progenies of emp17 (+/−) heterozygous plants produced about 1/4 kernels with an empty pericarp phenotype (Figure 2B), and no wild type Emp17 transcript was detected in the embryo and endosperm of the empty pericarp kernels ( Figure 2C). Linkage analysis in an F2 population showed that only plants segregating emp mutant phenotype carried Mu insertion in Emp17, indicating that the Mu insertion in Emp17 is either tightly linked to or the cause of the emp phenotype ( Figure 2B and Supplementary Figure 3).
To confirm whether Emp17 is the causal gene for the empty pericarp phenotype, we created transgenic plants over-expressing Emp17 (Emp17-OE) in the inbred line KN5585 by placing Emp17 under the Ubi-1 promoter. Eight independent transgenic lines (Emp17-OE1 to OE8) were obtained, and the expression level of Emp17 in these lines was 5-9.5 times higher than that in the wild type as detected by qRT-PCR ( Figure 2D). We crossed Emp17-OE1 with the emp17 (+/−) heterozygous plants and selfed the F1 to obtain the F2 progeny. The F2 seedlings were genotyped by PCR. To distinguish between the endogenous Emp17 and transgene Emp17, the Emp17-F1 primer was anchored to the 5'-UTR of the endogenous Emp17, ensuring that it could not anneal to the transgene that lacks this 5 -UTR sequence. The ubiquitin specific primer Ubi-F and Mu primer TIR8 were then used. Thus, the Emp17-F1/R1 primer 1 http://www.cbs.dtu.dk/services/TargetP 2 http://urgi.versailles.inra.fr/predotar/predotar.html pair amplified the endogenous Emp17, the Ubi-F/Emp17-R1 pair amplified the transgene, and the TIR8/Emp17-F1 pair detected the Mu3 insertion in Emp17. Genotyping 32 of the F2 seedlings identified 5 seedlings that were homozygous for emp17 harboring the Emp17 transgene (Figures 2E,F). The Emp17-F1/R1 primer pair amplified a larger fragment in these five emp17 (−/−) seedlings, which were proven to contain the Mu3 element. These five seedlings showed normal growth and development compared with the wild type (Figure 2E), indicating over-expression of Emp17 rescued the embryo-lethal phenotype of this mutant. These results demonstrate that Emp17 (GRMZM2G019689) is the causal gene for the empty pericarp phenotype of maize kernels in emp17.
The emp17 kernels were substantially smaller than the wild type siblings throughout kernel development. At 14 DAP, the emp17 embryo and endosperm was smaller than the wild type (Figures 3A,B). Paraffin sectioning indicated that the embryogenesis and endosperm development were severely arrested in emp17. At 10 DAP, leaf primordia (LP), shoot apical meristem (SAM), and root apical meristem (RAM) were clearly developed in the wild type embryo ( Figure 3C). Conversely, the mutant embryo only proceeded to the early transition stage (Figure 3F). At 14 DAP, the wild-type embryo entered the late embryogenesis stage (Figures 3D,E), while the mutant embryo remained at the transition stage and the endosperm was arrested at the cellularization stage (Figures 3G,H).
Loss-of-Function of EMP17 Abolishes the Editing at ccmF C -799 and nad2-677 Sites As previously reported, most of the known PPR-DYW proteins function in the C-to-U editing of organellar RNA (Guillaumot et al., 2017;Wang et al., 2019). The maize mitochondrial genome was predicted to encode 22 electron transport chain proteins, 11 ribosomal proteins, a maturase MatR, and a membrane transporter protein MttB (Clifton et al., 2004). These 35 gene transcripts were amplified from the emp17 mutant and wild type kernels, and directly sequenced. Comparison of the sequences revealed that the C-to-U editing at ccmF C -799 and nad2-677 sites was abolished in emp17, and completely edited in the wild type ( Figure 4A). In addition, the editing at ccmF C -906 and -966 sites was substantially decreased in emp17 in comparison with the wild type ( Figure 4A). In the Emp17-OE1:emp17 (−/−) seedlings, the editing at nad2-677, ccmF C -799, -906, and -966 sites was restored ( Figure 4A). These results indicate that the loss-of-function of Emp17 abolishes the editing at ccmF C -799 and nad2-677, and decreases the editing at ccmF C -906 and -966.
PPR proteins recognize the target RNA sequence in a "one-PPR motif: one-nucleotide" manner, in which the sixth amino acid residue in the first PPR-motif and the 1 amino acid residue in the next PPR-motif specifies the base of the RNA sequence (Barkan et al., 2012). Based on this code, the PPR motifs of EMP17 were aligned with the upstream sequence of nad2-677, ccmF C -799, -906, and -966 sites. The results showed that the codes aligned mostly with the ccmF C -799 and nad2-677 sites, but poorly with the ccmF C -906 and -966 sites (Figure 4B), suggesting that EMP17 probably binds strongly to the upstream sequences of the ccmF C -799 and nad2-677 sites, but weakly with the sequences of the ccmF C -906 and -966 sites.

The Ser 267 Residue in CcmF C and Phe 226 Residue in Nad2 Are Conserved in Plants
The deficient editing at ccmF C -799 and nad2-677 in the emp17 mutant led to a Ser-to-Pro change at CcmF C -267 and a Phe-to-Ser change at Nad2-226, respectively ( Figure 4A). Alignment of both gDNA and cDNA sequences of the ccmF C and nad2 orthologs indicated that these two amino acids were conserved in lower plants (Physcomitrella patens and Marchantia polymorpha), dicots (Glycine max, Nicotiana tabacum, Beta vulgaris, and Brassica napus), and monocots (Zea mays, Triticum aestivum, and Oryza sativa) (Figures 5A,B). The conservation of these two amino acid residues implies that these residues are probably important to the functional integrity of the CcmF C and Nad2 proteins.

Mitochondrial Complex III Biogenesis Is Severely Reduced in emp17
Nad2 is a subunit located in the hydrophobic arm of mitochondrial complex I (Braun et al., 2014), and CcmF C is involved in the maturation of Cytc and Cytc 1 , which are components of mitochondrial complex III . To assess the impact of the amino acid residue alteration on these complexes, we analyzed the assembly and activity of mitochondrial complexes and the complex proteins in emp17. Mitochondria were isolated from the embryo and endosperm of the emp17 mutant and wild type, respectively. The mitochondrial complexes were separated with BN-PAGE. Coomassie Brilliant Blue (CBB) staining showed that the level of complex I was comparable between emp17 and wild type, but the levels of complex III and supercomplex I + III 2 were remarkably decreased in emp17 ( Figure 6A). Furthermore, in-gel staining of the NADH dehydrogenase activity showed consistent results ( Figure 6B). Western blot analysis indicated that the level of complex III as detected by anti-Cytc 1 antibody was drastically decreased in emp17 (Figure 6C). Similarly, western blotting detection of complex V using anti-Atp1 antibody and in-gel staining of complex IV activity indicated that complex V and complex IV were increased in emp17 as compared with wild type (Figures 6D,E). Western blot assays using anti-Nad9, Cytc 1 , Cytc, Cox2, and Atp1 antibodies showed that the levels of Cytc and Cytc 1 was dramatically reduced. On the other hand, the level of Nad9 and Cox2 was substantially increased and the level of Atp1 was moderately increased in emp17 as compared with wild type (Figure 6F). These results indicate that the abolished editing at nad2-677 site causing the Phe-to-Ser change at Nad2-226 appears not to significantly affect the assembly of mitochondrial complex I. However, the Ser-to-Pro change at CcmF C -267 severely inhibits the maturation of Cytc 1 and Cytc and biogenesis of complex III, implying a critical role of Ser 267 to the function of CcmF C .

Emp17 Loss of Function Decreases Cytochrome Respiration Rate and Increases ZmAOX Expression
As previously reported, blocking of the cytochrome pathway often induces the alternative non-phosphorylating pathway in the respiratory chain (Yang et al., 2017;Wang et al., 2019;Xu et al., 2020). The maize genome hosts three ZmAOX genes (ZmAOX1, ZmAOX2, and ZmAOX3) (Karpova et al., 2002). As indicated by RT-PCR and qRT-PCR results, the level of ZmAOX2 and ZmAOX3 transcripts was dramatically increased in emp17, especially ZmAOX2 (Figures 7A,B). Western blotting analysis revealed that ZmAOX is expressed at low levels in the wild type kernels, but drastically enhanced in emp17 (Figure 6F).
To assess the physiological impact, we then measured the V t , V alt , and V cyt by using a Chlorolab II liquid-phase oxygen electrode, and specific inhibitors SHAM and KCN, respectively . The ratios of V alt to V t and V cyt to V t in emp17 and wild type kernels were calculated. The results indicated that V t and V cyt in emp17 were decreased to about 18 and 10% in comparison with that in wild type, respectively ( Table 1). The ratio of V alt /V t (75%) in emp17 is higher than that by about twofold (43%) in wild type. These data confirm that mitochondrial complex III dysfunction severely reduces the cytochrome respiration and induces the alternative respiration pathway in emp17. Co-evolution Between EMP17 and Editing Sites ccmF C -799 and nad2-677 Putative orthologous proteins of EMP17 can be found in many sequenced plant species in the NCBI and Uniport database ( Figure 1B). Phylogenetic analysis revealed a high degree of conservation of putative EMP17 orthologs in A. trichopoda and mono-and dicotyledonous species. EMP17 shared an over 50% sequence identity with most of its orthologous proteins (Figure 1B), i.e., 54.7% with that in A. trichopoda, the single living species of the sister lineage to all other extant flowering plants (Soltis et al., 2011). However, a putative EMP17 ortholog cannot be found in O. sativa, A. thaliana, B., and G. hirsutum. The most closely related protein in rice is Os12g0109300, sharing a 40.99% sequence identity with EMP17. In turn, Os12g0109300 shares an 81.94% sequence identity with maize protein GRMZM5G811022. As such, Os12g0109300 is unlikely to be an EMP17 ortholog in rice. The most closely related homologs of EMP17 are OTP82 in A. thaliana, hypothetical protein XP_013641116 in B. napus (accession number XP_013641116), and hypothetical protein (accession number KAG4215874) in G. hirsutum. All of these proteins share a less than 38% sequence identity with EMP17, much lower than that with the A. trichopoda homolog (54.7%).
The results suggest that these species may have lost the Emp17 gene in the genome.
This raises the question of why this protein is conserved in some species, but lost in others. The gDNA sequence of CcmF C shows that the ccmF C -799 site is "T" in A. trichopoda and eudicots. However, in monocots, both "T" and "C" present (Supplementary Figure 5A). In all the species harboring ccmF C -799C in mitochondrial DNA, putative orthologs of EMP17 can be identified in the nuclear genomes (Phoenix dactylifera, Z. mays, Sorghum bicolor, and T. aestivum) (Supplementary Figure 5A). In the species with ccmF C -799T in the mitochondrial gene, putative ortholog of EMP17 may be lost (Supplementary Figure 5A). A consistent result can be found between the existence of EMP17 orthologs in the nuclear genome and the nad2-677 site in mitochondria as well (Supplementary Figure 5B). These results suggest that EMP17 orthologs probably exist in the early flowering plants that do not require the editing function of EMP17 at the ccmF C -799 or nad2-677 sites because both sites are "T." But later in evolution, when the "T" was mutated to "C, " EMP17 was recruited for the editing function. For the species that maintain a "T" at this site, the Emp17 orthologs may be degenerated or lost. This notion offers a possible explanation for the disappearance of a clear EMP17   This study provides strong evidence that EMP17, a previously uncharacterized PPR-DYW protein, functions in the editing of mitochondrial transcripts in maize. Since only 6 out of the 82 maize PPR-DYW proteins have been fully characterized so far, the elucidation of the EMP17 function adds a new piece of information to the repertoire of this large protein family. Our results show that EMP17 is exclusively localized in mitochondria (Figure 1C), and loss-of-function in EMP17 abolishes the editing at ccmF C -799 and nad2-677 sites while reducing the editing at ccmF C -906 and ccmF C -966 sites ( Figure 4A). Conversely, the expression of Emp17 restores the editing defects in the emp17 mutant ( Figure 4A), demonstrating that EMP17 is required for the editing at these sites. In addition, based on the "one PPR motif: one nucleotide" recognition codes (Barkan et al., 2012), the 6,1'-amino acid residue combinations of EMP17 align well with the upstream sequences of ccmF C -799 and nad2-677 sites, and weakly with those of the ccmF C -906 and -966 sites (Figure 4B), suggesting that EMP17 may recognize its substrates specifically. Furthermore, the lack of editing is accompanied by a reduced mitochondrial complex III assembly, inhibition of the cytochrome pathway, elevated alternative pathway, and severely reduced respiration rates in the emp17 mutants. All of these results provide convincing evidence that EMP17 functions on the C-to-U editing at these sites in mitochondria and the loss of function of EMP17 impairs the cytochrome respiratory pathway.
The severely reduced assembly of complex III suggests that the impaired OXPHOS chain can likely be owed to the loss of editing at the ccmF C sites in emp17 ( Figure 4A). Except for the three ccmF C sites and the nad2-677 site, no other defects were found in the transcripts that are directly associated with complex III. In the four mitochondrion-encoded CCM pathway proteins (CcmB, CcmC, CcmF N , and CcmF C ) that are essential for the CCM and biogenesis of complex III , no defects were detected in the transcripts of ccmB, ccmC, and ccmF N . Additionally, expression levels of ccmB, ccmF C , and ccmF N in emp17 were either indistinguishable from that in wild type, or increased (ccmC) (Supplementary Figure 4). CcmF C has been shown to be important for plant growth and development, as the loss of ccmF C expression led to a deficiency of the c-type cytochromes and complex III in wtf9 (Francs- Small et al., 2012). Consistent with that, the Cytc 1 and Cytc proteins were barely detectable (Figure 6F), and the assembly of complex III was severely decreased in emp17 ( Figure 6C). These results suggest that the unedited ccmF C -799 disrupts the complex III assembly and leads to a dysfunction of CCM.
The editing deficiency at the ccmF C -799 site causing the Ser 267 -to-Pro 267 change in CcmF C is probably the major cause Respiration rate (nmol O 2 min −1 g −1 fresh weight) Mitochondrial V t , V cyt , and V alt were measured by the Chlorolab II liquid-phase oxygen electrode. The activity of cytochrome c oxidase and an alternative oxidase were inhibited by 2 mM KCN and 2 mM SHAM, respectively. Data are mean values SEs from three independent biological samples. Statistical differences (P < 0.01*) in the same row were based on student's t-test analysis. V t , total respiration rate; V cyt , cytochrome pathway capacity; V alt , alternative pathway capacity; KCN, potassium cyanide; SHAM, salicylhydroxamic acid; SEs, standard errors. for the inhibited kernel development in emp17. In plant mitochondria, Cytc and Cytc 1 are the essential components of mitochondrial complex III in the OXPHOS chain. Maturation of Cytc and Cytc 1 is crucial to mitochondrial functions, and hence, to plant growth and development. Impairment of c-type cytochrome maturation arrests seed development or plant growth. A deficiency of mature Cytc 1 and Cytc in the ccmh mutant causes embryo lethality in Arabidopsis (Meyer et al., 2005). Similarly, a deficiency of Cytc 1 and Cytc resulting from abolished editing at the ccmF N -1553 and -1357 sites in emp7 and ppr27, respectively, results in embryo lethality in maize (Sun F. et al., 2015;Liu et al., 2020). The lack of Cytc 1 and Cytc blocks the assembly of mitochondrial complex III, decreases the cytochrome respiration rate, and elevates the alternative nonphosphorylating pathway in the emp7 and ppr27 mutants. For the emp17 mutants, we found consistent results. The assembly of complex III was severely inhibited in emp17, and V t and V cyt were decreased to about 18 and 10% in comparison with that in wild type, respectively (Table 1), and the expression of ZmAOX2 and ZmAOX3 was dramatically increased in emp17 compared with wild type (Figures 6F, 7A,B). Thus, the lack of Cytc 1 and Cytc and severely reduced biogenesis of mitochondrial complex III blocks the cytochrome pathway and impairs the kernel development.
Our data suggest that the abolished editing at the nad2-677 site is probably not a major cause for the defective seed development in emp17. The editing deficiency at nad2-677 results in a Phe-to-Ser change at Nad2-226. As an essential component of the mitochondrial complex I, a deficiency of Nad2 impairs the assembly of complex I and arrests kernel development in maize (Xiu et al., 2016;Yang et al., 2020). However, the Phe-to-Ser change at Nad2-226 in emp17 does not significantly affect the assembly of mitochondrial complex I, or its activity as determined by in-gel NADH dehydrogenase activity ( Figure 6B). Although we cannot rule out the possibility that the detected assembled complex I is in fact non-functional in the electron transfer chain, it is likely that the amino acid residue change in Nad2-226 may have a less detrimental impact on the Nad2 function.
The Ser 267 Residue in CcmF C Is Essential for the CcmF C Function The lack of editing at ccmF C -799 as a result of the Emp17 mutation constitutes a surrogate mutation of Ser 267 to Pro 267 in the CcmF C protein. The severe impact of this mutation on CCM and complex III assembly illustrates the importance of the Ser residue at CcmFc-267 for the CcmF C function. In bacteria and plant mitochondria, CcmF and CcmH are proposed to take part in the final step of CCM, ligating heme delivered by CcmE to apo-cytc Rayapuram et al., 2008). In plant mitochondria, the ccmF gene has been split into multiple genes. For instance, ccmF is split into three genes (ccmF N 1 , ccmF N 2 , and ccmF C ) in Arabidopsis (Unseld et al., 1997), and two genes (ccmF N and ccmF C ) in maize (Clifton et al., 2004). Structure prediction and trypsin digestion experiments suggest that AtCcmF C has six transmembrane helices, three intermembrane space loops, and four mitochondrial matrix domains (Rayapuram et al., 2008). Transmembrane helice prediction by the TMHMM Server v. 2.0 3 indicates that the structure of ZmCcmF C and AtCcmF C is quite similar. The two proteins share a 78% sequence identity. Based on this structure, Ser 267 is located in the third intermembrane space loop named domain 6 (D6) of CcmF C (Figure 8A). Alignment of the protein sequences shows that the Ser 267 in various CcmF C proteins is highly conserved across species (Figure 5A), suggesting its importance to CcmF C function. It is known that Pro residue is a disruptor of protein α-helix and not favored in the β-sheet structures, as it has structure is limited and cannot complete the H-bonding network. For example, the Leu-to-Pro change in Nad7-279 and Atp6-213 and in the α-helix region is attributed to the destruction of the Nad7 and Atp6 function in the maize mutant emp18 and smk1 (Li et al., 2014(Li et al., , 2019. The Ser 267 in CcmF C was predicted in a β-sheet by the Swiss-model algorithm 4 (Figure 8B). Proline is not favored in β-sheet structures as it cannot complete the H-bonding network. It is possible that the Ser 267 -to-Pro 267 mutation in CcmF C in the β-sheet may negatively impact the structural stability of D6 that leads to nonfunctional CcmF C . The exact function of the D6 intermembrane space loop in CcmF C is unknown, however, these results imply that the Ser 267 residue in D6 plays a critical role for the function of CcmF C .

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
YW and B-CT designed the research and analyzed the data. YW and X-YL conducted most of the experiments. FS performed the BN gel assay. Z-QH, Y-YL, and Z-QG participated in the linkage and genetic complementarity analysis. XW performed the respiration rate assay. YW, X-YL, Y-ZY, AS, and B-CT wrote the article. All authors contributed to the article and approved the submitted version.