Skip to main content
Advertisement
  • Loading metrics

Vaccine Escape Recombinants Emerge after Pneumococcal Vaccination in the United States

Abstract

The heptavalent pneumococcal conjugate vaccine (PCV7) was introduced in the United States (US) in 2000 and has significantly reduced invasive pneumococcal disease; however, the incidence of nonvaccine serotype invasive disease, particularly due to serotype 19A, has increased. The serotype 19A increase can be explained in part by expansion of a genotype that has been circulating in the US prior to vaccine implementation (and other countries since at least 1990), but also by the emergence of a novel “vaccine escape recombinant” pneumococcal strain. This strain has a genotype that previously was only associated with vaccine serotype 4, but now expresses a nonvaccine serotype 19A capsule. Based on prior evidence for capsular switching by recombination at the capsular locus, the genetic event that resulted in this novel serotype/genotype combination might be identifiable from the DNA sequence of individual pneumococcal strains. Therefore, the aim of this study was to characterise the putative recombinational event(s) at the capsular locus that resulted in the change from a vaccine to a nonvaccine capsular type. Sequencing the capsular locus flanking regions of 51 vaccine escape (progeny), recipient, and putative donor pneumococci revealed a 39 kb recombinational fragment, which included the capsular locus, flanking regions, and two adjacent penicillin-binding proteins, and thus resulted in a capsular switch and penicillin nonsusceptibility in a single genetic event. Since 2003, 37 such vaccine escape strains have been detected, some of which had evolved further. Furthermore, two new types of serotype 19A vaccine escape strains emerged in 2005. To our knowledge, this is the first time a single recombinational event has been documented in vivo that resulted in both a change of serotype and penicillin nonsusceptibility. Vaccine escape by genetic recombination at the capsular locus has the potential to reduce PCV7 effectiveness in the longer term.

Author Summary

The 7-valent pneumococcal conjugate vaccine is a remarkable public health success story. It has significantly reduced invasive pneumococcal disease in the United States not only by protecting vaccinated children, but also by protecting unvaccinated older children and adults by herd immunity. However, there was always a concern that use of a limited-valency vaccine would result in an increase in disease due to nonvaccine serotypes, and this has now occurred in the US. The predominant nonvaccine serotype causing invasive disease is 19A, and this increase is partially explained by “vaccine escape” pneumococci, strains that have exchanged the genes that encode a vaccine serotype 4 capsule for genes that encode a nonvaccine serotype 19A capsule. These strains are then able to escape vaccine-induced immunity. Characterisation of the genetic event that resulted in these vaccine escape strains was the focus of our study and the results were surprising. The results of this study have important relevance to the long-term effectiveness of the current vaccine and to the development of future pneumococcal vaccines.

Introduction

Streptococcus pneumoniae (the “pneumococcus”) is one of the most important bacterial pathogens worldwide, especially among children. Pneumococcal pneumonia, meningitis, and septicemia result in 1 million deaths annually among children <5 y of age [1]. PCV7 protects against seven pneumococcal capsular types (serotypes)—4, 6B, 9V, 14, 18C, 19F, and 23F [2]—and has been used to vaccinate children in the United States since 2000. PCV7 has been remarkably effective in reducing disease among vaccinated children, and even among unvaccinated children and adults as a result of a striking herd immunity effect, due to the disruption of pneumococcal transmission from young children to older children and adults [26]. PCV7 is increasingly being used to vaccinate children in other countries.

There are 91 known pneumococcal serotypes, the last of which was discovered recently [7,8]. The ecological niche for the vast majority of the pneumococcal population is the nasopharynx of healthy children [9]; thus, any serotype-specific vaccine that is of limited valency and affects nasopharyngeal carriage will perturb the composition of the circulating pneumococcal population, with unknown consequences. The capsule is the principal known virulence factor with respect to invasive pneumococcal disease [10], and population biology studies indicated that certain serotypes have a greater potential to cause invasive disease than others [11,12]. Another important aspect of pneumococcal biology is that there is a strong association between genotype as defined by multilocus sequence typing (MLST), and serotype—that is, strains with the same MLST sequence type (ST) are usually of the same serotype [11,13,14].

Since the capsule is the principal invasive disease determinant, and is the target for serotype-specific prevention of disease by vaccination, there are two events that are especially important with respect to understanding vaccine effects: serotype replacement and capsular switching. At a population level, serotype replacement simply refers to a decrease in the prevalence of vaccine serotype pneumococci in the nasopharynx accompanied by a corresponding increase in nonvaccine serotype pneumococci, as they fill the newly vacant ecological niche. Serotype replacement in the nasopharynx of a healthy child may or may not be problematic; the public health concern is whether or not replacement serotypes also cause disease. Serotype replacement invasive disease among children and adults has significantly increased in the US post-vaccination [1518], and invasive pneumococcal disease among children <5 y of age is now predominantly due to serotype 19A [17]. A recent JAMA report described a significant overall increase in nonvaccine serotype disease among Native Alaskan children, most frequently due to serotype 19A [18]. Genotypic characterisation of these serotype 19A strains by MLST showed that most of the serotype 19A replacement disease can be explained by clonal expansion of one genotype, ST199, which existed prior to vaccination [15,18].

The second major vaccine-related concern is the possibility of capsular switching, when the genes encoding one type of capsule are exchanged, via transformation and recombination, for the genes encoding a different type of capsule. Capsular switches from one vaccine serotype to another were first described 16 y ago [1921], but it is the vaccine-to-nonvaccine serotype switch that is of primary concern, because it contributes to serotype replacement and allows for the possibility of vaccine escape. Acquisition of a nonvaccine capsule by a pneumococcal strain capable of causing invasive disease has been a serious concern related to the use of any serotype-specific vaccine [11,22].

The capsular locus of the pneumococcus is located between two genes, dexB and aliA ([7]; Figure 1). (Serotype 37 is one exception, as it has a defective capsular locus and the serotype is determined by the type 37 synthase gene [tts] located elsewhere in the pneumococcal genome [7,23].) Two of the six penicillin-binding proteins (PBPs) possessed by the pneumococcus are near the capsular locus: pbp2x is upstream of dexB and pbp1a is downstream of aliA. Alterations in PBPs confer penicillin resistance, and alterations in pbp2x, pbp1a, and pbp2b are the most important [2426]. Penicillin-resistant pneumococci are a major problem throughout the world [27]; prior to use of PCV7 in the US, nearly one-quarter of invasive pneumococci were penicillin-nonsusceptible [28]. Initially, use of PCV7 in the US reduced the incidence of invasive disease due to antimicrobial-resistant vaccine serotypes, but more recently, antimicrobial resistance has increased with the increase of serotype 19A disease [28,29].

thumbnail
Figure 1. Schematic of the Pneumococcal Genome

Housekeeping genes characterised in the MLST scheme (aroE, gdh, gki, recP, spi, xpt, and ddl), the genes for PBPs (pbp2x, -1a, -3, -2b, -2a, and -1b), the capsular locus flanked by the dexB and aliA genes, and the region around the capsular locus that was sequenced in this study (shown also in expanded version) are depicted. The capsular loci of serotypes 4 and 19A are 20.9 kb and 18.6 kb, respectively, in length.

https://doi.org/10.1371/journal.ppat.0030168.g001

The Centers for Disease Control and Prevention (CDC) has been monitoring invasive pneumococcal disease since 1995 through the Active Bacterial Core (ABC) surveillance program [6,14,30] and as a result, the post-vaccination increase in nonvaccine serotype 19A disease in the US was quickly detected. Serotype 19A strains collected by the CDC through 2005 were genotyped by MLST, which revealed that vaccine escape strains had begun to emerge in 2003 [14,15]. These strains possessed an MLST genotype, ST695, that had always been associated with vaccine serotype 4 (ST6954), but now expressed a serotype 19A capsule (ST69519A). These strains were detected only 3 y after vaccine implementation, but rapidly increased in prevalence. The first three strains were detected in 2003; two strains were detected in 2004; and 32 strains were detected in 2005, some of which had evolved further. Moreover, in 2005, two new types of serotype 19A vaccine escape strains emerged, ST236519A (n = 4) and ST89919A (n = 1); these appeared to represent new recombinational events that also occurred between serotype 4 recipients and serotype 19A donors. The aim of this study was to sequence the regions upstream and downstream of the capsular locus, including both PBPs, to identify the putative recombinational event(s) that resulted in these vaccine escape strains.

Results

Recipient, Donor, and Progeny Pneumococcal Strains Selected for Characterisation

Strains were selected based upon serotyping and MLST genotyping data and are described in Table 1. Between 1998 and 2005, 31,669 cases of invasive pneumococcal disease among all ages were identified from ABC sites, from which 28,363 pneumococcal strains were available for serotyping. 15,736 (55%) strains were of nonvaccine serotypes, and 1,935 of those were serotype 19A, 88% of which were recovered from 2000 to 2005 ([17] and ABC surveillance program, unpublished data). The following serotype 19A strains were genotyped by MLST: 82 of 113 (73%) strains from 1999 [14,31] and 779 of 1,597 (49%) strains from 2001 to 2005 ([14,15] and ABC surveillance program, unpublished data). No serotype 19A strains from 1998 or 2000 were genotyped by MLST. Possible serotype 19A capsular switches were identified among 57 of 779 (7%) of 2001–2005 genotyped strains: ST271, 690, 899, 1092, and 1790 (n = 1 strain each); ST236, 338, and 557 (n = 2 strains each); ST230 and 292 (n = 3 strains each); ST156 (n = 4 strains); and ST695 (n = 37 strains) [14,15,31]. Among these strains, the ST69519A plus related ST89919A strains were characterised in this study. All other possible capsular switch strains were present in very low numbers and/or were also found outside the US; vaccine and nonvaccine serotype variants of ST156 are recognised in the global database (http://www.mlst.net/).

thumbnail
Table 1.

Description of the Pneumococcal Strains Characterised in This Study

https://doi.org/10.1371/journal.ppat.0030168.t001

Recipient strains were those with the background genotype into which the 19A capsular locus recombined: ST6954 (n = 3), recovered in 1999, and ST8994 (n = 1), recovered in 2002. Based upon pulsed-field gel electrophoresis, selected MLST data, and eBURST analysis, it is likely that strains of ST8994 were also common among children <5 y of age before vaccine introduction [14]. Putative donor strains of the 19A capsular locus were selected from pneumococci recovered in 1999 and 2003. The most likely donor of the 19A capsular locus was ST19919A, because ST19919A existed as a prevalent strain pre-vaccine implementation and was the major clone responsible for serotype replacement disease post-vaccine implementation [15,17]. Thus, four pre- and post-vaccine implementation strains of ST19919A were selected, plus one strain of ST64519A, which is a single-locus MLST variant of ST19919A. All serotype 19A strains collected through 2005 that were ST695 or closely related MLST genotypes were characterised as progeny strains. This resulted in the characterisation of 42 unique strains: ST69519A (n = 36), ST236519A (n = 4), ST236319A (n = 1), and ST89919A (n = 1).

Pneumococci were recovered from patients aged <1 to 96 y with invasive pneumococcal disease (Table 1). All pneumococci were recovered from blood, apart from one strain (cdc46) that was recovered from cerebrospinal fluid. Among vaccine escape strains, 40 of 42 (95%) were collected in the northeastern US (note that 37% of all 20,196 available isolates since 2000 were from the Northeast; ABC surveillance program, unpublished data); and two strains were from Minnesota (cdc27) and Colorado (cdc21). Three putative donor strains were from Georgia (cdc8, cdc9, and cdc11); all other putative donor and recipient strains were collected in the Northeast. All vaccine escape strains were recovered from 2003 to 2005. Strains from 2003 to 2004 were recovered from children ≤4 y of age, although only pediatric isolates were genotyped [15]; 2005 strains were recovered from children and adults.

Sequence Homology within Regions Flanking the Capsular Locus

Regions upstream and downstream of the capsular locus were sequenced in all 51 pneumococci (Figure 1). Sequence alignment comparisons showed that recipient, donor, and progeny strains could be grouped based upon the relatedness of the regions flanking the capsular locus, considered as one joined sequence of 20.6 kb (Figure 1). The capsular locus has not yet been sequenced, but serotypes have been confirmed phenotypically, and the serotypes of a subset of the vaccine escape strains were confirmed by a PCR assay [32]. The capsular locus is considered here only as it encodes either a serotype 4 or 19A capsule.

Progeny strains consisted of three groups, P1, P2, and P3 (Table 1). Twenty-seven ST69519A progeny strains (P1) were identical to each other in the capsular locus flanking regions and were the major group of vaccine escape strains. There were five additional P1 variants, all of which differed from the P1 strains by only one to two single nucleotide polymorphisms (SNPs) over the entire region of sequence. Furthermore, one strain, P1var6, was a single-locus MLST variant of the P1 strains at spi, but the capsular locus flanking regions of P1 and P1var6 were identical. P2 progeny strains and strain P3 were 0.3% and 1%, respectively, divergent from P1 progeny over the entire region of sequence.

Among recipient strains, the flanking sequences of cdc1, cdc3, and cdc4 (R1) were identical; recipient strain cdc7 (R2) was 0.5% divergent from R1 strains (Table 1). Among 19A donor strains there were four unique sequences: cdc2 and cdc8 (D1), cdc9 (D2), cdc11 (D3), and cdc6 (D4).

Recombinational Crossover Points and Likely Serotype 19A Donors among Progeny Strains

Sequence alignments of the capsular locus flanking regions revealed the upstream and downstream recombinational crossover points in the ST69519A strains. Three SNPs (heterologous to donor and homologous to recipient) marked and confirmed the upstream crossover point, just upstream of pbp2x; similarly, two SNPs marked and confirmed the downstream crossover point, at the 3′ end and just beyond the 3′ end of pbp1a (Figures 2A and S1), resulting in a recombinational fragment of 38.6 kb. To our knowledge, this is the first reported in vivo example of a recombinational event in pneumococcus in which the capsular locus and both adjacent PBPs recombined in what appeared to be a single event. This was shown to be experimentally possible by Trzcinski and colleagues [33].

thumbnail
Figure 2. Sequencing Results among Progeny, Donor, and Recipient Pneumococci, including Differences in pbp2x and pbp1a Genes

P1 and P1variants, and the P2 and P2variant, are each represented by one P1 and P2 schematic, respectively, for simplicity. Unique regions of sequence are indicated by variations of shading and patterns, e.g. D1, P1, D4, P2 all have the same pbp2x sequence (black), whereas R1 has a different pbp2x sequence (speckled).

(A) Schematic of the recombinational crossover points revealed in the major group of progeny strains, ST69519A. Three SNPs marked and confirmed the upstream crossover region and two SNPs marked and confirmed the downstream crossover region in the progeny strains, as shown with black arrows. The length of the recombinational fragment was 38.6 kb (10.3 kb upstream + 18.6 kb capsular locus + 9.7 kb downstream).

(B) Illustration of the second recombinational event. ST236519A is identical to ST64519A in the entire sequenced region; no obvious crossover points have yet been revealed and the recipient is unknown.

https://doi.org/10.1371/journal.ppat.0030168.g002

The likely donor of the serotype 19A capsular locus and flanking regions in ST69519A was ST19919A. Although the actual strain that donated the DNA cannot be identified with certainty, it is represented by the pre- and post-vaccination D1 strains from 1999 and 2003, respectively (Figures 2A and S1).

Sequence alignments also revealed a second recombinational event around the capsular locus in four strains of P2 progeny, all of which were first identified in 2005. Three P2 strains were identical (or one SNP different, P2var1) to D4 in the entire flanking region sequenced (Figures 2B and S1), suggesting that D4 was the likely donor of the 19A capsule among P2 progeny strains.

The third type of 19A progeny strain identified in this study was P3 (ST89919A). ST899 appears to have circulated at least since 1999 as a serotype 4 strain [14], but in 2005 appeared as a 19A serotype strain (Table 1). None of the sequenced 19A donors in this study share distinctive regions of homology to the P3 progeny; thus, the donor of this 19A capsular locus is still unknown.

Relatedness Based on MLST Genotype

Figure 3 further clarifies the evolution of these vaccine escape strains by depicting the relatedness of the background MLST genotypes. The P1 progeny (including P1 variants) are all in the cluster with ST695 as the founder genotype. ST695 is a single-locus variant of ST899, differing only at xpt (allele 113) as a result of a 39-codon deletion in the middle of the xpt gene. This truncated xpt allele was first identified in 1999 among US strains [31] and was present in ST236319A and ST22844. ST2363 and ST2284 differ from ST695 at spi and recP, respectively. The truncated xpt allele is identical to xpt allele 10 apart from the large deletion, thus allele 113 probably derived from the full-length allele 10 (Table 1; http://www.mlst.net/).

thumbnail
Figure 3. eBURST Analysis of Closely Related STs Found in This Study and in the MLST Database (http://www.mlst.net/) as of March 2007

STs described in the current study are indicated with arrows and boxes. The size of each circle is proportional to the number of isolates of that ST, the blue circle indicates the founder, ST247, of the entire clonal complex, and the yellow circles indicate subgroup founders (STs that have at least two single-locus variants). Strains from the ST247 founder complex have been recovered throughout Europe, South America, and the US (ST2365 recovered only in the US), from 1978 to 2005. The allelic profile of ST247 is aroE 16, gdh 13, gki 4, recP 5, spi 6, xpt 10, and ddl 14. Strains from the ST246 subgroup have been detected in Spain, Australia, and the United Kingdom from 1997 to 2005; strains from the ST899 subgroup have been recovered in the US and the Czech Republic from 1999 to 2002; and strains of the ST695 subgroup have been detected in the US and South Africa from 1999 to 2005. Pneumococcal strains of all STs in this diagram are serotype 4 except where indicated; however, the serotype 19A strains have only been reported in the US to date, and only since 2003. P1, P1var6, P2, and P3 refer to the type of progeny strains (see Table 1).

https://doi.org/10.1371/journal.ppat.0030168.g003

The P2 progeny (including the P2variant) appeared to evolve from the founder ST2474 clonal complex. The oldest known strains of ST2474 were identified in The Netherlands in 1978, and ST2474 and the related single-locus variants in the ST2474 cluster have been detected in other European and South American countries through at least 2005 (Figure 3; http://www.mlst.net/). The P2 progeny were first identified in 2005 in the US. ST2474 strains have not yet been characterised, thus recombinational crossover points in P2 progeny have not yet been identified. The third serotype 19A capsular change was the P3 progeny (ST89919A); Figure 3 depicts the relatedness of this ST to other STs in the clonal complex.

Sequence Diversity among PBPs

Alterations in PBP genes pbp2x, pbp1a, and pbp2b confer penicillin resistance [2426]. Allele assignments were given to each unique pbp2x and pbp1a sequence among this collection of strains; the pbp2x and pbp1a sequences from TIGR4, a serotype 4 penicillin-susceptible strain, were used for comparison and designated allele 1 (Table 1; Figure 4). All R1, R2, and P3 strains possessed pbp2x and pbp1a alleles similar to those in TIGR4, which correspond to the basal penicillin minimum inhibitory concentration (MIC) of 0.03 μg/ml that is typical of all ST6954 isolates. Among D3, D4, and P2 strains, pbp2x was altered while pbp1a was largely conserved; both pbp2x and pbp1a were altered in the D1, D2, P1, and P1variant strains. Phenotypically, strains with altered pbp2x and pbp1a genes resulted in a MIC to penicillin of 0.06–0.12 μg/ml (Table 1). The breakpoints for penicillin resistance are ≤0.06 μg/ml (susceptible), 0.12–1.0 μg/ml (intermediate resistance), and ≥2 μg/ml (high-level resistance), but the variation in susceptibility testing is ±1 doubling dilution [34]. All donor strains and all progeny strains apart from the P3 strain had MIC values at the penicillin-susceptible breakpoint; however, based on the pbp2x and pbp1a sequences, D1, D2, P1, and P1variants had divergent PBPs typical of penicillin-nonsusceptible strains ([2426]; Figure 4). All pbp sequences in this study were unique to the National Center for Biotechnology Information database.

thumbnail
Figure 4. Sequence Alignments of pbp2x and pbp1a (2,253 and 2,160 Base Pairs, Respectively, in Total Length), Showing the Polymorphic Sites (Numbered Vertically) in Each of the PBP Alleles Present in This Dataset

Note that pbp1a is reverse complemented. Allele 1 for each gene corresponds to the penicillin-susceptible TIGR4 pneumococcal strain [38]. Pbp1a_allele 5 and pbp1a_allele 5rec are identical apart from base pairs 1784 and 2009, which mark the downstream recombinational crossover point in progeny strains ST69519A.

https://doi.org/10.1371/journal.ppat.0030168.g004

Resistance to Other Antimicrobials

Among all strains, four donor and five vaccine escape strains were erythromycin resistant, possessing the mef(A) efflux phenotype [29,35]; all five vaccine escape strains and two donors (cdc8 and cdc9) were confirmed to be mef(A)-positive by PCR ([36,37]; Lesley McGee, Emory University, personal communication). The other donor strains were not tested by PCR. One erythromycin-resistant donor strain (D4) was also clindamycin resistant (Table 1). D1 and D2 donor strains were co-trimoxazole resistant, and one vaccine escape strain (cdc13) was resistant to tetracycline. One P1var2 strain (cdc44) was resistant to chloramphenicol; R2 was resistant to chloramphenicol and levofloxacin.

Discussion

To our knowledge, this is the first report of vaccine escape events in the US subsequent to national pneumococcal vaccination. The main event resulted in nonvaccine serotype, penicillin-nonsusceptible ST69519A pneumococci as a result of recombination between the recipient ST6954 and donor ST19919A (D1). A second event resulted in the emergence of ST236519A from the well-established serotype 4 clonal complex of ST247, and these strains had a different serotype 19A donor, ST64519A. No recipient has yet been identified, but recipient strains of ST247 will be characterised to elucidate the recombinational crossover points in these progeny. Finally, a third event resulted in the emergence of ST89919A, but the serotype 19A donor is not yet known.

The most plausible explanation for the emergence of these strains is that one main recombinational event, ST69519A, occurred around 2003, and the fact that it was a nonvaccine serotype variant provided the selective advantage by which these strains could evade the immune pressures invoked by the vaccine. Furthermore, with the additional selective advantage of penicillin nonsusceptibility, this progeny strain increased in frequency over the next 2 y, disseminating through the northeastern US to become the fourth most common serotype 19A genetic complex in the US (ABC surveillance program, unpublished data). The recombinational event did not appear to result in a decrease in fitness among these strains, as strains increased over time within the population. More worryingly, all of these strains were recovered from patients with bacteremia and meningitis. We cannot be absolutely certain that ST69519A strains never existed pre-vaccination, but extensive surveillance pre- and post-vaccination in the US failed to reveal any such strains [14,31], and no such strains have been reported to the MLST database from other parts of the world. Hence, even if these strains did exist pre-vaccination they were likely to be very rare, and it could still be maintained that the immune pressure resulting from PCV7 use selected for the emergence of such strains.

The two additional progeny strain types, ST236519A and ST89919A, were identified in 2005, and this may be the very early detection of these new recombinants. Whether or not these strains will increase in prevalence remains to be seen. Since serotype 19A was the major replacement disease serotype to emerge post-vaccination [15,17], our efforts to understand the genetics of this event have been focused on serotype 19A strains. This will be expanded to include other putative capsular switching events to better understand the phenomenon of capsular switching. Although it has been known for many years that capsular switching can and does occur, it is not at all clear how often it occurs within the pneumococcal population.

An exciting component of the current study is that recombination in the pneumococcus has been revealed almost as it occurred. This presents a unique opportunity to measure recombination in nature as a result of vaccine-induced changes, and may shed light on how much recombination occurs in general within the pneumococcal genome. Studies to explore these strains in detail are ongoing.

What are the implications for vaccine escape? Clearly, vaccine escape by recombination at the capsular locus has the very real potential to reduce PCV7 effectiveness in the longer term. This will almost certainly be true for any serotype-specific pneumococcal vaccine, given the diversity and complexity of serotypes possessed by the pneumococcus. These US data reinforce two key points: i) the importance of surveillance pre- and post-vaccination in countries preparing to implement PCV7, to detect changes as and when they occur; and ii) the importance of understanding the genetic events that result in vaccine escape. Discerning the genetic event is crucial to understanding vaccine escape and pneumococcal recombination in general. Such knowledge will provide guidance about the design and use of future pneumococcal vaccines.

Materials and Methods

Invasive pneumococcal strains were collected as part of the ABC surveillance program at the CDC. Microbiological testing and molecular characterisation by MLST was performed as previously described [14,31]. Strains for this study were selected from pre- and post-vaccine implementation ABC collections and sent to the University of Oxford. Sequence alignments of the TIGR4 [38], R6 [39], and Spain23F-1 (http://www.sanger.ac.uk/Projects/S_pneumoniae/) pneumococcal genomes were used to design 51 sets of PCR primers (available upon request) in regions upstream and downstream of the capsular locus (Figure 1).

DNA was extracted from pneumococcal strains using the DNeasy Tissue Kit (Qiagen UK). PCR assays were identical for all combinations of PCR primers: 2.5 μl of PCR buffer (Qiagen), 0.5 μl of dNTPs (200 μM stock), 0.5 μl of forward primer (100 μM stock), 0.5 μl of reverse primer (100 μM stock), 0.3 μl of Taq polymerase (Qiagen UK), 19.7 μl of nuclease-free water (Invitrogen UK), and 1 μl (∼1 μg) of extracted DNA. Thirty-five cycles of PCR amplification were performed: denaturation at 94 °C for 30 s, annealing at 55 °C for 30 s, and elongation at 72 °C for 60 s. Ten-microliter sequencing reactions were prepared with 4 μl primer (1 μM stock of PCR primers), 1.75 μl sequencing buffer (Applied Biosystems), 0.5 μl Big Dye (Applied Biosystems), 1.75 μl nuclease-free water, and 2 μl cleaned PCR amplicons; followed by 25 sequencing cycles of 96 °C for 10 s, 50 °C for 5 s, and 60 °C for 2 min. PCR amplicons and sequencing products were cleaned as described previously [15].

Sequences were assembled and edited using Pregap4 and Gap4 (Staden Package, http://staden.sourceforge.net/). Consensus sequences were aligned using Multalin [40] and trimmed to the longest possible consensus sequence shared by all strains. Trimmed sequences were uploaded into MEGA version 3.1 [41], and compared in an iterative process to ascertain which strains were donors of the 19A capsule and reveal recombinational crossover points in progeny strains. The eBURST algorithm was used to compare all progeny strains to all pneumococcal strains represented in the MLST database [42].

Supporting Information

Figure S1. Sequence Alignments of Major Donor, Recipient, and Progeny Strains in This Study, Aligned with the Corresponding Region of the TIGR4 Sequence

Polymorphic sites are shown across the entire sequenced region upstream and downstream of the capsular locus. Black arrows highlight the SNPs in the upstream and downstream crossover regions; the locations of pbp2x, pbp1a, and the (missing) capsular locus are also shown. Strains are labeled A and B to correspond to the strains involved in the two main recombinational events shown in Figure 2.

https://doi.org/10.1371/journal.ppat.0030168.sg001

(95 KB PPT)

Accession Numbers

All unique pbp2x and pbp1a sequences in this study were submitted to GenBank (http://www.ncbi.nlm.nih.gov/Genbank/index.html) under accession numbers EU034013–EU034023.

Acknowledgments

The authors acknowledge the laboratory expertise of Robert E. Gertz, Jr. in the Respiratory Diseases Branch (RDB) Streptococcus Laboratory at the Centers for Disease Control and Prevention (CDC), Lynne Richardson, Becky Busby, and Rory Bowden for sequencing expertise at the University of Oxford, and the Computational Biology Research Group, University of Oxford, for bioinformatics support. We are indebted to all of the hospitals, laboratories, the RDB Streptococcus Laboratory and the RDB Epidemiology section, that participate in CDC's Emerging Infections Program Network and the Active Bacterial Core surveillance core program. We are grateful to Brian G. Spratt for insightful comments on the manuscript.

Author Contributions

ABB performed the molecular work, analyzed the data, and wrote the manuscript. RP and BB were responsible for phenotypic characterisation of the strains and MLST genotyping. All authors contributed to the design of the study, discussed the results, and commented on the manuscript.

References

  1. 1. World Health Organization (2007) Pneumococcal conjugate vaccine for childhood immunization–WHO position paper. Wkly Epidemiol Rec 82: 93–104.
  2. 2. Black S, Shinefield H, Fireman B, Lewis E, Ray P, et al. (2000) Efficacy, safety and immunogenicity of heptavalent pneumococcal conjugate vaccine in children. Pediatr Infect Dis J 19: 187–195.
  3. 3. O'Brien KL, Moulton LH, Reid R, Weatherholtz R, Oski J, et al. (2003) Efficacy and safety of seven-valent conjugate pneumococcal vaccine in American Indian children: group randomised trial. Lancet 362: 355–361.
  4. 4. Whitney CG, Farley MM, Hadler J, Harrison LH, Bennett NM, et al. (2003) Decline in invasive pneumococcal disease after the introduction of protein-polysaccharide conjugate vaccine. N Engl J Med 348: 1737–1746.
  5. 5. Lexau CA, Lynfield R, Danila R, Pilishvili T, Facklam R, et al. (2005) Changing epidemiology of invasive pneumococcal disease among older adults in the era of pediatric pneumococcal conjugate vaccine. JAMA 294: 2043–2051.
  6. 6. Whitney CG, Pilishvili T, Farley MM, Schaffner W, Craig AS, et al. (2006) Effectiveness of seven-valent pneumococcal conjugate vaccine against invasive pneumococcal disease: a matched case-control study. Lancet 368: 1495–1502.
  7. 7. Bentley SD, Aanensen DM, Mavroidi A, Saunders D, Rabbinowitsch E, et al. (2006) Genetic analysis of the capsular biosynthetic locus from all 90 pneumococcal serotypes. PLoS Genet 2: e31.
  8. 8. Park IH, Pritchard DG, Cartee R, Brandao A, Brandileone MC, et al. (2007) Discovery of a new capsular serotype (6C) within serogroup 6 of Streptococcus pneumoniae. J Clin Microbiol 45: 1225–1233.
  9. 9. Crook DW, Brueggemann AB, Sleeman K, Peto TEA (2004) Pneumococcal carriage. In: Tuomanen EI, Mitchell TJ, Morrison DA, Spratt BG, editors. The pneumococcus. Washington (D.C.): ASM Press. pp. 136–147.
  10. 10. Butler JC (2004) Epidemiology of pneumococcal disease. In: Tuomanen EI, Mitchell TJ, Morrison DA, Spratt BG, editors. The pneumococcus. Washington (D.C.): ASM Press. pp. 148–168.
  11. 11. Brueggemann AB, Griffiths DT, Meats E, Peto T, Crook DW, et al. (2003) Clonal relationships between invasive and carriage Streptococcus pneumoniae and serotype- and clone-specific differences in invasive disease potential. J Infect Dis 187: 1424–1432.
  12. 12. Brueggemann AB, Peto TE, Crook DW, Butler JC, Kristinsson KG, et al. (2004) Temporal and geographic stability of the serogroup-specific invasive disease potential of Streptococcus pneumoniae in children. J Infect Dis 190: 1203–1211.
  13. 13. Enright MC, Spratt BG (1998) A multilocus sequence typing scheme for Streptococcus pneumoniae: identification of clones associated with serious invasive disease. Microbiology 144: 3049–3060.
  14. 14. Beall B, McEllistrem MC, Gertz RE, Wedel S, Boxrud DJ, et al. (2006) Pre- and postvaccination clonal compositions of invasive pneumococcal serotypes for isolates collected in the United States in 1999, 2001 and 2002. J Clin Microbiol 44: 999–1017.
  15. 15. Pai R, Moore MR, Pilishvili T, Gertz RE, Whitney CG, et al. (2005) Postvaccine genetic structure of Streptococcus pneumoniae serotype 19A from children in the United States. J Infect Dis 192: 1988–1995.
  16. 16. Byington CL, Samore MH, Stoddard GJ, Barlow S, Daly J, et al. (2005) Temporal trends of invasive disease due to Streptococcus pneumoniae among children in the intermountain west: emergence of nonvaccine serogroups. Clin Infect Dis 41: 21–29.
  17. 17. Hicks LA, Harrison LH, Flannery B, Hadler JL, Schaffner W, et al. (2007) Incidence of pneumococcal disease due to non-pneumococcal conjugate vaccine (PCV7) serotypes in the United States during the era of widespread PCV7 vaccination, 1998–2004. J Infect Dis 196: 1346–54.
  18. 18. Singleton RJ, Hennessy TW, Bulkow LR, Hammitt LL, Zulz T, et al. (2007) Invasive pneumococcal disease caused by nonvaccine serotypes among Alaska Native children with high levels of 7-valent pneumococcal vaccine coverage. JAMA 297: 1784–1792.
  19. 19. Coffey TJ, Dowson CG, Daniels M, Zhou J, Martin C, et al. (1991) Horizontal transfer of multiple penicillin-binding protein genes, and capsular biosynthetic genes, in natural populations of Streptococcus pneumoniae. Mol Microbiol 5: 2255–2260.
  20. 20. Coffey TJ, Enright MC, Daniels M, Morona JK, Morona R, et al. (1998) Recombinational exchanges at the capsular polysaccharide biosynthetic locus lead to frequent serotype changes among natural isolates of Streptococcus pneumoniae. Mol Microbiol 27: 73–83.
  21. 21. Coffey TJ, Daniels M, Enright MC, Spratt BG (1999) Serotype 14 variants of the Spanish penicillin-resistant serotype 9V clone of Streptococcus pneumoniae arose by large recombinational replacements of the cpsA-pbp1a region. Microbiology 145: 2023–2031.
  22. 22. Spratt BG, Greenwood BM (2000) Prevention of pneumococcal disease by vaccination: does serotype replacement matter? Lancet 356: 1210–1211.
  23. 23. Llull D, Muñoz R, López R, García E (1999) A single gene (tts) located outside the cap locus directs the formation of Streptococcus pneumoniae type 37 capsular polysaccharide. Type 37 pneumococci are natural, genetically binary strains. J Exp Med 190: 241–251.
  24. 24. Coffey TJ, Daniels M, McDougal LK, Dowson CG, Tenover FC, et al. (1995) Genetic analysis of clinical isolates of Streptococcus pneumoniae with high-level resistance to expanded-spectrum cephalosporins. Antimicrob Agents Chemother 39: 1306–1313.
  25. 25. Grebe T, Hakenbeck R (1996) Penicillin-binding proteins 2b and 2x of Streptococcus pneumoniae are primary resistance determinants for different classes of β-lactam antibiotics. Antimicrob Agents Chemother 40: 829–834.
  26. 26. Muñoz R, Dowson CG, Daniels M, Coffey TJ, Martin C, et al. (1992) Genetics of resistance to third-generation cephalosporins in clinical isolates of Streptococcus pneumoniae. Mol Microbiol 6: 2461–2465.
  27. 27. Hoban D, Baquero F, Reed V, Felmingham D (2005) Demographic analysis of antimicrobial resistance among Streptococcus pneumoniae: worldwide results from PROTEKT 1999–2000. Int J Infect Dis 9: 262–273.
  28. 28. Kyaw MH, Lynfield R, Schaffner W, Craig AS, Hadler J, et al. (2006) Effect of introduction of the pneumococcal conjugate vaccine on drug-resistant Streptococcus pneumoniae. N Engl J Med 354: 1455–1463.
  29. 29. Farrell DJ, Klugman KP, Pichichero M (2007) Increased antimicrobial resistance among nonvaccine serotypes of Streptococcus pneumoniae in the pediatric population after the introduction of 7-valent pneumococcal vaccine in the United States. Pediatr Infect Dis J 26: 123–128.
  30. 30. Schuchat A, Hilger T, Zell E, Farley MM, Reingold A, et al. (2001) Active Bacterial Core surveillance of the Emerging Infections Program network. Emerg Infect Dis 7: 92–99.
  31. 31. Gertz RE, McEllistrem MC, Boxrud DJ, Li Z, Sakota V, et al. (2003) Clonal distribution of invasive pneumococcal isolates from children and selected adults in the United States prior to 7-valent conjugate vaccine introduction. J Clin Microbiol 41: 4194–4216.
  32. 32. Pai R, Gertz RE, Beall B (2006) Sequential multiplex PCR approach for determining capsular serotypes of Streptococcus pneumoniae isolates. J Clin Microbiol 44: 124–131.
  33. 33. Trzcinski K, Thompson CM, Lipsitch M (2004) Single-step capsular transformation and acquisition of penicillin resistance in Streptococcus pneumoniae. J Bacteriol 186: 3447–3452.
  34. 34. Clinical and Laboratory Standards Institute (2007) Performance standards for antimicrobial susceptibility testing. 17th informational supplement, CSLI document M100-S17E. Wayne (Pennsylvania): Clinical and Laboratory Standards Institute. pp. 64–68.
  35. 35. Brueggemann AB (2006) Antibiotic resistance mechanisms among pediatric respiratory and enteric pathogens: a current update. Pediatr Infect Dis J 25: 969–973.
  36. 36. Sutcliffe J, Grebe T, Tait-Kamradt A, Wondrack L (1996) Detection of erythromycin-resistant determinants by PCR. Antimicrob Agents Chemother 40: 2562–2566.
  37. 37. Widdowson C, Klugman KP (1998) Emergence of the M-phenotype of erythromycin-resistant pneumococci in South Africa. Emerg Infect Dis 4: 277–281.
  38. 38. Tettelin H, Nelson KE, Paulsen IT, Eisen JA, Read TD, et al. (2001) Complete genome sequence of a virulent isolate of Streptococcus pneumoniae. Science 293: 498–506.
  39. 39. Hoskins J, Alborn WE, Arnold J, Blaszczak LC, Burgett S, et al. (2001) The genome of the bacterium Streptococcus pneumoniae strain R6. J Bacteriol 183: 5709–5717.
  40. 40. Corpet F (1988) Multiple sequence alignment with hierarchical clustering. Nucl Acids Res 16: 10881–10890.
  41. 41. Kumar S, Tamura K, Nei M (2004) MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform 5: 150–163.
  42. 42. Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG (2004) eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol 186: 1518–1530.