Research Article

Probable synonymy of the nitrogen-fixing genus Azotobacter and the genus Pseudomonas

International Journal of Systematic and Evolutionary Microbiology 2007; 57(12):2894 · https://doi.org/10.1099/ijs.0.64969-0

View at publisher PubMed

Abstract

The relationships of the genus Azotobacter, Azomonas macrocytogenes and the genus Pseudomonas were revealed by comparative analysis of partial 16S rRNA and atpD, carA and recA gene sequences and as concatenated nucleotide and peptide sequences. Sequence similarities of Azotobacter species and Azomonas macrocytogenes indicated that these may be considered to be synonyms at the molecular level. In addition, these species show an intimate relationship with species of Pseudomonas, especially P. aeruginosa (the type species of the genus). In terms of the current circumscription of the genus Pseudomonas, Azotobacter and Azomonas macrocytogenes should be considered for amalgamation with Pseudomonas. Azotobacter and Azomonas comprise nitrogen-fixing strains with large pleomorphic cells that form cysts, and peritrichous flagella insertion; characteristics not included in the current circumscription of Pseudomonas. The data are discussed in the light of whether lateral transfer of genes could be involved in the determination of significant morphological characteristics, thus leading to a problem that may be encountered more frequently: how to resolve classification of taxa based on conserved sequences with those based on their phenotype. More fundamentally, the results illuminate problems that will increasingly be encountered: by what criteria can taxa be delineated, what are the most appropriate methods for classification, and what are the proper assumptions of bacterial classification?
Abbreviations: ML, maximum-likelihood; MLSA, multilocus sequence analysis; NJ, neighbour-joining

The GenBank/EMBL/DDBJ accession numbers for the 16S rRNA, atpD, carA and recA gene sequences determined in this study are given in Fig. 1 and Supplementary Figs S1 and S2.

All authority references for names can be found at Euzéby (1997).

Figures showing the comparative analysis of individual atpD, carA and recA nucleotide and peptide sequences of Azomonas, Azotobacter and Pseudomonas species using the maximum-likelihood and neighbour-joining algorithms are available as supplementary material with the online version of this paper.

The genus Pseudomonas has been circumscribed as a heterogeneous taxon (Palleroni, 1984). Revisions of the genus based on comparative analyses of 16S rRNA gene sequence data have resulted in the transfer of many species previously included in Pseudomonas, usually to novel genera proposed in the Alpha-, Beta- and Gammaproteobacteria (Kersters et al., 1996). Pseudomonas now includes species in a clade in the Gammaproteobacteria. Relocation of taxa within Pseudomonas (sensu stricto) has been indicated by comparative analyses of 16S rRNA gene sequence data (Moore et al., 1996; Anzai et al., 2000) and on combined data of gyrase B and the rho factor (Yamamoto et al., 2000). Anzai et al. (2000) noted a close relationship between Pseudomonas, Azotobacter and Azomonas. Azotobacter and Azomonas are specifically circumscribed as genera containing species capable of fixing nitrogen in aerobic conditions (Kennedy & Rudnick, 2005; Kennedy et al., 2005). Pseudomonas is not recorded as having this capability (Palleroni, 2005) although this has been reported for a limited number of species, some of which were of doubtful identification (Line, 1997; Vermeiren et al., 1999; Kulakov et al., 2002).

A preliminary comparison of 16S rRNA gene sequences suggested a much closer relationship between Azotobacter beijerinckii and Pseudomonas species than was previously indicated by their inclusion in separate families of the order Pseudomonadales (Palleroni, 1984), now as separate genera in the family Pseudomonadaceae (Palleroni, 2005). These considerations led to a re-evaluation of the relationships using partial sequences of the catalytic subunit of the ATP synthase (atpD), the small subunit of carbamoyl phosphate synthase (carA), recombinase A (recA) and the 16S rRNA gene. Sequences from Azomonas macrocytogenes (the type species Azomonas agilis was not available to us) and Azotobacter species were compared with those representing Pseudomonas. Representative strains of Azospirillum species were included to confirm the classification of this genus also.

Strains used.
Strains of Azomonas, Azospirillum, Azotobacter and Pseudomonas species were obtained from the International Collection of Microorganisms from Plants (ICMP), Landcare Research, Auckland. Information on these strains is given at the ICMP website:

DNA extraction.
Strains were cultured on nutrient agar (Difco) incubated at 27 °C for 2–3 days. Genomic DNA was extracted directly from colony growth on nutrient agar using phenol and CTAB (hexadecyltrimethylammonium bromide) (Ausubel et al., 1987).

Amplification and sequencing.
The primers and conditions used for PCR amplification of atpD, carA, recA and 16S rRNA gene sequences as given by Young & Park (2007) were followed. PCR products were column-purified with a Roche High Pure PCR Product Purification kit. In cases where more than one DNA band occurred on the gel, the DNA fragment of the correct size was extracted using the PCR purification kit. Purified products were cycle sequenced with the appropriate primers using BigDye Terminator Ready Reaction Mix v3.1 (ABI) and sequences were obtained in both directions with an ABI 3100 Avant Genetic Analyzer. Sequences were assembled and edited with Sequencher 4.5 (Gene Codes Corp.).

Sequence analysis.
All 16S rRNA, atpD, carA and recA gene sequences for Azomonas, Azospirillum, Azotobacter and Pseudomonas species included in this study have been registered in GenBank. GenBank accession numbers for the 16S rRNA genes are shown in Fig. 1. GenBank accession numbers for individual gene sequences are shown in Supplementary Figs S1 and S2 in IJSEM Online.



(32K):

Fig. 1. Comparative analysis of 16S rRNA gene sequences of Azomonas (Azm.), Azospirillum (Azs.), Azotobacter (Azt.) and Pseudomonas species, using the NJ algorithm. Bar, number of substitutions per site. Accession numbers of sequences are held in GenBank. Nucleotide sequences marked with an asterix (*) are derived from type strains. Sequences AY324110*, EF100149 and EF100150* refer to Azospirillum brasilense (ATCC 29145T) and Azospirillum lipoferum (8672 and ATCC 29707T), respectively.

The outgroups for alignments of the 16S rRNA, atpD, carA and recA genes were the appropriate sequence for the type strain of Escherichia coli (ICMP 15663T) and Xylella fastidiosa clone 9a5c, obtained directly from GenBank.

Nucleotide sequences were aligned with CLUSTAL_X 1.83 with default parameters and both ends of each alignment were trimmed to the following sizes: 16S rRNA gene, 1386 positions; atpD, 793 positions; carA, 630 positions; recA, 598 positions.

Peptide sequences for atpD, carA and recA were determined using GeneDoc, aligned using CLUSTAL_X and edited manually, again using GeneDoc, to ensure that alignment gaps in protein-coding genes did not cause errors in the amino acid sequence. Calculated peptide sequence lengths for atpD, carA and recA were 263, 209 and 198 amino acids respectively.

Individual nucleotide and peptide sequences were compared using PAUP* version 4.0b10 with the neighbour-joining (NJ) algorithm. For the distance measure, an uncorrected p distance and Jukes–Cantor model were used.

The atpD, carA and recA sequences were concatenated and trimmed to give a total length of 2070 positions, which resulted in a concatenated peptide length of 670 amino acids.

Concatenated sequences were treated using Gblocks 9.1b () (Castresana, 2000), which eliminates poorly aligned positions and divergent regions of an alignment of DNA sequences; positions which may not be homologous or may have been saturated by multiple substitutions. Because the atpD sequences for Pseudomonas putida and Pseudomonas stutzeri were not available, a non-stringent option was used to allow inclusion of the sequences from other strains. Concatenated nucleotide and peptide sequences were also analysed using the NJ algorithm as described above and maximum-likelihood (ML) using PHYML Online (Guindon et al., 2005) with the default models for nucleotides (HYK85) and for peptides (Blosum62).

16S rRNA gene sequences
The comparative analysis of 16S rRNA gene sequences using the NJ algorithm produced a tree (Fig. 1) with sequences representing Azospirillum (Azospirillum brasilense and Azospirillum lipoferum) as outliers to a complex cluster representing Pseudomonas, Azomonas and Azotobacter. Sequences representing Pseudomonas formed two groups: group 1, comprising P. aeruginosa, P. flavescens, P. mendocina, P. putida, P. resinovorans and P. stutzeri; and group 2, comprising P. agarici, P. aureofaciens, P. aurantiaca, P. chlororaphis, P. cichorii, P. fluorescens, P. syringae and P. tolaasii. Sequences representing Azotobacter and Azomonas macrocytogenes were grouped together within rRNA group 1 of Pseudomonas, being closest to P. aeruginosa and P. resinovorans.

Individual atpD, carA and recA gene and peptide sequences
Nucleotide sequences of atpD, carA and recA were obtained for all strains except those of Azospirillum species. The online version of this paper contains supplementary figures showing results of analyses of individual nucleotide and peptide sequences based on ML (Supplementary Fig. S1) and NJ (Supplementary Fig. S2) algorithms.

atpD sequences.
Nucleotide sequences divided into two discrete groups. All sequences representing Azotobacter and Azomonas macrocytogenes were present in one group together with P. aeruginosa and P. resinovorans. The remaining sequences of Pseudomonas species formed a second group. Peptide sequences produced a tree in which sequences of Azotobacter and Azomonas macrocytogenes were outliers to Pseudomonas.

carA sequences.
Nucleotide sequences of Azomonas macrocytogenes were outliers to a complex AzomonasAzotobacterPseudomonas structure in which Azotobacter sequences formed a discrete group embedded among Pseudomonas sequences. Peptide sequences produced a tree in which Pseudomonas was divided into the two groups expressed by the 16S rRNA gene, with sequences representing Azotobacter and Azomonas macrocytogenes included occurring separately in rRNA group 1.

recA sequences.
Nucleotide sequences of Azomonas macrocytogenes were outliers to a complex AzotobacterPseudomonas structure in which Azotobacter species sequences formed a discrete group embedded within the Pseudomonas complex, adjacent to P. aeruginosa. The peptide tree has a different structure, with Azomonas macrocytogenes included with Azotobacter sequences and P. stutzeri as nearest neighbour. These are embedded in sequences identified with rRNA group 1.

There was considerable minor variability between tree structures of individual Pseudomonas atpD, recA and carA nucleotide and peptide sequences, whether analysed using NJ or ML, but in general terms, individual nucleotide and peptide trees were similar (see Supplementary Figs S1 and S2 in IJSEM Online). For Azomonas macrocytogenes, the atpD sequence was similar to those of Azotobacter, but carA and recA nucleotides were outliers to those of AzotobacterPseudomonas and the carA peptide was nearest neighbour to P. stutzeri and P. flavescens sequences. The recA peptide was grouped with all Azotobacter sequences.

Removal of Azomonas macrocytogenes nucleotide and peptide sequences from the analyses or removal of all atpD data resulted in trees with nitrogen-fixing species embedded in Pseudomonas, with P. aeruginosa as nearest neighbour (data not shown).

Concatenated sequences
Analyses of concatenated nucleotide and peptide data treated using Gblocks indicated that <2 % of the sequences might be poorly aligned or hypervariable regions. These regions were considered to be so small that all sequence data were included. Analyses of concatenated sequences are shown in Fig. 2. The concatenated nucleotide and peptide sequences analysed using the NJ algorithm both produced trees in which Azotobacter species and Azomonas macrocytogenes were outliers to Pseudomonas species. The ML algorithm applied to nucleotide and peptide sequences resulted in trees in which Azotobacter and Azomonas macrocytogenes lay within the Pseudomonas aeruginosa cluster. The tree comparing nucleotide sequences was not very well resolved, but Azotobacter and Azomonas macrocytogenes sequences adjoined P. aeruginosa and P. resinovorans. The peptide tree reproduced the aeruginosa and fluorescens clusters, with Azotobacter and Azomonas macrocytogenes sequences adjoining P. aeruginosa and P. resinovorans.



(38K):

Fig. 2. Comparative analysis of concatenated atpD, carA and recA nucleotide and peptide sequences of Azomonas (Azm.), Azotobacter (Azt.) and Pseudomonas species, using the NJ algorithm (a) and the ML algorithm (b). Accession numbers of strains held in the International Collection of Microorganisms from Plants are specified. Bars, number of substitutions per site.
Three 16S rRNA gene sequences from species representing the genus Azospirillum, Azospirillum brasilense and Azospirillum lipoferum, were outliers to all other sequences, including E. coli. This is consistent with the placement of Azospirillum in the family Rhodospirillaceae in the Alphaproteobacteria (Pfennig & Trüper, 2005), therefore unrelated to Pseudomonas and the other nitrogen-fixing genera under consideration, and the genus is not considered further in this study.

The analysis of 16S rRNA gene sequences of Pseudomonas produced the two groups corresponding to the aeruginosa cluster and fluorescens cluster reported previously (Hilario et al., 2004; Moore et al., 1996). [As well as comparing 16S rRNA gene sequences, Hilario et al. (2004) concatenated the rRNA gene together with three genes, without appreciating that the rRNA gene has regions that are continuously conserved, whereas degeneracy in codons means that genes are more variable and therefore these different forms of sequence make incommensurable contributions when they are concatenated.] The nucleotide and peptide sequence analyses reported here also gave general support for this subdivision. However, these studies have included relatively small representations of the approximately 110 species now included in the genus (Euzéby, 1997). Where larger numbers have been compared using different nucleotides for comparison, this differentiation is not so obvious and may be obliterated or result in subdivision into more groups (Anzai et al., 2000; Ait Tayeb et al., 2005; Yamamoto et al., 2000). The subgeneric groups of Moore et al. (1996) are used here only as points of reference for discussion of analyses of the other sequences (Young et al., 1992).

Sequences of Azotobacter species consistently formed a single heterogeneous group within the aeruginosa cluster, usually very close to P. aeruginosa and P. resinovorans. The nucleotide carA and recA sequences of Azomonas macrocytogenes were outliers to Pseudomonas but the derived peptide sequences were embedded in Pseudomonas. This can be attributed to the similarity between individual sequences of these genes and significant differences between nucleotides and peptides due to degeneracy in the second and third bases of codons. The influence of Azomonas macrocytogenes sequences appears to account for the different tree structures of the concatenated sequences in which Azotobacter and Azomonas macrocytogenes sequences are outliers to those of Pseudomonas. In general, however, individual and concatenated nucleotide and peptide sequence data with the 16S rRNA gene indicate a close relationship between the three genera, with Azotobacter species and Azomonas macrocytogenes being members of Pseudomonas.

Relationship of Azotobacter and Azomonas
The description of the genus Azotobacter (Kennedy et al., 2005) is of coccoid to rod-shaped cells that are up to 2.0 µm in diameter and 4 µm in length, but can be filamentous, up to 60 µm in length; if motile then by peritrichous flagella, and are capable of fixing nitrogen in aerobic conditions. Cells can form cysts. Thompson & Skerman (1981) proposed a separate genus for Azotobacter paspali, as Azorhizophilus paspali, but Kennedy et al. (2005) did not recognize this proposal and included the species in Azotobacter. The description of Azomonas (Kennedy & Rudnick, 2005) is of ellipsoidal to rod-shaped cells that are more than 2 µm in diameter, usually 2.5–3.5 µm in length, motile by peritrichous, lophotrichous or polar flagella, produce water-soluble fluorescent pigments, and are capable of fixing nitrogen under aerobic conditions. A consideration of the phenotypic determinative characteristics of the genera supports the synonymy of these genera as indicated by the sequence data. An investigation of these two genera, with a wider selection of strains, including the type species of Azomonas (Azomonas agilis), and perhaps based on other genes, might support their synonymy. If they are established as a single genus, Azotobacter would take priority over Azomonas.

Synonymy of Azotobacter with Pseudomonas?
At the molecular level, the similarity of Azotobacter species housekeeping gene and 16S rRNA gene sequences with those of Pseudomonas point to an intimate relationship between these taxa. Rediers et al. (2004), who studied a large selection of genes, concluded that Azotobacter vinelandii was very closely related to Pseudomonas. Kennedy et al. (2005) also noted the close relationship of these genera. Insofar as the peptide trees give expression to phenotypic relationships these also confirm the intimate relationship of the genera although these are not always congruent with the nucleotide trees. Although Palleroni (2005) did not describe nitrogen fixation as a characteristic in Pseudomonas, Chan et al. (1994) suggested that this capability could be widespread in the genus. That it has not been confirmed widely is perhaps because there have been few comprehensive investigations. It was reported for P. stutzeri (Chan et al., 1994; Vermeiren et al., 1999). Barraquio et al. (1986) and Kulakov et al. (2002) attributed a nitrogen-fixing capability to P. saccharophila, but this species has been transferred to a new and unrelated genus, Pelomonas (Xie & Yokota, 2005). Line (1997) recorded P. acidovorans, P. alcaligenes and P. aureofaciens as nitrogen-fixing species, but identifications were made using determinative methods that might not have given reliable identifications. Jenni et al. (1989) claimed nitrogen-fixing capability for strains of P. pseudoflava, but their identification should also be confirmed.

The consideration of the question of synonymy of the two genera hinges on the structure of Pseudomonas itself. As presently defined, the sequence data reported here could be considered strong circumstantial evidence for synonymy, although comparative analyses of the 16S rRNA gene (Anzai et al., 2000), 16S rRNA gene and rpoB (Ait Tayeb et al., 2005), 16S rRNA gene (Moore et al., 1996) and gyrB and rpoD (Yamamoto et al., 2000) hint at subgeneric species groupings that could imply subdivision of the genus. If further studies showed that Pseudomonas should be represented more properly by a small number of species closely associated with P. aeruginosa then Azotobacter might be retained as a sister genus. While it is sometimes claimed that phenotypic data correlate with sequence data, a careful reading suggests that this is generally not the case. Comparisons have been made using a number of phenotypic datasets: ribotyping (Brosch et al., 1996), fatty acid content of whole-cell hydrolysates and phospholipid fractions (Vancanneyt et al., 1996a), SDS-PAGE of whole-cell protein (Vancanneyt et al., 1996b), westprinting (Tesar et al., 1996) and fatty acid profiling (Stead, 1992). None of these datasets indicate coherent relationships between themselves, and none indicate relationships that correspond to those indicated by sequence data. There is no substantive evidence at present to suggest that Pseudomonas should be represented more properly by a small number of species closely associated with P. aeruginosa, the condition necessary to support retention of Azotobacter as a separate genus.

If Azotobacter species were transferred to Pseudomonas as presently defined, they would be included in an aeruginosa cluster (Anzai et al., 2000; Moore et al., 1996; Yamamoto et al., 2000). Because of their very close proximity to P. aeruginosa, the type species, they would be considered to be authentic members of the genus. Key generic phenotypic determinative characteristics of Pseudomonas recorded by Palleroni (2005) that differentiate them from Azotobacter are their size (0.5–1.0x1.5–5.0 µm), cell morphology and motility (by one or more polar flagella; rarely non-motile). Apparently, no studies have been made to discover whether the diffusible fluorescent pigments produced by Pseudomonas and Azotobacter are chemically similar. Transfer would require emendation of the generic description to include peritrichous flagella insertion, cell pleomorphism and nitrogen-fixing capability.

The diversity of these phenotypic characteristics, morphology, flagella insertion and nitrogen-fixing capability, associated with a population of bacteria sharing core housekeeping genes, invites the consideration that they are peripheral to the core molecular structures of these bacteria. The form of flagella insertion and cell size are generally regarded as distinct generic, determinative characteristics, but these, and the capacity to fix nitrogen, may be adaptive characteristics that are stably integrated but nevertheless arise from transmissible elements. It has been clear for some time that the transfer of genetic material between unrelated bacteria is a widespread phenomenon, though less clear is how regularly such transfer leads to phenotypic expression (Young, 2001). If expression is common, as indicated here, the implications for bacterial systematics are far-reaching. The assumption of current classification is still that there will be some correspondence between bacterial taxa established on the basis of comparative analyses of conserved sequences, and their phenotypic circumscriptions (Murray et al., 1990; Stackebrandt et al., 2002). However, reticulated evolution (Karlin et al., 1997; Rivera et al., 1998; Doolittle, 1999) may be so pervasive that such correspondence fails too often to give support to this concept.

The genes were chosen for this study when such selections were largely arbitrary, on the basis that they were involved in some deep-seated metabolic activity and therefore were unlikely to be hypervariable; influenced by trivial adaptive or random mutation. In a study of bacterial species for which whole chromosomes have been sequenced, Zeigler (2003) compared 32 individual gene sequences as representatives of available chromosomes. Of those tested, genes such as recN, lig, dnaX, glyA and cysS gave best results, with atpD (23rd), recA (29th) and 16S rRNA (32nd) being relatively unrepresentative of the chromosome. atpD, carA and recA were used because primers could be developed from the limited data available at the time (Gaunt et al., 2001; Weir et al., 2004). In future, it should be possible to choose genes that better represent the taxa under consideration. Notwithstanding the caution indicated here, a conclusion of this study, that there is a significant lack of correspondence between relationships indicated by conserved metabolic genes and key phenotypic determinative characteristics, is not likely to be contradicted by a refinement of gene selection.

These data also reinforce a point that has been noted in the past, that tree structures are highly dependent on sequence selection and on the choice of tree-building algorithm (Ludwig & Schleifer, 1999). A consideration of recent issues of IJSEM shows that NJ is still the preferred algorithm for the interpretation of sequence data. More rapid computational methods now permit the application of other analyses. With this done, however, it is not clear what are appropriate criteria for deciding which algorithm is to be preferred. The tendency is probably to choose algorithms that produce trees that conform to the taxonomic inclinations of authors.

Natural selection acts at the molecular level to refine protein efficiency. Identical peptides can be encoded by different nucleotide sequences because of the redundancy in coding of the second and third bases of codons. This is shown by the greater variability shown in nucleotide trees compared with peptide trees (see Supplementary Figs S1 and S2 in IJSEM Online). Peptide sequences are the units of selection. Nucleotide sequences are the vehicles for transfer of the functional structure between generations. This leaves open the question as to which sequences, nucleotides or peptides, are most relevant in the consideration of the relationships of organisms.

Although one or more combined gene analyses may give what appear to be robust separations into distinct groups, there is no objective method of calibration to determine whether such groups merit recognition as distinct taxa (Young, 2001; Gevers et al., 2005). For comparisons using DNA–DNA reassociation, the arbitrary value of 70 % reassociation was set to discriminate individual species. An attempt to calibrate 16S rRNA gene sequence data with DNA–DNA reassociation data was made by Stackebrandt & Goebel (1994). Using data available then, they demonstrated that strains that had 16S rRNA gene sequence similarity values lower than 97 % were invariably from different species based on reassociation data. However, this co-relationship is of limited use because many species based on reassociation have 16S rRNA gene sequence similarities greater than 99 %. Using other genes or concatenated sequences does not give even this limited calibration. Determination of taxonomic divisions can only be made by some arbitrary criterion, such as a specified DNA–DNA reassociation value, or by the subjective decision of an author. This problem applies to all datasets gathered by any means.

General taxonomic considerations
In recent times, bacterial classification, especially at the generic level, has been based largely on comparative analyses of a single molecular sequence (usually the 16S rRNA gene) using a single algorithm (usually NJ). This approach has been considered to provide a key to the systematization of relationships at higher taxonomic levels in classifications that aim to demonstrate inferred historical, ancestral relationships; viz., phylogenies (Stackebrandt & Woese, 1984; Woese, 1987; Vandamme et al., 1996). Underpinning such phylogenetic classification is the assumption that change in these, necessarily conserved, sequences occurs with clock-like regularity (Mindell & Thacker, 1996). It has been clear for some time, however, that such interpretations are oversimple and sometimes misleading. Although the rRNA gene is ubiquitous, this is not of itself the basis for asserting that comparative analyses of this sequence alone indicate reliable relationships of organisms (Young, 2001; Young et al., 1992) or that any other gene can be assumed, a priori, to reflect relationships with accuracy.

Increasingly it has come to be understood that molecular analyses can be refined by using some form of consensus based on several genes. Such analyses have been restricted in the past by the technical difficulties of developing primers and of producing large numbers of sequences, which have now largely been overcome. The concept of multilocus sequence typing (MLST) (Maiden et al., 1998) (now called multilocus sequence analysis; MLSA) involves systematic selection of a number of genes to represent the chromosome, and analysis of these genes or their peptide products as concatenated sequences (Gevers et al., 2005). Clearly, studies using properly conducted MLSA will give improved taxonomic resolution in future systematics sequence studies.

The goal of MLSA is to reflect the overall similarities and differences in molecular structures of organisms. It is not often noticed that this approach is overtly phenetic (Sneath, 1989) and is explicitly in contrast with the phylogenetic approach. Phenetic classification is based on the overall similarities and differences of organisms. In practice, until now investigation has been confined to a relatively small proportion of the totality of an organism, but the intention is clear. Phenetic classification is natural and evolutionary in the sense that it gives expression to bacterial groupings determined by the interplay of natural selection on the phenotype as determined by the genome. It explicitly does not aim to show ancestral relationships, accepting the possibility that apparently diverse taxa may be the product of rapid changes between phylogenetically closely related organisms, or that the ancestrally distantly related organisms may appear to be closely related by convergent events that have left no molecular trace. Any dataset can be interpreted as phenetic or phylogenetic classification. Applied as phylogenetic classification, it makes the additional, unverifiable and historical, claim that the analyses indicate reliable ancestral relationships. Unverifiable, because without a fossil record there is no method by which phylogenies can be calibrated (Wilkins, 2002). Attempts at phylogenetic classification based on macro-morphological structures were abandoned as it came to be understood that such classifications were not adequately robust and were open to the logical contradiction that structures could not be construed as responding to change with time when they were necessarily highly conserved by natural selection (Davis & Hayward, 1963). Phylogenetic interpretations based on molecular morphological structures, with the assumed correlation between molecular structural change and time, have the same logical objection.

The Marsden Fund of the Royal Society of New Zealand under contract 97-LAN-LFS-0011 and the New Zealand Foundation for Research, Science and Technology provided financial support. H. M. Harman and N. Waipara are thanked for critically reading the manuscript.

References

Ait Tayeb, L., Ageron, E., Grimont, F. & Grimont, P. A. D. (2005). Molecular phylogeny of the genus Pseudomonas based on rpoB sequences and application for identification of isolates. Res Microbiol 156, 763–773.[Medline]

Anzai, Y., Kim, H., Park, J.-Y., Wakabayashi, H. & Oyaizu, H. (2000). Phylogenetic affiliation of the pseudomonads based on 16S rRNA sequence. Int J Syst Evol Microbiol 50, 1563–1589.[Abstract]

Ausubel, F. M., Brent, B., Kingston, R. E., Moore, D. D., Seidman, J. G., Smith, J. A. & Struhl, K. (1987). Current Protocols in Molecular Biology. New York: John Wiley.

Barraquio, W. L., Parde, B. C., Jr, Watanabe, I. & Knowles, R. (1986). Nitrogen fixation by Pseudomonas saccharophila Doudoroff ATCC 15946. J Gen Microbiol 132, 237–241.[Abstract/Free Full Text]

Brosch, R., Lefèvre, M., Grimont, F. & Grimont, P. A. D. (1996). Taxonomic diversity of pseudomonads revealed by computer interpretation of ribotyping data. Syst Appl Microbiol 19, 541–555.

Castresana, J. (2000). Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol 17, 540–552.[Abstract/Free Full Text]

Chan, Y.-K., Barraquio, W. L. & Knowles, R. (1994). N2-fixing pseudomonads and related soil bacteria. FEMS Microbiol Rev 13, 93–118.

Davis, P. H. & Hayward, V. H. (1963). Principles of Angiosperm Taxonomy. Edinburgh: Oliver & Boyd.

Doolittle, W. F. (1999). Phylogenetic classification and the universal tree. Science 284, 2124–2128.[CrossRef][Medline]

Euzéby, J. P. (1997). List of Bacterial Names with Standing in Nomenclature: a folder available on the Internet (updated: March, 2007; revised URL: ). Int J Syst Bacteriol 47, 590–592.[Abstract/Free Full Text]

Gaunt, M. W., Turner, S. L., Rigottier-Gois, L., Lloyd-Macgilp, S. A. & Young, J. P. W. (2001). Phylogenies of atpD and recA support the small subunit rRNA-based classification of rhizobia. Int J Syst Evol Microbiol 51, 2037–2048.[Abstract]

Gevers, D., Cohan, F. M., Lawrence, J. G., Spratt, B. G., Coenye, T., Feil, E. J., Stackebrandt, E., Van de Peer, Y., Vandamme, P. & other authors (2005). Re-evaluating prokaryotic species. Nat Rev Microbiol 3, 733–739

Guindon, S., Lethiec, F., Duroux, P. & Gascuel, O. (2005). PHYML Online – a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res 33, W557–W559

Hilario, E. M., Buckley, T. R. & Young, J. M. (2004). Improved resolution on the phylogenetic relationships among Pseudomonas species by the combined analysis of atpD, carA, recA and 16S rDNA. Antonie Van Leeuwenhoek 86, 51–64.[CrossRef][Medline]

Jenni, B., Isch, C. & Aragno, M. (1989). Nitrogen fixation by new strains of Pseudomonas pseudoflava and related bacteria. J Gen Microbiol 135, 461–467.[Abstract/Free Full Text]

Karlin, S., Mrazek, J. & Campbell, A. M. (1997). Compositional biases of bacterial genomes and evolutionary implications. J Bacteriol 179, 3899–3913.[Abstract/Free Full Text]

Kennedy, C. & Rudnick, P. (2005). Genus II. Azomonas Winogradsky 1938, 391AL. In Bergey's Manual of Systematic BacteriologyThe Proteobacteria, 2nd edn, vol. 2, part B, pp. 379–384. Edited by D. J. Brenner, N. R. Krieg, J. T. Staley & G. M. Garrity. New York: Springer.

Kennedy, C., Rudnick, P., MacDonald, M. L. & Melton, T. (2005). Genus III. Azotobacter Beijerinck 1901, 567AL. In Bergey's Manual of Systematic Bacteriology – The Proteobacteria, 2nd edn, vol. 2, part B, pp. 384–402. Edited by D. J. Brenner, N. R. Krieg, J. T. Staley & G. M. Garrity. New York: Springer.

Kersters, K., Ludwig, W., Vancanneyt, M., De Vos, P., Gillis, M. & Schleifer, K. H. (1996). Recent changes in the classification of the pseudomonads: an overview. Syst Appl Microbiol 19, 465–477.

Kulakov, L. A., McAlister, M. B., Kimberley, L. O., Larkin, M. J. & O'Hanlon, J. F. (2002). Analysis of bacteria contaminating ultrapure water in industrial systems. Appl Environ Microbiol 68, 1548–1555.[Abstract/Free Full Text]

Line, M. A. (1997). A nitrogen-fixing consortia associated with the bacterial decay of a wooden pipeline. Lett Appl Microbiol 25, 220–224.[CrossRef]

Ludwig, W. & Schleifer, K. H. (1999). Phylogeny of bacteria beyond the 16S rRNA standard. ASM News 65, 752–757.

Maiden, M. C. J., Bygraves, J. A., Feil, E., Morelli, G., Russell, J. E., Urwin, R., Zhang, Q., Zhou, J., Zurth, K. & other authors (1998). Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci U S A 95, 3140–3145.[Abstract/Free Full Text]

Mindell, D. P. & Thacker, C. E. (1996). Rates of molecular evolution: phylogenetic issues and applications. Annu Rev Ecol Syst 27, 279–303.[CrossRef]

Moore, E. R. B., Mau, M., Arnscheidt, A., Böttger, E. C., Hutson, R. A., Collins, M. D., van de Peer, Y., de Wachter, R. & Timmis, K. N. (1996). The determination and comparison of the 16S rRNA gene sequences of species of the genus Pseudomonas (sensu stricto) and estimation of the natural intrageneric relationships. Syst Appl Microbiol 19, 478–492.

Murray, R. G. E., Brenner, D. J., Colwell, R. R., De Vos, P., Goodfellow, M., Grimont, P. A. D., Pfennig, N., Stackebrandt, E. & Zavarzin, G. A. (1990). Report of the ad hoc committee on approaches to taxonomy within the Proteobacteria. Int J Syst Bacteriol 40, 213–215.[Free Full Text]

Palleroni, N. J. (1984). Genus I. Pseudomonas Migula 1894, 237AL. In Bergey's Manual of Systematic Bacteriology, vol. 1, pp. 141–199. Edited by N. R. Krieg & J. G. Holt. Baltimore: Williams & Wilkins.

Palleroni, N. J. (2005). Genus I. Pseudomonas Migula 1894, 237AL. In Bergey's Manual of Systematic Bacteriology – The Proteobacteria, 2nd edn, vol. 2, part B, pp. 323–379. Edited by D. J. Brenner, N. R. Krieg, J. T. Staley & G. M. Garrity. New York: Springer.

Pfennig, N. & Trüper, H. (2005). Family 1. Rhodospirillaceae Pfennig and Trüper 1971, 17AL. In Bergey's Manual of Systematic Bacteriology – The Proteobacteria, 2nd edn, vol. 2, part C, p. 17. Edited by D. J. Brenner, N. R. Krieg, J. T. Staley & G. M. Garrity. New York: Springer.

Rediers, H., Vanderleyden, J. & De Mot, R. (2004). Azotobacter vinelandii: a Pseudomonas in disguise? Microbiology 150, 1117–1119.[Free Full Text]

Rivera, M. C., Jain, R., Moore, J. E. & Lake, J. A. (1998). Genomic evidence for two functionally distinct genome classes. Proc Natl Acad Sci U S A 95, 6239–6244.[Abstract/Free Full Text]

Sneath, P. H. A. (1989). Analysis and interpretation of sequence data for bacterial systematics: the view of a numerical taxonomist. Syst Appl Microbiol 12, 15–31.

Stackebrandt, E. & Goebel, G. M. (1994). Taxonomic note: a place for DNA-DNA reassociation and 16S rRNA sequence analysis in the present species definition in bacteriology. Int J Syst Bacteriol 44, 846–849.[Abstract/Free Full Text]

Stackebrandt, E. & Woese, C. R. (1984). The phylogeny of the prokaryotes. Microbiol Sci 1, 117–122.[Medline]

Stackebrandt, E., Frederiksen, W., Garrity, G. M., Grimont, P. A. D., Kämpfer, P., Maiden, M. C. J., Nesme, X., Rosselló-Mora, R., Swings, J. & 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, 1043–1047.[Abstract]

Stead, D. E. (1992). Grouping of plant-pathogenic and some other Pseudomonas spp. by using cellular fatty acid profiles. Int J Syst Bacteriol 42, 281–295.[Abstract/Free Full Text]

Tesar, M., Hoch, C., Moore, E. R. B. & Timmis, K. N. (1996). Westprinting: development of a rapid immunochemical identification for species within the genus Pseudomonas sensu stricto. Syst Appl Microbiol 19, 577–588.

Thompson, J. P. & Skerman, V. B. D. (1981). Azotobacteraceae: the Taxonomy and Ecology of the Aerobic Nitrogen-fixing Bacteria. London: Academic Press.

Vancanneyt, M., Witt, S., Abraham, W.-R., Kersters, K. & Fredrickson, H. L. (1996a). Fatty acid content in whole-cell hydrolysates and phospholipid fractions of pseudomonads: a taxonomic evaluation. Syst Appl Microbiol 19, 528–540.

Vancanneyt, M., Torck, U., Dewettinck, D., Vaerewijck, M. & Kersters, K. (1996b). Grouping of pseudomonads by SDS-PAGE of whole-cell proteins. Syst Appl Microbiol 19, 556–568.

Vandamme, P., Pot, B., Gillis, M., De Vos, P., Kersters, K. & Swings, J. (1996). Polyphasic taxonomy, a consensus approach to bacterial systematics. Microbiol Rev 60, 407–438.[Abstract/Free Full Text]

Vermeiren, H., Willems, A., Schoofs, G., de Mot, R., Keijers, V. & Hai, W. (1999). The rice inoculant strain Alcaligenes faecalis A15 is a nitrogen-fixing Pseudomonas stutzeri. Syst Appl Microbiol 22, 215–224.[Medline]

Weir, B. S., Turner, S. J., Silvester, W. B. & Young, J. M. (2004). Unexpectedly diverse Mesorhizobium strains and Rhizobium leguminosarum nodulate native legume genera of New Zealand, while introduced legume weeds are nodulated by Bradyrhizobium species. Appl Environ Microbiol 70, 5980–5987.[Abstract/Free Full Text]

Wilkins, A. S. (2002). The Evolution of Developmental Pathways. Sunderland, MA: Sinauer Associates.

Woese, C. R. (1987). Bacterial evolution. Microbiol Rev 51, 221–271.[Free Full Text]

Xie, C.-H. & Yokota, A. (2005). Reclassification of Alcaligenes latus strains IAM 12599T and IAM 12664 and Pseudomonas saccharophila as Azohydromonas lata gen. nov., comb. nov., Azohydromonas australica sp. nov. and Pelomonas saccharophila gen. nov., comb. nov., respectively. Int J Syst Evol Microbiol 55, 2419–2425.[Abstract/Free Full Text]

Yamamoto, S., Kasai, H., Arnold, D. L., Jackson, R. W., Vivian, A. & Harayama, S. (2000). Phylogeny of the genus Pseudomonas: intrageneric structure reconstructed from the nucleotide sequences of gyrB and rpoD genes. Microbiology 146, 2385–2394.[Abstract/Free Full Text]

Young, J. M. (2001). Implications of alternative classifications and horizontal gene transfer for bacterial taxonomy. Int J Syst Evol Microbiol 51, 945–953.[Abstract]

Young, J. M. & Park, D.-C. (2007). Relationships of plant pathogenic enterobacteria based on partial atpD, carA, and recA as individual and concatenated nucleotide and peptide sequences. Syst Appl Microbiol 30, 343–354.[CrossRef][Medline]

Young, J. M., Takikawa, Y., Gardan, L. & Stead, D. E. (1992). Changing concepts in the taxonomy of plant pathogenic bacteria. Annu Rev Phytopathol 30, 67–105.

Zeigler, D. R. (2003). Gene sequences useful for predicting relatedness of whole genomes in bacteria. Int J Syst Evol Microbiol 53, 1893–1900.[Abstract/Free Full Text]