Cinnamaldehyde inhibits the growth of Phytophthora capsici through disturbing metabolic homoeostasis

Strain and reagents

P. capsici susceptible to conventional treatments was isolated and purified from a stem base of the diseased pepper in Henan Province of China. Single-spore isolate was obtained and cultured for 3 days at 25 °C on V8 agar plates in full darkness for subsequent experiments, in which V8 agar medium was composed of 100 mL filtered V8 juice, 900 mL deionized water, 1.5 g CaCO₃ and 20 g agar. V8 medium, which was consisted of 100 mL filtered V8 juice, 0.02 g CaCO₃ and 900 mL deionized water, was prepared for a large number of mycelia of P. capsici.

Technical-grade CA was purchased from Sigma and its assay was ≥95%. It was dissolved in acetone and Tween 80 to 3 mg/mL for subsequent experiments, in which acetone and Tween 80 were used as solvents. The medical solution was stored at 4 °C in the dark. qRT-PCR Kit and fluorescence quantitative PCR kit were purchased from TaKaRa [TaKaRa Biotechnology (Dalian) Co., Ltd., China]. Q5 Hot Start High-Fidelity2 × Master Mix was purchased from New England BioLabs (United States).

In vitro effect of CA on mycelial growth of P. capsici

Mycelial growth inhibition was used to determine the inhibitory effect of CA on P. capsici according to Khan et al.’s method (2011). CA solution was serially diluted by sterile water with 0.4% (vol/vol) acetone and Tween 80 and the desired concentration of CA emulsions were prepared. one mL of CA emulsion was blended with nine mL of V8 agar medium at 40–45 °C to obtain the agar plates amended with CA at 50, 80, 120, 150 and 250 mg/L. The same volume of solvents only was considered as the control.

A fresh mycelial disc (five mm in diameter) was obtained from culturing 3 days P. capsici dish and placed on treated plates (90 mm in diameter) after this medium solidification. After six days, each colony diameter was measured by the cross method. And the average value of the colony diameter was used for later data analysis. Each test group was repeated three replicates. According to the formula, the inhibitory rate of CA on the growth of P. capsici was calculated. By analyzing the correlation between the logarithm of the CA concentration and the probability value of the P. capsici colony growth inhibition rate, the virulence regression equation was obtained. The effective concentration for 50% and 75% inhibition of P. capsici mycelial growth (EC₅₀, EC₇₅) were estimated by the toxicity regression equation (Hou et al., 2020). ${\displaystyle \text{Inhibition}\left(\text{\%}\right)=\left(\mathrm{C}-\mathrm{T}\right)∕\left(\mathrm{C}-5\right)\times 100}$

Among them, C and T are colony diameter in the control sample and the treated sample respectively.

Effect of CA on P. capsici mycelium morphology

Mycelium morphology was observed by scanning electron microscope. A mycelial disc was placed into Petri dish containing 10 mL of V8 agar medium with a final concentration of 140 mg/L CA. The same volume of solvents only was considered as the control. After 3 days, mycelia were cut into 0.5 cm*0.5 cm pieces and fixed in 2.5% (w/v) glutaraldehyde at room temperature for 24 h. In order to remove the glutaraldehyde solution, the samples were washed with 0.1M PBS (pH = 7.2). The samples were fixed with 1% O_SO₄ for 1–2 h. The fixed solution in the sample was washed away by sterile water. The samples were dehydrated by serial rinses with ethanol in water solutions with 30% ethanol for 15 min, 50% ethanol for 20 min, 70%, 80%, 90%, 95%, 100% ethanol for 20 min respectively. After samples were dried and sprayed with gold, P. capsici mycelium morphology was observed under Scanning electron microscope.

Sample preparation for proteomic analysis

P. capsici was incubated on V8 agar medium at 25 °C in the dark for 3 days, and then mycelial plugs were obtained at 1/3 of the edge of the colony. Five mycelial plugs were put into Petri dish with 10 mL V8 medium for 3 days, and the culture medium was complemented with CA to produce a final concentration 140 mg/L. Petri dish with no CA was made as control. The mycelia were collected by centrifugation (4 °C, 4 000g, 15 min) and then discarded the supernatant. The mycelia were washed for 3 times with PBS (pH = 7.4) buffer and stored at −80 °C.

Protein extraction and hydrolysis

Total protein was extracted by acetone/ trichloroacetic acid (TCA) precipitation and SDT lysis method (Méchin, Damerval & Zivy, 2007). Mycelia were pulverized with liquid nitrogen in a mortar and poured into five-time-volume of TCA/ acetone and mixed by vortex. Then they were placed at −20 °C for 5 h. After the supernatant was removed by centrifugation (6000 g for 40 min at 4 °C), the samples were dried in a fume hood to obtain dry powder. The samples were lysed by putting into SDT lysate. The supernatant was collected by centrifugation (14000 g for 40 min) and filtered with a 0.22 M filter. The protein concentration was determined by the BCA protein assay. The total protein concentration was diluted to 100 mM by adding DTT. The mixture was eluted with UA buffer, IAA buffer and dissolution buffer at room temperature and hydrolyzed with 0.1 g/L trypsin buffer at 37 °C for 16 h. Then the peptide fragment was desalted with C18 Cartridge and the lyophilized peptide was dissolved with 40 µL dissolution buffer.

iTRAQ sample labeling and LC-MS/MS analysis

Six samples (three biological replicates) were labeled with the iTRAQ tags as follows: solvent control (113, 114 and 115), CA treatments (116, 117 and 118). These peptides were incubated and centrifuged in vacuo until dried. The labeled samples were pooled and then graded by using AKTA Purifier 100. The chromatographic column was balanced with loading buffer A (10 mM KH₂PO₄ in 25% ACN). The samples were eluted by gradient elution in liquid chromatography. During elution process, the absorbance was detected at 214 nm, and each fraction was separated for one minute. About 30 fractions were obtained and desalted with C18 cartridge after being lyophilized.

Each sample was separated using Nano-HPLC system Easy-nLC. The column was equilibrated with 95% loading buffer A. After the sample was loaded by the autosampler, it was separated by the analytical column. A linear gradient was performed at a velocity of flow of 300 nL/min for 1 h. The loading buffer B (10 mM KH₂PO₄ and 500 mM KCl in 25%ACN) was also gradient elution in chromatography.

Proteomics data processing

This study referred to the uniprot Phytophthora protein database. The data was analyzed as Ross et al.’s (2004) described. In short, a normal distribution of iTRAQ ratios was obtained by log2 conversion. And then logarithm, population mean and standard deviation were normalized respectively. Proteins with 1.2-fold change between the two treatments (P < 0.05) were identified as differentially expressed proteins, in which were analyzed by GO and KEGG analysis. An average linkage hierarchical clustering analysis was conducted on interested proteins using Multiple Experiment Viewer (MEV v 4.9.0).

Parallel reaction monitoring (PRM) Analysis

After the same mass protein of each sample was digested by trypsin, the peptide was desalted on the C18 column. The TripleTOF 5600+ LC-MS/MS system (AB SCIEX) was used for targeted MS analysis in PRM. The mass resolution of MS1 scan and MS/MS scan were ∼35000 and ∼15000, respectively. MS data acquisition was carried out in DDA mode to obtain proteins and peptide precursor ions, which were identified by ProteinPilot software. PRM acquisition included a MS1 scan (250 ms), and then a target MS/MS scan, with a period of 1.3–3.3 s. Finally, target proteins were introduced into the software Skyline for processing.

RNA extraction and qRT-PCR analysis

The mycelia cultured in V8 medium for 3 days (solvent control, CA treatment) were collected in 1.5 mL centrifuge tubes, centrifuged (8000 g, 2 min), and then pulverized with liquid nitrogen in a mortar. Total RNA was extracted with TRIZOL. The RNA was disposed using RNase free DNAse (TaKaRa) at 42 °C for 2 min for removing the DNA contamination. The synthetic cDNA was synthesized from 1 µg of total RNA by a reverse transcription kit. Reaction mixtures (12 µL) contained 6 µL of SYBR Green (TaKaRa), 0.5 μL of forward gene-specific primer, 0.5 µL of reverse gene-specific primer, 1 µL of cDNA, and 4 µL of ddH₂O. Primers were listed in Table S1. The G3PDH gene was used as the reference.

Data analysis

All experiments were performed thrice. The data were analyzed using one-way analysis of variance (ANOVA) with PASW Statistics 20.0. Differences with P values <0.05 were considered statistically significant.

Inhibitory effect of CA on mycelial growth of P. capsici

The CA obviously effect the mycelial radical elongation of P. capsici (Fig. 2). CA of higher concentration showed stronger inhibitory effect against P. capsica. The toxic regression equation of the CA for P. capsici was Y =  − 6.1485 + 5.5080X. Moreover, the EC₅₀ and EC₇₅ values of CA against P. capsici were 105.69 and 140.12 mg/L respectively.

Effect of CA on the mycelial morphology

After treated with CA (140 mg/L) for 3 days in V8 agar medium, the results of scanning electron microscopy (SEM) observations revealed that the mycelium growth of P. capsici was significantly limited (Fig. 3). In addition, the mycelium was distorted and shriveled. The embranchment increased. While, the mycelium of the control group was stretched, uniform and full. And the surface of mycelium was smooth and the diameter of mycelium was constant (Fig. S1). Since the CA of EC₇₅ concentration has obvious inhibition on the mycelium growth of P. capsici, the EC₇₅ value was selected as the test concentration in this study.

Overview of the quantitative proteomics analysis of P. capsici treated with CA

By the iTARQ-LC-MS/MS analysis, a total of 4370 proteins were identified in P. capsici treated with CA (140 mg/L) for 3 days. Most proteins have a molecular weight between 10–70 kDa, in which the number of proteins with molecular weight between 30–40 kDa was in the largest proportion (Fig. 4).

Effect of CA on protein levels

With the cut-off values P ≤ 0.05 and fold changes >1.2 or ≤0.83, proteins between the two samples were considered differentially expressed. Compared CA treatment with control group, a total of 1502 differentially expressed proteins were identified of which 647 proteins were up-regulated and 855 down- regulated based on the cut-off value (Table S2).

P. capsici proteins related to the mode of action of CA

Go analysis was performed for level 2, and the differentially expressed proteins were categorized according to their biological process, molecular function and cellular component. The biological process annotation revealed that the differentially expressed proteins of P. capsici treated with CA were involved in response to metabolic process (27%), cellular process (25%) and single-organism process (20%) (Fig. 5A). In the molecular function, as shown in Fig. 5B, the number of differentially expressed proteins involved in catalytic activity and binding, which accounted for 46% and 44%, respectively. In the category of cell component (Fig. 5C), the differentially expressed proteins involved in cell and cell part were the most, accounting for 23% respectively. So, CA mainly caused significant changes in protein expression related to metabolism, cell process, catalytic activity, protein binding and cell and cell separation in P. capsici, indicating that these differential proteins may be involved in the main biochemical pathway of CA acting on the target pathogen.

KEGG pathway analysis was used to analyze possible links between a protein and biochemical pathways. The differentially expressed proteins of P. capsici treated with CA were related to degradation pathways, including caprolactam degradation, other glycan degradation and valine, leucine and isoleucine degradation; associated with fatty acid pathways, including fatty acid metabolism and fatty acid elongation and involved with other pathways, including glyoxylate and dicarboxylate metabolism, hypertrophic cardiomyopathy, neomycin, kanamycin and gentamicin biosynthesis, propanoate metabolism (Fig. 6). These results indicated that fatty acid metabolism and the degradation of valine, leucine and isoleucine might play an important role in the inhibition of P. capsica by CA.

Functional classification of interested differentially expressed proteins

To determine the antimicrobial action of CA against P. capsici, we used the GO terms data and KEGG data to analyze the biochemical pathways involved in differentially expressed proteins. Analysis demonstrated that the 12 interested proteins involved in fatty acid metabolism (3 down-regulated, 3 up-regulated), polysaccharide metabolic (4 down-regulated) and leucine metabolism (2 down-regulated) were also significantly altered (Table 1). As shown in Fig. 7, the heat map of hierarchical clustering verified the rationality and accuracy of the selected differentially expressed proteins.

PRM verification of interested proteins

PRM technology based on mass spectrometry was used to quantitatively verify protein. Five interested proteins including glucan 1,3-beta-glucosidase, 1,3-beta-glucanosyltransferase, CAMK/CAMK1 protein kinase, methylcrotonoyl-CoA carboxylase subunit alpha and isovaleryl-CoA dehydrogenase were selected for initial validation. The PRM results showed that the five proteins of P. capsici treated with CA were down-regulated (Table 2), which was consistent with the results of proteomics analysis.

Quantitative of qRT-PCR analysis

The results of qRT-PCR showed that the gene expression level of the interested proteins was consistent with the protein expression level, except for CAMK/CAMK1 protein kinase, acetyl-CoA carboxylase and fatty acid synthase subunit alpha (Fig. 8). Moreover, the protein expression of cellulose synthase 3 decreased 1.43-fold while its RNA expression decreased 15.73-fold. The discrepancy between protein expression and RNA expression may be caused by the time difference between transcription and translation.