Previous Article | Next Article ![]()
Journal of Clinical Microbiology, January 2004, p. 250-256, Vol. 42, No. 1
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.1.250-256.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Infectious Diseases Epidemiology Research Unit, Division of Infectious Diseases, School of Medicine,1 Center for Biomedical Informatics,3 Intelligent Systems Program,4 Department of Infectious Diseases and Microbiology, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania,2 Department of International Health, Bloomberg School of Hygiene and Public Health, Johns Hopkins University, Baltimore, Maryland5
Received 6 April 2003/ Returned for modification 19 July 2003/ Accepted 22 September 2003
|
|
|---|
1% (0 to 5 of 476 bp) from each other and that were
16% (78 to 83 of 476 bp) divergent from that of the 9V-3 clone. Allowing for a 2-bp difference in the cpsB sequence resulted in the highest correlation between the cpsB gene and serotype. Overall, 95% (84 of 88) of the strains were classified correctly by serotype with the cpsB sequence. The distal recombination site of the Baltimore serotype 14 variants of the 9V-3 clone was not identical to that of the European serotype 14 variants. The cpsB gene was serotype specific regardless of whether capsular switching occurred. Although the correlation between serotype and the cpsB sequence was high, the overall diversity of the cpsB gene within a serotype likely will limit the role of this gene in a sequence-based serotyping method. |
|
|---|
The France 9V-3 pneumococcal clone is 1 of approximately 20 international clones and has been detected as serotype 14, serotype 9A, and serogroup 19 variants (18). These international clones account for a significant proportion of non-penicillin-susceptible Streptococcus pneumoniae (PNSP) strains in the United States (26). The proximal recombination site for the serotype 14 variants of the 9V-3 clone from Uruguay, Denmark, and Spain, but not Poland, was detected within the cpsA gene (2). The distal site of recombination occurred within the pbp1a gene, which has an orientation divergent from that of the capsular genes. The serotype 14 variants had a pbp1a sequence that was a combination of that of the non-penicillin-susceptible 9V-3 clone and that of the nonencapsulated penicillin-susceptible reference strain R6. The pbp1a gene sequence of the variants was identical to that of the 9V-3 clone until position 1854 or 1922. Thereafter, the gene sequence diverged from that of clone 9V-3 but was identical to that of strain R6.
Among all of the available non-penicillin-susceptible pneumococcal strains collected over a 2-year period from the Baltimore, Md., metropolitan area, 7.7% (11 of 143) were serotype 14 strains; 27.2% (3 of 11) of the serotype 14 strains were highly genetically related to the France 9V-3 clone rather than to the other serotype 14 strains (17). In this study, we characterized all serotype 14 PNSP strains and compared them to the France 9V-3 clone and to a genetically related serotype 9V strain by using PFGE, MLST, and phylogenetic trees of the concatenated housekeeping gene sequences. We also compared the PCR pbp1a restriction profiles and sequenced internal fragments of the pbp1a and cpsB genes from all serotype 14 strains to determine whether serotype 14 variants of the 9V-3 clone detected in Baltimore were similar to the European variants. Finally, we compared the internal fragments of the cpsB genes from all serotype 9V and 14 strains from Baltimore, 15 pneumococcal clones, and 23 strains from GenBank to assess the correlation between serotype and the cpsB sequence.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. cpsB sequences of 88 strains of selected serotypes
|
80% relatedness on the dendrogram. Due to an error that occurred while making plugs for PFGE, strain VIII was represented twice, while strain XII was not included. The correction of this error creates six, not five, clonal groups (16, 17). The MIC for penicillin-intermediate strains is 0.12 to 1 µg/ml, while that for penicillin-resistant strains is
2 µg/ml. PCR. Pneumococcal strains were incubated overnight on Trypticase soy agar containing 5% sheep blood. Genomic DNA was isolated by using Prepman Ultra in accordance with the manufacturer's instructions (Applied Biosystems, Foster City, Calif.). PCR primers for restriction profiles were pbp1a F (GGCATTCGATTTGATTCGCTTCTATCAT) and pbp1a R (CTGAGAAGATGTCTTCTCAGGCTTTTG) (8). The 30-µl reaction mixture contained 1.5 mM MgCl2, 0.33 µM each primer, 25 µM each deoxyribonucleotide, 1.5 U of the thermostable DNA Taq polymerase mixture, 3 µl of 10x buffer (Invitrogen Corporation, Carlsbad, Calif.) and 20 ng of DNA template. PCR was performed with a model 9700 thermal cycler (Perkin-Elmer, Wellesley, Mass.). The samples were subjected to an initial denaturation at 95°C for 1 min, followed by 10 cycles of denaturation at 94°C for 15 s, annealing at 58°C for 30 s, and elongation at 72°C for 1 min 50 s. This procedure was followed by 20 cycles with the same parameters but with sequential 10-s increments in the elongation cycle. A 7-min extension at 72°C followed the final cycle (8). PCR products were purified with Multiscreen PCR plates (Millipore, Bedford, Mass.). Five microliters of purified PCR products were restricted with 3 U each of MseI and DdeI for 1 to 2 h at 37°C (24). The restriction digests were electrophoresed for 2.5 to 3.0 h at 125 to 150 V through 4% Nusieve 3:1 agarose gels and captured with GelDoc 2000 (Bio-Rad, Hercules, Calif.). Restriction patterns of bands between 300 to 10,000 bp were classified as patterns 1 to 4.
DNA sequencing. (i) MLST.
MLST was performed as described at http://www.mlst.net with internal fragments of the following seven housekeeping genes (protein products are shown in parentheses): aroE (shikimate dehydrogenase), gdh (glucose-6-phosphate dehydrogenase), gki (glucose kinase), recP (transketolase), spi (signal peptidase I), xpt (xanthine phosphoribosyltransferase), and ddl (D-alanine-D-alanine ligase). The primer sets were elongated for all except gki as described previously (7). Identical sequence types (ST), single-locus variants (SLV), and double-locus variants were determined, and strains with
5 of 7 alleles were denoted as a single complex.
(ii) Phylogenetic analysis. Phylogenetic analysis was performed on the nucleotide and amino acid sequences obtained by concatenating the seven genes in order of occurrence in a serotype 4 strain: spi, gki, gdh, aroE, recP, ddl, and xpt (28). The 3,199 nucleotides were translated to 1,064 amino acids with Sequencher. The sequences were edited with Sequencher and aligned with ClustalX. Neighbor-joining (NJ), maximum-likelihood (ML), and maximum-parsimony (MP) trees were created with Paup 4.0, version beta10 (Sinauer Associates, Sunderland, Mass) (10, 27). Five thousand bootstrap replicates were done to assess the statistical significance of the phylogenetic tree generated by the NJ algorithm. Analysis based upon a related ST (BURST analysis) was performed with the START program (http://outbreak.ceid.ox.ac.uk/software.shtml) to determine the MLST lineages.
(iii) pbp1a and cpsB. A 485-bp fragment from positions 1498 to 1982 of the pbp1a gene was sequenced with the following primers: pbp1a (F), 5'-GCAAGTAGTGAAAARATGGCTGCTGC, and pbp1a (R), 5'-GACTGTGAAGTTGAACTWTCTGATG-3'. A 476-bp fragment from positions 72 to 547 of the cpsB gene was analyzed with published primers (3). The cpsB and pbp1a genes were amplified from the genomic DNA with the following PCR parameters: initial denaturation at 94°C for 4 min, followed by 35 cycles of denaturation at 94°C for 45 s, annealing at 50°C for 30 s, and elongation at 72°C for 1 min. A 7-min extension at 72°C followed the final cycle. PCR products were sequenced with a Big Dye terminator cycle sequencing ready reaction kit (Applied Biosystems) and run on a model 3700 DNA sequencer (Applied Biosystems). Both the forward and the reverse strands were sequenced. Raw sequences were aligned with Sequencher and ClustalX. The cpsB genes from 23 strains representing 13 serotypes, obtained from GenBank, also were analyzed. For the serotype 37 strain, a base-pair insertion was present at position 212; thus, a 477-bp fragment was analyzed.
|
|
|---|
![]() View larger version (48K): [in a new window] |
FIG. 1. PFGE and MLST for serotype 9V and 14 strains and the 9V-3 clone.
|
|
View this table: [in a new window] |
TABLE 2. PFGE-based clonal groups, MLST alleles, Pbp1a profiles, and penicillin susceptibility patterns for serotype 14 and 9V strains and the 9V-3 clone
|
![]() View larger version (40K): [in a new window] |
FIG. 2. Polymorphic sites within the pbp1a genes of strains R6, VII to IX, I to III, and XII and the 9V-3 clone. Identical sites are represented by dots. The polymorphic sites are numbered in vertical format. From positions 1498 to 1710, the serotype 14 variants had 37-bp differences from the 9V-3 clone and 45-bp differences from strain R6. From positions 1711 to 1874, the serotype 14 variants were identical to the 9V-3 clone and had 41-bp differences from strain R6. From positions 1875 to 1982, the serotype 14 variants had 33-bp differences from the 9V-3 clone and were identical to strain R6. Strain XII was identical to the 9V-3 clone. Strains VII to IX were nearly identical to strain R6.
|
The NJ, ML, and MP trees of the concatenated MLST gene sequences produced clusters similar to those produced by BURST analysis. The ML and NJ methods each produced one tree, while the MP method produced two trees, both with the same structure but slightly different branch lengths. NJ tree bootstrapping with 5,000 replications revealed that strains I to III clustered with the 9V-3 clone 65% of the time. The remainder of the time, strain III clustered with either strain VII (18%) or both strain VII and the 9V-3 clone (17%). When amino acid sequences were used to create an NJ tree with 5,000 bootstrapping replications, strain III always clustered with the 9V-3 clone.
All 41 serotype 9V strains had identical cpsB gene sequences. The cpsB sequences of all of the serotype 14 strains, including the Maryland strains, GenBank strains, and international clones, differed by
5 bp. All serotype 14 strains were
16% (78 to 83 of 476 bp) divergent from the 9V-3 clone. The serotype 14 variants could be differentiated from the 9V-3 clone with the cpsB and the pbp1a genes but not by PFGE or MLST.
A bootstrapped NJ tree of the 15 international clones and the 23 strains from GenBank further demonstrated a correlation between serotype and the cpsB sequence. Moreover, the tree graphically displayed the most common sequence for each serotype (Fig. 3). The sensitivity and specificity of using the cpsB sequence as a surrogate for serotype designation were determined. For each serotype, the most common cpsB sequence was chosen as the reference sequence. For example, the identical cpsB sequence among serotype 6B strains occurred in 75% (three of four) of the strains. The sensitivity was defined as the ability of this cpsB sequence to correctly classify 6B strains as belonging to serotype 6B. The specificity was defined as the ability of this cpsB sequence to correctly classify strains from other serotypes as not belonging to serotype 6B. Among serotypes represented by multiple strains, allowing for a
2-bp difference in the cpsB gene resulted in the highest specificity, and allowing for a >2-bp difference resulted in the highest sensitivity overall. The cpsB sequence had a sensitivity of 75 to 100% and a specificity of 94 to 100% for the following serotypes: 6B, 8, 9V, 14, 19A, 19F, and 23F (Table 3). Overall, 95% (84 of 88) of the strains were classified correctly by serotype with the cpsB sequence.
![]() View larger version (24K): [in a new window] |
FIG. 3. NJ tree of 15 international clones and 23 GenBank strains. Strains which clustered 50% of the time with bootstrapping are indicated.
|
|
View this table: [in a new window] |
TABLE 3. Sensitivity and specificity for cpsB sequences at various thresholds of base-pair differences
|
|
|
|---|
Unlike the European strains, all three serotype 14 variants of the 9V-3 clone had a unique pbp1a restriction profile and pbp1a DNA sequence from positions 1498 to 1710 compared to all other serotype 14 strains, strain R6, and the 9V-3 clone. Thereafter, the Baltimore variants were related to the 9V-3 clone and then to strain R6. Although the exact crossover point between 9V-3 and R6 was different from that in the European variants, the pattern was conserved. These data suggest two possibilities. Perhaps additional recombination events in the pbp1a gene occurred after capsular transformation with the 9V-3 clone, causing increased diversity in our strains compared to the European strains. Second, perhaps our variants arose through a recombination event different from that of the European strains. Previous studies demonstrated that variants have arisen on multiple occasions through unique recombination events (2, 3).
The variability that we noted within the pbp1a gene simply may reflect the mosaic patterns associated with the penicillin binding proteins. For example, some of the genetically unrelated serotype 14 strains had pbp1a sequences which were identical to or nearly identical to those of the 9V-3 clone or strain R6, respectively. Compared to the Spanish 23F-1 clone, strain XII has the same ST and a highly related PFGE pattern. The 23F-1 clone has the same pbp1a DNA sequence as the 9V-3 clone (21). Therefore, it is not unexpected that strain XII has the same pbp1a sequence as the 9V-3 clone. Strain R6 was not related to any of our serotype 14 strains by MLST but had a pbp1a sequence highly related to those of five strains. A search of the pbp1a sequences of the serotype 14 variants in GenBank demonstrated that this sequence was 99% related to that of a previously described 23F penicillin-resistant strain (6).
An inherent limitation of all studies focusing on recombination is the uncertainty of the original donor and recipient strains. An earlier study assumed that the pbp1a sequence of the original donor was identical to that of strain R6 and that the original recipient had a sequence identical to that of the 9V-3 clone (2). In contrast, we used all serotype 14 PNSP strains from a population base over a 2-year period and the 9V-3 clone as the donors and the recipient, respectively. Perhaps another suitable comparison group would be serotype 14 penicillin-susceptible strains, since the pbp1a gene may have undergone further recombination events after the capsular recombination event.
MLST assigns equal weights to one nucleotide difference and multiple nucleotide differences within the same locus, thereby minimizing differences within a single locus. However, since recombination events are typically several kilobases in size, this approach appears to be an accurate means for assessing the genetic relatedness of strains (4). In our study, the concatenated gene sequence analysis provided results different from those provided by MLST in only one instance. Strain III, found to be an SLV of the 9V-3 clone by MLST, had 32 base-pair differences within the ddl gene. While the ML and MP trees clustered strains I to III with the 9V-3 clone, the bootstrapped NJ tree of the DNA sequence revealed that these strains clustered together only 65% of the time. On the surface, these data suggest that strain III may not be a serotype 14 variant of the 9V-3 clone. However, evidence clearly indicates that the variation seen within the ddl gene can be due to a hitchhiking phenomenon related to a single recombination event associated with the penicillin binding protein 2b gene (4). Our data support this theory, since strain III is the only penicillin-intermediate serotype 14 variant. In contrast, strains I and II and the 9V-3 clone are penicillin resistant.
We found that nearly a third of the serotype 14 PNSP strains collected over a 2-year period from the Baltimore metropolitan area were variants of the 9V-3 clone that differed from the European variants. Since they were detected over an 8-month period from persons residing in three different counties, it appears that this clone is widespread within the Baltimore metropolitan area. As shown previously, MLST and PFGE had a high degree of correlation and concurred with the concatenated gene sequence the majority of the time. The recombination history for the serotype 14 variants from Baltimore was different from that for the European variants.
As noted previously (2), the cpsB genes of all serotype 14 strains were highly related but markedly different from that of the 9V-3 clone, even among strains that had undergone capsular transformation. Allowing for a 2-bp difference, the cpsB sequence usually correlated with the serotype designation among this small sample of strains. The key limitation of this study was the choice of the reference sequence for each serotype; the reference sequence for each serotype was the most common sequence. Since the compositions of strains in other databases may differ markedly, the validity of the cpsB sequence as a surrogate for serotype could decrease dramatically. Thus, although the cpsB sequence was correlated with the serotype among this small sample of strains, the role of this gene as part of a sequence-based serotyping method may be limited.
This work was supported in part by the State of Maryland, the Centers for Disease Control and Prevention, and career development awards from the National Institutes of Health to M. C. McEllistrem (K23 AI01788-03) and to L. H. Harrison (K24 AI52788).
|
|
|---|
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»