ABSTRACT
We present the results of a study using high-throughput whole-transcriptome sequencing (RNA-seq) and vibrational spectroscopy to characterize and fingerprint pathogenic-bacterium injury under conditions of unfavorable stress. Two garlic-derived organosulfur compounds were found to be highly effective antimicrobial compounds against Cronobacter sakazakii, a leading pathogen associated with invasive infection of infants and causing meningitis, necrotizing entercolitis, and bacteremia. RNA-seq shows changes in gene expression patterns and transcriptomic response, while confocal micro-Raman spectroscopy characterizes macromolecular changes in the bacterial cell resulting from this chemical stress. RNA-seq analyses showed that the bacterial response to ajoene differed from the response to diallyl sulfide. Specifically, ajoene caused downregulation of motility-related genes, while diallyl sulfide treatment caused an increased expression of cell wall synthesis genes. Confocal micro-Raman spectroscopy revealed that the two compounds appear to have the same phase I antimicrobial mechanism of binding to thiol-containing proteins/enzymes in bacterial cells generating a disulfide stretching band but different phase II antimicrobial mechanisms, showing alterations in the secondary structures of proteins in two different ways. Diallyl sulfide primarily altered the α-helix and β-sheet, as reflected in changes in amide I, while ajoene altered the structures containing phenylalanine and tyrosine. Bayesian probability analysis validated the ability of principal component analysis to differentiate treated and control C. sakazakii cells. Scanning electron microscopy confirmed cell injury, showing significant morphological variations in cells following treatments by these two compounds. Findings from this study aid in the development of effective intervention strategies to reduce the risk of C. sakazakii contamination in the food production environment and on food contact surfaces, reducing the risks to susceptible consumers.
INTRODUCTION
Cronobacter sakazakii is a motile, nonsporeforming Gram-negative facultative anaerobe (1). It was noted to be an opportunistic pathogen responsible for life-threatening forms of meningitis, sepsis, necrotizing colitis, bacteremia, and meningoencephalitis in preterm neonates and infants (2–6). The mortality rate for neonatal infections has been reported to be as high as 80% (7). C. sakazakii infections may also result in severe neurological sequelae such as hydrocephalus, quadriplegia, and retarded neural development in survivors (8). Epidemiological studies implicated powdered infant formula as a primary source of transmission (6); indeed, this bacterium was detected in 2.4% to 14.2% of powdered infant formula products (5). Antibiotic therapy with a combination of ampicillin and gentamicin is an efficient method to treat C. sakazakii infections (4). However, the emergence of strains resistant to antibiotics has led the design and use of novel antimicrobials. In recent years, plant-derived biologically active compounds, such as trans-cinnamaldehyde (2, 3), phenolic benzaldehyde (9), curcumin (10), and anthocyanin (11), have attracted more attention due to strong antimicrobial activity, particularly for effective inhibition of the survival and growth of multidrug-resistant bacteria (12).
Garlic (Allium sativum L.), besides its function as spice and food, has been used as medicinal plant for over 4,000 years. It contains more than 100 biologically active secondary metabolites (13). Previous research demonstrated that garlic may be effective against cardiovascular disease because of its antithrombotic, hypolipidemic, hypocholesterolemic, antihypertensive, antidiabetic, and antihyperhomocysteinemia effects (14, 15). Further, garlic possesses many other biological activities, including antimicrobial, antioxidant, anticancer, antimutagenic, antiasthmatic, and anti-inflammatory immune modulation activities and prebiotic activities (14–22). Evidence suggests that the biological and medical functions of garlic are mainly due to its high level of organosulfur compounds (15), while proteins derived from garlic do not have antimicrobial activity (23) and phenolic compounds make only a minor contribution (24).
The primary sulfur-containing constituents in garlic are S-alk(en)yl-l-cysteine sulfoxides, such as alliin [(+)-(S)-allyl-l-cysteine-sulfoxide] and γ-glutamyl cysteines (15, 16). When a raw garlic bulb is crushed, alliin and alliinase, which are stored in separate cellular compartments of the garlic clove, are mixed together. In the role of alliinase, alliin is further subjected to dehydration by pyridoxal phosphate and transformed to allyl sulfenic acid, pyruvic acid, and ammonia (16). Allyl sulfenic acid is unstable and very reactive at room temperature. With the elimination of water, two molecules of allyl sulfenic acid condense spontaneously to allicin (see Fig. S1A in the supplemental material), validated as an important substance for the medicinal properties of garlic (13, 15–17). However, allicin is an unstable compound which is rapidly decomposed into other stable sulfur-containing constituents. The manufacturing procedure and garlic varieties are two important factors that affect the categories and compositions of organosulfur compounds decomposed from allicin (15, 25). Specifically, ajoene [(E)- and (Z)-4,5,9-trithiadodeca-1,6,11-triene-9-oxides] is commonly found in oil-macerated/ether-extracted garlic oil, with E-ajoene usually present at levels twice as high as those of Z-ajoene, while steam-distilled garlic oil is mainly composed of lipid-soluble sulfur compounds, including diallyl sulfide, diallyl disulfide, and diallyl trisulfide, and of minor amounts of many other volatile compounds, such as allyl methyl disulfide and diallyl tetrasulfide (see Fig. S1B) (15, 16, 23).
Both ajoene and lipid-soluble sulfur compounds (e.g., diallyl sulfides) have been reported to display strong antimicrobial activity against a number of Gram-positive and Gram-negative bacteria, yeasts, and viruses (15, 16, 19, 22–30). We recently validated the hypothesis that diallyl sulfides could effectively inactivate planktonic cells of Listeria monocytogenes, Escherichia coli O157:H7 (31), and Campylobacter jejuni (24) and sessile C. jejuni cells in a biofilm (18). Velliyagounder et al. also confirmed the in vitro antimicrobial efficacy against Aggregatibacter actinomycetemcomitans (23). Further, Yoshida et al. confirmed that Z-ajoene has stronger antimicrobial activity than E-ajoene, indicating the impact of the conformation of the compound on its biological activity (22, 32). Previous to this study, researchers had not examined the effect of the antimicrobial activity of ajoene and diallyl sulfide on C. sakazakii.
In the current study, we examined the in vitro antimicrobial effect of ajoene and diallyl sulfide on C. sakazakii. Additionally, the mechanism of sublethal injury of C. sakazakii under conditions of stress from these two organosulfur compounds was characterized and fingerprinted on both the genotypic and phenotypic levels using whole-transcriptome sequencing (RNA-seq) and Raman spectroscopy, respectively. Raman spectroscopy is based on vibrational spectroscopy and can be applied to characterize the biochemical components of bacterial cell membranes (33) and bacterial metabolites (34); thus, it is a reasonable technological tool to determine bacterial stress and injury in response to unfavorable environmental conditions. In fact, we have recently conducted vibrational spectroscopy coupled with multivariate analyses to determine bacterial stress and injury induced by cold and freezing storage (35) and by the presence of metal oxide nanoparticles (36) and of bioactive compounds derived from plants (18, 24, 31).
Several studies have validated the idea that DNA microarray analysis, a method relying upon hybridization of targeted oligonucleotides to particular loci, could be applied to investigate the global gene expression and transcriptomic response of food-borne pathogens to specific stress, such as chlorine treatment (37, 38) and hydrogen peroxide treatment (39). However, the relatively low detection sensitivity of DNA-DNA hybridization has limited microarrays to detection of only the trace level of bacterial transcriptomic response to the unfavorable stress. For example, it is very challenging to apply the DNA microarray technique to detect the rare transcripts and noncoding regulatory RNAs (40). Recently, RNA-seq has been applied to study the bacterial transcriptome (41–43) and further characterize the changes in pathogenic bacterial gene expression during infection of host cells (44). Further, transcriptome sequencing of Salmonella enterica serovar Enteritidis under conditions of desiccation and starvation stress in oil was studied using RNA-seq (40).
To the best of our knowledge, this is the first report of a study that used both high-throughput genotypic (i.e., RNA-seq) and phenotypic (i.e., Raman spectroscopy) tools to systematically study bacterial injury under unfavorable conditions. We used C. sakazakii as an exemplary microorganism for the study. Our innovative analytical approaches were employed to develop a model for illustrating the mechanism for stress and injury of C. sakazakii under unfavorable conditions.
MATERIALS AND METHODS
Chemicals and reagents.Diallyl sulfide (purity, 97%), crude diallyl disulfide (purity, 80%), and l-cysteine (purity, 97%) were purchased from Sigma-Aldrich (St. Louis, MO). Dimethyl sulfoxide (DMSO; purity, >99%) was obtained from Sangon Biotech Co., Ltd. (Shanghai, China).
Synthesis of ajoene.Ajoene was synthesized from diallyl disulfide (compound 1) via 2 steps (20) as follows:
Peracetic acid (35%; 2.58 g, 11.87 mmol) was added dropwise via an addition funnel to a solution of commercially available diallyl disulfide (Sigma-Aldrich; 80%) (2.1 ml, 11.3 mmol) in chloroform (18.0 ml) at 0°C. The reaction mixture was stirred at 0°C for 30 min, and anhydrous sodium carbonate (4.0 g) was added in small portions during 30 min with vigorous stirring. The mixture was stirred for an additional 30 min at 0°C and then filtered through a pad of Celite and anhydrous magnesium sulfate. The filtrate was concentrated under a vacuum and purified by the use of a chromatography column to get allicin (compound 2) (1.2 g, 71% yield) as a light yellow liquid. 1H nuclear magnetic resonance (NMR) (CDCl3) levels of 6.30 to 5.04 (m, 6 H) and 3.94 to 3.65 (m, 4 H) were determined as follows; the NMR data are in accordance with previous publication (20):
Allicin (compound 2) (1.0 g, 6.2 mmol) was dissolved in 40% aqueous acetone (10.0 ml), and the homogeneous solution was heated at 65°C for 4 h. The reaction mixture was diluted with 50% aqueous methanol (40.0 ml) and washed with hexane (5×; 20.0 ml). The aqueous methanolic layer was then saturated with ammonium sulfate and extracted with methylene chloride (20.0×; 2 ml). The methylene chloride extract was dried using magnesium sulfate, concentrated in vacuo, and purified by flash chromatography (silica gel/ethyl acetate), affording 0.10 g (21% yield) of Z-4,5,9-trithiadodeca-1,6,1l-triene (compound 3) (Z/E > 9/1). 1H NMR (CDCl3) levels of 6.55 (d, J = 9.0 Hz, 1 H), 6.14 to 5.36 (m, 7 H), and 3.62 to 3.23 (m, 6 H), which are in accordance with a previous publication (20), were determined.
Bacterial strains and culture methods.Five strains of C. sakazakii (ATCC BAA-894, ATCC 12868, ATCC 29004, IQCC 10423, and IQCC 10449) were used in this study. All strains were separately cultivated overnight in 5 ml Luria-Bertani (LB) broth (BD Difco) at 37°C without shaking to achieve a concentration of ca. 9 log CFU/ml, and 1 ml of each bacterial culture was centrifuged at 8,000 × g for 10 min at 4°C. The supernatant was discarded, and the bacterial pellets were washed twice with sterile phosphate-buffered saline (PBS) (pH ∼7.0 to ∼7.2), containing NaCl (9 g), Na2HPO4 (13.76 g), and NaH2PO4 (1.794 g) (per liter of distilled water). Equal volumes of the resultant culture suspensions, each containing one of the five strains, were then combined to form a cocktail, and the cocktail was subsequently added into sterile PBS and tryptic soy broth (TSB; BD Bacto), respectively. An initial concentration of approximately 5 log CFU/ml was achieved by serial dilution in PBS or TSB.
Antibacterial effects of garlic-derived organosulfur compounds on C. sakazakii.The synthesized Z-ajoene was dissolved in (filter-sterilized) DMSO at the initial molar concentration of 0.2 M and stored at 4°C without light exposure. Different volumes of Z-ajoene and diallyl sulfide were added into the bacterial culture systems, resulting in concentrations of Z-ajoene of 0.0777, 0.777, 3.88, and 7.77 mM and of diallyl sulfide of 0.777, 3.88, 7.77, and 77.7 mM. Control groups included bacterial cultures without addition of organosulfur compounds and bacterial cultures with addition of DMSO (0.55 M) only. Each inoculated sample was mixed well by a vortex procedure and incubated statically at 4, 22, and 37°C for 0, 1, 2, 4, 8, and 12 h (both PBS and broth samples). At each testing time, the bacterial samples were serially diluted in sterile PBS and spirally plated on tryptic soy agar (TSA; BD Difco). After incubation at 37°C for 24 h, viable cells were enumerated. All the tests were repeated in duplicate for one experiment, and the experiments were individually repeated at least three times.
Antibacterial effects of garlic-derived organosulfur compounds supplemented with l-cysteine on C. sakazakii.The 3.88 mM concentration was selected as the initial concentration of both ajoene and diallyl sulfide in the bacterial culture systems, and equal amounts of l-cysteine were immediately supplemented into the two systems to counteract the effects of the organosulfur compounds. The control groups consisted of bacterial suspensions with equal molar concentrations of DMSO and without any additives. C. sakazakii cultures were incubated at 37°C for 0, 1, 2, 4, 8, and 12 h in PBS. At each sampling time, the cultures were serially diluted with sterile PBS and spiral plated on TSA, and the numbers of viable cells were determined after incubation at 37°C for 24 h. All experiments were repeated in triplicate.
RNA-seq and real-time qPCR.Triplicate 7-ml samples of C. sakazakii at an optical density at 540 nm of 0.3 (OD540) were treated with 3.88 mM ajoene, 3.88 mM diallyl sulfide, or vehicle (i.e., DMSO) for 30 min at 22°C. The bacterial cultures were centrifuged at 8,000 × g for 5 min at 4°C, and the supernatants were discarded. Total RNA was extracted using an Ambion Bacterial RiboExpress RNA extraction kit (Life Technologies, Grand Island, NY). A portion of the extracted RNA was used to generate cDNA (Thermoscript III; Life Technologies), and quantitative PCR (qPCR) was performed on triplicate cDNA samples using Power SYBR green PCR Master Mix (Applied Biosystems, Warrington, United Kingdom) and an ABI Prism 7000 Fast instrument (Life Technologies). The rRNA in the sample was removed using a Ribo-Zero rRNA removal kit (Gram-negative bacteria) (Epicentre, Madison, WI). The purified mRNA was submitted to the Genomics Core of the Laboratory for Biotechnology and Bioanalysis, at Washington State University. The mRNA was sequenced using an Ion Torrent sequencing system (Life Technologies), and data were compiled and analyzed using CLC genomics workbench software (CLCBio, Cambridge, MA). The analyzed transcriptomes were sorted by false discovery rate (FDR) P values of less than 0.05 with a relative expression change of greater than 2-fold. The identified gene lists were submitted to DAVID for cluster analysis.
Raman instrumentation.A confocal micro-Raman spectroscopic system was used in this study. This system includes a Raman spectrometer (Renishaw, Gloucestershire, United Kingdom), a Leica microscope (Leica Biosystems, Wetzlar, Germany) and a frequency-doubled, Nd:YAG (λ = 532 nm) laser (Coherent, Santa Clara, CA). The spectrometer has an entrance aperture of 50 μm and a focal length of 300 mm and is equipped with 1,200-line/mm grating. The Rayleigh scattering can be eliminated by the filters, and only the Raman scattered light photons were collected and dispersed by a diffraction grating and finally recorded as a spectrum by a 576-by-384-pixel charge-coupled-device (CCD) array detector. The size of each pixel is 16 by 16 μm. Gold-coated microarray chips covered with C. sakazakii samples left untreated or treated with 3.88 mM garlic-derived organosulfur compounds (i.e., ajoene and diallyl sulfide) for 30 min at 22°C were mounted on a standard stage of the microscope, which was focused under the collection assembly, and spectra were recorded using a 20× Nikon objective (numerical aperture [NA] = 0.4, working distance [WD] = 3.9 mm) with a wave number range of 4,000 to 100 cm−1 in an extended mode. The measurement was conducted over 60 s of exposure time (3-s integration time, 20 accumulations) for 10 random spots on each sample with ∼2 mW of incident laser power. Prior to acquisition of each Raman spectrum, the bacterial sample was exposed to incident laser illumination for 30 s to reduce fluorescence background and achieve the highest signal-to-noise ratio (SNR). A Thorlabs PM100A optical power meter was used to measure and calibrate the incident laser power. WiRE 3.0 software (Renishaw, Gloucestershire, United Kingdom) was employed for instrumental control and spectral collection. The experiment was conducted in triplicate.
Spectral preprocessing.The preprocessing of the raw Raman spectra can separate and subsequently eliminate the side effects which may influence the quality of spectroscope-based chemometric models and multivariate analyses. We first conducted a polynomial background fit (45) combined with baseline subtraction using identification and discrimination of minima via adaptive and least-squares thresholding (46) to remove fluorescence background derived from C. sakazakii cells on gold-coated microarray slides, Gaussian noise, white noise, CCD background noise, and cosmic spikes (47–49). Besides fluorescence background, most of the spectral interference is contributed by CCD background noise, which is generated due to the thermal fluctuations on the CCD detector. Spectral binning (2 cm−1) and smoothing (9-point Savitzky-Golay algorithm) were subsequently applied, followed by normalization of the spectra based upon the intensity of the C-H band in the wave numbers of 3,100 to 2,950 cm−1, the indication of the total biomass of C. sakazakii cells. Recent studies conducting Raman spectral normalization of C. jejuni (48) and E. coli and Staphylococcus spp. (50) validated that using this C-H vibrational band as the standard for normalization effectively removed the spectral fluctuation derived from the small focal volume of bacterial cells and yielded the best results for spectral baseline correction.
Spectral reproducibility, sensitivity, and specificity.Spectral reproducibility from three independent experiments was investigated by calculating the differentiation index (Dy1y2) value using the following equations:
Spectral sensitivity and specificity from three independent experiments were determined using the Wards cluster algorithm at the cutoff value established at 99% similarity for the classification model (i.e., principal component analysis [PCA]) (49).
Second-derivative transformation analysis of spectra.The second-derivative transformed Raman spectrum is the result of applying a derivative transform to the raw Raman spectrum, and this algorithm function may swing with greater amplitude than that of the raw spectrum to magnify the minor spectral variation and separate out overlapping bands (51). In addition, it could be an ideal noise filter because the changes in the spectral baseline have an insignificant impact on second derivatives (49). It was calculated with the following formula (52):
Classification chemometric models.An unsupervised principal component analysis (PCA) model was conducted to investigate the variations in the samples without a priori knowledge about C. sakazakii with or without treatment of selective garlic-derived organosulfur compounds (i.e., ajoene and diallyl sulfide). The Mahalanobis distance in the PCA model is defined as the distance between groups in units of within-group standard deviations and was calculated using the following equation (24, 31):
A naive Bayesian classifier was employed to classify spectrum based upon the features defined on selected principal components (PCs) from PCA. The posterior of a category given features X of an observation is obtained by the Baye's theorem:
Here, p(X|wi) is assumed as a multivariate normal distribution with means and covariance estimated from samples in each category. The mean is estimated as follows:
The prior of each category is proportional to the number of samples in each category as follows:
SEM analyses.Scanning electron microscopy (SEM) was performed to examine morphological changes of C. sakazakii cells before and after treatment with a 3.88 mM concentration of the garlic-derived organosulfur compounds (i.e., ajoene and diallyl sulfide) in sterilized PBS for 1 h at 22°C. Untreated and treated C. sakazakii cells were harvested by centrifugation at 8,000 × g for 10 min at 4°C. C. sakazakii cells were fixed with 2.5% glutaraldehyde overnight at 4°C. The samples were then rinsed twice with 0.1 M phosphate buffer. Postfixation was carried out using 1% osmium tetroxide for 1 h prior to dehydration in an ethanol series (25%, 50%, 70%, 80%, 90% [10 min for each concentration], and 100% ethanol dehydrate [twice for 15 min each time]), and the reaction mixture was freeze-dried in a Christ alpha 1-4 lyophilizer (Christ, Osterode, Germany). The samples were mounted onto SEM stubs and sputter coated with a thin layer of gold. The coated samples were examined under a Leo 1530 Gemini field emission scanning electron microscope using an accelerating voltage of 5 kV (Zeiss, Jena, Germany).
Statistical analysis.Experiments were performed in at least three replicate trials. The results are expressed as the means of the results of three independent replicates ± the standard deviations. The differences between samples were shown to be significant (P < 0.05) by one-way analysis of variance (ANOVA) using Matlab.
RESULTS
Antibacterial effects of garlic-derived organosulfur compounds on C. sakazakii.Both ajoene and diallyl sulfide showed obvious antimicrobial effects on C. sakazakii. In general, ajoene has a higher antimicrobial effect than diallyl sulfide (an example is shown in Fig. S2A in the supplemental material). As shown in Table S1 in the supplemental material, the bactericidal effects were enhanced with an increase in molar concentrations of garlic-derived organosulfur compounds (i.e., ajoene and diallyl sulfide), and the antibacterial effect was also time dependent. When 0.777 mM ajoene was applied to inactivate bacterial cocktails in PBS at 37°C for 8 h, viable cells were not recovered. Along with the increase in the molar concentration of ajoene, the time for bacterial inactivation was reduced significantly (i.e., viable cells were reduced to nondetectable levels at 2 h for 3.88 mM and 1 h for 7.77 mM in PBS at 37°C). There was still a significant bactericidal effect even though the concentration of ajoene dropped to 0.0777 mM (see Table S1A). The bactericidal effect of ajoene at 22°C and 4°C was lower than that at 37°C, indicating that the volatility may contribute to the antimicrobial effect of ajoene. Although lower temperatures (e.g., 4°C) could alter the bactericidal effect of diallyl sulfide, the variation of the effect was not as remarkable as that seen with ajoene (see Table S1B). C. sakazakii cells could be completely inactivated in several minutes by the treatment with 77.7 mM diallyl sulfide at 4, 22, and 37°C.
When nutritional ingredients were present, the bactericidal effects of the garlic-derived organosulfur compounds (i.e., ajoene and diallyl sulfide) decreased significantly compared to incubation in PBS alone, indicating a protective effect from nutrients in broth medium. This protective effect was more obvious with the treatment using ajoene than with that using diallyl sulfide, resulting in a larger difference in the bacterial survival numbers in the two different matrices (i.e., PBS and TSB). The suppressive (bacteriostatic) effects of garlic-derived organosulfur compounds on a cocktail of C. sakazakii cells in TSB are also shown in Table S1 in the supplemental material, and a representative example is shown in Fig. S2B. C. sakazakii cells were inactivated within the first several hours of treatment with 7.77 mM diallyl sulfide in TSB at 37°C; however, incubation for 8 h allowed trace (i.e., nondetectable [ND; see Table S1 in the supplemental material]) levels of the injured and viable bacterial cells to recover and rebound to a certain high level (∼8 log CFU/ml at 12 h). Previous work has shown that only very high concentrations of allicin prevent bacterial cells from overcoming the inhibition and resuming growth, indicating that they were able to metabolize the allicin to a noninhibitory compound (56). This could be a plausible explanation for our current observation of the recovery of C. sakazakii cells treated by diallyl sulfide, a decomposed derivative of allicin.
Antibacterial effects of garlic-derived organosulfur compounds together with l-cysteine on C. sakazakii.Equal amounts of l-cysteine was added into the bacterial suspensions treated by garlic-derived organosulfur compounds (i.e., 3.88 mM ajoene and diallyl sulfide) in PBS at 37°C. Due to the addition of l-cysteine, the amounts of ajoene were decreased significantly because the sulfhydryl (−SH) group in l-cysteine is highly reactive with disulfide bonds in ajoene (see Fig. S3 in the supplemental material). Consequently, the majority of the antimicrobial activity of ajoene on C. sakazakii was diminished; however, the antimicrobial effect of diallyl sulfide did not change significantly (P > 0.05) (see Table S2 in the supplemental material).
Transcriptomic responses of C. sakazakii to ajoene and diallyl sulfide.To better understand the effect of ajoene and diallyl sulfide treatment, we examined changes in gene expression due to sublethal injury. Samples of C. sakazakii were added to TSB containing 3.88 mM ajoene or diallyl sulfide and incubated for 30 min at 22°C. Bacteria treated with vehicle (DMSO) served as a negative control. The RNA was extracted, purified, and depleted of rRNA. Patterns in gene expression were determined using RNA sequencing.
The RNA-seq results provided an adequate depth of genome coverage with an average of 10 reads per gene and a median read number of 437. The RNA isolated and sequenced from untreated bacteria and that from bacteria treated with ajoene and diallyl sulfide produced similar qualities of reads. Collectively, reads from all treatment conditions resulted in detection of all 4,392 chromosomal genes, as well as both plasmids, comprising an additional 165 genes. The accuracy of the sequencing results was validated with selected use of qPCR, targeting groEL, dnaK, oxyR, and ribH. The responses from the individual biological experiments tested by qPCR agreed with the pooled sequencing results (data not shown).
To identify patterns of differential gene regulation, we determined the genes that were differentially regulated at least 1.5-fold with an FDR-corrected P value of less than 0.05 (Table S3). Using DAVID analysis with the default settings, the differentially regulated genes were functionally categorized (57). Annotation clusters with enrichment scores of greater than 2.0 and GO terms with a Benjamini FDR of less than 0.05 are reported (Table 1 and Table 2). Under these conditions, we found that treatment with either ajoene or diallyl sulfide resulted in large changes in gene expression profiles of C. sakazakii (Fig. 1). Treatments using each of the two garlic-derived organosulfur compounds lead to similar directional changes in gene regulation. Ajoene-treated C. sakazakii was found to contain 637 genes with increased expression and 373 downregulated genes. Treatment of C. sakazakii with diallyl sulfide resulted in 1,434 genes with increased expression and 268 genes with decreased expression.
Regulation of genes in C. sakazakii treated by ajoene
Regulation of genes in C. sakazakii treated by diallyl sulfide
Treatment with either ajoene (A) or diallyl sulfide (B) resulted in a global change in gene expression in C. sakazakii.
Treatment with ajoene resulted in no significant pattern of increased gene expression as judged by Benjamini FDR (Table 1A). The downregulated genes were found to form four annotation clusters with an enrichment score of greater than 2.0. These clusters were NADH expression factors, flagellum and motility proteins, and nitrate reductases (Table 1B). The downregulation of both the nitrate reductases and the NADH expression factors seems to be related to a reaction to oxidative species. The downregulation of the flagellum is important as it contributes to the formation of the extracellular polymeric substances and serves as glue for the formation of biofilm (58).
Diallyl sulfide treatment results in the upregulation of coenzymes and quinones, cell shape and cell wall maintenance pathways, and lipopolysaccharide synthesis (Table 2A). These functional groups are all related to the bacterial membrane and may indicate an attempt to repair those structures. The genes that are downregulated in response to diallyl sulfide are RNA and amino acid biosynthesis biological clusters (Table 2B). This indicates a general decrease in the metabolic processes of this microorganism.
The difference between the two compounds in the patterns of gene expression is informative. Both treatments resulted in the downregulation of a set of 113 common genes, leaving 260 and 155 unique downregulated genes with ajoene and diallyl sulfide treatment, respectively. Further, 443 common genes were upregulated in both ajoene and diallyl sulfide treatments. Ajoene treatment resulted in the upregulation of 194 unique genes, while diallyl sulfide treatment resulted in the upregulation of 991 unique genes. The consistency of the regulation was seen when comparing the upregulated genes of one treatment condition with the downregulated genes of the other. However, 11 genes were upregulated following ajoene treatment and downregulated following diallyl sulfide treatment. Likewise, 55 genes were downregulated following ajoene treatment but upregulated by diallyl sulfide; among those genes are members of the glutamine synthesis pathway. Analysis of the 55 genes using DAVID showed a significant difference in the regulation results for two-component systems following the two treatments. Ajoene treatment resulted in the downregulation of five genes (ESA_01517, ESA_01523, glnG, glnA, and glnL), whereas diallyl sulfide treatment resulted in the upregulation of the same five genes. This indicates that ajoene reduces the ability of the bacteria to detect environmental changes whereas diallyl sulfide increases elements of environmental sensitivity, highlighting the effective differences between the two compounds.
Raman spectroscopy depicted C. sakazakii injury caused by ajoene and diallyl sulfide.Raman spectra of untreated C. sakazakii and C. sakazakii treated by the same molar concentration (i.e., 3.88 mM) of ajoene and diallyl sulfide were recorded and then converted to the second-derivative transformed patterns. The second-derivative transformed Raman spectral features were compared to study the mechanism of C. sakazakii injury under conditions of unfavorable stress (Fig. 2A, B, and C). The area and height of each second-derivative transformed Raman band were calculated using Matlab, and only the significant (P < 0.05) band variations were characterized and analyzed (Table 3). For diallyl sulfide treatment, the band at 1,638 cm−1 is assigned to amide I (both α-helix and β-sheet) (33) and the band at 1,297 cm−1 is derived from CH2 deformation of fatty acids (59). The band at the wave number of 511 cm−1 originated from S-S disulfide stretching (18, 24, 31). For ajoene treatment, the band at 2,914 cm−1 is from CH of lipids and proteins (59) and the band at 1,605 cm−1 is assigned to phenylalanine, tyrosine, and C = C groups of proteins (60). The band at the wave number of 675 cm−1 is assigned to S-S disulfide stretching and ν(C-S) of amino acid methionine (60). Compared to those seen after the treatments by organosulfur compounds, the spectra derived from untreated C. sakazakii showed significant band variations at the wave numbers of 1,522 cm−1 and 1,160 cm−1 assigned to carotenoids on cell membrane (60), 1,129 cm−1 derived from acyl backbone in lipids (59), and 750 cm−1 assigned to tryptophan (60).
C. sakazakii injury caused by diallyl sulfide and ajoene. C. sakazakii cells cultured overnight were statically treated with 3.88 mM ajoene (0.77 μl/ml) and 3.88 mM diallyl sulfide (0.5 μl/ml) in tryptic soy broth (TSB) at 22°C for 0.5 h. Panels A, B, and C show second-derivative transformations of Raman spectral features of untreated C. sakazakii (black) and of C. sakazakii treated with diallyl sulfide (red) and ajoene (blue), respectively. (D) Two-dimensional principal component analysis (PCA) for segregation of untreated C. sakazakii (C) and C. sakazakii treated with diallyl sulfide (DS) and ajoene (AJ). PCA was established using wave numbers of 1,800 to 400 cm−1 and 3,200 to 2,700 cm−1.
Band assignment of significant (P < 0.05) variations of second-derivative transforms of Raman spectral features of untreated C. sakazakii and C. sakazakii treated with diallyl sulfide and ajoene
Therefore, ajoene and diallyl sulfide demonstrated some differences in the antimicrobial mechanisms against C. sakazakii planktonic cells on the basis of biochemical features. Initially, for the first phase, both diallyl sulfide and ajoene bind to thiol-containing proteins/enzymes to generate a disulfide stretching band, resulting in denaturation of proteins/enzymes in bacterial cells. This was in accordance with our previous study results showing that Raman spectroscopy could monitor the transmembrane transfer of sulfur-containing compounds into bacterial cells and that diallyl sulfides could freely penetrate through bacterial cell membrane and selectively bind to thiol-containing proteins, altering their structures (24, 31). Kyung also demonstrated that most sulfur-containing compounds derived from cysteine sulfoxides in Allium species are believed to inhibit microorganisms through an identical mechanism by reacting with the −SH groups of cellular proteins to generate mixed disulfides, as their antimicrobial activities are diminished and finally abolished by continuously adding cysteine, except that seen with allyl alcohol (28). That finding is supported by our current study, which showed that the significant “sulfur” band variations at wave numbers of 511 cm−1 and 675 cm−1 for C. sakazakii treated by diallyl sulfide and ajoene, respectively, and by l-cysteine could significantly reduce or eliminate the antibacterial effect of ajoene on C. sakazakii (see Table S2 in the supplemental material). However, different variations in protein bands of C. sakazakii treated by diallyl sulfide and ajoene indicated different second phases to alter the secondary structures of proteins. Specifically, diallyl sulfide mainly altered the α-helix and β-sheet of amide I (band at 1,638 cm−1) whereas ajoene mainly changed the structures of phenylalanine and tyrosine (band at 1,605 cm−1) (Table 3). Diallyl sulfide is derived from steam-distilled garlic oil, and ajoene originates from oil-macerated/ether-extracted garlic oil, and these two compounds have variations in chemical structures which may be the critical factor contributing to their different antimicrobial mechanisms (in phase II) and effects against bacteria (i.e., C. sakazakii).
The result of the PCA classification determined using the Raman spectral features in the wave numbers of 3,200 to 2,700 cm−1 and 1,800 to 400 cm−1 is shown in Fig. 2D. The spectra collected from the first two experiments were employed to construct the model, and the spectra recorded from the third experiment were used to validate the model. Clear segregation of PCA models indicated a significant (P < 0.05) difference between the groups of samples (i.e., untreated C. sakazakii and C. sakazakii treated with ajoene and diallyl sulfide), with the interclass distance ranging from 4.01 to 12.83 on the basis of Mahalanobis distance measurements calculated between the centroids of classes. Clusters with interclass distance values greater than 3 are considered to be significantly different from each other (24, 55). This PCA cluster model achieved a classification sensitivity of 93.2% and a specificity of 95.1% using Wards cluster algorithm at the cutoff value established at 99% similarity. Bayesian probability was then conducted to compare the first 10 significant features with PCs determined by using PCA. Good agreement in the results of these two methods was noted, indicating the robust nature of the PCA model for classification of different samples on the basis of Raman spectral features.
Electron microscope-examined morphological variations of C. sakazakii treated with ajoene and diallyl sulfide.To correlate RNA-seq and Raman spectroscopic data with structural changes of C. sakazakii cells caused by treatment with garlic-derived organosulfur compounds, SEM data were collected for untreated and treated samples with 3.88 mM ajoene–PBS or diallyl sulfide–PBS for 1 h at 22°C. Figure S4 in the supplemental material shows that untreated bacterial cells had an intact cellular structure with well-defined membranes. Exposure to ajoene or diallyl sulfide resulted in morphological damage, such as the loss of cell structural integrity, cell deformation, and breakage of cell walls and membranes.
DISCUSSION
Garlic-derived organosulfur compounds, such as thiosulfinates and their derivatives, have been found to have a strong antimicrobial effect on a variety of pathogenic bacteria (28, 61). In the current study, two stable organosulfur compounds derived from decomposed thiosulfinates were tested for their antibacterial effect on C. sakazakii. The purity and stability of ajoene and diallyl sulfide were monitored using reverse-phase high-performance liquid chromatography (HPLC) (24, 25) during the course of the study. Both organosulfur compounds were stable for at least a month when kept in the dark at 4°C without light exposure. Both ajoene and diallyl sulfide have shown significant antimicrobial activity against C. sakazakii, and the antibacterial effect of ajoene was significantly higher than that of diallyl sulfide (see Table S1 and Fig. S2 in the supplemental material). A previous study confirmed that Z-ajoene has higher antimicrobial activity than E-ajoene (19). The MIC value of ajoene for C. sakazakii in the current study was slightly higher than 100 μg/ml, which was similar to the MICs for Gram-negative bacteria reported by Yoshida et al. (22) but higher than the ones for Helicobacter pylori reported by Ohta et al. (19). Naganawa et al. reported that ajoene was more effective in inactivating Gram-positive bacteria than Gram-negative bacteria (26). The antibacterial effects of ajoene on C. sakazakii cells in two different matrices (i.e., broth and PBS) and at different temperatures (i.e., 4, 22, and 37°C) were significantly (P < 0.05) different (see Table S1). This may have been due to the differences in solubility of ajoene in different solvents and at different temperatures affecting the antimicrobial activity.
Diallyl sulfide is present in garlic oil at a concentration ∼112 μg/g (23). In the current study, the MIC value of diallyl sulfide for C. sakazakii was about 500 μg/ml. Ross et al. (62) and O'Gara et al. (63) reported the significant antimicrobial effect of diallyl constituents in garlic oil on a variety of bacteria and also confirmed that the antimicrobial activity of the diallyl sulfides increased with the number of sulfur atoms, which was also recently validated by our group (24). A recent study has shown that the heat-stable components of garlic extract play a major role in the antimicrobial activity of garlic extract and that both diallyl sulfide and ajoene are thermally stable (23). Those authors have confirmed that diallyl sulfide could significantly inhibit the growth of Aggregatibacter actinomycetemcomitans even at a low concentration (i.e., 0.01 μg/ml). Allicin in the garlic extract lost all biological activity after heat treatment, but diallyl sulfide showed no significant difference of the antimicrobial effects before and after heat treatment.
When l-cysteine was added to C. sakazakii cultures with ajoene, the antimicrobial effect of ajoene was abolished rapidly. In contrast, l-cysteine did not affect the antibacterial effect of diallyl sulfide on C. sakazakii (see Table S2 in the supplemental material). This result indicated that the effective functional groups for the antimicrobial activity of ajoene and diallyl sulfide are different. Although diallyl sulfide and ajoene are similar in structure (e.g., consisting of allyl groups), the results suggest that the disulfide bond is important for the antimicrobial activity of ajoene. It used to be believed that thiosulfinates (i.e., ajoene) inhibit microorganisms because of their disulfide bond, which reacts generally with the −SH groups of cellular proteins. For example, Feldberg et al. validated that the in vitro antimicrobial mechanism of allicin was associated with immediate and total inhibition of RNA synthesis and general inhibition of −SH proteins/enzymes (56). When ajoene and l-cysteine exist simultaneously at neutral pH, they react quickly and produce a 1 molar equivalent of S-allyl mercaptocysteine (S-AMC) and E,Z-2-amino-9-oxo-4,5,9-trithaidodeca-6,11-dienoic acid (ajocysteine) (see Fig. S3) (64). This clearly highlighted that the vast majority of antimicrobial properties of ajoene were rapidly eliminated but that diallyl sulfide was not when l-cysteine was present. Cavallito et al. also reported that allicin could react with a −SH group of cysteine and forms S-allylmercaptocysteine (allyl-ss-cysteine) by cleaving the disulfide bond and replacing the allyl group with cysteine (65). l-Cysteine is a critical part of glutathione S-transferase (GST), which is one of the most common mechanisms that bacteria use to antagonize oxidative stresses (23). Jakobsen et al. (25) have found that synthesized ajoene loses activity in in vivo infectious models but that fresh garlic extract showed a much more pronounced effect, which may have been due to the presence of l-cysteine or GST.
Ajoene also contains a sulfinyl group (see Fig. S1B in the supplemental material), which was reported for antibacterial agents such as allicin. Therefore, the antimicrobial activity of ajoene is derived from the presence of both the disulfide bond and the sulfinyl group (26). There was no clear evidence to validate which structures contribute to the antimicrobial activity of diallyl sulfide. The antimicrobial effects of diallyl sulfide in two different matrices and at different temperatures, especially the antimicrobial effects of diallyl sulfide in TSB, were significantly different from those of ajoene, confirming that there were significant differences in the antimicrobial mechanisms of ajoene and diallyl sulfide.
Previous work has revealed that thiosulfinates and diallyl sulfide can easily penetrate the bacterial cell membrane (66) and bind to thiol-containing proteins in bacterial cells and form disulfide bonds (24, 31). An interesting observation in the current study was that both ajoene and diallyl sulfide could bind to −SH in proteins in bacterial cells (described in detail below). These results seemed to be contradictory with respect to the fact that cysteine could react only with ajoene rather than with diallyl sulfide, resulting in the variations in the effects of the antimicrobial activity of ajoene and diallyl sulfide on bacteria with the addition of cysteine. A plausible explanation is that the size of the molecules attached to −SH affects the binding of this functional group with organosulfur compounds to form a disulfide stretching bond.
Here we have discussed in detail the mechanism of C. sakazakii stress and injury resulting from treatment with diallyl sulfide and ajoene (Fig. 3). Once the cells were treated with these garlic-derived organosulfur compounds, bacterial stress occurred.
The mechanism of bacterial stress and injury by the treatment of diallyl sulfide and ajoene.
Although both diallyl sulfide treatment and ajoene treatment resulted in moderate upregulation of gene expression, the profiles were very different. There were fewer upregulated genes in ajoene-treated bacteria (637) than in diallyl sulfide-treated bacteria (1,434). Ajoene-treated bacteria were also found to have no significantly enriched annotation clusters that were upregulated. However, treatment with ajoene resulted in a significant decrease in expression of flagellum-related genes. This may indicate a specific bacterial response to ajoene. Diallyl sulfide is known to cause cell wall disruption (18). Therefore, it was not surprising to find that many of the genes that were upregulated belonged to biological clusters that are involved in cell wall synthesis. The downregulation of housekeeping genes may imply a state of decreased metabolic activity in response to the stress. Taken the results together, ajoene and diallyl sulfide result in very different bacterial responses, suggesting different mechanisms for the garlic-derived compounds.
A critical observation in this study was that two distinct phases of interaction between garlic-derived organosulfur compounds (i.e., ajoene and/or diallyl sulfide) and thiol-containing proteins were illustrated (Fig. 3). By using a confocal Raman microscope, we observed that ajoene and diallyl sulfide could penetrate the bacterial cell membrane and then combine with the −SH group of cellular proteins in bacteria to form an S-S stretching bond (bands at 511 cm−1 and 675 cm−1 in Fig. 2C), indicated as phase I. This subsequently altered the conformational structures of the relative enzymes/proteins (bands at 1,638 cm−1 and 1,605 cm−1 in Fig. 2B), resulting in malfunction, indicated as phase II. Ajoene and diallyl sulfide showed similar phase I results but different phase II results during the interaction with thiol-containing proteins in bacterial cells. Ajoene mainly altered the secondary structural features of phenylalanine and tyrosine (band at 1,605 cm−1), while diallyl sulfide changed the structural features of amide I for both α-helix and β-sheet (band at 1,638 cm−1).
It should be emphasized that the strain whose genome was sequenced (i.e., C. sakazakii ATCC 29004) and used in this study has undergone a series of phenotypic changes over time and may not be completely representative of the wild type. Further investigation using several different C. sakazakii strains is required to further test the generality of the current hypotheses.
In the current study, we systematically investigated bacterial injury using various analytical tools. Both RNA-seq and Raman spectroscopy are high-throughput techniques to independently characterize bacterial stress and injury under unfavorable conditions. Moen et al. used an explorative multifactor approach for investigating global survival mechanisms of C. jejuni under unfavorable environmental conditions (67). In that study, data were generated with DNA microarrays for information about gene expression patterns and with Fourier transform infrared (FT-IR) spectroscopy to study global macromolecular changes in cell. However, FT-IR spectroscopy was applied only to the wave numbers of 1,200 cm−1 to 900 cm−1 to specifically determine the variations in the amounts of polysaccharides and oligosaccharides on C. jejuni cell membranes under the nongrowth survival conditions. In contrast, RNA-seq has a lower detection limit and can be used to analyze the stresses more completely than DNA microarrays. Raman spectroscopy can investigate all the major macromolecular variations in the bacterial cell, such as polysaccharides, proteins, and lipids (Fig. 2 and Table 3). Further, the confocal technique of Raman microscope allows us to clearly demonstrate the disulfide bond formation during the interaction between garlic-derived organosulfur compound and thiol-containing proteins in cells, which could not be monitored using FT-IR spectroscopy. Therefore, complementary RNA-seq and confocal micro-Raman spectroscopy are ideal to study bacterial stress and injury under conditions of unfavorable treatments, such as those using plant-derived bioactive compounds as novel antimicrobials.
In conclusion, we identified two garlic-derived organosulfur compounds (i.e., ajoene and diallyl sulfide) that have significant antimicrobial activity against C. sakazakii. RNA-seq was applied to determine the transcriptomic response of C. sakazakii treated by sublethal concentrations of ajoene and diallyl sulfide. The two compounds differed significantly in the patterns of induced gene expression. Confocal micro-Raman spectroscopy was used to fingerprint macromolecular variations in bacterial cells, and two distinct phases of interaction between organosulfur compounds and thiol-containing proteins in bacterial cells could be monitored. This research demonstrates the feasibility of combining different high-throughput analytical tools to study the mechanism of bacterial stress and injury under unfavorable conditions. Novel antimicrobial agents discovered from plants are potential resources to inactivate pathogens in the environment and on food contact surfaces to reduce food-borne diseases.
ACKNOWLEDGMENTS
This work was supported by funds awarded to X.L. by the UBC new faculty start-up grants; funds awarded to S.W. by the Ministry of Science and Technology of China (2012CB720803); and funds awarded to M.E.K. by the National Institutes of Health (R56 AI088518-01A1) and the United States Department of Agriculture, National Institute of Food and Agriculture (2011-67015-30772). T.P.E. was supported by the National Institute of General Medical Sciencesy (T32GM008336).
The content is solely our responsibility and does not necessarily represent the official views of NIH or USDA.
We thank Chunrong Liu and Ming Xian of the Department of Chemistry, Washington State University, for the instruction of synthesis of ajoene and Mingyuan Zhong at Applied Mathematics, University of Washington, for critical discussion about Bayesian's modeling for spectroscopic analysis applied in the current study.
FOOTNOTES
- Received 19 October 2013.
- Accepted 18 November 2013.
- Accepted manuscript posted online 22 November 2013.
Supplemental material for this article may be found at http://dx.doi.org/10.1128/AEM.03460-13.
- Copyright © 2014, American Society for Microbiology. All Rights Reserved.