Previous Article | Next Article ![]()
Journal of Clinical Microbiology, July 2002, p. 2653-2658, Vol. 40, No. 7
0095-1137/02/$04.00+0 DOI: 10.1128/JCM.40.7.2653-2658.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Department of Microbiology and Cancer Research Institute, Institute of Endemic Diseases, SNUMRC, Seoul National University College of Medicine, and Clinical Research Institute, Seoul National University Hospital, Seoul 110-799,1 Laboratory of Rickettsial and Zoonotic Disease, Department of Microbiology, Korean National Institute of Health, Seoul 122-701, Korea,2 Department of Microbiology, School of Medicine, University of Occupational and Environmental Health, Kitakyushu 807, Japan3
Received 21 December 2001/ Returned for modification 4 March 2002/ Accepted 23 April 2002
|
|
|---|
|
|
|---|
The rpoB gene encodes the ß subunit of DNA-dependent RNA polymerase (24), and rifampin resistance is related to mutations in a particular region of rpoB (10). Recently, rpoB sequences were used as an alternative tool either for determining the phylogeny of or for identifying enteric bacteria (16), Mycobacterium (10), spirochetes including Borrelia (14, 21), and Bartonella (22). In this study, partial rpoB sequences (300 bp), containing a region that though highly conserved still has a remarkable ability for species differentiation (10, 12, 14, 22), were determined for genotypic classification of Legionella species and the results were compared with those results from the 16S rRNA gene and mip. In addition, culture isolates of Legionella species were identified based on rpoB sequences.
Fifty-five reference strains of 38 Legionella species were used in this study (Table 1). Among the 15 serogroups of L. pneumophila, three strains of serogroups, i.e., 4, 5, and 15, were L. pneumophila subsp. fraseri, while the others were L. pneumophila subsp. pneumophila. One hundred and five culture isolates, of which 10 strains were isolated from pneumonia patients and the others from cooling water, were identified by rpoB sequence analysis. These isolates were grown on buffered charcoal yeast extract agar and identified on the basis of cysteine requirement, autofluorescence, and biochemical tests such as those for gelatinase, urease, catalase, oxidase, peroxidase, and ß-lactamase activities; hippurate hydrolysis; nitrate reduction; and browning of tyrosine-supplemented agar (8). They were provided by Mi-Yeoun Park at the Korean National Institute of Health and H. Miyamoto at the University of Occupational and Environmental Health for blinded identification by rpoB sequence analysis.
|
View this table: [in a new window] |
TABLE 1. Legionella strains used to determine the rpoB sequences in this study
|
Sequences of the purified PCR products were directly determined with forward and reverse primers using an Applied Biosystems automated sequencer (model 377) and a BigDye Terminator Cycle Sequencing kit (Perkin-Elmer Applied Biosystems, Warrington, United Kingdom). For the sequencing reaction, 30 ng of purified PCR products, 2.5 pmol of each primer, and 4 µl of BigDye Terminator RR mix (part no. 4303153; Perkin-Elmer Applied Biosystems) were mixed and adjusted with distilled water to a final volume of 10 µl. The reaction was run with 5% (vol/vol) dimethyl sulfoxide for 30 cycles of 15 s at 95°C, 5 s at 50°C, and 4 min at 60°C. To confirm the identity of the reference strains used in this study, we also amplified and sequenced about 400 bp from the 3' end of the 16S rRNA gene (8). In addition, mip sequence analysis (402 bp) for 79 Korean isolates was performed to compare with the result of rpoB sequence analysis. For the amplification and sequencing of mip DNA, a new primer set, ML1 (5'-GAT AAG TTG TCT TAT AGC ATT GG-3') and ML2 (5'-TCT GTC CAT CCT GGG ATA ACT TG-3'), was used.
The partial rpoB sequences (300 bp) were aligned using the multiple alignment program CLUSTAL X (29) and the sequences determined were submitted to GenBank (Table 1). Phylogenetic trees of Legionella were constructed by the neighbor-joining (NJ) method and the parsimony method in the PAUP program (28). In the NJ method, pairwise distances were calculated using the maximum likelihood option, and in parsimony analysis, heuristic searches were conducted with the option of tree bisection reconstruction branch swapping. C. burnetii, the etiologic agent of Q fever, was used as an outgroup. Branch supporting values were evaluated with 1,000 bootstrap replications. The 16S rRNA gene (1,396 bp) and mip (525 bp) sequences of Legionella were retrieved from the GenBank database and used to infer the phylogenetic relationships by the same method. Incongruence length difference (ILD) tests (partition homogeneity test in the PAUP program) (3) were conducted to determine whether rpoB, the 16S rRNA gene, and mip data sets were coalescent together or not. The degree of incongruence between trees was assessed by comparing the log likelihood values (11) of three NJ trees with those of the phylogenetic topologies obtained from the other genes (7). All analyses were performed using the PAUP program (28).
Similarities between the partial rpoB DNA sequences (300 bp) of Legionella species were lower than 95% in all cases except for L. jamestownensis-L. londiniensis (98.7%). The sequence divergence in rpoB was 3.5 times greater than that for the 16S rRNA gene, indicating the greater variance of rpoB. The mip sequences were more divergent than those of rpoB by a factor of 1.5. However, the relationship between the pairwise dissimilarities of rpoB and other genes was not linear, which suggested that the degree of sequence divergences was not proportional to the time needed for speciation. In contrast, the pairwise distances between the mip and 16S rRNA genes exhibited a close linear correlation (data not shown). The homoplasy index of the rpoB sequences (0.779), calculated by parsimony analysis using the PAUP program, was higher than those of the 16S rRNA gene and mip sequences (0.558 and 0.693, respectively).
The phylogeny inferred from rpoB sequences has long terminal branches with the exception of L. jamestownensis and L. londiniensis. However, bootstrap values supporting each branch in the rpoB phylogeny were relatively lower (Fig. 1) than those of the 16S rRNA gene and mip phylogenies (Fig. 2). Although rpoB, 16S rRNA gene, and mip phylogenies demonstrated partially similar relationships, they represented different topologies in many respects (Fig. 1 and 2). According to the results of ILD tests (3) of PAUP (28), three of the gene data sets were incongruent with each other (P < 0.001), indicating that three gene sequences could not be combined. Log likelihood tests (11) also indicated that the tree topology of one gene differed significantly from those of the other two genes (P < 0.0001).
![]() View larger version (40K): [in a new window] |
FIG. 1. Phylogenetic relationships of Legionella species inferred from partial rpoB DNA sequences. This tree was constructed by the NJ method in PAUP (28). The species C. burnetii was used as an outgroup. The bootstrap values presented at corresponding branches were evaluated from 1,000 replications. Values below 50% are not indicated. The autofluorescence (BW, blue-white; YG, yellow-green; R, red), ubiquinone (A to F), and fatty acid (I to IV) groups of each Legionella species (8, 9, 13, 20, 30) are indicated at the right. Biochemical traits identical to phylogenetic grouping are represented in boxes. The scale bar represents 5 substitutions per 100 nucleotides.
|
![]() View larger version (32K): [in a new window] |
FIG. 2. Gene trees based on 16S rRNA gene sequences (A) and mip sequences (B). These were constructed using the same method as described above for Fig. 1 for the rpoB sequences. The mip sequence of L. geestiana was not available in GenBank. B, ubiquinone group B; R, red autofluorescence group; and B-BW, ubiquinone group B and blue-white autofluorescence group. The scale bars represent the number of inferred substitutions per site.
|
There are eight blue-white autofluorescent species in the genus Legionella: L. anisa, L. parisiensis, L. bozemanii, L. cherrii, L. steigerwaltii, L. gormanii, L. dumoffii, and L. tucsonensis (8, 20). All of these species belong to ubiquinone group B (13, 20). In the mip tree (Fig. 2B), the clade formation of these eight blue-white autofluorescent species with L. wadsworthii was robustly supported by bootstrap analysis (97%), but there was no equivalent single clade in the rpoB (Fig. 1) and 16S rRNA gene phylogenies (Fig. 2A).
The other species belonging to ubiquinone group B, L. cincinnatiensis, L. santicrucis, L. longbeachae, L. sainthelensi, and L. gratiana (13, 20), all grouped well into a single clade in both the rpoB (Fig. 1) and mip (Fig. 2B) phylogenies. While their relationships were supported only moderately by the rpoB phylogeny (bootstrap value of 73%), they were fully supported by the mip phylogeny. In the 16S rRNA gene tree (Fig. 2A), L. gratiana displayed a close relationship with one of the subgroups in the blue-white autofluorescence group, though with very low bootstrap support. The other four species formed a clade that was fully supported (bootstrap value of 100%) in the 16S rRNA gene tree (Fig. 2A).
In the rpoB tree, L. maceachernii, L. micdadei, and L. nautarum formed one clade (Fig. 1), which was identical in the mip phylogeny (Fig. 2B). These three species belong to the same ubiquinone group (D) and fatty acid group III (9, 13, 20). While their relationships were poorly supported by a very low bootstrap value in the rpoB, the support was much stronger as determined by the mip phylogeny. However, in the 16S rRNA gene phylogeny, rather than clustering with L. maceachernii and L. micdadei, L. nautarum formed a distinct group along with L. israelensis (Fig. 2A).
The red autofluorescent species, L. erythra and L. rubrilucens, formed a distinct clade in all analyses (R group in Fig. 2), showing a particularly high similarity in terms of 16S rRNA gene sequences (99.4%), though the sequence similarities in rpoB and mip were not as high (88.7 and 89.0%, respectively). The yellow-green autofluorescent species, L. birminghamensis and L. wadsworthii (30), were clustered into one group in the rpoB phylogeny (Fig. 1) but not in the other two trees. These two species with yellow-green autofluorescence showed a sequence similarity of 80.7%.
L. geestiana, which contains unique isoprenoid quinone (Q-14) (8, 13) and fatty acid compositions (5), was placed at a basal position with respect to all Legionella species in both rpoB and 16S rRNA gene phylogenies (Fig. 1 and 2A). For rpoB sequences, the sequence similarities between L. geestiana and all other species were below 80%.
Species identification of the 105 strains was accomplished using rpoB sequence analysis of reference strains to measure the similarities and thereby infer the phylogenetic reconstruction (data not shown). Of these, 98 isolates belonged to L. pneumophila and exhibited 99.3 to 100% rpoB sequence similarities with the reference strains. The others were identified as other Legionella species (two L. bozemanii, two L. dumoffii, two L. feeleii, and one L. micdadei species), which showed 98.3 to 100% sequence similarities with each reference strain. It was interesting that 87 strains could be identified as L. pneumophila subsp. pneumophila and 11 as L. pneumophila subsp. fraseri. These strains of L. pneumophila subsp. pneumophila and subsp. fraseri showed 97.0 to 100% and 98.3 to 100% sequence similarities within each subspecies, respectively. The identified strains that belong to each subspecies showed 87.3 to 89.7% sequence similarity. However, no correlation was observed between the rpoB sequence and those of serogroups. For example, rpoB sequences of four reference strains belonging to serogroup 1 of L. pneumophila were not identical. On the other hand, mip sequences of 97 L. pneumophila strains, including 18 reference strains and 79 isolates, showed 93.4 to 100% similarities. However, unlike in rpoB analysis, no clear differentiation between two subspecies of L. pneumophila was observed in mip sequence analysis.
Recently, it was demonstrated that the rpoB-based approach to microbial community analysis or identification could practically overcome the inherent limitations of 16S rRNA gene intraspecies heterogeneity (4). In this study, we have shown that rpoB sequences are as useful as 16S rRNA gene and mip sequences in accessing and evaluating the relationships among the Legionella species. Such a simple genetic analysis promises to provide a practical scheme for the classification of Legionella species and the identification of culture isolates in the diagnostic or reference clinical laboratory.
In contrast to analysis of other genes, rpoB analysis has several advantages. In spite of the sequence's shortness (300 bp), rpoB sequence analysis was able to clearly differentiate among the Legionella species. One such example was the clear differentiation of L. erythra and L. rubrilucens by rpoB analysis. Because of their high similarity in 16S rRNA gene sequences (99.4%) and their DNA hybridization test results, the latter representing a value of just below 70% (23), they had long been considered subspecies of the same species (8). Another example was the differentiation of blue-white autofluorescent species. While their 16S rRNA gene sequences indicated similarities of at least 97.8% (8), their rpoB sequences represented 81.7 to 95.3% similarities. A further example, which might be the most significant one, was the differentiation of subspecies. L. pneumophila includes two subspecies (L. pneumophila subsp. pneumophila and subsp. fraseri), which featured 99.2% similarity (8) with 16S rRNA gene sequence and were not differentiated by mip analysis (19). However, rpoB analysis was able to clearly distinguish them (12). Coupled with the high nucleotide similarity within each subspecies (97.0 to 100% and 98.3 to 100%), the significant difference of nucleotides between the two subspecies (87.3 to 89.7%) indicates that rpoB can be a better marker in differentiation of L. pneumophila than can mip, which could not discriminate two subspecies of L. pneumophila. These results confirm the usefulness of rpoB sequence analysis in population genetic and epidemiological studies of L. pneumophila, including molecular typing.
The rpoB tree showed quite a different topology from those of the 16S rRNA gene and mip sequences (Fig. 1 and 2). The ILD and log likelihood tests suggested that phylogenetic relationships inferred from the three genes were statistically significantly different. The lack of congruence between the rpoB, 16S rRNA gene, and mip trees did not result only from frequent synonymous base substitutions. Otherwise, the observed incongruence among the three gene trees is most likely the result of horizontal gene transfer between Legionella species (7, 15, 25, 26). It is known that such horizontal gene transfer disrupts the treelike branching pattern, thus complicating the phylogenetic relationships of species (7).
Due to these discrepancies among the three gene trees, none was able to exactly explain the phylogenetic relationships within the Legionella species and the evolution of this species. This suggests that the use of several markers, such as rpoB in combination with 16S rRNA gene or mip, may be necessary for the reliable identification and phylogenetic study of Legionella. Such an approach will reduce the risk of error in molecular typing or identification.
We thank H. S. Jung (School of Biological Sciences, Seoul National University, Seoul, Korea) for his help with the phylogenetic analysis and critical reading of the manuscript. We also thank B. S. Field, R. F. Benson, and E. Brown (National Center for Infectious Diseases, Centers for Disease Control and Prevention, Atlanta, Ga.) and M. J. Kim (College of Medicine, Korea University, Seoul, Korea) for providing reference strains and clinical isolates.
This work was supported by a grant of the Korea Health 21 R&D Project, Ministry of Health and Welfare, Seoul, Republic of Korea (01-PJ10-PG6-01GM03-0002), and in part by the BK21 project for Medicine, Dentistry, and Pharmacy.
|
|
|---|
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»