Bacterial strains and growth conditions.
S. pneumoniae serotype 4 strain TIGR4 [47], ΔpotABCD [11], and the complement strain (pABG5-potABCD) [15] were used in this study. All strains were grown in either chemically defined medium (CDM) [48] or Todd- Hewitt broth supplemented with 0.5% yeast extract (THY) or on 5% sheep blood agar plates (BAP) in 5% CO2 at 37°C. All assays were performed in triplicate in three independent experiments.
RNA Sequencing.
Total RNA was isolated and purified from mid-log phase cultures (OD600 0.4) of TIGR4 and ΔpotABCD (n = 4) grown in THY (a complete medium that mimics nutrients in the host milieu) using the RNeasy® Mini Kit (Qiagen, Valencia, CA, USA). RNA quality was checked with an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). RNA-Seq analysis was performed as described earlier [15]. Briefly, libraries for RNA-Seq were prepared with the KAPA RNA Hyper Kit with RiboErase (KAPA Biosystem, Wilmington, MA, USA) with 5 µg RNA as input. The concentration and quality of libraries were determined by the Qubit ds DNA HS Assay Kit (Life Technologies, Carlsbad, CA, USA) and Agilent Tapestation (Agilent Technologies, Santa Clara, CA, USA). Sequencing was done on Illumina Hiseq 3000, the quality of the data was checked with Illumina SAV and de-multiplexing was performed with Illumina Bcl2fastq2 v 2.17. Removal of failed reads, mapping of the short sequence reads to the TIGR4 reference genome, and identification of differentially expressed genes were performed with CLC Genomic Workbench 20.0.3 (Qiagen, Valencia, CA, USA).
Paired end reads of both WT and ΔpotABCD were mapped to the TIGR4 genome using CLC proprietary read mapper and read counts were estimated by EM estimation algorithm [49] and DEGs were identified based on the fold change generated by the edgeR algorithm. Changes in gene expression with a fold change of ± 1.3 at a false discovery rate (FDR) of ≤ 0.05 were considered significant. Functions and pathways represented by DEGs were identified utilizing multiple bioinformatics resources such as MetaCyc [50], Gene Ontology [51], KEGG [52], UniProt [53], and STRING [54]. RNA-Seq raw data and metadata are available at NCBI GEO with the accession number XXXXXXX.
Quantitative real time PCR.
To validate RNA-Seq results, we measured expression of selected genes by quantitative reverse transcription PCR (qRT-PCR). The primers used for qRT-PCR are listed in supplementary material (Supp. Table 1). All primers were validated by performing a melt curve analysis with SYBR Green (Thermo Fisher Scientific Waltham, MA, USA). In brief, total RNA was purified from mid-log phase cultures (OD600 0.4) of TIGR4 and ΔpotABCD grown in THY (n = 3). Purified RNA (7.5 ng/reaction) was reverse-transcribed into cDNA and PCR was performed using the SuperScript III Platinum SYBR Green One-Step qRT-PCR Kit (Thermo Fisher Scientific, Waltham, MA, USA) as previously described [55]. Relative quantification of gene expression was determined by the Stratagene Mx3005P qPCR System (Agilent, Santa Clara, CA, USA). Expression of selected genes was normalized to the expression of gapDH and fold changes determined by the comparative CT method.
Measurement of intracellular pH.
The intracellular pH (pHi) was determined based on the method described by Clementi et al., [56] with slight modifications. Briefly, mid-log phase cultures (OD600 0.4) of TIGR4 and ΔpotABCD grown in THY (n = 3) were collected by centrifugation, washed, and suspended in phosphate buffered saline (PBS). Cells (108 CFU/mL) were loaded with 5 mM BCECF/AM dye (Millipore-Sigma, St. Louis, MO, USA) and incubated for 30 min at 30°C in the dark. Cells were then pelleted, washed, and reenergized with 10 mM glucose in PBS. To obtain the in vivo calibration curve for each strain, 400 µL of energized cells were pelleted and suspended in potassium buffers ranging from pH 6.5 to 8.0. Nigericin (1 mM) (Thermo Fisher Scientific, Waltham, MA, USA) was added to the cells (to equilibrate the pHi of the cells to the pH of the surrounding buffer) and incubated at 37°C for 5 min. Fluorescence was then measured by a Synergy H1 plate reader (BioTek, Winooski, VT, USA), and a calibration curve was obtained by plotting fluorescence against the pH of the buffers. To measure the pH of individual samples, 200 µL (108 CFU/mL) of the loaded and energized cells was added to the wells of a 96-well plate in duplicate and fluorescence detected using a plate reader for 5 min. 10 µM carbonyl cyanide 3-chlorophenylhydrazone (CCCP) was added to one well (to serve as control) and the reading taken for another 5 min. CCCP is a protonophore that uncouples proton motive force and causes a rapid decrease in pHi (Millipore-Sigma, St. Louis, MO, USA). Nigericin was added to both CCCP-treated control and untreated sample (to a final concentration of 1 mM) to equilibrate the pHi to the pH of the buffer and fluorescence was read for an additional 5 min. Fluorescence was calculated and the pHi was interpolated from the calibration curve.
Measurement of intracellular NADPH.
The intracellular concentration of NADPH was determined using the NADP/NADPH Assay Kit (Abcam, Cambridge, MA, USA). Mid-log phase cultures (n = 3) were harvested at 5,000 × g for 10 min at 4°C, suspended in PBS and transferred to beadbeater tubes (MP Biomedicals, Irvine, CA, USA). Cell suspensions were lysed with a FastPrep-24™ Classic benchtop homogenizer (MP Biomedicals, Irvine, CA, USA) and centrifuged at 6,000 × g for 5 min at 4°C. The cells were processed according to the manufacturer’s instructions. NADPH concentrations were determined with a SpectraMax® M5 multi-mode microplate reader (Molecular Devices, Sunnyvale, CA, USA). The concentration of the protein extracts was determined with the Pierce BCA Protein Assay Kit (Thermo Fisher Scientific, Waltham, MA, USA) and used to normalize NADPH concentrations.
Measurement of intracellular glutathione.
The ratios of reduced (GSH) to oxidized (GSSG) intracellular glutathione concentrations were determined using the GSH/GSSG-Glo™ Assay Kit (Promega, Madison, WI, USA). The cells were processed, and protein concentration determined as was done for NADPH quantification above. Luminescence was measured with a Cytation™ 5 cell imaging multi-mode reader (BioTek, Winooski, VT, USA) and used to calculate glutathione concentrations. The concentration of the protein extracts was determined with the Pierce BCA Protein Assay Kit (Thermo Fisher Scientific, Waltham, MA, USA) and used to normalize glutathione concentrations. GSH/GSSG (reduced/oxidized glutathione) ratios were calculated from the normalized glutathione concentrations according to the manufacturer’s instructions.
UPLC–HRMS untargeted metabolomics.
Approximately 109 CFU/mL of cells from mid-log phase cultures of TIGR4 and ΔpotABCD grown in THY (n = 5) were transferred onto a 0.2 µm Whatman polycarbonate membrane by vacuum filtration. The membranes were snap-frozen in liquid nitrogen and stored at − 80°C. Metabolites were extracted from bacteria on the membranes with extraction solvent (40:40:20 methanol, acetonitrile, and water with 0.1% formic acid) at 4°C (16). The extracts were transferred to 2.0 mL tubes, centrifuged for 5 min (16,100 × g) at 4°C and the supernatant transferred to new 2.0 mL tubes. Tubes containing ∼1.7 mL of the total supernatant were dried under a stream of N2 and solid residue was suspended in 300 µL of sterile water and transferred to autosampler vials for mass spectrometric analysis. A 10 µL aliquot was injected through a Synergi 2.5-micron reverse-phase Hydro-RP 100, 100 × 2.00 mM LC column (Phenomenex, Torrance, CA, USA) kept at 25°C. The eluent was introduced into the MS via an electrospray ionization source conjoined to an Exactive™ Plus Orbitrap Mass Spectrometer (Thermo Scientific, Waltham, MA, USA) through a 0.1 mm internal diameter fused silica capillary tube. The mass spectrometer was run in full scan mode with negative ionization mode with a window from 85 to 1000 m/z with a method adapted from [57]. Samples were run with a spray voltage of 3 kV. The nitrogen sheath gas was set to a flow rate of 10 psi with a capillary temperature of 320°C. Automatic gain control target was set to 3e6. The samples were analyzed with a resolution of 140,000 and a scan window of 85–800 m/z from 0 to 9 min and 110–1000 m/z from 9 to 25 min. Files generated by Xcalibur (RAW) were converted to the open source mzML format [49] via the open source msconvert software as part of the ProteoWizard package [49]. Maven (mzroll) software, Princeton University [58, 59] was used to automatically correct the total ion chromatograms based on the retention times for each sample [58, 59]. Metabolites were manually identified and integrated using known masses (± 5 ppm mass tolerance) and retention times (1 ≤ 1.5 min). Unknown peaks were automatically selected via Maven’s automated peak detection algorithms. A database of 275 metabolites verified using exact m/z and known retention times, expanded from the original database [57] was used. The statistical analysis on metabolite peak intensity post CFU normalization was done by MetaboAnalyst 4.0 [60]. Quantile normalization which is highly efficient in normalizing metabolite variations from mass spectrometry [61] was used to normalize the data. Significant differences in metabolite peak intensity between ΔpotABCD and TIGR4 were identified by a T-test at an adjusted FDR of ≤ 0.05. .
Hydrogen peroxide production.
H2O2 generated from mid-log phase cultures of TIGR4 and ΔpotABCD (n = 3) was compared using a quantitative peroxide assay (Pierce, Thermo Fisher Scientific Waltham, MA, USA). Briefly, 1 mL of bacterial culture (108 CFU/mL) grown in THY was centrifuged at 4 C for 2 min at 10,000 x g and the supernatant filtered with a 0.22 µm filter. The concentration of H2O2 was measured in the filtrate following the manufacturer’s instructions.
Hydrogen peroxide survival.
Mid-log phase cultures of TIGR4, ΔpotABCD and the complement pABG5-potABCD strain grown in either THY or CDM (OD600 0.4–0.5) were centrifuged at 10,000 x g for 2 min and suspended in PBS. The cells (108 CFU/mL) (in 1 mL PBS) were then supplemented with final concentrations of hydrogen peroxide 2.5 mM and incubated at 37°C with 5% CO2 for 15 min. CDM is devoid of polyamines but has the amino acid precursors for polyamine synthesis. To determine the effects of polyamines on pneumococcal H2O2 stress, ΔpotABCD challenged with 2.5 mM H2O2, was supplemented with cadaverine, putrescine or spermidine (½MIC, ¼MIC, ⅛MIC) and incubated for 15 min. MICs for cadaverine, spermidine or putrescine are 10.0, 4.8 and 5.7 mM, respectively. Control reactions had untreated bacteria, and CFUs were determined by serial dilution in PBS and plating on BAP. Results from three independent experiments were expressed as percent survival of treated bacteria relative to the untreated bacteria.
S -nitrosoglutathione susceptibility.
Mid-log phase cultures of TIGR4, ΔpotABCD and the complement pABG5-potABCD strain grown in CDM were centrifuged at 10,000 x g for 2 min and cells suspended in PBS. The cells (107 CFU/mL) in 100 µl were supplemented with a final concentration of 2.5 mM GSNO (Sigma-Aldrich, Israel), a nitric oxide producer, and incubated at 37°C in 5% CO2 for 60 min. In addition, ΔpotABCD challenged with 2.5 mM GSNO was supplemented with cadaverine, putrescine or spermidine (½MIC, ¼MIC, ⅛MIC). Control reactions had untreated bacteria, and CFUs were determined by serial dilution in PBS and plating on BAP every after 15 min. Results from three independent experiments were expressed as percent survival of treated bacteria relative to the untreated bacteria.
Statistical analysis.
Significant differences between the susceptibility of TIGR4, ΔpotABCD and the complement pABG5-potABCD strain to the different stressors, changes in pHi, production of endogenous H2O2, levels of NADPH and GSH/GSSG ratio, as well as changes in gene expression measured by qRT-PCR were determined by a Student’s T-test at p ≤ 0.05.