Comparative Proteomic Analysis Reveals Elevated Capacity for Photosynthesis in Polyphenol Oxidase Expression-Silenced Clematis terniflora DC. Leaves

Polyphenol oxidase (PPO) catalyzes the o-hydroxylation of monophenols and oxidation of o-diphenols to quinones. Although the effects of PPO on plant physiology were recently proposed, little has been done to explore the inherent molecular mechanisms. To explore the in vivo physiological functions of PPO, a model with decreased PPO expression and enzymatic activity was constructed on Clematis terniflora DC. using virus-induced gene silencing (VIGS) technology. Proteomics was performed to identify the differentially expressed proteins (DEPs) in the model (VC) and empty vector-carrying plants (VV) untreated or exposed to high levels of UV-B and dark (HUV-B+D). Following integration, it was concluded that the DEPs mainly functioned in photosynthesis, glycolysis, and redox in the PPO silence plants. Mapman analysis showed that the DEPs were mainly involved in light reaction and Calvin cycle in photosynthesis. Further analysis illustrated that the expression level of adenosine triphosphate (ATP) synthase, the content of chlorophyll, and the photosynthesis rate were increased in VC plants compared to VV plants pre- and post HUV-B+D. These results indicate that the silence of PPO elevated the plant photosynthesis by activating the glycolysis process, regulating Calvin cycle and providing ATP for energy metabolism. This study provides a prospective approach for increasing crop yield in agricultural production.


Introduction
Polyphenol oxidase (PPO) is an oxidoreductase that catalyzes the oxidation of monophenols and/or o-diphenols to o-quinones, which form brown melanin pigments in fruits and vegetables by covalently modifying and cross-linking proteins [1,2]. However, PPO is not only involved in the formation of pigments, but also plays a crucial role in the biosynthesis of secondary metabolites such as aurones [3] and betalins [4]. PPO solitarily catalyzes the hydroxylation and oxidative cyclization of chalcones, leading to the formation of aurone [3]. In addition, PPO, a tyrosinase can also hydroxylate tyramine into dopamine, which in the presence of betalamic acid yields dopamine-betaxanthin that of genome-walking technologies ( Figure S2a). Sequence analysis of CtPPO cDNA revealed a 166-bp 5 -untranslated region (UTR) located upstream of a start codon and an open reading frame (ORF) of 1681 bp encoding 586 amino acids with a calculated molecular mass of 66.11 kDa and an isoelectric point of 7.18. Blast analysis of the predicted CtPPO amino acid sequence revealed 62% similarity with Nelumbo nucifera. To further understand the evolutionary relationships between CtPPO and PPO in other plants, a phylogenetic tree was constructed based on amino acid sequence. Phylogenetic analysis showed CtPPO has a close evolutionary relationship with the PPO gene in N. nucifera, but not other plant species ( Figure S2b).

Virus-Induced Gene Silencing of CtPPO in C. terniflora DC
VIGS provides an alternative approach for functional analysis of genes. Here, to determine the role of CtPPO, we silenced CtPPO using tobacco rattle virus (TRV)-mediated gene silencing. We constructed two modified TRV vectors: TRV-CtPPO and TRV-PDS. To determine the extent of silencing, we performed qRT-PCR to analysis the transcript levels of the CtPPO gene in inoculated plants. RNA extracted from C. terniflora agro-infiltrated with empty vector TRV or TRV-CtPPO or control leaves 20 days post-infection were used for qRT-PCR. As shown in Figure 1a, CtPPO expression levels were reduced in plants infiltrated with TRV-CtPPO compared to the control plants (VIGS-vector and control). To further confirm CtPPO expression, we measured CtPPO activity in control, TRV-vector and TRV-CtPPO plants. In VIGS-CtPPO (VC), PPO activity was significantly lower than in VIGS-vector (VV) and control, illustrating that a successful VC plant model was constructed (Figure 1b). Despite these obvious changes in CtPPO in plants, however, there were no macroscopic phenotypic differences in the VC, VV, and control leaves at the starting point ( Figure 1c). application of genome-walking technologies ( Figure S2a). Sequence analysis of CtPPO cDNA revealed a 166-bp 5'-untranslated region (UTR) located upstream of a start codon and an open reading frame (ORF) of 1681 bp encoding 586 amino acids with a calculated molecular mass of 66.11 kDa and an isoelectric point of 7.18. Blast analysis of the predicted CtPPO amino acid sequence revealed 62% similarity with Nelumbo nucifera. To further understand the evolutionary relationships between CtPPO and PPO in other plants, a phylogenetic tree was constructed based on amino acid sequence. Phylogenetic analysis showed CtPPO has a close evolutionary relationship with the PPO gene in N. nucifera, but not other plant species ( Figure S2b).

Virus-Induced Gene Silencing of CtPPO in C. terniflora DC.
VIGS provides an alternative approach for functional analysis of genes. Here, to determine the role of CtPPO, we silenced CtPPO using tobacco rattle virus (TRV)-mediated gene silencing. We constructed two modified TRV vectors: TRV-CtPPO and TRV-PDS. To determine the extent of silencing, we performed qRT-PCR to analysis the transcript levels of the CtPPO gene in inoculated plants. RNA extracted from C. terniflora agro-infiltrated with empty vector TRV or TRV-CtPPO or control leaves 20 days post-infection were used for qRT-PCR. As shown in Figure 1a, CtPPO expression levels were reduced in plants infiltrated with TRV-CtPPO compared to the control plants (VIGS-vector and control). To further confirm CtPPO expression, we measured CtPPO activity in control, TRV-vector and TRV-CtPPO plants. In VIGS-CtPPO (VC), PPO activity was significantly lower than in VIGS-vector (VV) and control, illustrating that a successful VC plant model was constructed (Figure 1b). Despite these obvious changes in CtPPO in plants, however, there were no macroscopic phenotypic differences in the VC, VV, and control leaves at the starting point ( Figure  1c). Data are shown as mean ± SD of the independent biological replicates. Asterisks indicate significant changes as measured by Student's t-test (* p < 0.05, ** p < 0.01, and *** p < 0.001); (b) PPO activity was assayed in C. terniflora leaves of control, VIGS-vector and VIGS-CtPPO plants. Data are shown as mean ± SD of independent biological replicates. Asterisks indicate significant changes as measured by Student's t-test (* p < 0.05, ** p < 0.01, and *** p < 0.001); (c) Morphology of control, VIGS-vector, and VIGS-CtPPO C. terniflora plants before and after HUV-B+D treatment. Photograph showing C. terniflora growth under different treatments. The damaged areas were marked with red arrows.

Effects of VIGS-CtPPO and VIGS-Vector on Leaf Proteins in C. terniflora DC.
To investigate the effects of CtPPO on C. terniflora leaf proteins, a gel-free/label-free proteomic technique was used. Proteins extracted from VC-and VV-infected C. terniflora leaves were reduced, alkylated, digested, and analyzed by nano-LC-MS/MS. The obtained proteomic data were analyzed using the UniProtKB/Swiss-Prot database and protein content was estimated by mol%. The functions of identified proteins were predicted based on comparisons with functional annotations of the Arabidopsis genome and classified using MapMan bin codes. The results of the functional analyses demonstrated that the proteins with increased expression were mainly enriched in photosynthesis, glycolysis, redox, protein metabolism, and secondary metabolism, while the proteins with decreased expression were mainly enriched in protein metabolism and photosynthesis. Notably, expression of proteins related to secondary metabolism, development, and glycolysis were increased in the VC plants, whereas the expression of proteins related to N-metabolism, amino acid metabolism, and cell were decreased (Figure 2a; Tables S8 and S9). and VIGS-CtPPO plants. Data are shown as mean ± SD of the independent biological replicates. Asterisks indicate significant changes as measured by Student's t-test (* p < 0.05, **p < 0.01, and ***p < 0.001); (b) PPO activity was assayed in C. terniflora leaves of control, VIGS-vector and VIGS-CtPPO plants. Data are shown as mean ± SD of independent biological replicates. Asterisks indicate significant changes as measured by Student's t-test (* p < 0.05, ** p < 0.01, and *** p < 0.001); (c) Morphology of control, VIGS-vector, and VIGS-CtPPO C. terniflora plants before and after HUV-B+D treatment. Photograph showing C. terniflora growth under different treatments. The damaged areas were marked with red arrows.

Effects of VIGS-CtPPO and VIGS-Vector on Leaf Proteins in C. terniflora DC.
To investigate the effects of CtPPO on C. terniflora leaf proteins, a gel-free/label-free proteomic technique was used. Proteins extracted from VC-and VV-infected C. terniflora leaves were reduced, alkylated, digested, and analyzed by nano-LC-MS/MS. The obtained proteomic data were analyzed using the UniProtKB/Swiss-Prot database and protein content was estimated by mol%. The functions of identified proteins were predicted based on comparisons with functional annotations of the Arabidopsis genome and classified using MapMan bin codes. The results of the functional analyses demonstrated that the proteins with increased expression were mainly enriched in photosynthesis, glycolysis, redox, protein metabolism, and secondary metabolism, while the proteins with decreased expression were mainly enriched in protein metabolism and photosynthesis. Notably, expression of proteins related to secondary metabolism, development, and glycolysis were increased in the VC plants, whereas the expression of proteins related to N-metabolism, amino acid metabolism, and cell were decreased (Figure 2a; Tables S8 and S9).

Effects of High-Level UV-B and Dark Treatment on Proteins in VIGS-CtPPO and VIGS-Vector C. terniflora DC. Leaves
HUV-B+D dramatically impacts plant physiology [24]. To further investigate the effects of CtPPO on C. terniflora, we exposed VC, VV, and control plants to HUV-B+D. Subsequent phenotypic analysis revealed the leaves to be wilted with burned patches and a high degree of crispation in control plants. The leaves of VV plants were similar to the control plants, while the leaves from the VC plants appeared significantly less damaged than in the control and VV plants (Figure 1c). To determine the expression profiles of leaf proteins in VC and VV after HUV-B+D, a gel-free/label-free proteomic approach was used. Proteins extracted from VC and VV plant leaves after HUV-B+D were reduced, alkylated, digested, and analyzed by nano-LC-MS/MS. Obtained proteomic data were analyzed using the UniProtKB/Swiss-Prot database and the protein content was estimated by mol%. The functions of identified proteins were predicted based on comparisons with functional annotations of the Arabidopsis genome and classified using MapMan bin codes. The results of the functional analyses demonstrated that HUV-B+D increased expression of proteins involved in most of the functional classes, including photosynthesis, stress, glycolysis, protein metabolism, and redox, in VC plants compared to VV plants. However, expression of proteins involved in N-metabolism was decreased in VC plants compared to VV plants (Figure 2b). For proteins differentially expressed in response to HUV-B+D in VV and VC plants, further integration analysis was performed. Figure 3 shows that the number of proteins involved in photosynthesis displaying increased expression was decreased in VC plants compared with that in VV plants. The number of increased proteins related to stress, OPP, secondary metabolism, and transport was increased in VC plants compared to VV plants. The number of decreased proteins associated with stress, glycolysis, amino acid metabolism, and carbon-one(C1)-metabolism was decreased in VC plants compared to VV plants. Furthermore, proteins involved in glycolysis, amino acid metabolism, cell, mitochondrial ETC, TCA, and signaling displaying increased expression were only detected in VC plants ( Figure 3, Tables S6 and S7).

MapMan Analysis of the VIGS-CtPPO and VIGS-Vector C. terniflora Leaf Proteomic Data at the Starting Point and after HUV-B+D
To examine changes in the levels of the identified proteins in-depth in VV and VC plants, the DEPs were analyzed using MapMan software. The analysis identified the main functional categories of the proteins displaying significant changes in expression as involved in the light reaction, photorespiration, and Calvin cycle. HUV-B+D conditions induced dramatic decrease in Calvin cycle and light reaction related proteins in VV, however, very interestingly, it was significantly rescued by the silencing of CtPPO in VC. Further analysis showed that the silencing of CtPPO led to the increase

MapMan Analysis of the VIGS-CtPPO and VIGS-Vector C. terniflora Leaf Proteomic Data at the Starting Point and after HUV-B+D
To examine changes in the levels of the identified proteins in-depth in VV and VC plants, the DEPs were analyzed using MapMan software. The analysis identified the main functional categories of the proteins displaying significant changes in expression as involved in the light reaction, photorespiration, and Calvin cycle. HUV-B+D conditions induced dramatic decrease in Calvin cycle and light reaction related proteins in VV, however, very interestingly, it was significantly rescued by the silencing of CtPPO in VC. Further analysis showed that the silencing of CtPPO led to the increase of Calvin cycle related proteins in VC than VV at starting point (Figure 4a,b). After exploring HUV-B+D, the silencing of CtPPO gave a further increase of Calvin cycle related proteins. Additionally, it also resulted in the synchronously increase on light reaction and photorespiration related proteins in VC than VV at HUV-B+D (Figure 4c,d).  Following this preliminary analysis, we further studied the differences in enzymes in the light reaction and Calvin cycle. The proteins involved in the light reaction and Calvin cycle was mostly decreased in VV but there was no obvious change in VC after HUV-B+D treatment. At starting point, 1 and 3 proteins in photosystem I and II of light reaction was increased in VC compared to VV. In Calvin cycle, 3 and 2 proteins involved in carboxylation and reduction was increased in VC compared with VV. However, after HUV-B+D treatment, 19 proteins related to NADP + -dependent aldehyde Following this preliminary analysis, we further studied the differences in enzymes in the light reaction and Calvin cycle. The proteins involved in the light reaction and Calvin cycle was mostly decreased in VV but there was no obvious change in VC after HUV-B+D treatment. At starting point, 1 and 3 proteins in photosystem I and II of light reaction was increased in VC compared to VV. In Calvin cycle, 3 and 2 proteins involved in carboxylation and reduction was increased in VC compared with VV. However, after HUV-B+D treatment, 19 proteins related to NADP + -dependent aldehyde reductase, ATP synthase, and glycerol phosphate dehydrogenase in light reactions was increased and 11 proteins related to the enzymes in reduction and regeneration of Calvin cycle was increased in comparison of VC and VV ( Figure S4).

Expression Profile Analysis of Glycolysis-Related Proteins in C. terniflora DC. Leaves
Glycolysis related proteins were significantly changed by the silencing of CtPPO. To further characterize the relationship between glycolysis and CtPPO, expression profile analysis of glycolysis-related proteins in C. terniclora leaves was conducted. The abundances of glyceraldehyde-3-phosphate dehydrogenase (GAPDH), UTP-glucose-1-phosphate uridylytransferase (GPUT), phosphoglycerate kinase (PGK), phosphoenolpyruvate carboxylase (PEPC), and enolase were increased in VC plants compared to VV plants at the starting point. Moreover, after HUV-B+D, the abundances of GPUT, fructose-1,6-bisphosphate aldolase (FBA), GAPDH, PGK, BPGM, enolase, and PEPC were all increased in VC plants compared to VV plants ( Figure 5).

Effects of VIGS-CtPPO on Photosynthesis Characteristics in C. terniflora DC. Leaves
To confirm the relationship between CtPPO and photosynthesis, the chlorophyll content was measured first. The results showed that the chlorophyll a, chlorophyll b, and total chlorophyll contents were all higher in VC plants than in VV and control plants ( Figure 6). Then the photosynthesis characteristics were measured by using a Li-6400 Portable Photosynthesis System.

Effects of VIGS-CtPPO on Photosynthesis Characteristics in C. terniflora DC. Leaves
To confirm the relationship between CtPPO and photosynthesis, the chlorophyll content was measured first. The results showed that the chlorophyll a, chlorophyll b, and total chlorophyll contents were all higher in VC plants than in VV and control plants ( Figure 6). Then the photosynthesis characteristics were measured by using a Li-6400 Portable Photosynthesis System. The photosynthesis rate, intercellular CO2, stomatal conductance, and transpiration rate were further measured in the study. The photosynthesis rate, transpiration rate, and stomatal conductance in VC plants were higher than in VV plants, while the intercellular CO2 was lower than that in VV ( Figure 7). After HUV-B+D, the damage of photosynthesis rate in VC plants was significantly lower than in VV and control plants, while the stomatal conductance and transpiration rate were low in control, VV, and VC plants. The intercellular CO2 was still lower in VC and control plants than in VV plants, but both were higher than that at the starting point ( Figure S5). The photosynthesis rate was dramatically decreased in by HUV-B+D in C. terniflora, however, it seems that the silencing of CtPPO has a mitigating effect to the damage degree in VC (Figure 7 and Figure S5). Data are shown as mean ± SD for independent biological replicates. Asterisks indicate significant changes as measured by Student's t-test (* p < 0.05, ** p < 0.01). fw: fresh weight.
The photosynthesis rate, intercellular CO 2 , stomatal conductance, and transpiration rate were further measured in the study. The photosynthesis rate, transpiration rate, and stomatal conductance in VC plants were higher than in VV plants, while the intercellular CO 2 was lower than that in VV ( Figure 7). After HUV-B+D, the damage of photosynthesis rate in VC plants was significantly lower than in VV and control plants, while the stomatal conductance and transpiration rate were low in control, VV, and VC plants. The intercellular CO 2 was still lower in VC and control plants than in VV plants, but both were higher than that at the starting point ( Figure S5). The photosynthesis rate was dramatically decreased in by HUV-B+D in C. terniflora, however, it seems that the silencing of CtPPO has a mitigating effect to the damage degree in VC (Figure 7 and Figure S5). Analysis of photosynthesis characteristics in control, VIGS-vector, and VIGS-CtPPO at starting point. The photosynthesis rate, intercellular CO2, stomatal conductance, and transpiration rate were measured using an open gas-exchange system. Data are shown as mean ± SD of independent biological replicates. Asterisks indicate significant changes as measured by Student's t t-test (* p < 0.05, ** p < 0.01, and *** p < 0.001).

Effects of VIGS-CtPPO on ATP Synthase in C. terniflora DC. Leaves
To confirm the relationship between CtPPO and energy metabolism, qRT-PCR-based analysis of ATP synthase was performed. Leaves of VV and VC were collected at starting point and after HUV-B+D. At the starting point, the expression of ATP synthase was increased without statistical significance in VC compared with control and VV ( Figure 8). However, after HUV-B+D, ATP synthase was increased by 2-folds in VC plants compared with VV and control plants (Figure 8). These results indicate that the silencing of CtPPO might activate the energy metabolism by upregulating the ATP synthase in response to HUV-B+D. Figure 7. Analysis of photosynthesis characteristics in control, VIGS-vector, and VIGS-CtPPO at starting point. The photosynthesis rate, intercellular CO 2 , stomatal conductance, and transpiration rate were measured using an open gas-exchange system. Data are shown as mean ± SD of independent biological replicates. Asterisks indicate significant changes as measured by Student's t t-test (* p < 0.05, ** p < 0.01, and *** p < 0.001).

Effects of VIGS-CtPPO on ATP Synthase in C. terniflora DC. Leaves
To confirm the relationship between CtPPO and energy metabolism, qRT-PCR-based analysis of ATP synthase was performed. Leaves of VV and VC were collected at starting point and after HUV-B+D. At the starting point, the expression of ATP synthase was increased without statistical significance in VC compared with control and VV ( Figure 8). However, after HUV-B+D, ATP synthase was increased by 2-folds in VC plants compared with VV and control plants (Figure 8). These results indicate that the silencing of CtPPO might activate the energy metabolism by up-regulating the ATP synthase in response to HUV-B+D. Data are shown as mean ± SD of independent biological replicates. Asterisks indicate significant changes as measured by Student's t-test (** p < 0.01, and *** p < 0.001).

Silencing of PPO Promoted the Light Reaction in C. terniflora
The photosynthetic process in plants can operate through linear or cyclic electron flow, which involves three major complexes of the electron transfer chain: photosystem II (PSII), photosystem I (PSI), and cytochrome b6f complex [29,30]. Light reactions are an important step in photosynthesis and can transfer light energy to ATP and NADPH [31]. The initial step in this process is the absorption of light energy by the chlorophyll molecule. Ferredoxins were extensively employed as electron shuttles by anaerobes long before the advent of oxygenic photosynthesis [32], while the ferredoxin-NADP + reductase catalyzes an electron-hydride exchange between reduced ferredoxin and NADP + to yield NADPH [33]. The activation state of chloroplast ATP synthase is regulated by proton-motive forces generated by photosynthetic electron transfer reactions and reduction of disulphides in the γ subunit by TRX [34,35]. A lower ATP content resulting from the loss of ATP synthesis can inhibit the synthesis of ribulose biphosphate and influence the photosynthetic assimilation of CO2 [36]. In our study, the amounts of chlorophyll a and chlorophyll b were increased in VC plants compared with the VV and control plants, indicating the silencing of CtPPO could upregulate the PSs system in photosynthesis. In addition, the expression of ATP synthase, and ferredoxin-NADP + reductase were decreased in response to HUV-B+D (Table S6), indicating the inhibition effect of HUV-B+D on C. terniflora. However, the expression of ATP synthase ( Figure 8) and ferredoxin-NADP + reductase (Table 1)   Data are shown as mean ± SD of independent biological replicates. Asterisks indicate significant changes as measured by Student's t-test (** p < 0.01, and *** p < 0.001).

Silencing of PPO Promoted the Light Reaction in C. terniflora
The photosynthetic process in plants can operate through linear or cyclic electron flow, which involves three major complexes of the electron transfer chain: photosystem II (PSII), photosystem I (PSI), and cytochrome b6f complex [29,30]. Light reactions are an important step in photosynthesis and can transfer light energy to ATP and NADPH [31]. The initial step in this process is the absorption of light energy by the chlorophyll molecule. Ferredoxins were extensively employed as electron shuttles by anaerobes long before the advent of oxygenic photosynthesis [32], while the ferredoxin-NADP + reductase catalyzes an electron-hydride exchange between reduced ferredoxin and NADP + to yield NADPH [33]. The activation state of chloroplast ATP synthase is regulated by proton-motive forces generated by photosynthetic electron transfer reactions and reduction of disulphides in the γ subunit by TRX [34,35]. A lower ATP content resulting from the loss of ATP synthesis can inhibit the synthesis of ribulose biphosphate and influence the photosynthetic assimilation of CO 2 [36]. In our study, the amounts of chlorophyll a and chlorophyll b were increased in VC plants compared with the VV and control plants, indicating the silencing of CtPPO could upregulate the PSs system in photosynthesis. In addition, the expression of ATP synthase, and ferredoxin-NADP + reductase were decreased in response to HUV-B+D (Table S6), indicating the inhibition effect of HUV-B+D on C. terniflora. However, the expression of ATP synthase ( Figure 8) and ferredoxin-NADP + reductase (Table 1)

Silencing of PPO Activated the Calvin Cycle in C. terniflora
The Calvin cycle is an essential carbon fixation process in photosynthesis that converts carbon dioxide to glucose, accompanied by reduction reactions and ribulose 1,5-bisphosphate (RuBP) regeneration [37]. The Calvin cycle is redox-activated process, and the redox homeostasis is mediated by reducing power from photosynthetic electron transport to ferredoxin (Fd) and NADPH, via Fd-thioredoxin (TRX) reductase (FTR) and NADPH-dependent TRX reductase (NTRC) [38]. FTR and multiple TRXs consist the Fd-TRX system, while the NTRC contains a complete TRX system in a single polypeptide. Both systems have an impact on activation of Calvin cycle-related enzymes, including fructose-1,6-bisphosphatase (FBPase), sedoheptulose-1,7-bisphosphatase (SBPase), phosphoribulokinase (PRK), and GAPDH [39][40][41]. In our current study, Calvin cycle-related PGK, GAPDH, PRK, TK, and RuBP were decreased in response to HUV-B+D, indicating HUV-B+D significantly inhibited the Calvin cycle. However, the drop in CtPPO activity led to a dramatic increase in Calvin cycle-related proteins in C. terniflora at the starting point and after HUV-B+D. The redox-related proteins were further analyzed and TRX was found to be increased in VC plants compared with VV plants at the starting point (Table 2) and Fd-NADP reductase was increased in VC plants compared to VV plants after HUV-B+D ( Table 1). As PPO is located in chloroplast [42,43], these results indicate VIGS of CtPPO can protect the Calvin cycle by activating the FTR and NTRC system. Chlorophyll biosynthetic enzymes glutamyl-tRNA reductase, Mg-protoporphyrin IX monomethylester cyclase, and plastidic 2-Cys PRXs are direct substrates of NTRC in Arabidopsis [44]. As a result, NTRC protects the formation of chlorophyll in Arabidopsis [45]. The NTRC, together with 2-Cys PRXs, forms a two component peroxide detoxifying system that acts as a reductant under stress conditions [44]. Overexpression of NTRC raises the CO 2 fixation rate and lowers non-photochemical quenching by enhancing the activation of TRX-regulated enzymes and chloroplast ATP synthase in the Calvin cycle [46]. In this study, the amounts of total chlorophyll, chlorophyll a, and chlorophyll b were increased in VC plants compared with VV and control plants ( Figure 6). The amount of 2-Cys peroxiredoxin was decreased in VV plants in response to HUV-B+D (Table S6), however, it was increased in VC plants compared with VV plants after HUV-B+D (Table 1). Furthermore, chloroplast ATP synthase was increased in VC plants compared with VV and control plants at the starting point and were increased more significantly in VC plants than in VV and control plants following HUV-B+D (Figure 8). The intercellular CO 2 was significantly lower in VC than VV and control (Figure 7) which illustrated the silencing of CtPPO increased the utilization of carbon resource. These results help to increase understanding that CtPPO might influence the carbon fixation by regulating NTRC system.

Silencing of PPO Enhanced the glycolysis in C. terniflora
Glycolysis is a central metabolic pathway in plants that oxidizes hexoses to provide ATP, reduces power and pyruvate, and produces precursors for anabolism [47]. Studies have proved that the increase of PGK, and BPGM in Arabidopsis thaliana provide energy supply by ATP production [48,49]. Arabidopsis double mutants lacking BPGM enzyme activity exhibited defects in blue light, low CO 2 , and abscisic acid-regulated stomatal movements [49], which are responsible for all gaseous diffusion and can control photosynthetic CO 2 uptake to influence photosynthesis [50]. The increased glycolytic proteins in root of oat enhanced ATP production and promoted adaptation to anaerobic conditions [51]. It has been suggested that the stomatal red-light response signal may be provided by the redox state of photosynthetic electron transport chain components, such as the redox state of plastoquinone and production of ATP [52] In our study, the increases in GPUT, FBA, GAPDH, PGK, Enolase, and PEPC in CtPPO-silenced plants demonstrated the activation of glycolysis, which could provide more ATP for energy metabolism in C. terniflora. Further analysis also proved the enhanced stomatal conductance by the silencing of CtPPO in VC (Figure 7). The interaction of glyceraldehyde-3-phosphate dehydrogenase and phospholipase D might provide a direct connection between signal transduction, energy metabolism, and growth control in plant response to stress conditions [53]. Yasmeen's study suggest that Cu nanoparticles might enhance the tolerance of wheat to drought and salinity by increasing the glycolysis related protein abundance [54]. These results indicate the decrease in CtPPO activity in C. terniflora might elevate its stress tolerance by activating glycolysis metabolism.

Artificial Interference with PPO Activity Has Potential Applications in Agricultural Production
PPOs catalyze the oxidation of phenols to quinones that subsequently react with amino acids or proteins to form brown and black pigments, greatly reducing the appearance of quality of wheat products [55]. The role of plant PPO in postharvest browning has been the primary focus of research [56]. Although there is no evidence that high PPO activity is associated with depressed nutritional value, the darkened color still negatively affects consumer choice [57]. Studies on wheat concluded that the development of wheat cultivars with low grain PPO activity is one of the main objectives in wheat breeding programs [58]. Photosynthesis is the major trait for any further increase in the yield potential of crops [59]. Various genetic engineering approaches for enhancing C3 plant photosynthesis have been consistently proposed over the past years to improve the crop productivity [60][61][62]. Our study clearly suggests the potential for suppression of PPO activity on photosynthesis in C. terniflora, which indicates the feasibility of artificially mediating photosynthesis through the control of PPO gene expression. Nevertheless, it is undeniable that PPO has positive effects on plant defense in response to abiotic or biotic stresses as many PPO genes have been shown to be upregulated by wounding, pathogens, and hormones [63]. However, latent PPO enzymes can be activated by interactions with their substrates [64], occurring in great measure when plants are under stress. All in all, our current work provides a prospective approach for increasing crop yield in agricultural production.

Plant Materials and Growth Conditions
C. terniflora DC. seeds sprouted in incubators and the sprouts were sown into seedbeds. The seedlings were then transplanted into pots and placed in a greenhouse, which was controlled at 28-30 • C, 70-80% relative humidity, and 160 µmol m −2 s −1 of white light irradiance. In the garden, conditions ranging from the soil to microclimate were equivalent among all plant samples. After six weeks, the plants were used for experiments [28].

High Level UVB and Dark Treatment
For HUV-B+D, 6-week-old plants were exposed to 104.4 kJ m −2 d −1 of UV-B irradiation at conditions of 25-30 • C and 80% relative humidity in a cabinet. Plants before irradiated were regarded as starting point samples. Intensity of UV-B irradiation was determined by a UV Light Meter (Beijing Normal University, Beijing, China). After irradiated by HUV-B for 5 h, plants were incubated in dark for 48 h. Leaves in the basal 10 to 60 cm of each experimental plant were collected as simples for further proteomic analyses. Five leaves were collected for each replicate and 3 independent biological replicates were assessed [28].

RNA Extraction and Cloning of CtPPO Gene
An RNeasy Plant Mini kit (Qiagen, Hilden, Germany) was used to obtain total RNA. Quantity and quality of obtained RNA were determined using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). The extracted total RNA was used as template to synthesize first-strand complementary DNA (cDNA) using an Oligo (dT) and OneScript TM Reverse Transcriptase OneScrpt TM cDNA Synthesis Kit (Applied Biological Materials Inc., Vancouver, BC, Canada). The first-strand cDNA was utilized as template for PCR using the primers listed in Table S5, which was designed based on the single EST sequence of a suspected CtPPO gene in the transcriptome of C. terniflora and amplification was performed as described [28]. The amplified PCR fragments, which were approximately 264 bp in length, were isolated, inserted into the pMD TM 18-T Vector (pMD TM 18-T Vector Cloning Kit, Takara, Kyoto, Japan), and sequenced. The complete coding sequence of CtPPO was obtained using 5and 3 -rapid amplification of cDNA ends (RACE) with internal specific primers (Table S3) using the SMARTer RACE cDNA Amplification Kit (Takara) according to the manufacturer's instructions.
The downstream sequence of CtPPO was obtained using a PCR-based genome walking approach with the Genome Walking Kit (Takara). Three designed primers (gene specific primer (GSP); Table S2) and 4 internal short degenerate arbitrarily primed (AP) primers included in the kit (AP1, AP2, AP3, and AP4) were used to perform genome-walking PCR. The PCR conditions and designed primer Tm were based on the manufacturer's instructions. Three-step nested PCR was performed to increase specificity. The nested PCR products were excised, purified, and sequenced.

Construction of CtPPO Virus-Derived Vectors
A 439-bp fragment of the obtained CtPPO gene was amplified using Permix TaqTM (Ex TaqTM Version 2.0 plus dry; Takara) and inserted into a pMD TM 18-T vector (pMD TM 18-T Vector Cloning Kit, Takara) using primers TXPPO-F and TXPPO-R (Table S4) containing XbaI and BamHI restriction enzyme sites. The CtPPO sequence was harvested from the pMD TM 18-T vector using the EcoRI restriction enzyme and inserted into the EcoRI site between the CaMV 35S promoter (2 × 35S) and NOS terminator of the pTRV2 vector. The orientation of the plasmid with pTRV2-CtPPO ( Figure S1) was verified by sequencing. Assembly of the pTRV vector was performed as described [65]. The plasmids were sequenced to verify correct insertion of the fragment and transformed into Agrobacterium tumefaciens GV3101.

Virus-Induced Gene Silencing
VIGS was performed according to the method described by Salim et al. [66]. The GV3101 strains of A. tumefaciens carrying pTRV1, pTRV2, and pTRV2-PPO were stored at −80 • C. Agrobacteria were cultured in 300 mL of Luria-Bertani liquid medium containing 10 mM 2-(N-morpholino) ethanesulfonic acid (MES), 50 µg mL −1 kanamycin, and 20 µM acetosyringone. Then the cultures were centrifuged at the speed of 5000× g for 10 min and the resulting bacterial pellets were resuspended in 5 mL infiltration buffer containing 10 mM MgCl 2 , 10 mM MES and 200 µM acetosyringone, and incubated at 28 • C with shaking for 3 h. The mixed suspensions were used for surface infiltration of three-week-old Nicotiana benthamiana seedlings using a needleless syringe and 4-6 leaves were infected. The PDS (phytoene desaturase) gene was used as a marker for the evaluation of the efficiency of VIGS in CtPPO genes silencing in N. benthamiana. The VIGS N. benthamiana plants were grown in a greenhouse and the TRV virus was detected by PCR to confirm infection of the N. benthamiana leaves. The successfully infected N. benthamiana leaves were used for the abrasion inoculation of C. terniflora leaves as described below.
Virus-infected N. benthamiana leaves (1 g) were ground with 2 mL of 10 mM potassium phosphate buffer (pH 7) containing 1-2% w/v Celite 545 AW abrasive (Sigma-Aldrich, Milwaukee, WI, USA) using a mortar and pestle. The resulting homogenate was used for inoculation by the abrasion method on all tender leaves of 3-week-old C. terniflora plants. Two groups of plants, one containing the pTRV virus vector and the other the pTRV-PPO virus vector, were prepared and each was inoculated, where six plants were used as biological replicates.

Protein Extraction, Enrichment, and Digestion for Proteomics Analysis
The tissues were ground in liquid nitrogen. For protein extraction, the phenol extraction method was performed according to a previously published protocol [67]. Briefly, 6 g of frozen powder was suspended in 15 mL of homogenization medium containing 50% (v/v) phenol, 100 mM KCl, 50 mM EDTA, 0.2% (v/v) 2-mercaptoethanol, 1 mM PMSF, 700 mM sucrose, and 500 mM Tris-HCl at pH 7.5, shaken at 37 • C for 15 min, and then centrifuged at the speed of 2500× g at 4 • C for 20 min. Four volumes of precipitation solution (0.1 M ammonium acetate in 100% methanol) were added in to the phenolic phase collected to precipitate the proteins. The precipitate was collected by centrifugation for 10 min at the speed of 2500× g. The supernatant was discarded and washed with precipitation solution. Precipitation was allowed to occur at −20 • C for 3 h and the precipitate was obtained by centrifugation, which was repeated twice. The collected proteins were dissolved in 0.5 mL of 50 mM ammonium bicarbonate (ABC). Protein concentrations were determined using the Bradford assay with bovine serum albumin as the standard.
The protein (100 µg) and an additional 50 mM ABC were transferred into a new tube to create a final volume of 100 µL. The protein was reduced with 10 mM dithiothreitol for 1 h at 55 • C and alkylated with 50 mM iodoacetamide at 37 • C for 30 min in darkness, then digested with sequence-grade modified trypsin (Promega, Fitchburg, WI, USA). Trypsin and 1% lysyl endopeptidase were added to the protein and then the solutions were incubated in a water bath at 37 • C for 4 h and another 1% enzyme was added. The solutions were incubated at 37 • C for 12 h, centrifuged at 10,000× g, and then dried in a vacuum. Trypsinized peptides in 200 µL of 0.1% trifluoroacetic acid (TFA) and 0.5% acetonitrile were loaded onto activated and balanced C18 SPE columns (Sep-Pak C18; Waters, Milford, MA, USA). After washing the column with 200 µL of 0.1% TFA and 0.5% acetonitrile, 300 µL of 0.1% TFA and 60 % acetonitrile were added and then the elution solution was collected and dried in a vacuum. Obtained supernatant was collected for subsequent analysis.

Nano-HPLC-MS/MS Analysis
The peptides were resuspended in 100 µL of solvent A (water with 0.1% formic acid; B: ACN with 0.1% formic acid), separated by nano-UPLC, and analyzed by online electrospray tandem MS, which were performed on a Nano Aquity UPLC system (Waters Corporation, Milford, MA, USA) connected to a quadrupole-Orbitrap mass spectrometer (Q-Exactive; Thermo Fisher Scientific, Bremen, Germany) equipped with an online nano-electrospray ion source. Peptide samples (2 µL) were loaded onto the trap column (Thermo Scientific Acclaim PepMap C18, 100 µm × 2 cm) with a flow of 10 µL/min for 3 min and subsequently separated on an analytical column (Acclaim PepMap C18, 75 µm × 25 cm) with a linear gradient of 1 to 30% solvent B (acetonitrile with 0.1% formic acid) in 95 min. The column was re-equilibrated at the initial conditions for 15 min. The column flow rate was maintained at 300 nL/min. The electrospray voltage of 2.0 kV versus the inlet of the mass spectrometer was used.
The Q-Exactive mass spectrometer was set at the data-dependent mode allowing the machine switch between MS and MS/MS acquisition automatically. Survey full-scan MS spectra (m/z 350-1600) were acquired with a mass resolution of 70 K followed by fifteen sequential high-energy collisional dissociation MS/MS scans with a resolution of 17.5 K. One microscan was recorded using a dynamic exclusion of 30 s in all cases. The MS/MS fixed first mass was set at 100.

Protein Identification Based on Mass Spectrometry Data
Tandem mass spectra were extracted by Proteome Discoverer software (Thermo Fisher Scientific, version 1.4.0.288, Bremen, Germany). Charge state deconvolution and deisotoping were not performed. All MS/MS samples were analyzed using Mascot (Matrix Science, London, UK; version 2.3). Mascot was set up to search the Uniprot database (Taxonomy: Viridiplantae, 5716577 entries, https://www. uniprot.org/) assuming the digestion enzyme trypsin. Mascot was searched with a fragment ion mass tolerance of 0.050 Da and a parent ion tolerance of 10.0 PPM. Carbamidomethyl of cysteine was specified in Mascot as fixed modifications. Oxidation of methionine was specified in Mascot as a variable modification. Minimum peptide length was set to six amino acids the false discovery rate was set to 0.01 for peptide identifications. Identified peptides shared between two proteins were combined and reported as one protein group. The minimum requirement for the identification of a protein was a minimum of 2 matched peptides and a p-value < 0.05. Spectral counting was used to estimate relative protein abundance [68]. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [69] partner repository with the dataset identifier PXD011439.

Functional Annotation
Protein functions were categorized using MapMan bin codes (available online: http://mapman. gabipd.org/) [70]. Annotations were transferred to the Arabidopsis genome with consideration for orthologous genes to predict functions of identified proteins from C. terniflora. Pathway mapping of identified proteins was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (available online: http://www.genome.jp/kegg/) [71].

qRT-PCR
Total RNA of leaves was extracted using an RNA isolation kit (Huayueyang, Beijing, China) and served as the template for reverse transcription using 5× All-In-One RT MasterMix with an AccuRT Genomic DNA Removal Kit (Applied Biological Materials, Richmond, BC, Canada) to obtain cDNA. Next, qRT-PCR was performed on a Bio-Rad IQ2 Multicolor Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA) with EvaGreen 2× qPCR Master Mix-iCycler (Applied Biological Materials) as the fluorescent dye. GAPDH, a housekeeping gene, was used as a standard for relative quantification of target genes. The gene-specific primers are listed in Additional file 1: Table S1. Three biological replicates were assessed and the relative expression levels were calculated using the 2 −∆∆Ct method [26].

PPO Enzymatic Activity
Enzymatic activity of PPO was assayed as described by Marko et al. [75]. Frozen leaf (0.5 g) was homogenized with 1 mL of 50 mM sodium phosphate buffer (pH 7.5) containing 1% (w/v) polyvinyl polypyrrolidone and 0.1 mM EDTA using a pre-chilled mortar and pestle. After extraction in an end-over-end shaker at 4 • C for 20 mins, the mixture was centrifuged at 16000× g for 20 min at 4 • C and the resulting supernatant was used in enzyme extraction.
Then, 800 µL of 0.1 M sodium phosphate buffer (pH 6.0) containing 50 mM (+)-catechin (Carl Roth, Karlsruhe, Germany) was mixed with 150 µL of enzyme extraction sample and incubated for 15 min at 37 • C. The reaction was terminated by adding 150 µL of 6 N HCl and the absorbance at 420 nm was recorded and used in further calculations. Protein concentration was determined according to Bradford using bovine serum albumin as the reference protein. Specific PPO activity was calculated as units (U) per mg sample protein with 1 U corresponding to an absorbance change of 0.01/min. Since levels of PPO activity in the control plants were somewhat variable between independent experiments, data were normalized to the specific PPO activity in leaves of healthy wild-type controls and the results are shown as fold induction (with 1 corresponding to a specific activity of 3.2 to 23 U/mg).

Analysis of Photosynthesis Characteristics
An Li-6400 Portable Photosynthesis System (Li-Cor, Lincoln, NE, USA) equipped with a red/blue LED light source were used to measure the photosynthesis characteristics of the leaves, and was operated using air from a large volume with a stable CO 2 partial pressure. All measurements were carried out at a photon flux density of 500 µmol m −2 s −1 at 28 • C and made after stable reading was achieved. The photosynthesis rate (µmol m −2 s −1 ), stomatal conductance (mol m −2 s −1 ), and transpiration rate (mmol m −2 s −1 ) of each leaf were recorded. Each experiment was repeated for different plants three times and three independent biological replicates were assessed.

Analysis of Chlorophyll Content
C. terniflora leaves (0.2 g) were cut up and blended with 10 mL of 95% alcohol in a tube. The tube was covered and the leaves soaked in the dark until they had turned completely white.

Statistical Analysis
Statistical significance comparing two groups was determined using Student's t-test and multiple groups using one-way ANOVA followed by Tukey's test. Statistical evaluation was performed on SPSS statistical software version 22.0 (IBM, Armonk, NY, USA). A p-value <0.05 was considered statistically significant.

Accession Codes
The nucleotide sequence of CtPPO has been deposited in the GenBank™ database under the accession number MK070494.

Conclusions
PPO is an oxidoreductase that plays a crucial role in the biosynthesis of secondary metabolites. C. terniflora is a Chinese folk medical resource with potential pharmaceutical value for the treatment of various diseases [21]. High-intensity UV-B and dark incubation (HUV-B+D), which are helpful stressors for medicinal plants [28]. Omics technologies were used on C. terniflora to prospectively understand the inherent mechanism underlying its medicinal quality [24,26]. To clarify the in vivo physiological functions of PPO, we performed a comparative proteomic analysis on the leaves of C. terniflora with down-regulated PPO activity by VIGS. The main findings are as follows: (i) The differentially expressed proteins were mainly functioned in photosynthesis, glycolysis and redox in the VC plants; (ii) the differentially expressed proteins related to photosynthesis were mainly involved in light reaction and Calvin cycle; (iii) the expression level of adenosine triphosphate (ATP) synthase, the content of chlorophyll, and the photosynthesis rate were increased in VC plants compared to VV plants pre-and post HUV-B+D. Taken together, these results indicate that the silence of PPO can activate the glycolysis process, regulate Calvin cycle and provide ATP for energy metabolism to elevate the plant photosynthesis. And this study provides a prospective approach for increasing crop yield in agricultural production.