ORIGINAL RESEARCH article

Front. Microbiol., 12 October 2017

Sec. Food Microbiology

Volume 8 - 2017 | https://doi.org/10.3389/fmicb.2017.01989

iTRAQ Proteomic Analysis Reveals That Metabolic Pathways Involving Energy Metabolism Are Affected by Tea Tree Oil in Botrytis cinerea

  • Department of Food Science and Engineering, Ningbo University, Ningbo, China

Abstract

Tea tree oil (TTO) is a volatile essential oil obtained from the leaves of the Australian tree Melaleuca alternifolia by vapor distillation. Previously, we demonstrated that TTO has a strong inhibitory effect on Botrytis cinerea. This study investigates the underlying antifungal mechanisms at the molecular level. A proteomics approach using isobaric tags for relative and absolute quantification (iTRAQ) was adopted to investigate the effects of TTO on B. cinerea. A total of 718 differentially expression proteins (DEPs) were identified in TTO-treated samples, 17 were markedly up-regulated and 701 were significantly down-regulated. Among the 718 DEPs, 562 were annotated and classified into 30 functional groups by GO (gene ontology) analysis. KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis linked 562 DEPs to 133 different biochemical pathways, involving glycolysis, the tricarboxylic acid cycle (TCA cycle), and purine metabolism. Additional experiments indicated that TTO destroys cell membranes and decreases the activities of three enzymes related to the TCA cycle. Our results suggest that TTO treatment inhibits glycolysis, disrupts the TCA cycle, and induces mitochondrial dysfunction, thereby disrupting energy metabolism. This study provides new insights into the mechanisms underlying the antifungal activity of essential oils.

Introduction

Botrytis cinerea, one of the most destructive fungal pathogens, causing gray mold rot in a wide range of fresh fruits and vegetables. The resulting reduction in shelf life is responsible for enormous economic losses in the produce industry. Although chemical fungicides are widely used to control the incidence of the disease, this practice potentially introduces harmful substances into the food chain, and also selects for B. cinerea strains with increased drug resistance (Brul and Coote, 1999; Leroux et al., 2002). These limitations provide a strong stimulus to explore safer and more effective antifungal agents. Essential oils are promising natural substitutes that offer disease control by inhibiting pathogen growth (Prakash et al., 2012). For example, the essential oils of Angelica archangelica L. (Apiaceae) roots and Solidago canadensis L. have been characterized and tested in vitro as antifungal agents against B. cinerea (Fraternale et al., 2014; Liu et al., 2016). Lemongrass essential oil significantly reduces the incidence of B. cinerea and prolongs the shelf-life and sensory properties of frozen mussels and vegetables (Abdulazeez et al., 2016). Essential oils of aromatic plants, which belong to the Lamiacea family such as origanum (Origanum syriacum L. var. bevanii), lavender (Lavandula stoechas L. var. stoechas) and rosemary (Rosmarinus officinalis L.), have been reported to cause considerable morphological degenerations of the fungal hyphae of B. cinerea and suppress in vivo disease development on tomato against B. cinerea (Soylu et al., 2010).

Tea tree oil (TTO) is a volatile natural plant essential oil obtained from the leaves of the Australian tree Melaleuca alternifolia by vapor distillation (Homer et al., 2000). The oil exhibits a broad spectrum of antimicrobial activities against a variety of bacteria, fungi, and virus (Carson et al., 2006; Miao et al., 2016). Growth and metabolic activity of Escherichia coli and Candida albicans are inhibited after treatment with TTO (Gustafson et al., 1998; Bona et al., 2016). Our previous studies showed that TTO treatment effectively inhibits spore germination and mycelial growth of B. cinerea, modifies its morphology and cellular ultrastructure, and controls gray mold on strawberry and cherry fruits (Shao et al., 2013a; Li et al., 2017a). TTO's antifungal mechanism in B. cinerea involves the loss of membrane integrity and the subsequent release of intracellular compounds, probably due in part to changes in membrane fatty acid and ergosterol composition (Shao et al., 2013b; Li et al., 2017a). TTO also causes mitochondrial damage in B. cinerea, disrupting the tricarboxylic acid (TCA) cycle and leading to the accumulation of reactive oxygen species (ROS) (Li et al., 2017b). Metabolomic analysis by quadrupole time-of-flight mass spectrometer was consistent with these results (Xu et al., 2017). However, the molecular mechanisms underlying the effects of TTO against B. cinerea have not yet been associated with specific proteins.

Proteomics can be used to study the changes in protein levels under stress conditions in great detail (Franco et al., 2013), and has been applied to investigate the mode of action of the antimicrobial agent apidaecin IB against membrane proteins in E. coli cells (Zhou and Chen, 2011). Other studies have revealed that proteins related to energy and DNA metabolism, and amino acid biosynthesis are down-regulated in E. coli JK-17 in the presence of rose flower extract (Cho and Oh, 2011). Syzygium aromaticum essential oil perturbs the expression of virulence-related genes involved in the synthesis of serine protease, flagella, and lipopolysaccharide in Campylobacter jejuni (Kovács et al., 2016). In this study, we conducted a proteomics analysis using isobaric tags for relative and absolute quantification (iTRAQ) to study B. cinerea to identify proteins and potential mechanisms underlying the antifungal activity of TTO.

Materials and methods

B. cinerea growth and exposure to TTO

Highly virulent B. cinerea (ACCC 36028) was purchased from the Agricultural Culture Collection of China and grown at 25°C on potato dextrose agar (PDA, containing 1 L potato liquid, 20 g/L glucose, and 15 g/L agar) before use. TTO was purchased from Fuzhou Merlot Lotus Biological Technology Company (Fujian Province, China). The primary components of TTO are terpinen-4-ol (37.11%), γ-terpinene (20.65%), α-terpinene (10.05%), 1, 8-cineole (4.97%), terpinolene (3.55%), ρ-cymene (2.14%), and α-terpineol (3.82%), as specified by the supplier. B. cinerea cultures were maintained on PDA at 25°C for 3 days. Spore suspensions were harvested by adding 10 mL sterile 0.9% NaCl solution to each petri dish and then gently scraping the mycelial surface three times with a sterile L-shaped spreader to free the spores. The spore suspension was adjusted using a hemocytometer to 5 × 106 spores/mL. One milliliter suspension was inoculated into 250 mL flasks containing 150 mL sterile potato dextrose broth medium and cultured at 25°C on a rotary shaker at 150 revolutions per minute for 3 days. Before mycelia were harvested, TTO was added to the medium to a final concentration of 5 mL/L, and cultures incubated for another 2 h (Xu et al., 2017). Mycelia were collected and rinsed three times with 0.1 M phosphate buffered saline (PBS) (pH 7.4). Samples were stored at −80°C. Cultures without TTO were used as a control. Three samples were prepared in parallel for each condition.

Protein extraction

Approximately 200 mg of frozen mixed mycelium from control or TTO treated cultures was ground into powder in liquid nitrogen and suspended in 25 mL 10% (v/v) trichloroacetic acid in acetone containing 65 mM dithiothreitol (DTT). The suspension was vortexed and incubated at −20°C for 2 h, centrifuged at 12,000 × g for 45 min at 4°C, and the supernatant discarded. The precipitate was rinsed three times with chilled acetone. The pellet was vacuum dried and dissolved in lysis buffer (4% SDS, 100 mM Tris-HCl, 100 mM DTT, pH 8.0). After incubation for 5 min in boiling water, the suspension was sonicated on ice at 50 W for 5 min. The crude extract was incubated in boiling water again for 5 min, and clarified by centrifugation at 14,000 × g for 40 min at 20°C. To digest protein in the supernatant, 200 μL UA buffer (8 M urea, 150 mM Tris-HCl, pH 8.5) was added and the mixture was centrifuged at 14,000 × g for 30 min at room temperature. This step was repeated three times. Subsequently, 100 μL 50 mM iodoacetamide (IAM) was added, the samples were incubated for 30 min in darkness, and then centrifuged at 14,000 × g for 30 min at room temperature. The precipitate was resuspended in 100 μL UA buffer and samples were centrifuged at 14,000 × g for 30 min at room temperature. 100 μL dissolution buffer was added, followed by centrifugation at 14,000 × g for 30 min at room temperature. This step was repeated three times. The supernatant was removed, the pellet was dissolved in 40 μL trypsin buffer, incubated at 37°C for 18 h, and clarified by centrifugation at 14,000 × g for 30 min at room temperature. Finally, 40 μL 25 mM dissolution buffer was added and samples were centrifuged at 14,000 × g for 30 min at room temperature. The supernatant was transferred to a new tube and quantified with the Bradford assay using BSA as the standard, and SDS-PAGE was performed to verify protein quality.

iTRAQ labeling and strong cation exchange (SCX) fractionation

iTRAQ labeling was performed according to the manufacturer's instructions. Peptides were prepared using the 8-plex iTRAQ labeling kit (AB Sciex, CA, USA). Control replicates were labeled with reagents 113, 114, and 115, and the TTO treatment replicates were labeled with reagents 116, 117, and 118. The labeled peptide mixtures were pooled and dried by vacuum centrifugation.

The labeled peptide mixtures were dissolved in 3 mL buffer A (10 mM KH2PO4 in 25% acetonitrile, pH 3.0) and loaded onto a polysulfoethyl 4.6 × 100 mm column (5 μm, 200 Å, PolyLC, Inc., Maryland, USA). The peptides were eluted at a flow rate of 1 mL/min with a gradient of buffer A for 30 min, 5–70% buffer B (10 mM KH2PO4, 500 mM KCl in 25% acetonitrile, pH 3.0) for 65 min, and 70–100% buffer B for 80 min. The eluted peptides were pooled into 10 fractions, desalted on C18 cartridges (Sigma), and vacuum-dried.

LC-MS/MS analysis

For nano LC–MS/MS analysis, 10 μL of supernatant from each fraction was injected into an Obitrap-Elite (ThermoFinnigan) equipped with an Easy nLC (Proxeon Biosystems, now Thermo Fisher Scientific). The mobile phase was a mixture of water containing 0.1% formic acid and acetonitrile with 0.1% formic acid isocratically delivered by a pump at a flowrate of 250 nL/min. The elution gradient was: 0–105 min, 0–50% B; 105–110 min, 50–100% B; 110–120 min, 100% B. The MS scanning range was 300–1,800 m/z, MS resolution was 70,000, the number of scans range was 1, and the dynamic exclusion time was 40 s. The MS/MS activation type was HCD, the isolation window was 2 m/z, the MS/MS resolution was 17,500, the normalized collision energy was 30 eV, and the underfill ratio was 0.1%.

Analysis of differentially expression proteins

For protein quantitation, one protein was required to contain at least two unique peptides. The quantitative protein ratios were weighted and normalized by the median ratio in Mascot (http://www.matrixscience.com). When differences in protein expression between TTO-treated and control groups were >1.5-fold or <0.67-fold, with p < 0.05, the protein was considered to be differentially expressed.

Bioinformatic analysis

Gene Ontology (GO) is a standardized gene function classification system that describes the properties of proteins using three attributes: biological process, molecular function, and cellular components (Ashburner et al., 2000). A GO analysis (http://www.geneontology.org) was conducted to assign functional annotations for differentially expression proteins (DEPs), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg) was used to predict the primary metabolic and signal transduction pathways in which the identified DEPs are involved.

Confocal laser scanning microscopy

To assess the effects of TTO on the cytoplasmic membranes of B. cinerea, confocal laser scanning microscopy (LSM 880, Carl Zeiss, Germany) was performed, using the fluorescent indicator propidium iodide (PI) (Sigma-Aldrich, USA) and a modified protocol (Lee and Kim, 2017). B. cinerea cells containing 4 × 106 spores/ml were added to each glass tube and incubated with TTO (final concentration 5 mL/L) with shaking at 200 rpm at 25°C for 2 h. The cells were washed and resuspended in 0.5 mL PBS (pH 7.4), stained with PI (10 μM final concentration) for 30 min at room temperature in the dark, and then washed twice with PBS. Images were acquired using confocal laser scanning microscopy. The experiment was repeated three times.

Measurement of enzyme activities related to TCA cycle

Using the protocol described above (see Protein Extraction), ground mycelium was suspended in PBS (pH 7.4) and centrifuged at 10,000 × g for 10 min at 4°C. Enzyme activities were measured in the supernatant for malate dehydrogenase (MDH), citrate synthase (CS), and oxoglutarate dehydrogenase (OGDH), using kits purchased from Nanjing Jiancheng Bioengineering Institute (Nanjing, Jiangsu, China), following the manufacturer's instructions. Protein concentration was determined using a method based on the (Bradford, 1976) assay. MDH activity was calculated as μmol of NAD reduced per minute per mg of protein (U/mg protein). One unit of CS activity was defined as the amount of enzyme that produces 1 μmol of citric acid per minute (U/mg protein). OGDH activity was defined as the amount of enzyme that produces 1 nmol of NADH per minute (U/mg protein). Measurements were performed at 595 nm using three replicates for each sample.

Statistical analysis

All experiments were repeated three times. Mean values and standard deviations were calculated using Excel 2010 (Microsoft Inc., Seattle, WA, USA). Statistical analyses were performed using one-way ANOVA with SPSS Statistics 17.0 (SPSS Inc., Chicago, USA).

Results

Identification of B. cinerea proteins by iTRAQ

A total of 204,639 spectra were generated by iTRAQ proteomic analysis using control and TTO-treated B. cinerea and were analyzed using the Mascot search engine. As shown in Figure 1A, 17,337 spectra matched known spectra, comprising 10,001 peptides, 9,720 unique peptides, and 2,397 proteins from control and TTO-treated samples. The distribution of the number of peptides, predicted molecular weights, and isoelectric points, and peptide sequence coverage are shown in Figures 1B–D, respectively. Over 87% of the proteins were represented by at least two peptides. Molecular weights ranged from 20 to 200 kDa, and isoelectric points ranged from 5.0 and 7.0. Approximately 51% of identified proteins had more than 10% peptide sequence coverage.

Figure 1

Identification of differentially expressed proteins using iTRAQ

The threshold for differential expression (TTO-treated vs. control) was a protein level difference >1.5 or < 0.67, with a p < 0.05. 718 differentially expressed proteins were identified in the TTO sample, of which 17 were up-regulated and 701 were down-regulated. Details for each protein are provided in Table 1.

Table 1

AccessionProtein nameScoreSequence coverage (%)Foldap-value
gi|154691848cytochrome c96.337.90.3280.007
gi|347441783citrate synthase133.18.01.8190.028
gi|472236008malate dehydrogenase protein957.755.42.1200.017
gi|472241505oxoglutarate dehydrogenase protein698.327.21.6110.037
gi|347827327pyruvate carboxylase2, 263.638.71.7510.027
gi|347833674phosphoenolpyruvate carboxykinase548.730.21.6250.044
gi|347839725succinyl-CoA ligase subunit alpha420.324.31.6120.040
gi|347826865fructose-1,6-bisphosphatase308.139.11.6400.031
gi|154323902enolase2, 009.946.61.6210.008
gi|472238209glucose-6-phosphate isomerase protein574.229.91.9800.032
gi|472246374phosphoglycerate mutase protein54.32.61.5760.021
gi|4722404356-phosphofructokinase protein539.928.11.7750.022
gi|472237248bisphosphoglycerate-independent phosphoglycerate mutase protein823.044.12.1640.018
gi|347841748fructose-bisphosphate aldolase1, 045.242.21.7250.027
gi|536718572phosphoglycerate kinase 1587.540.21.7230.040
gi|347833674phospho-2-dehydro-3-deoxyheptonate aldolase548.730.21.8700.029
gi|347835540phosphoglycerate mutase family protein36.04.71.7920.015
gi|4722409746-phosphofructo-2-kinase fructose bisphosphatase protein98.89.41.8510.037
gi|347441437inosine 5-monophosphate dehydrogenase581.819.91.6060.020
gi|347841600adenine phosphoribosyltransferase182.437.81.7770.022
gi|347829189adenosine kinase465.931.31.9560.016
gi|347441679adenosylhomocysteinase1, 287.461.71.8810.027
gi|347837737S-adenosylmethionine synthetase423.130.12.0040.008
gi|347831618AMP deaminase 3111.14.51.6730.029
gi|347828730adenylosuccinate synthetase333.030.91.6020.036
gi|347837737S-adenosylmethionine synthetase423.130.12.0040.008
gi|347837845adenylyl cyclase-associated protein417.920.71.8100.022
gi|472242224guanyl-nucleotide exchange factor protein65.41.51.6740.004
gi|154691052uracil phosphoribosyltransferase90.69.41.7960.046
gi|154697015nucleoside diphosphate kinase522.442.81.9350.010
gi|347840376UTP-glucose-1-phosphate uridylyltransferase1, 333.645.71.6230.038
gi|347832865ribulose-phosphate 3-epimerase38.67.92.2040.031
gi|154300519alcohol dehydrogenase protein167.716.51.9600.026
gi|347836330alcohol dehydrogenase (NADP dependent)281.124.42.0190.020
gi|347441899zinc-containing alcohol dehydrogenase636.544.81.6560.032
gi|347440923aldehyde dehydrogenase1, 070.948.01.8650.021
gi|154703069ATP synthase D chain, mitochondrial252.126.41.9240.050
gi|563298521ATP synthase subunit e, mitochondrial60.29.91.7570.033
gi|347839842ATP citrate lyase subunit549.037.51.5890.023
gi|154703371vacuolar ATP synthase subunit E93.312.72.3820.013
gi|154692979vacuolar ATP synthase subunit D74.819.51.7150.024
gi|347441643vacuolar ATP synthase subunit H307.622.31.7610.028
gi|472245494vacuolar ATP synthase catalytic subunit a protein577.727.81.5800.012
gi|347835157v-type proton ATPase subunit B274.117.62.0410.019
gi|507414597mitochondrial import protein 131.18.61.8720.043
gi|472243251mitochondrial pyruvate dehydrogenase kinase protein61.43.42.6320.009
gi|229891130amino-acid acetyltransferase, mitochondrial44.22.12.1150.022
gi|3282211isocitrate lyase 1, partial27.82.51.8740.029
gi|347832197malate synthase46.45.71.8750.048
gi|347840647acetyl-CoA carboxylase2, 370.733.81.6220.039
gi|347842358acetyl-CoA acetyltransferase449.446.31.9820.018
gi|347841050fatty acid synthase1, 414.525.21.6930.042
gi|472245418fatty acid synthase beta subunit dehydratase protein1, 668.624.81.5670.045
gi|347841364NADP-specific glutamate dehydrogenase1, 138.846.91.8400.021
gi|347827914homocitrate synthase454.539.01.5010.031
gi|347837008homoserine kinase190.128.71.9200.042
gi|347836521GABA transaminase483.927.71.5440.018
gi|472242205aspartate aminotransferase protein385.826.11.8370.048
gi|347841990tryptophan synthase611.028.31.5420.024
gi|347832506threonine synthase348.416.41.5600.047
gi|154692095cysteine synthase292.825.01.5890.028
gi|347833148glutamine synthetase484.026.91.7780.015
gi|347839014histidine biosynthesis protein184.69.31.8400.027
gi|347828253dihydrodipicolinate synthetase family protein518.728.01.8690.013
gi|347836881D-3-phosphoglycerate dehydrogenase656.425.51.7580.018
gi|472242394saccharopine dehydrogenase protein338.736.21.7430.039
gi|347441047glycine dehydrogenase286.912.31.7080.029
gi|507414630C-1-tetrahydrofolate synthase905.431.21.7370.031
gi|347831191glutamate carboxypeptidase protein298.023.21.9770.020
gi|347841903methionine aminopeptidase 1221.220.32.0400.021
gi|332313356methionine aminopeptidase 273.110.32.0440.027
gi|347829817serine/threonine protein kinase32.64.41.6930.037
gi|472244536glutamate-cysteine ligase protein61.63.61.6980.037
gi|3478294875-methyltetrahydropteroyltriglutamate-homocysteine methyltransferase2, 013.241.12.5050.004
gi|347836712glycine cleavage system H protein116.022.01.9820.025
gi|472236211amino acid permease protein39.74.01.9990.031
gi|347830997peptide methionine sulfoxide reductase82.820.51.9190.029
gi|472243795aromatic-l-amino-acid decarboxylase protein287.612.11.9360.015
gi|347833024lysine decarboxylase-like protein79.78.61.5850.033
gi|472246546glutathione-dependent formaldehyde dehydrogenase587.548.71.6660.043
gi|347840830NADH-cytochrome b5 reductase305.523.01.5450.029
gi|347827019cytochrome P450 monooxygenase31.42.41.7220.042
gi|125949746calcineurin194.312.41.7770.023
gi|154289817chitin synthase129.24.71.5550.023
gi|347840218sorbitol dehydrogenase28.32.91.7060.028
gi|347440923aldehyde dehydrogenase1, 070.948.01.8650.021
gi|347833737mitochondrial peroxiredoxin Prx142.87.61.8560.044
gi|347828993antioxidant129.533.12.1270.028
gi|347839043superoxide dismutase163.117.01.7170.012
gi|166408944flavohemoglobin294.735.71.9940.009
gi|347828340oxidoreductase305.114.932.1190.045
gi|347841065nuclear control of ATPase protein84.74.70.2190.001
gi|347836808heat shock protein 703, 060.853.21.7500.014
gi|47224275330 kda heat shock protein296.747.51.9590.019
gi|347827157heat shock protein 901, 603.737.71.6500.032
gi|347830903heat shock protein STI1689.135.82.4510.011
gi|347830415heat shock protein Hsp881, 199.534.31.8170.020
gi|347833633heat shock protein748.334.31.9990.020
gi|154288804short chain dehydrogenase105.920.72.1420.005
gi|347840162translation initiation factor 3284.646.81.9050.031
gi|472245156eukaryotic translation initiation factor 3 subunit749.418.91.8900.015
gi|229463757eukaryotic translation initiation factor 3 subunit H195.820.71.8510.013
gi|229501208eukaryotic translation initiation factor 3 subunit K232.733.51.7510.044
gi|347841080eukaryotic translation initiation factor 2 subunit alpha193.517.11.5740.030
gi|347830243eukaryotic translation initiation factor 4e151.712.01.7980.044
gi|347840917actin-depolymerizing factor 1519.953.61.9590.018
gi|3182891actin1, 055.452.81.5550.035
gi|347831507actin binding protein276.916.61.9420.003
gi|347840551actin related protein 2/3 complex217.421.91.8350.013
gi|347838304F-actin capping protein beta subunit isoforms 1 and 2156.027.71.5950.044
gi|205716451actin cytoskeleton-regulatory complex protein end 3109.410.41.8270.022
gi|347827283actin lateral binding protein691.250.32.6210.002
gi|347441258myosin regulatory light chain cdc4327.643.91.7750.049
gi|347838471survival factor 1321.928.41.6080.038
gi|347441690transcription factor HMG78.821.83.5650.004
gi|347838526transcription factor CCAAT39.13.54.9700.001
gi|374093884transcription regulator PAC1, partial42.13.22.5010.023
gi|472235708cp2 transcription factor protein92.26.21.7480.040
gi|347826783transcription initiation factor subunit28.97.42.0830.024
gi|347837746transcription factor CBF/NF-Y46.16.11.8690.021
gi|347840266transcription factor Zn, C2H250.51.73.4070.003
gi|347837101EF-hand calcium-binding domain protein42.83.70.0310.001
gi|472246130cell division control protein cdc48 protein1, 298.440.71.6540.021
gi|472235945cell lysis protein103.720.51.9300.025
gi|206558271cell division cycle protein 12338.63.91.8090.050
gi|347828695apoptosis-inducing factor 3267.717.22.2900.003
gi|472242094thioredoxin protein388.651.42.6340.003
gi|472244889sulfate adenylyltransferase protein328.925.81.8580.011
gi|347839319protein disulfide-isomerase542.339.11.8620.031
gi|347442007transaldolase1, 216.650.21.9840.022
gi|154703303elongation factor 1-alpha2, 637.450.01.8310.034
gi|347830450elongation factor 21, 896.644.61.6880.020
gi|472244387elongation factor 1-beta protein597.240.02.0060.024
gi|347841449NAD-dependent formate dehydrogenase1, 663.050.11.9310.042
gi|34783578526S protease regulatory subunit 6A355.127.61.8480.017
gi|472242788proteasome component pre3 protein101.723.91.9420.023
gi|347841691arp2/3 complex subunit Arc16249.241.71.7290.020
gi|15431920726S protease regulatory subunit 7221.919.42.0090.026
gi|347833025proteasome subunit alpha type 1133.216.91.7060.025
gi|347441407protein kinase C substrate282.518.11.7030.028
gi|347827686sec14 cytosolic factor240.141.41.7110.030
gi|347840528peptidyl-prolyl cis-trans isomerase D431.339.92.0700.019
gi|563298153inorganic pyrophosphatase317.829.71.7140.015
gi|347830035aldose 1-epimerase338.429.62.1140.040
gi|347831189carbohydrate-Binding Module family 48 protein330.427.13.7440.014
gi|347839149carbohydrate-Binding Module family 50 protein196.525.32.2760.047
gi|347841295cystathionine beta-synthase416.026.01.7900.031
gi|347842143diphosphomevalonate decarboxylase303.625.91.7880.022
gi|347836348protein phosphatase PP2A regulatory subunit A414.121.11.5760.045
gi|347838932class I/II aminotransferase340.323.91.8440.015
gi|347831623amidophosphoribosyltransferase1, 467.620.81.5730.025
gi|472236449enoyl- hydratase isomerase protein101.119.11.8490.026
gi|472237246tubulin-specific chaperone c protein222.720.71.6210.044
gi|347826898trans-2-enoyl-CoA reductase31.91.90.0310.001
gi|3478378641,3,8-naphthalenetriol reductase89.019.62.2130.029
gi|472243905casein kinase i protein148.319.81.5910.043
gi|347831955acetate kinase193.118.91.7260.015
gi|347839614aspartyl aminopeptidase293.318.81.5640.036
gi|4722385383-hydroxybutyryl-dehydrogenase protein133.317.21.6450.025
gi|347441025arf gtpase-activating protein249.917.22.0740.008
gi|347828551phosphatidyl synthase72.69.41.9670.029
gi|154294387mitogen-activated protein kinase101.917.11.6640.039
gi|472240101alpha beta hydrolase fold-3 domain protein45.39.01.8120.020
gi|347827703BAR domain protein271.643.41.7510.037
gi|347830570ThiJ/PfpI family protein645.537.01.7030.016
gi|347832713DUF1688 domain-containing protein437.727.61.7260.034
gi|472245392DUF718 domain-containing protein75.627.31.8030.019
gi|347836108C2 domain-containing protein286.023.61.9470.021
gi|347833490DUF757 domain-containing protein74.822.41.8400.045
gi|472245612c6 finger domain protein248.422.41.7820.029
gi|347838618UBX domain-containing protein101.116.22.2520.036
gi|472236354yip1 domain-containing protein66.011.12.0520.033
gi|347836200FAD binding domain-containing protein117.410.72.0720.015
gi|347836441DUF89 domain-containing protein69.46.01.6380.027
gi|472240877bar domain-containing protein69.25.91.7840.040
gi|347832303acyl-CoA dehydrogenase domain protein202.219.92.0100.042
gi|472237107saff domain-containing protein94.88.51.9330.015
gi|347828586CUE domain-containing protein53.83.13.8330.008
gi|472244807calponin domain protein79.32.92.0670.033
gi|563296966KH domain protein31.21.71.9000.011
gi|347829378R3H domain-containing protein32.31.61.9380.001
gi|347836748pumilio domain-containing protein37.91.42.3130.007
gi|347836261methyltransferase domain-containing protein27.92.90.0310.001
gi|154691472eukaryotic peptide chain release factor subunit 1426.930.81.9120.036
gi|347837479glia maturation factor gamma102.730.61.7030.028
gi|347837628CORD and CS domain-containing protein134.329.81.7870.013
gi|347828828ruvB-like helicase 1417.530.41.5020.035
gi|347442085CND899.46.30.4050.001
gi|15605143040S ribosomal protein S31, 591.360.81.6380.040
gi|34782780540S ribosomal protein S5418.338.51.5310.044
gi|34783512040S ribosomal protein S6332.834.31.7630.046
gi|34783642940S ribosomal protein S7276.130.41.8570.007
gi|15604347140S ribosomal protein S8688.840.21.5840.026
gi|15429114540S ribosomal protein S10106.225.41.8910.016
gi|15606167940S ribosomal protein S13404.433.81.8670.035
gi|47223738440S ribosomal protein S18546.842.31.9020.018
gi|34783725040S ribosomal protein S19363.851.02.7150.025
gi|34744146740S ribosomal protein S21157.163.62.7620.018
gi|34782932640S ribosomal protein S23190.620.01.8980.048
gi|15606588140S ribosomal protein S24348.132.61.8610.040
gi|15606563340S ribosomal protein S25174.326.82.0730.037
gi|34783233340S ribosomal protein S27322.437.81.8230.028
gi|34782811840S ribosomal protein S29126.942.92.5080.013
gi|34782751340S ribosomal protein S3063.116.10.1990.002
gi|34782877160S ribosomal protein L4497.813.22.9190.014
gi|15606208460S ribosomal protein L91, 053.663.41.5710.031
gi|22989153654S ribosomal protein L4, mitochondrial54.36.80.3750.024
gi|15603753060S ribosomal protein L12608.940.01.5620.010
gi|34783240160S ribosomal protein L13444.733.01.6620.032
gi|34783580560S ribosomal protein L6611.833.01.6700.023
gi|34783624860S ribosomal protein L10126.511.32.3360.030
gi|34783976660S ribosomal protein L16271.729.72.0550.039
gi|15431625760S ribosomal protein L17563.930.52.1360.011
gi|15431024860S ribosomal protein L19409.829.42.6520.009
gi|34784017860S ribosomal protein L21247.935.61.9770.029
gi|34783098560S ribosomal protein L23425.648.91.9360.030
gi|34783553460S ribosomal protein L24274.029.02.2910.015
gi|34783134860S ribosomal protein L26236.536.82.1740.030
gi|34744154960S ribosomal protein L27a708.848.31.6030.018
gi|34784147460S ribosomal protein L28236.952.73.5930.010
gi|47224583160S ribosomal protein L31295.448.02.2300.019
gi|34782664860S ribosomal protein L33274.237.61.9090.034
gi|15431503960S ribosomal protein L35140.218.92.7440.024
gi|15603647460S ribosomal protein L36166.035.91.6480.038
gi|15429764860S acidic ribosomal protein P01, 277.741.71.8960.029
gi|34783523760S acidic ribosomal protein P1553.241.22.3790.011
gi|34783855860S acidic ribosomal protein P2500.455.92.1780.012
gi|347441053ribosome associated DnaJ chaperone Zuotin635.225.31.8630.029
gi|156044830ribosome biogenesis protein Nhp2106.99.81.5940.024
gi|229485392ribosome biogenesis protein erb156.14.21.6360.045
gi|347837666nuclear transport factor 2236.228.22.2490.020
gi|472246396nuclear segregation protein466.527.03.040.013
gi|347835094leucyl-tRNA synthetase722.125.91.8090.016
gi|347835240methionyl-tRNA synthetase183.718.91.9310.029
gi|347828755tryptophanyl-tRNA synthetase283.923.41.8640.037
gi|563295297histidyl-tRNA synthetase286.921.91.6910.027
gi|347835339glutamyl-tRNA synthetase353.021.51.7830.032
gi|347841257threonyl-tRNA synthetase522.018.21.9180.013
gi|347840344valyl-trna synthetase535.313.71.6810.046
gi|347833265aspartyl-tRNA synthetase271.915.11.8610.017
gi|347836347phenylalanyl-tRNA synthetase beta chain159.913.52.1480.003
gi|347842507tRNA methyltransferase31.72.91.7350.006
gi|347837080polyadenylate-binding protein621.119.81.7550.039
gi|563292520histone H1-binding protein84.17.01.8940.025
gi|472237673oxysterol-binding protein154.56.53.3780.014
gi|154692219glycogen synthase204.911.11.8840.038
gi|154308576glucose-6-phosphate 1-dehydrogenase365.525.11.9860.023
gi|3478330531,3-beta-glucan biosynthesis protein131.710.62.1310.033
gi|347841047plasma membrane stress response protein34.62.03.1950.009
gi|347830640methylenetetrahydrofolate reductase196.213.41.5520.019
gi|154309515ca/CaM-dependent kinase-1141.718.41.5660.036
gi|347829911GTP-binding nuclear protein Ran301.838.11.7320.025
gi|472236275tRNA splicing endonuclease subunit protein96.814.52.0130.007
gi|347831289RNA binding effector protein Scp160853.422.11.5680.050
gi|347839263DNA-directed RNA polymerase I subunit49.614.12.6620.041
gi|347441996HAD superfamily hydrolase203.132.51.5990.041
gi|347840552ubiquitin carboxyl-terminal hydrolase362.927.11.9760.026
gi|347837756ubiquitin-like protein SMT334.918.82.3010.030
gi|472238757ubiquitin-activating enzyme e1 1 protein489.317.31.6650.016
gi|154695558ubiquitin-conjugating enzyme E236.37.51.5790.042
gi|472241717ubiquitin thioesterase protein56.48.31.7490.027
gi|347440894translocon beta subunit Sbh1225.344.61.7530.042
gi|472236180minor allergen alt a 7 protein282.347.82.8440.005
gi|472235513anthranilate synthase component 2 protein392.720.71.5900.029
gi|347833273nipsnap family protein154.319.91.6330.026
gi|347832071phosphoglucomutase1, 936.253.11.8960.017
gi|347829895phosphomannomutase182.721.51.8540.028
gi|347832016N-acetylglucosamine-phosphate mutase436.926.41.8530.011
gi|347841616UDP-galactopyranose mutase549.033.12.1490.020
gi|347841593UDP-N-acetylglucosamine pyrophosphorylase519.935.01.9220.008
gi|472237006UDP-glucose 4-epimerase gal10 protein191.120.51.8670.009
gi|347441001mannose-1-phosphate guanyltransferase alpha-a584.136.31.6310.033
gi|472241485nad h-dependent d-xylose reductase xyl1 protein247.928.61.5410.046
gi|347828612transketolase1, 284.841.22.0200.013
gi|154321267phosphoketolase883.524.41.8360.042
gi|347842358acetyl-CoA acetyltransferase449.446.31.9820.018
gi|347830285phospho-2-dehydro-3-deoxyheptonate aldolase460.236.11.9500.027
gi|3478407153-isopropylmalate dehydratase593.029.81.5190.019
gi|347440697cyanide hydratase/nitrilase353.717.02.5510.012
gi|347832595aldo/keto reductase family oxidoreductase497.642.51.9990.018
gi|154322845aldo/keto reductase327.828.91.7240.044
gi|347838695nitroreductase family protein228.332.71.8930.018
gi|154293270glucose 1-dehydrogenase263.427.81.6360.043

The main differentially expressed proteins in B. cinerea after treatment with TTO.

a

Fold: the average ratio (control/TTO-treated) of protein levels from three biological replicates as determined by iTRAQ approach. A protein was considered a differential expression protein as it exhibited a >1.5-fold or < 0.67-fold change and P < 0.05.

GO analysis of DEPs

GO analysis was conducted to identify significantly enriched GO functional groups. DEPs were categorized by biological process, cellular component, and molecular function. Of the 718 DEPs, 562 were annotated and classified into 30 functional groups (Figure 2). Biological processes accounted for 12 GO terms (with “metabolic process” accounting for 44.11% of these, and “cellular process” 34.32%). Cellular components accounted for 7 GO terms, dominated by “cell” (31.60%) and “cell part” (31.60%). Molecular functions accounted for 11 GO terms, the most abundant being “catalytic” (44.72%) and “binding” (43.61%).

Figure 2

The agriGO analysis tool was used to detect and visualize significantly enriched GO terms associated with the 562 annotated proteins, with an adjusted p-value cutoff of 0.05. Significant functions included “regulation of biological quality” (GO:0065008, p = 0.033) and “primary metabolic process” (GO:0044238, p = 0.016). There are 5 DEPs, accounting for about 45.45% of the total protein in regulation of biological quality. And 189 DEPs, accounting for about 73.82% of the total protein in primary metabolic process.

KEGG analysis of DEPs

Proteins typically do not exercise their functions independently, but coordinate with each other to complete a series of biochemical reactions. Pathway analysis can help reveal cellular processes involved in disease mechanisms or drug action. Using the KEGG database as a reference, 562 DEPs were linked to 133 different pathways. Glycolysis, the TCA cycle, and purine metabolism were among the pathways most significantly altered by exposure to TTO.

Confocal microscopy

Confocal laser scanning microscopy was used to investigate B. cinerea cell membrane integrity after TTO treatment. PI easily penetrates a membrane-damaged cell and binds to DNA, resulting in red fluorescence. B. cinerea cells were examined by both bright-field microscopy (Figures 3A,C) and fluorescence microscopy (Figures 3B,D). Control cells have no detectable red fluorescence (Figure 3B), indicating that they have intact cell membranes. In contrast, red fluorescence was observed after cells were treated for 2 h with TTO at 5 mL/L (Figure 3D). These results suggest that TTO compromises the integrity of the B. cinerea cell membrane, potentially causing cell death.

Figure 3

Enzyme activities related to TCA cycle

Because the iTRAQ analysis clearly implicated the TCA cycle as a possible TTO target, we investigated the activities of MDH, CS, and OGDH, three key enzymes related to the TCA cycle (Figure 4). The results indicate that activities for these enzymes decreased significantly in TTO-treated cells (87.4, 53.3, and 40.3%, respectively), consistent with our observation that the MDH, CS, and OGDH proteins are significantly down-regulated in TTO-treated cells.

Figure 4

Discussion

The antifungal activity of essential oils is probably based on their ability to significantly reduce total lipid and ergosterol content, thereby disrupting membrane permeability and resulting in leakage of cell components such as ATP, DNA, and potassium ions (Tian et al., 2011; Tao et al., 2014; Cui et al., 2015). Our previous study demonstrated that TTO considerably increases membrane permeability, causing extrusion of abundant material (Shao et al., 2013b; Yu et al., 2015) and decreasing intracellular ATP in B. cinerea (Li et al., 2017b). In this study, observations using confocal laser scanning microscopy indicate that TTO damages the B. cinerea cell membrane, potentially causing the release of internal material such as ATP.

Levels for many DEPs related to glycolysis metabolism, such as glucose-6-phosphate isomerase, 6-phosphofructokinase, phosphoenolpyruvate carboxykinase, fructose-1, 6-bisphosphatase, and enolase, are decreased by TTO treatment (Table 1). Glucose-6-phosphate isomerase catalyzes the conversion of glucose-6-phosphate into fructose 6-phosphate in the second step of glycolysis (Achari et al., 1981). 6-phosphofructokinase is a key enzyme in the control of the glycolytic pathway in nearly all cells (Wang et al., 2016). The activity of this enzyme is controlled by several metabolites, most notably its two substrates, fructose 6-phosphate and ATP. Glycolysis is also an important pathway for energy production in the cytosol of plant cells. Our results suggest that TTO inhibits glycolysis and may affect energy supply in B. cinerea.

Mitochondria are the primary sites of aerobic respiration in eukaryotic cells. They generate energy for cellular functions through oxidative phosphorylation and the TCA cycle, and also play a crucial role in regulating the apoptosis (Shaughnessy et al., 2014). In this study, several proteins associated with the mitochondrial respiratory chain and TCA cycle, such as ATP synthase D chain, ATP synthase subunit e, MDH, CS, and OGDH, were significantly down-regulated in cells treated with TTO (Table 1). ATP synthase D chain and ATP synthase subunit e are involved in the biosynthesis of ATP. Dill oil inhibits mitochondrial ATPase activity and dehydrogenase activities, and affects mitochondrial function in Aspergillus flavus (Tian et al., 2012). Mustard essential oils decrease intracellular ATP and increase extracellular ATP in E. coli O157:H7 and Salmonella typhi (Turgis et al., 2009). Citral decreases intracellular ATP content, increases extracellular ATP content, inhibits the TCA pathway, and decreases the activities of CS and α-ketoglutarate dehydrogenase in Penicillium digitatum (Zheng et al., 2015). Our additional study demonstrates that TTO treatment significantly inhibits the activities of MDH, CS, and OGDH (Figure 4). In our previous study, we found that TTO decreases intracellular ATP and the activities of MDH, succinate dehydrogenase, ATPase, CS, isocitrate dehydrogenase, and α-ketoglutarate dehydrogenase, disrupting the TCA cycle in B. cinerea (Li et al., 2017b). The down-regulation of two MDHs suggests that the Krebs cycle is not completely functional in Paracoccidioides lutzii upon exposure to argentilactone (Prado et al., 2014). Together, these results imply that TTO affects proteins in B. cinerea involved in glycolysis, the TCA cycle, and ATP synthesis, thereby disrupting the TCA cycle, interrupting energy metabolism, and inducing mitochondrial dysfunction.

Cytochrome c (cyt c) is a hemoglobin located in the inner mitochondrial membrane, and is responsible for transferring electrons between mitochondrial electron transport chain complexes III and IV (Reed, 1997; Lo et al., 2017). ATP is produced by the aerobic mitochondrial respiratory chain. Abnormal cyt c disrupts the mitochondrial respiratory chain and impacts ATP production (Zhou et al., 2015). Our study shows that cyt c is up-regulated in B. cinerea after TTO treatment at 5 mL/L (Table 1). The increase in cyt c levels may improve the performance of the oxidative respiratory chain, perhaps as a protective response to TTO toxicity.

Purines are one of the building blocks for nucleic acids. Their synthesis pathways generate many kinds of energy molecules (Qian et al., 2014). Inosine 5′-monophosphate dehydrogenase (IMPDH) is a rate-controlling enzyme in the de novo synthesis of the guanine nucleotide, and plays crucial roles in cell growth and proliferation (Fotie, 2016). IMPDH inhibition reduces guanine nucleotide pools and interrupts cellular functions such as DNA replication, RNA synthesis, and signal transduction (Weber, 1983; Weber et al., 1996). These effects are associated with cell cycle disruption, cellular differentiation, and apoptosis (Vitale et al., 1997; Yalowitz and Jayaram, 2000). Nucleoside diphosphate kinases (NDPK) are critical enzymes related to the maintenance of intracellular nucleotide levels, and catalyze the conversion of nucleoside triphosphates to nucleoside diphosphates in all living organisms (Véron et al., 1994). Both NDPK and AK can mediate the conversion of adenosine into ATP by ADP and AMP (Senft and Crabtree, 1983). In our study, TTO treatment decreased IMPDH levels (Table 1). Furthermore, levels of adenosine kinase AK and NDPK were also reduced after TTO treatment (Table 1). From these results, we can conclude that TTO may block the accumulation of energy and disrupt the cell cycle, ultimately inducing apoptosis.

Conclusion

The effect of TTO treatment on proteins in B. cinerea is summarized in Figure 5. We found that important metabolic pathways, including glycolysis, the TCA cycle, and purine metabolism, were compromised by TTO treatment, while cyt c increased. We conclude that the disruption of energy metabolism by TTO contributes to its antifungal activity against B. cinerea.

Figure 5

Statements

Author contributions

JX and XS designed the experiments. JX and YW performed the experiments. FX and HW analyzed the data. JX, XS, and HW drafted the manuscript. All authors read and approved the final manuscript.

Acknowledgments

This study was funded by the National Science Foundation of China (No. 31371860), the Public Welfare Applied Research Project of Zhejiang Province (No. 2017C32010), the Science and Technology Program of Ningbo City (2017C10065), the School Research Project (XYL17014), and the K.C. Wong Magna Fund in Ningbo University.

Conflict of interest

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.

References

  • 1

    AbdulazeezM. A.AbdullahiA. S.JamesB. D. (2016). Chapter 58 – Lemongrass (Cymbopogon spp.) Oils, in Essential Oils in Food Preservation Flavor and Safet (London, UK: Academic Press Books-Elsevier), 509516.

  • 2

    AchariA.MarshallS. E.MuirheadH.PalmieriR. H.NoltmannE. A. (1981). Glucose-6-phosphate isomerase. Philos. Trans. R. Soc. London. 293, 145157. 10.1098/rstb.1981.0068

  • 3

    AshburnerM.BallC. A.BlakeJ. A.BotsteinD.ButlerH.CherryJ. M.et al. (2000). Gene Ontology: tool for the unification of biology. Nat. Genet.25, 2529. 10.1038/75556

  • 4

    BonaE.CantamessaS.PavanM.NovelloG.MassaN.RocchettiA.et al. (2016). Sensitivity of Candida albicans to essential oils: are they an alternative to antifungal agents?J. Appl. Microbiol. 121, 15301545. 10.1111/jam.13282

  • 5

    BradfordM. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein–dye binding. Anal. Biochem.6, 31713188. 10.1016/0003-2697(76)90527-3

  • 6

    BrulS.CooteP. (1999). Preservative agents in foods: mode of action and microbial resistance mechanisms. Int. J. Food Microbiol. 50, 117. 10.1016/S0168-1605(99)00072-0

  • 7

    CarsonC. F.HammerK. A.RileyT. V. (2006). Melaleuca alternifolia (Tea Tree) oil: a review of antimicrobial and other medicinal properties. Clin. Microbiol. Rev. 19, 5062. 10.1128/CMR.19.1.50-62.2006

  • 8

    ChoY. S.OhK. H. (2011). Cellular and proteomic responses of Escherichia coli JK-17 exposed to the Rosa hybrida flower extract. Biotechnol. Bioprocess Eng. 16, 885893. 10.1007/s12257-011-0051-5

  • 9

    CuiH.ZhangX.ZhouH.ZhaoC.LinL. (2015). Antimicrobial activity and mechanisms of Salvia sclarea essential oil. Bot. Stud. 56, 16. 10.1186/s40529-015-0096-4

  • 10

    FotieJ. (2016). Inosine 5'-Monophosphate Dehydrogenase (IMPDH) as a potential target for the development of a new generation of antiprotozoan agents. Mini Rev. Med. Chem. 10.2174/1389557516666160620065558 [Epub ahead of print].

  • 11

    FrancoC.SoaresR.PiresE.KociK.AlmeidaA. M.SantosR.et al. (2013). Understanding regeneration through proteomics. Proteomics13, 686709. 10.1002/pmic.201200397

  • 12

    FraternaleD.FlaminiG.RicciD. (2014). Essential oil composition of Angelica archangelica L. (Apiaceae) roots and its antifungal activity against plant pathogenic fungi. Plant Biosyst. 150, 558563. 10.1080/11263504.2014.988190

  • 13

    GustafsonJ.LiewY.ChewS.MarkhamJ.BellH.WyllieS.et al. (1998). Effects of tea tree oil on Escherichia coli. Lett. Appl. Microbiol. 26, 194198. 10.1046/j.1472-765X.1998.00317.x

  • 14

    HomerL. E.LeachD. N.LeaD.LeeL. S.HenryR. J.BaverstockP. R. (2000). Natural variation in the essential oil content of Melaleuca alternifolia Cheel (Myrtaceae). Biochem. Syst. Ecol. 28, 367382. 10.1016/S0305-1978(99)00071-X

  • 15

    KovácsJ. K.FelsoP.MakszinL.PápaiZ.HorváthG. Á. H.et al. (2016). Antimicrobial and virulence modulating effect of clove essential oil on the food-borne pathogen Campylobacter jejuni. Appl. Environ. Microbiol. 82, 61586166. 10.1128/AEM.01221-16

  • 16

    LeeH. S.KimY. (2017). Paeonia lactiflora inhibits cell wall synthesis and triggers membrane depolarization in Candida albicans. J. Microbiol. Biotechnol. 27, 395404. 10.4014/jmb.1611.11064

  • 17

    LerouxP.FritzR.DebieuD.AlbertiniC.LanenC.BachJ.et al. (2002). Mechanisms of resistance to fungicides in field strains of Botrytis cinerea. Pest Manag. Sci. 58, 876888. 10.1002/ps.566

  • 18

    LiY.ShaoX.XuJ.WeiY.XuF.WangH. (2017a). Effects and possible mechanism of tea tree oil against Botrytis cinerea and Penicillium expansum in vitro and in vivo test. Can. J. Microbiol. 63, 219227. 10.1139/cjm-2016-0553

  • 19

    LiY.ShaoX.XuJ.WeiY.XuF.WangH. (2017b). Tea tree oil exhibits antifungal activity against Botrytis cinerea by affecting mitochondria. Food Chem. 234, 6267. 10.1016/j.foodchem.2017.04.172

  • 20

    LiuS.ShaoX.WeiY.LiY.XuF.WangH. (2016). Solidago canadensis L. essential oil vapor effectively inhibits Botrytis cinerea growth and preserves postharvest quality of strawberry as a food model system. Front. Microbiol. 7:1179. 10.3389/fmicb.2016.01179

  • 21

    LoY. T.HuangH. W.HuangY. C.ChanJ. F.HsuY. H. (2017). Elucidation of tRNA-cytochrome c interactions through hydrogen/deuterium exchange mass spectrometry. BBA Proteins Proteom. 1865, 539546. 10.1016/j.bbapap.2017.02.015

  • 22

    MiaoL.ZhuL.LiuB.DuL.JiaX.LiH.et al. (2016). Tea tree oil nanoemulsions for inhalation therapies of bacterial and fungal pneumonia. Colloids. Surface. B. Biointerfaces141, 408416. 10.1016/j.colsurfb.2016.02.017

  • 23

    PradoR. S.BailãoA. M.SilvaL. C.OliveiraC. M. A. D.MarquesM. F.SilvaL. P.et al. (2014). Proteomic profile response of Paracoccidioides lutzii to the antifungal argentilactone. Front. Microbiol. 6:616. 10.3389/fmicb.2015.00616

  • 24

    PrakashB.SinghP.KediaA.DubeyN. K. (2012). Assessment of some essential oils as food preservatives based on antifungal, antiaflatoxin, antioxidant activities and in vivo efficacy in food system. Food Res. Int. 49, 201208. 10.1016/j.foodres.2012.08.020

  • 25

    QianZ.LiuZ. L.KangN.WangA.ZengX.JianX. (2014). Genomic and transcriptome analyses reveal that MAPK- and phosphatidylinositol-signaling pathways mediate tolerance to 5-hydroxymethyl-2-furaldehyde for industrial yeast Saccharomyces cerevisiae. Sci. Rep. 4:6556. 10.1038/srep06556

  • 26

    ReedJ. C. (1997). Cytochrome c: can't live with it–can't live without it. Cell91, 559562. 10.1016/S0092-8674(00)80442-0

  • 27

    SenftA. W.CrabtreeG. W. (1983). Purine metabolism in the schistosomes: potential targets for chemotherapy. Pharmacol. Ther. 20, 341356. 10.1016/0163-7258(83)90031-1

  • 28

    ShaoX.ChengS.WangH.YuD.MungaiC. (2013b). The possible mechanism of antifungal action of tea tree oil on Botrytis cinerea. J. Appl. Microbiol. 114, 16421649. 10.1111/jam.12193

  • 29

    ShaoX.WangH.XuF.ChengS. (2013a). Effects and possible mechanisms of tea tree oil vapor treatment on the main disease in postharvest strawberry fruit. Postharvest Biol. Tec. 77, 94101. 10.1016/j.postharvbio.2012.11.010

  • 30

    ShaughnessyD. T.McAllisterK.WorthL.HaugenA. C.MeyerJ. N.DomannF. E.et al. (2014). Mitochondria, energetics, epigenetics, and cellular responses to stress. Environ. Health. Perspect. 122, 12711278. 10.1289/ehp.1408418

  • 31

    SoyluE. M.KurtS.SoyluS. (2010). In vitro and in vivo antifungal activities of the essential oils of various plants against tomato grey mould disease agent Botrytis cinerea. Int. J. Food Microbiol. 143, 183189. 10.1016/j.ijfoodmicro.2010.08.015

  • 32

    TaoN.OuYangQ.JiaL. (2014). Citral inhibits mycelial growth of Penicillium italicum by a membrane damage mechanism. Food Control41, 116121. 10.1016/j.foodcont.2014.01.010

  • 33

    TianJ.BanX.ZengH.HeJ.ChenY.WangY. (2012). The mechanism of antifungal action of essential oil from dill (Anethum graveolens L.) on Aspergillus flavus. PLoS ONE7:e30147. 10.1371/journal.pone.0030147

  • 34

    TianJ.BanX.ZengH.HuangB.HeJ.WangY. (2011). In vitro and in vivo activity of essential oil from dill (Anethum graveolens L.) against fungal spoilage of cherry tomatoes. Food Control22, 19921999. 10.1016/j.foodcont.2011.05.018

  • 35

    TurgisM.HanJ. J.CailletS.LacroixM. (2009). Antimicrobial activity of mustard essential oil against Escherichia coli O157:H7 and Salmonella typhi. Food Control20, 10731079. 10.1016/j.foodcont.2009.02.001

  • 36

    VéronM.TepperA.HildebrandtM.LascuI.LacombeM. L.JaninJ.et al. (1994). Nucleoside diphosphate kinase: an old enzyme with new functions?Adv. Exp. Med. Biol. 370, 607611. 10.1007/978-1-4615-2584-4_126

  • 37

    VitaleM.ZamaiL.FalcieriE.ZauliG.GobbiP.SantiS.et al. (1997). IMP dehydrogenase inhibitor, tiazofurin, induces apoptosis in K562 human erythroleukemia cells. Cytometry30, 6166. 10.1002/(SICI)1097-0320(19970215)30:1<61::AID-CYTO9>3.0.CO;2-I

  • 38

    WangJ.ZhangP.ZhongJ.TanM.GeJ.TaoL.et al. (2016). The platelet isoform of phosphofructokinase contributes to metabolic reprogramming and maintains cell proliferation in clear cell renal cell carcinoma. Oncotarget7, 2714227157. 10.18632/oncotarget.8382

  • 39

    WeberG. (1983). Enzymes of purine metabolism in cancer. Clin. Biochem. 16, 5763. 10.1016/S0009-9120(83)94432-6

  • 40

    WeberG.PrajdaN.YangH.YehY. A.ShenF.SinghalR. L.et al. (1996). Current issues in the regulation of signal transduction. Adv. Enzyme Regul. 36, 3355. 10.1016/0065-2571(96)00003-9

  • 41

    XuJ.ShaoX.LiY.WeiY.XuF.WangH. (2017). Metabolomic analysis and mode of action of metabolites of tea tree oil involved in the suppression of Botrytis cinerea. Front. Microbiol. 8:1017. 10.3389/fmicb.2017.01017

  • 42

    YalowitzJ. A.JayaramH. N. (2000). Molecular targets of guanine nucleotides in differentiation, proliferation and apoptosis. Anticancer Res. 20, 23292338.

  • 43

    YuD.WangJ.ShaoX.XuF.WangH. (2015). Antifungal modes of action of tea tree oil and its two characteristic components against Botrytis cinerea. J. Appl. Microbiol. 119, 12531262. 10.1111/jam.12939

  • 44

    ZhengS.JingG.WangX.OuyangQ.JiaL.TaoN. (2015). Citral exerts its antifungal activity against Penicillium digitatum by affecting the mitochondrial morphology and function. Food Chem. 178, 7681. 10.1016/j.foodchem.2015.01.077

  • 45

    ZhouM.LiY.HuQ.BaiX. C.HuangW.YanC.et al. (2015). Atomic structure of the apoptosome: mechanism of cytochrome c- and dATP-mediated activation of Apaf-1. Genes Dev. 29, 23492361. 10.1101/gad.272278.115

  • 46

    ZhouY.ChenW. N. (2011). iTRAQ-coupled 2-D LC-MS/MS analysis of membrane protein profile in Escherichia coli incubated with apidaecin IB. PLoS ONE75:e20442. 10.1371/journal.pone.0020442

Summary

Keywords

iTRAQ, proteomics, essential oil, Botrytis cinerea, antifungal

Citation

Xu J, Shao X, Wei Y, Xu F and Wang H (2017) iTRAQ Proteomic Analysis Reveals That Metabolic Pathways Involving Energy Metabolism Are Affected by Tea Tree Oil in Botrytis cinerea. Front. Microbiol. 8:1989. doi: 10.3389/fmicb.2017.01989

Received

06 September 2017

Accepted

27 September 2017

Published

12 October 2017

Volume

8 - 2017

Edited by

Boqiang Li, Institute of Botany (CAS), China

Reviewed by

Jun Tian, Jiangsu Normal University, China; Soner Soylu, Mustafa Kemal University, Turkey

Updates

Copyright

*Correspondence: Xingfeng Shao

This article was submitted to Food Microbiology, a section of the journal Frontiers in Microbiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics