Previous Article | Next Article ![]()
Journal of Clinical Microbiology, April 2009, p. 932-939, Vol. 47, No. 4
0095-1137/09/$08.00+0 doi:10.1128/JCM.01449-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Department of Bacteriology, Mycology and Parasitology, Statens Serum Institut, Copenhagen,1 Department of Clinical Microbiology, Rigshospitalet, Copenhagen,2 Department of Clinical Microbiology, Odense University Hospital, Odense, Denmark3
Received 29 July 2008/ Returned for modification 30 July 2008/ Accepted 23 January 2009
|
|
|---|
|
|
|---|
There are reports describing associations between some specific species and clinical manifestations. Streptococcus gallolyticus and Streptococcus lutetiensis in the bovis group are reported to have a strong association with colon cancer (31-33). Streptococcus anginosus and Streptococcus intermedius in the anginosus group are associated with abscess formation (9, 26), and Streptococcus mitis and other normal oral commensals have a strong association with infective endocarditis (6, 10, 22). The precise species-level identification of NHS from relevant clinical specimens is crucial to making the right diagnosis and understanding the pathogenesis of the infection.
The conventional phenotypic tests do not always allow accurate identification. The automated systems, such as the Rapid ID 32 Strep and the Vitek 2 GP systems (BioMerieux, France), which are based on phenotypic tests, are widely applied in the clinical microbiology laboratory. The large number of species relative to the number of biochemical traits examined, the variability of several traits within species, the poor reproducibility of some tests, and the lack of sufficient phenotypic data for more recently described species in the underlying databases often result in shortcomings with regard to the exact species designations (14, 23, 29, 41). In a study by Hoshino et al. (23) with 148 strains consisting of 115 clinical isolates and 33 reference strains, the rate of correct identification by commercial kits was below 50% but varied significantly between species. The most significant problems were observed with S. mitis, S. oralis, and the 11 Streptococcus species that have been described since 1991. These inherent problems call for alternative means of identification.
Early and effective antimicrobial treatment can result in negative cultures with important clinical specimens, e.g., heart valve tissue or brain abscess material. This stresses the need for the possibility of performing non-culture-based molecular biology examinations (43).
Gene sequence analysis has been applied in an attempt to make an accurate species-level identification of streptococci. The target sequences have included genes encoding functional RNA (the 16S rRNA gene [3], rnpB [45]), protein-coding genes (sodA [23], tuf [37], groESL [42], rpoB [11], gdh and ddl [12]), and noncoding spacer regions like the intergenic spacer (ITS) region (7). However, it seems that the analysis of no single gene is sufficient for this very heterogeneous group (23, 40).
Alternatively, the performance of multilocus sequence analysis, as suggested by some studies (23), is not always realistic in the local clinical microbiology laboratory or even in reference laboratories.
The ITS region has been reported to be used for the strain typing of staphylococci and Clostridium difficile (15, 16). It is also suggested to be a good candidate for species identification (17, 20, 36). According to Chen et al., the intraspecies similarity scores of 11 viridans group streptococcal species were high and ranged from 0.97 to 1.0, and pairwise comparison of two species of viridans group streptococci revealed a lower level of sequence similarity between their ITS regions than between their 16S rRNA gene sequences (7). These characteristics indicate that the ITS region might constitute a more discriminative target sequence than the 16S rRNA gene for the differentiation of closely related species of NHS.
The gdh gene, which is about 1,500 bp in length, encodes a 45-kDa glutamate-6-phosphate dehydrogenase. It is used together with six other housekeeping genes for multilocus sequence typing of pneumococci (23). The gdh gene is reportedly highly conserved, as it exhibits an extremely small number of point mutations relative to the numbers in other genes (30). It has been shown that partial gdh sequences could be used to unambiguously differentiate S. pneumoniae from S. mitis and S. oralis (25, 30).
The purpose of this study was to establish a method based on ITS and partial gdh sequence analysis that is capable of unambiguously identifying clinically relevant NHS to the species level. At the same time, the method should be able to be easily applied in a reference laboratory.
|
|
|---|
-galactosidase, alkaline phosphatase, arginine decarboxylase, urease, extracellular polysaccharide (dextran and levan), and esculetin from esculin; performance of the peroxide test; and the detection of acid production from inulin, salicin, raffinose, amygdalin, and glycogen. The reference strains from SSI, ATCC, NCTC, and CCUG were also identified by 16S rRNA gene sequence analysis at SSI (8). Most of the SK strains were also well identified by phylogenetic analysis of the nucleotide sequences of four housekeeping genes, ddl, sodA, gdh, and rpoB, at the Institute of Medical Microbiology and Immunology, University of Aarhus (23). |
View this table: [in a new window] |
TABLE 1. The 68 strains of the genus Streptococcus used in the study
|
PCR primers. To amplify the ITS region, we designed a forward primer, primer Strep16S-1471F (5'-GTG GGA TAG ATG ATT GGG GTG AAG T-3'), the 5' end of which is located at position 1471 of the 16S rRNA gene (Escherichia coli numbering). Reverse primer 6R-IGS (5'-GGG TTC CCC CAT TCG GAH AT-3') was adapted and improved from the reverse primer of Chen et al. (7). The 5' end of primer 6R-IGS is located at position 108 downstream of the 5' end of the 23S rRNA gene (E. coli numbering).
To amplify the partial gdh gene, we used two primers, primer Strep-gdhF (5'-ATGGACAAACCAGCNAGYTT-3') and primer Strep-gdhR (5'-GCT TGA GGT CCC ATR CTN CC-3'), which amplify a 660-bp amplicon.
PCR analysis of ITS region and partial gdh gene. PCR of the ITS region and subsequent sequence analysis were performed with all 68 streptococcal strains. PCR of the partial gdh gene and subsequent sequence analysis were performed only with the 39 strains belonging to the mitis group (11 S. oralis strains, 11 S. mitis strains, and 17 S. pneumoniae strains).
The PCR was performed with 50-µl reaction volumes consisting of 1x HotStarTaq master mix (containing final concentrations of 2.5 U HotStarTaq DNA polymerase, 1x PCR buffer, 200 µM each deoxynucleoside triphosphate; Qiagen) and 0.4 µM (final concentration) each primer. The conditions of the PCR with the primers for the ITS region (primers Strep16S-1471F and 6R-IGS) were as follows: 95°C for 15 min and 40 cycles of 94°C for 30 s, 61°C for 30 s, and 72°C for 30 s. The conditions of the PCR with the primers for the partial gdh region were as follows: 95°C for 15 min and 35 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 30 s. Both PCRs were performed with a DNA Engine Dyad cycler (Bio-Rad). The PCR products were analyzed by 2% gel electrophoresis after staining of the gels with ethidium bromide. The PCR products were purified with spin columns (Microcon YM-100 filter units; Millipore, Billerica, MA).
Sequencing of PCR products. (i) sequencing of ITS region and sequence editing. Both DNA strands of the amplicons were sequenced on an ABI Prism 3100 Avant genetic analyzer (Applied Biosystems, Foster City, CA) with a BigDye (version 3.1) kit (Applied Biosystems). Primers Strep16S-1471F and 6R-IGS were used as sequencing primers.
After sequencing of the PCR product, the sequence had to be edited, as described by Chen et al. (7). Portions of the 16S and 23S rRNA gene regions were removed to obtain the full sequences of the ITS region with CTAAGG at the 5' ends and AATAA at the 3' ends of the sequences of the ITS region.
(ii) Sequencing of partial gdh sequence and sequence editing. Primers Strep-gdhF and Strep-gdhR were used for the sequencing PCR. The sequences were manually edited as described by M. Kilian (Institute of Medical Microbiology and Immunology, University of Aarhus) so that the sequences were more comparable in the NCBI BLAST search engine. For all the strains evaluated in this study, the 5' ends of the gdh sequences were TTTAAAAACCT, whereas the 3' ends of the gdh sequences were cut just before the sequence AGA ACC ATA C, so that the 3' end of the edited sequence was TGC TTC/A TCC.
The edited sequences of the ITS region and the partial gdh gene were compared to sequences deposited in the NCBI database by using the BLAST search engine and by taking into consideration the percentage and number of identities, the maximum score, and E values for the best and the next best taxon matches.
Phylogenetic analysis. DNA sequences were aligned by using the ClustalW program built into the MEGA (version 4.0) program package. Phylogenetic analysis on the basis of the sequences of the ITS region and the partial gdh gene for the 39 strains belonging to the mitis group were performed by the neighbor-joining and minimal evolution methods in the MEGA (version 4.0) program package (downloaded from http://www.megasoftware.net). The distance between the sequences was calculated by using the Kimura two-parameter model.
Nucleotide sequence accession numbers. The GenBank accession numbers of the sequences of the ITS region and the partial gdh gene obtained from our study are listed in Table 1.
|
|
|---|
PCR of the partial gdh gene yielded a single band at 660 bp, with the size being constant independent of the species.
Identification of Streptococcus strains (HS and NHS) on the basis of sequence analysis of ITS region. Sequence analysis of the ITS region was performed with all 68 streptococcal strains. The sizes of the edited ITS region sequences varied from 248 to 498 bp. All species except for the S. mitis and S. pneumoniae strains achieved the expected taxon as the best taxon match. S. anginosus had only a 96% sequence identity. Low maximum score differences from the second best taxon match (difference, <10) were obtained for the S. salivarius and S. vestibularis strains, although the percent sequence identities were high (99 to 100%), and the next best taxon matches belonged to the salivarius group. Among the 39 strains in the mitis group, all 11 S. oralis strains achieved the expected taxon as the best taxon match, although maximum score differences to the next best taxon match (S. mitis or S. pneumoniae) were relatively low (differences, 7 to 17). Four of the 11 strains of S. mitis (strains SK612, SK614, SK648, and SK661) had S. pneumoniae as the best taxon match and S. mitis as the second best taxon match. However, the maximum score differences between the two taxons were very low (differences, 5 to 6). Six of the 11 S. mitis strains and the 17 S. pneumoniae strains had the expected taxon as the best taxon match, and the next best taxon matches were S. pneumoniae and S. mitis, respectively. The maximum score differences between the two best taxon matches were low (differences, 0 to 11) (Table 2).
|
View this table: [in a new window] |
TABLE 2. Identification of HS and NHS strains on the basis of ITS sequence dataa
|
|
View this table: [in a new window] |
TABLE 3. Identifications of the 39 strains of S. mitis, S. pneumoniae, and S. oralis on the basis of the partial gdh sequencesa
|
Phylogenetic analysis on the basis of the sequences of the ITS region and partial gdh gene. The phylogenetic relationship of S. mitis, S. oralis, and S. pneumoniae derived from the sequences of the ITS region is presented in Fig. 1. The 39 strains failed to form any distinct clusters according to the taxon. The evolutionary distances for 38 strains were within 0.002. The distance between the last S. oralis strain (strain SK152) to the other strains was only 0.005. Therefore, sequence analysis of the ITS region was insufficient for discrimination of the three species S. oralis, S. pneumoniae, and S. mitis from each other.
![]() View larger version (13K): [in a new window] |
FIG. 1. Phylogenetic tree determined on the basis of the sequences of the ITS regions of 11. S. mitis (Smit) strains, 11 S. oralis (Soral) strains, and 17 S. pneumoniae (Spneu) strains obtained by the unrooted neighbor-joining method in the MEGA (version 4.0) program package. The scale bar indicates the evolutionary distance between the sequences determined by calculation of the percent sequence divergence. It clearly demonstrates that S. mitis, S. pneumoniae, and S. oralis are genetically closely related species and cannot be discriminated from each other on the basis of the sequences of their ITS regions.
|
![]() View larger version (4K): [in a new window] |
FIG. 2. Minimal evolution algorithm (suppressed) obtained by using the MEGA (version 4.0) program and based on the partial gdh gene sequences of 11 S. oralis strains, 17 S. pneumoniae strains, and 13 S. mitis strains. It shows that the three species form three distinct clusters. The S. oralis cluster has a longer distance to the two other clusters, indicating that S. pneumoniae and S. mitis are genetically more closely related on the basis of gdh gene evolution. There are three subclusters within the S. mitis cluster, indicating that the species S. mitis contains a heterogeneous group of strains.
|
|
|
|---|
In our study, sequence analysis of the ITS region was used as a first-line tool for the identification of species in the genus Streptococcus. All 11 strains (100%) belonging to the five HS species achieved the correct species as the best taxon match on the basis of sequence analysis of the ITS region. The maximum score differences from the next best taxon match varied from 18 to 422. These differences are large enough for the differentiation of species. Of the strains belonging to the 19 NHS species, 53/57 strains (94%) achieved the correct species designation as the best taxon match.
Members of the mutans, sanguinis, and anginosus groups achieved unambiguous identifications with high identity scores, and the differences in the maximum scores from the best to the next best taxon match were significant (232 to 543, 127 to 397, and 38 to 394, respectively).
A high degree of heterogeneity within the S. anginosus species has been reported previously (27). It was suggested that the species contains several subspecies or new species. In our study, we included only one strain of S. anginosus. Therefore, the low percent identity (96%) probably reflects the heterogeneity in this species. Further molecular taxonomic studies are needed to explore this heterogeneity within strains belonging to this species.
The type strain of S. gallolyticus subsp. gallolyticus achieved the expected taxon as the best taxon match, although the maximum score difference from the next best taxon match, S. macedonicus, was only 6. The taxonomy of the bovis group has undergone dynamic changes in the last two decades. Recent studies based on DNA-DNA hybridization, 16S rRNA gene sequencing, and sodA gene sequencing revealed that the S. bovis/S. equinus complex consists of five clusters. S. gallolyticus and S. macedonicus belong to one cluster. The DNA-DNA hybridization data from the same study revealed that the genomes of S. macedonicus and S. gallolyticus display >70% homology, which supports the hypothesis that S. macedonicus and S. gallolyticus are a single species (38). The four other strains belonging to the bovis group included in this study, S. infantarius subsp. infantarius, S. lutetiensis, and S. equinus, achieved unambiguous species identification on the basis of sequence analysis of the ITS region. Strain ATCC 15351 was previously named S. bovis. The high degrees of similarity by both DNA-DNA hybridization and 16S rRNA gene sequencing brought the conclusion that the species S. equinus and S. bovis belong to a single species. The name S. equinus has nomenclatural priority. Therefore, S. bovis is no longer a recognized taxon (39).
Among the members of the salivarius group, only S. salivarius is commonly identified from a variety of human infections (4, 13, 21). S. thermophilus has been isolated only from dairy products. S. vestibularis was identified from the human oral cavity, and its association with human infections has not been confirmed. In our study, both S. salivarius and S. vestibularis achieved the correct taxon as the best taxon match, although the maximum score distance to the next taxon was only 8. The second best taxon for S. salivarius was S. vestibularis and vice versa. The failure of S. salivarius to produce extracellular polysaccharides on sucrose-containing agar is helpful in securing a correct distinction between these two species.
The differences in maximum scores between the best and the next best taxon matches were very small for the strains belonging to the mitis group (range of differences, 0 to 17), often making it impossible to allocate the strain examined to a specific species. Four strains of S. mitis had S. pneumoniae as the best taxon match, followed very closely by S. mitis (maximum score differences, 5 to 6). The phylogenetic analysis of the results of sequence analysis of the ITS region revealed that S. oralis did not form a distinct cluster in relation to S. mitis and S. pneumoniae. It is well known that members of the mitis group are also closely related on the basis of their 16S rRNA gene sequences (18, 28).
On the basis of these results, we concluded that the ITS region can be used for the species-level identification of strains belonging to the hemolytic, anginosus, mutans, and sanguinis groups. Strains belonging to the salivarius group and the mitis group can be identified only to the group level.
Sequence analysis of the partial gdh gene proved to be useful in separating S. mitis, S. oralis, and S. pneumoniae, which are otherwise hard to differentiate from each other. In this study, 17/17 S. pneumoniae strains (100%), 11/11 S. oralis strains (100%), and 10/11 S. mitis strain (90.9%) were correctly identified. Only one S. mitis strain could not be discriminated from S. pseudopneumoniae. The maximum score difference between these two taxons was 0. This is probably because of the taxonomic changes that have been made in recent years. S. pseudopneumoniae is a relatively new taxon that was first described in 2004 and that is closely related to S. mitis and S. pneumoniae (1). Phylogenetic analysis based on the sodA gene sequence showed that the strains assigned to the species S. pseudopneumoniae were more closely aligned with S. mitis than with S. pneumoniae (1). The phylogenetic analysis based on concatenated partial sequences of the ddl, gdh, rpoB, and sodA genes showed that S. pseudopneumoniae is included within the pneumoniae-mitis-pseudopneumoniae cluster (29). This species was, however, not included in our study.
The three distinct clusters of S. mitis, S. pneumoniae, and S. oralis in the phylogenetic algorithm based on partial gdh sequences proved that gdh sequence analysis is capable of discriminating these three genetically closely related species from each other. As a result, we have a tool for the species-level identification of these species.
The phylogenetic algorithm also gave us information about the genetic relationship between the three species. There was a much longer distance from the S. oralis cluster to the S. mitis and S. pneumoniae clusters than between the S. mitis and S. pneumoniae clusters. This suggests that S. mitis and S. pneumoniae are genetically more closely related to each other than they are to S. oralis. The S. mitis cluster was much more scattered than the clusters of the other two species, and it formed several subclusters. This suggests that S. mitis is genetically more heterogeneous. This is in accordance with the recent observations of Bek-Thomsen et al., who observed that the range of interstrain gdh sequence distances was significantly larger for S. mitis than for what was found among S. pneumoniae strains (2).
In this study, we present a reliable method for the identification of clinically relevant streptococci, with a focus on NHS, to the species level. The method is easy to perform in a laboratory that has sequencing facilities. By sequence analysis of the ITS region, all HS strains and most NHS strains could be identified to the species level. S. mitis, S. pneumoniae, and S. oralis could not be unambiguously discriminated from each other by sequence analysis of the ITS region. A second sequence analysis based on the partial gdh gene distinguished these three species from each other. Only one S. mitis strain could not be unambiguously discriminated from S. pseudopneumoniae, which is probably because of the new nomenclature change. The phylogenetic tree based on the gdh gene sequences clearly shows that S. oralis, S. mitis, and S. pneumoniae form three distinct clusters. If colonies are available, sequencing of the ITS region and the partial gdh gene can both be completed within 24 to 72 h.
On the basis of the results from this study, we conclude that the combination of sequence analysis of the ITS region and sequence analysis of the partial gdh gene is a potential tool for the identification of clinically relevant streptococci in a clinical microbiology reference laboratory.
Published ahead of print on 4 February 2009. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»