Previous Article | Next Article ![]()
Journal of Clinical Microbiology, July 2003, p. 2952-2960, Vol. 41, No. 7
0095-1137/03/$08.00+0 DOI: 10.1128/JCM.41.7.2952-2960.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Department of Microbiology,1 Department of Rehabilitation Medicine, College of Medicine, Konkuk University, Chungju, Chungchongbuk-Do 380-701,2 Department of Microbiology, National Institute of Health, Seoul 122-701,3 Department of Microbiology, College of Medicine, Hallym University, Chunchon, Kangwon-Do 200-702,4 Department of Microbiology, Seoul National University College of Medicine, Seoul 110-799, Korea5
Received 26 December 2002/ Returned for modification 12 March 2003/ Accepted 16 April 2003
|
|
|---|
|
|
|---|
Rickettsiae have been characterized conventionally and identified by serotyping and protein analysis by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (15, 17). Since rickettsiae are difficult to grow because of their strictly intracellular nature, these methods are time-consuming and expensive. Also, cross-reactions with different species and different genera by serotyping have been reported (8, 14, 18).
Various DNA-based techniques have been developed for the species identification of rickettsiae. Phylogenetic analyses based on the 16S rRNA gene sequence have frequently been used, but since the sequences are highly conserved, significant inferences about intragenus phylogeny are not possible (15, 16). Phylogenetic analyses based on the citrate synthase-encoding gene (gltA) and the outer membrane protein-encoding gene (ompA) have been developed to find a more sensitive and significant phylogenetic relationship among the rickettsiae (4, 16). Comparison of gltA gene sequences is more sensitive than comparison of 16S rRNA gene sequences, but reliable phylogenetic interrelationships could be drawn only among the rickettsiae that diverged early from the common ancestor of SFG (4). The interspecies differences of the ompA gene are greater than those of both the 16S rRNA and the gltA genes among the SFG, but the gene sequences of some species (R. akari, R. australis, R. helvetica, R. belii, R. canada, R. typhi, R. prowazekii) were not available for testing (4, 15). Therefore, a new identification method that counteracts the disadvantages of these methods is needed.
Recently, the groEL gene was used to identify Streptococcus (5), Enterococcus (6), Staphylococcus (7), Bartonella (12), Mycobacterium (20), and Ehrlichia (21). The groEL genes, which encode the 60-kDa heat shock protein GroEL, are ubiquitous in both prokaryotes and eukaryotes and encode highly conserved housekeeping proteins which are essential for the survival of these cells. The groEL genes provide the defining evolutionary relationships among the members of the eubacterial lineage (7, 12).
In this study, we developed a new method for the identification of the rickettsiae based on groEL gene analysis. The nucleotide sequences (534 to 546 bp) of the groEL genes from 15 rickettsial strains were determined and compared. Korean isolates were identified by comparison of their groEL gene sequences with those of reference strains. Also, we developed a rapid identification method using PCR-restriction fragment length polymorphism (RFLP) analysis and primer sets specific for species of the STG, SFG, and TG rickettsiae based on their groEL gene sequences.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Rickettsial strains used in this study
|
DNA was extracted with a High Pure PCR Template Preparation kit (Roche Diagnostics Co., Indianapolis, Ind.).
Nucleotide sequencing. PCR was performed with a set of primers (primer RF [5'-GTTGAAGTT/AGTTAAAGG-3'; positions 650 to 666 in the R. rickettsii R numbering] and primer RR [5'-TTTTTCTTTT/ATCATAATC-3'; positions 1183 to 1166]) to amplify the groEL gene (534 to 546 bp). Template DNA (50 ng) and 20 pmol of each primer were added to a PCR mixture tube (AccuPower PCR PreMix; Bioneer, Daejeon, Korea), which contained 1 U of Taq DNA polymerase, 250 µM each deoxynucleoside triphosphate, 50 mM Tris-HCl (pH 8.3), 40 mM KCl, 1.5 mM MgCl2, and gel loading dye. The volume was adjusted to 20 µl with distilled water. The reaction mixture was subjected to 30 cycles of amplification (30 s at 94°C, 45 s at 44°C, and 45 s at 72°C), followed by a 5-min extension at 72°C (model 9600 thermocycler; Perkin-Elmer Cetus). The PCR products were electrophoresed on a 1.5% agarose gel and purified with a QIAEX II Gel Extraction kit (QIAGEN, Hilden, Germany). The purified PCR product (20 ng) was ligated with 50 ng of the pGEM-T vector (Promega) at 16°C for 18 h and used for the transformation of Escherichia coli XL-1 Blue. Recombinant DNA was extracted with a High Pure Plasmid Isolation kit (Boehringer Mannheim, Indianapolis, Ind.). At least two clones of each strain were used for nucleotide sequence determination. The purified DNA was sequenced with the CEQ L DNA Analysis system and the CEQ 2000 Dye Terminator Cycle Sequencing kit (Beckman Coulter Inc., Fullerton, Calif.) with primers RF and RR and a forward and reverse sequencing primer (M13). For the sequencing reaction, 20 ng of purified DNA, 100 pmol of primer, 2 µl of 10x sequencing reaction buffer, 1 µl of the deoxynucleoside triphosphate mixture, 2 µl each of ddUTP, ddGTP, ddCTP, and ddATP, and 1 µl of DNA polymerase enzyme were mixed together; and the final volume was adjusted to 20 µl by adding distilled water. The reaction was performed for 30 cycles of 20 s at 96°C, 20 s at 50°C, and 4 min at 60°C.
Sequence analysis. Sequences were aligned by using the multiple-alignment algorithm in the MegAlign software package (Windows version 3.12e; DNASTAR, Madison, Wis.). Phylogenetic trees were constructed by the neighbor-joining method and the unweighted pair group method with arithmetic averages (UPGMA) by using the MEGA program (11). A bootstrap analysis (100 repeats) was performed to evaluate the topology of the phylogenetic tree.
PCR-RFLP analysis with DdeI digestion. PCR was performed with the same set of primers described above (primers RF and RR) to amplify the groEL DNA (534 to 546 bp), as described above. Endonuclease DdeI (New England Biolabs, Beverly, Mass.) was used to cleave the PCR products, as recommended by the manufacturer. The restriction fragments were electrophoresed on a 3% Amplisize agarose gel (Bio-Rad Lab., Hercules, Calif.) for 15 min at 100 V.
PCR for species-specific identification of rickettsiae. Primer sets specific for species of the TG, the STG, and the genus Rickettsia (SFG and TG) were designed on the basis of the groEL DNA sequences determined. Template DNA (50 ng) and 20 pmol of each primer were added to a PCR mixture tube (AccuPower PCR PreMix; Bioneer), and the volume was adjusted to 20 µl with distilled water. See Table 3 for the PCR amplification conditions.
|
View this table: [in a new window] |
TABLE 3. Intragroup variations of groEL and 16S rRNA gene sequences of rickettsial strains
|
|
|
|---|
More than 91.6% similarity was observed among the groEL gene sequences of strains of the genus Rickettsia (SFG and TG). The groEL gene sequence similarities of strains of the genus Rickettsia (SFG and TG) to STG strains and Ehrlichia strains were 55.6 to 61.8 and 52.4 to 55.6%, respectively (Table 2). SFG strains showed 91.6 to 93.3% similarities to TG strains, but they showed lower levels of similarity (55.6 to 60.7%) to STG strains. The sequence similarities of STG strains to SFG strains, TG strains, and Ehrlichia strains were 55.6 to 60.7, 57.5 to 61.8, and 44.1 to 47.0%, respectively.
|
View this table: [in a new window] |
TABLE 2. Levels of similarity between groEL sequences and between 16S rRNA gene sequences
|
The similarities of the groEL and 16S rRNA gene sequences of the SFG, TG, and STG rickettsiae are shown in Table 3. The groEL gene sequence similarities among the SFG, TG, and STG rickettsiae were 95.3 to 99.6, 97.0, and 95.4 to 99.1%, respectively; and the 16S rRNA gene sequence similarities among the SFG, TG, and STG rickettsiae were 97.4 to 99.4, 99.4, and 99.2 to 99.7%, respectively.
Comparison of phylogenetic trees determined from the 16S rRNA and groEL sequences. Phylogenetic trees based on the groEL gene sequences were constructed by UPGMA and the neighbor-joining method (data not shown) to investigate the relationships among the rickettsial strains (Fig. 1). The rickettsial strains could be distinguished from Ehrlichia strains. Five SFG strains (R. rickettsii, R. japonica, R. conorii, R. sibirica, and R. akari) and two TG strains (R. typhi and R. prowazekii) formed a cluster separate from the STG strains. Also, four STG strains (O. tsutsugamushi Kato, Karp, Kawasaki, and Gilliam) formed a cluster separate from strains of the genus Rickettsia (SFG and TG). SFG strains were closer to TG strains than to STG and Ehrlichia strains. Phylogenetic trees based on 16S rRNA gene sequences were constructed by UPGMA and the neighbor-joining method (data not shown) to investigate the relationships among the rickettsial strains (Fig. 2). Results similar to those found with the phylogenetic tree based on the groEL gene sequences were obtained with the phylogenetic tree based on the 16S rRNA gene sequences. However, R. akari was branched independently and formed a distinct cluster clearly separate from the other SFG strains (Fig. 2), as was shown in a previous report (16).
![]() View larger version (19K): [in a new window] |
FIG. 1. Phylogenetic tree based on groEL gene sequences of the rickettsial strains tested. The phylogenetic tree was constructed by UPGMA with MEGA software. Bootstrap analysis was performed with 100 replicates.
|
![]() View larger version (20K): [in a new window] |
FIG. 2. Phylogenetic tree based on 16S rRNA gene sequences of the rickettsial strains tested. The phylogenetic tree was constructed by UPGMA with MEGA software. Bootstrap analysis was performed with 100 replicates.
|
|
View this table: [in a new window] |
TABLE 4. Levels of similarity between groEL sequences and between GroEL amino acid sequences
|
Phylogenetic analysis with groEL genes from Korean isolates. Phylogenetic trees were constructed by the neighbor-joining method (Fig. 3) and UPGMA (data not shown). In the phylogenetic tree, isolates Youngworl and Hwasung formed a cluster with the STG strains separate from the other rickettsial strains. Isolates 87-91 and 87-100 formed a cluster with R. typhi Wilmington, were closer to R. prowazekii Breinl, and were separate from the other rickettsial strains tested (Fig. 3).
![]() View larger version (19K): [in a new window] |
FIG. 3. Phylogenetic tree based on the groEL gene sequences of Korean isolates. The phylogenetic tree was constructed by the neighbor-joining method with MEGA software. Bootstrap analysis was performed with 100 replicates.
|
|
View this table: [in a new window] |
TABLE 5. DdeI restriction fragments of amplified groEL DNAa
|
![]() View larger version (43K): [in a new window] |
FIG. 4. DdeI restriction patterns of the amplified groEL DNA (534 to 546 bp) from rickettsial strains. Lanes: M, marker DNA (25- and 100-bp mixed DNA ladder); 1, O. tsutsugamushi Karp; 2, O. tsutsugamushi Kato; 3, O. tsutsugamushi Boryong; 4, O. tsutsugamushi Kawasaki; 5, O. tsutsugamushi Youngworl; 6, O. tsutsugamushi Hwasung; 7, R. prowazekii Breinl; 8, R. akari MK; 9, R. typhi Wilmington; 10, R. typhi 87-91; 11, R. typhi 87-100; 12, R. sibirica 246; 13, R. conorii Indian Tick Typhus; 14, R. japonica YH.
|
|
View this table: [in a new window] |
TABLE 6. Oligonucleotide sequences of species-specific primer sets
|
![]() View larger version (102K): [in a new window] |
FIG. 5. Amplification of groEL DNA from rickettsial strains with the primers specific for the species of STG (A), TG (B), and the genus Rickettsia (C). Lanes M, marker DNA (25- and 100-bp mixed DNA ladder). (A) Lanes: 1, O. tsutsugamushi Karp; 2, O. tsutsugamushi Kato; 3, O. tsutsugamushi Boryong; 4, O. tsutsugamushi Kawasaki; 5, O. tsutsugamushi Youngworl; 6, O. tsutsugamushi Hwasung; 7, R. prowazekii Breinl; 8, R. typhi Wilmington; 9, R. sibirica 246; 10, R. conorii Indian Tick Typhus; 11, R. japonica YH; 12, R. akari MK; 13, Coxiella burnetii Henzerling; 14, Ehrlichia sennetsu Miyayama. (B) Lanes: 1, R. prowazekii Breinl; 2, R. typhi Wilmington; 3, R. typhi 87-91; 4, R. typhi 87-100; 5, O. tsutsugamushi Karp; 6, O. tsutsugamushi Kato; 7, O. tsutsugamushi Boryong; 8, O. tsutsugamushi Kawasaki; 9, R. sibirica 246; 10, R. conorii Indian Tick Typhus; 11, R. japonica YH; 12, R. akari MK; 13, C. burnetii Henzerling; 14, E. sennetsu Miyayama. (C) 1, R. sibirica 246; 2, R. conorii Indian Tick Typhus; 3, R. japonica YH; 4, R. akari MK; 5, R. prowazekii Breinl; 6, R. typhi Wilmington; 7, R. typhi 87-91; 8, R. typhi 87-100; 9, O. tsutsugamushi Karp; 10, O. tsutsugamushi Kato; 11, O. tsutsugamushi Boryong; 12, O. tsutsugamushi Kawasaki; 13, C. burnetii Henzerling; 14, E. sennetsu Miyayama.
|
|
|
|---|
GroEL is a highly conserved heat shock chaperonin system of proteins in bacteria, thus establishing a phylogenetic relationship between rickettsiae and bacteria (7, 12, 20). The groEL gene provides better topological and statistical support than the gltA and 16S rRNA genes for the determination of similarity within Bartonella species (12). The groEL gene can also be used as a tool for the rapid identification of Streptococcus (5), Enterococcus (6), Staphylococcus (7), and Mycobacterium (20) by dot blot hybridization, PCR-RFLP analysis, etc.
To date, the groEL sequences of rickettsiae available from GenBank are those of R. rickettsii R (GenBank accession no. U96733), R. prowazekii Breinl (GenBank accession no. Y15783), and O. tsutsugamushi Karp (GenBank accession no. M31887). We designed and synthesized a set of primers whose sequences were based on regions of groEL gene sequences, which were previously shown to be conserved between the rickettsiae and other, different microorganisms, and that could successfully amplify groEL DNA (534 to 546 bp).
The interspecies differences of the groEL gene sequences among the SFG rickettsial strains could not be compared with those of the gltA and ompA gene sequences because the groEL gene sequences of only five SFG strains were compared.
Analysis of the groEL gene has several advantages over the analysis of other genes. Differentiation between members of the genus Rickettsia and the STG rickettsiae by gltA and ompA gene sequence analyses was not possible because the sequences of the gltA and ompA genes of the STG rickettsiae had not been available until recently. Analysis of the groEL genes of the STG rickettsiae could be done with the same primers used for amplification of the groEL genes of members of the genus Rickettsia. In comparison with the 16S rRNA gene sequences, the groEL gene sequences have higher degrees of divergence among the rickettsiae. More than 91.6% similarity was observed among the groEL gene sequences of strains of the genus Rickettsia (SFG and TG). More than 97.1% similarity was observed among the 16S rRNA gene sequences of strains of the genus Rickettsia (SFG and TG). The similarities of the groEL gene sequences among SFG, TG, and STG strains were 95.3 to 99.6, 97.0, and 95.4 to 99.1%, respectively, whereas the similarities of the 16S rRNA gene sequences among SFG, TG, and STG strains were 97.4 to 99.4, 99.4, and 99.2 to 99.7%, respectively. Using the higher level of divergence of the groEL gene sequence, we successfully developed rapid differentiation methods: PCR-RFLP analysis and a species-specific PCR based on the groEL gene sequences. These methods will be especially useful in areas of the Asia-Pacific region where scrub typhus is endemic. These results show that groEL gene analysis is useful for the differentiation of STG strains and strains of the genus Rickettsia.
The sequences of some species (R. akari, R. australis, R. helvetica, R. belii, R. canada, R. typhi, and R. prowazekii) were not available for ompA gene analyses. However, the sequences of R. akari, R. typhi, and R. prowazekii were available for groEL sequence analysis, although other species (R. australis, R. helvetica, R. belii, and R. canada) could not be used for groEL sequence analysis. R. akari clustered with the SFG rickettsiae in the tree based on the groEL sequence, although it was considered an outlying member of SFG (16). R. akari also clustered with the SFG rickettsiae in the tree based on the gltA sequence (16) and was readily distinguished from the STG and TG rickettsiae by the species-specific PCR based on groEL sequences. These results show that groEL sequence analysis is also very useful for the differentiation of the SFG and TG rickettsiae.
We characterized four Korean isolates (isolates Youngworl, Hwasung, 87-91, and 87-100). The four isolates were identified as O. tsutsugamushi and R. typhi only by their reactivities with monoclonal antibodies. In the phylogenetic tree, Youngworl and Hwasung formed a cluster with STG strains separate from the other rickettsial strains (Fig. 3). Isolates 87-91 and 87-100 formed a cluster with R. typhi Wilmington and were closer to R. prowazekii Breinl and were separate from other rickettsial strains (Fig. 3). Youngworl and Hwasung had the same DdeI restriction pattern (407 and 139 bp) as the other O. tsutsugamushi strains and were readily distinguished from the other rickettsial strains tested. Isolates 87-91 and 87-100 had the same DdeI restriction pattern (322, 114, and 98 bp) as R. typhi Wilmington and were readily distinguished from the other rickettsial strains tested.
The results suggest that groEL DNA is useful for the identification and characterization of rickettsiae. groEL gene analysis is able to facilitate the rapid diagnosis of rickettsial diseases and to differentiate rickettsial diseases from other acute febrile diseases.
Scrub typhus is one of the most prevalent febrile illnesses in Korea (19). R. typhi and R. prowazekii have also been isolated in Korea, although there has been no recent documentation of this (1). Although no case of SFG rickettsiosis had been detected in Korea since the isolation of R. akari from Korean voles in 1957 (9), the possibility is strong that the SFG rickettsiae exist in Korea. There are two reasons for this. First, new SFG rickettsiae have been isolated in all parts of the world (15). Since the isolation of R. japonica in Japan, new SFG rickettsiae have been reported in China, Thailand, and the eastern part of Russia (3, 13). Second, serological surveys of wild rodents suggest that SFG rickettsiae exist in Korea (1, 10). groEL gene sequence analysis is able to facilitate the rapid detection and isolation of SFG rickettsiae in Korea.
|
|
|---|
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»