- Research article
- Open Access
A comprehensive genomic, transcriptomic and proteomic analysis of a hyperosmotic stress sensitive α-proteobacterium
© Kohler et al.; licensee BioMed Central. 2015
- Received: 4 December 2014
- Accepted: 11 March 2015
- Published: 26 March 2015
With the aim of remaining viable, bacteria must deal with changes in environmental conditions, including increases in external osmolarity. While studies concerning bacterial response to this stress condition have focused on soil, marine and enteric species, this report is about Caulobacter crescentus, a species inhabiting freshwater oligotrophic habitats.
A genomic analysis reported in this study shows that most of the classical genes known to be involved in intracellular solute accumulation under osmotic adaptation are missing in C. crescentus. Consistent with this observation, growth assays revealed a restricted capability of the bacterium to propagate under hyperosmotic stress, and addition of the compatible solute glycine betaine did not improve bacterial resistance. A combination of transcriptomic and proteomic analyses indicated quite similar changes triggered by the presence of either salt or sucrose, including down-regulation of many housekeeping processes and up-regulation of functions related to environmental adaptation. Furthermore, a GC-MS analysis revealed some metabolites at slightly increased levels in stressed cells, but none of them corresponding to well-established compatible solutes.
Despite a clear response to hyperosmotic stress, it seems that the restricted capability of C. crescentus to tolerate this unfavorable condition is probably a consequence of the inability to accumulate intracellular solutes. This finding is consistent with the ecology of the bacterium, which inhabits aquatic environments with low nutrient concentration.
- Osmotic stress adaptation
- Compatible solute
- Integrative omic analysis
- Caulobacter crescentus
An increase in the extracellular concentration of either ionic or non-ionic solutes is one of the several conditions to which microorganisms must cope with in order to remain viable. High external osmolarity leads to dehydration of the cytoplasm and reduction of the cell turgor [1,2]. The consequent increase in the cytoplasmic ion concentration causes disturbances on several cellular processes, such as central and secondary metabolic pathways, protein metabolism, and cell motility and development. As no active transport mechanism for water is present, microorganisms must employ one of two basic strategies to counteract the deleterious effects of the hyperosmotic stress. In members of the archaeal group Halobacteriaceae as well as in obligately halophilic bacterial anaerobes, all adapted to high-saline environments, potassium ions (K+) are accumulated at very high intracellular concentrations [3,4]. This is related to an extensive protein structure adaptation such as an increased proportion of acidic amino acid residues on the surface, attracting K+ and water molecules.
Most other microorganisms adapt to high osmolarity by accumulating a restricted range of highly soluble molecules with no net charge at physiological pH, referred to as compatible solutes owing to their compatibility with cellular processes at high internal concentrations [5,6]. In the Gram-negative bacterium Escherichia coli, where the hyperosmotic response has been extensively studied, the initial phase of osmoadaptation is characterized by a fast intracellular accumulation of K+ and the counter-ion glutamate [7,8]. While K+ is rapidly taken up from the environment within minutes upon the osmotic shock, glutamate is mainly synthesized and accumulates at a much slower rate. K+ uptake in E. coli is performed by the combined action of three specialized transport systems: the P-type ATPase Kdp, which is energized by ATP hydrolysis, and the Trk and Kup systems, that take up K+ in symport with a proton [7,8].
It is assumed that K+/glutamate act as a second messenger to trigger and coordinate the subsequent osmotic response, characterized by the uptake and/or synthesis of compatible solutes [9,10]. Glycine betaine is the preferred compatible solute for the majority of bacteria. Although de novo synthesis of glycine betaine is restricted to some phototrophic bacteria, many microorganisms, including E. coli, are able to take up choline from the surrounding environment and convert it in glycine betaine by the Bet system . In addition, glycine betaine released in the bathing medium by either the primary microbial producers or the decaying activity can be transported into the cells by the transmembrane protein ProP and the multi-component binding-protein-dependent transport system ProU. Both transport systems are also involved in the uptake of the compatible solute proline. Despite the widespread capability of synthesizing proline, Gram-negative bacteria accumulate this solute mainly by the transport systems ProP and ProU. Like proline, the compatible solute carnitine plays a role in osmoadaptation, accumulating into the cells mainly by transport from the external environment. Furthermore, endogenous osmoprotection in Gram-negative bacteria is achieved by synthesis of the nonreducing glucose disaccharide trehalose by using three distinct substrates: activated glucose (OstAB system), glucose oligosaccharides (TreYZ) and the reducing disaccharide maltose (TreS) [12-14].
Caulobacter crescentus is an oligotrophic free-living α-proteobacterium, normally found in habitats with very low concentrations of organic substances [15,16]. Genome analysis reveals a large number of genes that could enable an adaptation of the bacterium to dilute aquatic conditions . Besides continuous starvation, C. crescentus can also experience fluctuations in other environmental parameters, including increases in the external osmotic pressure. To date only a few genes were reported to be involved in C. crescentus response to hyperosmotic stress, including the paralogous extracytoplasmic function (ECF) sigma factors σT and σU [18-20]. Therefore, the capability of this bacterium to cope with hyperosmotic stress and the molecular systems involved in this response are still poorly characterized. Here, we list C. crescentus genes presumably involved in the synthesis or uptake of classical solutes related to osmoadaptation, analyze the sensitivity of cells to high extracellular concentrations of an ionic (sodium chloride) and a non-ionic (sucrose) solute, and describe changes in global gene expression, protein synthesis and metabolites to characterize the response of C. crescentus to hyperosmotic stress.
C. crescentus lacks most of the classical genes involved in hyperosmotic stress
With the aim of comparing the expected capability of C. crescentus to transport and synthesize classical solutes with that of other α-proteobacteria, we performed similar searches in the genomes of representative species of this bacterial group. Like C. crescentus, Bradyrhizobium japonicum, Phenylobacterium zucineum, Sphingomonas sp., Zymomonas mobilis, and Acetobacter pasteurianus possess genes related to potassium uptake, but none belonging to the TrK system (Figure 1A and Additional file 1: Table S1). However, genes presumably involved in the synthesis or transport of compatible solutes other than proline were identified in the genome of these α-proteobacterial species, suggesting a more pronounced capacity of these bacteria to undergo the secondary compatible solute based response to increased extracellular osmotic pressure.
The trK genes were identified only in Sinorhizobium meliloti, Rhodobacter sphaeroides, Roseobacter denitrificans, and Rhodospirillum rubrum (Figure 1A and Additional file 1: Table S1). Curiously, genes related to the uptake and synthesis of compatible solutes were found to be well represented in these α-proteobacteria, suggesting a higher capacity in coping with hyperosmotic stress. Considering the presumable ability to accumulate solute in response to hyperosmotic stress, Rickettsia rickettsii and Wolbachia sp. were the more compromised of the species analyzed, as these bacteria seem to possess only genes related to proline transport.
Taken together, the in silico analysis implies that members of the α-Proteobacteria group are very heterogeneous with respect to the expected capability to transport and synthesize classical solutes following osmotic upshock and that C. crescentus is one of the bacteria whose genome is deficient in genes involved in such biological process.
Sucrose exposure causes a more drastic impact on C. crescentus than salt stress
Sucrose treatment (from 50 to 110 mM) also caused an impact on C. crescentus growth in a dose-dependent manner (Figure 2A). However, in contrast to the effect of salt stress, exposure to high sucrose concentrations (110 mM or higher) resulted in decreased viability of the cells. Furthermore, cells at these high sucrose concentrations lost their characteristic morphology throughout the experiments, shifting to a spherical shape (Figure 2B) and becoming more susceptible to lysis after suspension in distilled water. Thus, sucrose stress produces a more pronounced effect upon C. crescentus.
Although glycine betaine is known to be a compatible solute that is taken up and/or synthesized in response to osmotic upshift , we found no protective effect of this compound upon growth and viability of C. crescentus cells exposed to either NaCl or sucrose at concentrations that decreased the rate of growth but did not lead to a complete growth arrest (data not shown). Combined, these results are in accordance to data from in silico analysis showing a very restricted or no capability of C. crescentus cells in accumulating classical solutes under hyperosmotic stress.
Global changes in gene expression and protein synthesis
To get a comprehensive picture of the adaptation process to an increase in osmolarity of the medium, we monitored gene expression changes in C. crescentus after an osmotic upshift using DNA microarrays. In addition, changes in protein levels of the cytosolic proteome were monitored by two-dimensional gel electrophoresis followed by protein identification by mass spectrometry. As treatment of cells with 110 mM sucrose resulted in a slight growth similar to that observed in the presence of 60 mM NaCl (Figure 2A), and a higher incorporation of L-[35S]-methionine into protein was measured under these conditions relative to the complete growth arrested cell cultures, 110 mM sucrose and 60 mM NaCl were chosen for the transcriptomic and proteomic analyses. To minimize secondary effects of the hyperosmotic stress, samples were taken just 30 min after the stress challenge.
By analyzing global transcriptome, 393 genes were identified as differentially expressed following osmotic upshift (Additional file 2: Table S2). In the proteome experiments, 125 protein spots displayed significant change in intensity after stress, and proteins corresponding to 72 spots could be identified (Additional file 3: Figure S1 and Additional file 2: Table S2). About two-thirds of the proteins identified correspond to genes also differentially expressed according to transcriptomic analyses. The remaining proteins were identified only by proteomic analyses, suggesting that both transcriptional and posttranscriptional controls are involved in the regulation of genes under hyperosmotic stress in C. crescentus. Together, our experimental procedures allowed us to identify a total of 425 open reading frames whose expression levels changed under at least one of the conditions tested (NaCl or sucrose stress), which are about 10% of all genes predicted in C. crescentus genome. Among those open reading frames identified in both experimental approaches as affected by osmotic stress, we could observe a perfect correlation between the transcriptome and the proteome data.
High salinity stress led to up-regulation of 143 genes and resulted in an increased synthesis of 27 proteins, whereas 209 genes and 33 proteins were found to be down-regulated (Additional file 2: Table S2). Sucrose stress increased the expression of 100 genes and 33 proteins and led to down-regulation of 113 genes and 33 proteins. A total of 75 genes and 25 proteins were up-regulated, and 97 genes and 29 proteins were down-regulated under both stress conditions. At a first glance, this seemed to imply that about half of the differentially regulated genes in the presence of NaCl were not regulated after sucrose stress and only few genes were differentially expressed exclusively in the presence of sucrose. However, most of the up- or down-regulated genes observed after high salinity exposure were regulated in the same way after sucrose stress, although many of them did not reach the selected threshold (log2 ratio of +1 or −1) (data not shown). Therefore, in general our data suggest that the two different osmotic stressors produce similar changes in gene expression. This is more evident in the protein levels, as just few proteins were differentially synthesized under only one stress condition (Additional file 3: Figure S1 and Additional file 4: Figure S2 and Additional file 2: Table S2).
Cellular processes affected by osmotic stress in C. crescentus
Interestingly, a close inspection in the relative abundance of transcript corresponding to the repertory of genes presumably involved in K+ uptake (kup and kdpABCDE), glutamate synthesis (glnA, gltAB and gdhA) and proline syntheis (proABC) showed no significant up-regulation. Instead, glnA and gltD were down-regulated after both sucrose and salt stresses (Figure 1B). Accordingly, three protein spots whose intensities decreased following hyperosmotic stress corresponded to GlnA, GltB and GltD. Therefore, besides the absence in C. crescentus of most of the genes related to classical solutes accumulation, increase in the corresponding transcript and protein under hyperosmotic stress was not ascribed to any of those that are present.
Data showed a general down-regulation of genes/proteins involved in translation, including almost all ribosomal proteins, several translation factors and some proteins required for tRNA synthesis, aminoacylation and ribosome biogenesis. The unique gene/protein acting on translation with increased expression is related to the hibernation factor YfiA, a protein responsible for inactivation of 70S ribosomes by direct binding to the 30S ribosomal subunit . A general decrease in expression was also observed for genes/proteins related to amino acid and folate biosynthesis. Several genes/proteins involved in chaperoning and secretion processes were also down-regulated in the presence of salt or sucrose. Thus, all these changes might suggest a strong reduction in protein synthesis following an increase in the external osmotic pressure.
Carbon and energy metabolism
We observed down-regulation of several genes/proteins participating in the central metabolism (glycolysis, pyruvate dehydrogenase complex, citric acid cycle, electron transport chain and ATP synthesis). On the other hand, some genes/proteins presumably involved in the Entner Doudoroff pathway were found up-regulated in osmotically stressed cells. Thus, the glucose metabolism seems to be redirected under hyperosmotic stress. Very striking was the observation of the increased expression of genes/proteins involved in the fatty acid metabolism and polyhydroxyalkanoate (PHA) synthesis. This observation points out to the possibility that synthesis of PHA under hyperosmotic stress relies on the fatty acid de novo synthesis pathway, as previously shown [22,23]. In addition, genes/proteins involved in the synthesis of polyhydroxybutyrate, a very common kind of PHA, were found up-regulated under the tested stress conditions. Therefore, more than one of the routes for the PHA synthesis might be activated when external osmotic pressure is increased.
Stress related metabolism and adaptations
Several of the up-regulated genes/proteins have been described previously as stress-responsive genes. In this category are a number of proteins required for degradation of denatured and misfolded proteins and those with a possible detoxification role, such as catalase, glutathione S-transferases, aldehyde dehydrogenases and cytochrome P450 family proteins. Up-regulation of genes/proteins required for detoxification of H2O2 and toxic organic compounds might suggest that C. crescentus suffers oxidative stress under hyperosmotic stress. Curiously, expression of enzymes that detoxify H2O2 with the expense of reduced coenzyme (glutathione peroxidase and alkyl hydroperoxide reductase) was reduced in the presence of salt or sucrose. Other up-regulated targets include structural porins and enzymes involved in synthesis and degradation of surface sugars. This observation might be indicative of substantial modification in bacterial envelope during the process of adaptation to hyperosmotic stress.
Transcription and regulatory function
Our analyses revealed a decrease in expression of some genes/proteins related to RNA metabolism, including proteins necessary for RNA unwinding, degradation and synthesis. Interestingly, genes/proteins related to the RNA polymerase core subunits and the principal RNA polymerase sigma factor σ73 (RpoD), which is essential for expression of housekeeping genes , were down-regulated. Among the regulatory proteins involved in transcription of selected genes under specific conditions, the extracytoplasmic function (ECF) sigma factor σT, σU and σM, and those involved in the control of expression of the two first sigma factors were found up-regulated [19,20]. In contrast, other regulatory proteins, including the ones involved in phosphate and nitrogen metabolism, were down-regulated. In sum, the changes in expression of these regulatory genes might be responsible for up- and down-regulation of genes following hyperosmotic stress.
Genes affected specifically or to a greater extent by one of the osmotic stressors
By comparing data of this analysis with that obtained when comparisons of transcriptome and proteome were carried out between osmotic shock (NaCl or sucrose exposure) and the non-stress condition, several of the genes/proteins were identified at higher levels under sucrose stress because of their stronger down-regulation following NaCl treatment with respect to sucrose exposure (Figure 4 and Additional file 4: Figure S2). Most of these genes code for ribosomal proteins, suggesting that protein synthesis is more negatively affected under salt stress than following sucrose exposure.
Only a few genes/proteins were affected exclusively or mainly by sucrose stress. A number of these genes are organized in two operons distantly located in the chromosome of C. crescentus. One of these operons is CC_1133 - CC_1135, whose genes code for proteins involved in transport across the bacterial envelope and enzymes acting upon carbohydrates. Upstream of this operon is another gene (CC_1136) involved in transport, which was also up-regulated only under sucrose stress (Figure 4). The genes CC_1631 and CC_1632 are part of the operon CC_1628 - CC_1632 and were expressed at higher levels under sucrose stress relative to the salt exposure (Figure 4). Although CC_1628 and CC_1630 were not differentially expressed by comparing transcripts of cells exposed to either sucrose or salt, both genes were up-regulated in the presence of sucrose with respect to the no stress condition. Therefore, all genes belonging to this transcriptional unit are expected to be specially up-regulated under sucrose stress. Genes of this operon also encode proteins with transport or catalytic function.
Many genes expressed at higher levels under salt stress relative to sucrose treatment displayed stronger up-regulation in the presence of NaCl than under sucrose exposure. Most of genes have no putative function assigned yet, and those functionally annotated are predicted to play distinct roles. Therefore, no functional category was enriched with genes expressed at higher levels under high salinity condition. Only one protein (a TonB-dependent receptor) was synthesized at low levels under sucrose stress relative to the salt treatment due to a more pronounced down-regulation after sucrose exposure. Therefore, although the similar general effect of NaCl and sucrose on genes and proteins expression, each stressor leads to small specific changes.
Analysis of cellular extracts indicates the absence of typical compatible solutes
Metabolites significantly changed under hyperosmotic stress in Caulobacter crescentus
NaCl vs control
Sucrose vs control
In sum, we could show that C. crescentus is not able to synthesize the compatible solutes sorbitol and trehalose following treatment with either salt or sucrose, ruling out their involvement as counteracting molecules during hyperosmotic stress. However, a validated protocol for intracellular metabolomic analysis of C. crescentus would be needed to clarify all other changes on the metabolome in a deeper detail.
This report describes the capability of C. crescentus in tolerating and responding to osmotic upshift triggered by either sugar or salt. Our data showed a very restricted tolerance of this bacterium to hyperosmotic stress, as even a relatively low increase in extracellular concentration of an ionic or non-ionic solute impairs the bacterial propagation in culture. Despite the similar consequence upon growth, salt and sucrose were found to be different with respect to the effect on survival, as cell viability was kept in the presence of high NaCl concentrations, but decreased under more expressive sucrose stress. This indicates that C. crescentus is more sensitive to sugar stress than to high salinity condition. This bacterium is well adapted to survive in environments characterized by low concentrations of carbon sources [15,16], where the scarce organic nutrients are scavenged in order to keep essential biological processes and so allow cells retain viability. The low probability of encountering abundant organic nutrients in such environments could have selected for organisms with restricted or no capability to tolerate high concentration of carbohydrates, as the absence of the molecular systems required for the response to these solutes confers the advantage of energy economy related to the maintenance of genetic information and expression of these systems. Accordingly, the distribution of C. crescentus in freshwater habitats might be related to the sensitivity to osmotic stress on the whole. Therefore, the ecology of C. crescentus suggests that this bacterium does not experience conditions of high external osmolarity in its natural habitats, or this condition is very sporadically encountered.
Data from genomic analysis is in accordance with the low tolerance of C. crescentus to high external osmotic pressure, as only some genes related to intracellular accumulation of solutes under this stress condition were identified. As these genes are involved in potassium uptake and synthesis of the amino acids glutamate and proline, their presence could just reflect the requirement of the ion and amino acids for the normal physiology of C. crescentus. The down-regulation (glutamate synthesis) or no change (potassium transport and proline synthesis) in expression of genes required for accumulation of osmoprotective solutes observed in transcriptomic and proteomic experiments could favor this hypothesis. The slight increase of glutamate and proline in C. crescentus after saline or sucrose stress might be linked with the increased protein degradation during these stresses rather than with activation of biosynthetic pathways. This assumption is supported by i) the up-regulation of various peptidases; ii) the general down-regulation of amino acid biosynthesis, especially those involved in the glutamate synthesis; iii) the increase in the levels of unrelated amino acids (phenylalanine, leucine and isoleucine). Thus, it could be speculated that increased protein degradation is important to supply cells with free amino acids for the biosynthesis of stress related proteins or just means a mechanism to avoid accumulation of misfolded and aggregated proteins. Besides the apparent inability to synthesize glycine betaine, the main compatible solute used by bacteria, C. crescentus cells were not beneficiated by the addition of this molecule in the bathing medium, probably because glycine betaine is not transported into the cells. In fact, genes corresponding to uptake systems for this compatible solute (BetT, ProP and ProU) are absent from the genome of C. crescentus.
The combined data of transcriptomic, proteomic and metabolomic analyses reported here provide reasons for the striking reduction of C. crescentus propagation in the presence of salt or sucrose. Down-regulation of transcripts/proteins related to amino acid biosynthesis, transcription, translation and protein folding suggests that protein synthesis rate is decreased in C. crescentus under hyperosmotic stress. In addition, energy conversion seems to be affected as judged by the decreased expression of transcripts/proteins corresponding to metabolic pathways involved in glucose degradation, the carbohydrate used by C. crescentus cells during experiments performed in this study. Glucose accumulation in intracellular extracts of stressed cells, as revealed by GC-MS analysis, corroborates this assumption. As protein synthesis and energy conversion are absolutely necessary for bacterial propagation and considering the extension to which they were affected, reduction of both processes might be the main cause of the decrease in growth rate after osmotic upshift. Down-regulation of transcripts/proteins involved in these processes might be caused by decreased expression of the principal sigma factor σ73. Conversely, we could rule out the involvement of cell cycle regulator genes in growth arrest during hyperosmotic stress, as no change in the corresponding transcript/proteins was observed.
Whereas an increase in the external osmotic pressure affects essential processes and so reduces bacterial propagation, C. crescentus cells are also able to respond to this unfavorable condition by increasing expression of transcripts/proteins related to environmental adaptation. More specifically, these transcripts/proteins could contribute to the modification of the cell envelope and elimination of reactive oxygen species, toxic organic compound and denatured and misfolded proteins that might accumulate in cells under the stress condition. In addition, production of polyhydroxyalkanoate in stressed cells is expected based on the transcriptomic and proteomic data. As the key role of polyhydroxyalkanoate for the hyperosmotic stress resistance in several bacterial species has been reported [27-29], our data could also suggest a possible involvement of accumulation of this class of organic compounds in the C. crescentus adaptation following osmotic upshift.
Induction of these functions under hyperosmotic stress appears to be triggered by the proteins supposed or already demonstrated to act in directing RNA polymerase to transcribe specific regulons. One of the up-regulated genes/proteins with this function, the ECF sigma factor σT, was recently reported to be absolutely essential for the osmotic stress response in C. crescentus, and 49 genes were identified as belonging to its regulon, including σU [19,20]. A total of 27 σT-dependent genes were found up-regulated under at least one of the hyperosmotic stress conditions tested here. The use of lower concentration of sucrose in our experiments (110 mM instead of 150 mM) could explain the absence of the other σT-dependent genes. Whereas the function and regulation of σT and σU are well described, nothing is known about the role of the ECF sigma factor σM, which was also significantly induced under salt stress condition.
In general, sucrose and NaCl triggered similar responses in C. crescentus. Just few genes/proteins were affected specifically or in a greater extent in the presence of one stressor relative to the other. Of particular interest were those transcripts/proteins up-regulated under sucrose stress when compared to cells in the presence of NaCl, as almost all are expected to be involved in transport or metabolism of carbohydrates. Therefore, these genes/proteins could function in carbohydrate utilization. Their induction when the extracellular concentration of carbohydrates increases could allow cells to quickly degrade or store these molecules, which are scarce in the natural environments of the bacterium. Alternatively, these genes and proteins could indirectly contribute to the adaptation to high concentrations of this sugar by degrading sucrose, thus causing a decrease in the osmotic pressure of the medium.
Interestingly, most of the genes/proteins affected by osmotic upshift were also differentially expressed upon carbon starvation in C. crescentus (see Additional file 2: Table S2 for comparison) , including those involved in essential processes such as transcription, translation, energy conversion and transport across membrane, besides a number of proteins of unknown function. Again, changes in expression of the regulatory proteins σD, σT, σU and σM appear to be responsible for the effect of carbon starvation upon the physiological functions. The role of σT and σU in this response is expected as these sigma factors mediate the general stress response in C. crescentus. Instead, the involvement of σM in hyperosmotic stress and carbon starvation is somewhat surprising, reinforcing the requirement of further characterization of these regulatory genes with the aim of better understanding their contribution for stress adaptation.
Among the free-living α-proteobacteria analyzed with respect to the ability to cope with an increase in the extracellular osmotic pressure, C. crescentus is the species with the more reduced set of genes related to accumulation of solute under hyperosmotic stress, and so with the lower expected tolerance to this stress condition. The sensitivity of C. crescentus to NaCl is quite similar to that of B. japonicum, with both bacteria growing at reduced rate in the presence of about 50 mM NaCl, but not at the concentration of 100 mM . According to our genomic analysis, these bacteria differ only with respect to the occurrence of trehalose related genes, which were found just in B. japonicum. The trehalose biosynthetic genes are induced under hyperosmotic stress, allowing accumulation of this compatible solute . However, since we could not detect trehalose by GC-MS analysis, we assume that trehalose does not have an impact as a compatible solute in C. crescentus.
Another interesting comparison that deserves attention is with Z. mobilis. The genomic analysis shown here in principle suggests a restricted ability to accumulate classical compatible solutes, similar to C. crescentus. However, Z. mobilis tolerates uncommonly high sucrose concentrations, due to its ability to synthesize and accumulate the unusual compatible solute sorbitol . From our metabolic data, we could exclude sorbitol as a compatible solute for C. crescentus. Several characterized α-proteobacterial species, such as S. meliloti, R. sphaeroides and R. rubrum are also very resistant to osmotic upshift, but for these bacteria this feature could be related to their intrinsic capability to accumulate classical compatible solutes by both uptake and synthesis. Finally, in intracellular species such as R. rickettsii and Wolbachia sp. that are encountered in more controlled environments the complete absence of genes related to accumulation of solutes is justified. In general, the capacity to transport and/or synthesize solutes is correlated with the versatility in occupying natural environments.
Even though the ability to accumulate solutes and the related resistance to hyperosmotic stress is very different among α-proteobacteria, an increase in the external osmotic pressure triggers changes in global gene expression of C. crescentus quite similar to those observed for other representatives of this phylogenetic group, including down-regulation of genes involved in protein synthesis and energy metabolism and up-regulation of genes acting upon cell envelope and detoxification functions [33-36]. All these changes are related to a stress condition, in which transcription and translation become reduced and more specific in order to favor physiological modifications required under the harmful condition. Therefore, the striking variation in the capability of α-proteobacteria in tolerating osmotic stress is mainly linked with their ability to accumulate solute.
Integrative data presented in this study contribute to expand our knowledge about the capacity of C. crescentus to tolerate an increase in the external osmotic pressure, as well as the strategies employed by this bacterium to respond to this stress condition. It was possible to correlate the low genomic representation of genes required for intracellular accumulation of solutes with the restricted capability of C. crescentus cells to propagate following an osmotic upshift. This correlation prompted us to suppose that the bacterium does not experience hyperosmotic stress in its natural environment or that the increase in external osmotic pressure is very weak, allowing cells to respond to this stress condition and keep viable. Results are therefore consistent with the fact that C. crescentus is very well adapted to aquatic habitats containing low amount of inorganic ions and organic compounds.
Bacterial strains and growth conditions
C. crescentus strain NA1000 was cultured at 30°C in M2 minimal salts medium plus glucose . Overnight culture of strain NA1000 was used to inoculate fresh M2 minimal salts medium and culture was grown to the exponential growth phase (OD600 = 0.7). For osmotic/salt stress experiments, cells were incubated in the presence of sodium chloride (NaCl) or sucrose at final concentrations ranging from 40/50 to 200 mM. Growth curves were constructed by incubating cultures for additional 4 hours. To measure cell viability, aliquots were taken, serially diluted and plated 4 and 8 hours after stress on M2 minimal medium for counting colony-forming units.
For isolation of RNA used in DNA microarray experiments, cultures of C. crescentus NA1000 were grown to exponential growth phase (OD600 = 0.7) and submitted for 30 minutes to stress (60 mM NaCl or 110 mM sucrose) or kept under control conditions (no stress). Cells (four aliquots of 2 ml from each treatment) were collected by centrifugation in a microcentrifuge for 1 min and suspended in 1 ml of Trizol Reagent (Invitrogen). After the extraction procedure according to manufacturer’s instructions, the integrity of the RNA was checked by agarose gel electrophoresis and samples were tested for the absence of DNA contamination by PCR.
Three distinct biological RNA samples isolated from control cells and from each stress condition were reverse transcribed and labeled by using the FairPlay III Microarray Labeling system (Agilent) as previously described [19,38]. Briefly, cDNA was synthesized from total RNA (20 μg) in the presence of amino allyl modified dUTP and random primer. After purification, the resulting amino-modified cDNA was chemically labeled by incorporation of the dyes Alexa Fluor 555 (Cy3) or Alexa Fluor 647 (Cy5). Labeled cDNAs were combined, mixed with Agilent hybridization buffer, and competitively hybridized to custom-designed Agilent microarrays according to the manufacturer’s instructions (Agilent). Data extraction and normalization was performed using Agilent Feature Extraction Software 188.8.131.52 (Agilent). The custom-designed arrays contain 9–11 probes covering a region around predicted translational start sites (−300 to +200 relative to the translational start site +1) of each gene. Only those probes downstream of the translational start site were considered for estimating the fold change of gene expression. Ratios obtained for probes corresponding to the same gene were averaged and genes showing a ratio log2 (stress - NaCl or sucrose/control) < −1 or log2 (stress - NaCl or sucrose /control) > 1 in all three biological replicates were considered as differentially expressed. Microarray datasets were deposited in GEO (GSE49654).
2-DE and image generation
Two distinct biological protein samples from each stress and control conditions were isolated. For protein analyses, cultures of C. crescentus NA1000 were grown to the exponential growth phase (OD600 = 0.7) and submitted for 30 minutes to stress (60 mM NaCl or 110 mM sucrose) or kept under control conditions. Afterwards, pulse-labeling was performed for 5 min by using 25 μCi of L-[35S]-methionine per mL. Then, 5 min later, L-[35S]-methionine incorporation was stopped by adding an excess of unlabeled methionine (1 mM) and chloramphenicol (100 mg/mL) and subsequently by transferring the sample to ice. After disruption of the harvested cells by sonication, the protein solution was separated from cell debris by centrifugation, and the protein amount and incorporated radioactivity were determined. Crude protein extracts (100 μg of protein) were loaded onto Pharmacia ready-made IPG strips (pH range 3–10) for the IEF of 2-DE as recommended. The separation in the second dimension was carried out as suggested by the manufacturer. After fixing with ethanol, acetic acid and water (40/10/50% vol/vol), the wet gels were dried on a chromatography paper backing by using a heated vacuum dryer. For autoradiography of the incorporated 35S methionine, dried gels were exposed to storage phosphor screens (Molecular Dynamics Storage Phosphor Screen, 20 by 25 cm) for a time span ensuring usage of the whole dynamic range by the strongest spot and corresponding to the amount of radioactivity separated on the gel. Screens were scanned using a Typhoon 9400 Variable Mode Imager (Amersham Biosciences) at 65 536 (16 bit) graylevels) and a resolution of 200 micron.
2-D gel analysis
To visualize changes in the protein pattern we used dual channel imaging as supported by Delta2D 4.3 (DECODON Greifswald). This technique involves overlaying the protein synthesis images of the control (pseudocolored green) and the protein synthesis image of the stress- NaCl or sucrose (pseudocolored red) onto each other . To ensure that results were not influenced by spot mismatches, distorted gels were adjusted by using the program DECODON Delta2D 4.3 . After pseudocolor overlaying by DECODON Delta2D, the histograms of the autoradiograms were normalized according their gray level integrated overall quantity. This procedure generates a dual channel image with an equal representation of detectable quantities of each sub-image. For spot detection on autoradiograms the images were fused using the union image fusion algorithm of Delta2D. After image fusion, spot positions and boundaries were detected on the fused image (proteome map). Their shapes and positions were transferred to each single autoradiogram and were subsequently quantified by using grey level integration inside spot boundaries and normalized by using 100% detected overall gel quantity as reference. An RSD of 30–40% is considered as sufficient in parallels of 2-D gel experiments. The reproducibility of our data was high, with more than 90% of all data of the identified proteins falling below an RSD of 35%.
For the extraction of differentially expressed spots, relative spot volumes were standardized (mean centering and division by standard deviation). The significance of changes in the protein synthesis of detected spots was determined by one-factorial ANOVA for each spot considering both the stress types (α = 0.1, P-values based on F-distribution) supported by Mev 4.1 . Low abundant spots with a spot volume less than 0.02% on all images were eliminated from these analyses. Finally, changes of the average values exceeding two-fold up- or down-regulation were considered as substantial.
For identification of proteins by mass spectrometry, preparative 2D PAGE was performed. Some 300 μg of protein extracts was loaded onto IPG strips (GE-Healthcare) in the pH range of 3–10. The resulting 2D gels were fixed in the presence of Sypro Ruby. The stained gels were scanned with the Typhoon 9400 Variable Mode Imager (Amersham Biosciences). For identification of proteins by MALDI-TOF-MS, Sypro Ruby stained protein spots were cut out from gels using a spot cutter (Proteome WorkTM) with a picker head of 2 mm and transferred into 96-well microtiter plates. Digestion with trypsin and subsequent spotting of peptide solutions onto the MALDI targets were performed automatically in the Ettan Spot Handling Workstation (GE-Healthcare) using a modified standard protocol . MALDI-TOF-MS analyses of spotted peptide solutions were carried out on a Proteome-Analyzer 4700 (Applied Biosystems, Foster City, CA, USA). The spectra were recorded in reflector mode in a mass range from 900 to 3700 Da. Automatic or manual calibration was performed as described previously . After calibration, peak lists were created using the ‘peak to mascot’ script of the 4700 ExplorerTM software. The resulting peak lists were analyzed by using the mascot search engine (Matrix Science, London, UK), GPMAW 4.1 (Lighthouse data). The annotation of C. crescentus CB15 was used for protein identification and denotation. Peptide mixtures that yielded at least twice a Mowes score of at least 50 and sequence coverage of at least 30% was regarded as positive identifications. Proteins that failed to exceed the 30% sequence coverage cut-off were subjected to MALDI-MS/MS .
Integrated transcriptome and proteome data representation
Treemaps of gene functional categories of C. crescentus have been constructed according to the TIGRFam classification . For this purpose the 2D plane was iteratively partitioned according to the TIGR main roles, sub-roles, C. crescentus predicted operons and proteins . The treemap cells on the deepest hierarchy level represent the C. crescentus proteins. Treemap cells were reordered to ensure that consecutive genes in the genome are symbolized by consecutive cells. Circular nodes represent the genes, the connectors in between link operon members. Gene expression was encoded by using a color ramp starting at blue (decreased expression) via grey (no change) and ending at orange (high expression). Messenger RNA data were used for coloring the gene circles and connectors, protein data for coloring the treemap cells. Dark grey cells have no assigned data.
Cell sampling for GC-MS analysis
C. crescentus were cultivated in M2 to an OD600 of 0.5 and stressed with either 60 mM NaCl or 110 mM sucrose. After 3 hours, equal amount of control and stressed cells were harvested by centrifugation at 12000 × g for 3 minutes at 2°C. Supernatants were discarded and the cell pellets were overlayed with 75% (w/v) ethanol −20°C. The tubes were imminently frozen in liquid nitrogen until the final sample preparation.
Sample preparation for GC-MS analysis
All samples were thawed on ice and the cell pellets were disrupted using a Ribolyser (MagNa Lyser, Roche, Germany), at 5000 m × s for 30 seconds. Samples were transferred on 50 ml falcon tubes and the ribolyser tubes were washed with a total of 5 ml 75% (w/v) ethanol −20°C until all cell rests are washed out and were transferred to the falcon tube. Next, the internal standards were added (20 nmol ribitol, 20 nmol norvaline and 2.5 nmol camphorsulfonic acid; Sigma-Aldrich) and the extraction solution was mixed and incubated for 10 min on ice. Afterwards, the Falcon tubes were filled with destilated water to the final volume of 35 ml and stored at – 70°C. Finally, the samples were lyophilized for GC-MS analysis.
Lyophilized samples were derivatized using a two-step method with MeOX (Sigma-Aldrich) and MSTFA (Chromatographie Service GmbH) as described . For identification and quantification of metabolites a GC-MS method was applied as described previously . Qualitative and quantitative analysis were performed using the ChromaTOF software (LECO Corporation). Identification of peaks was carried out by comparison of mass-spectra and retention time with an in-house database and the Fiehn GC/MS metabolomics RTL library (Agilent Technologies).
In silico search for genes related to accumulation of typical osmoprotective solutes
The Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used in a search for the annotated open reading frame displaying the highest identity to each gene known in E. coli to be involved in potassium transport (kup, trkAHG and kdpABCDEF), glutamate synthesis (glnA, gltBD and gdhA), glycine betaine synthesis (betABCI), glycine betaine transport (betS), glycine betaine, proline and carnitine transport (proP and proVWZ) proline synthesis (proABC) and trehalose synthesis (ostAB). This analysis resulted in an initial list of genes in representative α-proteobacterial species (Sinorhizobium meliloti 1021, Bradyrhizobium japonicum USDA110, Caulobacter crescentus CB15, Phenylobacterium zucineum, Rhodobacter sphaeroides 2.4.1, Roseobacter denitrificans, Sphingomonas sp. MM-1, Zymomonas mobilis subsp. mobilis ZM4, Rhodospirillum rubrum ATCC11170, Acetobacter pasteurianus IFO3283-01 Rickettsia rickettsii Sheila Smith and Wolbachia sp. wMel). Each gene from this list was used in a similar search for the annotated open reading frame displaying the highest identity in E. coli and the other α-proteobacterial species. In both searches, the threshold for Smith-Waterman score was 100 and the minimal overlap in local sequence alignment considered was 75% of the hit’s length. Genes identified in this best hit reciprocal search were assumed to be present in the α-proteobacteria analyzed. All genes uncovered by this in silico analysis are shown in Additional file 1: Table S1. Protein domains were ascribed according to the Pfam protein families database .
This work was supported by a grant to S.L.G. from Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP). R.F.L. and C.K. were postdoctoral fellows from FAPESP, and S.L.G. was partially supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq-Brazil). We thank Michael T. Laub for assistance with microarray analysis, Philipp Gierok and Michael Lalk for GC-MS analysis, and Anne Kohler, Luci D. Navarro and Sandra M. Fernandes for expert technical assistance.
- Yancey P. Compatible and counteracting solutes. In: Strange K, editor. Cellular and molecular physiology of cell volume regulation. Boca Raton: CRC Press; 1994. p. 81–109.Google Scholar
- Record Jr MT, Courtenay ES, Cayley DS, Guttman HJ. Responses of E. coli to osmotic stress: large changes in amounts of cytoplasmic solutes and water. Trends Biochem Sci. 1998;23(4):143–8.View ArticlePubMedGoogle Scholar
- Galinski E, Trüper HG. Microbial behaviour in salt-stressed ecosystems. FEMS Microbiol Rev. 1994;15:95–108.View ArticleGoogle Scholar
- Ventosa A, Nieto JJ, Oren A. Biology of moderately halophilic aerobic bacteria. Microbiol Mol Biol Rev. 1998;62(2):504–44.PubMed CentralPubMedGoogle Scholar
- Kempf B, Bremer E. Uptake and synthesis of compatible solutes as microbial stress responses to high-osmolality environments. Arch Microbiol. 1998;170(5):319–30.View ArticlePubMedGoogle Scholar
- Empadinhas N, da Costa MS. Osmoadaptation mechanisms in prokaryotes: distribution of compatible solutes. Int Microbiol. 2008;11(3):151–61.PubMedGoogle Scholar
- Csonka LNEW. Osmoregulation. In: Neidhardt FC et al., editors. Escherichia coli and Salmonella. Washington, DC: Cellular and molecular biology. ASM Press; 1996. p. 1210–23.Google Scholar
- Stumpe SSA, Schleyer M, Bakker EP. K+ circulation across the prokaryotic cell membrane: K+-uptake systems. In: Konings WN, Kaback HR, Lolkema JS, editors. Transport processes in eukaryotic and prokaryotic organisms. Amsterdam: Elsevier; 1996. p. 473–99.View ArticleGoogle Scholar
- Higgins CFCJ, Stirling DA, Sutherland L, Booth IR. Osmotic regulation of gene expression: ionic strength as an intracellular signal? Trends Biochem Sci. 1987;12:339–44.View ArticleGoogle Scholar
- Epstein W. Osmoregulation by potassium transport in Escherichia coli. FEMS Microbiol Rev. 1986;39:73–8.View ArticleGoogle Scholar
- Sleator RD, Hill C. Bacterial osmoadaptation: the role of osmolytes in bacterial stress and virulence. FEMS Microbiol Rev. 2002;26(1):49–71.View ArticlePubMedGoogle Scholar
- Giaever HM, Styrvold OB, Kaasen I, Strom AR. Biochemical and genetic characterization of osmoregulatory trehalose synthesis in Escherichia coli. J Bacteriol. 1988;170(6):2841–9.PubMed CentralPubMedGoogle Scholar
- Strom AR, Kaasen I. Trehalose metabolism in Escherichia coli: stress protection and stress regulation of gene expression. Mol Microbiol. 1993;8(2):205–10.View ArticlePubMedGoogle Scholar
- Empadinhas N, da Costa MS. Diversity and biosynthesis of compatible solutes in hyper/thermophiles. Int Microbiol. 2006;9(3):199–206.PubMedGoogle Scholar
- Lapteva NA. Ecological features of distribution of bacteria of the genus Caulobacter in freshwater bodies. Mikrobiologiya. 1987;56:537–43.Google Scholar
- Poindexter JS. The caulobacters: ubiquitous unusual bacteria. Microbiol Rev. 1981;45(1):123–79.PubMed CentralPubMedGoogle Scholar
- Nierman WC, Feldblyum TV, Laub MT, Paulsen IT, Nelson KE, Eisen JA, et al. Complete genome sequence of Caulobacter crescentus. Proc Natl Acad Sci U S A. 2001;98(7):4136–41.View ArticlePubMed CentralPubMedGoogle Scholar
- Zuleta LF, Italiani VC, Marques MV. Isolation and characterization of NaCl-sensitive mutants of Caulobacter crescentus. Appl Environ Microbiol. 2003;69(6):3029–35.View ArticlePubMed CentralPubMedGoogle Scholar
- Lourenco RF, Kohler C, Gomes SL. A two-component system, an anti-sigma factor and two paralogous ECF sigma factors are involved in the control of general stress response in Caulobacter crescentus. Mol Microbiol. 2011;80(6):1598–612.View ArticlePubMedGoogle Scholar
- Alvarez-Martinez CE, Lourenco RF, Baldini RL, Laub MT, Gomes SL. The ECF sigma factor sigma(T) is involved in osmotic and oxidative stress responses in Caulobacter crescentus. Mol Microbiol. 2007;66(5):1240–55.View ArticlePubMedGoogle Scholar
- Polikanov YS, Blaha GM, Steitz TA. How hibernation factors RMF, HPF, and YfiA turn off protein synthesis. Science. 2012;336(6083):915–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Rehm BH, Kruger N, Steinbuchel A. A new metabolic link between fatty acid de novo synthesis and polyhydroxyalkanoic acid synthesis. The PHAG gene from Pseudomonas putida KT2440 encodes a 3-hydroxyacyl-acyl carrier protein-coenzyme a transferase. J Biol Chem. 1998;273(37):24044–51.View ArticlePubMedGoogle Scholar
- Fiedler S, Steinbuchel A, Rehm BH. PhaG-mediated synthesis of Poly(3-hydroxyalkanoates) consisting of medium-chain-length constituents from nonrelated carbon sources in recombinant Pseudomonas fragi. Appl Environ Microbiol. 2000;66(5):2117–24.View ArticlePubMed CentralPubMedGoogle Scholar
- Fujita M, Tanaka K, Takahashi H, Amemura A. Transcription of the principal sigma-factor genes, rpoD and rpoS, in Pseudomonas aeruginosa is controlled according to the growth phase. Mol Microbiol. 1994;13(6):1071–7.View ArticlePubMedGoogle Scholar
- Sootsuwan K, Thanonkeo P, Keeratirakha N, Thanonkeo S, Jaisil P, Yamada M. Sorbitol required for cell growth and ethanol production by Zymomonas mobilis under heat, ethanol, and osmotic stresses. Biotechnol Biofuels. 2013;6(1):180.View ArticlePubMed CentralPubMedGoogle Scholar
- Loos H, Kramer R, Sahm H, Sprenger GA. Sorbitol promotes growth of Zymomonas mobilis in environments with high concentrations of sugar: evidence for a physiological function of glucose-fructose oxidoreductase in osmoprotection. J Bacteriol. 1994;176(24):7688–93.PubMed CentralPubMedGoogle Scholar
- Kadouri D, Jurkevitch E, Okon Y. Poly beta-hydroxybutyrate depolymerase (PhaZ) in Azospirillum brasilense and characterization of a phaZ mutant. Arch Microbiol. 2003;180(5):309–18.View ArticlePubMedGoogle Scholar
- Ayub ND, Pettinari MJ, Ruiz JA, Lopez NI. A polyhydroxybutyrate-producing Pseudomonas sp. isolated from Antarctic environments with high stress resistance. Curr Microbiol. 2004;49(3):170–4.View ArticlePubMedGoogle Scholar
- Langenbach S, Rehm BH, Steinbuchel A. Functional expression of the PHA synthase gene phaC1 from Pseudomonas aeruginosa in Escherichia coli results in poly(3-hydroxyalkanoate) synthesis. FEMS Microbiol Lett. 1997;150(2):303–9.View ArticlePubMedGoogle Scholar
- Britos L, Abeliuk E, Taverner T, Lipton M, McAdams H, Shapiro L. Regulatory response to carbon starvation in Caulobacter crescentus. PLoS One. 2011;6(4):e18179.View ArticlePubMed CentralPubMedGoogle Scholar
- Boncompagni E, Osteras M, Poggi MC, le Rudulier D. Occurrence of choline and glycine betaine uptake and metabolism in the family rhizobiaceae and their roles in osmoprotection. Appl Environ Microbiol. 1999;65(5):2072–7.PubMed CentralPubMedGoogle Scholar
- Sugawara M, Cytryn EJ, Sadowsky MJ. Functional role of Bradyrhizobium japonicum trehalose biosynthesis and metabolism genes during physiological stress and nodulation. Appl Environ Microbiol. 2010;76(4):1071–81.View ArticlePubMed CentralPubMedGoogle Scholar
- Chang WS, Franck WL, Cytryn E, Jeong S, Joshi T, Emerich DW, et al. An oligonucleotide microarray resource for transcriptional profiling of Bradyrhizobium japonicum. Mol Plant Microbe Interact. 2007;20(10):1298–307.View ArticlePubMedGoogle Scholar
- Fida TT, Breugelmans P, Lavigne R, Coronado E, Johnson DR, van der Meer JR, et al. Exposure to solute stress affects genome-wide expression but not the polycyclic aromatic hydrocarbon-degrading activity of Sphingomonas sp. strain LH128 in biofilms. Appl Environ Microbiol. 2012;78(23):8311–20.View ArticlePubMed CentralPubMedGoogle Scholar
- Tsuzuki M, Moskvin OV, Kuribayashi M, Sato K, Retamal S, Abo M, et al. Salt stress-induced changes in the transcriptome, compatible solutes, and membrane lipids in the facultatively phototrophic bacterium Rhodobacter sphaeroides. Appl Environ Microbiol. 2011;77(21):7551–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Dominguez-Ferreras A, Perez-Arnedo R, Becker A, Olivares J, Soto MJ, Sanjuan J. Transcriptome profiling reveals the importance of plasmid pSymB for osmoadaptation of Sinorhizobium meliloti. J Bacteriol. 2006;188(21):7617–25.View ArticlePubMed CentralPubMedGoogle Scholar
- Ely B. Genetics of Caulobacter crescentus. Methods Enzymol. 1991;204:372–84.View ArticlePubMedGoogle Scholar
- Kohler C, Lourenco RF, Avelar GM, Gomes SL. Extracytoplasmic function (ECF) sigma factor sigmaF is involved in Caulobacter crescentus response to heavy metal stress. BMC Microbiol. 2012;12:210.View ArticlePubMed CentralPubMedGoogle Scholar
- Bernhardt J, Buttner K, Scharf C, Hecker M. Dual channel imaging of two-dimensional electropherograms in Bacillus subtilis. Electrophoresis. 1999;20(11):2225–40.View ArticlePubMedGoogle Scholar
- Berth M, Moser FM, Kolbe M, Bernhardt J. The state of the art in the analysis of two-dimensional gel electrophoresis images. Appl Microbiol Biotechnol. 2007;76(6):1223–43.View ArticlePubMed CentralPubMedGoogle Scholar
- Saeed AI, Bhagabati NK, Braisted JC, Liang W, Sharov V, Howe EA, et al. TM4 microarray software suite. Methods Enzymol. 2006;411:134–93.View ArticlePubMedGoogle Scholar
- Eymann C, Dreisbach A, Albrecht D, Bernhardt J, Becher D, Gentner S, et al. A comprehensive proteome map of growing Bacillus subtilis cells. Proteomics. 2004;4(10):2849–76.View ArticlePubMedGoogle Scholar
- Haft DH, Selengut JD, White O. The TIGRFAMs database of protein families. Nucleic Acids Res. 2003;31(1):371–3.View ArticlePubMed CentralPubMedGoogle Scholar
- Bernhardt J, Funke S, Hecker M, Siebourg J. Visualizing gene expression data via voronoi treemaps, Proceedings of the 2009 Sixth International Symposium on Voronoi Diagrams. IEEE Comput Soc. 2009:233–241.Google Scholar
- Strelkov S, von Elstermann M, Schomburg D. Comprehensive analysis of metabolites in Corynebacterium glutamicum by gas chromatography/mass spectrometry. Biol Chem. 2004;385(9):853–61.View ArticlePubMedGoogle Scholar
- Gierok P, Harms M, Richter E, Hildebrandt JP, Lalk M, Mostertz J, et al. Staphylococcus aureus alpha-toxin mediates general and cell type-specific changes in metabolite concentrations of immortalized human airway epithelial cells. PLoS One. 2014;9(4):e94818.View ArticlePubMed CentralPubMedGoogle Scholar
- Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42(Database issue):D222–30.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.