Previous Article | Next Article ![]()
Journal of Clinical Microbiology, August 2004, p. 3686-3695, Vol. 42, No. 8
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.8.3686-3695.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Dominique K. Boudreau,1 Maurice Boissinot,1,2 Ann Huletsky,1,2 Dave Richard,1,2 Marc Ouellette,1,2 Paul H. Roy,1,3 and Michel G. Bergeron1,2*
Centre de Recherche en Infectiologie de l'Université Laval, Centre Hospitalier Universitaire de Québec (Pavillon CHUL), Sainte-Foy, Québec, Canada G1V 4G2,1 Division de Microbiologie, Faculté de Medicine,2 Département de Biochimie et de Microbiologie, Faculté des Sciences et de Génie, Université Laval, Sainte-Foy, Québec, Canada G1K 7P43
Received 6 October 2003/ Returned for modification 15 January 2004/ Accepted 22 April 2004
|
|
|---|
|
|
|---|
Current systems for identification of clinically relevant streptococcal species largely depend on an array of culture-based biochemical tests (12). Some simple and rapid presumptive physiological tests or serological tests are available (12, 15, 39). The most clinically significant pathogens among streptococci can be rapidly identified using phenotypic and immunological tests. However, additional phenotypic tests may be required to confirm identification at the species level. In fact, identification of streptococci to the species level may require up to 7 days because these bacteria grow slowly and because identification may rely on a cumbersome classification system that does not always correlate with phylogenetic analysis (13, 22, 46).
Many DNA-based methods have been applied for the identification and detection of clinically important streptococcal species. Hybridization-based assays for the specific detection of S. pneumoniae (8, 14), S. pyogenes (19, 37), S. agalactiae (3), and S. bovis (47) have been developed. However, these probe-based tests are prone to a lack of sensitivity. Consequently, a number of PCR-based approaches having increased analytical sensitivities and allowing detection of streptococci directly from clinical specimens have been developed (1, 6, 10, 17, 21, 25, 34, 35). However, there is no published study on the development of a Streptococcus-specific PCR assay. The use of genus- or group-specific PCR assays coupled with species-specific internal probes should allow researchers to substantially decrease the number of primers used for bacterial identification, thereby simplifying the development of molecular assays for bacteria (2).
Phylogenetic analyses of streptococci conducted using several conserved genes, including those coding for 16S rRNA, heat shock proteins, glucose pyrophosphorylase, and superoxide dismutase, have been reported (4, 13, 23, 33, 43). All of these phylogenetic studies demonstrated the usefulness of genetic approaches to improve the accuracy of streptococcal species identification.
Genus-specific PCR-based assays, targeting the tuf gene encoding elongation factor Tu, for the specific detection of enterococci (26) and staphylococci (31) have been previously described by our group. Similarly, newly generated tuf streptococcal sequences were used in the present study for PCR detection and extensive phylogenetic analysis of streptococci.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Reference bacterial strains used to validate the Streptococcus-specific assay
|
|
View this table: [in a new window] |
TABLE 2. Validation of the PCR assay with DNA from clinical isolates of streptococci
|
Oligonucleotides. The partial tuf gene sequences obtained in this study as well as those available from public databases were analyzed using GCG Wisconsin software (version 10.3; Accelrys Inc., San Diego, Calif.). PCR primers were analyzed using Oligo primer analysis software (version 5.0; Molecular Biology Insights, Cascade, Col.). Oligonucleotides were synthesized with a model 394 DNA/RNA synthesizer (Applied Biosystems).
Streptococcus-specific PCR. For all bacterial species tested, PCR amplifications using the Streptococcus-specific primers were performed from 1 µl of a genomic DNA preparation at 1 ng/µl which was transferred directly to a 19-µl PCR mixture containing 50 mM KCl, 10 mM Tris-HCl (pH 9.0), 0.1% Triton X-100, 2.5 mM MgCl2, 0.4 µM concentrations of each of the Streptococcus-specific primers (Str1 [5'-GTACAGTTGCTTCAGGACGTATC-3'] and Str2 [5'-ACGTTCGATTTCATCACGTTG-3']), 200 µM (each) deoxynucleoside triphosphate (Amersham Biosciences, Piscataway, N.J.), 3.3 µg of bovine serum albumin (Sigma-Aldrich Canada Ltd., Oakville, Ontario, Canada) per µl, and 0.5 U of Taq DNA polymerase (Promega, Madison, Wis.) combined with the TaqStart antibody (BD Biosciences Clontech, Palo Alto, Calif.). Thermal cycling for PCR amplification and agarose gel analysis of the amplified products were performed as previously described (32). The analytical sensitivity (i.e., the minimal number of genome copies detected per PCR) of the PCR assay was determined using serial twofold dilutions of quantitated genomic DNA purified from bacterial strains representing 28 streptococcal species (Table 1). Strict precautions to prevent carry-over of amplified DNA were used (28). Pre- and post-PCR manipulations were conducted in separate areas. Aerosol-resistant pipette tips were used to handle all reagents and samples. Control reactions to which no DNA was added were routinely performed to verify the absence of DNA carry-over.
Phylogenetic analysis. Multiple sequence alignments were performed using PILEUP (GCG Wisconsin package, version 10.3) and/or CLUSTAL W software (version 1.83) (44) and checked manually with a GCG SeqLab editor to verify the quality of the alignments. The SeqLab editor was also used to identify regions containing gaps, indels, or ambiguities to be excluded for phylogenetic analysis. This edition process yielded a 761-bp tuf sequence suitable for phylogenetic analysis and a 1,260-bp sequence for 16S rDNA analysis. Distance phylogenetic trees were generated using a neighbor-joining or heuristic method with MEGA2 software (version 2.1) (27). Evolutionary distance values were calculated by using Kimura's two-parameter substitution model (18). Bootstrap values were obtained for 1,000 randomly generated trees. This number of replicates was sufficient to produce stable tree topologies. Maximum parsimony analyses were performed using the heuristic method of PAUP software (version 4.0b10; Sinauer Associates Inc., Sunderland, Mass.) with general search parameters (41). A sequence of Enterococcus faecalis V583 was used as an outgroup, because this species is phylogenetically close to streptococci.
Nucleotide sequence accession numbers. GenBank accession numbers for the 761-bp tuf sequences determined in this study are as follows: AY266992 for S. acidominimus; AY266993, AY266994, AY266995, AY266996, and AF276256 for S. agalactiae; AF276257 for S. anginosus; AY266997 and AF276258 for S. bovis; AF276259 for S. constellatus subsp. constellatus; AF276260 for S. criceti; AF276261 for S. cristatus; AF276262 for S. downei; AF276263, AY582541, and AY582542 for S. dysgalactiae; AF276264 for S. equi subsp. equi; AF276265 for S. ferus; AF276266 and AY267005 for S. gordonii; AF276267 for S. intermedius; AF276268 for S. macacae; AF276269 and AY582543 for S. mitis; AF274741 for S. mutans; AF276270 and AY582544 for S. oralis; AF276271 for S. parasanguinis; AY267004 for S. parauberis; AF274742, AY267000, AY267001, AY267002, and AY267003 for S. pneumoniae; AF274743 for S. pyogenes; AF276272 for S. ratti; AF276273 for S. salivarius; AF276274 for S. sanguinis; AF276275 for S. sobrinus; AF274744 for S. suis; AF276276 for S. uberis; and AF276277 for S. vestibularis. GenBank accession numbers for the 1,260-bp 16S rDNA sequences are as follows: AY584476 for S. cristatus; AY584477 for S. parauberis; AY584478 for S. dysgalactiae; and AY584479 for S. ferus.
|
|
|---|
|
View this table: [in a new window] |
TABLE 3. DNA sequence identities between the streptococcal tuf sequencesa
|
A multiple alignment of the tuf sequences from streptococcal species as well as those from staphylococci, enterococci, and lactococci revealed regions conserved among streptococci but distinct from those of other bacteria. The regions for which the sequence mismatches for the nonstreptococcal bacterial species were mainly located at the 3' end were chosen as targets for the Streptococcus-specific PCR primers (Str1 and Str2). This strategy allowed effective discriminatory PCR, since mismatches at the 3' end of primers are the most detrimental to PCR. The selected Streptococcus-specific primers generated 197-bp amplicons. Direct sequencing of these amplicons provided sequence information for the 153-bp sequence between the two Streptococcus-specific primers. Comparisons of these 153-bp sequences obtained for the 28 streptococcal species tested revealed sequence identities ranging from 79.7 to 100% (Table 3). For these sequences, there were two pairs of streptococcal species showing more than 98% identity (S. mitis versus S. pneumoniae [98.7%] and S. pyogenes versus S. dysgalactiae [100%]).
The tuf sequences either determined in our laboratory or available in public databases for different strains of the same streptococcal species allowed the analysis of intraspecies sequence variations. As shown in Fig. 1, the species for which two or more strains have been sequenced are S. pneumoniae (n = 8), S. agalactiae (n = 7), S. pyogenes (n = 6), S. gordonii (n = 3), S. dysgalactiae (n = 3), S. mitis (n = 2), S. oralis (n = 2), S. bovis (n = 2), S. mutans (n = 2), and S. uberis (n = 2). We compared the region corresponding to the 153-bp sequence for all of these available sequences and found that the level of intraspecies sequence variation ranged from 0 to 2.6%, depending on the species. More specifically, intraspecies variations were (i) 0% for S. agalactiae, S. dysgalactiae, S. mutans, and S. uberis, (ii) 0 to 0.7% for S. pneumoniae, (iii) 0.7% for S. mitis, (iv) 0 to 1.3% for S. pyogenes, (v) 0.7 to 2.6% for S. gordonii, (vi) 2% for S. oralis, and (vii) 2.6% for S. bovis.
![]() ![]() View larger version (83K): [in a new window] |
FIG. 1. Phylogenetic relationships among 28 streptococcal species. (A) Phylogenetic tree based on a 761-bp portion of tuf. (B) Phylogenetic tree based on a 1,260-bp portion of 16S rDNA. The trees were generated using the MEGA2 heuristic method, and evolutionary distance values were calculated by Kimura's two-parameter substitution model. The value on each branch represents the percentage of bootstrap replications supporting the branch. A total of 1,000 bootstrap replications were calculated. Bootstrap values lower than 50% are not shown. GenBank accession numbers are given in parentheses. The tuf and 16S rDNA portions correspond to nucleotide positions 340 to 1,100 of the complete tuf gene of S. pneumoniae R6 (AE008504) and 93 to 1,382 of the complete 16S rRNA gene of S. pneumoniae R6 (AE008546). All sequences used for these phylogenetic analysis were obtained either from this study or from the following sources: GenBank (http://www.ncbi.nlm.nih.gov) (A), TIGR ongoing genome projects (http://www.tigr.org) (B), Sanger ongoing genome projects (http://www.sanger.ac.uk) (C), and our group (for previously determined sequences) (24) (D); sequences from these four sources are indicated with A, B, C, and D, respectively.
|
Phylogenetic analysis. The entire Streptococcus group is monophyletic as determined on the basis of tuf and 16S rDNA phylogenies and forms a phylum distinct from L. lactis (Fig. 1). tuf-based phylogenetic relationships between the 28 streptococcal species revealed clusters that are generally in agreement with those observed with 16S rDNA phylogeny (Fig. 1 and Table 4). Incidentally, the 16S rDNA tree, which we constructed mostly using available database sequences for the species selected for the present study, was similar in terms of branching to the 16S rDNA tree recently reported by Facklam (13) for 55 streptococcal species. Comparison of tuf and 16S rDNA phylogenetic trees revealed a different ancestor for the pyogenes group (Fig. 1). As determined on the basis of tuf phylogeny, that streptococcal group is linked with the mitis and salivarius groups whereas 16S rDNA results show that the pyogenes group branches separately from the other streptococcal groups. For the other major deep branches, the bootstrap values were too low with both tuf- and 16S rDNA-based phylogenetic trees to permit reliable interpretations (Fig. 1). Phylogenetic analysis of available database sequences for sodA, groEL, atpD (data not shown), and rnpB (42) showed a similar lack of resolution for basal branches of the streptococcal tree. However, branches between closely related taxa are well supported and can be used to establish their grouping (Fig. 1 and Table 4).
|
View this table: [in a new window] |
TABLE 4. Phylogenetic clusters observed with tuf and 16S rRNA genes within phenotypic streptococcal species groups
|
Distance analysis using neighbor-joining or the MEGA2 heuristic method produced similar trees except for minor differences present in the topology of the basal branches for which bootstrap support is weak. Distance and parsimony phylogenetic analysis of tuf yielded trees showing similar end-branching structures (data not shown). However, many basal branches are still poorly resolved.
|
|
|---|
The use of genus- or group-specific PCR assays can substantially decrease the number of primers used for bacterial identification. Indeed, identification of the most frequently encountered species at the genus level is often sufficient to permit selection of an appropriately targeted antibiotic (2). Furthermore, the use of post-PCR hybridization with species-specific internal probes bound onto a solid support (e.g., oligonucleotide arrays) would allow identification at the species level. Fluorescent probes (e.g., TaqMan probes) may also be used for the development of species-specific real-time PCR assays for clinically important streptococci. We chose the tuf gene as a genetic target to develop a PCR-based assay for the detection of streptococci, because it has both conserved and variable regions suitable for the design of genus- and species-specific probes, respectively (26, 31). Our assay efficiently detected all 28 streptococcal species tested, with an analytical sensitivity of 1 to 10 genome copies per PCR. It also detected DNA purified from two phylogenetically closely related species (i.e., E. durans and L. lactis). However, the 15-to-31% nucleotide sequence divergence in the tuf gene of these two species compared to corresponding sequences for the 28 streptococcal species provides much flexibility for the development of internal probes specific to streptococci.
Based on 16S rDNA sequence analysis, the genus Lactococcus is phylogenetically very closely related to the genus Streptococcus (40). It is therefore not surprising that L. lactis DNA was detected by the Streptococcus-specific PCR assay. Indeed, analysis of the primer binding sites for this species revealed that there was no mismatch in the Str1 primer and only one at the ninth nucleotide of the 3' end for the Str2 primer. Phylogenetic analysis of tuf sequences from different lactococcal species confirmed that they form a distinct phylum closely related to streptococci (data not shown).
Analysis of tuf sequences from a variety of streptococci and lactococci revealed that only one copy of tuf is present in their genome. The nonspecific amplification of E. durans by the Streptococcus-specific PCR assay can be explained by the finding that E. durans and some closely related enterococcal species carry two divergent tuf genes, one of which shares a common ancestor with the tuf gene of streptococci and lactococci (24). Indeed, it is the horizontally transferred tufB gene of E. durans that was amplified by the Streptococcus-specific PCR-based assay. However, no other enterococcal tufB genes were amplified by this assay.
The 153 clinical isolates representing a variety of streptococci which were obtained from the Quebec City region were all detected by our assay except for one S. bovis strain which was not amplified. The two other clinical strains of S. bovis as well as the reference S. bovis strain ATCC 33317 were amplified efficiently. The nondetectable S. bovis strain was of biotype II. S. bovis of this biotype has been shown to be more genetically polymorphic and frequently associated with animal hosts (13, 43, 45). Sequence data for this S. bovis strain revealed the presence of a single mismatch at the 3' end (first nucleotide) of the Streptococcus-specific primer Str1 which probably explains the absence of PCR amplification.
S. oralis and S. mitis are phylogenetically very closely related to S. pneumoniae (22). These three species form well-supported distinct groups in tuf phylogeny but not with 16S rDNA phylogeny, where S. mitis is present in two different branches (Fig. 1). The findings with 16S rDNA phylogeny regarding S. mitis are consistent with previous studies showing that the genes coding for the streptococcal pneumolysin (ply), autolysin (lytA), and superoxide dismutase (sodA) are more polymorphic in S. mitis (23, 46). The genetic heterogeneity of S. mitis may be associated with the horizontal transfer of genes between streptococci (9, 20, 29). The present study also revealed a very low intraspecies tuf sequence divergence in S. pneumoniae (0 to 0.3% for 8 strains). The homogeneity of tuf sequences observed in S. pneumoniae could be explained by the clonal spread of a limited number of variants of this bacterial pathogen (9, 11).
Sequence analysis of the 761-bp portion of tuf revealed that sequence variations between S. oralis, S. mitis, and S. pneumoniae ranged from 1.3 to 2.9% (Table 3) compared to 0.3 to 0.6% for the 1,260-bp portion of 16S rDNA (data not shown). Hence, the higher level of sequence variations in tuf compared to 16S rDNA for these very closely related species provides more potential for the development of primers and probes allowing them to be distinguished. These three species can also be distinguished by using less-conserved genes like ddl (coding for D-alanine-D-alanine ligase), lytA, and sodA (16, 17, 23, 38).
According to our tuf sequence data for the 153-bp amplicon sequences, species-specific sequence variations are present for all streptococcal species except for S. pyogenes and S. dysgalactiae, for which the tuf amplicon sequences are identical. However, other regions of tuf may be appropriate to distinguish these two species, as suggested by the 1.4% divergence in their sequences for the 761-bp tuf portion. Surprisingly, the 16S rDNA sequences for S. pyogenes and S. dysgalactiae revealed a significantly higher level of divergence (i.e., 3.2%). Consequently, S. dysgalactiae branched together with S. agalactiae according to 16S rDNA phylogeny while it appeared to be more related to S. pyogenes on the basis of tuf phylogeny (Table 4 and Fig. 1). Phylogenetic studies performed with sodA, groEL, atpD (data not shown), and rnpB (42) support this relationship between S. dysgalactiae and S. pyogenes.
The comparison of tuf and 16S rDNA phylogenetic trees revealed differences in branching topologies and clustering for some streptococcal species. For S. criceti, S. cristatus, S. mutans, S. ratti, S. sanguinis, and S. uberis, the clustering differences for these species could be explained by the lower resolution of the 16S rDNA phylogenetic tree. In the case of S. dysgalactiae, S ferus, and S. suis, apparent discrepancies between the phylogenetic trees are not statistically well supported. It is possible that these differences are associated with different evolutionary rates for tuf and 16S rDNA. It can also be linked to the heterogeneity of 16S rDNA operons in bacteria. Indeed, it has been reported that phylogenetic studies are severely limited by 16S rDNA heterogeneity and that analysis of distinct rRNA operons within the same microbial strain may lead to different results (5, 7, 36). Finally, for S. suis and S. mitis, differential branching could be the result of previously reported greater genetic variation within these two species (4, 13, 23, 43, 46).
Phylogenetic analysis using either tuf or 16S rDNA sequences revealed that three clinically important streptococcal phenotypic species groups (i.e., mitis, pyogenes, and mutans) are composed of three to five phylogenetic clusters and that the groups anginosus and salivarius were monophyletic. This suggests that the mitis, pyogenes, and mutans groups are much more genetically heterogeneous than the anginosus and salivarius groups. It also indicates that the standard streptococcal species grouping, which is based on physiological and biochemical characteristics determined by conventional methods, does not correlate well with the level of genetic diversity within each group.
Phylogenetic analysis of multiple strains of the same streptococcal species revealed that the level of intraspecies tuf sequence variations depends on the streptococcal species. We analyzed the 761-bp tuf sequences from a total of six to eight strains for each of the three most clinically important streptococcal species (i.e., S. pyogenes, S. pneumoniae, and S. agalactiae). This analysis clearly demonstrated that S. pyogenes has more polymorphisms (intraspecies variation of 0.1 to 3.5%) than S. pneumoniae and S. agalactiae (intraspecies variation of 0 to 0.3% and 0 to 0.4%, respectively). By contrast, analysis of 16S rDNA sequences suggests that S. pyogenes has slightly fewer intraspecies sequence variations than S. agalactiae or S. pneumoniae (0 to 0.2% versus 0 to 0.5%). These observations indicate that evolutionary rates for these two genes involved in different components of the protein synthesis machinery may differ for streptococcal species.
In conclusion, we have performed an extensive sequence analysis of streptococci showing that tuf generally offers a better discrimination power than 16S rDNA to distinguish streptococcal species. tuf and 16S rDNA phylogenetic trees were generally in agreement, although different clustering of some closely related streptococcal species was observed. However, these phylogenetic clusters revealed that classical streptococcal phenotypic groups may comprise different genetic subgroups. We have used tuf sequences to develop a PCR-based approach for the detection of streptococci. Future developments will seek to combine this genus-specific assay with detection of species-specific tuf sequence polymorphisms by using internal hybridization probes to provide a molecular diagnostic tool for rapid and accurate diagnosis of streptococcal infections.
We thank Louise Côté, director of the microbiology laboratory of CHUQ (Pavillon CHUL) for access to the laboratory and for providing clinical isolates of streptococci. We also thank Gisèle Chassé and Ève Bérubé for their contribution in microbial species culture and identification. Sequence data for S. agalactiae A909, S. gordonii NCTC 7868, and S. pneumoniae 670 were obtained from The Institute for Genomic Research website at http://www.tigr.org. Sequencing of S. agalactiae A909 was accomplished with support from the National Institute of Allergy and Infectious Diseases (NIAID). Sequencing of S. pneumoniae 670 was accomplished with support of the National Institute of Allergy and Infectious Diseases and the University of Alabama. Sequencing of S. gordonii NCTC 7868 was accomplished with support from National Institute of Dental and Craniofacial Research. Sequence data for S. pyogenes Manfredo and S. uberis 0140J were produced by the sequencing groups at the Sanger Institute and can be obtained from ftp://ftp.sanger.ac.uk/pub/pathogens/sp/ and ftp://ftp.sanger.ac.uk/pub/pathogens/su/.
Present address: Bioniche Therapeutics Research Center, Montreal, Quebec, Canada H4P 2R2. ![]()
|
|
|---|
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»