Morgan L Wright1, Ken Pendarvis2, Bindu Nanduri3, Mariola J Edelmann4, Haley N Jenkins1, Joseph S Reddy3, Jessica G Wilson1, Xuan Ding5, Paul R Broadway6, Mais G Ammari2, Oindrila Paul1, Brandy Roberts1 and Janet R Donaldson1*#
Received Date: January 22, 2016; Accepted Date: March 29, 2016; Published Date: April 05, 2016
Citation: Wright ML, Pendarvis K, Nanduri B, Edelmann MJ, Jenkins HN, et al. (2016) The Effect of Oxygen on Bile Resistance in Listeria monocytogenes. J Proteomics Bioinform 9:107-119. doi:10.4172/jpb.1000396
Copyright: © 2016 Wright ML, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Proteomics & Bioinformatics
Listeria monocytogenes is a Gram-positive facultative anaerobe that is the causative agent of the disease listeriosis. The infectious ability of this bacterium is dependent upon resistance to stressors encountered within the gastrointestinal tract, including bile. Previous studies have indicated bile salt hydrolase activity increases under anaerobic conditions, suggesting anaerobic conditions influence stress responses. Therefore, the goal of this study was to determine if reduced oxygen availability increased bile resistance of L. monocytogenes. Four strains representing three serovars were evaluated for changes in viability and proteome expression following exposure to bile in aerobic or anaerobic conditions. Viability for F2365 (serovar 4b), EGD-e (serovar 1/2a), and 10403S (serovar 1/2a) increased following exposure to 10% porcine bile under anaerobic conditions (P < 0.05). However, HCC23 (serovar 4a) exhibited no difference (P > 0.05) in bile resistance between aerobic and anaerobic conditions, indicating that oxygen availability does not influence resistance in this strain. The proteomic analysis indicated F2365 and EGD-e had an increased expression of proteins associated with cell envelope and membrane bioenergetics under anaerobic conditions, including thioredoxin-disulfide reductase and cell division proteins. Interestingly, HCC23 had an increase in several dehydrogenases following exposure to bile under aerobic conditions, suggesting that the NADH:NAD+ is altered and may impact bile resistance. Variations were observed in the expression of the cell shape proteins between strains, which corresponded to morphological differences observed by scanning electron microscopy. These data indicate that oxygen availability influences bile resistance. Further research is needed to decipher how these changes in metabolism impact pathogenicity in vivo and also the impact that this has on susceptibility of a host to listeriosis.
ACN: Acetonitrile; BHI: Bran Heart Infusion; FA: Formic Acid; PBS: Phosphate Buffered Saline; SEM: Scanning Electron Microscope
Listeria monocytogenes is a Gram-positive facultative anaerobe and the causative agent of listeriosis . With a nearly 20% mortality rate in the United States, L. monocytogenes typically manifests disease in the young, elderly, pregnant women and the immunocompromised through the consumption of contaminated foods [2,3]. Upon consumption L. monocytogenes must resist the multiple stressors encountered within the gastrointestinal tract, including bile, variations in pH, and oxygen availability [4,5]. Bile is specific to the gastrointestinal tract and consists of multiple components such as ions, cholesterol, proteins, bile salts, and pigments . Of these, bile salts have been shown to possess antimicrobial activity through the induction of DNA damage and degradation of viral and bacterial membranes [7,8]. Listeria monocytogenes is able to grow in the gall bladder , indicating this bacterium is resistant to the highest concentration of bile salts encountered within the body.
Bile resistance mechanisms have been extensively studied, including the bile salt hydrolase bsh [9,10], the general stress response sigma factor sigB [11,12], the bile exclusion system bilE , and virulence regulator prfA . However, information is lacking in regards to the response of L. monocytogenes to bile under physiologically relevant anaerobic and microaerophilic conditions. The expression of genes needed for survival under acidic conditions was found to increase under anaerobic conditions [14,15] and an increase in branch-chain fatty acids in the cell membrane was observed when L. monocytogenes was cultured under elevated carbon dioxide and anaerobic conditions . Additionally, the activity of the bile salt hydrolase has been found to increase under anaerobic conditions . Oxygen restriction also enhances growth at lower temperatures (~ 19°C). Together, these data suggest that oxygen availability influences resistance to stressors, which could potentially impact the virulence capability of L. monocytogenes.
Although much research has been conducted to identify and characterize proliferation and pathogenesis, few studies have analyzed the effect of oxygen availability on the bile resistance properties of L. monocytogenes including the influence of oxygen on regulation of the proteome expressed in response to bile. Therefore, this study focused on comparing the influence that reduced oxygen has on bile resistance in strains representing three serovars of L. monocytogenes. Two of these serovars (1/2a and 4b) represent nearly 90% of all listeriosis cases, whereas serovar 4a is rarely associated with listeriosis . Here, we report that anaerobiosis increased bile resistance with the serovar 1/2a and 4b strains tested and the resistance may be linked to differential expression patterns of metabolic and membrane bioenergetic proteins.
Bacterial strains and culture conditions
Strains used in this study were F2365 (serotype 4b), EGD-e (serotype 1/2a), 10403S (serotype 1/2a), and HCC23 (serotype 4a) [17-19]. All strains of L. monocytogenes were cultured in Brain Heart Infusion (BHI) broth and incubated in a shaker incubator at 37°C at 250 rpm.
Anaerobic conditions: Overnight cultures were diluted 1:100 into 10 mL of fresh BHI in an anaerobic chamber (Coy Laboratories) with a gas mix of 95% N2, 5% H2 in Wheaton serum vials capped with rubber stoppers and sealed with aluminum; the redox indicator resazurin (5 μM) was added in order to visually monitor anaerobiosis. Inoculated vials were grown to mid-logarithmic phase (OD600 ~ 0.4) at 37ºC with agitation, at which time cells were divided into 4 separate 2 mL aliquots, pelleted immediately at 8,000 × g for 5 min, then resuspended in 2 mL of BHI supplemented with 0%, 1% (0.02 g), 5% (0.1 g), or 10% (0.2 g) porcine bile extract (Sigma B8631, Sigma Aldrich) and 0.1% methanol (for solubility of bile) in the anaerobic chamber. Samples (100 μL) were collected using a syringe needle at 0, 1, 2, 3, 4, 5 and 6 h post exposure to porcine bile extract. Samples were serially diluted in phosphate buffered saline (PBS) and plated on BHI agar. Plates were incubated under anaerobic conditions using an AnaeroPack System (Mitsubishi Gas Chemical) at 37ºC for 18 h prior to viable plate count analysis.
Aerobic conditions: Bacterial cultivation was performed as described for the anaerobic cultivation assay, but without the use of sealed vials or the addition of resazurin. All cultures were incubated under normal atmospheric conditions at 37ºC. A minimum of three independent experiments was performed for each strain under each condition tested.
Protein sample preparation
Cultures were grown to mid-log in either aerobic or anaerobic conditions or subsequently exposed to either 0% or 5% porcine bile for 1 h as described for the survival assays. At 1 h post exposure, 10 mL were collected and immediately centrifuged at 8,000 × g for 10 min at 10ºC. Cell pellets were resuspended in 700 μL NP-40 lysis buffer (150 mM NaCl, 1.0% IGEPAL, 50 mM Tris, pH 8.0) supplemented with protease inhibitors (Sigma P2714, Sigma Aldrich) and lysed for 60 sec through sonication using Covaris model S220. Lysate was plated to ensure complete lysis. Lysates were centrifuged at 20,000 × g for 15 min to remove cell debris. A Fluka Protein Quantification Kit was used to determine protein concentration using the manufacturer’s protocol (Sigma Aldrich). Proteins (100 μg - 200 μg) were then purified by chloroform-methanol precipitation in order to remove any residual detergent, which would have interfered with mass spectrometry analysis. Briefly, samples were mixed with 600 μL of methanol, 150 μL of chloroform, and 450 μL of sterile deionized water, vortexed, and centrifuged at 13,000 × g for one min. Supernatant was removed, 450 μL of methanol was added to the interface layer, and centrifuged again at 13,000 × g for two min. Protein pellets were dried in a vacuum centrifuge at room temperature, denatured and alkylated, and digested with 2 μg of trypsin overnight as previously described  and desalted using a macrotrap (Michrome Bioresources). Cleaned samples were resuspended in 20 μL of 0.1% formic acid (FA) and 5% acetonitrile (ACN) for proteomic analysis.
Proteomic analysis by mass spectrometry
Peptides were analyzed using a Dionex UltiMate 3000 (Thermo Scientific) high performance liquid chromatography machine (HPLC) coupled with an LTQ-OrbiTrap Velos (Thermo Scientific) tandem mass spectrometer. The linear trap of the LTQ was used for precursor and fragment scans, as the resolution of the OrbiTrap was not required for these analyses . The HPLC was configured for reverse phase chromatography using a C18 Acclaim PepMap RSLC column (Thermo Scientific) with a flow rate of 300 nL/min. Peptides were separated for mass spectrometry analysis using an acetonitrile gradient starting at 2% ACN, 0.1% FA and reaching 50% ACN, 0.1% FA in 120 min, followed by a 15 min wash of 95% ACN, 0.1% FA. Column equilibration was handled automatically using the Dionex UltiMate 3000. The eluate from the HPLC was fed directly to the LTQ for nanospray ionization followed by MS / MS analysis of detected peptides by Collision-induced dissociation (CID). The LTQ was configured to perform 1 ms scan followed by 17 MS / MS scans of the 17 most intense peaks repeatedly over the 135 min duration of each HPLC run. Dynamic exclusion was enabled with duration of 3 min, repeat count of 3, and a list length of 500. For fragment mass analysis, activation time was 40 ms and normalized collision energy was 35. Ion trap mass spectrometer (IT-MS) was used for both full mass analysis, as well as fragment mass analysis. Raw files from the LTQ were converted to mgf format using the MSConvert GUI software from the ProteoWizard toolset [21,22]. L. monocytogenes strain-specific protein FASTA databases for were downloaded from the National Center for Biotechnology Information reference sequences (NCBI RefSeq 10403S gi|386042347; EGD-e gi|16802048; F2365 gi|85700163; and HCC23 gi|217963303). The X!tandem [PMID: 14976030] and OMSSA algorithms [PMID: 15473683] were used to match ms/ms spectra to the FASTA databases. Precursor and fragment mass tolerances were set to 1 Da and 0.5 Da respectively for X!tandem. Both precursor and fragment mass tolerances were set to 1 Da for OMSSA. Tryptic cleavage rules were used when calculating in-silico peptide precursor masses; b and y ions were used for fragment m/z matching. Amino acid modifications that were included in the database searches were single and double oxidation of methionine and both carboxymethylation and carbamidomethylation of cysteine, phosphorylation of serine, threonine and tyrosine, as well as water loss from serine and threonine. Randomized versions of each protein FASTA were concatenated to the originals as a means to calculate false discovery rates (FDR). Peptide-spectrum matches with e-values < 0.05 were accepted for down-stream analysis. Peptide matches were organized by protein using Perl, at which time non-unique peptide sequences and proteins identified by a single peptide sequence were removed. The maximum acceptable FDR was 1%; if the FDR was above 1%, the peptide e-value cutoff was lowered incrementally by 0.005 until the FDR fell below 1%. After each iteration, proteins reduced to a single peptide sequence were removed from the results. The mass spectrometry proteomics data has been deposited to the ProteomeXchange Consortium  via the PRIDE partner repository with the dataset identifier PXD002243 and 10.6019/PXD002243.
Differential expression of proteins between porcine bile treated and non-treated samples as well as between aerobic and anaerobic conditions upon bile treatment was performed pairwise using peptide elution profiles. Precursor mass spectra were extracted from the raw data in MS1 format using the MSConvert GUI software from the ProteoWizard toolset [21,22]. Peptide precursor m/z values were extracted from the previously compiled protein identifications using Perl. Elution profiles for peptide-spectrum matches were calculated by parsing each corresponding MS1 file and summing the ion current for that match’s m/z value within a 0.25 Da tolerance, effectively integrating the elution profiles. Each trace started at the scan number of the peptide-spectrum match and preceded both forward and backward until the chromatogram noise level, or a distance of 250 scans, was reached. Multiple peptide-spectrum matches with the same precursor m/z were only counted once, ensuring the same integral was not included multiple times. Once all peptide-spectrum matches were processed, intensities were summed for each protein on a per-replicate basis. Proteins not identified in a replicate were represented with the average noise level of the replicate’s chromatogram for further calculations. The reasoning behind this is two-fold: 1) peptides not identified in a replicate could be present at levels at or below the noise level of the chromatogram, causing the mass spectrometer to ignore them and 2) for calculating expression ratios between lines, zero cannot be in the denominator. Data were normalized using a mode-based technique. First, the mode of the protein intensities for each replicate was calculated, representing the most commonly occurring protein intensity. Next, for each identified protein, the intensity per replicate was divided by the mode of the same replicate. This ensures that normalization is not affected by the minimum and maximum intensities, which can vary tremendously between replicates. A permutation analysis was performed for each protein by evaluating the difference in means of the replicates of both conditions. From this permutation, a p-value was calculated to indicate the significance of the difference in means. Two additional permutations were performed for each protein, comparing both conditions to their own baselines. These baseline permutations provided a mechanism to further reduce false positives introduced by differences in chromatogram ion current as electron multiplier performance deceases. Proteins were considered to be differentially expressed if the difference in means between conditions resulted in a P < 0.05 and the difference in means between one of the conditions and its baseline was P < 0.05.
Protein datasets identified as significantly differentially expressed were subjected to KOBAS 2.0 analysis (http://kobas.cbi.pku.edu.cn/home.do). KOBAS 2.0 was applied to first annotate all of the entries with Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and to then identify significantly enriched pathways with P < 0.05.
Gene expression assays
All strains were cultured in 5 mL BHI under aerobic or anaerobic conditions and exposed to either 0% or 5% porcine bile extract for 1 h. Cells were pelleted by centrifugation at 10,000 × g for 2 min, washed in ice cold PBS, then treated with RNAprotect Bacteria Reagent (Qiagen) according to the manufacturer’s protocol. RNA was isolated using the RNeasy mini kit (Qiagen) following an altered manufacturer’s protocol as followed. Briefly, 400 μL of RLT buffer with 2-β-mercaptoethanol was added to cell pellets, followed by lysis using a bead beater for two 2 min intervals. Ethanol was then directly added to lysed cells at a 60% final concentration to the samples prior to being homogenized using a QIAshredder (Qiagen). Samples were treated with RNase free DNase. RNA and DNA quantitation was determined using Qubit RNA BR assay and DNA HS assay, respectively, and analyzed using a Qubit 2.0 fluorometer (Life Technologies) following manufacturer’s protocol. RNA was normalized to 100ng and converted to cDNA using an Applied Biosystems High Capacity cDNA Reverse Transcriptase Kit following manufacturer’s protocol. The RT-PCR protocol was as follows: 25°C for 10 min, 37°C for 120 min and 85°C for 5 sec. Concentrations were determined using a Nanodrop ND-100.
Expression of bsh was determined in relation to expression of 16S rRNA gene using an Applied Biosystems Step One Plus System. For the bsh target, the forward primer was 5’- CCTGTTGGCGTGTTAACAAATAA-3’, the reverse primer was 5’- CCATCCCACGACTATAAGCATC -3’, and the probe was 5’-FAM-TCGCGTTCT/ZEN/TTCGAGTGAAACTCCA-IowaBlackFQ. For the 16S target, the forward primer was 5’-CTTGTCCCTTGACGGTATCTAAC-3’ the reverse primer was 5’-GCGCTTTACGCCCAATAAATC-3’ and the probe was 5’-FAM/CGGTAATAC/ZEN/GTAGGTGGCAAGCGT-IowaBlackFQ. Standard curves were generated to verify primer efficiency. Each reaction received 10 μL of 2x Taqman Gene Expression Master Mix (Life Technologies), 2 μL of 10x PrimeTime qPCR assay mix designed for either bsh or 16S rRNA (Integrated DNA Technologies) and 5 μL of 1:100 diluted cDNA template, with volume adjusted to 20 μL with RNAse free water. The qPCR reaction followed two stages; first stage: 50°C for 2 min and 95°C for 10min; second stage: 95°C for 15 s and 60°C for 1 min, for 40 cycles. Fold changes in expression of bsh were calculated based on expression levels of 16S rRNA from three independent experiments as previously described .
Scanning electron microscopy
F2365, EGD-e, 10403S, and HCC23 were cultured aerobically or anaerobically in 5 mL BHI to mid-log phase as described above, at which time cultures were split to two 2 mL aliquots and treated with either 0% or 5% bile for 1 h at 37°C. Cells were then processed for scanning electron microscopy (SEM) as previously described . Briefly, cells were pelleted by centrifugation at 8,000 × g for 5 min, fixed in 2.5% glutaraldehyde in 0.1M cacodylate buffer, washed in 0.1 M cacodylate buffer, post-fixed in 2% osmium tetraoxide in 0.1 M cacodylate buffer, dehydrated in an ethanol series, and dried in hexamethyldisilazane series. Samples were sputter coated with platinum prior to observation on a JOEL 6500F SEM. Cell size was determined using the ImageJ software.
Survival assays were analyzed using the Proc Glimmix procedure of SAS (v. 9.4, SAS Institute, Cary, NC). Fixed effects included bacterial strain, treatment, oxygen availability, time, and interactions between these parameters. All concentrations were converted to Log10 to achieve normality. Means were separated using LSMeans with an alpha of 0.05 using Tukey-Kramer adjustment. For SEM analyses, the length and width of both control and treatment samples under either aerobic or anaerobic conditions were analyzed using a T-tailed type II T-test, with P < 0.05 deemed as significant.
Oxygen availability influences the survival of L. monocytogenes in bile
The aim of this study was to determine if oxygen availability influenced bile resistance among strains of L. monocytogenes. The viability of strains 10403S (serovar 1/2a), EGD-e (serovar 1/2a), F2365 (serovar 4b) and HCC23 (serovar 4a) following exposure to porcine bile extract was first assessed under both aerobic and anaerobic conditions (Figure 1). Exposure to 1% or 5% bile under aerobic conditions did not impact viability of 10403S. The 10% concentration resulted in a 2.4 Log10 decrease in viability within the first hour post exposure (P = 0.001). EGD-e decreased in viability following exposure to 5% and 10% bile, with a nearly 1 Log10 and 3 Log10 reductions respectively (P < 0.01 and P < 0.001, respectively). F2365 was able to continue replicating in media supplemented with either 1% or 5% porcine bile extract, but not 10%. Exposure to 10% bile resulted in a nearly 3 Log10 decrease in viability within the first hour post exposure (P < 0.001), which was similar to what was observed in the other strains analyzed (P < 0.05). Interestingly, the avirulent strain HCC23 had a slight increase in viability immediately following exposure to 1% bile, but then remained static for the remainder of the times analyzed. Exposure to 5% and 10% resulted in a decrease of 1.2 Log10 and 2.7 Log10 within the first hour, respectively (P < 0.001).
The virulent strain F2365 was able to survive and continue to replicate following exposure to all concentrations of bile tested under anaerobic conditions (Figure 2). A slight increase of 0.2 Log10 was observed within 1 h exposure to 1% bile (P = 0.03). Also, an increase of 0.2 Log10 after 1 h exposure to 5% and 10% bile was observed in F2365 under anaerobic conditions.
An increase of approximately 1 Log10 was observed after 1 h exposure to 1% bile in both 10403S and EGD-e (P = 0.01). The viability of 10403S increased until 3 h, at which point no changes in viability were noted (Figure 2). In 10% bile, EGD-e decreased by approximately 1 Log10 within 1 h post exposure (P = 0.02) and recovered at 2 h. Following an initial decrease in viable colonies, the avirulent strain HCC23 had relatively no change in viability in all bile concentrations analyzed, except for a significant decrease in viability following 6 h of exposure to 10% bile HCC23 (4.69 Log10, P < 0.001).
Statistical comparisons of the viability of L. monocytogenes exposed to bile under aerobic versus anaerobic conditions revealed a difference in F2365 in 10% bile for each time point tested (P < 0.001; Table 1); exposure to 5% bile only showed significance at 1 h (P = 0.017). 10403S also had an increase in survival under anaerobic conditions in 10% bile for up to 4 h post exposure (P < 0.05). Interestingly, no significant change was observed in HCC23 upon exposure to bile in relation to oxygen availability except for 1% at 1 h (P = 0.011) and 2 h (P = 0.002).
|0% Bile||1% Bile||5% Bile||10% Bile|
Comparisons between aerobic and anaerobic conditions for each strain at each time point revealed significant survival difference between both conditions under various concentrations of bile salt stress
Table 1: Statistical analysis between growth in aerobic and anaerobic conditions.
Bile salt hydrolase activity varies between strains under aerobic and anaerobic conditions
The expression of bsh has been shown to increase under anaerobic conditions . To determine whether the differences in the bile survival between the four strains tested was due to variations in the expression of bsh, real-time PCR was utilized to quantify the expression following a 1 h exposure to bile under aerobic or anaerobic conditions (Table 2). Expression of bsh increased for 10403S, EGD-e, and HCC23 under anaerobic conditions (P < 0.001). The expression also increased slightly for F2365 under anaerobic conditions, but the change was not significant (P > 0.05). This suggests that the expression of bsh alone does not impact the increased resistance observed under anaerobic conditions.
|Aerobic Fold Change bsh (±StDev)||Anaerobic Fold Change bsh (±StDev)||P-value (aerobic v. anaerobic)|
|10403S||0.57 (0.05)||4.19 (0.42)||<0.001|
|EGD-e||1.50 (0.02)||4.38 (0.39)||<0.001|
|F2365||1.22 (0.38)||1.62 (0.15)||0.175|
|HCC23||1.78 (0.39)||2.65 (0.08)||<0.001|
Table 2: Fold changes in expression of bsh in aerobic and anaerobic conditions.
Proteomes vary between strains in response to bile
Since variations were observed in the viability of the four strains tested that could not be solely attributed to the variations in bsh expression, the proteomes of F2365, 10403S, EGD-e, and HCC23 were analyzed following exposure to bile in either aerobic or anaerobic conditions. Strains were treated for 1 h with either 0% or 5% porcine bile extract in aerobic or anaerobic conditions and proteins were isolated from whole cell lysates. Table 3 represents a summary of the number of proteins identified as significantly differentially expressed. Proteomics data quality was assessed by including a randomized version of the protein FASTA database when performing spectrum matching. The maximum FDR reported from these analyses was 0.993%, with a mean FDR of 0.375%.
|Aerobic proteins||Anaerobic proteins|
|Strain||Total identified||1h sig different||1h bile sig different||Total identified||1h sig different||1h bile sig different|
Table 3: Summary of proteins identified by mass spectrometry.
Proteins identified were analyzed by KOBAS to identify whether trends were evident in expression of specific pathways. Table 4 represents a summary of the pathways that were significantly differentially expressed following exposure to bile under either aerobic or anaerobic pathways. Alterations in purine metabolism were identified in all four strains analyzed. 10403S, F2365, and HCC23 all had an alteration in pyruvate metabolism under anaerobic conditions in the presence of bile; EGD-e only had an alteration in this pathway under aerobic conditions.
|KEGG pathway name||Background number||1h sig different||1h bile sig different||1h sig different||1h bile sig different|
|na||P value||na||P value||na||P value||na||P value|
|Nucleotide excision repair||7||3||0.043||-||-||-||-||-||-|
|Citrate cycle (TCA cycle)||11||-||-||10||0.031||-||-||-||-|
|Fatty acid biosynthesis||14||-||-||-||-||-||-||9||0.030|
|Fatty acid metabolism||15||-||-||-||-||-||-||10||0.019|
|Glycine, serine and threonine metabolism||31||6||0.045||-||-||6||0.022||-||-|
|Fatty acid metabolism||15||-||-||-||-||-||-||7||0.033|
|Glyoxylate and dicarboxylate metabolism||10||-||-||4||0.041||-||-||-||-|
|Nucleotide excision repair||7||-||-||5||0.004||-||-||-||-|
|One carbon pool by folate||10||-||-||-||-||4||0.020||6||0.021|
|Valine, leucine and isoleucine degradation||10||-||-||4||0.041||4||0.020||-||-|
aNumber of proteins in the input file associated with each pathway
Table 4: Significant pathways identified using KOBAS in significantly differentially expressed protein datasets.
Proteins associated with cell envelope and cellular processes are differentially expressed under aerobic and anaerobic conditions upon exposure to bile
Proteomic analyses revealed a change in the expression of proteins involved in redox reactions, invasion, and cell division following bile exposure (Table 5). The MreB protein, which is involved in determining cell shape, decreased upon bile treatment under aerobic and anaerobic conditions for virulent strains 10403S, EGD-e, and F2365 (Table 5). Expression of MreC only decreased following bile exposure under anaerobic conditions; expression increased under aerobic conditions. A decrease in several cell division proteins was detected under both aerobic and anaerobic conditions following a 1h exposure to bile. FtsZ, FtsH, and FtsA decreased in all strains tested following exposure to bile under both aerobic and anaerobic conditions.
|Protein||ListiList||0%||5%||0%||5%||0% vs. 5%||0% vs.5%|
|Peptide/nickel transport, ATP binding||1.2||0.001||0.976||6.666||3.459||Up||Down|
|PTS mannose-specific, factor IIAB||1.2||1.589||0.168||0.904||0.001||Down||Down|
|Cell division protein FtsZ||1.7||1.558||0.293||3.718||0.259||Down||Down|
|ATP synthase, epsilon chain||1.4||0.001||0.585||0.840||0.001||Up||Down|
|Cell division protein FtsZ||1.7||2.090||0.001||2.980||0.001||Down||Down|
|PTS mannose-specific, factor IIAB||1.2||1.798||0.001||2.011||0.001||Down||Down|
|Cell wall surface anchor family protein||1.8||0.003||1.802||0.161||0.001||Up||Down|
|Manganese ABC transporter protein||1.2||0.001||1.550||1.659||0.673||Up||Down|
|PTS mannose-specific, factor IIAB||1.2||0.769||0.328||2.446||0.001||Down||Down|
|ATP synthase epsilon chain||1.4||0.001||0.594||-||-||Up||-|
Table 5: Select proteins associated with the cell envelope with a significant change in expression following bile exposure under aerobic and anaerobic conditions.
Membrane bioenergetics proteins were also differentially expressed. Thioredoxin reductase decreased following exposure to bile under both aerobic and anaerobic conditions for all strains tested (Table 5). NADH dehydrogenase decreased following exposure to bile in both aerobic and anaerobic conditions for 10403S and HCC23. Several proteins related to invasion were differentially expressed following bile treatment under aerobic and anaerobic conditions. Flavocytochrome c, which is involved in intracellular survival , increased after exposure to bile in HCC23 under aerobic conditions.
Additional cell envelope associated proteins detected that are worth noting include several flagellar associated proteins and cell wall anchor proteins. Expression of the flagellar motor protein and flagellin increased in HCC23 under aerobic conditions following exposure to bile. Flagellin also increased in F2365 under aerobic conditions following bile treatment. A cell wall anchor protein increased in expression under aerobic conditions only in F2365. Lipoproteins were increased in expression in HCC23 following exposure to bile under both aerobic and anaerobic conditions. However, the lipoprotein adhesin B increased following exposure to bile only under aerobic conditions in HCC23.
Morphological changes occur following exposure to bile under aerobic, but not anaerobic conditions
Proteomic analysis revealed changes in expression of cell shape determining proteins and proteins involved in cell division, including FtsZ. Therefore, scanning electron microscopy (SEM) was utilized to observe the morphological changes that occurred between aerobic and anaerobic bile treatment conditions (Figure 3). Under aerobic conditions, an increase in cell length and width was observed in 10403S, EGD-e, and F2355 (P < 0.001) following bile exposure. An increase in elongation of HCC23 was seen under aerobic conditions as well (P = 0.04), which correlated with the increase in MreC that was observed by the proteomic analysis. In contrast, under anaerobic conditions, the only change in length was observed in EGD-e (P = 0.003). Changes in width were only observed in 10403S and HCC23 (P < 0.001).
Proteins associated with metabolism are differentially expressed under aerobic and anaerobic conditions upon exposure to bile
Many proteins involved in metabolism were differentially expressed following bile treatment under aerobic and anaerobic conditions (Table 6). Enzymes involved in the entrance into the TCA cycle were differentially expressed between aerobic and anaerobic conditions. For instance, pyruvate dehydrogenase, which is involved in the decarboxylation of pyruvate into acetyl CoA, decreased following bile exposure in both 10403S and F2365. However, expression of pyruvate dehydrogenase increased in HCC23 following bile exposure in aerobic conditions. The expression of pyruvate carboxylase, which converts pyruvate to oxaloacetate, decreased following bile exposure in aerobic conditions for HCC23, 10403S, and F2365. Tagatose-diphosphate aldolase, which is involved in galactose metabolism, decreased in expression for 10403S, F2365, and HCC23 under aerobic and anaerobic conditions. Several dehydrogenases were increased in expression in HCC23 following exposure to bile under aerobic conditions, including aspartate dehydrogenase, alanine dehydrogenase, glutamate dehydrogenase, and naloxone dehydrogenase.
|Protein||ListiList||0%||5%||0%||5%||0% vs.5%||0% vs.5%|
|Ribulose P 3- epimerase||2.1||0.456||0.992||2.617||0.001||Up||Down|
|Tagatose 1,6-P aldolase||2.1||2.264||0.805||1.719||0.001||Down||Down|
|Purine nucleoside phosphorylase||2.3||1.866||2.915||2.199||0.459||Up||Down|
|Tagatose 1,6-P aldolase||2.1||2.186||0.001||2.419||0.001||Down||Down|
|Propanediol utilization: dioldehydratase||2.2||0.001||0.652||-||-||Up||-|
|Purine nucleoside phosphorylase||2.3||0.987||0.134||1.106||0.001||Down||Down|
|Tagatose 1,6-P aldolase||2.1||3.001||0.475||3.588||0.001||Down||Down|
|Purine nucleoside phosphorylase||2.3||1.904||2.534||2.271||1.248||Up||Down|
Table 6: Select proteins associated with intermediary metabolism with a significant change in expression following bile exposure under aerobic and anaerobic conditions.
The osmotic stress response protein alanine dehydrogenase had altered expression following bile exposure in either aerobic or anaerobic conditions. Expression increased in HCC23 and 10403S following exposure to bile under aerobic conditions. Interestingly, expression decreased in F2365 and EGD-e following bile exposure under aerobic conditions. Metabolism protein GuaB, an inosine-monophosphate dehydrogenase associated with protection from stress through protein folding, decreased in all strains tested under both aerobic and anaerobic conditions.
Proteins associated with information pathways are differentially expressed under aerobic and anaerobic conditions upon exposure to bile
There were several significant differences in the expression of proteins related to stress responses and repair mechanisms between aerobic and anaerobic conditions (Table 7). Expression of the chaperone proteins DnaK and DnaJ decreased in HCC23, EGD-e, F2365, and 10403S following bile exposure in both aerobic and anaerobic conditions. Multiple proteins involved in the Clp operon, which is transcribed under stress conditions, were found to decrease in expression under both aerobic and anaerobic conditions following bile exposure in EGD-e, 10403S, and F2365. Interestingly, ClpP, which is involved in intracellular replication, increased in expression in HCC23 following bile exposure in aerobic conditions .
|Protein||ListiList||0%||5%||0%||5%||0% vs.5%||0% vs.5%|
|GTP-sensing repressor CodY||3.5||1.718||0.508||-||-||Down||-|
|GTP-sensing repressor CodY||3.5||1.882||0.704||2.204||0.001||Down||Down|
|GTP-sensing repressor CodY||3.5||2.077||0.001||1.951||0.001||Down||Down|
|GTP-sensing repressor CodY||3.5||1.486||0.547||1.979||0.001||Down||Down|
Table 7: Select proteins associated with information pathways with a significant change in expression following bile exposure under aerobic and anaerobic conditions.
The nucleotide excision repair protein UvrA decreased following bile exposure in EGD-e, 10403S, and F2365. Expression of RecA decreased following bile treatment in EGD-e, F2365, and HCC23 under both aerobic and anaerobic conditions. Expression of the recombination repair protein RecA increased in 10403S following bile exposure under aerobic conditions. The SigB negative regulator RsbW decreased in all strains tested following bile exposure in aerobic conditions.
OsmC, which is involved in the osmotic stress response, increased in HCC23 and 10403S under aerobic conditions following bile treatment, but not anaerobic conditions. Expression of the pore-forming toxin listeriolysin O decreased following exposure bile salts under anaerobic conditions for F2365 and EGD-e. Catalase, which is involved in degradation of hydrogen peroxide, was detected in all strains. However, expression in catalase decreased following exposure to bile under both aerobic and anaerobic conditions for all strains tested.
Bile salts, the bactericidal component of bile, induce oxidative damage in the DNA and membrane . Therefore, the ability of enteric pathogens to colonize the GI tract is dependent on their ability to survive in the presence of bile salts. Though many studies have elucidated the bile salt stress response of Listeria, information is lacking in regards to this response under physiologically relevant anaerobic conditions. A previous study conducted by our group analyzed the anaerobic bile response, but no information is available as to the comparative analysis between bile survival under aerobic and anaerobic conditions . Therefore, the purpose of this study was to determine if reduced oxygen affects the bile survival of L. monocytogenes.
In this study resistance of L. monocytogenes to various concentrations of bile salts was shown to be strain specific, as well as influenced by oxygen availability. A significant increase in bile resistance was observed for strains F2365, 10403S, and EGD-e under anaerobic conditions. Interestingly, HCC23 was the only strain tested that had an increase in bile resistance under aerobic conditions. This suggests that the mechanism by which L. monocytogenes responds to bile will differ depending upon oxygen availability in a strain dependent manner. This is of particular interest as the gastrointestinal tract is an anaerobic / microaerophilic environment.
Listeria monocytogenes possesses many proteins that assist in invasion and dispersal throughout the body . Even though these bacteria have the ability to evade stressors and the immune system by invading phagocytic and non-phagocytic cells, L. monocytogenes has been found to remain extracellular within the gallbladder, where the concentration of bile salts is the greatest . In the current study, lipoproteins involved in the invasion of cells increased upon exposure to porcine bile extract in HCC23, but not other strains tested. As lipoproteins are involved in cell invasion and intracellular survival , the differences in expression between strains suggests that the mechanism by which L. monocytogenes replicates within the lumen of the gallbladder is mediated by bile exposure . Additionally, exposure to bile has been found to increase biofilm formation . This finding supports previous data reported by our group that indicated variations in biofilm production occurs following exposure to bile  and provides additional support that oxygen availability may contribute to the ability of L. monocytogenes to invade.
Following bile salt stress, changes in morphology were observed through SEM analysis. Significant elongation and width of L. monocytogenes cells were seen following treatment. However, bile treatment under anaerobic conditions resulted in very little morphological changes in comparison to controls. This suggests that the mechanism utilized by L. monocytogenes to adequately respond to bile differs based on oxygen availability.
Previous studies have suggested a link between the phosphotransferase system (PTS) and PrfA (virulence regulator), in which PrfA-regulated genes are repressed by sugars such as mannose, glucose, and fructose that are transported through the PTS . Interestingly, PrfA has been shown to interfere with glucose uptake, resulting in a decrease in expression of PTS [29,30]. Expression of PTS associated proteins, including glucose and mannose-specific systems, decreased following bile exposure. Interestingly, HCC23 exhibited an increase in PrfA under aerobic conditions, but this was not observed in the other strains tested. This suggests that a link between carbon metabolism and PrfA expression may reflect the invasive capability of L. monocytogenes. Also, since proteins associated with the PTS system were detected under aerobic conditions, this could provide a way of reducing the activation of virulence genes until appropriate conditions are encountered.
Several studies have demonstrated the importance of repair mechanisms to overcome damage induced by bile salts. In Salmonella enterica, bile salts were found to induce oxidative DNA damage primarily in the form of transitions, specifically GC–AT . However, the type of damage induced under anaerobic conditions has not been characterized. To provide an assessment of the response to bile, the stress response proteins expressed by bile resistant strains 10403S, EGD-e, and F2365 and the bile sensitive strain HCC23 were analyzed. HCC23 had several stress response proteins that were differentially expressed in response to bile in comparison to the other strains tested. For instance, OsmC, which is an enzyme involved in osmotic and oxidative stress , was increased under aerobic conditions for HCC23 and 10403S and decreased in expression following bile treatment under anaerobic conditions. UvrA, which is associated with DNA repair, decreased following bile exposure in all strains tested. Catalase, which is needed for detoxifying reactive oxygen species that can induce oxidative damage, was detected in all strains tested. However, the expression of this protein decreased following exposure to bile regardless of whether exposure was aerobic or anaerobic. Superoxide dismutase was also detected; the expression increased in 10403S, though never to concentrations detected under aerobic conditions. Recombinational repair protein RecA was also detected in all strains, though expression only increased following bile treatment in 10403S. Together, the increased expression of DNA repair proteins under anaerobic conditions suggests oxygen influences the expression of stress response genes required for DNA repair. However, oxidative damage may be a secondary effect of bile as the increase in expression was not observed in all strains tested. The dependency of repair proteins in bile survival is most likely a multi-faceted mechanism that relies on strain specific mechanisms.
Phosphogluconate dehydrogenase is an enzyme involved in reducing NADP+ through the Pentose Phosphate Pathway. This protein was decreased in all strains analyzed following bile exposure, suggesting that bile shunts metabolism away from this pathway. Therefore, in order to recycle NAD+, the expression of additional dehydrogenases was assessed following bile exposure. HCC23 had an increase in expression of several dehydrogenases, following bile exposure under aerobic conditions. This was interesting, as F2365 and 10403S only had an increase in expression of aspartate dehydrogenase and alanine dehydrogenase, respectively. This could indicate an imbalance of NADH:NAD+ that could impact the membrane integrity and therefore damage induced by bile. Further research is needed to characterize these aspects of L. monocytogenes in relation to oxygen availability.
In summary, bile resistance for all strains studied was found to be strain dependent and this resistance is influenced by a reduction in available oxygen. Several changes in expression of proteins associated with the cell morphology, DNA repair, invasion, and metabolism were identified to be involved in the resistance under this stress. Proteins expressed under anaerobic conditions, including virulence factors and the SOS response, suggest that limited oxygen is needed for regulation of these proteins to overcome the damaging affects of bile salts. Oxygen availability influences bile resistance and this may be dependent on regulation of general stress responses by available oxygen. Further research is needed to characterize the mechanisms at which oxygen is detected and how it is used to regulate genes needed for survival under stress.
We would like to thank Ms. Lindsey Brown, Dr. Steven Ricke, and Dr. Justin Thornton for their helpful insights into this project. We would also like to thank Ms. Amanda Lawrence and Dr. Tibor Pechan for their technical assistance. This work was funded through the NIH funded Pathogen-Host Interactions COBRE at MSU, P20GM103646.
The authors have declared no conflict of interest.