|
|
||||||||
1 The Royal Veterinary College, Royal College Street, London NW1 0TU, UK
2 Max-Planck-Institute for Infection Biology, Core Facility Protein Analysis, Campus Charité Mitte, Schumannstr. 21/22, D-10117 Berlin, Germany
Correspondence
Neil G. Stoker
nstoker{at}rvc.ac.uk
| ABSTRACT |
|---|
|
|
|---|
Sample outputs of the programs AlternaStart.pl, proteogest and ParseProteogest.pl are available as supplementary data with the online version of this paper.
| INTRODUCTION |
|---|
|
|
|---|
Determining the TSS accurately is important, not only because this defines the amino acid sequence of the protein (as the stop codon is unambiguous), but also because this defines the upstream region in the DNA. Genome-wide and focused studies of promoter structure and regulatory motifs depend on the intergenic regions defined by the gene-finding process (Edwards et al., 2005
; Salgado et al., 2000
). The TSS assignment therefore affects the analysis of both protein function and transcriptional regulation.
Traditionally, TSS identification has been achieved using N-terminal sequencing by Edman degradation (Edman, 1950
). This is often technically demanding, and requires large (picomole) quantities of protein. Furthermore, some proteins are blocked through N-terminal modifications and cannot be sequenced. Proteomic approaches using sensitive MS methods (Jungblut et al., 1999
, 2000
) have revolutionized protein analysis because of the speed and sensitivity of the technology, but are not used to identify TSSs.
In this paper we apply rapid proteomic analysis methods to the issue of identifying TSSs. As the same data generated for protein identification could simultaneously identify protein starts, this is therefore a highly efficient approach. We were able to confirm the predicted N termini of 11 proteins, and to correct the predictions for two proteins from Mycobacterium tuberculosis. We also demonstrated that in contrast to the predicted situation in Escherichia coli, a high proportion (23 %) of the proteins showed N-terminal acetylation.
| METHODS |
|---|
|
|
|---|
|
Protein test set selection.
In order to assess the feasibility of our strategy, we decided to focus initially on the proteins most likely to have their definitive translation start identified by MS. Thus, the alternative protein sequences were screened according to a number of technical criteria. First we discarded any protein not in a set of 289 proteins for which a 2D-gel spot had been identified in a previous large-scale proteomics analysis of M. tuberculosis (http://web.mpiib-berlin.mpg.de/cgi-bin/pdbs/2d-page/extern/index.cgi). We next eliminated any protein for which the predicted mass of the first N-terminal tryptic fragment of the p0 protein was not between 800 and 2000 Da (the best resolution range for our MS equipment), and proteins without an arginine at the C-terminal end of the p0 tryptic peptide [as ionization in MALDI-MS is less effective in lysine-containing peptides than in arginine-containing ones (Krause et al., 1999
)]. The remaining 76 proteins were then ranked according to the minimal number of alternative starts (i.e. least total number of m and p variants) and least number of fragments assuming a p0 start. This selection and classification task was performed using the ParseProteogest.pl program, which also calculates N-terminal tryptic fragment weights, including possible modifications such as methionine formylation and fragment acetylation. The output of the ParseProteogest.pl program is available as Supplementary Table S3. The data for M. tuberculosis protein GlpX (Rv1099c) were added to the list, because although the protein does not meet all our selection criteria, we had strong previous evidence of incorrect translational start prediction for it (Movahedzadeh et al., 2004
). From these proteins, the 15 with the best available spectral data were selected for further investigation (Table 1
). All these proteins were tested for the presence of a signal peptide using the online SignalP 3.0 resource available at http://www.cbs.dtu.dk/services/SignalP/ (Bendtsen et al., 2004
), and no signal peptide was detected in any of them.
|
For the digestion of spots, the gel pieces were destained in 500 µl destaining buffer [200 mM NH4HCO3, 50 % acetonitrile (ACN)], equilibrated in 500 µl digestion buffer (50 mM NH4HCO3, 5 % ACN). The supernatant was removed and the gel piece was dried within 30 min in a Microconcentrator 5301 (Eppendorf) at 30 °C, and then digested in 25 µl digestion buffer containing 0.1 µg trypsin (sequencing-grade modified trypsin, Promega). After digestion overnight at 37 °C, the reaction tube was centrifuged and the supernatant (S1) transferred into another tube. To the gel piece, 25 µl 60 % ACN, 0.3 % trifluoroacetic acid (TFA) was added to stop the trypsin reaction and to shrink and wash the gel piece. After 10 min the supernatant (S2) was added to S1. The gel piece was washed and shrunk with 25 µl 100 % ACN. This supernatant was then added to S1+S2 and dried in the Microconcentrator 5301 at 65 °C. For the MS analysis, the dried peptides were dissolved in 1 µl 33 % ACN and 0.1 % TFA, and 0.25 µl of this solution was mixed with 0.5 µl alpha-cyano-4-hydroxycinnamic acid (CHCA) solubilized in 50 % ACN, 0.3 % TFA on Parafilm, and the resulting mixture was transferred to the template of the mass spectrometer and analysed.
The peptide masses were obtained using the following parameters: reflectron mode, 20 kV accelerating voltage, a low mass gate of 500 Da, and a mass range between 500 and 4000 Da. MS/MS spectra were obtained with and without collision gas. The database searches were performed with Mascot (http://www.matrixscience.com). Search parameters were: 30 p.p.m. peptide mass tolerance for peptide mass fingerprints, and 0.3 Da for MS/MS spectra. MS/MS data were compared with the theoretical sequence data of the protein under investigation. Search criteria were: one and two missed cleavages allowed and possible oxidation of methionine; N-terminal acetylation of the protein; and pyroglutamic acid formation from N-terminal glutamine, propionamide and sodium adducts. MS/MS spectra were manually evaluated.
The MS/MS peaklists can be found within the 2D-PAGE database (http://www.mpiib-berlin.mpg.de/2D-PAGE/), linked to the relevant protein.
| RESULTS |
|---|
|
|
|---|
Our predicted reassignment of the start of GlpX (Rv1099c) was based mainly on extensive comparative genomic analysis (Movahedzadeh et al., 2004
). We observed that amino acid homology extended beyond the annotated start, and we proposed that translation actually started 34 residues earlier (Fig. 1
). Crucially, some experimental N-terminal sequence data had been reported for the Corynebacterium glutamicum GlpX orthologue, Fbp (Rittmann et al., 2003
), which allowed us to be more confident about our new prediction for the M. tuberculosis GlpX. The extension of GlpX led to a conflict with the adjacent and divergently transcribed gene Rv1100, and we proposed an alternative start for this as well. This apparent misannotation therefore completely altered the predicted shared promoter region for two genes.
As GlpX had been identified proteomically in M. tuberculosis extracts (http://web.mpiib-berlin.mpg.de/cgi-bin/pdbs/2d-page/extern/index.cgi), we reasoned that existing MS data might be able to confirm or refute our reannotation. We therefore predicted peptide masses that would distinguish the predictions (Fig. 1
, Table 1
). Analysis of the existing data was inconclusive, but generation of new MS data confirmed (a) the presence of a peptide (NLAMELVR) that crossed the originally predicted start site and (b) the presence of a peptide that corresponded to our new prediction (Fig. 2
). No peak corresponding to the original start prediction was seen.
|
Developing a higher-throughput approach
The success with GlpX showed the potential of peptide mass mapping to experimentally verify protein starts on a large scale. We therefore calculated the masses of start peptides for all predicted proteins in the M. tuberculosis genome. We called the currently annotated start codon p0, and identification of tryptic peptides starting at this codon would support the current annotation.
We then identified potential alternative start codons upstream (m1, m2, etc. for minus) and downstream (p1, p2, etc. for plus) (see Fig. 1
); for each of these we recalculated the predicted masses of start peptides. Use of an upstream start codon would result in a longer protein, with the p0 codon-encoded residue contained within another peptide, the identification of which would be evidence against p0 (as occurred with GlpX). We therefore also calculated the masses of such trans-p0 peptides. The algorithm used is described in Methods, and shown in Fig. 1
.
To test our approach more rigorously, a shortlist of 15 proteins was drawn up. This was based on (a) the presence of existing MS data, (b) no obvious signal sequence being present, (c) the mass of the p0 peptide falling within the range of ideal sensitivity for MS, and (d) the peptide being preferably cleaved by trypsin at arginine rather than lysine residues (see Methods).
Table 1
shows the data that we obtained. We initially analysed existing data, which were adequate in some instances but not in others, requiring repeat MALDI-MS analyses. Of the 15 proteins, we had difficulties with two. In one case (Rv1017c; PrsA), the predicted N-terminal peptide coincided with a common contamination peak. Another (Rv2557) came from a spot containing more than one protein and was therefore unusable for the present investigation. The other 13 proteins were successfully analysed, and are briefly presented below.
Peptide start analyses
Confirmed starts.
We confirmed the starts for 11 proteins (see Table 1
). In most cases, a peptide lacking fMet was detected, although in one case (Rv0738), deformylated methionine was present. In two cases (ArgD and PrcA), both N-terminally acetylated peptides were found in addition to the non-acetylated form.
Revised starts
GlpX (Rv1099c): m2.
As described above, our hypothesis that GlpX is 34 aa residues larger than originally predicted was confirmed (Fig. 1
). The fMet is cleaved and, in addition to an unmodified protein species, there is also a less-abundant spot with an acetylated N-terminal peptide (Fig. 2
).
RibH (Rv1416): m1.
The ribH gene lies downstream of ribA2 in what appears to be an operon (Fig. 3
). There is a 14 bp gap between the stop codon of ribA2 and the predicted start valine of ribH. No peak corresponding to this start was seen; instead, there was a clear peak corresponding to the peptide crossing the originally predicted start. The m1 start lies 2 aa upstream of the peptide identified, and as the second residue is a lysine, should be cleaved in the procedure used. We identified both the predicted cleaved peptide (GGA...) and a peptide in which the trypsin had not cleaved (MKGGA...). Where there are alternative start codons (GTG, TTG), fMet is inserted rather than the usual amino acid that these would code for (Laursen et al., 2005
). The finding therefore of a methionine encoded by a GTG, which would normally code for valine, is in itself proof that this is the translational start. This sequence was confirmed by MS/MS. The peptide starting with GGA was confirmed by eight y ions (1, 2, 3, 5, 6, 9, 10, 12) and three b ions (4, 5, 8) and the peptide starting with MKGGA by eight y ions (1, 3, 4, 5, 6, 9, 10, 12). The expected ions y5 and y10, both with a cleavage after D and a mass loss of 64, further confirmed the sequence. This new ribH start would overlap with the ribA2 stop codon (GTGA; Fig. 3
) at position 4. Such an arrangement is known to be the most common for prokaryotic genes thought to be part of the same operon (Salgado et al., 2000
), as is the case with ribA2 and ribH.
|
| DISCUSSION |
|---|
|
|
|---|
Our approach is simple in concept, and has the great advantage that most of the data are routinely generated during standard proteomic procedures. TSS identification has previously been a relatively neglected area; in the past, N-terminal sequencing was the main method for protein identification, but this has been largely supplanted by more sensitive MALDI-MS technology. There have been large-scale N-terminal sequencing projects; for example, Link et al. (1997)
studied the N termini of 295 E. coli proteins. Of these, 72 failed because of signal peptide cleavage or N-terminal blockage. Of 223 N-terminal starts characterized, 10 (4.5 %) were reassigned. However, they used Edman degradation, and the need to use different chemistries makes this less practical as a routine method. We suggest that application of the method described here in proteomics laboratories would dramatically increase the number of experimentally determined protein TSSs, with minimal extra resources. There is also the possibility of retrospective examination using pre-existing MS datasets, although we preferred to carry out confirmatory analyses. This emphasizes the importance of storing MS data in proteome databases.
It would be possible to identify some of the predicted N-terminal peptides using a standard Mascot search, and the fact that this does not happen, or that the data do not feed into the wider community, may reflect different priorities (protein identification being the normal aim). However, our method goes further, and evidence for a start peptide is supported (a) by the absence of alternatives and (b) by the search for peptides spanning potential starts. We accept that as with all peptide mass mapping, while the identification of a single peak with a particular molecular mass provides supporting evidence for the presence of a particular peptide, additional MS/MS analysis to confirm peptide identity is required to be confident about the assignment.
More recent large-scale proteome projects have used proteomics to help genome annotation (Jaffe et al., 2004
; Lipton et al., 2002
). In a reanalysis of the Mycoplasma pneumoniae genome, over 81 % of the predicted ORFs were identified proteomically, and 16 new ORFs were detected. The algorithm used resulted in the extension of the predicted N termini of 19 proteins, on the basis that peptides were identified that extended beyond the annotated start. The translational start was then reassigned by looking upstream for a start codon. However, these groups did not actually identify start peptides.
We looked for proteins starting at ATG, GTG or TTG codons. ATG is the classical initiation codon, but in GC-rich genomes such as that of M. tuberculosis, GTG is extremely common. Our method conclusively confirms the starts in which a methionine residue is seen with a GTG or TTG codon, as this is the only scenario in which such a substitution would take place. Thus, methionine residues were detected with both RibH and Rv0738, even though they start with GTG (valine) codons.
There are limitations to our method. TSSs of proteins that are cleaved (for example removal of signal sequences) will not be resolved. However, it would be possible to adapt our method to identify such cleaved N termini, as signal sequence cleavage sites can be predicted. We deliberately chose a subset of proteins that would maximize the chance of success, suggesting that our success rate (84 %) would be less likely with a random set of proteins. One criterion we used was to choose proteins for which the predicted N-terminal peptide was readily resolvable using our experimental conditions and mass spectrometer. If a protein had an N-terminal peptide that was not suitable, then alternative proteases could be used, or an orthologue in another species might be more amenable. Also, the presence or absence of a peptide that crosses the predicted start might be detectable, even if the start peptide itself was not. We encountered another problem with Rv2908c, in which a predicted start (p2) corresponded to a tryptic peptide (because it lay next to an arginine or lysine residue); in such cases it would be necessary to include additional evidence before confidently assigning this as a start. Finally, our approach is directly applicable to all prokaryote genomes, but will need adapting for eukaryotes, in which exons need to be considered.
The applied method described here can use data from 2-DE/MS and LC/MS identifications. 2-DE/MS data have the advantage of high sequence coverage and therefore have a better chance to have access to the N terminus. The mass analysis of a 2-DE spot with peptide mass fingerprinting investigates only one protein, because 2-DE separates proteins. Therefore, the N-terminal peptide can be searched in a relatively small dataset of the proteolytically obtained peptides of this single protein. There is also access to MS and MS/MS data in already existing proteome databases (e.g. http://www.mpiib-berlin.mpg.de/2D-PAGE/). LC/MS methods have an advantage for sensitive identification of proteins, because only two or three peptides are sufficient for identification. This advantage turns into a limitation for the problem of the identification of posttranslational modifications and N termini (Schmidt et al., 2004
), because the probability that the N terminus or the modification is within the two or three identified peptides is very low. Nevertheless, our strategy for N-terminus confirmation or correction has the potential to be applied to the LC/MS approach by changing the experimental design to focus on the N-terminal peptide, as for example by the Cofradic approach (Gevaert et al., 2003
).
Both of the misannotated proteins that we identified were longer than the predicted proteins. Although our sample was small, we might expect this to be the case with a carefully annotated genome, as alignments to other homologues will tend to give a conservative indication of length. With genomes for which gene start assignment uses other assumptions, the balance of m and p reassignments might change.
One benefit of proteomic analyses is the ability to identify posttranslational modifications. A common modification is removal of fMet, and we observed this for 11 of the 13 proteins. In the two cases in which there was no cleavage, the residue had been deformylated. In some cases, methionine residues (both N-terminal and internal) had been oxidized (Schmidt et al., 2006
); we cannot tell if this occurred intracellularly or during the experimental procedure. Most interestingly, N-terminal peptides from three proteins (GlpX, ArgD and PrcA) were at least partially acetylated. MS/MS sequencing showed the acetyl group to be on the N-terminal residue; although N-terminal acetylation is the norm in eukaryotic proteins, it is reported to be rare in prokaryotes (Polevoda & Sherman, 2002
). It is amino acid dependent (Persson et al., 1985
), and our observations of acetylated serine or threonine residues are in line with previous work. N-terminal acetylation may not always alter function, but in some cases it does; thus, the important M. tuberculosis antigen ESAT-6, which normally interacts with the protein CFP-10, fails to do so when acetylated (Okkels et al., 2004
). We show here that three out of 13 proteins (23 %) have an acetylated form. This raises the question of whether the assumption that few prokaryotic proteins are acetylated is true for all bacteria. Much of the existing dogma has come from the analysis of E. coli proteins, of which in a survey of 810 proteins (Polevoda & Sherman, 2003
), only three were acetylated, all three being ribosomal subunits (S5, S18 and L12); our data suggest that the mycobacteria, at least, may be very different. A larger survey of proteins from these and other bacteria is needed to confirm this.
In this paper we suggest a simple strategy for the experimental identification of TSSs. The correct identification of these sites affects not only the length and possibly the function of the encoded proteins, but also key regulatory factors, such as the presence or absence of upstream regulatory sequences and the identification of genes likely to be co-expressed as part of an operon. We believe that the proposed nomenclature used herein for describing alternative start sites could be a useful standard for TSS identification in genome annotation. Since the data are often available and the strategy is therefore not particularly onerous in time or materials, we would recommend that researchers working on a particular gene attempt to mine proteomics databases for TSS identification, and that researchers performing large-scale proteomics investigations consider making TSS identification part of their protein-characterization procedure. The experimental identification of the TSS of a protein should be of use in annotating orthologues in other species, so that work in different laboratories would support a general effort to improve ORF annotation.
| ACKNOWLEDGEMENTS |
|---|
Edited by: J. Parkhill
| REFERENCES |
|---|
|
|
|---|
Brent, M. R. (2005). Genome annotation past, present, and future: how to define an ORF at each locus. Genome Res 15, 17771786.
Cagney, G., Amiri, S., Premawaradena, T., Lindo, M. & Emili, A. (2003). In silico proteome analysis to facilitate proteomics experiments using mass spectrometry. Proteome Sci 1, 5.[CrossRef][Medline]
Camus, J. C., Pryor, M. J., Medigue, C. & Cole, S. T. (2002). Re-annotation of the genome sequence of Mycobacterium tuberculosis H37Rv. Microbiology 148, 29672973.
Cole, S. T., Brosch, R., Parkhill, J., Garnier, T., Churcher, C., Harris, D., Gordon, S. V., Eiglmeier, K., Gas, S. & other authors (1998). Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature 393, 537544.[CrossRef][Medline]
Doherty, N. S., Littman, B. H., Reilly, K., Swindell, A. C., Buss, J. M. & Anderson, N. L. (1998). Analysis of changes in acute-phase plasma proteins in an acute inflammatory response and in rheumatoid arthritis using two-dimensional gel electrophoresis. Electrophoresis 19, 355363.[CrossRef][Medline]
Edman, P. J. (1950). Method for determination of the amino acid sequence in peptides. Acta Chem Scand 4, 283293.
Edwards, M. T., Rison, S. C., Stoker, N. G. & Wernisch, L. (2005). A universally applicable method of operon map prediction on minimally annotated genomes using conserved genomic context. Nucleic Acids Res 33, 32533262.
Gevaert, K., Goethals, M., Martens, L., Van Damme, J., Staes, A., Thomas, G. R. & Vandekerckhove, J. (2003). Exploring proteomes and analyzing protein processing by mass spectrometric identification of sorted N-terminal peptides. Nat Biotechnol 21, 566569.[CrossRef][Medline]
Jaffe, J. D., Berg, H. C. & Church, G. M. (2004). Proteogenomic mapping as a complementary method to perform genome annotation. Proteomics 4, 5977.[CrossRef][Medline]
Jungblut, P. R., Schaible, U. E., Mollenkopf, H. J., Zimny-Arndt, U., Raupach, B., Mattow, J., Halada, P., Lamer, S., Hagens, K. & Kaufmann, S. H. (1999). Comparative proteome analysis of Mycobacterium tuberculosis and Mycobacterium bovis BCG strains: towards functional genomics of microbial pathogens. Mol Microbiol 33, 11031117.[CrossRef][Medline]
Jungblut, P. R., Bumann, D., Haas, G., Zimny-Arndt, U., Holland, P., Lamer, S., Siejak, F., Aebischer, A. & Meyer, T. F. (2000). Comparative proteome analysis of Helicobacter pylori. Mol Microbiol 36, 710725.[CrossRef][Medline]
Jungblut, P. R., Muller, E. C., Mattow, J. & Kaufmann, S. H. (2001). Proteomics reveals open reading frames in Mycobacterium tuberculosis H37Rv not predicted by genomics. Infect Immun 69, 59055907.
Krause, E., Wenschuh, H. & Jungblut, P. R. (1999). The dominance of arginine-containing peptides in MALDI-derived tryptic mass fingerprints of proteins. Anal Chem 71, 41604165.[Medline]
Lamer, S. & Jungblut, P. R. (2001). Matrix-assisted laser desorption-ionization mass spectrometry peptide mass fingerprinting for proteome analysis: identification efficiency after on-blot or in-gel digestion with and without desalting procedures. J Chromatogr B Biomed Sci Appl 752, 311322.[CrossRef][Medline]
Laursen, B. S., Sorensen, H. P., Mortensen, K. K. & Sperling-Petersen, H. U. (2005). Initiation of protein synthesis in bacteria. Microbiol Mol Biol Rev 69, 101123.
Link, A. J., Robison, K. & Church, G. M. (1997). Comparing the predicted and observed properties of proteins encoded in the genome of Escherichia coli K-12. Electrophoresis 18, 12591313.[CrossRef][Medline]
Lipton, M. S., Pasa-Tolic, L., Anderson, G. A. & other authors (2002). Global analysis of the Deinococcus radiodurans proteome by using accurate mass tags. Proc Natl Acad Sci U S A 99, 1104911054.
Movahedzadeh, F., Rison, S. C., Wheeler, P. R., Kendall, S. L., Larson, T. J. & Stoker, N. G. (2004). The Mycobacterium tuberculosis Rv1099c gene encodes a GlpX-like class II fructose 1,6-bisphosphatase. Microbiology 150, 34993505.
Okkels, L. M., Muller, E. C., Schmid, M., Rosenkrands, I., Kaufmann, S. H., Andersen, P. & Jungblut, P. R. (2004). CFP10 discriminates between nonacetylated and acetylated ESAT-6 of Mycobacterium tuberculosis by differential interaction. Proteomics 4, 29542960.[CrossRef][Medline]
Persson, B., Flinta, C., von Heijne, G. & Jornvall, H. (1985). Structures of N-terminally acetylated proteins. Eur J Biochem 152, 523527.[Medline]
Polevoda, B. & Sherman, F. (2002). The diversity of acetylated proteins. Genome Biol 3, reviews0006 (http://genomebiology.com/2002/3/5/reviews/0006).
Polevoda, B. & Sherman, F. (2003). N-terminal acetyltransferases and sequence requirements for N-terminal acetylation of eukaryotic proteins. J Mol Biol 325, 595622.[CrossRef][Medline]
Rittmann, D., Schaffer, S., Wendisch, V. F. & Sahm, H. (2003). Fructose-1,6-bisphosphatase from Corynebacterium glutamicum: expression and deletion of the fbp gene and biochemical characterization of the enzyme. Arch Microbiol 180, 285292.[CrossRef][Medline]
Salgado, H., Moreno-Hagelsieb, G., Smith, T. F. & Collado-Vides, J. (2000). Operons in Escherichia coli: genomic analyses and predictions. Proc Natl Acad Sci U S A 97, 66526657.
Schmidt, F., Donahoe, S., Hagens, K., Mattow, J., Schaible, U. E., Kaufmann, S. H., Aebersold, R. & Jungblut, P. R. (2004). Complementary analysis of the Mycobacterium tuberculosis proteome by two-dimensional electrophoresis and isotope-coded affinity tag technology. Mol Cell Proteomics 3, 2442.
Schmidt, F., Krah, A., Schmid, M., Jungblut, P. R. & Thiede, B. (2006). Distinctive mass losses of tryptic peptides generated by matrix-assisted laser desorption/ionization time-of-flight/time-of-flight. Rapid Commun Mass Spectrom 20, 933936.[CrossRef][Medline]
Received 30 August 2006;
revised 14 November 2006;
accepted 14 November 2006.
This article has been cited by other articles:
![]() |
C. Ansong, S. O. Purvine, J. N. Adkins, M. S. Lipton, and R. D. Smith Proteogenomics: needs and roles to be filled by proteomics in genome annotation Brief Funct Genomic Proteomic, March 10, 2008; (2008) eln010v1. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |