Previous Article | Next Article ![]()
Journal of Clinical Microbiology, September 2003, p. 4194-4216, Vol. 41, No. 9
0095-1137/03/$08.00+0 DOI: 10.1128/JCM.41.9.4194-4216.2003
Division of Bacterial and Mycotic Diseases, Centers for Disease Control and Prevention, Atlanta, Georgia,1 Division of Infectious Diseases, University of Pittsburgh, Pittsburgh, Pennsylvania,2 Acute Disease Epidemiology Section and Division of Public Health Laboratories, Minnesota Department of Health, Minneapolis, Minnesota,3 Department of International Health, Johns Hopkins University Bloomberg School of Public Health, Baltimore, Maryland4
Received 7 February 2003/ Returned for modification 7 April 2003/ Accepted 17 April 2003
|
|
|---|
|
|
|---|
Pulsed-field gel electrophoresis (PFGE) of chromosomal macrorestriction digests and multilocus sequence typing (MLST) are complementary techniques that are well suited for the analysis of large numbers of isolates. PFGE is ideal for the rapid clustering of genetically related isolates. MLST can be used to assign nonsubjective and electronically portable clone identifiers to clusters of related isolates.
We have used PFGE and MLST to clearly delineate strain distribution within a large set of population-based invasive pneumococcal isolates recovered prior to the use of PCV7. Isolates recovered from adults with selected serotypes that are rare in children, or in which reduced antibiotic susceptibility had been previously reported (34), were also included in this analysis.
|
|
|---|
![]() View larger version (28K): [in a new window] |
FIG. 1. PFGEdendrograms within serotypes included in this study accounting for
multiple invasive isolates. Serotypes 1 and 3 are shown. The numbers of
isolates represented in each dendrogram correspond to the total values
for the corresponding serotypes given in Tables
3 and
4. The serotype is
indicated beside each dendrogram. All MLST results that have been
determined are shown with the ST designation and antibiotic resistance
profile (if any), followed by seven numbers specifying the alleles.
Dice coefficient values (percentages) are indicated above each PFGE
dendrogram. The vertical black lines extending various lengths of some
of the dendrograms indicate penicillin-resistant (penR) isolates (MIC,
2.0 µg/ml). Intermediate penicillin resistance (penI)
is depicted by vertical gray lines. All antibiotic resistances among
the isolates are indicated for serotypes that do not include
penicillin-resistant isolates. Divergent type B ddl alleles,
found only among certain penicillin-resistant strains, are indicated by
asterisks. cmR, chloramphenicol resistance; tmpR,
trimethoprim-sulfamethoxazole resistance; ermR, erythromycin
resistance; tetR, tetracycline resistance; cliR, clindamycin
resistance; cefR, cefotaxime resistance; lfxR, levofloxine
resistance.
|
![]() View larger version (31K): [in a new window] |
FIG. 13. Serotypes
33F, 35B, and 38. For details, see the legend to Fig.
1.
|
|
View this table: [in a new window] |
TABLE 1. Pneumococcal
strains of known STs used as PFGE references
|
|
View this table: [in a new window] |
TABLE 2. Clonal
complex or single ST assignment of 1,042 study isolates based on PFGE
and MLST data
|
PFGE
analysis.
The Bionumerics
software (Applied-Maths, Kortrijk, Belgium) was used to make UPGMA
(unweighted pair group method with arithmetic means) dendrograms of
fragment patterns with the dice coefficient. The dice similarity
coefficient was used with optimization and position tolerance settings
of 1.0 and 1.5%, respectively. PFGE-based clusters were defined
as isolates with
80% genetic relatedness on the
dendrogram.
MLST analysis. MLST alleles from www.mlst.net were downloaded for convenient FASTA screening with the Wisconsin version 10.2 package. New alleles were screened at www.mlst.net to ensure that they were still not included in the pneumococcal MLST database. For the 28 new alleles discovered during this work, sequence traces of both strands were electronically submitted to the pneumococcal database curator for approval prior to their inclusion in the database. Clonal complex (CC) or lineage assignment (BURST) was carried out with software downloaded from http://outbreak.ceid.ox.ac.uk/software.shtml and was used as directed, requiring the parameter of five identical targets with at least one other isolate in the set. MLST data were also used for UPGMA-generated dendrograms with the software downloaded from this web site.
|
|
|---|
The PFGE-MLST associations reported here are documented by serotype (Fig. 1 to 13 and Table 2). All of the serotypes that included 20 or more isolates were represented in at least six of the eight states included in the surveillance (data not shown). All MLST results are superimposed on the PFGE dendrograms shown in Fig. 1 to 13.
We attempted to
relate the PFGE clusters (individual patterns differing bya UPGMA dice coefficient of
20%) within specific
serotypes with specific allelic profiles or clusters of related allelic
profiles by applying MLST to representative isolates. When we observed
similar PFGE profiles shared by isolates of differing serotypes, we
subjected representative isolates to MLST to determine if they were
genetically related.
Isolates that had five or more identical
alleles almost always also had similar PFGE patterns with UPGMA dice
coefficients of
80% (intraserotype relationships shown
in Fig. 1 to
13; data not shown for
isolates of differing serotypes). This relationship was
tested as follows. In 19 of the dendrograms including 40 PFGE clusters
(individual PFGE results differing by
20% similarity)
of 2 or more isolates, 121 isolates were subjected to MLST to explore
the relationship between the PFGE cluster and the allelic profile. In
37 of the 40 PFGE clusters, all isolates subjected to MLST had five to
seven allelic identities. An exception was seen within a related PFGE
cluster in serotype 6B, where ST641 and ST384 differed at four of seven
MLST loci (Fig. 4). Within
serotype 4, ST695 and ST205 were both found within a related cluster
and differed at five of seven MLST loci (however, note that the
xpt-113 allele [GenBank accession
no.
AF474331]
is a deletion derivative of xpt-10, found in ST205)
(Fig. 2). In serotype 19A,
ST172 and ST75 displayed three allelic differences (Fig.
10). Two other exceptions
were seen. One large PFGE cluster within serotype 9V included ST156
(found twice), ST166 (found once), ST642 (found once), and ST643 (found
twice) (Fig. 6). Although
all four of these STs lie within BURST CC156 (Table
3; BURST analysis is explained in the next section), ST642 differed by
four loci from ST156 and by three loci from ST166. The
final exception was seen within a large serotype 18C cluster and
concerned ST113 (found four times), ST637 (found once),
ST638 (found twice), and ST648 (found once) (Fig.
9). Although these four
STs lie within CC113 (Table
3), ST648
differs by three loci from both ST113 and ST637. In contrast, there
were no exceptions to the general
80% PFGE relatedness;
five to seven allelic identity relationship were seen within all 32
PFGE clusters tested within serotypes 1, 3, 6A, 7F, 9N, 11A, 12F, 14,
19F, 22F, 23F, 33F, 35B, and 38 and serogroup 15 (Fig.
1,
3,
5,
7,
8,
9,
11,
12, and
13).
![]() View larger version (35K): [in a new window] |
FIG. 4. Serotype
6B. Reference strains representing ST384, ST270, ST185, and ST90 are
underlined. cefR, cefotaxime resistance. For details, see the legend to
Fig.
1.
|
![]() View larger version (31K): [in a new window] |
FIG. 2. Serotype
4. For details, see the legend to Fig.
1.
|
![]() View larger version (40K): [in a new window] |
FIG. 10. Serotype
19A. Reference strains representing ST175, ST41, ST75, and ST226 are
underlined. cliR, clindamycin resistance. For details, see the legend
to Fig.
1.
|
![]() View larger version (31K): [in a new window] |
FIG. 6. Serotype
9V. The reference strain representing ST156 is underlined. For details,
see the legend to Fig.
1.
|
|
View this table: [in a new window] |
TABLE 3. BURST
analysis of reference STs combined with STs obtained from year 1999
invasive study isolates
|
![]() View larger version (28K): [in a new window] |
FIG. 9. Serogroup
15, serotype 16F, and serogroup 18. For details, see the legend to Fig.
1.
|
![]() View larger version (37K): [in a new window] |
FIG. 3. Serotype
6A. For details, see the legend to Fig.
1.
|
![]() View larger version (25K): [in a new window] |
FIG. 5. Serotypes
7F and 9N. For details, see the legend to Fig.
1.
|
![]() View larger version (21K): [in a new window] |
FIG. 7. Serotypes
11A and 12F. For details, see the legend to Fig.
1.
|
![]() View larger version (37K): [in a new window] |
FIG. 8. Serotype
14. Reference strains representing to ST20, ST67, ST9, and ST14 are
underlined. cliR, clindamycin resistance. For details, see the legend
to Fig.
1.
|
![]() View larger version (37K): [in a new window] |
FIG. 11. Serotype
19F. Reference strains representing ST646, ST236, ST271, and ST650 are
underlined. lfxR, levofloxacin resistance. For details, see the legend
to Fig.
1.
|
![]() View larger version (32K): [in a new window] |
FIG. 12. Serotypes
22F and 23F. Serotype 23F reference strains corresponding to ST37,
ST173, ST242, and ST81 are underlined. cliR, clindamycin resistance.
For details, see the legend to Fig.
1.
|
80%
relatedness, only 12 STs representing seven pairs of ST comparisons
(ST641-ST384, ST695-ST205, ST642-ST156, ST642-ST166, ST648-ST113,
ST648-ST647, and ST172-ST75) did not exhibit five or more allelic
identities (Fig. 2,
4,
6, and
9). Note that only two of
these seven pairs represented pairs of isolates that were not found
within the same CC (ST641-ST384 and ST695-ST205) (Tables
3 and
4). Except for these two instances, all of the isolates tested within the
same PFGE cluster were found to be within the same CC, even when they
only had four or fewer allelic identities (Fig.
1 to
13; Tables
3 and
4). Typically
representative of the PFGE cluster-MLST relationship, for example, was
the serotype 14 dendrogram, which revealed seven PFGE clusters of two
or more isolates (Fig. 8).
Within each cluster, the isolates tested (15 in total) invariably had
five to seven allelic identities. Thus, selected isolates and reference
strains within the same serotype exhibiting PFGE profiles with
80% similarity nearly always had five to seven allelic
identities. On the basis of this observation, combined with a BURST
analysis of 187 MLST results (105 different STs) representing the study
isolates and reference strains, we were able to assign MLST-based CCs
or lineages to 1,042 (89.2%) of the 1,168 isolates subjected to
PFGE (Tables 2,
3, and
4). |
View this table: [in a new window] |
TABLE 4. STs
differing by three or more alleles from all other STs found among year
1999 study isolates and reference strainsa
|
Although the relatively limited numbers of STs (from the study isolates plus those of the reference strains in Table 1) used to construct Tables 3 and 4 seriously limited the robustness of BURST analysis for the purpose of determining ancestral STs, we felt that these data could only be presented clearly by limiting this algorithm as much as possible to the STs used in this study. When this was done, in several instances, the ancestral STs could not be assigned from the data set and instead assignments were predicted from BURST analysis of all known STs (CC81, CC156, CC218, CC124, CC473, CC62, and CC177 in Table 3). Despite this limited number of ST samples, reference STs predicted for CC199, CC236, CC662, and CC113 coincided with the ancestral ST generated by using all known allelic profiles (667 STs at the time of this study). Reference STs for CC146 and CC663 were chosen because of an unsurpassed number of SLVs or double-locus variants (DLVs) within their respective complexes; however, these STs were not actually identifiable as ancestral within a BURST analysis expanded to all known strains. Finally, a single CC (CC66) depicted an ancestral ST (based on a BURST analysis of all known STs) that was not actually encountered in this study (ST66), since the three STs in this study were DLVs of each other.
By using the cumulative PFGE and MLST data, we assigned CCs or lineages consisting of single STs to 67 to 100% of the isolates within each individual serotype (Table 2). For each serotype representing multiple isolates, 2 to 20 isolates were subjected to MLST. Single isolates that were sole representatives of serotypes (types 5 and 7B) were also subjected to MLST. A total of 50 CCs and single ST lineages were assigned on the basis of a BURST analysis of the 105 STs from the study isolates and reference strains (Tables 3 and 4). Seventeen CCs were represented by 2 to 10 related allelic profiles, while 33 lineages consisted of single STs with four or fewer allelic identities to other STs found in the sample set of study isolates and reference strains.
Upon performance of BURST analysis on all known allelic profiles (49 new STs plus the 618 STs recorded at www.mlst.net as of 27 November 2002), the number of single allelic profiles depicted as unique lineages in Table 4 was reduced from 33 to 8 (ST376, ST452, ST660, ST665, ST659, ST41, ST175, and ST448 [the latter 3 were associated with reference strains described in Table 1 only]). In addition, the number of CCs consisting of two or more STs was decreased from 17 to 13 when all known STs were included in the BURST analysis. In this manner, CC13, CC663, ST18, ST43, ST173, and ST226 (together constituting 20 STs) were merged into one large complex that included 88 additional STs in the MLST database (data not shown). Similarly, CC37, CC146, CC124, ST90, ST185, and ST384 (together constituting 13 STs) were merged together into one large complex that also included 88 additional STs). Finally, CC81 and CC66 (together constituting 7 STs) were merged into one complex that included 23 additional STs from the MLST database. Since such global relationships based on BURST analysis would be very difficult to depict graphically, and our primary interest was the description of this population-based set of invasive isolates within our surveillance, we present the BURST data based primarily on the 105 STs represented by our isolate set and reference strains (Tables 3 and 4). SLVs and DLVs of indicated ancestral or reference STs are depicted, as are relationships between other STs included within given CCs (solid lines and dotted lines, respectively). We do not depict an UPGMA-generated dendrogram of the MLST results because the method generated easily discernible artifacts (data not shown; see http://outbreak.ceid.ox.ac.uk/software.shtml for details of why this occurs).
Allelic profile associations with specific serotypes. With few exceptions, serotype-ST associations were remarkably consistent between our data set and the extensive database available at www.mlst.net, indicating strong global associations of given STs with specific serotypes. We did find several additional serotype associations within individual STs; however, in the majority of cases, these results were still consistent with the observed serotype associations within the respective CCs. For example, within ST199, serotype 15C had not been previously reported at the database; however, serotypes 15B and 15C are known to interconvert (33) and the association of type 15C isolates with ST199 has recently been reported (24). Although we did not see serotype 19F ST9 and ST485 isolates, there are two other STs within CC13 that were found in type 19F isolates (Table 3). Within CC218, we observed a novel serotype 7F isolate with ST218, which was formerly only associated with serotype 12F (Table 3). Within lineages consisting of single STs (defined by the presence of four or fewer identities with other study isolates and reference strains), we found only one example of our results differing in serotype associations from those documented at the MLST database (Table 4). ST90, formerly found only within type 6B isolates, was found in this study among type 6A and 6B isolates.
New MLST allelic profiles and alleles. During this work, we discovered 28 new MLST alleles not present at www.mlst.net and deposited them in that database. These included the alleles ddl-75 to ddl-78, ddl-82 to ddl-84, ddl-95, ddl-98, ddl-99, ddl-101, xpt-72, xpt-75 to xpt-78, xpt-85, xpt-86, xpt-96, xpt-113, recP44, recP45, aroE41, aroE42, gdh-46, gki-64, and spi-55. Twelve of these new alleles were independently corroborated by others, since they subsequently appeared in STs entered at www.mlst.net after we deposited them.
Seven of the new ddl alleles from ß-lactam-nonsusceptible isolates were relatively highly divergent (indicated by asterisks in Fig. 3 [serotype 6A], 4 [serotype 6B], 6 [serotype 9V], 8 [serotype 14], 10 [serotype 19A], 11 [serotype 19F], 12 [serotype 23F], and 13 [serotype 35B]) from each other and from pre-antibiotic era type A pneumococcal ddl alleles. Type B alleles have about 90 to 97% sequence identity to type A alleles, while type A alleles generally have >99% sequence identity to each other. Our results are in agreement with previous work documenting the presence of divergent ddl sequences in a subset of ß-lactam-resistant isolates because of cotransformation with pbp2b from nonpneumococcal donors (9). We did not encounter type B ddl alleles in any of the ß-lactam-sensitive isolates subjected to MLST, which is also consistent with the "hitchhiking" effect of type B ddl alleles linked with divergent pbp2b sequences (9).
We present 49 new STs, 47 of which were assigned numbers ranging between 622 and 695 and 2 of which are the recently described STs ST451 and ST452 (1), which were originally discovered during this survey of year 1999 ABCs isolates (Fig. 13). In 11 additional instances, unique STs were corroborated by independent researchers, since we encountered these profiles (three of which contained new alleles that we supplied to the database) prior to their becoming publicly available at www.mlst.net (ST376 [xpt-72], ST393, ST395, ST423, ST433, ST460, ST473, ST485, ST496 [aroE42], ST448, and ST570 [ddl-76]). With 10 of these 11 STs, the ST-serotype associations that we observed were in agreement with the independently listed documentation at www.mlst.net for these STs. The single exception was an invasive ST485 isolate of type 19A, while the carriage isolate described at www.mlst.net was type 19F.
Antibiotic-resistant
isolates.
Penicillin-resistant (MIC,
0.12 µg/ml) isolates were restricted among serogroups
6, 9, 14, 19, and 23, with the exception of serotypes 35B (12 isolates)
and 7B (1 isolate) (Fig. 5
and 13 and Table
2). An additional
nontypeable isolate (ST668, described below within the description of
serotype 12F results) displayed intermediate penicillin resistance
(Fig. 7). The majority of
penicillin-resistant isolates (depicted by vertical broad lines in Fig.
3 to
5,
6 to
8, and
10 to
13) were resistant to
various combinations of other antibiotics (depicted only for isolates
subjected to MLST). For this reason, resistance to additional
antibiotics is only indicated in Fig.
1 to
13 for serotypes in which
resistance to any antibiotics occurs
infrequently.
CC or lineage distribution
within individual serotypes. (i) Serotype 1.
Two lineages were apparent within
serotype 1, each represented by a single ST (ST217 and ST227). Six of
the seven year 1999 invasive isolates were available from children
5 years of age in our 1999 invasive disease surveillance. In
addition, we included 15 serotype 1 isolates from individuals ranging
from 7 to 63 years of age. All serotype 1 isolates, except one
trimethoprim-sulfamethoxazole-resistant isolate (Fig.
1), were sensitive to the
eight antibiotics included in this analysis. Two isolates from within
the major PFGE cluster of 20 isolates, and the single PFGE variant,
were subjected to MLST, which revealed two nonrelated STs previously
documented at
www.mlst.net
as being from serotype 1 invasive isolates in other countries (ST217
from Denmark and Sweden and ST227 from the United Kingdom and Denmark)
(Fig. 1). The 20 isolates
within the PFGE cluster that included 2 ST227 isolates were recovered
in eight different states, while the sole unrelated ST217 variant was
recovered from a child in Minnesota. Typical serotype 3 PFGE profiles
highly associated with allelic profile ST180 had fortuitous similarity
(UPGMA dice coefficients of up to 95%) to the single serotype 1
PFGE profile that corresponded to ST217; however, MLST analysis
indicated no relatedness (compare ST217 to ST180 from the serotype 1
and 3 dendrograms shown in Fig.
1). Thus, we found no
evidence of genetic relatedness between serotype 1 isolates and
isolates of other serotypes.
(ii) Serotype
3.
From the PFGE and MLST
data compiled from the 56 serotype 3 isolates from CDC surveillance (6
from children <5 years old and 50 from individuals
5
years old representing all eight survey states), it is apparent that at
least the majority of serotype 3 invasive isolates within the United
States are represented by a single clone with ST180 (Fig.
1 and Table
2), which has been
documented in the global MLST database as being from type 3 isolates in
at least six different countries. Four of the isolates depicted
exhibited resistance to combinations of tetracycline, erythromycin, and
chloramphenicol (Fig.
1).
(iii) Serotype 4. Approximately 85% (61 of 71 pediatric isolates) of type 4 isolates were projected to represent three distinct lineages represented by single STs (ST205, ST244, and ST695), each with four or fewer identical alleles in common with other STs seen in this study. BURST analysis using all known STs clustered ST244 and ST695 within the same lineage. ST695 (three-locus variant of ST244) was characterized by the novel xpt-113 allele with an in-frame deletion of 39 codons. ST695 was obtained from five independent isolates and represented approximately 72% of the type 4 isolates (three PFGE clusters and an additional isolate with a divergent PFGE profile) (Fig. 2 and Table 2). With the exception of the 39-codon deletion, xpt-113 is identical to xpt-10, which is common to both ST244 (two isolates) and ST205 (eight isolates). ST695 is an SLV (xpt-10 to xpt-113) of ST246 recorded at www.mlst.net as being from type 4 invasive isolates recovered in the United Kingdom during 1997.
One ST695 isolate was erythromycin resistant, and another isolate within the ST205 PFGE cluster was trimethoprim-sulfamethoxazole resistant. Two additional isolates of undetermined genetic types were trimethoprim-sulfamethoxazole resistant and tetracycline resistant.
(iv) Serotype 5. A single trimethoprim-sulfamethoxazole-resistant serotype 5 pediatric isolate was ST653, a DLV of a type 5 isolate recovered in Colombia in 1994 (www.mlst.net).
(v) Serotype 6A. Fifty-three percent (31 isolates) of the 58 serotype 6A isolates were penicillin resistant. Fifty-seven of the 58 isolates were recovered from pediatric patients. The sole isolate from an adult patient was the result of an isolate mistakenly initially serotyped as 9N. This was a genetically heterogeneous isolate set represented by at least six distinct genetic lineages deduced from 11 MLST results (Fig. 3 and Tables 3 and 4). We discovered four of these STs independently (ST376, ST460, ST473, and ST395) prior to their appearance at www.mlst.net, which documented serotype 6A isolates recovered between 1994 and 1998 (ST376 [from a case of meningitis; blood isolate], ST460 [unspecified diagnosis; sputum isolate], ST473, and ST395 [both carriage isolates]). ST376 and ST460 putatively accounted for more than 50% of serotype 6A isolates (Fig. 3 and Table 2). ST376, circumstantially representing 19 isolates (Fig. 3), was characterized by the new type B ddl-77 allele and was obtained from 4 independent isolates. ST460 was obtained from 2 isolates representing a PFGE cluster of 13 penicillin-sensitive isolates. CC473 isolates (ST473 and its SLV ST647) represented five penicillin-resistant isolates. ST395 represented two penicillin-resistant isolates. ST660 and ST90, both with four or fewer similarities to all of the other STs in this study, were found in single isolates that had PFGE profiles dissimilar to other type 6A isolates. ST90 was additionally found among serotype 6B isolates.
(vi) Serotype 6B. One hundred three serotype 6B isolates (all pediatric) were analyzed, constituting at least six CCs or single ST lineages (nine STs). Forty-two percent (n = 43) of these isolates were penicillin resistant. CC663 ST138 (documented as being from invasive serotype 6B isolates from three different countries in the MLST database) and ST639 represented 24 (23.3%) of these isolates, all but 1 of which were penicillin sensitive.
Isolates with PFGE patterns and/or STs similar or identical to those of reference strains Spain6B-2 (ST90; 13 isolates), MD6B-17 (ST384; 19 isolates), and S. Africa6B-8 (ST185; 1 isolate) were also identified (Fig. 4 and Table 2).
ST644 (CC156) was found within single serotype 6B and 19F isolates (Table 3).
(vii) Serotype 7B. Only one serotype 7B isolate was analyzed, which was mistakenly included in this study as a serotype 7F isolate (bottom of Fig. 5 serotype 7F dendrogram). This isolate was intermediately penicillin resistant and was resistant to three additional antibiotics. The ST of this isolate (ST644) was found among CC663 isolates that include serotypes 19A and 6B (Table 3).
(viii) Serotype 7F. Thirty-one serotype 7F isolates were analyzed (seven from pediatric patients), all of which were antibiotic susceptible (Fig. 5). Twenty-nine (93.5%) of these isolates putatively have STs identical or highly similar to ST191, which was previously found among invasive serotype 7F isolates that have been recovered from multiple countries (see www.mlst.net). One serotype 7F isolate, with a PFGE profile markedly different from those of all other serotype 7F isolates, was ST218 (within CC218), which was otherwise found only in serotype 12F isolates (Table 3).
(ix) Serotype 9N. Only 2 of the 13 serotype 9N isolates were from pediatric patients. These isolates and nine other isolates from adult patients made up the major PFGE cluster, which appeared to be composed of CC66 isolates (ST631 and ST632) (Fig. 5).
(x) Serotype 9V. Forty-two of the 98 serotype 9V isolates were initially mistakenly tested as belonging to other serogroup 9 serotypes (9N and 9A). In this manner, we fortuitously included 26 adult isolates in this analysis (Fig. 6). About 92% of these isolates (both pediatric and adult) were found to be within CC156 (Table 3,ST156, ST166, ST642, and ST643) and closely related to internationally disseminated clone Spain9V-3 (ST156). More than 63% of the CC156 isolates examined were penicillin resistant (Table 2). Two independent penicillin-resistant isolates were found to have in common ST644, a DLV of ST156 (Table 3) characterized by the previously undocumented ddl-98 and xpt-76 alleles. A single ST166 isolate was found. ST166, an SLV of ST156, was initially found in a multiresistant isolate recovered in Taiwan (see www.mlst.net). ST642 (which also carried new allele xpt-76) differed from the other STs found in CC156 isolates by two to four alleles.
(xi) Serotype 11A. Twenty-two of the 25 serotype 11A isolates (88%, all but one from adults) were putatively within CC62 (Table 3), which consists of ST62 (recorded at www.mlst.net as being from type 11A meningitis isolates recovered in Spain) and an SLV of ST62 characterized by the new recP45 allele (ST630) (Fig. 7). Antibiotic resistance within CC62 was restricted to one trimethoprim-sulfamethoxazole-resistant isolate. A single tetracycline-resistant isolate displaying a divergent PFGE profile was typed as ST662, which lies within CC662 (Table 3). This isolate, although originally typed as 11A, was subsequently found to be nontypeable. It is interesting that CC662 appears to consist primarily of type 18C and 33F isolates (along with a single type 23F isolate) and that, on the basis of BURST analysis, ST662 is the ancestral profile.
(xii) Serotype 12F. The 60 serotype 12F isolates (11 from pediatric patients) were primarily within CC218 (Table 3), consisting of ST218 and ST220 (Fig. 7 and Table 2). Both of these profiles are associated solely with invasive type 12F isolates at www.mlst.net (ST218 reported as being from several different countries and ST220 reported as being from a single isolate recovered in Denmark). One isolate initially typed as 12F had a divergent PFGE profile and divergent ST668, found within CC199 (Table 3). Subsequent reserotyping revealed this strain to be nontypeable (Fig. 7).
(xiii) Serotype 14. Serotype 14 is most often associated with invasive isolates recovered from children and consisted of 208 isolates recovered from pediatric patients (63.5% of the year 1999 type 14 invasive isolates recovered from children in our survey area).
On the basis of MLST results from seven serotype 14 isolates within a major related PFGE cluster (also including the ST9 reference strain England14-9), the single largest serotype 14 CC (96 isolates = 46.2%) was CC13, in which 81.2% of the strains (78 isolates) were penicillin resistant (Fig. 8 and Tables 2 and 3). Besides ST9 (also found in penicillin-sensitive, internationally dispersed clone England14-9), these CC13 isolates were characterized by ST13 (characterized by both penicillin-resistant and penicillin-sensitive isolates in the MLST database), ST16 (listed at www.mlst.net as being from a single pneumonia-causing serotype 14 isolate in Canada), ST622 (SLV of ST13 characterized by new allele xpt-78), and ST623 (DLV of ST13 associated with two isolates on the basis of PFGE) (Table 3).
The majority of penicillin-susceptible isolates (78 [71%] of the 110 susceptible isolates) appeared to be within CC124, apparently primarily with ST124 in common (seven independent isolates had ST124). ST124 is documented in the MLST database as being from invasive type 14 isolates recovered from several different countries. An SLV of ST124 (ST656), characterized by the new aroE41 allele, was found within a single multiresistant isolate that displayed a divergent PFGE profile (Fig. 8, top of type 14 dendrogram).
Five multiresistant isolates were represented by ST671. This was the only ST representing a type 14 strain within the type 9V-dominated CC156 group (Table 3). This was not an unexpected result, since ST156 (the ancestor ST found in Spain9V-3) has also been documented among type 14 isolates (7, 25). Fourmultiresistant isolates putatively had in common ST67 (CC66 in Table 3) on the basis of two isolates with PFGE profiles highly similar to that of reference strain Tennessee14-10 (ST67), and another pair of isolates had highly similar PFGE profiles (one of which was directly determined to be ST67) (Fig. 8).
CC199 serotype 14 isolates were composed of an ST200 tetracycline-resistant isolate with a unique PFGE profile and the SLV of ST200, ST667. ST200 has been previously documented in serotype 14 invasive isolates recovered in Denmark and Taiwan. One multiresistant isolate had a PFGE profile identical to that of reference strain Spain14-5 (ST18) and was directly typed as ST18, which has four or fewer allelic identities with other STs in the sample set (Table 4). A single multiresistant ST633 isolate (CC81) was also observed.
(xiv) Serotypes 15B and 15C. The six serogroup 15 isolates from pediatric patients that were available to us (from patients in Georgia, Minnesota, and Connecticut) had highly similar PFGE profiles and ST199 (the ancestral ST of CC199 in Table 3 on the basis of study STs and a cumulative BURST analysis of all known STs) (Fig. 9). ST199 was also found among serotype 19A and 19F study isolates and is documented at www.mlst.net as being from serotype 15B and 19A invasive isolates recovered in the United Kingdom and The Netherlands. As recently reported, ST199 is apparently commonly found within both serotype 15B and 15C pneumococci recovered from nasopharyngeal specimens of children in the United Kingdom (24).
(xv) Serotype 16F. One of the two serotype 16F isolates was from an adult and was originally mistakenly identified as serotype 9A. Reserotyping this isolate revealed that it was serotype 16F. The ST of this isolate, which was characterized by our newly discovered ddl-76 allele, was independently corroborated on 21 October 2002 at www.mlst.net as being from a serogroup 16 isolate recovered in Scotland. The other serotype 16F isolate (from a child) was the unique ST659 (Fig. 9). The closest known matches to ST659 have only three alleles in common and are from a variety of serotypes, none of which are type 16F.
(xvi) Serotype 18C. Ninety-three serotype 18C isolates from children were analyzed, five of which were mistakenly included and subsequently found to be serotype 18B (three isolates) or 18F (two isolates). These five isolates had PFGE profiles very similar to those of the majority of type 18C isolates. The major PFGE cluster, consisting of 74 isolates (79.6% of the total) represents CC113 (Fig. 9, all 9 isolates tested). On the basis of a BURST analysis of the study STs and all known STs, ST113 is the ancestor ST within this CC, which also includes SLVs ST637 and ST638 (three independent isolates) and satellite ST648 (three-locus variant of ST113) (Table 3). This observation is consistent with ST113 being observed in four of the nine isolates within this PFGE cluster that were sequence typed. ST113 is documented as being from isolates dating back to 1980 that were recovered from meningitis isolates in the United Kingdom, The Netherlands, and Spain. ST496 (within CC662, which additionally represents type 33F, 11A, and 23F study isolates) was characterized by our new aroE42 allele and was found in four isolates within a total cluster of eight isolates that had highly similar PFGE profiles. This ST was independently entered at www.mlst.net on 10 July 2002 as being from a type 18C carriage isolate recovered in 1996 in Finland.
Single isolates with divergent PFGE profiles were also found. A unique trimethoprim-sulfamethoxazole-resistant ST666 isolate was characterized by the new xpt-96 allele and is a DLV of STs found in serotypes 14, 18F, and 21. An ST632 isolate (CC66 shared with type 14 and 9N isolates), characterized by the new xpt-75 allele, was additionally found in two serotype 9N isolates in this study.
(xvii) Serotype 19A.
Eighty-two serotype 19A isolates were
analyzed, 18 of which were recovered from pediatric patients. For 54
(65.8%) of these isolates, the penicillin MICs were
0.12 µg/ml. On the basis of PFGE and MLST data shown
in Fig. 10, 72%
(n = 59) of the type 19A isolates were within CC199
(see serogroup 15 above) and typed as ST199 (two independent isolates)
and ST645 (an ST199 SLV) (Table
3). The lower 30%
of the serotype 19A dendrogram shown in Fig.
10 represents a diverse
set of reference strains and isolates representing at least three
distinct CCs or single ST lineages (ST663, ST172, and ST75 representing
CC663; ST665, which has four or fewer allelic identities to all known
STs; and ST485 within CC13). ST172 and ST485 are listed at
www.mlst.net
as being from serotype 19A meningitis and serotype 19F carriage
isolates, respectively.
(xviii) Serotype 19F. A set of 110 serotype 19F pediatric isolates were analyzed (Fig. 11) revealing a diverse population of at least 10 different CCs or single ST lineages (17 different STs from study isolates recorded in Table 2). Approximately 21% of the isolates were categorized as CC236 (23 multiresistant isolates), consisting of ST236 (ancestral sequence, associated with major resistant clone Taiwan19F-14), ST271 (SLV of ST236 documented at www.mlst.net as being from invasive multiresistant strains in Korea and the United States), ST652 (SLV of ST236 with new allele xpt-86), ST651 (SLV of ST236), and ST635 (DLV of ST271 characterized by new allele spi-55 and new type B allele ddl-78) (Table 3).
Another 20% of serotype 19F isolates (22 isolates, 21 of which were penicillin sensitive) were assigned to CC251, which consists of ST251 (found in three isolates, one of which was intermediately penicillin resistant) and ST654 (SLV of ST251 found in one isolate). ST251 is documented as being from an antibiotic-sensitive type 19F meningitis isolate recovered in the United Kingdom in 1995.
About 12% (13 multiresistant isolates) of serotype 19F isolates were within CC81 since two isolates within this cluster had ST81 in common. ST81 is found within reference strain Spain23F-1 and was also presumed to be within multiple serotype 23F study isolates on the basis of PFGE profiles (Fig. 11).
Ten penicillin-sensitive isolates (9.1% of serotype 19F isolates) were within CC199 (ancestral ST199 and DLV ST649, Tables 2 and 3).
Multipleserotype 19F isolates were also represented within CC13 (ST425 [four isolates] and ST655 [two isolates]) and CC177 (two isolates; one isolate was typed as ST177 and, with another isolate, had a PFGE profile similar to that of ST646 [SLV of ST177], which is found in the major Chile multiresistant type 19F reference strain) (Fig. 11). ST177 is recorded in the MLST database as coming from type 19F isolates recovered in Spain (associated with meningitis) and Portugal.
A single trimethoprim-sulfamethoxazole-resistant, intermediately penicillin-resistant serotype 19F representative of CC156 was detected (ST644). ST644 was also found within serotype 6B (Tables 3 and 4 and Fig. 11). ST43 (from a tetracycline-resistant study isolate and recorded at the MLST database as coming from a 19F isolate causing meningitis in the United Kingdom) and ST423 (a lineage 2 isolate also independently discovered in a recent United Kingdom carriage isolate) were also detected in single isolates.
(xix) Serotype 22F. Thirty-seven serotype 22F isolates were analyzed, six of which were from pediatric patients. All but one isolate was predicted to be ST433 on the basis of MLST of six isolates and the observation of related PFGE profiles (Fig. 12). ST433 was also independently discovered and associated with a carriage isolate recently recovered in England (see www.mlst.net).
(xx) Serotype 23F. Forty-six (71.9%) of the 64 serotype 23F pediatric isolates were predicted to be within CC37, consisting of ST37 (ancestral type); SLVs ST33, ST626 (characterized by new allele ddl-83), and ST627; and ST37 DLVs ST628 (new allele ddl-95), ST629, and ST625 (Fig. 12 and Table 3). More than half of these lineage 4 type 23F isolates were penicillin resistant (Table 2). Eleven multiresistant isolates were presumed to be in CC81 (ST81) (Table 3) on the basis of PFGE and resistance profiles related to reference strain Spain23F-1 (Fig. 12). A set of four ST242 isolates (three of which were sequence typed) were identified that displayed PFGE and resistance profiles highly related to multiresistant, invasive reference strain Taiwan23F-15 (30). A single isolate with a divergent PFGE profile was determined to have new allelic profile ST661, within CC662, which is composed primarily of serotype 33F and 18C isolates (Table 3).
(xxi) Serotype 33F. All seven serotype 33F pediatric isolates were predicted to be within CC662 (ST100, which was associated with meningitis in Spain, and DLV ST636) (Fig. 13). CC662 was common to several serotype 18C (ST496) isolates and a single 23F isolate (Table 3).
(xxii) Serotype 35B. The 16 serotype 35B isolates were divided into two newly discovered STs that were divided into two corresponding PFGE clusters (Fig. 13). All 12 penicillin-nonsusceptible isolates, including one pediatric isolate, were ST451, which was characterized by the new recP44, xpt-77, and ddl-82 (type B) alleles. This clone has been discovered among invasive U.S. isolates recovered during other years prior and subsequent to 1999 and is further described elsewhere (1). It is interesting that SLVs of ST451 from type 35B penicillin-resistant isolates have appeared at the MLST database. One of these is additionally documented within a serotype 14 isolate. The four penicillin-sensitive isolates were predicted to be ST452, which has also been detected among U.S. invasive isolates during other years (1).
(xxiii) Serotype 38. Four of the five serotype 38 pediatric isolates were predicted to have in common unique ST393 (Fig. 13), which was also discovered independently in at least one recent type 38 United Kingdom carriage isolate (www.mlst.net).
|
|
|---|
A second, unanticipated, purpose of this analysis was the use of MLST to detect serotyping errors. Although there were relatively few instances within most serogroups, we did uncover a major error involving a large number of year 1999 serotype 9A and 9N isolates that were actually found to be type 9V. We were alerted to this after first noting that several of these type 9A and 9N isolates were ST156 (Spain9v-3 type strain) or ST156 derivatives and after noting that ST156 has never been found previously within types 9A and 9N.
The increasing characterization of pneumococci collected worldwide through the use of MLST became increasingly evident to us during this study. For example, four STs (now designated ST376, ST460, ST473, and ST395) appeared to represent more than 60% of the serotype 6A isolates and were discovered during our work prior to their appearance at www.mlst.net. Several other, similar, instances were observed during the course of this project. We discovered a penicillin-resistant type 35B clone disseminated across the United States that carried three new MLST alleles that we submitted to the MLST database. Within a few months, an SLV was reported at the MLST database as being from Sweden and two SLVs have appeared in other areas within the United States. These observations reflect the fact that, through the use of MLST, there is effective documentation of a high percentage of circulating clonal types. This nonsubjective electronic transfer of precise genetic information is a necessary first step in understanding clone-specific genetic traits differentially affecting the epidemiology and virulence of various globally disseminated strains. The addition of 49 new STs found during this work, combined with 618 STs already recorded at www.mlst.net, brought the total number of known STs to 667 as of 27 November 2002.
It is known that multiple penicillin-resistant strains within individual serotypes arose during the past 25 years through steps requiring multiple horizontal DNA transfer events and point mutations. It is also known that there are multiple existing examples of strains within indistinguishable genotypes expressing distinct serotypes that appear to have arisen through a simpler mechanism of horizontal transfer requiring a single double-crossover event. Such strains have only been examined extensively within PCV7 serotypes. It is plausible that successful clones currently targeted by potent vaccines could gradually shift to alternate nontargeted serotypes through genetic transfer events, analogous to the major shift in the frequency of individual successful penicillin-resistant strains that occurred during the past 25 years. For example, a penicillin-nonsusceptible type 24F strain has recently been characterized within an ST previously only documented in invasive, intermediately penicillin-resistant type 14 isolates (29).
Vaccination tends to cause a shift toward carriage of strains of nonvaccine serotypes. Vaccine trial results indicate that PCV7 will reduce the carriage of vaccine serotype strains (5, 6). Among children who were immunized with PCV7, a marked reduction in the number of episodes of acute otitis media occurred; however, the proportion of isolates within nonvaccine serotypes significantly increased (11). The degree to which transformation of present vaccine serotype strains to nonvaccine serotypes will contribute to this shift is unknown. The only means by which to detect such occurrences is precise genetic surveillance. We have provided a broad database of genetic types that will greatly facilitate our future invasive pneumococcal surveillance within the United States. This database will allow us to compare CC (or single ST lineage) and serotype associations from invasive population-based isolate sets from subsequent years (and possibly prior years) to assess the potential selective effects of PCV7.
We and others have found several examples of multiple serotypes sharing the same ST; however, these examples have been primarily limited to involve serogroups 6, 9, 19, 14, and 23. This is still a reason for concern, since PCV7 only includes types 6B, 9V, 14, 19F, and 23F, and the potential protective effects of PCV7 against other serotypes included within serogroups (i.e., 6A, 9A, 9N, 19A, and 23A) have not been well studied. It is believed that the common basis for clinical isolates of varied serotypes having the same genotype is horizontal transfer of serotype-specific genes; however, this has only been shown at the molecular level in a limited number of instances involving interserogroup transformation events (3, 4). There are no data indicating that intrachromosomal changes could shift serotypes; however, shifts changing the serotype without changing the serogroup have not been examined at the DNA sequence level. This appears likely to be the case with the documented interconversion of types 15B and 15C (24, 33). In the near future, sequence information from more of the capsular biosynthetic operons will help to resolve this issue.
CC199, CC13, CC81, CC37, CC66, CC663, CC662, and ST90 represented sets of isolates with both PCV7 and non-PCV7 serotypes, indicating the possibility that horizontal transfers involving an extensive array of both PCV7 and non-PCV7 serotype-specific genes have occurred in the past (Table 3).
ST199 (15B, 15C, 19A, and 19F) and ST90 (6A and 6B) are two examples of STs shared between study set isolates with PVC7 serotypes and one or more non-PCV7 serotypes. Many more examples exist in the MLST database, including ST66 (9N and 19F) and ST81 (19F, 23F, 19A, 9V, and 14) (Table 3). There were only two examples of STs representing non-PCV7 types outside of PCV7 serogroups in combination with other serotypes, but this may be due to the fact that fewer non-PCV7 serotype isolates were included in this analysis. These included ST199 (15B, 15C, 19A, and 19F) and ST218 (12F and 7F) (CC199 and CC218 in Table 3).
It has been estimated that MLST alleles are changed 10-fold more frequently through recombination than by mutation in pneumococci (12); however, the ranges of ST decay rates among pneumococcal populations in nature are not known. The data shown (Fig. 1 to 13 and Tables 2, 3, and 4) provide a pre-PCV7 baseline for determining ST decay rates and for determining serotype shifts within these CCs and single ST lineages.
In summary, this work is necessary to properly monitor horizontal transfer events resulting in serotype changes within virulent clones and the generation of new antibiotic-resistant strains. Currently, pneumococcal clone expansion (and decay) and the emergence of new strains can only be accurately assessed through a system based on MLST. For future proper vaccine formulation strategies within the United States, it is necessary to put into place the most precise surveillance possible. The present effort has provided a reference strain network that will greatly decrease the effort involved in subsequent large genetic surveillance efforts in the United States.
We are very grateful to State Laboratory members for providing the isolates and demographic information concerning the isolate set. We especially thank Carolyn Wright and other CDC ABCs-EIP members for providing isolate information. We greatly appreciate the efforts of Angela Brueggemann in verifying our new allele sequences and providing designations for new alleles and STs, and we thank all of the others who are involved in the maintenance of the global pneumococcal MLST database at www.mlst.net. We thank Brian Spratt for providing unpublished ST data. We are grateful to Delois Jackson and Alma Ruth Franklin in the CDC streptococcal serology laboratory for serotyping and MIC results. We thank James Jorgensen and his laboratory personnel for some of the MIC data and Elizabeth Blanton for her summer of work in the CDC laboratory.
|
|
|---|
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»