|
|
||||||||
1 Department of Chemical Engineering, University of Washington, Seattle, WA 98195, USA
2 Department of Microbiology, University of Washington, Seattle, WA 98195, USA
Correspondence
Ludmila Chistoserdova
milachis{at}u.washington.edu
| ABSTRACT |
|---|
|
|
|---|
A supplementary table is available with the online version of this paper.
| INTRODUCTION |
|---|
|
|
|---|
M. mobilis JLW8 is the only formally described representative of Methylotenera, a new genus within the family Methylophilaceae (Kalyuzhnaya et al., 2006
). All known members of this family are obligate methylotrophs, i.e. they are specialized in degradation of organic compounds containing no carbon–carbon bonds (C1 compounds), most prominently methanol and methylamine (Lidstrom, 2006
). Based on culture-independent surveys, Methylotenera species appear to be ubiquitous, having been detected in marine, freshwater, acid mine drainage, polluted soil, sedimentary rock and glacier environments (according to the entries in the non-redundant database, http://www.ncbi.nlm.nih.gov/), consistent with an important role in global carbon cycling. However, the exact function of these bacteria remains unknown. M. mobilis JLW8 was isolated from a 63 m deep sampling site in Lake Washington sediment (Kalyuzhnaya et al., 2006
). Samples used for stable isotope probing-based enrichment of methylotrophic communities that were also used for metagenomic sequencing were collected from the same site, but at a later date (Kalyuzhnaya et al., 2008b
). From metagenomic analysis conducted on these samples, the methylamine microcosm was dominated by M. mobilis strains, which allowed for extracting a composite genome of M. mobilis and conducting metabolic reconstruction for this species (Kalyuzhnaya et al., 2008b
). However, the dominant strains of M. mobilis represented in the metagenome appear to be distinct from M. mobilis JLW8, based on 16S rRNA sequence comparisons. The study undertaken here pursued three main objectives: (1) to validate the metagenome-based sequence assembly and binning of M. mobilis sequences conducted in the previous study, (2) to test whether this metagenomic sequence provides a quality protein database for peptide identification in a strain without perfect matches in the metagenome and test the limitation of such a database, and (3) to obtain new insights into the physiology of M. mobilis.
| METHODS |
|---|
|
|
|---|
Protein extraction and digestion.
Frozen cell pellets were resuspended in 500 µl hot resuspension buffer (20 mM Tris/HCl pH 8.0, 5 mM DTT) and lysed by boiling for 2 min in a water bath, followed by cooling on ice for 10 min. For nucleic acid digestion, 10 units benzonase nuclease (Roche) was added after adjusting the suspension to 2 mM MgCl2 (final concentration). After 15 min incubation at room temperature, 200 µl of ethanol- and buffer-washed glass beads (150 µm) was added, the sample volume was adjusted to 1.5 ml with resuspension buffer, and bead beating was performed for 4 min at 48 r.p.m. in a MiniBeadBeater (BioSpec Products). Total protein concentration was determined by Bradford protein assay (Bio-Rad); this was approximately 1 mg ml–1. The homogenate was lyophilized to dryness. After resuspension in 500 µl digestion buffer (2 M urea, 5 % acetonitrile, 5 mM dithiothreitol, 0.1 % Rapigest), the sample was digested as previously described (Bosch et al., 2008
).
HPLC pre-fractionation and linear ion trap mass spectrometry.
The soluble fraction after digestion was lyophilized to a volume of approximately 150 µl and centrifuged for 10 min at 14 000 r.p.m. in a tabletop centrifuge (Eppendorf). A 10–20 µl sample of the supernatant was applied to a PLRP-S reversed-phase column (2.1 mm i.d.x150 mm, 300 Å, 5 µm; Polymer Laboratories) with mobile phases of 0.1 % trifluoroacetic acid in H2O and acetonitrile. Peptides were eluted at 0.2 ml min–1 with a gradient of 2–60 % acetonitrile in 60 min and 60–90 % acetonitrile in 20 min and collected as five separate fractions. The fractions were each lyophilized to 20 µl, and after reconstitution to volumes of 120 µl with acetic acid and acetonitrile at final concentrations of 0.5 % and 5 % (v/v), respectively, subjected to LC/LC-MS/MS (two-dimensional capillary HPLC/linear ion trap tandem mass spectrometry) analysis as previously described (Bosch et al., 2008
).
SEQUEST database searching, DTAselect filtering and estimation of random false positive identification rate.
Using the SEQUEST algorithm, raw data were searched against a concatenated FASTA database that consisted of polypeptide sequences derived from the M. mobilis composite genome (Kalyuzhnaya et al., 2008b
), the human subset of the non-redundant database (nrdb) depleted of all virus sequences, and reversed sequence (decoy) versions of both databases appended. The concatenated database comprised a total of 214 148 protein sequences, 12 951 of which represented the M. mobilis specific meta-database. The search results were filtered in DTAselect applying the parameters described by Bosch et al. (2008)
. Sequest
Cn/Xcorr values of 0.08/1.9, 0.08/2.2 and 0.08/3.3 were used for singly, doubly and triply charged peptides, respectively. Three peptides unique to a particular ORF were required for positive identification. Random false positive identifications were determined at the peptide level as a ratio of two times the number of reversed peptide identifications in the concatenated database to the total number of identified peptides (Elias & Gygi, 2007
). By this means, the qualitative random false positive rate was estimated to be 1.5 % for the first biological replicate and 2.1 % for the second biological replicate.
Data processing and normalization.
The in-house methods developed for estimating protein abundances from spectral counting values have been described by Xia et al. (2007a
, b
). In order to normalize the two biological replicates, the total sum of the average spectral counts for the second biological replicate was divided by the total sum of the average spectral counts for the first biological replicate, resulting in a normalization factor of 1.077. Each average spectral count value in the first biological replicate was multiplied by this normalization factor.
Reproducibility of biological replicates.
Linear regression analysis demonstrated good data reproducibility for the two biological replicates (R=0.86 for the 1924 proteins common to both replicates; data not shown).
DNA amplification, sequencing and analysis.
DNA of M. mobilis JLW8 was isolated using a QIAamp DNA minikit (Qiagen). Semi-randomly selected fragments of DNA, as shown in Table 1
, were amplified by PCR, cloned and sequenced essentially as described by Kalyuzhnaya et al. (2008a)
. Oligonucleotide primers for amplification were designed based on the metagenomic sequence, as follows. First, a gene to be amplified was chosen (genes involved in methylotrophy, central pathways or hypothetical genes). All the homologues representing different M. mobilis strains in the metagenome were then aligned using tools available as part of the Joint Genome Institute IMG/M package (http://img.jgi.doe.gov/cgi-bin/m/main.cgi), and consensus sequences were chosen with no or few mismatches.
|
Enzyme assays.
Cells were grown on methylamine as above, to OD600 0.8, pelleted by centrifugation, as above, resuspended in Tris/HCl pH 8.0 buffer and disrupted by passage through a French pressure cell at 1.2x108 Pa. Cell extracts were centrifuged at 14 000 r.p.m. for 25 min at 4 °C to remove cell debris. The activity of glutamine synthetase was assayed as described by Murrell & Dalton (1983)
, the activity of glutamate synthase was assayed as described by Bravo & Mora (1988)
, and the activity of glutamate dehydrogenase was assayed as described by Yarrison et al. (1972)
.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
Qualitative proteome coverage
Spectral counts as used here are summations of the number of redundant peptides (i.e. including repeated measurements of the same peptide at different stages of the HPLC separation) observed that map uniquely to a particular protein or protein grouping as defined by the user (Bosch et al., 2008
). While the precise relationship between observed spectral counts and abundance within a given proteome is complex, for microbial systems high counts usually correlate with high abundance, as validated by global transcription measurements and other means (Bosch et al., 2008
). However, low counts do not necessarily imply low abundance, due to detectability differences among peptides. Data for two biological replicates were collected and analysed. The M. mobilis specific protein database utilized for peptide identification was represented by polypeptide sequences translated from the 12 951 gene sequences that constituted the M. mobilis composite genome (covering approximately five closely related M. mobilis genomes; Kalyuzhnaya et al., 2008b
), downloaded from the IMG/M database, along with their automated annotations. These annotations were manually curated for proteins with high spectral counts and thus inferred high relative abundance (see Supplementary Table S1, available with the online version of this paper). A total of 2353 and 2336 polypeptides in this database were matched with the peptides identified in the first and the second biological replicate samples, respectively (Table S1), with a positive detection cutoff of at least three peptide matches per protein (at a cutoff of two peptide matches, 2767 and 2787 peptides were identified, respectively; data not shown). This identification rate (approx. 20 % of total proteins in the database) is, as expected, in a lower range compared to the identification rate typical of proteins matched to a database specific for a given organism (up to 60 %, for example, in M. extorquens AM1; Bosch et al., 2008
), but within the range of the previously conducted metagenome-based proteomic analyses (Ram et al., 2005
) and in agreement with the recent models for protein detection efficiency as a function of sequence divergence (Denef et al., 2007
). An additional factor affecting the peptide detection was the partial nature of many of the polypeptides translated from the M. mobilis metagenome, due to low sequence assembly (Kalyuzhnaya et al., 2008b
).
Protein-based reconstruction of methylotrophy-specific metabolic pathways
In most cases, peptides of M. mobilis JLW8 were matched to a number of homologous proteins in the database, representing protein clusters, similarly to the results represented in Table 1
. Those homologues with positive peptide matches were grouped manually into clusters, and peptide spectral counts within each cluster were summed. Some of the inferred proteins with the highest summed spectral counts observed in both biological replicates were the ones usually found highly expressed in bacteria, such as the translation elongation factor, chaperonins GroEL and DnaK, subunits of ATP synthase, DNA-directed RNA polymerase, ribosomal proteins, etc. (Table S1) that can be defined as core proteins (Callister et al., 2008
). However, of most interest to us were the peptides representing major enzymes and pathways specifically related to methylamine metabolism in M. mobilis. One of the most highly expressed metabolic enzymes was hexulosephosphate synthase (Hps, phylotype 1; Table 2
, Table S1 and Fig. 1
), the first enzyme in the ribulose monophosphate (RMP) cycle operating in formaldehyde assimilation/oxidation in Methylophilaceae (Lidstrom, 2006
). The gene for this enzyme is part of a methylotrophy gene island in M. mobilis (Kalyuzhnaya et al., 2008b
) that is very similar to an island also found in Methylobacillus flagellatus (Chistoserdova et al., 2007b
). Interestingly, both Methylotenera mobilis and Methylobacillus flagellatus encode a second copy of Hps (phylotype 2, approximately 89 % identity with phylotype 1 in M. mobilis). This enzyme is also detectable in the M. mobilis proteome, but at a much lower level of total counts (Table 1
). Other enzymes that are parts of the RMP cycle, both assimilatory and dissimilatory branches, were also identified at relatively high spectral counts, the most abundant being transaldolase and transketolase, which are involved in RMP regeneration (Fig. 1
). Peptides for the catalytic subunits of methylamine dehydrogenase (MauAB) were found at high counts, along with the proposed novel electron-accepting cytochrome (MauO; Kalyuzhnaya et al., 2008b
), suggesting that this is the main primary oxidation system. Interestingly, MauD, whose function has been proposed as being in small-subunit maturation (van der Palen et al., 1997
), was found at relatively high counts, unlike the proteins involved in cofactor biosynthesis (for example, MauG; Li et al., 2008
), suggesting a possible role as a chaperone. Polypeptides were detected at high counts forming the glutamine synthetase/glutamate synthase (GS/GOGAT) pathway. While not typically considered as part of the methylamine utilization pathway, the presence of these enzymes suggests a mechanism for utilizing the ammonia that results from methylamine oxidation. In this respect, the incomplete tricarboxylic acid cycle (Kalyuzhnaya et al., 2008b
), in addition to its anaplerotic function in providing central metabolism intermediates (Lidstrom, 2006
), must play another important role, in providing the
-ketoglutarate that is an essential intermediate in the GS/GOGAT cycle. The presence of the activities of methylamine dehydrogenase and the hexulosephosphate synthase/isomerase pair in M. mobilis JLW8 has been previously reported, supporting the analysis conducted here (Kalyuzhnaya et al., 2006
). To confirm identification of the GS/GOGAT cycle, we measured the activities of the respective enzymes and found them at high levels (95±3 and 174±5 mU, respectively). Glutamate dehydrogenase (GDH), another enzyme known to assimilate ammonia, was undetectable in cell extracts, in agreement with the lack of a gene homologue encoding GDH in the M. mobilis genomic database.
|
|
One metabolic feature that distinguishes Methylotenera mobilis from Methylobacillus flagellatus, as predicted from genome–genome comparisons, is the presence of the genes for a complete methylcitric acid (MCA) cycle. All the proteins constituting this pathway were positively identified in this study (Table 2
, Supplementary Table S1 and Fig. 1
). The role of this pathway in M. mobilis remains unknown. Its involvement in metabolism of methylamine is not likely as all the essential metabolites are produced in the RMP cycle in combination with the enzymes interconverting C3 metabolites and the partial tricarboxylic acid cycle (Fig. 1
). More likely, the presence of the specific MCA cycle proteins during growth on methylamine is due to the peculiarities of regulation of this metabolic pathway. It has been previously suggested that this pathway may play a role in utilization of methylated compounds whose degradation results in propionate as a product (Kalyuzhnaya et al., 2008b
). While this hypothesis is a subject for experimental testing, the present study demonstrates expression of the entire set of the genes involved, thus supporting functionality of this pathway.
Curiously, XoxF and XoxG, homologues of, respectively, the large subunit of methanol dehydrogenase and an associated cytochrome, were found at high abundance. The potential role of the Xox system has been attracting attention recently, implicating it in C1 metabolism, but so far without a well-defined function (Wilson et al., 2008
; Kalyuzhnaya et al., 2008a
). Another highly abundant enzyme was a putative aldehyde dehydrogenase (Adh) whose function also remains unknown.
Conclusions
We utilized a metagenomic sequence of M. mobilis for a strain-specific proteomic analysis of the closely related M. mobilis JLW8 grown on methylamine. To address the first specific objective (see Introduction), we provided an experimental validation of sequence assembly and M. mobilis-specific sequence binning conducted in our previous study (Kalyuzhnaya et al., 2008b
), by detecting proteins involved in major metabolic pathways predicted to operate in metabolism of methylamine. To address the second specific objective, we demonstrated that a metagenomic dataset of inherently lower quality compared to completed genomic sequences is usable for cross-strain protein detection, by identifying a large portion of polypeptides predicted to be encoded in the genome of M. mobilis JLW8, with high confidence. To address the third specific objective, we obtained new insights into the metabolism of methylamine by M. mobilis, including expression of multiple pathways for formaldehyde oxidation (linear versus cyclic) and multiple enzyme paralogues (Hps, Fae), suggesting redundancy of key methylotrophy functions previously noted for other organisms (Chistoserdova et al., 2000
, 2007a
). In addition, we established a potential link between methylamine oxidation and ammonia utilization via the GS/GOGAT pathway. This pathway is likely to play a role in the removal of ammonia produced by the methylamine dehydrogenase reaction to prevent its accumulation to toxic levels. We also demonstrated that all the proteins involved in the proposed MCA cycle were expressed, suggesting a central role for this cycle in M. mobilis metabolism. This work is a prelude to future applications of this metagenomic dataset for protein identification in other M. mobilis strains as well as in natural microbial populations (metagenome-based metaproteomics).
| ACKNOWLEDGEMENTS |
|---|
Edited by: J. A. Vorholt
| REFERENCES |
|---|
|
|
|---|
Bravo, A. & Mora, J. (1988). Ammonium assimilation in Rhizobium phaseoli by the glutamine synthetase-glutamate synthase pathway. J Bacteriol 170, 980–984.
Callister, S. J., McCue, L. A., Turse, J. E., Monroe, M. E., Auberry, K. J., Smith, R. D., Adkins, J. N. & Lipton, M. S. (2008). Comparative bacterial proteomics: analysis of the core genome concept. PLoS One 3, e1542[CrossRef][Medline]
Chistoserdova, L., Vorholt, J. A., Thauer, R. K. & Lidstrom, M. E. (1998). C1 transfer enzymes and coenzymes linking methylotrophic bacteria and methanogenic Archaea. Science 281, 99–102.
Chistoserdova, L., Gomelsky, L., Vorholt, J. A., Gomelsky, M., Tsygankov, Y. D., Thauer, R. K. & Lidstrom, M. E. (2000). Analysis of two formaldehyde oxidation pathways in Methylobacillus flagellatus KT, a ribulose monophosphate cycle methylotroph. Microbiology 146, 233–238.
Chistoserdova, L., Chen, S.-W., Lapidus, A. & Lidstrom, M. E. (2003). Methylotrophy in Methylobacterium extorquens AM1 from a genomic point of view. J Bacteriol 185, 2980–2987.
Chistoserdova, L., Crowther, G. J., Vorholt, J. A., Skovran, B., Portais, J.-C. & Lidstrom, M. E. (2007a.). Identification of a fourth formate dehydrogenase in Methylobacterium extorquens AM1 and confirmation of the essential role of formate oxidation in methylotrophy. J Bacteriol 189, 9076–9081.
Chistoserdova, L., Lapidus, A., Han, C., Goodwin, L., Saunders, L., Brettin, T., Tapia, R., Gilna, P., Lucas, S. & other authors (2007b). The genome of Methylobacillus flagellatus, the molecular basis for obligate methylotrophy, and the polyphyletic origin of methylotrophy. J Bacteriol 189, 4020–4027.
Denef, V. J., Shah, M. B., Verberkmoes, N. C., Hettich, R. L. & Banfield, J. F. (2007). Implications of strain- and species-level sequence divergence for community and isolate shotgun proteomic analysis. J Proteome Res 6, 3152–3161.[CrossRef][Medline]
Elias, J. E. & Gygi, S. P. (2007). Target-decoy search strategy for increased confidence in large-scale protein identifications by mass spectrometry. Nat Methods 4, 207–214.[CrossRef][Medline]
Frias-Lopez, J., Shi, Y., Tyson, G. W., Coleman, M. L., Schuster, S. C., Chisholm, S. W. & DeLong, E. F. (2008). Microbial community gene expression in ocean surface waters. Proc Natl Acad Sci U S A 105, 3805–3810.
Gasteiger, E., Hoogland, C., Gattiker, A., Duvaud, S., Wilkins, M. R., Appel, R. D. & Bairoch, A. (2005). Protein identification and analysis tools on the ExPASy server. In The Proteomics Protocols Handbook, pp. 571–607. Edited by John M. Walker. Totowa, NJ: Humana Press
Gilbert, J. A., Field, D., Huang, Y., Edwards, R., Li, W., Gilna, P. & Joint, I. (2008). Detection of large numbers of novel sequences in the metatranscriptomes of complex marine microbial communities. PLoS One 3, e3042[CrossRef][Medline]
Kalyuzhnaya, M. G., Korotkova, N., Crowther, G. J., Marx, C. J., Lidstrom, M. & Chistoserdova, L. (2005). Analysis of gene islands involved in methanopterin-linked C1 transfer reactions reveals new functions and provides evolutionary insights. J Bacteriol 187, 4607–4614.
Kalyuzhnaya, M. G., Bowerman, S., Lara, J. C., Lidstrom, M. E. & Chistoserdova, L. (2006). Methylotenera mobilis gen. nov., sp. nov, an obligately methylamine-utilizing bacterium within the family Methylophilaceae. Int J Syst Evol Microbiol 56, 2819–2823.
Kalyuzhnaya, M. G., Hristova, K. R., Lidstrom, M. E. & Chistoserdova, L. (2008a). Characterization of a novel methanol dehydrogenase in representatives of Burkholderiales: implications for environmental detection of methylotrophy and evidence for convergent evolution. J Bacteriol 190, 3817–3823.
Kalyuzhnaya, M. G., Lapidus, A., Ivanova, N., Copeland, A. C., McHardy, A., Szeto, E., Salamov, A., Grigoriev, I. V., Suciu, D. & other authors (2008b). High resolution metagenomics targets major functional types in complex microbial communities. Nat Biotechnol 26, 1029–1034.[CrossRef][Medline]
Konstantinidis, K. T. & Tiedje, J. M. (2005). Towards a genome-based taxonomy for prokaryotes. J Bacteriol 187, 6258–6264.
Li, X., Fu, R., Liu, A. & Davidson, V. L. (2008). Kinetic and physical evidence that the diheme enzyme MauG tightly binds to a biosynthetic precursor of methylamine dehydrogenase with incompletely formed tryptophan tryptophylquinone. Biochemistry 47, 2908–2912.[CrossRef][Medline]
Lidstrom, M. E. (2006). Aerobic methylotrophic prokaryotes. In The Prokaryotes. Edited by A. Balows, H. G. Truper, M. Dworkin, W. Harder & K.-H. Schleifer. New York: Springer.
Murrell, J. C. & Dalton, H. (1983). Purification and properties of glutamine synthetase from Methylococcus capsulatus (Bath). J Gen Microbiol 129, 1187–1196.
Ram, R. J., Verberkmoes, N. C., Thelen, M. P., Tyson, G. W., Baker, B. J., Blake, R. C., II, Shah, M., Hettich, R. L. & Banfield, J. F. (2005). Community proteomics of a natural microbial biofilm. Science 308, 1915–1920.
Rappé, M. S. & Giovannoni, S. J. (2003). The uncultured microbial majority. Annu Rev Microbiol 57, 369–394.[CrossRef][Medline]
Rusch, D. B., Halpern, A. L., Sutton, G., Heidelberg, K. B., Williamson, S., Yooseph, S., Wu, D., Eisen, J. A., Hoffman, J. M. & other authors (2007). The Sorcerer II Global Ocean Sampling Expedition: Northwest Atlantic through Eastern Tropical Pacific. PLoS Biol 5, e77[CrossRef][Medline]
Tringe, S. G., von Mering, C., Kobayashi, A., Salamov, A. A., Chen, A., Chang, H. W., Podar, M., Short, J. M., Mathur, A. J. & other authors (2005). Comparative metagenomics of microbial communities. Science 308, 554–557.
van der Palen, C. J., Reijnders, W. N., de Vries, S., Duine, J. A. & Spaning, R. J. (1997). MauE and MauD proteins are essential in methylamine metabolism of Paracoccus denitrificans. Antonie Van Leeuwenhoek 72, 219–228.[CrossRef][Medline]
Vorholt, J. A., Marx, C. J., Lidstrom, M. E. & Thauer, R. K. (2000). Novel formaldehyde-activating enzyme in Methylobacterium extorquens AM1 required for growth on methanol. J Bacteriol 182, 6645–6650.
Wilson, S. M., Gleisten, M. P. & Donohue, T. J. (2008). Identification of proteins involved in formaldehyde metabolism by Rhodobacter sphaeroides. Microbiology 154, 296–305.
Xia, Q., Wang, T., Taub, F., Park, Y., Capestany, C. A., Lamont, R. J. & Hackett, M. (2007a). Quantitative proteomics of intracellular Porphyromonas gingivalis. Proteomics 7, 4323–4337.[CrossRef][Medline]
Xia, Q., Wang, T., Park, Y., Lamont, R. J. & Hackett, M. (2007b). Differential quantitative proteomics of Porphyromonas gingivalis by linear ion trap mass spectrometry: non-label methods comparison, q-values and LOWESS curve fitting. Int J Mass Spectrom 259, 105–116.[CrossRef][Medline]
Yarrison, G., Young, D. W. & Choules, G. L. (1972). Glutamate dehydrogenase from Mycoplasma laidlawii. J Bacteriol 110, 494–503.
Received 9 October 2008;
revised 2 January 2009;
accepted 5 January 2009.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |