Metabolomic Characterization of the Salt Stress Response in Streptomyces coelicolor Kol, Stefan; Merlo, M. Elena; Scheltema, Richard A.; de Vries, Marcel; Vonk, Roelf; Kikkert, Niels A.; Dijkhuizen, Lubbert; Breitling, Rainer; Takano, Eriko Published in: Applied and Environmental Microbiology DOI: 10.1128/AEM.01992-09 IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check the document version below. Document Version Publisher's PDF, also known as Version of record
Publication date: 2010 Link to publication in University of Groningen/UMCG research database
Citation for published version (APA): Kol, S., Merlo, M. E., Scheltema, R. A., de Vries, M., Vonk, R. J., Kikkert, N. A., ... Takano, E. (2010). Metabolomic Characterization of the Salt Stress Response in Streptomyces coelicolor. Applied and Environmental Microbiology, 76(8), 2574-2581. DOI: 10.1128/AEM.01992-09
Copyright Other than for strictly personal use, it is not permitted to download or to forward/distribute the text or part of it without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license (like Creative Commons). Take-down policy If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. Downloaded from the University of Groningen/UMCG research database (Pure): http://www.rug.nl/research/portal. For technical reasons the number of authors shown on this cover page is limited to 10 maximum.
Download date: 16-01-2017
Vol. 76, No. 8
Metabolomic Characterization of the Salt Stress Response in Streptomyces coelicolor£‡ Stefan Kol,1† M. Elena Merlo,1,2† Richard A. Scheltema,2† Marcel de Vries,3 Roel J. Vonk,3 Niels A. Kikkert,2 Lubbert Dijkhuizen,1 Rainer Breitling,2* and Eriko Takano1* Microbial Physiology, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Kerklaan 30, 9751 NN Haren, Netherlands1; Groningen Bioinformatics Centre, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Kerklaan 30, 9751 NN Haren, Netherlands2; and Centre for Medical Biomics, University Medical Centre Groningen, 9713 AV Groningen, Netherlands3 Received 19 August 2009/Accepted 18 February 2010
The humicolous actinomycete Streptomyces coelicolor routinely adapts to a wide variety of habitats and rapidly changing environments. Upon salt stress, the organism is also known to increase the levels of various compatible solutes. Here we report the results of the first high-resolution metabolomics time series analysis of various strains of S. coelicolor exposed to salt stress: the wild type, mutants with progressive knockouts of the ectoine biosynthesis pathway, and two stress regulator mutants (with disruptions of the sigB and osaB genes). Samples were taken from cultures at 0, 4, 8, and 24 h after salt stress treatment and analyzed by liquid chromatography-mass spectrometry with an LTQ Orbitrap XL mass spectrometer. The results suggest that a large fraction of amino acids is upregulated in response to the salt stress, as are proline/glycine-containing diand tripeptides. Additionally we found that 5*-methylthioadenosine, a known inhibitor of polyamine biosynthesis, is downregulated upon salt stress. Strikingly, no major differences between the wild-type cultures and the two stress regulator mutants were found, indicating a considerable robustness of the metabolomic response to salt stress, compared to the more volatile changes in transcript abundance reported earlier. phy (HILIC) (23). The accurate mass determination and high resolving power of the Orbitrap (22, 40) have made the analysis of complex metabolite mixtures feasible, making the Orbitrap an attractive option for untargeted metabolomics screens (5). The use of HILIC ensures that the resolving power of the mass spectrometer is further extended by fractionation of the sample, resulting in a less complex mixture being injected in the mass spectrometer at each moment in time. In addition to wild-type S. coelicolor, we studied strains with disruptions in various steps of the biosynthesis of ectoine and hydroxyectoine (Fig. 1), the two best-characterized osmoprotectants of Streptomyces (6, 7), as well as mutants with mutations in two stress regulators (OsaB and sB) under continuous salt stress and during a time course following salt shock. The two stress regulators were previously shown to be involved in the salt shock response (3, 31). This study provides the first attempt at a global metabolomic characterization of the salt adaptation process in bacteria.
Salt stress conditions are known to cause major changes in the primary and secondary metabolism of bacterial cells (26, 34). This is also the case for the Gram-positive soil bacteria of the genus Streptomyces, which are well known for their complex secondary metabolism and production of a wide range of industrially relevant metabolites, including various antibiotics (21). Upon salt stress, these species are known to increase the levels of various compatible solutes, including ectoine, alanine, glutamine, and proline (6, 11, 16, 28). The genome sequence of Streptomyces indicates a largely unexplored capacity for production of secondary metabolites (2), some of which could also be involved in the salt stress response. Here we conducted an untargeted metabolomics analysis of salt-stressed Streptomyces coelicolor, the genome-sequenced model species of the genus (2), to investigate the complexity of the metabolite changes associated with the salt stress response in more detail and to increase our understanding of the mechanisms of osmoadaptation by metabolomic rearrangement. For this purpose, we used a high-resolution LTQ Orbitrap mass spectrometer coupled with hydrophobic-interaction liquid chromatogra-
MATERIALS AND METHODS Strains and medium conditions. Escherichia coli strain ET12567/pUZ8002 was used for conjugation of cosmids into S. coelicolor (14). Salt shock or osmoadaptation experiments using the S. coelicolor M145 parent strain (27), LW35 (ectA::Tn [where Tn is transposon 5062]) (this study), LW36 (ectC::Tn) (this study), LW37 (ectD::Tn) (this study), the osaB::Tn strain (3), and the sigB::aac(3)IV strain (DsigB) (30) were performed with Difco nutrient broth (DNB), whereas solid growth experiments were performed with supplemented minimal medium (SMMS) with various amounts of NaCl or ectoines added (20). For the liquid cultures, a high inoculum of spores was used to reduce the formation of clumps. Ectoine mutant construction and characterization. Three cosmids containing Tn5062 transposon insertions in ectA (SCO1864), ectC (SCO1866), and ectD (SCO1867) were used to disrupt the ectoine biosynthesis cluster genes in S. coelicolor. Conjugation to S. coelicolor was performed as previously reported (3). Allelic replacements were screened by selecting for Aprr and Kans clones. To
* Corresponding author. Mailing address for E. Takano (Streptomyces physiology): Microbial Physiology, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Kerklaan 30, 9751 NN Haren, Netherlands. Phone: 31-50-3632143. Fax: 31-503632154. E-mail:
[email protected]. Mailing address for R. Breitling (metabolomics): Groningen Bioinformatics Centre, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Kerklaan 30, 9751 NN Haren, Netherlands. Phone: 31-503638088. Fax: 31-50-3637976. E-mail:
[email protected]. † Equally contributing authors. ‡ Supplemental material for this article may be found at http://aem .asm.org/. £ Published ahead of print on 26 February 2010. 2574
DOWNLOADED FROM aem.ASM.ORG - at University of Groningen July 13, 2010
APPLIED AND ENVIRONMENTAL MICROBIOLOGY, Apr. 2010, p. 2574–2581 0099-2240/10/$12.00 doi:10.1128/AEM.01992-09 Copyright © 2010, American Society for Microbiology. All Rights Reserved.
METABOLOMICS OF S. COELICOLOR
2575
FIG. 1. Biosynthesis pathway for ectoine. L-Aspartate-b-semialdehyde (ASA) is converted into L-2,4-diaminobutyrate (DABA) by the enzyme DABA aminotransferase, encoded by the ectB gene. DABA is converted into N-g-acetyl-L-2,4-diaminobutyrate (ADABA) by DABA acetyltransferase, encoded by the ectA gene, and L-ectoine (Ect) is formed from ADABA by ectoine synthetase, encoded by the ectC gene (8, 33, 35). Ectoine can then be hydroxylated to form 5-hydroxyectoine (EctOH) by the enzyme ectoine hydroxylase, encoded by the ectD gene (7, 17). CoA, coenzyme A.
confirm the genomic replacement of the ectA, ectC, and ectD genes with the Tn5062-containing copies, genomic DNA from each mutant was isolated and used as a template for PCR using one of the transposon primers EZR1 and EZL2 (3) and ect gene primers which anneal to the start and stop codons of ectA, ectC, and ectD (for primer sequences, see Table SB in the supplemental material). PCR conducted with wild-type strain M145 total DNA as a template gave rise to amplified products which corresponded in size with the predicted products of 633 (ectA), 1,250 (ectB), 414 (ectC), and 888 (ectD) bp by use of the ect gene primers (see Fig. SA.b, lanes 1 to 4, in the supplemental material). Due to the insertion of the large transposon, no PCR product was amplified when only the ect gene primers were used with genomic DNA from the mutants (see Fig. SA.b, lanes 7, 10, and 13). However, the transposon primers in combination with the ect gene primers yielded products corresponding in size with the predicted distance from the ect start and stop codons to the location of the transposon (see Fig. SA.b, lanes 5, 6, 8, 9, 11, and 12). These data confirmed the genomic disruption of the ectA, ectC, and ectD genes. This was further validated by a complementation analysis using a drop dilution assay. Cells were grown on SMMS with 1 M NaCl, with and without an additional 20 mM ectoine and hydroxyectoine, and monitored for growth (Fig. 2). Growth conditions and metabolite extraction. Flasks containing 50 ml DNB were inoculated with 1 3 108 S. coelicolor spores of the M145, ectA::Tn, ectC::Tn, ectD::Tn, osaB::Tn, or DsigB strain, and cultures were grown for 24 h, at which point 5 ml of 5 M NaCl or 5 ml of MilliQ water was added for salt shock samples. Cells exposed to continuous salt stress were inoculated and grown in 50 ml DNB supplemented with 0.5 M NaCl. Samples of 5 ml (optical density at 450 nm [OD450] of 1.0) or an equivalent thereof needed to obtain an identical cell mass were harvested at 0, 4, 8, and 24 h after addition of salt for the salt shock cultures or at 24 and 48 h for the continuously salt-stressed cultures. Cells were collected by centrifugation for 10 min at 4,500 rpm and 4°C, the supernatant was removed completely, and the cell pellet was resuspended in 500 ml methanol. After incubation on ice for 10 min with occasional mixing, cells were spun down at 4°C, 13,000 rpm, and the supernatant was collected and stored at 280°C. LC-MS. Data were acquired with a Luna 3-mm-particle-size HILIC 200A high-pressure liquid chromatograph (HPLC) (150 by 2.0 mm; Bester, Rotterdam, Netherlands) coupled to an LTQ Orbitrap XL mass spectrometer (Thermo Fisher Scientific, Bremen, Germany) in positive-ionization mode at a resolution of 30,000. The liquid chromatography-mass spectrometry (LC-MS) system was run in binary gradient mode with a flow rate of 400 ml/min; a volume of 20-ml methanol samples was injected. The solvents used were 90/10 acetonitrile with 5 mM ammonium acetate (solvent A) and 50/50 acetonitrile with 5 mM ammonium acetate (solvent B). The gradient used was set to hold 90% acetonitrile isocratic for 2.5 min, after which it was moved to 50% for 7.5 min and finally held isocratic for 2.5 min. The same conditions were maintained for the tandem MS (MS-MS) measurements. Analysis of the fragmentation patterns was done with XCalibur 2.0.5. Extracts were supplemented before injection with solvent A at a ratio of 1:10, improving the LC peak shape. Even though molecules were diluted 10-fold, the peak intensity levels of compounds of interest were not affected negatively, due to sharper peaks. Ectoine and hydroxyectoine standards (Fluka, Zwijndrecht, Netherlands) were used as controls. Data processing. Data processing was performed using a configurable software pipeline capable of extracting and automatically matching peaks from multiple measurements into peak groups (40). Retention time alignment of the extracted mass chromatograms was done using the COW-CoDA algorithm (9). Noise resulting from the greedy peak picking approach was reduced by application of a relative standard deviation (RSD) filter set to 35% for matched replicates, as
FIG. 2. Ectoine biosynthesis disruption mutants display a salt-sensitive phenotype, which is complemented by the extracellular addition of ectoines. S. coelicolor spore stocks of the M145, ectA::Tn, ectC::Tn, ectD::Tn, DsigB, and osaB::Tn strains were diluted so that a spot of 2.5 ml contained approximately 105, 104, 103, 102, 101, or 100 spores. Spots were plated on SMMS without any addition (A, left) or with the addition of 1 M NaCl (A, right), 1 M NaCl and 20 mM ectoine (B, left), 1 M NaCl and 20 mM hydroxyectoine (B, right), or 1 M NaCl, 20 mM ectoine, and 20 mM hydroxyectoine (C). Plates were incubated at 30°C for 3 days in the absence of salt and for 6 days in the presence of salt.
DOWNLOADED FROM aem.ASM.ORG - at University of Groningen July 13, 2010
VOL. 76, 2010
KOL ET AL.
APPL. ENVIRON. MICROBIOL.
FIG. 3. Levels of ectoine and related compounds detected under continuous salt stress reflect the ectoine mutants. The presence of DABA, ADABA, ectoine, and hydroxyectoine was determined in wild-type cells and in mutants lacking ectA, ectC, or ectD grown under continuous salt stress. DABA was detected only in the etcA::Tn strain, which is expected to be deficient in the conversion of DABA into ADABA and therefore accumulates the precursor. ADABA, another intermediate compound, accumulates in the ectC::Tn strain and to a lesser degree in the ectD::Tn strain, as expected. The etcD::Tn strain, which is deficient in the conversion of ectoine to hydroxyectoine, accumulated large amounts of ectoine, while hydroxyectoine accumulated in the M145 strain, as expected.
quantification is expected to be at least 20% accurate over multiple runs (41). The selected mass chromatograms were putatively identified by matching the masses progressively to those from metabolite-specific databases. First, ScoCyc (2), LIPID MAPS (12), and a contaminant database (25) were used. Unidentified peak groups were then matched to KEGG (24) and the remaining unidentified ones to METLIN (42) and the Human Metabolome Database (47). This iterative process was used in order to restrict the number of potential matches to the most likely (39). As a last step, all peak groups related to another peak group (e.g., isotopes, ion adducts, and fragments) were removed from the set based on a correlation analysis of the intensity patterns and peak shape (46). Hierarchical clustering was performed to identify metabolites that showed similar dynamic changes after salt shock with various mutants. To achieve a robust clustering, the distance score was calculated for data with sufficient variation (RSD of .0.2) based on the Pearson correlation of the discretized slopes of the z-score-normalized time courses, rather than on the absolute values. By this approach, instead of determining whether the original signal intensities correlate, we determined the correlation of the slopes, i.e., the pattern of changes over time, and in order to reduce the influence of noise on the slopes, they were discretized to three possible values (21, 0, or 1, respectively) depending on whether signal intensities were decreasing, stable, or increasing. The resulting hierarchical tree was consequently cut into 4 clusters. For each cluster, a quality score was calculated based on the average distance of each observation for each strain to the average profile of the cluster; this values indicates the coherence of expression profiles within the cluster (smaller values indicate better quality). For convenient visualization of the time courses for individual metabolites, the intensity values are displayed as fold changes relative to the initial time point (i.e., for each strain, the values are divided by the intensity value of the first time point for that strain).
RESULTS AND DISCUSSION Validation of metabolite quantification. The relative quantification accuracy of the mass spectrometry data was assessed by quantifying the S. coelicolor antibiotic undecylprodigiosins, which are visibly present in the extracts (the extracts are dark red when large amounts of antibiotics are produced). The traditional method of measuring the abundance of undecylprodigiosins by OD533 measurement (44) was used as a second quantification method. Figure SB in the supplemental material shows the total detected signal for undecylprodigiosin (mass of 391.26236 Da for butylcycloheptylprodigiosin, metacycloprodiginine, and methylcyclo-decyl-prodiginine; mass of 393.27801 Da for undecylprodiginine) compared to the OD533 values. The Pearson correlation is 0.95, indicating that the relative quantification achieved by the HILIC-Orbitrap combination is surprisingly precise, especially considering that these compounds do not bind well to the HILIC column and are flushed out with a large number of other compounds (mostly lipids), making exact quantification challenging due to ionization suppression effects (23). For the large number of more polar compounds in our metabolome screen, which are better separated by the HILIC column, we can therefore expect similarly reliable quantitation.
Metabolomics of continuous salt exposure. To validate our methodology, we first focused on a targeted metabolomics screen of the wild type and the ectoine biosynthesis mutants of S. coelicolor grown under continuous salt stress. Under these conditions, we expect very characteristic abundance profiles for ectoine, hydroxyectoine, and their precursors; these critical osmoprotectants should be absent in the mutants that show a clear growth defect under salt stress conditions (Fig. 2). This was indeed the case in our measurements, as the intensity levels that we found closely matched the expected phenotype of the knockouts (Fig. 3). (i) The early precursor metabolite L-2,4-diaminobutyrate (DABA) was found to accumulate in the ectA::Tn strain, where its conversion to N-g-acetyl-L-2,4diaminobutyrate (ADABA) is deficient, and was detected nowhere else. (ii) The later precursor metabolite ADABA was found to accumulate strongly in the ectC::Tn strain and, to a lesser degree, in the ectD::Tn strain. The presence of this metabolite in the ectD::Tn strain can be explained by the reversibility of its conversion to ectoine. (iii) Ectoine was found to accumulate strongly in the ectD::Tn strain, where its conversion to hydroxyectoine is deficient. Low intensity levels of ectoine were also detected for the M145 strain, suggesting that ectoine acts as an intermediate in S. coelicolor and is converted almost quantitatively to hydroxyectoine. Very low levels of the same mass and similar retention times for ectoine were also detected for the ectA::Tn and ectC::Tn strains; this could indicate the existence of an alternative, low-yield biosynthesis pathway for ectoine or the presence of an unrelated isomer of low abundance. In both cases, the detected levels were so low that we do not expect this to be of biological relevance for osmoprotection. (iv) Finally, hydroxyectoine was found to accumulate strongly in the M145 strain, indicating that it acts as the prime osmoprotectant. Retention time and tandem MS fragmentation patterns in comparison with the standards confirmed the identities of the detected compounds as ectoine and hydroxyectoine in the ectD::Tn and M145 strains, respectively (see Fig. SC and SD in the supplemental material). Metabolomics of the salt shock response. (i) Global metabolome screen. To obtain a global metabolomic characterization of the salt shock response, we measured time series for metabolite profiles of salt-shocked cells, using the same S. coelicolor parent strain (M145) and mutant strains (ectA::Tn, ectC::Tn, ectD::Tn, osaB::Tn, and DsigB strains). A total of 1,247 distinct peak groups were quantified and were combined into 363 peak groups of related peaks (each corresponding to a potential metabolite). Of the peak groups, 229 could be
DOWNLOADED FROM aem.ASM.ORG - at University of Groningen July 13, 2010
2576
METABOLOMICS OF S. COELICOLOR
2577
FIG. 4. Overview of the putatively identified metabolites detected in the global metabolome screen. A total of 229 metabolites could be assigned putative identities based on database matching. By far the largest class of detected metabolites were lipids (107 metabolites) and in particular a large number of short- and medium-chain fatty acids. The second-largest class contained amino acids (AA) and their derivatives (66 metabolites), including a large number of di- and tripeptides, many of which contained proline and glycine residues.
assigned putative identities based on exact mass matching (for an overview of the identifications, see Fig. 4 and Table SA in the supplemental material). This compares well to the predicted amount of metabolites based on genome annotations reported by Borodina et al. (4). (ii) Dynamic responses to salt exposure. The dynamic response to salt exposure was visualized by unsupervised clustering of the metabolite levels for the time series of the saltstressed cultures. Two coherent clusters of metabolites were revealed, one with consistently increasing abundance and one with consistently decreasing abundance (Fig. 5, clusters 1 and 4, respectively; for the detected masses, see Table SC in the supplemental material). All other clusters were much less coherent, as indicated by their quality scores. This coherence shows that the global metabolomic response to salt stress is
clearly dominant in our samples and involves a reasonably large number of 52 peak groups (33 increasing and 19 decreasing). Comparison of the time trends for all putatively identified metabolites from the salt shock cultures to those for the control, non-salt-shocked cultures revealed reproducible results for 15 metabolites accumulating in response to the salt shock and 3 metabolites with the reverse response (Table 1). The remaining metabolites either did not give reproducible results or showed equal behaviors in the salt shock and control cultures. Figure 6 shows exemplary time courses with and without salt shock, including the putative molecular identities of the compounds involved. A large number of the upregulated compounds were also contained in the medium used, and further analysis is required to determine whether they are synthesized or taken up by the organism. The most important accumulating
FIG. 5. Clustering analysis uncovers the salt stress response. Clustering analysis has been used to group metabolites showing similar temporal behaviors into clusters. For each of the 4 clusters, the mean behaviors are shown, with error bars depicting the standard deviations. Cluster 1 contains compounds that accumulate after the salt shock, including the well-known osmoprotectant proline. Cluster 4 contains compounds that decrease after the salt shock. Interestingly, these clusters are the most homogenous of all clusters, as indicated by the quality score, showing that the salt-responsive metabolites react as a single group.
DOWNLOADED FROM aem.ASM.ORG - at University of Groningen July 13, 2010
VOL. 76, 2010
KOL ET AL.
APPL. ENVIRON. MICROBIOL.
TABLE 1. Responses of putatively identified metabolites showing reproducible behavior in response to salt stress Response and class
Upregulation in response to salt shock Amino acids and derivatives
Peptides
Nucleotide derivatives
Downregulation in response to salt shock Amino acids and derivatives Nucleotide derivative No response to salt shock Amino acids and derivatives
Putative identification
Observed mass (Da)
Name
Formula
115.0633285 131.0946287 174.1116757 165.0789786 149.0510493 204.0898776 172.0847923 328.1859033 300.1797553 301.1273854 243.1219061 136.0385108 135.0544952 267.0967539 283.0916686
Proline (Iso)Leucine Arginine Phenylalanine Methionine Tryptophan GlyPro ArgGlyPro LysGlyPro GluGlyPro AlaGlyPro Hypoxanthine Adenine Adenosine/deoxyguanosine Guanosine
C5H9NO2 C6H13NO2 C6H14N4O2 C9H11NO2 C5H11NO2S C11H12N2O2 C7H12N2O3 C13H24N6O4 C13H24N4O4 C12H19N3O6 C10H17N3O4 C5H4N4O C5H5N5 C10H13N5O4 C10H13N5O5
117.0789786 147.0531578 297.0895601
Valine Glutamate 59-Methylthioadenosine
C5H11NO2 C5H9NO4 C11H15N5O3S
155.0694765 138.0429274 89.04767847
Histidine Urocanate Alanine
C6H9N3O2 C6H6N2O2 C3H7NO2
metabolite was proline, which is a well-known osmoprotectant that was previously reported for Streptomyces (11). The strong levels of accumulation of proline after the salt shock were the same for all of the strains, while in the nonshocked cultures proline levels remained stable or decreased slightly over time. The strong and immediate change indicates that proline is used for acute osmoprotection. Interestingly, this response is independent of the major regulators OsaB and sB. In addition to this classical salt-responsive metabolite, the global metabolome screen identified other potential osmoprotectants. Arginine, phenylalanine, methionine, tryprophan, and (iso)leucine showed the same strong accumulation as proline and thus could also play a supplementary role in acute osmoprotection. Like all amino acids, the listed compounds are zwitterions that could plausibly act as osmoprotectants. A similar dynamic, seen for guanosine, adenosine, adenine, and hypoxanthine nucleotide derivatives, is less obvious to explain. While affecting a large number of amino acids, the accumulation in response to the salt shock is not a general phenomenon for all amino acids, as can be seen in the time courses for glutamate and valine. In addition to the amino acids, we detected five di- and tripeptides that accumulated during salt stress (Fig. 4). All of these peptides contained a proline and a glycine residue. These proline/glycine-containing peptides were significantly enriched compared to the protein compositions for both S. coelicolor and Saccharomyces cerevisiae (P value of ,0.01, based on random sampling of peptide sequences). Glycine, glutamine, and alanine individually were observed at similar significance levels, while other amino acids were represented at the same frequencies in proteins and the accumulating peptides. This indicates that proline/glycine-containing peptides are either produced by specific proteolysis or taken up preferentially
from the peptone-containing medium. A possible candidate for a di- and tripeptide uptake system is SCO3064, which has homology with DtpT from Lactococcus lactis (18). This is in agreement with earlier studies of other organisms that also showed an involvement of proline-containing peptides in the salt stress response (1, 32, 37). One of the few putatively identified metabolites that showed a slight but consistent decrease in abundance during salt shock but not under control conditions was 59-methylthioadenosine. This metabolite acts as an inhibitor of polyamine biosynthesis in vitro and in vivo (19, 36). As polyamines are known to protect against salt stress in plants (48), active degradation of the inhibitor 59methylthioadenosine could be part of the salt response in Streptomyces as well. The expression of the pfs gene, encoding 59methylthioadenosine/S-adenosylhomocysteine nucleosidase (10), is increased by salt stress in E. coli (29). Alternatively, the decreased levels of 59-methylthioadenosine could indicate a decreased activity of the polyamine biosynthesis pathway under salt stress. Distinguishing these two alternatives requires targeted measurements of polyamines, which were not detected in our experimental setup. Other metabolites that decreased in abundance under highsalt conditions also did so under control conditions. For instance, the amino acid histidine and its catabolite urocanate decreased in both groups of cultures; these were the only amino acid-related compounds which showed this behavior. Not surprisingly, ectoine and hydroxyectoine were absent from the core salt shock response cluster, due to the disruption of the ectoine pathway in several of the strains. The detected compounds are thus those that show a consistent immediate response to salt shock independent of genotype, indicating that the salt shock metabolome extends far beyond the ectoines.
DOWNLOADED FROM aem.ASM.ORG - at University of Groningen July 13, 2010
2578
METABOLOMICS OF S. COELICOLOR
2579
FIG. 6. Comparison of time trends for putatively identified salt stress-responsive metabolites. The level of the major osmoprotectant proline, as well as those of several other putatively identified amino acids, shows an increase in response to salt shock. The detected proline/glycinecontaining di- and tripeptides also accumulate in response to salt shock (ArgGlyPro is shown as a representative example), as do a few nucleotide-related compounds, including adenosine. 59-Methylthioadenosine is one of the few compounds reproducibly showing the opposite pattern, an accumulation under low-salt conditions only. Levels of histidine and its catabolite urocanate decrease in all of the cultures, showing no response to the salt shock. For an overview of all 20 reproducible metabolites, see Fig. SE in the supplemental material. Results are shown for samples taken at 0, 4, 8, and 24 h after salt shock.
The only putative osmoprotectant that we detected but which showed no response to salt stress was alanine, which was previously described as an osmoprotectant for Streptomyces (11). The time series for all strains showed stable behavior of this metabolite for the salt-shocked samples and the non-saltshocked control cultures, indicating that it is in fact not used as an osmoprotectant under the conditions tested. Conclusions. Streptomyces species are known to be able to withstand considerable levels of osmotic stress (34) and to accumulate ectoine and hydroxyectoine under salt adaptation (7). We have studied the global metabolic response of S. coelicolor to salt stress by using a high-resolution LTQ Orbitrap
mass spectrometer. We did not aim for a complete coverage of the entire metabolome but rather intended to achieve an unbiased, comprehensive assessment of the moderately polar osmolytes most similar to ectoine, in particular, amino acids, nucleotides, and their derivatives. Sugars and polyols, the remaining major group of osmoprotectants (which were not detectable under the liquid chromatography-mass spectrometry conditions used), will therefore be the target of future studies using complementary analytical technologies with higher selectivity toward these compounds. It is very likely that among the described compounds, we have detected all that are osmotically relevant, as osmoprotectants need to accumulate in large
DOWNLOADED FROM aem.ASM.ORG - at University of Groningen July 13, 2010
VOL. 76, 2010
APPL. ENVIRON. MICROBIOL.
KOL ET AL.
amounts to be effective. However, the mass spectrometry data do not provide absolute quantification of the compounds involved. Instead, our results focus on the dynamic changes of the metabolome. This relative quantification, which is standard practice in metabolomics studies, clearly identifies the temporal trends for a large number of potentially osmoprotective compounds. The current study leads to the putative identification of potential novel osmoprotectants in S. coelicolor. Strikingly, we did not observe a major difference between the metabolomic responses to salt stress in wild-type cultures and in the two mutants with mutations in the major osmotic stress regulatory genes osaB and sigB (31). This indicates a considerable robustness of the metabolome, compared to the much more volatile transcriptome, in agreement with recent studies of plants (15). Robustness of metabolic flux under environmental and genetic perturbations has been found previously (13, 45) and seems to be a general feature of microbial metabolic networks (43). Our findings extend this observation to the robustness of steadystate metabolite levels. This robustness could be due to redundancy in the regulatory input leading to metabolic adaptation and has obvious evolutionary benefits. Elucidating the cellular networks that enable this metabolic robustness will require additional experimentation at multiple molecular levels, including transcriptomics, proteomics, and, most importantly, the salt-dependent kinetics of the enzymes involved. Furthermore, the strong accumulation of ectoine in the ectD::Tn strain is an indication that S. coelicolor could potentially be used as a producer of ectoine in (large-scale) bacterial milking (38), but the absolute levels of production would need to be quantified to determine whether this is an attractive option. ACKNOWLEDGMENTS We thank Erhard Bremer for his critical reading of the manuscript. The work was supported by grants from EU-FP6 ActinoGEN (to L.D. and S.K.), the Northern Netherlands Collaboration Initiative (SNN EZ/KOMPAS RM119 to L.D.), the University of Groningen (an Ubbo Emmius fellowship to M.E.M. and a Rosalind Franklin fellowship to E.T.), and the Netherlands Organization for Scientific Research NWO (a Vidi grant to R.B. and a medium investment grant to R.J.V.). REFERENCES 1. Amezaga, M. R., I. Davidson, D. McLaggan, A. Verheul, T. Abee, and I. R. Booth. 1995. The role of peptide metabolism in the growth of Listeria monocytogenes ATCC 23074 at high osmolarity. Microbiology 141:41–49. 2. Bentley, S. D., K. F. Chater, A. M. Cerdeno-Tarraga, G. L. Challis, N. R. Thomson, K. D. James, D. E. Harris, M. A. Quail, H. Kieser, D. Harper, A. Bateman, S. Brown, G. Chandra, C. W. Chen, M. Collins, A. Cronin, A. Fraser, A. Goble, J. Hidalgo, T. Hornsby, S. Howarth, C. H. Huang, T. Kieser, L. Larke, L. Murphy, K. Oliver, S. O’Neil, E. Rabbinowitsch, M. A. Rajandream, K. Rutherford, S. Rutter, K. Seeger, D. Saunders, S. Sharp, R. Squares, S. Squares, K. Taylor, T. Warren, A. Wietzorrek, J. Woodward, B. G. Barrell, J. Parkhill, and D. A. Hopwood. 2002. Complete genome sequence of the model actinomycete Streptomyces coelicolor A3(2). Nature 417:141–147. 3. Bishop, A., S. Fielding, P. Dyson, and P. Herron. 2004. Systematic insertional mutagenesis of a streptomycete genome: a link between osmoadaptation and antibiotic production. Genome Res. 14:893–900. 4. Borodina, I., P. Krabben, and J. Nielsen. 2005. Genome-scale analysis of Streptomyces coelicolor A3(2) metabolism. Genome Res. 15:820–829. 5. Breitling, R., A. R. Pitt, and M. P. Barrett. 2006. Precision mapping of the metabolome. Trends Biotechnol. 24:543–548. 6. Bursy, J., A. U. Kuhlmann, M. Pittelkow, H. Hartmann, M. Jebbar, A. J. Pierik, and E. Bremer. 2008. Synthesis and uptake of the compatible solutes ectoine and 5-hydroxyectoine by Streptomyces coelicolor A3(2) in response to salt and heat stresses. Appl. Environ. Microbiol. 74:7286–7296.
7. Bursy, J., A. J. Pierik, N. Pica, and E. Bremer. 2007. Osmotically induced synthesis of the compatible solute hydroxyectoine is mediated by an evolutionarily conserved ectoine hydroxylase. J. Biol. Chem. 282:31147–31155. 8. Canovas, D., C. Vargas, M. I. Calderon, A. Ventosa, and J. J. Nieto. 1998. Characterization of the genes for the biosynthesis of the compatible solute ectoine in the moderately halophilic bacterium Halomonas elongata DSM 3043. Syst. Appl. Microbiol. 21:487–497. 9. Christin, C., A. K. Smilde, H. C. Hoefsloot, F. Suits, R. Bischoff, and P. L. Horvatovich. 2008. Optimized time alignment algorithm for LC-MS data: correlation optimized warping using component detection algorithm-selected mass chromatograms. Anal. Chem. 80:7012–7021. 10. Cornell, K. A., and M. K. Riscoe. 1998. Cloning and expression of Escherichia coli 59-methylthioadenosine/S-adenosylhomocysteine nucleosidase: identification of the pfs gene product. Biochim. Biophys. Acta 1396:8–14. 11. Empadinhas, N., and M. S. da Costa. 2008. Osmoadaptation mechanisms in prokaryotes: distribution of compatible solutes. Int. Microbiol. 11:151–161. 12. Fahy, E., M. Sud, D. Cotter, and S. Subramaniam. 2007. LIPID MAPS online tools for lipid research. Nucleic Acids Res. 35:W606–W612. 13. Fischer, E., and U. Sauer. 2005. Large-scale in vivo flux analysis shows rigidity and suboptimal performance of Bacillus subtilis metabolism. Nat. Genet. 37:636–640. 14. Flett, F., V. Mersinias, and C. P. Smith. 1997. High efficiency intergeneric conjugal transfer of plasmid DNA from Escherichia coli to methyl DNArestricting streptomycetes. FEMS Microbiol. Lett. 155:223–229. 15. Fu, J., J. J. Keurentjes, H. Bouwmeester, T. America, F. W. Verstappen, J. L. Ward, M. H. Beale, R. C. de Vos, M. Dijkstra, R. A. Scheltema, F. Johannes, M. Koornneef, D. Vreugdenhil, R. Breitling, and R. C. Jansen. 2009. Systemwide molecular evidence for phenotypic buffering in Arabidopsis. Nat. Genet. 41:166–167. 16. Galinski, E. A. 1995. Osmoadaptation in bacteria. Adv. Microb. Physiol. 37:272–328. 17. Garcia-Estepa, R., M. Argandona, M. Reina-Bueno, N. Capote, F. IglesiasGuerra, J. J. Nieto, and C. Vargas. 2006. The ectD gene, which is involved in the synthesis of the compatible solute hydroxyectoine, is essential for thermoprotection of the halophilic bacterium Chromohalobacter salexigens. J. Bacteriol. 188:3774–3784. 18. Hagting, A., J. v. d. Velde, B. Poolman, and W. N. Konings. 1997. Membrane topology of the di- and tripeptide transport protein of Lactococcus lactis. Biochemistry 36:6777–6785. 19. Hibasami, H., and A. E. Pegg. 1978. Differential inhibition of mammalian aminopropyltransferase activities. Biochem. Biophys. Res. Commun. 81: 1398–1405. 20. Hopwood, D. A. 1967. Genetic analysis and genome structure in Streptomyces coelicolor. Bacteriol. Rev. 31:373–403. 21. Hopwood, D. A. 2007. Streptomyces in nature and medicine: the antibiotic makers. Oxford University Press, Oxford, United Kingdom. 22. Hu, Q., R. J. Noll, H. Li, A. Makarov, M. Hardman, and R. Graham Cooks. 2005. The Orbitrap: a new mass spectrometer. J. Mass Spectrom. 40:430– 443. 23. Kamleh, M. A., J. A. Dow, and D. G. Watson. 2009. Applications of mass spectrometry in metabolomic studies of animal model and invertebrate systems. Brief. Funct. Genomic. Proteomic. 8:28–48. 24. Kanehisa, M., and S. Goto. 2000. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28:27–30. 25. Keller, B. O., J. Sui, A. B. Young, and R. M. Whittal. 2008. Interferences and contaminants encountered in modern mass spectrometry. Anal. Chim. Acta 627:71–81. 26. Kempf, B., and E. Bremer. 1998. Uptake and synthesis of compatible solutes as microbial stress responses to high-osmolality environments. Arch. Microbiol. 170:319–330. 27. Kieser, T., M. J. Bibb, M. J. Butter, K. F. Chater, and D. A. Hopwood. 2000. Practical Streptomyces genetics. The John Innes Foundation, Norwich, United Kingdom. 28. Killham, K., and M. K. Firestone. 1984. Salt stress control of intracellular solutes in streptomycetes indigenous to saline soils. Appl. Environ. Microbiol. 47:301–306. 29. Kim, Y., C. M. Lew, and J. D. Gralla. 2006. Escherichia coli pfs transcription: regulation and proposed roles in autoinducer-2 synthesis and purine excretion. J. Bacteriol. 188:7457–7463. 30. Lee, E. J., Y. H. Cho, H. S. Kim, B. E. Ahn, and J. H. Roe. 2004. Regulation of sigmaB by an anti- and an anti-anti-sigma factor in Streptomyces coelicolor in response to osmotic stress. J. Bacteriol. 186:8490–8498. 31. Lee, E. J., N. Karoonuthaisiri, H. S. Kim, J. H. Park, C. J. Cha, C. M. Kao, and J. H. Roe. 2005. A master regulator sigmaB governs osmotic and oxidative response as well as differentiation via a network of sigma factors in Streptomyces coelicolor. Mol. Microbiol. 57:1252–1264. 32. Le Marrec, C., E. Bon, and A. Lonvaud-Funel. 2007. Tolerance to high osmolality of the lactic acid bacterium Oenococcus oeni and identification of potential osmoprotectants. Int. J. Food Microbiol. 115:335–342. 33. Louis, P., and E. A. Galinski. 1997. Characterization of genes for the biosynthesis of the compatible solute ectoine from Marinococcus halophilus and osmoregulated expression in Escherichia coli. Microbiology 143:1141–1149.
DOWNLOADED FROM aem.ASM.ORG - at University of Groningen July 13, 2010
2580
34. Malin, G., and A. Lapidot. 1996. Induction of synthesis of tetrahydropyrimidine derivatives in Streptomyces strains and their effect on Escherichia coli in response to osmotic and heat stress. J. Bacteriol. 178:385–395. 35. Ono, H., K. Sawada, N. Khunajakr, T. Tao, M. Yamamoto, M. Hiramoto, A. Shinmyo, M. Takano, and Y. Murooka. 1999. Characterization of biosynthetic enzymes for ectoine as a compatible solute in a moderately halophilic eubacterium, Halomonas elongata. J. Bacteriol. 181:91–99. 36. Pegg, A. E., R. T. Borchardt, and J. K. Coward. 1981. Effects of inhibitors of spermidine and spermine synthesis on polyamine concentrations and growth of transformed mouse fibroblasts. Biochem. J. 194:79–89. 37. Piuri, M., C. Sanchez-Rivas, and S. M. Ruzal. 2003. Adaptation to high salt in Lactobacillus: role of peptides and proteolytic enzymes. J. Appl. Microbiol. 95:372–379. 38. Sauer, T., and E. A. Galinski. 1998. Bacterial milking: a novel bioprocess for production of compatible solutes. Biotechnol. Bioeng. 57:306–313. 39. Scalbert, A., L. Brennan, O. Fiehn, T. Hankemeier, B. Kristal, B. van Ommen, E. Pujos-Guillot, E. Verheij, D. Wishart, and S. Wopereis. 2009. Mass-spectrometry-based metabolomics: limitations and recommendations for future progress with particular focus on nutrition research. Metabolomics 5:435–458. 40. Scheltema, R. A., A. Kamleh, D. Wildridge, C. Ebikeme, D. G. Watson, M. P. Barrett, R. C. Jansen, and R. Breitling. 2008. Increasing the mass accuracy of high-resolution LC-MS data using background ions: a case study on the LTQ-Orbitrap. Proteomics 8:4647–4656. 41. Shah, V. P., K. K. Midha, J. W. Findlay, H. M. Hill, J. D. Hulse, I. J. McGilveray, G. McKay, K. J. Miller, R. N. Patnaik, M. L. Powell, A. Tonelli, C. T. Viswanathan, and A. Yacobi. 2000. Bioanalytical method validation—a revisit with a decade of progress. Pharm. Res. 17:1551–1557.
METABOLOMICS OF S. COELICOLOR
2581
42. Smith, C. A., G. O’Maille, E. J. Want, C. Qin, S. A. Trauger, T. R. Brandon, D. E. Custodio, R. Abagyan, and G. Siuzdak. 2005. METLIN: a metabolite mass spectral database. Ther. Drug Monit. 27:747–751. 43. Stelling, J., U. Sauer, Z. Szallasi, F. J. Doyle III, and J. Doyle. 2004. Robustness of cellular functions. Cell 118:675–685. 44. Strauch, E., E. Takano, H. A. Baylis, and M. J. Bibb. 1991. The stringent response in Streptomyces coelicolor A3(2). Mol. Microbiol. 5:289–298. 45. Tang, Y. J., H. G. Martin, P. S. Dehal, A. Deutschbauer, X. Llora, A. Meadows, A. Arkin, and J. D. Keasling. 2009. Metabolic flux analysis of Shewanella spp. reveals evolutionary robustness in central carbon metabolism. Biotechnol. Bioeng. 102:1161–1169. 46. Tautenhahn, R., C. Bo ¨ttcher, and S. Neumann. 2007. Annotation of LC/ ESI-MS mass signals, p. 371–380. In S. Hochreiter and R. Wagner (ed.), Bioinformatics research and development. Springer-Verlag, Berlin, Germany. 47. Wishart, D. S., D. Tzur, C. Knox, R. Eisner, A. C. Guo, N. Young, D. Cheng, K. Jewell, D. Arndt, S. Sawhney, C. Fung, L. Nikolai, M. Lewis, M. A. Coutouly, I. Forsythe, P. Tang, S. Shrivastava, K. Jeroncic, P. Stothard, G. Amegbey, D. Block, D. D. Hau, J. Wagner, J. Miniaci, M. Clements, M. Gebremedhin, N. Guo, Y. Zhang, G. E. Duggan, G. D. Macinnis, A. M. Weljie, R. Dowlatabadi, F. Bamforth, D. Clive, R. Greiner, L. Li, T. Marrie, B. D. Sykes, H. J. Vogel, and L. Querengesser. 2007. HMDB: the Human Metabolome Database. Nucleic Acids Res. 35:D521–D526. 48. Yamaguchi, K., Y. Takahashi, T. Berberich, A. Imai, A. Miyazaki, T. Takahashi, A. Michael, and T. Kusano. 2006. The polyamine spermine protects against high salt stress in Arabidopsis thaliana. FEBS Lett. 580: 6783–6788.
DOWNLOADED FROM aem.ASM.ORG - at University of Groningen July 13, 2010
VOL. 76, 2010