Previous Article | Next Article ![]()
Journal of Clinical Microbiology, April 2006, p. 1359-1366, Vol. 44, No. 4
0095-1137/06/$08.00+0 doi:10.1128/JCM.44.4.1359-1366.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
and
E. C. Böttger
Institut für Medizinische Mikrobiologie, Universität Zürich, 8006 Zürich, Switzerland
Received 13 December 2005/ Returned for modification 9 January 2006/ Accepted 24 January 2006
|
|
|---|
|
|
|---|
Except for Pseudomonas aeruginosa, accurate identification of gram-negative nonfermenters in the clinical laboratory mainly relies on commercially available phenotypic identification systems. Commercial systems for bacterial identification have been offered for approximately 30 years. Of the original manual systems, only the API 20 NE (bioMérieux, Marcy l'Etoile, France) has remained available, while fully or partly automated identifications systems, such as VITEK (bioMérieux, Marcy l'Etoile, France), PHOENIX (BD, Sparks, Md.), and MicroScan (Dade Behring, West Sacramento, Calif.), have now been implemented in many laboratories. These systems have contributed to a more effective management of patients by enabling microbiologists to identify corresponding bacteria of clinical relevance more rapidly and accurately (12). However, all phenotypic test systems have potential inherent problems, e.g., (i) not all strains within a given species may exhibit a particular characteristic (6), (ii) the same strain may give different results upon repeated testing (6), and (iii) the corresponding databases are limited (2). Specifically, nonfermenters recovered from cystic fibrosis patients may pose problems due to phenotypic variations and slower growth rates because of the significant antimicrobial pressure that these organisms face in the lungs of these patients (17).
Genotypic identification methods are emerging as an alternative or complement to established phenotypic identification procedures. For bacteria, 16S rRNA gene sequence analysis is a widely accepted tool for molecular identification (5, 18, 21). Public databases (GenBank, Nucleotide Sequence Database at the European Molecular Biology Laboratory, DNA Data Bank of Japan, RDP II) contain a vast number of bacterial 16S rRNA sequences, allowing for rapid analysis and providing phylogenetically meaningful information. Few studies so far have systematically compared molecular and phenotypic identification procedures to determine the usefulness of sequence-based methods for the diagnostic laboratory (2, 3, 8, 10, 11, 16, 22, 24, 25, 27, 28, 30); available studies mainly focused on mycobacteria (8, 16, 22, 25) and gram-positive microorganisms (2, 3, 28). Several limitations are present in those rare studies investigating aerobic gram-negative rods: one study was restricted to isolates which phenotypically gave problematic results (11), one study was restricted to isolates from cystic fibrosis patients (30), and one study was restricted to the commercially available MicroSeq system (27). Given these limitations, and with a view toward developing a diagnostic algorithm for implementation into the microbiology laboratory, we here compared in a prospective fashion phenotypic systems (API 20 NE, VITEK 2) with 16S rRNA gene sequencing for identification of clinically relevant isolates of aerobic nonfermenting gram-negative bacilli.
|
|
|---|
Identification using the API 20 NE system.
The API 20 NE system covers 61 nonenterobacterial gram-negative taxa. Testing was performed according to the instructions of the manufacturer (bioMérieux, Marcy l'Etoile, France). Substrate assimilations were read after 24 and 48 h. Interpretation of the results was done after 48 h using the identification software version 6.0. Strains were classified into one of the following 3 groups: (i) identification at species level, (ii) identification at genus level, (iii) no identification (i.e., low discrimination). According to the manufacturer's instructions, strain identification at the species level was divided into 4 subgroups: (i) excellent species identification, percent identification of
99.9% and T value of
0.75; (ii) very good species identification, percent identification of
99.0% and T value of
0.5; (iii) good species identification, percent identification of
90.0% and T value of
0.25; (iv) acceptable species identification, percent identification of
80.0% and T value of
0.0.
Identification using the VITEK 2 fluorescent system (ID-GNB card).
The VITEK 2 fluorescent system (ID-GNB card) includes 43 nonenterobacterial gram-negative taxa. Testing was performed according to the instructions of the manufacturer. Briefly, strains were cultured on Columbia sheep blood agar (Difco 279020) or MacConkey agar (BBL 4311391) for 18 to 24 h at 37°C before the isolate was subjected to analysis. Strains which had been stored at 70°C were subcultured twice before analysis. A bacterial suspension was adjusted to a McFarland standard of 0.50 to 0.63 in a solution of 0.45% sodium chloride using the VITEK 2 DensiCheck instrument (bioMérieux). The time between preparation of the solution and filling of the card was always less than 1 h. Analysis was done using the identification card for gram-negative bacteria (ID-GNB card) containing 41 fluorescent biochemical tests (13). Cards are automatically read every 15 min. Data were analyzed using the VITEK 2 software version VT2-R03.1. Isolates which were not identified were retested with a fresh subculture. According to the manufacturer's instructions, strain identification at the species level was divided into 4 subgroups: (i) excellent species identification, T value of
0.75; (ii) very good species identification, T value of
0.5 and <0.75; (iii) good species identification, T value of
0.25 and <0.5; (iv) acceptable species identification, T value
0.0 and <0.25.
Sequencing of 16S rRNA gene. Sequencing was performed as described previously (2). In brief, a loopful of bacterial cells was digested using lysozyme (Sigma-Aldrich Chemie GmbH, Schnelldorf, Germany) and alkaline lysis. Following nucleic acid purification, the 5' part of the 16S rRNA gene (corresponding to Escherichia coli positions 10 to 806) was amplified using primers BAK11w [5'-AGTTTGATC(A/C)TGGCTCAG] and BAK2 [5'-GGACTAC(C/T/A)AGGGTATCTAAT] (4). Cycling parameters included an initial denaturation for 5 min at 95°C; 40 cycles of 1 min at 94°C, 1 min at 48°C, and 1 min at 72°C; and a final extension for 10 min at 72°C. Five microliters of the DNA extract was used for amplification in a total volume of 50 µl containing 1.25 U of AmpliTaq DNA polymerase LD (Applied Biosystems, Rotkreuz, Switzerland). Amplicons were purified and sequenced with forward primer BAK11w (4), and the fragments were analyzed using an automatic DNA sequencer (ABI Prism 3100 Genetic Analyzer; Applied Biosystems).
Sequence analysis. 16S rRNA gene sequences were compared with those available in the GenBank, EMBL, and DDBJ databases using a two-step procedure. A first search was performed with the FASTA algorithm of the Wisconsin GCG program package (9). All positions showing differences to the best-scoring reference sequence were visually inspected in the electropherogram, and the sequence was corrected manually if necessary. Thereafter, a second search was done using BLASTN. Undetermined nucleotides (designated N) in either the determined sequence or the reference sequence were counted as matches. The mean length of the sequences after manual editing was 508 ± 104 nucleotides (range, 315 to 692 nucleotides) containing 0.8 ± 1.5 undetermined (N) positions (range, 0 to 9 N).
For identification, the following criteria were used: (i) when the determined sequence yielded a similarity score of
99% with a reference sequence of a classified species, the unknown isolate was assigned to this species; (ii) when the score was <99% and
95%, the unknown isolate was assigned to the corresponding genus; and (iii) when the score was <95%, the unknown isolate was assigned to a family. If the unknown isolate was assigned to a species and the second classified species in the scoring list showed less than 0.5% additional sequence divergence, this was categorized as a "species with low demarcation to next species."
Discrepant analysis. According to Stackebrandt and Goebel (26), 16S rRNA gene similarities of less than 97% indicate that strains belong to different species. If the isolate was assigned to a different species by phenotypic and molecular criteria, the procedure was as follows: (i) if the isolate's sequence showed less than 97% similarity to the sequence of the species assigned phenotypically, it was concluded that the isolate does not belong to the species identified by phenotypic means and the sequencing result was considered correct; (ii) if the isolate's sequence showed 97% to 99% similarity to the sequence of the species assigned by phenotypic means, the isolate was categorized as unresolved.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Molecular versus phenotypic identification for 107 isolates (unresolved data)
|
|
View this table: [in a new window] |
TABLE 2. Eighty-eight isolates for which molecular and phenotypic identification were concordant
|
|
View this table: [in a new window] |
TABLE 3. Molecular and phenotypic identifications: discrepant results (n = 19)
|
Identification using VITEK 2. VITEK 2 identified 57 of the 107 isolates to the species level, yielding excellent, very good, good, and acceptable species identifications in 37, 12, 4, and 4 cases, respectively. In 1 case, identification at the genus level was achieved; 49 isolates could not be identified (Table 1). In 12 cases identified at the species level, a mixed taxon was reported, with all of these representing Acinetobacter baumannii or A. calcoaceticus or Acinetobacter genomic species 3 (belonging to the A. baumannii-A. calcoaceticus complex).
Comparing 16S rRNA gene sequencing with API 20 NE. For 45 of the 58 strains identified to the species level by API 20 NE (for raw data, see Tables 2 and 3; for resolved data, see Table 4), molecular identification assigned the isolate to the same taxon (18/20, 15/18, 10/15, and 2/5 cases with excellent, very good, good, and acceptable species identification by API 20 NE, respectively). In 1/58 isolates, sequencing gave genus level assignment (Wautersia sp.) and the API system resulted in species level identification (Wautersia paucula).
|
View this table: [in a new window] |
TABLE 4. Molecular identification versus identification by API 20 NE (resolved data)
|
In 7 of the 12 discrepant cases, the isolates' 16S rRNA gene sequences exhibited
97% similarity to the 16S rRNA gene sequence of the species assigned by API 20 NE. According to Stackebrandt and Goebel (26), 16S rRNA gene similarities of less than 97% indicate that strains belong to different species. Although only partial sequences were used here, it was thus concluded that these isolates do not belong to the species identified by the API system. For example, one strain was identified as Burkholderia cepacia by API 20 NE; the strain's 16S rRNA gene sequence showed 100.0% similarity with Herbaspirillum huttiense (formerly Pseudomonas huttiensis) but less than 95% similarity with B. cepacia. This isolate clearly does not represent B. cepacia but is H. huttiense. In 4 of these 7 cases, the isolate belongs to a species not covered by the API 20 NE database.
In 5/12 discrepant cases, the 16S rRNA gene sequence of the isolate exhibited between 97% and 99% sequence similarity with the 16S rRNA gene sequence of the taxon assigned by API 20 NE (for further details, see Table 3). We can thus not exclude the possibility that the isolates are representatives of the species identified by the API system (unresolved cases). In one of these 5 cases, sequence analysis assigned the isolate to a species not included in the API 20 NE database.
Using API 20 NE, 7/107 strains investigated were identified at the genus level and 42/107 isolates could not be identified at any taxonomic level (Table 1). 16S rRNA gene sequencing allowed species assignment in 5 of the former and 37 of the latter cases (for raw data, see Tables 2 and 3; for resolved data, see Table 4). The remaining 7 strains were identified at the genus level by molecular means.
Comparing 16S rRNA gene sequencing with VITEK 2. Molecular identification assigned the isolate to the same taxon in 46 of 57 strains identified to the species level by VITEK 2 (33/37, 6/12, 3/4, and 4/4 cases with excellent, very good, good, and acceptable species identification by the VITEK 2, respectively) (for raw data, see Tables 2 and 3; for resolved data, see Table 5). Eleven of the 57 results were discrepant (Table 3). In 3/11 discrepant results, the isolates were identified merely to genus level by sequencing (the sequences showed less than 99.0% similarity to the best-scoring reference sequence); in 8/11 discrepancies, sequencing assigned the strains to a different species compared to VITEK 2.
|
View this table: [in a new window] |
TABLE 5. Molecular identification versus identification by VITEK 2 (resolved data)
|
Using VITEK 2, one strain was identified at the genus level; 16S rRNA gene sequencing allowed species identification of this isolate. Forty-nine of 107 isolates could not be identified at any taxonomic level using VITEK 2; sequencing allowed species assignment in 42 and genus assignment in 7 of these 49 isolates (for raw data, see Tables 2 and 3; for resolved data, see Table 5).
|
|
|---|
16S rRNA gene sequencing is more accurate for the identification of gram-negative nonfermenters than API 20 NE and VITEK 2: 92% of isolates were assigned at the species level by sequencing, compared to 54% and 53% by API and VITEK 2, respectively. In 12 of 45 (API 20 NE) and 11 of 57 (VITEK 2) isolates with identification to the species level by phenotypic procedures, sequence analysis yielded discrepant results. For API 20 NE, we concluded that sequencing was correct for 7/12 strains, and 5/12 strains remained unresolved. For VITEK 2, it was concluded that the molecular approach delivered the correct result for 9/11 strains, and 2 strains remained unresolved.
For the API 20 NE, a result seems reliable when a species assignment is reported as an excellent or very good identification; 18 of 20 and 15 of 18 cases in these categories were correctly identified when compared to the results of sequencing. In the categories "good species identification" and "acceptable species identification," 10 of 15 cases and 2 of 5 cases were correct, respectively. Thus, if identification is of lower quality, the result is less reliable and the isolate should be subjected to sequence analysis if accurate identification is of concern. For the VITEK 2 fluorescent system, excellent species identifications are reliable; results of 33 of 37 cases in this category were consistent with molecular analysis. For the other categories, the numbers of isolates were too small to allow for firm conclusions.
One major problem of phenotypic test systems is that the available databases are limited. Of the 107 strains analyzed in the present study, 16 and 46 of the isolates corresponded to species not included in the API 20 NE and VITEK 2 fluorescent card databases, respectively. In other words, while the API 20 NE covers the majority of nonfermenters isolated, the VITEK 2 database needs to be expanded. More recently, efforts have been undertaken to enlarge the VITEK 2 database. The revised VITEK 2 database is associated with a new colorimetric detection card and covers 159 taxa (versus 101 for the ID-GNB card). Twenty of the 46 strains not identified in our study due to noninclusion of the corresponding species in the database are now included in the expanded VITEK 2 database. Evaluation with a large collection of 655 gram-negative strains (including 144 nonfermenters) gave encouraging results (12). Further studies, testing the system under routine laboratory conditions, and comparing the results with genotypic methods are required before firm conclusions can be drawn on the quality of the enlarged VITEK 2 database.
16S rRNA gene sequencing has distinct benefits compared to phenotypic identification procedures: (i) it is not restricted to a specific group of bacteria as public databases such as GenBank cover the whole spectrum of phylogenetic diversity; (ii) novel, not yet described species can be assigned to a group of related bacteria; (iii) results are in general unambiguous and not dependent on strain variation or individual interpretation. We have recently demonstrated the potential of partial 16S rRNA gene sequencing for identification of gram-positive microorganisms and we have proposed algorithms for integration of molecular identification procedures into the diagnostic work flow (2, 3). In the present study, we have extended our efforts to include 16S rRNA gene sequence analysis for accurate identification of gram-negative nonfermenting bacilli.
16S rRNA gene sequencing has some limitations. It is in part compromised by a low phylogenetic power at the species or subspecies level (26). In the present study, 35% of isolates could not be unambiguously assigned to a single species. This was particularly true for members of the Burkholderia cepacia complex, the Acinetobacter baumannii-A. calcoaceticus complex and close relatives, some members of the genus Pseudomonas, the genus Achromobacter, and the genus Ralstonia. For the Burkholderia cepacia complex, the recA gene seems more appropriate for recognizing the different genomovars (19). For Acinetobacter, complete 16S rRNA gene sequencing has been proposed for identification (15). More recently, the 16S-23S intergenic spacer region has been shown to differentiate closely related bacteria within the Acinetobacter baumannii-A. calcoaceticus complex (7).
The quality of public databases, such as GenBank, is critical. Sequences are deposited independently of their quality, e.g., regardless of the correct assignment, the length of the sequence, or the number of ambiguous nucleotides. Of particular concern is that sequences in public databases may be assigned to a naming which possibly no longer is valid due to taxonomic changes or which never has been validly published before. As an example, several sequences are annotated under Acinetobacter calcoaceticus subsp. anitratus in GenBank; however, this is an old denomination which phenotypically corresponds to members of the Acinetobacter baumannii-A. calcoaceticus complex (1). Correct interpretation of sequence data also requires some familiarity with taxonomy and recent taxonomic changes. As an example, a sequence of an isolate with 100% homology to Burkholderia multivorans but 12 mismatches to Burkholderia cepacia should not be reported as B. multivorans per se. Rather, the result reported should reflect that B. multivorans is regarded as a member of the B. cepacia complex; B. multivorans is the second most common species of the B. cepacia complex in cystic fibrosis infections (20).
Despite all these minor limitations, 16S rRNA gene sequencing does not lead to false identification, and with some knowledge about taxonomy, a sequence can unambiguously be assigned. In this study, partial 16S rRNA gene sequencing has been compared with two commercially available systems (API 20 NE, VITEK 2) for identification of nonfermenting gram-negative rods. The majority of strains could not be accurately identified by phenotypic profiling, as species assignment was found to be reliable only when an excellent (or very good) species identification according to the system's criteria was achieved; this was the case in 35% of the isolates. In practical terms, molecular identification is more laborious than phenotypic identification; results of sequencing are usually available within one to two working days. Given these considerations, we have developed an algorithm for the effective and proper identification of gram-negative nonfermenters in the diagnostic laboratory (Fig. 1). Thus, when API 20 NE or VITEK 2 do not yield an excellent (or very good) species identification, nonfermenters should be subjected to 16S rRNA gene sequencing if adequate species assignment is of concern.
![]() View larger version (31K): [in a new window] |
FIG. 1. Algorithm for the identification of nonfermenting gram-negative bacilli.
|
This study was supported by the University of Zürich.
Present address: Bio-Analytica AG, 6000 Luzern 6, Switzerland. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»