Original Research ARTICLE
Proteomic Analysis Reveals the Positive Roles of the Plant-Growth-Promoting Rhizobacterium NSY50 in the Response of Cucumber Roots to Fusarium oxysporum f. sp. cucumerinum Inoculation
- 1Key Laboratory of Southern Vegetable Crop Genetic Improvement in Ministry of Agriculture, College of Horticulture, Nanjing Agricultural University, Nanjing, China
- 2Department of Horticulture, Shanxi Agricultural University, Taigu, China
- 3Suqian Academy of Protected Horticulture, Nanjing Agricultural University, Suqian, China
Plant-growth-promoting rhizobacteria (PGPR) can both improve plant growth and enhance plant resistance against a variety of environmental stresses. To investigate the mechanisms that PGPR use to protect plants under pathogenic attack, transmission electron microscopy analysis and a proteomic approach were designed to test the effects of the new potential PGPR strain Paenibacillus polymyxa NSY50 on cucumber seedling roots after they were inoculated with the destructive phytopathogen Fusarium oxysporum f. sp. cucumerinum (FOC). NSY50 could apparently mitigate the injury caused by the FOC infection and maintain the stability of cell structures. The two-dimensional electrophoresis (2-DE) approach in conjunction with MALDI-TOF/TOF analysis revealed a total of 56 proteins that were differentially expressed in response to NSY50 and/or FOC. The application of NSY50 up-regulated most of the identified proteins that were involved in carbohydrate metabolism and amino acid metabolism under normal conditions, which implied that both energy generation and the production of amino acids were enhanced, thereby ensuring an adequate supply of amino acids for the synthesis of new proteins in cucumber seedlings to promote plant growth. Inoculation with FOC inhibited most of the proteins related to carbohydrate and energy metabolism and to protein metabolism. The combined inoculation treatment (NSY50+FOC) accumulated abundant proteins involved in defense mechanisms against oxidation and detoxification as well as carbohydrate metabolism, which might play important roles in preventing pathogens from attacking. Meanwhile, western blotting was used to analyze the accumulation of enolase (ENO) and S-adenosylmethionine synthase (SAMs). NSY50 further increased the expression of ENO and SAMs under FOC stress. In addition, NSY50 adjusted the transcription levels of genes related to those proteins. Taken together, these results suggest that P. polymyxa NSY50 may promote plant growth and alleviate FOC-induced damage by improving the metabolism and activation of defense-related proteins in cucumber roots.
Cucumber (Cucumis sativus L.) is an important and popular vegetable cash crop that is consumed worldwide. However, the production of this plant is severely threatened by cucumber Fusarium wilt, which is caused by the soilborne fungal pathogen Fusarium oxysporum f. sp. cucumerinum (FOC), a destructive vascular disease that can cause plant death and serious economic loss (Ahn et al., 1998; Huang et al., 2012). Because there are no available effective chemical products or resistant varieties and because grafting is time consuming (Chung et al., 2008; Cao et al., 2012), other effective, economical, and environmentally friendly cultivation methods have been developed and tested.
To date, alternative methods, such as biological control agents, have been shown to be effective and have been increasingly applied in the field (Yang et al., 2010; Li et al., 2012; Lin et al., 2014; Xu et al., 2014). Among these methods, some beneficial bacteria inhabiting the plant rhizosphere, called plant growth-promoting rhizobacteria (PGPR), are directly or indirectly involved in promoting plant growth and the biological control of plant diseases (Kloepper and Metting, 1992). A variety of mechanisms have been confirmed to improve plant health and increase crop productivity for those strains, such as synthesizing different antimicrobial compounds (Tamehiro et al., 2002; Ongena et al., 2007; Chen et al., 2014), various hormones (Bottini et al., 2004; Martínez-Viveros et al., 2010; Kochar et al., 2011), increasing the availability of plant nutrients in the rhizosphere (Idriss et al., 2002; Palacios et al., 2014), and inducing systemic resistance in host plants (Niu et al., 2011; Jiang et al., 2015).
Numerous technical methods, such as DNA microarray technology and proteomic analysis, have been applied to improve our understanding of plant-PGPR interactions. A study in the Arabidopsis root transcriptome showed that extensive changes were induced in hormone- and defense-related genes after exposure to Azospirillum brasilense Sp245 (Spaepen et al., 2014). In another study, the ISR- and iron acquisition-related transcription factors in Arabidopsis roots were activated by Pseudomonas fluorescens WCS417 and its volatiles (Zamioudis et al., 2015). A dynamic protein network was designed to explain the mechanism of P. putida UW4 to release hypoxic stress and promote cucumber growth (Li et al., 2013b). Proteins involved in growth promotion and defense reactions were significantly induced in Arabidopsis following inoculation with Paenibacillus polymyxa E681 (Kwon et al., 2016). In addition, many other researchers have used proteomics to investigate plant-pathogen interactions. Wang et al. (2011) found that a switch from glycolysis to the pentose phosphate pathway may be an active response of the cotton plant against V. dahliae infection to enhance wilt tolerance or resistance. Total root protein was isolated from infected cucumber roots of the susceptible bulk and resistant bulk of cucumber generation F2, and the result showed that jasmonic acid and redox signaling components occurred in response to F. oxysporum infection in resistant plants (Zhang et al., 2016). Eighty-six differentially expressed proteins were identified from the leaves of resistant tomato cultivar “Zheza-301” and susceptible cultivar “Jinpeng-1” after TYLCV infection, which demonstrated that an interaction network between tomato leaves and TYLCV infection was established (Huang et al., 2016). To date, although several papers have addressed the use of proteomics to investigate the interactions of plants and pathogens (Konishi et al., 2001; Wang et al., 2011; Carrillo et al., 2014) or plants and PGPR (Kandasamy et al., 2009; Wang et al., 2013), little research is available on the effects of PGPB on cucumber plants under pathogenic attack.
In this present study, a promising bio-control agent Paenibacillus polymyxa-NSY50, which was originally isolated from the high suppression capacity of compost (Du et al., 2015; Shi et al., 2016), was used to better understand direct or indirect interactions among P. polymyxa NSY50, the pathogen FOC and cucumber roots. A total of 56 proteins were identified that were up- or down-regulated by the inoculation of P. polymyxa NSY50 and/or FOC, and the proteins could be classified into multiple biological functions. This work provides key insights into the promotion mechanisms of P. polymyxa NSY50 action to determine the optimal use of this beneficial PGPR in sustainable agricultural practices.
Material and Methods
Plant Material, Microbial Culture Conditions, and Treatments
Cucumber (Cucumis sativus L. cv. Jinchun No. 2) seeds were sown in quartz sand and cultivated in a greenhouse. Paenibacillus polymyxa-NSY50 was grown on LB medium at 28°C for 3 days, and the cucumber Fusarium wilt pathogen F. oxysporum f. sp. cucumerinum (FOC) was incubated in PDA liquid culture for 7 days.
At the two true leaves stage, the cucumber seedlings were transplanted to tanks containing half-strength Hoagland nutrient solution (He et al., 2015). After 3 days of pre-culture, the seedlings were treated as follows: (1) control, control plants were grown in Hoagland's solution; (2) NSY50, seedlings grown in Hoagland nutrient solution containing 500-mL 1.0 × 108 NSY50 cell suspension; (3) FOC: 500 mL of a cell suspension of FOC at 1 × 107 conidia/mL was poured into the tanks after 6 days of culture; and (4) NSY50+FOC: plants were inoculated with 500 mL of NSY50 (1 × 108 CFU/mL) after 3 days of culture, and 3 days later, the plants were challenge-inoculated with 500 mL of a cell suspension of FOC (1 × 107 conidia/mL). The experiment was arranged in a randomized complete block design with three replicates for each treatment, which resulted in a total of 36 seedlings per treatment. Seedlings were reared in a greenhouse at 25–30°C in the daytime and at 15–18°C at night under natural lighting at 60–75% relative humidity.
Plant Growth Analysis
After 9 days of FOC treatment, shoot, and root were separated from one another and measured after being washed with sterile distilled water. The plant height, biomass, and volume of root in each treatment were measured as previously described (Shi et al., 2016). Each data point was the average of measurements collected in triplicate.
For transmission electron microscopy, roots of seedlings from different treatments were excised and immediately cut into small segments (1~2 mm, including the tip); The samples were then fixed with 2.5% glutaraldehyde in 0.1 M phosphate buffer (pH 7.4) for 24 h (primary fixation) and immersed in 2% osmic acid in the same buffer for 2 h (second fixation). After dehydration in acetone and embedding in Durcupan ACM (Fluka), the resulting roots were cut to obtain ultra-thin sections, stained with uranium acetate and lead citrate in series and examined using a HITACHI transmission electron microscope (Carl Zeiss, Göttingen, Germany) at an accelerating voltage of 80 kV.
The cucumber roots were prepared for protein fractionation 3 days post-inoculation of FOC. Total protein extraction was performed using a trichloroacetic acid, i.e., the acetone precipitation method modified described by Hurkman and Tanaka (1986). Fresh root samples (2 g) were milled with liquid nitrogen and suspended in extraction buffer as described by An et al. (2016). The homogenate was centrifuged at 15,000 g at 4°C for 20 min. An aliquot (1 mL) of the resulting supernatant was transferred to a new tube and precipitated with ice-cold acetone (10% TCA and 0.07% β-mercaptoethanol) overnight at −20°C, and then the resulting protein sample was centrifuged at 20,000 g for 25 min. The pellet was washed three times with cold acetone (0.07% β-mercaptoethanol) and allowed to stand at −20°C for 2 h. Finally, the protein pellet was air-dried and used for 2-DE.
2-DE and Image Analysis
Isoelectric focusing (IEF) was performed using an 18-cm IPG linear gradient strip, pH 4–7 (GE Healthcare, San Francisco, USA). The dried protein pellet was rehydrated in a rehydration buffer that contained 7 M urea, 2 M thiourea, 4% 3-[(3-cholanidopropyl) dimethylammonio]-1-propanesulfonic acid (w/v), 40 mM DTT, 0.5% (v/v) IPG buffer 4–7 and 0.01% (w/v) bromophenol blue. The protein levels were quantified using the Bradford method (Bradford, 1976). The 350 μL protein samples were loaded on IPG strips containing 800 μg of protein and rehydrated for 12–16 h at 25°C. After rehydration, the IPG strips were run on an Ettan IPGphor 3 (GE Healthcare, San Francisco, USA) using the following protocol: the voltage for IEF was set at 100 V for 1 h, followed by 200 V for 1 h, 200 V for 1 h, 500 V for 1 h, 1000 V for 1 h, 4000 V for 1 h, a gradient of 10,000 V for 1 h, and then a 10,000 V rapid focus, which reached a total of 75,000 V h with a maximum electric current of 50 μA per strip. After accomplishing the first dimension, IEF strips were equilibrated for 15 min in a 2D equilibrium buffer that consisted of 6 M urea, 30% glycerol (v/v), 50 mM Tris-HCl (pH 8.8), and 2% sodium dodecyl sulfate (SDS) containing 1% DTT (w/v); the strips were then incubated in 2.5% (w/v) iodoacetamide instead of DTT for 15 min. Then, the strips were placed directly onto 12.5% polyacrylamide-SDS slab gels and sealed using 1% molten agarose solution. The second dimensional separation was conducted using the EttanDalySix electrophoresis system (GE Healthcare, USA). At 15 mA per gel, the electrophoresis was not stopped until the bromophenol blue dye reached approximately 1 cm from the bottom of the gel. Protein spots were visualized using Coomassie Brilliant Blue (CBB) R-250.
The CBB-stained 2-DE gels were scanned using an Image scanner III (GE Healthcare, USA). The digitized images of three independent experiments were analyzed with Imagemaster 2D Platinum version 5.0 (GE Healthcare, USA). The abundance of each protein spot was estimated by the percentage volume (vol.%), and the spot volumes were normalized as the ratio of the total volume for all of the spots that were present. Spots with a vol.% that represented at least a 1.5-fold change and was significant at the P < 0.05 level (Duncan's multiple range tests) were used for identification.
Protein Identification and Functional Classification
The differentially expressed protein spots were extracted from CBB-stained preparative polyacrylamide gels and then identified using an ABI 5800 Proteomics Analyzer MALDI-TOF/TOF (Applied Biosystems, Foster City, CA, USA). The data lists were used as a query to search in the NCBI (http://www.ncbi.nlm.nih.gov/) and cucumber genomics databases (http://cucumber.genomics.org.cn) using the software MASCOT version 2.2 (Matrix Science, London, UK). The following search parameter criteria were used: trypsin cleavage, one missed cleavage site allowed; carbamidomethyl set as a fixed modification; oxidation of methionines allowed as a variable modification; a peptide mass tolerance within 100 ppm; fragment tolerance set to ± 0.4 Da; and a minimum ion score confidence interval for MS/MS data set to 95%.
Identified proteins were categorized according to the biological processes with which they were involved based on Gene Ontology (GO) (http://www.geneontology.org/) and the Uniprot Protein Knowledgebase (http://www.uniprot.org/). Hierarchical clustering of the protein expression patterns was performed on the log-transformed (using log base 2) spot abundance ratios using the software Cluster version 3.0. A heat map was created using Java Treeview.
Western Blot Analysis
Root tissues were milled in a mortar with ice-cold extraction buffer containing 30 mM Tris—HCl (pH 8.7), 1 mM MgCl2, 0.7 M saccharose, 1 mM EDTA, 1 mM DTT, 1 mM PMSF, and 1 mM ascorbic acid. The extracted protein was quantified using the Bradford method (1976), denatured at 95°C for 3–5 min and then stored at −20°C until analysis.
For western blot analysis, SDS-polyacrylamide gel electrophoresis (SDS-PAGE) was conducted using the methods of Laemmli (1970). Briefly, 11 μg protein samples were transferred to a 0.45 μm PVDF membrane at 10 V for 1.5 h and were washed with TBST three times; then, the PVDF membrane was blocked with 5% nonfat dry milk for 2 h, washed with TBST three times, and incubated with monoclonal antibodies against S-adenosylmethionine synthase and enolase (produced in rabbit; Univ-bio, Shanghai, China) for 2 h. The membrane was washed with TBST and incubated at room temperature for 1 h with a Goat Anti-Rabbit IgG HRP-conjugate. The membrane was then washed with TBST three times and developed using diaminobenzidene (DAB) and H2O2 (He et al., 2012).
Total RNA Extraction and Quantitative Real-Time PCR (qRT-PCR) Analysis
Root tissues were harvested at 1, 3, 9 days after inoculation with FOC. The total RNA was extracted according to the TRI reagent protocol (Takara Bio Inc.) and then converted into cDNA following the manufacturer's instructions. Primers were designed based on the sequences obtained from NCBI using Beacon Designer 7.90 (Supplementary Table 1). qRT-PCR was performed using the SYBR PrimeScript™ RT-PCR Kit (Takara Bio Inc.) in accordance with the manufacturer's instructions and run on a StepOne™ real-time PCR system (Applied Biosystems, Singapore). The relative gene expression was calculated using the 2−ΔΔCt method (Livak and Schmittgen, 2001), and the methods were based on a previously described protocol (Shi et al., 2016). All reactions were conducted with three biological replicates.
All data were statistically analyzed using SPSS 20.0 for Windows, and significance was assigned at the P < 0.05 level using Duncan's multiple comparisons test.
The biomass of the cucumber seedlings was measured 9 days post-inoculation (dpi) with the pathogen. The plant height, shoot freshness and dry mass of NSY50-treated seedlings were significantly greater than those of plants grown under control conditions (Table 1). However, pathogen inoculation (treatment FOC) seriously decreased the growth of the cucumber seedlings, the plant height, shoot fresh weight, dry weight and root volume and the root fresh weight and dry weight, which were significantly suppressed compared to the control. Moreover, under pathogen stress, plant biomass was significantly higher in P. polymyxa NSY50 pre-treated cucumber plants, which showed a 1.2–1.5-fold up-regulation in plant growth compared to untreated plants. Thus, the results clearly indicated that the pathogen inoculation inhibited growth in cucumber seedlings and that P. polymyxa NSY50 can release at least a portion of the stress and promote plant growth.
Table 1. Effects of plant-growth-promoting bacteria on the growth of cucumber seedlings under Fusarium wilt stress.
Ultrastructure of the Roots
Under control conditions, the vacuoles had small sizes, and many nuclei and mitochondria could be seen clearly under the microscope (Figure 1). The cellular ultrastructure of cucumber roots seemed to show no difference between the control plants and plants treated with NSY50 alone (NSY50). However, the inoculation of FOC (FOC) caused remarkable structural changes. Specifically, the karyotheca appeared to have contracted, and the number of the mitochondria decreased compared to the control. Furthermore, the number of central vacuoles and autophagosomes appeared to indicate that the pathogen attack advanced the cell aging and death process. Pre-treatment with NSY50 (NSY50+FOC) apparently enhanced the stability of the cell structure.
Figure 1. Effects of NSY50 on the cellular ultrastructure of cucumber root tips with or without FOC inoculation 3 days after treatment. Samples comprised 3 mm root tips. Control, seedlings cultured in normal nutrient solution; NSY50, control+ P. polymyxa NSY50 (2.5 × 106 CFU/mL); FOC, control+ FOC (2.5 × 105 conidia/mL); NSY50+FOC, pre-treated with NSY50 (2.5 × 106 CFU/mL) for 3 days, then challenged with FOC (2.5 × 105 conidia/mL).M, Mitochondria; N, Nucleus; V, Vacuole; P, Plastid; A, Autophagosome.
Identification and Functional Classification of Proteins
Approximately 350 reproducible protein spots were detected on the 2-DE gels, and 56 differentially expressed protein spots (changes ≥ 1.5-fold) were successfully identified by MALDI-TOF/TOF MS (Figure 2; Supplementary Figure 1). These differentially expressed proteins are listed in Table 2.
Figure 2. Representative 2-DE gel images of total protein extracts from root samples inoculated with NSY50 and/or FOC. An equal amount (800 μg) of total proteins were separated by IEF/SDS-PAGE and then stained with Coomassie Brilliant Blue (R-250) and loaded on each 18 cm gel strip (pH 4–7, linear). The numbers for the 56 differentially expressed proteins are marked and annotated according to the numbering in Table 2. Control, seedlings cultured in normal nutrient solution; NSY50, control+ P. polymyxa NSY50 (2.5 × 106 CFU/mL); FOC, control+ FOC (2.5 × 105 conidia/mL); NSY50+FOC, pre-treated with NSY50 (2.5 × 106 CFU/mL) for 3 days then challenged with FOC (2.5 × 105 conidia/mL).
The identified proteins were grouped into seven different categories on the basis of biological processes to which they were related according to the Gene Ontology and Uniprot Protein Knowledgebase (Figure 3A). The categories with a high level of expression variation included the defense response (28.6%), protein metabolism (25.0%), carbohydrate and energy metabolism (19.6%) and amino acid metabolism (17.9%). Among the 56 differentially expressed spots, 44 protein spots were significantly regulated by the pathogen FOC compared to the control (Figure 3B), and there were 16 up-regulated spots and 28 down-regulated spots (Figure 3C). The defense-response-related proteins had the highest up-regulation rate, and the decreased proteins were involved in protein metabolism (10 spots), carbohydrate, and energy metabolism (7 spots), and amino acid metabolism (7 spots). Additionally, 32 protein spots were significantly induced by NSY50 inoculation compared to the control (Figure 3B). Of those spots, 25 were up-regulated, and seven were down-regulated. The defense-response-related proteins (10 spots) and amino acid metabolism (5 spots) were the most highly enriched categories (Figure 3C). However, compared to FOC, 33 proteins had changes, including 24 proteins that increased and 9 proteins that decreased after pre-treatment with NSY50 (NSY50+FOC). The category of carbohydrate and energy metabolism had the highest up-regulation rate, whereas the decreased proteins were distributed among the defense response, protein metabolism, secondary metabolism and cell-related protein categories (Figure 3C).
Figure 3. Distribution of differentially expressed proteins by NSY50 and/or FOC in cucumber roots. (A) Functional classification and distribution of all 56 differentially expressed proteins. (B) Venn diagram showing the number of overlapping proteins that were differentially regulated by NSY50, FOC, and NSY50+FOC compared to the control. (C) Functional protein distribution in the compared groups (changes ≥1.5-fold or ≤0.67-fold).
Clustering Analysis of Differentially Expressed Proteins
To acquire a comprehensive overview of the differentially expressed proteins that were affected by NSY50 under control conditions and FOC stress, hierarchical clustering was performed, and proteins with similar expression patterns were grouped together (Figure 4). Cluster A consisted of 4 proteins (spots 9, 6, 13, and 10) that were up-regulated by NSY50, down-regulated after inoculation with FOC compared to the control, but recovered by the application of NSY50 (NSY50+FOC). Cluster B was composed of 23 proteins that decreased considerably under FOC stress with/without NSY50 inoculation. Most of the identified proteins in this cluster were involved in protein metabolism, carbohydrate and energy metabolism. Cluster C involved 7 proteins that were up-regulated by NSY50 under both control conditions and FOC stress conditions, and these proteins corresponded to defense responses, carbohydrate and energy metabolism, amino acid metabolism, and protein metabolism. Cluster D included 10 proteins that were down-regulated by NSY50 but up-regulated by the inoculation of FOC. Most of these proteins were involved in defense responses. Cluster E contained 4 proteins that were down-regulated by NSY50 but up-regulated by the inoculation of FOC, regardless of whether they were pre-inoculated with NSY50. Cluster F included 8 proteins that increased in abundance when inoculated with NSY50 or FOC alone but were down-regulated in combined inoculation treatment (NSY50+FOC).
Figure 4. Hierarchical clustering analysis of the differentially expressed proteins responding to NSY50 and/or FOC. The fold changes of protein abundance among the four treatments were log2 transformed and delivered to the Cluster and Treeview software. Each row represents individual protein spots and spot numbers, and the protein names are labeled to the right of the corresponding heat maps. Red and green show the higher and lower expression levels, respectively.
Validation of Differentially Expressed Proteins
As shown in Table 2, five differentially expressed protein spots belonging to ENO and eight SAMs protein spots were successfully identified by MALDI-TOF/TOF MS, and they were all significantly regulated by inoculation with NSY50 and/or FOC. This outcome indicated that ENO and SAMs were considered important plant proteins for coping with pathogen invasion. Therefore, ENO and SAMs were analyzed by western blotting to verify the proteomic data. The result showed a similar tendency, in that the expression of ENO and SAMs were slightly more highly regulated by the inoculation of NSY50 but significantly down-regulated by FOC infection compared to the control. However, the expression of both of these proteins was recovered by pre-treatment with NSY50 (NSY50+FOC) compared to FOC, regardless of whether measurements were taken 1, 3, or 9 days post-inoculation with FOC (Figures 5A,C; Supplementary Figures 2–5). Additionally, the level of SAMs was seriously decreased by FOC inoculation on the first day after inoculation (Figure 5C).
Figure 5. Western blot and qRT-PCR analysis of the proteins and genes of ENO and SAMs expression levels for the four treatments in the roots of cucumber seedlings challenged with F. oxysporum at 1, 3, and 9 days post-inoculation (dpi). (A) Western blot analysis of the ENO expression level for the four treatments over 3 days; (B) qRT-PCR analysis the ENO gene expression level for the four treatments over 3 days; (C) Western blot analysis of the SAMs expression level for the four treatments over 3 days; (D) qRT-PCR analysis of the SAMs gene expression level for the four treatments over 3 days. Bars marked with dissimilar letters are significantly different according to Duncan's multiple range test (P < 0.05). Control, seedlings were cultured in normal nutrient solution; NSY50, control+ P. polymyxa NSY50 (2.5 × 106 CFU/mL); FOC, control+ FOC (2.5 × 105 conidia/mL); NSY50+FOC, pre-treated with NSY50 (2.5 × 106 CFU/mL) for 3 days, then challenged with FOC (2.5 × 105 conidia/mL).
Expression Analysis of Several Identified Protein-Related Genes
The ENO and SAMs gene expression levels were then analyzed 1, 3, and 9 days post-inoculation. FOC up-regulated both proteins at all time points compared to the control (Figures 5B,D). The gene expression of ENO after NSY50 treatment was significantly lower compared to the FOC treatment. However, SAMs showed a contrary tendency in the early stages of inoculation (Figure 5D), and the transcript levels in NSY50+FOC were higher and increased more rapidly compared to those of plants that were inoculated with only FOC (FO C) from 1 to 3 days post-inoculation.
Seven other genes related to the identified proteins, including GST, SAMDC, ACS, ACO1, ACO2, HSP70, and OPR1, were selected and subjected to expression pattern analysis. NSY50 up-regulated almost all of these genes (except ACS and ACO1) compared to the control (Figure 6). FOC stress up-regulated the expression of GST, SAMDC and ACO1 but decreased the ACS and ACO2 expression. The combined inoculation (NSY50+FOC) treatments up-regulated expression of almost all of these genes (except SAMDC and ACS) compared with the FOC treatment. We summarize the metabolism pathways that connect the expression of these genes, the proteins related to an antioxidant response, amino acid metabolism, and the EMP-TCA cycle in Figure 7.
Figure 6. Effects of NSY50 and/or FOC treatment on transcript analysis of ENO, GST, SAMs, SAMDC, ACS, ACO1, ACO2, HSP70, and OPR1 in the roots of cucumber seedlings. Each histogram represents a mean ± SE of three independent experiments (n = 3). Different letters indicate significant differences between the treatments (P < 0.05) according to Duncan's multiple range test. Abbreviations: ENO, enolase; GST, glutathione S-transferase; SAMs, S-adenosylmethionine synthase; SAMDC, S-adenosylmethionine decarboxylase; ACS, 1-aminocyclopropane-1- carboxylate synthase; ACO1, 1-aminocyclopropane-1-carboxylate oxidase 1; ACO2, 1-amin ocyclopropane-1-carboxylate oxidase 2; HSP70, heat shock 70 kDa protein; OPR1, oxophytodienoate reductase 1.
Figure 7. Schematic presentation of the effects of PGPB and FOC stress on metabolism in cucumber roots. Changes in protein abundance (marked in red ellipses) and gene expression (blue) were integrated. Red arrows on the life side indicate changes induced by NSY50 compared to the control, green arrows in the middle indicate changes induced by FOC compared to the control, and the purple arrows on the right side indicate changes induced by combined inoculation (NSY50+FOC) conditions. Arrows ↑ or ↓ represent up-regulation or down-regulation, and short lines indicate no change. Abbreviations: BFRUCT, acid beta-fructofuranosidase; FRK, fructokinase; FBA, fructose-bisphosphate aldolase; ENO, enolase; MDH, malate dehydrogenase; ATPase, ATP synthase; SOD, superoxide dismutase; POD, peroxidase; APX, ascorbate peroxidase; GST, glutathione S-transferase; DHAR, glutathione S-transferase DHAR; TRX, thioredoxin; GS, glutamine synthetase; SAMs, S-adenosylmethionine synthase; SAMDC, S-adenosylmethionine decarboxylase; ACS, 1-aminocyclopropane-1-carboxylate synthase; ACO,1-aminocyclopropane-1-carboxylate oxidase; EF2, elongation factor 2; NTF, nuclear transport factor; eIF5A, eukaryotic translation initiation factor 5A; PDI, protein disulfide-isomerase; CPN20, 20 kDa chaperonin; TPA, thiol protease aleurain-like; PRs, proteasome; UEV1C, ubiquitin-conjugating enzyme E2 variant 1C.
Biological control is one of the most beneficial methods because it is both efficient and environmentally friendly for plant protection (Bargabus et al., 2003; Tjamos et al., 2005). Currently, the commercialization of bio-control agents is still limited (Szewczyk et al., 2006) due to complicated mechanisms of bio-control that have yet to be clarified (Niu et al., 2011). Paenibacillus polymyxa has been considered a bio-control agent for a wide range of plant pathogens (Dijksterhuis et al., 1999; Ling et al., 2011; Hong et al., 2016). It is an important Bacillus species because of its ability to produce different bioactive compounds with a broad spectrum of activities or induce systemic resistance against a variety of environmental stresses (Kim et al., 2015; Zhou et al., 2016). With the goal of commercializing P. polymyxa NSY50 in the future, we analyzed the proteins in cucumber roots that responded to NSY50 under FOC stress conditions to investigate the promotion and protection mechanisms of NSY50. The results of our proteome analysis are discussed below.
Proteins Related to Defense
Reactive oxygen species (ROS) have been demonstrated to be involved in symbiotic interactions between plants and microorganisms (Scheler et al., 2013). During pathogen invasion, plants have always suffered from ROS and reactive aldehyde bursts, which often lead to cellular damage (Niu et al., 2011; Gupta et al., 2014; Sengupta et al., 2015; Camejo et al., 2016). Fortunately, plants have developed an antioxidant enzymatic system to scavenge for and collect these toxic compounds. In the present study, a total of 16 spots were identified as defense-related proteins, and most of them were regulated by the inoculation of NSY50 and/or FOC (Table 2, Figure 3C).
Two superoxide dismutases (SOD, spots 25, and 27) were identified and were significantly up-regulated by pre-treatment with NSY50 (NSY50+FOC) compared to FOC treatment. Three peroxidases (POD, spots 19, 24, and 54) with different subunits were identified and showed different accumulation patterns. Regardless of whether they were inoculated with NSY50 or FOC, the expression levels were both up-regulated compared to the control. Additionally, spots 19 and 24 showed expression levels in NSY50+FOC that were higher than the proteins inoculated with FOC alone. Ascorbate peroxidase (APX, spot 21) was also up-regulated by NSY50 (NSY50) compared to the control. It is known that the activation of antioxidant enzymes is a major plant cell method for protection (Youssef et al., 2016). Several studies confirmed the ability of bio-control agents to mitigate the effects of oxidative bursts by increasing the activity of antioxidant enzymes (Israr et al., 2016; Youssef et al., 2016). Thus, the results of the present study indicated that application of NSY50 played an active role in determining the abundance of antioxidant enzymes.
GSTs belong to a family of multifunctional enzymes that play critical roles in protecting tissues from oxidative damage by quenching reactive molecules (An et al., 2016). Previous studies have reported that GSTs can enhance resistance to the virulent bacterial pathogen Pseudomonas syringae pv. tomato DC3000 in Arabidopsis (Jones et al., 2006) and were increased in response to various hormones, such as salicylic acid, ethylene, methyl jasmonate, and auxin, and to biotic and abiotic stress (Lieberherr et al., 2003; Faltin et al., 2010). Additionally, Kwon et al. (2016) found that four glutathione S-transferases were significantly up-regulated by plant-growth-promoting rhizobacterium P. polymyxa E681 in Arabidopsis and suggested that the activation of defense-related proteins may the reason for resistance to the fungal pathogen. It was reported that the overexpression of GSTs induced by another PGPR, P. fluorescens, had an essential role in the ISR by priming rice plants and protecting cells from oxidative damage (Kandasamy et al., 2009). In the present study, two glutathione S-transferases (GST, spots 32, and 37) were identified. One of them (spot 32) was significantly up-regulated in NSY50, and the other (spot 37) showed expression levels in NSY50+FOC that were much higher than those in FOC. The qRT-PCR analysis also showed that the application of NSY50 significantly up-regulated this gene expression, regardless of whether the plants were inoculated with FOC (Figure 6). This result indicated that NSY50 could increase GST content to remove toxins in cucumber seedlings.
Plant thioredoxins appear to play a fundamental role in plant tolerance of oxidative stress by supplying power to reductases detoxifying lipid hydroperoxides or repairing oxidized proteins (Vieira Dos Santos and Rey, 2006). Thioredoxin H-type 1 isoform X2 (TRX, spot 52) belongs to the thioredoxins family and was up-regulated by the inoculation of NSY50 and/or FOC. This result was consistent with the results of Laloi et al. (2004), who reported that the Arabidopsis cytosolic thioredoxin h5 gene was induced by oxidative stress and its response to a pathogen elicitor. In addition, Sun et al. (2010) found that the overexpression of NtTRXh3 (An h-type thioredoxin in tobacco) conferred resistance to Tobacco mosaic virus and Cucumber mosaic virus, both of which had reduced multiplication and pathogenicity in NtTRXh3- overexpressing plants compared to controls. These reductions indicated that TRX participated in defense mechanisms linked to the oxidative bursts resulting from pathogen attack.
In this study, a probable aldo-keto reductase 4 (AKRs, spots 50) was successfully identified, which showed that inoculation with NSY50 significantly up-regulated the expression of aldo-keto reductase (AKRs), with or without FOC. The AKRs superfamily comprises a large number of primarily monomeric protein members that reduce a broad spectrum of substrates, which range from simple sugars to potentially toxic aldehydes (Sengupta et al., 2015). There are a number of interesting reports on the active involvement of plant AKRs in detoxification of stress-induced reactive carbonyls as plant defense systems against biotic stress factors, such as pathogenic attack (Tremblay et al., 2010; Santos et al., 2012; Xu et al., 2013). This may be the first report on a type of protein induced by the PGPR that actively contributes to defense against biotic stressors.
Other defense response proteins, such as heat shock protein 70 (HSP70, spots 4 and 12) and MLP-like protein 328 (MLP328, spot 26), were also identified in this study. HSP70 can be induced by various environmental stresses, such as cold, salt, high temperature, and hypoxic stress, and this protein often acts as a molecular chaperone (Xu et al., 2011; Li et al., 2013b; Yuan et al., 2016) to prevent the irreversible aggregation of other proteins and regulate the folding and accumulation of proteins (Feder and Hofmann, 1999). In this study, it was observed that NSY50 enhances the expression of HSP70 (spot 4) by 3.31-fold compared to the control. Whereas, FOC significantly decreased HSP70 expression both under control conditions and pre-treatment with NSY50 (Table 2), the combined inoculation (NSY50+FOC) showed much higher up-regulation than FOC, although there was no significant difference. The transcription of HSP70 showed a significant increase in NSY50+FOC compared to inoculation with FOC alone (Figure 6). This result indicates that the PGPG are capable of enhancing the expression of this important protein. This result was consistent with the research of Li et al. (2013b), who reported that PGPR P. putida UW4 significantly enhanced the level of plant heat shock proteins under hypoxic stress to improve stress resistance. The role of the MLP-like protein was unclear, and it was hypothesized to be a physiological defense protein that responded to stress (Fukao et al., 2009; An et al., 2016). However, in this study, MLP328 was significantly up-regulated, with a 3.13–3.68-fold increase following the inoculation of NSY50 and/or FOC compared to the control. The mechanism of this type of physiological defense protein needs further study.
Of the differentially accumulated proteins, 14 were identified as responding to FOC and/or NSY50 treatments; these responses were attributed to protein metabolism, and the proteins could be divided into five functional groups (Table 2). The first group includes those involved in protein folding and assembly, including four protein disulfide-isomerase (PDI, spots 14, 42, 44, and 45) and a 20 kDa chaperonin (CPN20, spot 23). Both the NSY50 and FOC treatments enhanced the expression of PDI (spot 14) by 1.90–2.42-fold compared to the control. Several studies have shown that the overexpression of PDI in plants that are challenged with pathogens has been thought to be related to the correct folding of defense proteins (Gruber et al., 2007; Palomares-Rius et al., 2011; Cipriano et al., 2015). However, the remaining three PDI proteins (spots 42, 44, and 45) were significantly decreased by pathogen attack (inoculated with FOC). This effect may be due to the stresses that were inhibited expression of this type of protein, which was consistent with previous reports (Yuan et al., 2016). However, another report showed that the application of PGPR recovered the expression level of this protein under stress conditions (Li et al., 2013b). The CPN20 (spot 23), which was among the most-represented proteins that emphasize the impact of stress on post-translation modification/ modification machinery (Zhang et al., 2015) showed similar results when inoculated with FOC.
The second group included three spots of eukaryotic translation initiation factor 5A (eIF-5A, spots 33, 38, and 39), one eukaryotic translation initiation factor 5A-2 (eIF-5A2, spot 56) and one elongation factor 2 (EF2, spot 47), which are involved in the initiation and elongation stage of mRNA translation and protein synthesis (Thompson et al., 2004; Yuan et al., 2016). Overall, the FOC inoculation significantly decreased the eIF-5A expression levels in cucumber seedling roots, whereas pre-treatment with NSY50 (NSY50+FOC) enhanced the expression of eIF-5A. However, the eIF-5A2 protein had a higher accumulation after FOC inoculation conditions both with and without NSY50 compared to non-stress conditions. The protein eIF-5A2 was reported to likely be involved in the regulation of apoptotic cell death (Feng et al., 2007). EF2 was significantly up-regulated, with a 2.70-fold increase after inoculation with NSY50 compared to the control, whereas inoculation with FOC showed only a slight increase. These results indicated that the PGPR could balance ribosomal particles and protein synthesis, which was consistent with previous reports (Li et al., 2013b).
The third group included proteasome subunit alpha type-2-A (PRs, spot 22) and thiol protease aleurain-like (TPA, spot 41), which are involved in protein degradation. The proteasome was assigned a precise role in the degradation of the oxidized proteins generated by many stresses, which avoided the inhibition of the plant's metabolic pathway (Polge et al., 2009). In this study, infection with FOC down-regulated the expression level of TPA, and the plants metabolic pathway may have been damaged by this destructive fungal pathogen, which was consistent with the TEM result. The PRs significantly decreased under FOC stress but were remarkably up-regulated by inoculation with NSY50, which indicated that the degradation of proteins with oxidative damage was enhanced (An et al., 2016). The fourth group was protein transport and included one protein called nuclear transport factor 2-like (NTF2, spot 28). The NTF2 facilitated protein transport into the nucleus (Zhao et al., 2006). Both inoculation with NSY50 and FOC were significantly up-regulated the expression of this protein. The ubiquitin-conjugating enzyme E2 variant 1C (UEV1C, spot 30) belonged to the fifth group, protein modification, and it participates in protein modification by catalyzing the covalent attachment of ubiquitin to proteins (Yuan et al., 2016). In this study, FOC significantly decreased the expression of UEV1C, whereas pre-treatment with NSY50 (NSY50+FOC) enhanced its expression. These results suggest that the application of NSY50 alleviated the damage to protein metabolism that was induced by F. oxysporum.
Carbohydrate and Energy Metabolism
The proteins related to carbohydrate and energy metabolism, particularly the proteins involved in the glycolysis (EMP) and tricarboxylic acid (TCA) cycle were significantly changed and are shown in Table 2. A total of 11 identified proteins in response to two strains were found and were divided into four functional groups. The first group focused on glycolysis and included five spots with enolase isoform (ENO, spots 1, 2, 17, 18, and 46), protein DJ-1 homolog D (DJ1D, spot 13), probable fructokinase-4 (FRK, spot 20) and fructose-bisphosphate aldolase (FBA, spot 55). The others were involved with the tricarboxylic acid cycle, including malate dehydrogenase (MDH, spot 36), sucrose metabolism protein acid beta-fructofuranosidase-like (BFRUCT, spot 43), and energy metabolism protein ATP synthase (ATPase, spot 3). Almost all of them were down-regulated by FOC and recovered by treatment with NSY50 (NSY50+FOC), except FBA and BFRUCT.
Enolases, known as 2-phospho-D-glycerate hydrolases, are one of the most important enzymes in glycolysis. They catalyze the dehydration of 2-phosphoglycerate into phosphoenolpyruvate. Five spots were identified as enolases in this study, and western blotting was used to verify proteomic data (Table 2, Figure 5A; Supplementary Figures 2, 3) and showed a similar tendency, regardless of whether measurements were taken 1, 3, or 9 days post-inoculation with the pathogen FOC. The gene expression showed a reverse trend, and inoculation with FOC significantly up-regulated its gene expression (Figure 5B). Interestingly, protein abundance significantly recovered along with inoculation time (Figure 5A). Several studies reported that both the ENO protein and the gene expression increased under salt stress and attempted to generate more energy to cope with the stress (Ndimba et al., 2005; Zhong et al., 2015; An et al., 2016; Yuan et al., 2016).
The DJ-1D protein (DJ1D, spot 13), which was reported to have several functions, including functioning as an antioxidant (Taira et al., 2004), is a redox-dependent chaperone (Krebiehl et al., 2010) and can convert glyoxals to carboxylic acids (Kwon et al., 2013). Glyoxals are reactive 2-oxoaldehydes that are found in cells under various stress conditions (Thornalley, 1996). In this study, DJ1D was significantly up-regulated, with 4.25- and 2.28-fold increases after inoculation with NSY50 and/or FOC, respectively.
Li et al. (2013b) reported that the application of PGPR P. putida UW4 adjusted the EMP-TCA metabolism under hypoxic stress. In contrast, Kwon et al. (2016) found that some of the proteins related to the glycolysis pathway, such as FBA and pyruvate dehydrogenase complex component E2, were significantly up-regulated by the PGPR strain P. polymyxa E681 in Arabidopsis roots. Carbohydrate metabolism provides the plant with required carbon, which is critical for the production of new tissues, and also can have profound effects on plant growth through modulation of cell division and expansion (Eveland and Jackson, 2012; Kwon et al., 2016). In this study, we found that the application of P. polymyxa NSY50 increased the abundance of FRK, ENO, BFRUCT, and ATPase compared to control conditions. Additionally, NSY50 partially recovered the ENO, FRK, DJ1D, MDH, and ATPase expression, which had been decreased by FOC. These results suggested that the application of NSY50 could promote EMP-TCA pathway activity. Presumably, these results imply there was an enhancement in energy generation and carbohydrate production to promote plant growth and resist stress damage.
Amino Acid Metabolism
Ten spots were identified as proteins involved in amino acid metabolism. Eight S-adenosylmethionine synthases (SAMs, spots 6, 7, 8, 9, 10, 15, 16, and 48) comprising different subunits were identified and markedly decreased (except spot 16) in response to FOC infection. SAMs catalyze the biosynthesis of S-adenosylmethionine (SAM) from methionine and ATP, and they are active in the biosynthesis of lignin, glycine betaine and polyamines (Fontecave et al., 2004). These proteins can be induced by various environmental stresses, such as cold, salt, drought stress and even pathogen infection (Li et al., 2013a; Huang et al., 2016), and most of them were down-regulated by stress (Wang et al., 2016; Yuan et al., 2016). This outcome was consistent with the results showing that the inoculation of FOC significantly decreased SAMs expression (Table 2, Figure 5C; Supplementary Figures 4, 5). However, pre-treatment with NSY50 up-regulated SAMs abundance both under control and pathogen infection conditions. The accumulation of SAMs in plants has been shown to be related to the enhance tolerance to various kinds of environmental stresses. Overexpression of suadea salsa SAMs gene promoted salt tolerance in transgenic tobacco (Qi et al., 2010). An increase in production of SAMs was observed in cold stressed rice (Cui et al., 2005), mechanically wounded Phaseolus lunatus (Arimura et al., 2002), salt-stressed barley (Witzel et al., 2009), and in response to cotton worm feeding in soybean (Fan et al., 2012). Hence, the remarkable expression of SAMs by NSY50 suggested that plant might promote SAM-dependant metabolism through synthesize more SAM to further enhance the synthesis of new proteins and stress tolerance.
Glutamine synthetase (GS) is a key enzyme in plant N assimilation that catalyzes the combination of ammonia and glutamate into glutamine (Nam et al., 2011; Rogić et al., 2015). Two types of GS (spots 11 and 35) were identified and showed different accumulation patterns. However, they were both up-regulated by NSY50, regardless of whether the plants were inoculated with FOC. This result suggested that the application of NSY50 could lead to a high accumulation of amino acid metabolism, and this was probably one way in which NSY50 promoted plant growth and intensified resistance to pathogen infection for cucumber seedlings.
Fatty Acid Metabolism Proteins
Jasmonic acid is a signaling compound that influences multiple cellular functions and plays crucial roles in the signaling network regulating the development of induced resistance, including systemic acquired resistance (SAR) and induced systemic resistance (Glazebrook, 2001). Niu et al. (2011) found that PGPR Bacillus cereus AR156 induced ISR to Pseudomonas syringae pv. tomato DC3000 in Arabidopsis via the simultaneous activation of salicylic acid-, jasmonic acid-, and ethylene-dependent signaling pathways. In the present study, 12-oxophytodienoate reductase 1 (spot 49), which is involved in jasmonate biosynthesis, was identified to be induced by both NSY50 and FOC inoculation. Additionally, NSY50 treatment showed a 2.02-fold increase in abundance compared to FOC. A good correlation was observed between this protein and transcript abundance, as etermined by qPCR analysis (Figure 6). This outcome indicated that P. polymyxa NSY50 may stimulate the jasmonate signal transduction network.
Enoyl-[acyl-carrier-protein] reductase [NADH] (ACP, spot 51), which is a subunit of the fatty acid synthase complex that catalyzes the de novo synthesis of fatty acids (Mou et al., 2000), was significantly up-regulated by inoculation with FOC. It was recently reported that ACP was sensitive to the diphenyl ether cyperin produced by certain pathogenic fungi in the tissues of infected plants, which likely contributed to the virulence of these disease agents (Dayan et al., 2008). The role of the ACP during pathogen infection remains unclear, and further study is required.
Secondary Metabolism and Cell-Related Proteins
The NSY50 was also found to regulate proteins involved in other metabolic pathways, such as secondary metabolism (spot 34) and cell-related protein (spots 5 and 29), which indicated that NSY50 regulated various metabolic pathways to counteract FOC infection. The mechanism of action will need to be studied in future investigations.
To the best of our knowledge, this is the first proteomic-based research to focus on the interactions between cucumber plants, the pathogen F. oxysporum, and PGPR strain P. polymyxa NSY50. A total of 56 protein spots were identified with a change more than 1.5-fold or less than the 0.67-fold in the protein abundance ratio in cucumber seedling roots after inoculation with P. polymyxa NSY50 and/or FOC. The majority of these proteins were related to defense responses, protein metabolism, carbohydrate and energy metabolism and amino acid metabolism. These proteins might work cooperatively to enhance resistance to FOC attack and keep plant growth on a mostly even keel. Our results showed that the improved plant growth and defenses by P. polymyxa NSY50 might be associated with the following processes: (i) activation of stress defense-related proteins by stress injuries that are alleviated by NSY50; (ii) NSY50 stimulating protein synthesis and degrading damaged proteins induced by FOC; (iii) up-regulation of proteins in the EMP-TCA route to provide more carbohydrates and energy to promote plant growth and resist debilitation; or (iv) adjusting the plant amino acid metabolism to increase biomass of the plant and regulate other metabolic pathways, such as jasmonic acid, to enhance resistance to FOC. This research enriched our understanding of the mechanisms of how PGPR mediates the stress response in plants and promotes protection from pathogen infection.
We thank the numerous individuals who participated in this research. SG designed the study and guided the research. ND and LS wrote the main manuscript text and performed the experiments. YY prepared all the figures and performed some of the experiments. BL contributed new reagents and analytical tools. SS and JS edited this manuscript. All authors reviewed and approved the manuscript.
This work was financially supported by the National Natural Science Foundation of China (No. 31672199, No. 31471869, and No. 31401919), the China Agriculture Research System (CARS-25-C-03), and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.01859/full#supplementary-material
Ahn, I. P., Chung, H. S., and Lee, Y. H. (1998). Vegetative compatibility groups and pathogenicity among isolates of Fusarium oxysporum f. sp. cucumerinum. Plant Dis. 82, 244–246. doi: 10.1094/PDIS.1918.104.22.168
An, Y., Zhou, H., Zhong, M., Sun, J., Shu, S., Shao, Q., et al. (2016). Root proteomics reveals cucumber 24-epibrassinolide responses under Ca(NO3)2 stress. Plant Cell Rep. 35, 1081–1101. doi: 10.1007/s00299-016-1940-z
Arimura, G., Ozawa, R., Nishioka, T., Boland, W., Koch, T., Kühnemann, F., et al. (2002). Herbivore-induced volatiles induce the emission of ethylene in neighboring lima bean plants. Plant J. 29, 87–98. doi: 10.1046/j.1365-313x.2002.01198.x
Bargabus, R. L., Zidack, N. K., and Sherwood Jacobsen, B. J. (2003). Oxidative burst elicited by Bacillus mycoides isolate Bac, J., a biological control agent, occurs independently of hypersensitive cell death in sugar beet. Mol. Plant Microbe Interact. 16, 1145–1153. doi: 10.1094/MPMI.2003.16.12.1145
Bottini, R., Cassán, F., and Piccoli, P. (2004). Gibberellin production by bacteria and its involvement in plant growth promotion and yield increase. Appl. Microbiol. Biotechnol. 65, 497–503. doi: 10.1007/s00253-004-1696-1
Bradford, M. M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248–254. doi: 10.1016/0003-2697(76)90527-3
Camejo, D., Guzmán-Cedeño, Á., and Moreno, A. (2016). Reactive oxygen species, essential molecules, during plant-pathogen interactions. Plant Physiol. Biochem. 103, 10–23. doi: 10.1016/j.plaphy.2016.02.035
Cao, Y., Xu, Z., Ling, N., Yuan, Y., Yang, X., Chen, L., et al. (2012). Isolation and identification of lipopeptides produced by B. subtilis SQR 9 for suppressing Fusarium wilt of cucumber. Sci. Hortic. 135, 32–39. doi: 10.1016/j.scienta.2011.12.002
Carrillo, E., Rubiales, D., and Castillejo, M. (2014). Proteomic analysis of pea (Pisum sativum L.) response during compatible and incompatible interactions with the pea aphid (Acyrthosiphon pisum H.). Plant Mol. Biol. Rep. 32, 697–718. doi: 10.1007/s11105-013-0677-x
Chen, D., Liu, X., Li, C., Tian, W., Shen, Q., and Shen, B. (2014). Isolation of Bacillus amyloliquefaciens S20 and its application in control of eggplant bacterial wilt. J. Environ. Manag. 137, 120–127. doi: 10.1016/j.jenvman.2014.01.043
Chung, S., Kong, H., Buyer, J. S., Lakshman, D. K., Lydon, J., Kim, S. D., et al. (2008). Isolation and partial characterization of Bacillus subtilis ME488 for suppression of soilborne pathogens of cucumber and pepper. Appl. Microbiol. Biotechnol. 80, 115–123. doi: 10.1007/s00253-008-1520-4
Cipriano, A. K., Gondim, D. M., Vasconcelos, I. M., Martins, J. A., Moura, A. A., Moreno, F. B., et al. (2015). Proteomic analysis of responsive stem proteins of resistant and susceptible cashew plants after Lasiodiplodia theobromae infection. J. Proteomics 113, 90–109. doi: 10.1016/j.jprot.2014.09.022
Dayan, F. E., Ferreira, D., Wang, Y. H., Khan, I. A., McInroy, J. A., and Pan, Z. (2008). A pathogenic fungi diphenyl ether phytotoxin targets plant enoyl (acyl carrier protein) reductase. Plant Physiol 147, 1062–1071. doi: 10.1104/pp.108.118372
Dijksterhuis, J., Sanders, M., Gorris, L. G., and Smid, E. J. (1999). Antibiosis plays a role in the context of direct interaction during antagonism of Paenibacillus polymyxa towards Fusarium oxysporum. J. Appl. Microbiol. 86, 13–21. doi: 10.1046/j.1365-2672.1999.t01-1-00600.x
Du, N., Shi, L., Du, L., Yuan, Y., Li, B., Sang, T., et al. (2015). Effect of vinegar residue compost amendments on cucumber growth and Fusarium wilt. Environ. Sci. Pollut. Res. Int. 22, 19133–19141. doi: 10.1007/s11356-015-4816-9
Faltin, Z., Holland, D., Velcheva, M., Tsapovetsky, M., RoeckelDrevet, P., Handa, A. K., et al. (2010). Glutathione peroxidase regulation of reactive oxygen species level is crucial for in vitro plant differentiation. Plant Cell Physiol. 51, 1151–1162. doi: 10.1093/pcp/pcq082
Fan, R., Wang, H., Wang, Y., and Yu, D. (2012). Proteomic analysis of soybean defense response induced by cotton worm (prodenia litura, fabricius) feeding. Proteome Sci. 10, 338–348. doi: 10.1186/1477-5956-10-16
Feng, H., Chen, Q., Feng, J., Zhang, J., Yang, X., and Zuo, J. (2007). Functional characterization of the Arabidopsis eukaryotic translation initiation factor 5A-2 that plays a crucial role in plant growth and development by regulating cell division, cell growth, and cell death. Plant Physiol. 144, 1531–1545. doi: 10.1104/pp.107.098079
Fukao, Y., Ferjani, A., Fujiwara, M., Nishimori, Y., and Ohtsu, I. (2009). Identification of zinc-responsive proteins in the roots of Arabidopsis thaliana using a highly improved method of two-dimensional electrophoresis. Plant Cell Physiol. 50, 2234–2239. doi: 10.1093/pcp/pcp154
Gruber, C. W., Cemazar, M., Clark, R. J., Horibe, T., Renda, R. F., Anderson, M. A., et al. (2007). A novel plant protein-disulfide isomerase involved in the oxidative folding of cystine knot defense proteins. J. Biol. Chem. 282, 20435–20446. doi: 10.1074/jbc.M700018200
Gupta, K. J., Mur, L. A., and Brotman, Y. (2014). Trichoderma asperelloides suppresses nitric oxide generation elicited by Fusarium oxysporum in Arabidopsis roots. Mol. Plant Microbe Interact. 27, 307–314. doi: 10.1094/MPMI-06-13-0160-R
He, L., Li, B., Lu, X., Yuan, L., Yang, Y., Yuan, Y., et al. (2015). The effect of exogenous calcium on mitochondria, respiratory metabolism enzymes and ion transport in cucumber roots under hypoxia. Sci. Rep. 5:11391. doi: 10.1038/srep11391
He, L., Lu, X., Tian, J., Yang, Y., Li, B., Li, J., et al. (2012). Proteomic analysis of the effects of exogenous calcium on hypoxic-responsive proteins in cucumber roots. Proteome Sci. 10, 1–15. doi: 10.1186/1477-5956-10-42
Hong, C. E., Kwon, S. Y., and Park, J. M. (2016). Biocontrol activity of Paenibacillus polymyxa AC-1 against Pseudomonas syringae and its interaction with Arabidopsis thaliana. Microbiol. Res. 185, 13–21. doi: 10.1016/j.micres.2016.01.004
Huang, X., Shi, D., Sun, F., Lu, H., Liu, J., and Wu, W. (2012). Efficacy of sludge and manure compost amendments against Fusarium wilt of cucumber. Environ. Sci. Pollut. Res. Int. 19, 3895–3905. doi: 10.1007/s11356-012-1025-7
Huang, Y., Ma, H. Y., Huang, W., Wang, F., Xu, Z. S., and Xiong, A. S. (2016). Comparative proteomic analysis provides novel insight into the interaction between resistant vs susceptible tomato cultivars and TYLCV infection. BMC Plant Biol. 16:162. doi: 10.1186/s12870-016-0819-z
Idriss, E. E., Makarewicz, O., Farouk, A., Rosner, K., Greiner, R., Bochow, H., et al. (2002). Extracellular phytase activity of Bacillus amyloliquefaciens FZB45 contributes to its plant-growth-promoting effect. Microbiology 148(Pt 7), 2097–2109. doi: 10.1099/00221287-148-7-2097
Israr, D., Mustafa, G., Khan, K. S., Shahzad, M., Ahmad, N., and Masood, S. (2016). Interactive effects of phosphorus and Pseudomonas putida on chickpea (Cicer arietinum L.) growth, nutrient uptake, antioxidant enzymes and organic acids exudation. Plant Physiol. Biochem. 108, 304–312. doi: 10.1016/j.plaphy.2016.07.023
Jiang, C. H., Wu, F., Yu, Z. Y., Xie, P., Ke, H. J., Li, H. W., et al. (2015). Study on screening and antagonistic mechanisms of Bacillus amyloliquefaciens 54 against bacterial fruit blotch (BFB) caused by Acidovorax avenae subsp. citrulli. Microbiol. Res. 170, 95–104. doi: 10.1016/j.micres.2014.08.009
Jones, A. M. E., Thomas, V., Bennett, M. H., Mansfield, J., and Grant, M. (2006). Modifications to the Arabidopsis defense proteome occur prior to significant transcriptional change in response to inoculation with Pseudomonas syringae. Plant Physiol. 142, 1603–1620. doi: 10.1104/pp.106.086231
Kandasamy, S., Loganathan, K., Muthuraj, R., Duraisamy, S., Seetharaman, S., Thiruvengadam, R., et al. (2009). Understanding the molecular basis of plant growth promotional effect of Pseudomonas fluorescens on rice through protein profiling. Proteome Sci. 7, 9143–9148. doi: 10.1186/1477-5956-7-47
Kim, Y. S., Kotnala, B., Kim, Y. H., and Jeon, Y. (2015). Biological characteristics of Paenibacillus polymyxa GBR-1 involved in root rot of stored Korean ginseng. J. Ginseng Res. 40, 453–461. doi: 10.1016/j.jgr.2015.09.003
Kloepper, J., and Metting, F. Jr. (1992). “Plant growth-promoting rhizobacteria as biological control agents,” in Soil Microbial Ecology: Applications in Agricultural and Environmental Management, ed F. B. Metting (New York, NY: Marcel Dekker), 255–274.
Kochar, M., Upadhyay, A., and Srivastava, S. (2011). Indole-3-acetic acid biosynthesis in the biocontrol strain Pseudomonas fluorescens Psd and plant growth regulation by hormone overexpression. Res. Microbiol. 162, 426–435. doi: 10.1016/j.resmic.2011.03.006
Konishi, H., Ishiguro, K., and Komatsu, S. (2001). A proteomics approach towards understanding blast fungus infection of rice grown under different levels of nitrogen fertilization. Proteomics 1, 1162–1171. doi: 10.1002/1615-9861(200109)1:9<1162::AID-PROT1162>3.0.CO;2-S
Krebiehl, G., Ruckerbauer, S., Burbulla, L. F., Kieper, N., Maurer, B., Waak, J., et al. (2010). Reduced basal autophagy and impaired mitochondrial dynamics due to loss of Parkinson's disease-associated protein DJ-1. PLoS ONE 5:e9367. doi: 10.1371/journal.pone.0009367
Kwon, Y. S., Lee, D. Y., Rakwal, R., Baek, S.-B., Lee, J. H., Kwak, Y.-S., et al. (2016). Proteomic analyses of the interaction between the plant-growth promoting rhizobacterium Paenibacillus polymyxa E681 and Arabidopsis thaliana. Proteomics 16, 122–135. doi: 10.1002/pmic.201500196
Laloi, C., Mestresortega, D., Marco, Y., Meyer, Y., and Reichheld, J. P. (2004). The Arabidopsis cytosolic thioredoxin h5 gene induction by oxidative stress and its W-box-mediated response to pathogen elicitor. Plant Physiol. 134, 1006–1016. doi: 10.1104/pp.103.035782
Li, B., He, L., Guo, S., Li, J., Yang, Y., Yan, B., et al. (2013a). Proteomics reveal cucumber Spd-responses under normal condition and salt stress. Plant Physiol. Biochem. 67, 7–14. doi: 10.1016/j.plaphy.2013.02.016
Li, J., McConkey, B. J., Cheng, Z., Guo, S., and Glick, B. R. (2013b). Identification of plant growth-promoting bacteria-responsive proteins in cucumber roots under hypoxic stress using a proteomic approach. J. Proteomics 84, 119–131. doi: 10.1016/j.jprot.2013.03.011
Li, L., Ma, J., Li, Y., Wang, Z., Gao, T., and Wang, Q. (2012). Screening and partial characterization of Bacillus with potential applications in biocontrol of cucumber Fusarium wilt. Crop Protect. 35, 29–35. doi: 10.1016/j.cropro.2011.12.004
Lieberherr, D., Wagner, U., Dubuis, P. H., Métraux, J. P., and Mauch, F. (2003). The rapid induction of glutathione S-transferases AtGSTF2 and AtGSTF6 by avirulent Pseudomonas syringae is the result of combined salicylic acid and ethylene signaling. Plant Cell Physiol. 44, 750–757. doi: 10.1093/pcp/pcg093
Lin, Y., Du, D., Si, C., Zhao, Q., Li, Z., and Li, P. (2014). Potential biocontrol Bacillus sp. strains isolated by an improved method from vinegar waste compost exhibit antibiosis against fungal pathogens and promote growth of cucumbers. Biol. Control 71, 7–15. doi: 10.1016/j.biocontrol.2013.12.010
Ling, N., Huang, Q., Guo, S., and Shen, Q. (2011). Paenibacillus polymyxa SQR-21 systemically affects root exudates of watermelon to decrease the conidial germination of Fusarium oxysporum f. sp. niveum. Plant Soil 341, 485–493. doi: 10.1007/s11104-010-0660-3
Martínez-Viveros, O., Jorquera, M. A., Crowley, D. E., Gajardo, G., and Mora, M. L. (2010). Mechanisms and practical considerations involved in plant growth promotion by rhizobacteria. J. Soil Sci. Plant Nutr. 10, 293–319. doi: 10.4067/s0718-95162010000100006
Mou, Z., He, Y., Dai, Y., Liu, X., and Li, J. (2000). Deficiency in fatty acid synthase leads to premature cell death and dramatic alterations in plant morphology. Plant Cell 12, 405–418. doi: 10.1105/tpc.12.3.405
Nam, M. H., Sun, M. H., Kim, K. M., Park, W. J., Seo, J. B., Cho, K., et al. (2011). Comparative proteomic analysis of early salt stress-responsive proteins in roots of SnRK2 transgenic rice. Proteome Sci. 10, 1–19. doi: 10.1186/1477-5956-10-25
Ndimba, B. K., Chivasa, S., Simon, W. J., and Slabas, A. R. (2005). Identification of Arabidopsis salt and osmotic stress responsive proteins using two-dimensional difference gel electrophoresis and mass spectrometry. Proteomics 5, 4185–4196. doi: 10.1002/pmic.200401282
Niu, D. D., Liu, H. X., Jiang, C. H., Wang, Y. P., Wang, Q. Y., Jin, H. L., et al. (2011). The plant growth-promoting rhizobacterium Bacillus cereus AR156 induces systemic resistance in Arabidopsis thaliana by simultaneously activating salicylate- and jasmonate/ethylene-dependent signaling pathways. Mol. Plant Microbe Interact. 24, 533–542. doi: 10.1094/MPMI-09-10-0213
Ongena, M., Jourdan, E., Adam, A., Paquot, M., Brans, A., Joris, B., et al. (2007). Surfactin and fengycin lipopeptides of Bacillus subtilis as elicitors of induced systemic resistance in plants. Environ. Microbiol. 9, 1084–1090. doi: 10.1111/j.1462-2920.2006.01202.x
Palacios, O. A., Bashan, Y., and De-Bashan, L. E. (2014). Proven and potential involvement of vitamins in interactions of plants with plant growth-promoting bacteria—an overview. Biol. Fertil. Soils 50, 415–432. doi: 10.1007/s00374-013-0894-3
Palomares-Rius, J. E., Castillo, P., Navas-Cortés, J. A., Jiménez-Díaz, R. M., and Tena, M. (2011). A proteomic study of in-root interactions between chickpea pathogens: the root-knot nematode Meloidogyne artiellia and the soil-borne fungus Fusarium oxysporum f. sp. ciceris race 5. J. Proteomics 74, 2034–2051. doi: 10.1016/j.jprot.2011.05.026
Polge, C., Jaquinod, M., Holzer, F., Bourguignon, J., Walling, L., and Brouquisse, R. (2009). Evidence for the existence in Arabidopsis thaliana of the proteasome proteolytic pathway: activation in response to cadmium. J. Biol. Chem. 284, 35412–35424. doi: 10.1074/jbc.M109.035394
Qi, Y. C., Wang, F. F., Hui, Z., and Liu, W. Q. (2010). Overexpression of suadea salsa S-adenosylmethionine synthetase gene promotes salt tolerance in transgenic tobacco. Acta Physiol. Plant. 32, 263–269. doi: 10.1007/s11738-009-0403-3
Rogić, T., Horvatić, A., Tkalec, M., Cindrić, M., and Balen, B. (2015). Proteomic analysis of Mammillaria gracilis Pfeiff. in vitro-grown cultures exposed to iso-osmotic NaCl and mannitol. Plant Cell Tissue Organ Cult. 122, 127–146. doi: 10.1007/s11240-015-0756-9
Santos, C. S., Pinheiro, M., Silva, A. I., Egas, C., and Vasconcelos, M. W. (2012). Searching for resistance genes to Bursaphelenchus xylophilus using high throughput screening. BMC Genomics 13, 1–16. doi: 10.1186/1471-2164-13-599
Sengupta, D., Naik, D., and Reddy, A. R. (2015). Plant aldo-keto reductases (AKRs) as multi-tasking soldiers involved in diverse plant metabolic processes and stress defense: a structure-function update. J. Plant Physiol. 179, 40–55. doi: 10.1016/j.jplph.2015.03.004
Shi, L., Du, N., Yuan, Y., Shu, S., Sun, J., and Guo, S. (2016). Vinegar residue compost as a growth substrate enhances cucumber resistance against the Fusarium wilt pathogen Fusarium oxysporum by regulating physiological and biochemical responses. Environ. Sci. Pollut. Res. Int. 23, 18277–18287. doi: 10.1007/s11356-016-6798-7
Spaepen, S., Bossuyt, S., Engelen, K., Marchal, K., and Vanderleyden, J. (2014). Phenotypical and molecular responses of Arabidopsis thaliana roots as a result of inoculation with the auxin-producing bacterium Azospirillum brasilense. New Phytol. 201, 850–861. doi: 10.1111/nph.12590
Sun, L., Ren, H., Liu, R., Li, B., Wu, T., Sun, F., et al. (2010). An h-type thioredoxin functions in tobacco defense responses to two species of viruses and an abiotic oxidative stress. Mol. Plant Microbe Interact. 23, 1470–1485. doi: 10.1094/MPMI-01-10-0029
Szewczyk, B., Hoyos-Carvajal, L., Paluszek, M., Skrzecz, I., and Lobo de Souza, M. (2006). Baculoviruses– re-emerging biopesticides. Biotechnol. Adv. 24, 143–160. doi: 10.1016/j.biotechadv.2005.09.001
Taira, T., Saito, Y., Niki, T., Iguchi-Ariga, S. M., Takahashi, K., and Ariga, H. (2004). DJ-1 has a role in antioxidative stress to prevent cell death. EMBO Rep. 5, 213–218. doi: 10.1038/sj.embor.7400074
Tamehiro, N., Okamoto-Hosoya, Y., Okamoto, S., Ubukata, M., Hamada, M., Naganawa, H., et al. (2002). Bacilysocin, a novel phospholipid antibiotic produced by Bacillus subtilis 168. Antimicrob. Agents Chemother. 46, 315–320. doi: 10.1128/AAC.46.2.315-320.2002
Thompson, J. E., Hopkins, M. T., Taylor, C., and Wang, T. W. (2004). Regulation of senescence by eukaryotic translation initiation factor 5A: implications for plant growth and development. Trends Plant Sci. 9, 174–179. doi: 10.1016/j.tplants.2004.02.008
Thornalley, P. J. (1996). Pharmacology of methylglyoxal: formation, modification of proteins and nucleic acids, and enzymatic detoxification-A role in pathogenesis and antiproliferative chemotherapy. Gen. Pharmacol. 27, 565–573. doi: 10.1016/0306-3623(95)02054-3
Tjamos, S. E., Flemetakis, E., Paplomatas, E. J., and Katinakis, P. (2005). Induction of resistance to Verticillium dahliae in Arabidopsis thaliana by the biocontrol agent K-165 and pathogenesis-related proteins gene expression. Mol. Plant Microbe. Interact. 18, 555–561. doi: 10.1094/MPMI-18-0555
Tremblay, A., Hosseini, P., Alkharouf, N. W., Li, S., and Matthews, B. F. (2010). Transcriptome analysis of a compatible response by Glycine max to Phakopsora pachyrhizi infection. Plant Sci. 179, 183–193. doi: 10.1016/j.plantsci.2010.04.011
Wang, F. X., Ma, Y. P., Yang, C. L., Zhao, P. M., Yuan, Y., Jian, G. L., et al. (2011). Proteomic analysis of the sea-island cotton roots infected by wilt pathogen Verticillium dahliae. Proteomics 11, 4296–4309. doi: 10.1002/pmic.201100062
Wang, W., Chen, L. N., Wu, H., Zang, H., Gao, S., Yang, Y., et al. (2013). Comparative proteomic analysis of rice seedlings in response to inoculation with Bacillus cereus. Lett. Appl. Microbiol. 56, 208–215. doi: 10.1111/lam.12035
Wang, X., Oh, M. W., and Komatsu, S. (2016). Characterization of S-adenosylmethionine synthetases in soybean under flooding and drought stresses. Biol. Plant 60, 269–278. doi: 10.1007/s10535-016-0586-6
Witzel, K., Weidner, A., Surabhi, G. K., Börner, A., and Mock, H. P. (2009). Salt stress-induced alterations in the root proteome of barley genotypes with contrasting response towards salinity. J. Exp. Bot. 60, 3545–3557. doi: 10.1093/jxb/erp198
Xu, L., Liu, Z. Y., Zhang, K., Lu, Q., Liang, J., and Zhang, X. Y. (2013). Characterization of the Pinus massoniana transcriptional response to Bursaphelenchus xylophilus infection using suppression subtractive hybridization. Int. J. Mol. Sci. 14, 11356–11375. doi: 10.3390/ijms140611356
Xu, Q., Zou, Q., Zheng, H., Zhang, F., Tang, B., and Wang, S. (2011). Three heat shock proteins from Spodoptera exigua: gene cloning, characterization and comparative stress response during heat and cold shocks. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 159, 92–102. doi: 10.1016/j.cbpb.2011.02.005
Xu, Z., Zhang, R., Wang, D., Qiu, M., Feng, H., Zhang, N., et al. (2014). Enhanced control of cucumber wilt disease by Bacillus amyloliquefaciens SQR9 by altering the regulation of its degu phosphorylation. Appl. Environ. Microbiol. 80, 2941–2950. doi: 10.1128/AEM.03943-13
Yang, X., Chen, L., Yong, X., and Shen, Q. (2010). Formulations can affect rhizosphere colonization and biocontrol efficiency of Trichoderma harzianum SQR-T037 against Fusarium wilt of cucumbers. Biol. Fertil. Soils 47, 239–248. doi: 10.1007/s00374-010-0527-z
Youssef, S. A., Tartoura, K. A., and Abdelraouf, G. A. (2016). Evaluation of Trichoderma harzianum and Serratia proteamaculans effect on disease suppression, stimulation of ROS-scavenging enzymes and improving tomato growth infected by Rhizoctonia solani. Biol. Control. 100, 79–86. doi: 10.1016/j.biocontrol.2016.06.001
Yuan, Y., Zhong, M., Shu, S., Du, N., Sun, J., and Guo, S. (2016). Proteomic and physiological analyses reveal putrescine responses in roots of cucumber stressed by NaCl. Front. Plant Sci. 7:1035. doi: 10.3389/fpls.2016.01035
Zamioudis, C., Korteland, J., Van Pelt, J. A., van Hamersveld, M., Dombrowski, N., Bai, Y., et al. (2015). Rhizobacterial volatiles and photosynthesis-related signals coordinate MYB72 expression in Arabidopsis roots during onset of induced systemic resistance and iron-deficiency responses. Plant J. 84, 309–322. doi: 10.1111/tpj.12995
Zhang, D., Meng, K. X., Hao, Y. H., Fan, H. Y., Cui, N., Wang, S. S., et al. (2016). Comparative proteomic analysis of cucumber roots infected by Fusarium oxysporum f. sp. Cucumerium Owen. Physiol. Mol. Plant Pathol. 96, 77–84. doi: 10.1016/j.pmpp.2016.09.002
Zhao, Q., Leung, S., Corbett, A. H., and Meier, I. (2006). Identification and characterization of the Arabidopsis orthologs of nuclear transport factor 2, the nuclear import factor of ran. Plant Physiol. 140, 869–878. doi: 10.1104/pp.105.075499
Zhong, M., Yuan, Y., Shu, S., Sun, J., Guo, S., Yuan, R., et al. (2015). Effects of exogenous putrescine on glycolysis and Krebs cycle metabolism in cucumber leaves subjected to salt stress. Plant Growth Regul. 79, 319–330. doi: 10.1007/s10725-015-0136-9
Zhou, C., Guo, J., Zhu, L., Xiao, X., Xie, Y., Zhu, J., et al. (2016). Paenibacillus polymyxa BFKC01 enhances plant iron absorption via improved root systems and activated iron acquisition mechanisms. Plant Physiol. Biochem. 105, 162–173. doi: 10.1016/j.plaphy.2016.04.025
Keywords: cucumber, Fusarium oxysporum f. sp. cucumerinum, plant-growth-promoting rhizobacteria, proteomics, root
Citation: Du N, Shi L, Yuan Y, Li B, Shu S, Sun J and Guo S (2016) Proteomic Analysis Reveals the Positive Roles of the Plant-Growth-Promoting Rhizobacterium NSY50 in the Response of Cucumber Roots to Fusarium oxysporum f. sp. cucumerinum Inoculation. Front. Plant Sci. 7:1859. doi: 10.3389/fpls.2016.01859
Received: 26 September 2016; Accepted: 25 November 2016;
Published: 14 December 2016.
Edited by:Jie Zhou, Zhejiang University, China
Reviewed by:Yuan Huang, Huazhong Agricultural University, China
Abu Hena Mostafa Kamal, University of Texas at Arlington, USA
Copyright © 2016 Du, Shi, Yuan, Li, Shu, Sun and Guo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shirong Guo, firstname.lastname@example.org
†These authors have contributed equally to this work.