Research Article

Genomovars 11 to 18 of Pseudomonas stutzeri, identified among isolates from soil and marine sediment

,, Jorge Lalucat2 and Wilfried Wackernagel1

1 Genetik, Institut für Biologie und Umweltwissenschaften, Universität Oldenburg, Germany
2 Microbiologia, Departament de Biologia, Institut Mediterrani d'Estudis Avancats (CSICUIB), Palma de Mallorca, Spain

Correspondence
Johannes Sikorski
sikorski{at}research.haifa.ac.il

International Journal of Systematic and Evolutionary Microbiology 2005; 55(5):1767 · https://doi.org/10.1099/ijs.0.63535-0

View at publisher PubMed

Abstract

Amongst 440 strains of Pseudomonas stutzeri isolated from soil and marine sediment for a population genetic study, eight strains were each presumed to represent a novel genomic group and were compared with each other and to reference strains of P. stutzeri genomovars 1 to 10 and other Pseudomonas species by DNADNA hybridization, 16S rRNA and internally transcribed 16S23S rRNA spacer region (ITS1) sequences and basic physiological properties defining the species. While 16S rRNA and ITS1 gene sequences positioned the eight strains within the phylogenetic branch of P. stutzeri, the DNADNA hybridizations with reference strains of the 10 described genomovars and among the novel strains were generally below 70 %, which is the threshold for species and genomovar differentiation. Since the physiological properties studied in the eight strains fitted the profile of P. stutzeri, eight new genomovars of P. stutzeri, numbered 11 to 18, are proposed, with strains 28a50, 28a39, 28a22, 28a3, 4C29, 24a13, 24a75 and MT-1 being the reference strains. The highly transformable reference strain 28a3 of genomovar 14 had a localized 16S rRNA gene sequence tag characteristic of genomovar strains 2 and 3, suggesting a possible horizontal gene transfer event involving part of the 16S rRNA gene.
Abbreviations: RAPD, random amplified polymorphic DNA

Published online ahead of print on 22 April 2005 as DOI 10.1099/ijs.0.63535-0.

The GenBank/EMBL/DDBJ accession numbers for the ITS1 rRNA gene sequences of genomovars 11 to 18 of Pseudomonas stutzeri are AY850023AY850031.

Two phylogenetic trees based on 16S rRNA and ITS1 gene sequences and tables detailing DNADNA relatedness values and physiological properties of the genomovars are available as supplementary material in IJSEM Online.

J. S. is grateful for a research stay at the Departament de Biologia of the Universitat de les Illes Balears (Spain). This work was supported by grants from the Federal Ministry of Research and Technology (Germany), the Fonds der Chemischen Industrie (Germany) and the CICYT (Spain) project BOS2001-0303. The help of Maria Valens and Christina Ramon in DNADNA hybridization experiments is highly acknowledged. Also, the help of Ramon Rosselló-Mora in DNADNA hybridization, phylogenetic reconstructions and many helpful discussions is greatly acknowledged.

Footnotes

Present address: Institute of Evolution, University of Haifa, Mount Carmel, Haifa 31905, Israel.



The species Pseudomonas stutzeri is a non-fluorescent member of the genus Pseudomonas (γ-Proteobacteria) displaying high genetic (Rius et al., 2001; Cladera et al., 2004) and physiological diversity (Rosselló et al., 1991). Strains of P. stutzeri have been isolated from a variety of environmental and clinical habitats (Sikorski et al., 2002a and references therein). Some strains received attention as model organisms because of their specific metabolic properties (Musarrat & Hashsham, 2003; Obradors & Aguilar, 1991; Rosselló-Mora et al., 1994; Zumft, 1997) and their ability for natural genetic transformation (Berndt et al., 2003; Meier & Wackernagel, 2003; Sikorski et al., 1998, 2002b).

Taxonomically, P. stutzeri strains have been grouped into 10 genomovars by DNADNA hybridization (García-Valdés et al., 2003; Rosselló et al., 1991; Rosselló-Mora et al., 1996; Sepúlveda-Torres et al., 2001; Ursing et al., 1995). In a recent population-genetic study, approximately 440 strains from soil and marine environments were studied by random amplified polymorphic DNA-PCR (RAPD-PCR) and the 16S rRNA gene sequence was determined (>1450 bp) (Sikorski et al., 2002a) for 34 of the strains (several being representatives of the main RAPD groups). The 16S rRNA gene sequences suggested that seven representative strains (28a50, 28a39, 28a22 and 28a3, from soil close to Tel Aviv airport, Israel; 4C29, from marine sediment on the shore of the North Sea coast, Germany; 24a13 and 24a75, from a soil contaminated with mineral oil, Germany) and strain MT-1 (from Mariana Trench, Japan; Tamegai et al., 1997) were members of new genomovars. Based on their 16S rRNA gene sequence dissimilarity values, these eight strains were as different from each other and from reference strains of the established genomovars 1 to 10 as the genomovar reference strains differed from each other (Table 1). Moreover, the lowest dissimilarity value of any of these strains to each other or to a reference strain of an established genomovar (0·41 %; Table 1) was larger than the maximum dissimilarity value observed among nine strains within three established genomovars (0·31 %; Table 1). In the following, we present further data to support the identification of eight new genomovars including the results of DNADNA hybridization studies, sequence analyses of 16S23S rRNA internally transcribed spacer regions (ITS1) and comparisons of physiological properties. Additionally, two strains (28a18 and 28a69) which were suggested by RAPD-PCR to belong to the new genomovars represented by strains 28a3 (genomovar 14) and 28a22 (genomovar 13), respectively (Sikorski et al., 2002a), were included in part of the studies.


Table 1. 16S rRNA gene sequence dissimilarities Strains: 1, Genomovar 1 reference strains CCUG 11256T, ATCC 17589, ATCC 17593; 2, genomovar 2 reference strains ATCC 17591, ATCC 14405, ATCC 17587; 3, genomovar 3 reference strains DSM 50227, AN10, AN11; 4, genomovar type strain CCUG 11256T (gv. 1) and reference strains ATCC 17591 (gv. 2), DSM 50227 (gv. 3), DSM 6084 (gv. 4), DSM 6082 (gv. 5), DSM 50238 (gv. 7), JM 300 (gv. 8), KC (gv. 9), CLN 100 (gv. 10).


Genomic DNA from the eight new genomovar reference strains and from strains 28a18 and 28a69 was isolated using the Genomic DNA kit from QIAGEN and from the reference strains of genomovars 110 by the method of Marmur (1961). DNADNA hybridizations were performed using a modification of the hydroxyapatite method as described by Ziemke et al. (1998). Reference DNAs were double-labelled with DIG-11-dUTP and biotin-16-dUTP, using a nick-translation kit (Boehringer Mannheim). Sequences of the 16S rRNA and ITS1 genes were determined as described by Guasp et al. (2000) and Sikorski et al. (2002a). Phylogenetic analysis was performed using the ARB package for 16S rRNA gene sequences (Ludwig et al., 2004) and MEGA version 2.1 (Kumar et al., 2001) for ITS1 sequences. Physiological tests were done according to Smibert & Krieg (1994).

The high DNADNA relatedness values of strain 28a18 with 28a3 (91 %) and 28a69 with 28a22 (93 %) indicate membership of the respective genomovar, as suggested by RAPD-PCR (Sikorski et al., 2002a). In contrast, the DNADNA relatedness values of the eight representative strains of the new genomovars to each other and to the reference strains of the genomovars 1 to 10 were at or below the threshold value of 70 % for species delineation (Table 2) except for two pairs (77 %, MT-1/28a22; 82 %, 28a3/24a13; Supplementary Table S1 available in IJSEM Online). However, the 16S rRNA and ITS1 gene sequence analyses group all of them into the same phylogenetic branch as the previously described strains of P. stutzeri (Supplementary Figs S1 and S2 available in IJSEM Online), indicating their phylogenetic affiliation with the P. stutzeri group. Additionally, physiological traits characteristic of the species P. stutzeri were present in the eight strains (see Supplementary Table S2 in IJSEM Online). Motility, denitrification and starch hydrolysis are characteristic of most P. stutzeri strains, but exceptions have been described (strain JM300 is amylase-negative, strain ZoBell is described as non-motile). New genomic species, as indicated here by DNADNA hybridization for the eight strains, should not be classified as novel species unless differentiating phenotypes are found (Rosselló-Mora & Amann, 2001; Stackebrandt et al., 2002; Ursing et al., 1995), which is not the case here. Thus, we propose eight new genomovars, 11 to 18, of P. stutzeri, with strains 28a50, 28a39, 28a22, 28a3, 4C29, 24a13, 24a75 and strain MT-1, respectively, as the reference strains for each genomovar. The strains have been deposited as CCUG 5053850545 (=DSM 1708217089).


Table 2. DNADNA relatedness values (%) The strains used for comparisons were the same as in Table 1.


It is remarkable that strain 28a18 (Sikorski et al., 2002a) was found to be highly similar to the genomovar 14 representative 28a3, with a DNADNA relatedness value of 91 % (Supplementary Table S1 in IJSEM Online), yet this strain had a position in the 16S rRNA tree different from that of 28a3 (Supplementary Fig. S1 available in IJSEM Online). The two strains were isolated from the same soil sample and are nearly identical in their ITS1 sequences (Supplementary Fig. S2 available in IJSEM Online), their RAPD patterns (Sikorski et al., 2002a) and their partial rpoB sequences [as determined by restriction enzyme profiling of 1·5 kb PCR products (nucleotides 532 to 2034 of the Escherichia coli rpoB sequence); J. Sikorski and W. Wackernagel, unpublished]. The divergence of their 16S rRNA gene sequences results from five nucleotide changes within a stretch of 15 nucleotides (E. coli positions 74 to 92; E. coli has an insert of four nucleotides compared with all P. stutzeri strains), which makes the 16S rRNA gene sequence of strain 28a3 in this region identical to that of genomovars 2 and 3 strains. This rare specific local sequence identity between a single member of one genomovar and the members of two far-distant genomovars may be explained by a horizontal gene transfer event involving part of the 16S rRNA gene. Strain 28a3, as the putative recipient, was shown to be highly transformable, which is generally not the case for strains of the RAPD group of which 28a3 is the representative (Sikorski et al., 2002b).

References

Bennasar, A., Rosselló-Mora, R., Lalucat, J. & Moore, E. R. (1996). 16S rRNA gene sequence analysis relative to genomovars of Pseudomonas stutzeri and proposal of Pseudomonas balearica sp. nov. Int J Syst Bacteriol 46, 200205.[Abstract/Free Full Text]

Berndt, C., Meier, P. & Wackernagel, W. (2003). DNA restriction is a barrier to natural transformation in Pseudomonas stutzeri JM300. Microbiology 149, 895901.[Abstract/Free Full Text]

Cladera, A., Bennasar, A., Barceló, M., Lalucat, J. & García-Valdés, E. (2004). Comparative genetic diversity of Pseudomonas stutzeri genomovars, clonal structure, and phylogeny of the species. J Bacteriol 186, 52395248.[Abstract/Free Full Text]

García-Valdés, E., Castillo, M. M., Bennasar, A., Guasp, C., Cladera, A. M., Bosch, R., Engesser, K. H. & Lalucat, J. (2003). Polyphasic characterization of Pseudomonas stutzeri CLN100 which simultaneously degrades chloro- and methylaromatics: a new genomovar within the species. Syst Appl Microbiol 26, 390403.[CrossRef][Medline]

Guasp, C., Moore, E. R., Lalucat, J. & Bennasar, A. (2000). Utility of internally transcribed 16S-23S rDNA spacer regions for the definition of Pseudomonas stutzeri genomovars and other Pseudomonas species. Int J Syst Evol Microbiol 50, 16291639.[Abstract]

Kumar, S., Tamura, K., Jakobsen, I. B. & Nei, M. (2001). MEGA2: molecular evolutionary genetics analysis software. Bioinformatics 17, 12441245.[Abstract/Free Full Text]

Ludwig, W., Strunk, O., Westram, R. & 29 other authors (2004). ARB: a software environment for sequence data. Nucleic Acids Res 32, 13631371.[Abstract/Free Full Text]

Marmur, J. (1961). A procedure for the isolation of deoxyribonucleic acid from microorganisms. J Mol Biol 3, 208218.

Meier, P. & Wackernagel, W. (2003). Mechanisms of homology-facilitated illegitimate recombination for foreign DNA acquisition in transformable Pseudomonas stutzeri. Mol Microbiol 48, 11071118.[CrossRef][Medline]

Musarrat, J. & Hashsham, S. A. (2003). Customized cDNA microarray for expression profiling of environmentally important genes of Pseudomonas stutzeri strain KC. Teratog Carcinog Mutagen (Suppl. 1), 283294.

Obradors, N. & Aguilar, J. (1991). Efficient biodegradation of high-molecular-weight polyethylene glycols by pure cultures of Pseudomonas stutzeri. Appl Environ Microbiol 57, 23832388.[Abstract/Free Full Text]

Rius, N., Fusté, M. C., Guasp, C., Lalucat, J. & Loren, J. G. (2001). Clonal population structure of Pseudomonas stutzeri, a species with exceptional genetic diversity. J Bacteriol 183, 736744.[Abstract/Free Full Text]

Rosselló, R. A., García-Valdés, E., Lalucat, J. & Ursing, J. (1991). Genotypic and phenotypic diversity of Pseudomonas stutzeri. Syst Appl Microbiol 14, 150157.

Rosselló-Mora, R. & Amann, R. (2001). The species concept for prokaryotes. FEMS Microbiol Rev 25, 3967.[Medline]

Rosselló-Mora, R. A., Lalucat, J. & García-Valdés, E. (1994). Comparative biochemical and genetic analysis of naphthalene degradation among Pseudomonas stutzeri strains. Appl Environ Microbiol 60, 966972.[Abstract/Free Full Text]

Rosselló-Mora, R. A., Lalucat, J. & Moore, E. R. B. (1996). Strain JM300 represents a new genomovar within Pseudomonas stutzeri. Syst Appl Microbiol 19, 596599.

Sepúlveda-Torres, L. C., Zhou, J., Guasp, C., Lalucat, J., Knaebel, D., Plank, J. L. & Criddle, C. S. (2001). Pseudomonas sp. strain KC represents a new genomovar within Pseudomonas stutzeri. Int J Syst Evol Microbiol 51, 20132019.[Abstract]

Sikorski, J., Graupner, S., Lorenz, M. G. & Wackernagel, W. (1998). Natural genetic transformation of Pseudomonas stutzeri in a non-sterile soil. Microbiology 144, 569576.[Abstract/Free Full Text]

Sikorski, J., Moehle, M. & Wackernagel, W. (2002a). Identification of complex composition, strong strain diversity and directional selection in local Pseudomonas stutzeri populations from marine sediment and soils. Environ Microbiol 4, 465476.[CrossRef][Medline]

Sikorski, J., Teschner, N. & Wackernagel, W. (2002b). Highly different levels of natural transformation are associated with genomic subgroups within a local population of Pseudomonas stutzeri from soil. Appl Environ Microbiol 68, 865873.[Abstract/Free Full Text]

Smibert, R. & Krieg, N. R. (1994). Phenotypic characterization. In Methods for General and Molecular Bacteriology, pp. 607654. Edited by P. Gerhardt, R. G. E. Murray, W. A. Wood & N. R. Krieg. Washington, DC: American Society for Microbiology.

Stackebrandt, E., Frederiksen, W., Garrity, G. M. & 10 other authors (2002). Report of the ad hoc committee for the re-evaluation of the species definition in bacteriology. Int J Syst Evol Microbiol 52, 10431047.[Abstract]

Tamegai, H., Li, L., Masui, N. & Kato, C. (1997). A denitrifying bacterium from the deep sea at 11,000-m depth. Extremophiles 1, 207211.[CrossRef][Medline]

Ursing, J. B., Rosselló-Mora, R. A., García-Valdés, E. & Lalucat, J. (1995). Taxonomic note: a pragmatic approach to the nomenclature of phenotypically similar genomic groups. Int J Syst Bacteriol 45, 604.[Abstract/Free Full Text]

Ziemke, F., Hofle, M. G., Lalucat, J. & Rosselló-Mora, R. (1998). Reclassification of Shewanella putrefaciens Owen's genomic group II as Shewanella baltica sp. nov. Int J Syst Bacteriol 48, 179186.[Abstract/Free Full Text]

Zumft, W. G. (1997). Cell biology and molecular basis of denitrification. Microbiol Rev 61, 533616.[Abstract/Free Full Text]