Previous Article | Next Article ![]()
Journal of Bacteriology, June 2006, p. 4068-4078, Vol. 188, No. 11
0021-9193/06/$08.00+0 doi:10.1128/JB.01921-05
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Zhili He,1,4
Eric J. Alm,1,2
Adam P. Arkin,1,2,7
Edward E. Baidoo,1,2
Sharon C. Borglin,1,3
Wenqiong Chen,1,6
Terry C. Hazen,1,3
Qiang He,1,4
Hoi-Ying Holman,1,3
Katherine Huang,1,2
Rick Huang,1,3
Dominique C. Joyner,1,3
Natalie Katz,1,3
Martin Keller,1,6
Paul Oeller,1,6
Alyssa Redding,1,7
Jun Sun,1,6
Judy Wall,1,5
Jing Wei,1,6
Zamin Yang,1,4
Huei-Che Yen,1,5
Jizhong Zhou,1,4 and
Jay D. Keasling1,2,7*
Virtual Institute of Microbial Stress and Survival,1 Physical Biosciences Division, Lawrence Berkeley National Laboratory, Berkeley, California,2 Earth Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California,3 Environmental Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee,4 Biochemistry and Molecular Microbiology & Immunology Departments, University of Missouri, Columbia, Missouri,5 Diversa Inc., San Diego, California,6 Departments of Chemical Engineering and Bioengineering, University of California, Berkeley, California7
Received 15 December 2005/ Accepted 14 March 2006
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The bacterial response to hyperionic stress includes a range of mechanisms, such as aquaporins for water intake and functions regulated by the stationary-phase sigma factor, RpoS (12, 53). A commonly encountered response is accumulation of neutral, polar, small molecules, such as glycine betaine (GB), proline, trehalose, or ectoine (20, 29, 33). These compatible solutes serve as osmoprotectants and are synthesized and/or imported into the cell. The genomic sequence of D. vulgaris provides insight into its potential responses to increased salt in the environment. While genes for trehalose synthesis (otsBA) are present in the closely related microbe Desulfovibrio sp. strain G20, homologues are not present in D. vulgaris. Furthermore, no genes for ectoine synthesis or transport can be identifed. D. vulgaris does contain genes for proline biosynthesis and an ABC transport system annotated for the uptake of GB, choline, and proline.
Compared to the responses to hyperionic stress, the responses specific to Na+ stress are less understood, since they are often observed in conjunction with responses to osmotic or alkaline stress (15). RpoS-regulated genes have also been found to be responsive to Na+ stress (12, 19). Also often implicated are Na+/H+ antiporters, such as those encoded by the nha genes in Escherichia coli (15, 56), orthologs of which are present in D. vulgaris.
A comprehensive understanding of various stress response mechanisms requires the use of multiple techniques to examine changes in different types of biomolecules. Integrated analyses involving a combination of microarray data and proteomics are increasingly being used for studies at a cellular level. Although mRNA and protein levels are linked, additional regulation during and after translation requires that expression at both stages be monitored for a more complete view of the physiological landscape (3, 24, 40). Although not addressed here, critical information concerning cellular proteomic responses also lies in posttransalational modifications. In this study, changes in the transcript and protein levels were examined. The data sets obtained were complemented by the results of a preliminary metabolite analysis. In addition, phospholipid fatty acid (PLFA) composition and lipid content, which are measures of changes in cell membrane properties, were also studied.
| MATERIALS AND METHODS |
|---|
|
|
|---|
MIC. MIC was defined as the stressor concentration that doubled the generation time and/or decreased the overall yield by 50%. Growth curve experiments were conducted in 96-well plates using an OmniLog instrument (Biolog Inc., Hayward, CA), which captured digital optical density images every 15 min for 150 h. Each well was inoculated with 10% mid-log-phase cells, six replicates of each stressor dilution were used, and the plates were placed in an anaerobic atmosphere and sealed in airtight Retain bags (Nasco) before they were placed in the OmniLog. OmniLog measurements were calibrated against D. vulgaris cell densities obtained by determining the OD600, using a Biolog plate reader (OD590), and determining direct cell counts by the acridine orange direct count method. All growth curves were comparable at 95% confidence intervals for the exponential growth phase. A kinetic plot of D. vulgaris growth was used to determine the generation time and cell yield. For NaCl and KCl stress, the MIC for D. vulgaris was determined to be 250 mM (in addition to the concentration present in LS4D medium) by testing concentrations ranging from 0 to 5,000 mM. For all experiments in this study, 250 mM NaCl or KCl was added to growth medium to establish stress conditions.
Biomass production. Highly controlled and reproducible conditions were used for simultaneous production of cell cultures for transcriptomics, proteomics, metabolite assay, PLFA, and synchrotron Fourier transform infrared (sFTIR) spectromicroscopy studies. Control and experimental production cultures were prepared in triplicate. Cultures were exposed to stress at the mid-log phase to minimize in-culture variability. NaCl or KCl (250 mM) was added to experimental cultures, and an equivalent volume of sterile distilled water was added to control cultures. The time of stressor addition was defined as zero time, and samples were taken at 30, 60, 120, and 240 min after exposure. The longest time (<5 h) was less than one generation time, so all samples were collected prior to the stationary phase. Samples were chilled to 4°C in <15 s during collection by pulling samples from production cultures through 7 m of capillary tubing immersed in an ice bath using a peristaltic pump. Chilled samples were centrifuged, and pellets were washed with 4°C degassed sterile phosphate-buffered saline and centrifuged again at 6,000 x g (10 min, 4°C). The final pellet was flash frozen in liquid N2 and stored at 80°C. For the production of each biomass sample, the following purity and growth characteristics were recorded: temperature, pH, OD600, acridine orange direct count, sFTIR spectrum, total protein content, anaerobic colony morphology, absence of aerobic colonies, and PLFA biosignatures.
Osmoprotection assays. D. vulgaris was cultured anaerobically in LS4D medium. Stressors, osmolytes, and other metabolites were added as indicated below, and growth was monitored by determining the OD600 and the total cell protein content.
Microarray analysis. Oligonucleotide probe design and microarray construction have been described previously (34). After RNA extraction, purification, and labeling, a reverse transcription reaction was used to generate labeled cDNA probes. Labeled genomic DNA (Cy3) was used as a control and as the common reference to cohybridize with labeled RNA (Cy5) samples for each slide. Each comparison was done in triplicate. Finally, since there were duplicate arrays on a given slide and three biological replicates, there were a total of 18 possible spots for each gene. Hybridized microarray slides were scanned using the ScanArray Express microarray analysis system (Perkin-Elmer, Massachusetts). Spot signals, spot quality, and background fluorescence intensities were quantified with ImaGene, version 5.5 (Biodiscovery Inc., Los Angeles, CA). For raw microarray data see NCBI GEO accession number GSE4447.
Computational. (i) Gene models. Models of The Institute for Genomic Research (NCBI) were used.
(ii) Microarray data analysis.
Log expression levels, including global normalization, were first computed for each microarray. Log expression levels obtained from replicate arrays were averaged. Each gene was represented by two spots on each microarray, and spots flagged by the scanning software were excluded. The net signal of each spot was calculated by subtracting the background signal and adding a pseudosignal of 100 to obtain a positive value. For resulting net signals of <50, a value of 50 was used. For each spot, the level of expression was the ratio of the two channels (ratio of mRNA to genomic DNA). For each replicate, the levels were normalized so that the total expression levels for the spots that were present on all replicates were identical. Finally, mean expression levels and standard deviations of each spot were calculated, which required n > 1. To estimate differential gene expression in control and treatment conditions, normalized log ratios were used. The log ratio was log2(treatment) log2(control). This log ratio was normalized using LOWESS on the difference versus the sum of the log expression level (16). Sector-based artifacts were observed; therefore, the log ratio was further normalized by subtracting the median of all spots within each sector. Up to this point, data were processed using spots instead of genes to allow sector-based normalization. Finally, the spots for each gene were averaged to obtain a final normalized log ratio. To assess the significance of the normalized log ratio, a Z score was calculated by using the following equation:
![]() |
Metabolite assays using CE-MS. Three identical 50-ml cultures were pooled anaerobically and centrifuged at 10,000 x g (10 min, 4°C). The pellet was resuspended in 5 ml cold methanol (prechilled on dry ice) by vortexing and allowed to stand for 5 min on dry ice. A mixture of 5 ml chloroform (Acros) and 2 ml water was added to the lysate, and the preparation was mixed thoroughly to remove free phospholipids. The extraction mixture was centrifuged at 7,000 x g for 5 min to separate the aqueous and nonaqueous layers. The top layer was transferred to a 5,000-molecular-weight-cutoff filter (Vivaspin6) and centrifuged at 7,000 x g (10 min, 4°C). The filtrate was lyophilized and reconstituted in 4 ml water. For desalting, solid-phase extraction (1g C18 cartridge; (Varian) was used. The C18 cartridge was conditioned using 4 ml methanol, followed by 4 ml water. The reconstituted sample (4 ml) was introduced into the cartridge and washed with 4 ml of water to remove salts. The sample was eluted with 4 ml methanol to which enough water was added for the sample to remain frozen at 80°C, and it was lyophilized (Labconco). Dried samples were reconstituted to produce metabolite mixtures in 100 µl buffer containing 115 µM methionine sulfone as the internal standard. Data were obtained in the positive ion mode by capillary electrophoresis (CE) (Agilent) and were analyzed using electrospray mass spectrometry (MS) (Agilent MSD). The metabolite mixture was resolved on a fused silica column (1 m by 50 µm; Polymicro Technologies) with 1 M formic acid as the electrolyte and 5 mM ammonium acetate as the sheath liquid. Metabolites were identified by comparing the nominal mass and retention time with the nominal masses and retention times of standards. Quantities of selected metabolites were estimated by comparison with known quantities of internal standards. All reagents used were high-performance liquid chromatography grade (Honey Well Burdick and Jackson).
ICAT and tandem liquid chromatography-mass spectrometry (LC-MS). To obtain total cell protein, cell pellets from 300-ml cultures at OD600 of 0.3 to 0.4 were lysed in 2 ml of 50 mM Tris-HCl (pH 8.0) by sonication on ice. The lysates were clarified by centrifugation at 14,000 x g (30 min, 4°C), and protein levels were determined using the bicinchoninic acid assay (Pierce). Samples (300 µg) of total cell protein from control and stressed samples (120 min) were labeled by isotope-coded affinity tagging (ICAT) with light and heavy tags, respectively, according to the manufacturer's instructions (Applied Biosystems). Trypsinized samples were desalted by using a C18 macro spin column (The Nest Group, Inc., Massachusetts) according to the manufacturer's instructions. Each sample was dried using a Speedvac (Thermo Savant), dissolved in 40 µl of 0.1% (vol/vol) 1 M formic acid, resolved by two-dimensional nano liquid chromatography (LC) using an inline strong cation-exchange cartridge into salt cuts (5, 20, 30, 40, 50, and 1,000 mM KCl), and resolved further by reverse-phase separation (15-cm 300-Å-diameter C18 RP column; Dionex). A 0 to 30% gradient in the organic phase (80% acetonitrile, 0.1% 1 M formic acid) was used. The m/z values of the resolved peptides were determined using electrospray ionization and a time of flight mass spectrometer (QSTAR hybrid quadrupole time of flight mass spectrometer; Applied Biosystems). Time of flight MS and MS-MS results were analyzed using the ProICAT software (Applied Biosystems). The acquired data were searched against a theoretical database created using the FASTA file containing the open reading frames of D. vulgaris Hildenborough, E. coli K-12, and Bacillus subtilis 168. A mass tolerance of 0.3 for MS and a mass tolerance of 0.3 for MS-MS were used to obtain a list of candidate proteins. The database of results generated by ProICAT computed the ID confidence for all peptides identified and the quantitation quality for ICAT peptide pairs for which a heavy/light ratio could be determined. Quantitation quality ranges from 0 to 100 and is based on how well the mass difference of the ICAT pair agrees with the expected value. Only proteins identified at >99% confidence were used. Of these proteins, only those with a stressed/control (heavy/light) quantitation quality of >75% were used. The ICAT strategy has been determined to be very accurate for small changes, detecting as little as a 30% change in protein levels (3, 61). Heavy/light ratios for peptides associated with the same protein were used to assess the internal error in a manner similar to the manner described previously (22). Internal error refers to the error between technical replicates of the same ICAT-labeled sample examined twice by LC-MS as described above. The error was determined to be 29%, and therefore, only changes greater than 29% were considered significant. The complete list is available at http://vimss.lbl.gov/SaltStress/.
Three-dimensional nano LC-MS proteomics. For the three-dimensional proteomics analysis, performed at Diversa Corp., we used fractionation of total proteins by three-dimensional LC, followed by MS-MS analysis to identify the proteins, as described previously (59). Replicate cultures from a control (time zero and 120 min) and a stressed sample (120 min) were used to obtain total protein. The relative abundance of proteins in each sample was estimated based on the hypothesis that the more abundant a peptide ion is in a mixture, the more likely the peptide ion is sampled during the course of an MS-MS experiment (35, 61). When this model is used, the total numbers of qualified spectral counts represent the relative abundance of each protein under a specific condition. To identify proteins for which there were significant changes under different conditions, the statistical "local-pooled-error" test (28) was used. Only changers with a P value of <0.05 were considered to be significantly changed. A total of 1,356 proteins were identified in all samples, and for 47 of these proteins there were reproducible changes between the control and the stressed sample (http://vimss.lbl.gov/SaltStress/).
PLFA assays. Pellets from 40-ml cultures (centrifuged at 1,015 x g for 15 min) were stored at 20°C. Total lipids were extracted from each pellet with a modified Bligh-Dyer solution (60). Phospholipids were separated from total lipids on a C18 silicic acid column (Unisil; Clarkson Chemical), methylated, and then analyzed with an Agilent 6890N gas chromatograph equipped with a flame ionization detector. Peaks were confirmed by MS (Agilent 5972A MSD), and double bond positions were confirmed by dimethyl disulfide derivatization. Peaks were quantified using an internal 19:0 fatty acid methyl ester standard (Sigma).
sFTIR spectromicroscopy. All sample handling was done anaerobically at 4°C. Cells were washed in ice-cold, oxygen-free buffer (0.1 M Tris-HCl [pH 7.5], 1.0 mM EDTA) and placed on a gold-coated microscope slide. The free-flowing buffer solution was removed, and cells were placed inside a sample holder, which maintained a relative humidity of 100% and an anaerobic atmosphere. The ZnSe window of the sample holder allowed acquisition of infrared spectra. Infrared measurements were obtained at the Advanced Light Source Facility, Lawrence Berkeley National Laboratory, using an Fourier transform infrared interferometer bench (Nicolet Magna 760) equipped with an infrared microscope (Nic-Plan IR microscope), as described previously (27). Signals in the spectral region from 4,000 to 250 cm1 were interpreted using previously described principles (9). Characteristic CO2 peaks and water vapor fingerprints were removed from the spectra, and the resulting spectra were normalized against the amide I absorption peak (1,648 cm1). Each spectrum was the average of 20 raw spectra, and each raw spectrum represented at least 20 cells. Differences in spectra, which revealed differences in cellular changes in control and stressed biomasses, were analyzed.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
Accumulation of compatible solutes. The D. vulgaris genome encodes a putative GB-choline-proline ABC transport system. The three-gene operon contains proV, encoding an ATP-binding protein; opuBB, encoding a permease protein; and proW, encoding the periplasmic binding protein that serves as the substrate receptor in the transporter complex. Microarray analysis showed that the mRNA levels of all three genes were highly up-regulated in response to both excess NaCl and KCl (Fig. 2B and Table 1). ICAT proteomics data for the same biomass supported the opuBB up-regulation (Fig. 3), and OpuBB was among the mostly highly up-regulated proteins.
|
|
|
In osmoprotection assays, very low concentrations of GB not only efficiently alleviated the growth inhibition due to excess salt but also prevented the late-log-phase cell elongation (Fig. 1). While this indicates that GB is an effective osmoprotectant, the up-regulated permease encoded by opuBB may not be specific to GB. Hence, additional substrates known to be effective osmoprotectants or precursors of GB were also tested. Our results showed that choline did not alleviate salt stress in D. vulgaris (Fig. 1). This may have been due to either the inability of D. vulgaris to import choline or the absence of specific dehydrogenases that convert choline to GB. GB has also been reported to be synthesized via sequential methylation of the glycine amino group by methyltransferases (44). Although several methyltransferases have been annotated in D. vulgaris, addition of glycine, methyl glycine (sarcosine), or dimethyl glycine to salt-stressed cells resulted in no significant improvement in growth (http://vimss.lbl.gov/SaltStress/). Organisms that use GB as an osmoprotectant, including E. coli and B. subtilis, often can synthesize it (11, 20). However, similar to Lactobacillus spp. (20), D. vulgaris may use GB exclusively via import.
The up-regulation of the GB-choline-proline ABC transport system also suggests that proline is a candidate for osmoprotection in D. vulgaris. Proline, a well-documented osmoprotectant (4, 11), may also accumulate due to up-regulation of biosynthesis, increased transport, or decreased proline oxidation to glutamate. Up-regulation of proline biosynthesis genes was not observed upon salt stress (http://vimss.lbl.gov/SaltStress/). Furthermore, CE-MS analysis indicated that cells did not accumulate proline during salt stress, while glutamate levels did increase (Fig. 3). In osmoprotection assays, proline did alleviate salt stress, but not as effectively as GB (Fig. 1). Therefore, despite the homology between the B. subtilis and D. vulgaris opuBB genes, proline appears not to be the preferred osmoprotectant for D. vulgaris.
In contrast, ectoine alleviated salt stress at concentrations comparable to the concentration of GB (Fig. 1). Interestingly, neither an ectoine biosynthesis pathway nor any ectoine-specific transport systems have been annotated in D. vulgaris, although it is evident that uptake and accumulation of small polar osmolytes are major salt stress responses in D. vulgaris.
Among other amino acids, the mRNA levels of tryptophan biosynthesis genes were markedly up-regulated during NaCl stress and, to a lesser extent, during KCl stress (http://vimss.lbl.gov/SaltStress/). We found no reports of tryptophan up-regulation in relation to salt or osmotic stress in other bacteria and have no explanation for this phenomenon. To investigate if the up-regulation at the transcript level indicated that there was a tryptophan requirement during salt stress, an osmoprotection assay was performed. However, addition of tryptophan did not alleviate the reduction in growth in the presence of 250 mM salt (http://vimss.lbl.gov/SaltStress/).
Efflux of excess salt ions. D. vulgaris contains more than 14 genes annotated for cation/multidrug efflux that may pump cations out of the cell using the same mechanism as Na+/H+ antiporters. Genes in a three-gene operon (DVU2815 to DVU2817), as well as genes in a two-gene multidrug efflux operon (DVU3326 and DVU3327), were highly up-regulated upon salt stress (Fig. 2). Transcripts of a putative acrAB system (DVU0058 to DVU0063) were highly up-regulated throughout the stress treatment in response to both NaCl and KCl (Fig. 2), a phenomenon also observed in E. coli in response to 0.5 M NaCl (39). The large number of up-regulated genes with cation/multidrug efflux domains suggests that efflux systems may play a major role in actively pumping excess salt ions out of the cell in D. vulgaris.
Several highly conserved Na+/H+ antiporters are also present in D. vulgaris, but for most of these antiporters the transcript and protein profiles indicated that there was no change in stressed cells (http://vimss.lbl.gov/SaltStress/). While this may indicate that under the conditions studied D. vulgaris does not recruit Na+/H+ antiporters to exclude excess salt ions, the transcripts of these open reading frames were abundant (Huang and Alm, unpublished observations), adequate levels of the antiporters may be present for salt efflux, and no further up-regulation may be required. One exception was the ubiquinone oxidoreductase subunit, encoded by echA (DVU0434), whose transcript levels were increased in the presence of both NaCl-stressed cells (at 120 min, Z = 2.12) and KCl-stressed cells (at 120 min, Z = 2.14). While no other gene in the operon showed a significant change (http://vimss.lbl.gov/SaltStress/), oxidoreductases such as EchA have been implicated in respiration-coupled Na+ efflux in several other bacteria (46, 49), and the observed DVU0434 up-regulation may be important in D. vulgaris salt stress.
The down-regulated genes included the genes involved in Na+ uptake, such as the genes encoding flagellar systems (Fig. 2), suggesting that D. vulgaris eventually becomes nonmotile under salt stress conditions. Similar down-regulation of the flagellar assembly genes has been observed in other bacteria, including Shewanella oneidensis (36) and B. subtilis (55). However, unlike S. oneidensis cells, the D. vulgaris cells were observed to be highly motile after salt stress (data not shown), and several key chemotaxis genes, such as cheY (DVU2073), were very highly and reproducibly up-regulated within 30 min of treatment with either NaCl or KCl (Fig. 2A). While cheY remained up-regulated throughout the later times monitored, other chemotaxis genes (cheA, cheW, and cheD) were most highly up-regulated only at earlier times (30 and 60 min) (Fig. 2A). The up-regulation of the transcript levels of chemotaxis genes may indicate an initial response in cells to move away from the stressful cations.
At the proteome level several periplasmic-binding proteins of ABC transport systems were significantly down-regulated (Table 2). While both proteomics techniques used in our study revealed this down-regulation, only one of the genes showed a similar change at the mRNA level (DVU1937; Z = 1.72). Whether the down-regulation of these periplasmic proteins is beneficial during salt stress requires further investigation.
|
Salt stress is known to affect the fluidity of the cell wall, and lipid regulation plays a role in cell wall fluidity (37). Eight PLFAs accounting for 80% of the total D. vulgaris PLFA composition (17) were used to study the effect of NaCl stress on the PLFA profile (Fig. 4). Since our study did not address the long-term adaptation to salt stress, no drastic changes were observed in either the total PLFA content or the net contribution of each individual PLFA (http://vimss.lbl.gov/SaltStress/). However, changes that were observed during salt stress were not observed in the control biomass over the same period of time (Fig. 4). The level of the branched, unsaturated PLFA i17:1w9c, which accounted for 14% of the total PLFA content and is characteristic of Desulfovibrio species (17), increased in stressed cells but not in control cells. In contrast, the levels of the unbranched, saturated PLFAs 16:0 (7%) and 18:0 (25%) decreased in salt-stressed cells but not in control cells. In general, the levels of branched PLFAs increased more in stressed cells than in unstressed cells.
|
The small increase in the level of lipids in the salt-stressed biomass was also observed in sFTIR experiments performed with whole cells (Fig. 4). No significant changes were found in the expression of genes encoding biosynthesis or modification functions of fatty acids or polysaccharides (http://vimss.lbl.gov/SaltStress/). However, ICAT proteomics did reveal increases in the levels of OmpH, which is putatively involved in lipopolysaccharide biosynthesis, and RfbA, which is involved cell wall biogenesis (Table 2).
Energy production. GB import, ion efflux, and helicases all require ATP hydrolysis, which may explain the significant up-regulation in F-type ATPases seen at both the transcript and protein levels (Fig. 2B and Table 2). Consistent with an increase in the energy requirement, the transcript levels for all genes (DVU0531 to DVU0536) in the operon encoding the high-molecular-weight cytochrome Hmc increased (Fig. 2B) in both NaCl and KCl stress conditions. Hmc is a potential transmembrane redox protein complex that has been reported to be involved in electron transport and energy production (25). Interestingly, expression of the Rrf12 operon (DVU0529 and DVU0530), which is located immediately downstream of the Hmc operon and is the predicted negative regulator for Hmc (32), also increased during salt stress (http://vimss.lbl.gov/SaltStress/).
Increases in transcript levels were observed for the members of the operon encoding part of the sulfate reduction pathway (DVU0846 to DVU0851) after treatment with NaCl (Table 1). Increases were also observed at the protein level for ApsA, the adenylyl sulfate reductase
subunit (DVU0847), and QmoB, the quinone-interacting membrane-bound oxidoreductase (DVU0849), although there was not significant up-regulation of QmoA (DVU0848) (Table 2). Sulfate is the terminal electron acceptor for anaerobic respiration of D. vulgaris. Taken together, these data indicate that there was an increase in sulfate reduction and suggest that stressed cells may require more energy than unstressed cells. In addition, the mRNA level for the sulfate permease (DVU0279) was up-regulated. This may cause a paradox for the salt-stressed cell, since this permease has been shown to function as an Na+ symporter (14).
An unexplained contradiction is the apparent down-regulation of several lactate permeases in both NaCl- and KCl-stressed cells (Fig. 2A and Table 1). Additional energy requirements would be expected to require greater lactate uptake. Preliminary data from enzymatic assays performed with cell culture supernatant did not reveal a change in the lactate consumed by salt-stressed cells (data not shown), but more experiments are required to obtain a better understanding of lactate metabolism. The pyruvate metabolic pathway also appeared to be down-regulated (Fig. 2A). However, there is sufficient redundancy in the pyruvate-to-acetyl coenzyme A pathway in the genome such that this metabolic flux could have remained unchanged during salt stress.
Comparison of stress responses to Na+ and K+. In many bacteria, K+ accumulation via the trk and kdp systems is a primary response to excess Na+ ions (7, 36, 41, 54). D. vulgaris has four genes (DVU0412, DVU0413, DVU1606, DVU2302) that are predicted to be involved in K+ uptake and transport and an operon homologous to the E. coli kdp K+ transport system (DVU3334 to DVU3339). However, microarray data did not indicate a significant change in any of these candidates except DVU3335, a histidine sensor kinase, whose transcript level decreased with KCl stress (Table 1). Most importantly, addition of KCl did not alleviate the detrimental effect of NaCl on growth. In fact, increasing K+ concentrations were found to be additive with Na+ for decreasing the growth rate (data not shown). Thus, at the NaCl concentration used here, uptake of K+ ions was not an important response mechanism in D. vulgaris.
The overall similarity of microarray data obtained after NaCl and KCl challenges to D. vulgaris was striking (Fig. 2 and Table 1). Very few differences were observed, and the differences were primarily differences for hypothetical proteins whose functions are unknown. Those with homologies to products of known genes included an ABC transport protein for branched amino acids (DVU2744), and several genes for methionine and leucine biosynthesis were up-regulated only in KCl stress conditions (http://vimss.lbl.gov/SaltStress/).
Regulatory systems. Several regulons have been identified in D. vulgaris (51), and the most striking differential expression involved increases in gene expression for members of the ferric uptake regulon (FUR) after salt stress (Fig. 2). FUR derepression was also observed in response to heat shock (13), oxygen exposure, and acid stress (He and Zhou, unpublished data) and might be a global stress response in D. vulgaris. In B. subtilis, high salinity causes iron deficiency, leading to derepression of FUR genes, particularly siderophore biosynthesis genes (26). Possibly because of its anaerobic lifestyle, D. vulgaris appears to lack a siderophore biosynthetic pathway. However, even though omission of Fe2+ from the LS4D growth medium severely affected growth, it did not affect the net sensitivity of D. vulgaris to salt or the ability of GB to alleviate growth deficiency (http://vimss.lbl.gov/SaltStress/). The ability of GB to protect against salt stress in the absence of Fe2+ suggests that FUR may not regulate the GB uptake operon.
In bacteria, two-component systems play an important role in initiating responses to external stimuli and stress. In the cyanobacterium Synechocystis, a histidine kinase involved in a salt stress response has been identified (47). Microarray data for D. vulgaris transcripts showed that there were changes in several histidine kinases during salt stress (http://vimss.lbl.gov/SaltStress/), but further experiments are needed to confirm their specific involvement. Other regulatory proteins, such as the stationary-phase RpoS
70 sigma factor, have also been reported to be involved in the salt stress response (12). Although some
-proteobacteria apparently have a well-recognized RpoS ortholog (43), D. vulgaris appears not to encode RpoS. Three putative
70 factors in D. vulgaris, rpoD (DVU1788), rpoH (DVU1584), and fliA (DVU3229), did not change significantly during salt stress, although
54, encoded by rpoN (DVU1628), appeared to be up-regulated at the transcript level (Table 1). Curiously, proteomics data also indicated that the level of RpoC, the RNA polymerase ß' subunit, increased in the salt-stressed cells (Tables 1 and 2).
Comparison of genomics data. The microarray and proteomics analyses performed in this study were targeted to study the stress response rather than adaptation. To ensure that stress was measured and also that the cells were not moribund, parallel studies were performed 2 h following salt addition. Although the maximum number of changers was observed at this time, most of the changes at the microarray level did not exceed fivefold. This is not surprising since the doubling time of D. vulgaris growing in LS4D medium at 30°C is 5 h. The data suggest that changes at the proteome level were even milder. Using the cysteine-tagging ICAT method, 167 proteins were identified, and 139 of these proteins had stressed/control ratios. Of these, 79 showed changes greater than 30% and were considered to be proteins that showed a significant change (Table 2) (http://vimss.lbl.gov/SaltStress/). Twenty-seven of the changing proteins showed trends similar to trends in mRNA levels (http://vimss.lbl.gov/SaltStress/). Differences between microarray and proteomics data have been observed in most studies which used both techniques to map cellular responses (3, 13, 21, 24). Such differences may arise because of the inherent lag between mRNA changes and changes in protein levels. Posttranslational modifications and regulation, such as turnover and protein secretion, could also contribute to the differences. However, further validation is required to confirm the cause or significance of any of these differences. Additionally, while a correlation between absolute protein quantity and mRNA value may be expected (35, 58), correlations between fold changes in microarray data and proteomics data appear to be more variable (3, 13, 21, 24).
Salt stress model for D. vulgaris Hildenborough. The different techniques used in this study allowed us to observe the initial effects and responses to salt stress and to formulate a more complete salt stress model for a single organism (Fig. 5). Uptake of available molecules, such as GB and ectoine, from the environment is a robust strategy to counter osmotic stress and appears to be the favored mechanism in D. vulgaris. Up-regulation at the transcript level also suggests that a variety of efflux mechanisms are involved. Cells may also utilize available protective mechanisms, such as Na+/H+ antiporters that are abundant but show no up-regulation upon stress. Other immediate responses include up-regulation of chemotaxis systems, apparent motility of the cells, and adjustments in cell wall composition. Increases in ATP production are supported at both the transcript and protein levels. Increases in the transcript levels of all members of the Hmc operon suggest that electron channeling may also increase during salt stress in D. vulgaris.
|
| ACKNOWLEDGMENTS |
|---|
We thank Keith Keller and Janet Jacobsen (Lawrence Berkeley National Laboratory) for the supplemental data website (http://vimss.lbl.gov/SaltStress/) and Sherry Seybold (Lawrence Berkeley National Laboratory) for graphics support.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles: