PLoS OnePLoS ONEplosplosonePLoS ONE1932-6203Public Library of ScienceSan Francisco, USA242826053839905PONE-D-13-3081810.1371/journal.pone.0081539Research ArticlePopulation Genetic Structure of Streptococcus pneumoniae in Kilifi, Kenya, Prior to the Introduction of Pneumococcal Conjugate VaccinePneumococcal Population Structure, KilifiBrueggemannAngela B.12*MurokiBeth Mbesu3KulohomaBenard W.3KaraniAngela3WanjiruEva3MorpethSusan23KamauTatu4SharifShahnaaz4ScottJ. Anthony G.23Department of Zoology, University of Oxford, Oxford, United KingdomNuffield Department of Medicine, University of Oxford, Oxford, United KingdomKEMRI-Wellcome Trust Research Programme, Kilifi, KenyaMinistry of Public Health and Sanitation, Nairobi, KenyaBeallBernardEditorCenters for Disease Control & Prevention, United States of America* E-mail: angela.brueggemann@ndm.ox.ac.uk

Competing Interests: ABB and JAGS each have vaccine-related research funding from GlaxoSmithKline Biologicals; JAGS has received conference support from Merck, Inc. Other authors declare no conflicts of interest. This does not alter the authors' adherence to all the PLOS ONE policies on sharing data and materials.

Conceived and designed the experiments: ABB JAGS. Performed the experiments: ABB BMM BWK AK EW SM. Analyzed the data: ABB BMM BWK AK SM JAGS. Contributed reagents/materials/analysis tools: TK SS. Wrote the manuscript: ABB JAGS.

201325112013811e815393172013231020132013Brueggemann et alThis is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.Background

The 10-valent pneumococcal conjugate vaccine (PCV10) was introduced in Kenya in 2011. Introduction of any PCV will perturb the existing pneumococcal population structure, thus the aim was to genotype pneumococci collected in Kilifi before PCV10.

Methods and Findings

Using multilocus sequence typing (MLST), we genotyped >1100 invasive and carriage pneumococci from children, the largest collection genotyped from a single resource-poor country and reported to date. Serotype 1 was the most common serotype causing invasive disease and was rarely detected in carriage; all serotype 1 isolates were members of clonal complex (CC) 217. There were temporal fluctuations in the major circulating sequence types (STs); and although 1-3 major serotype 1, 14 or 23F STs co-circulated annually, the two major serotype 5 STs mainly circulated independently. Major STs/CCs also included isolates of serotypes 3, 12F, 18C and 19A and each shared ≤2 MLST alleles with STs that circulate widely elsewhere. Major CCs associated with non-PCV10 serotypes were predominantly represented by carriage isolates, although serotype 19A and 12F CCs were largely invasive and a serotype 10A CC was equally represented by invasive and carriage isolates.

Conclusions

Understanding the pre-PCV10 population genetic structure in Kilifi will allow for the detection of changes in prevalence of the circulating genotypes and evidence for capsular switching post-vaccine implementation.

This project was funded by the Wellcome Trust (ref: 081835, awarded to JAGS) and the Robert Austrian Research Award in Pneumococcal Vaccinology (awarded to ABB). ABB and JAGS are funded by research fellowships from the Wellcome Trust (ref: 083511/Z/07/Z and 098532, respectively). The KEMRI-Wellcome Trust Research Programme is supported by core funding from the Wellcome Trust (ref: 092654/Z/10/A). Surveillance for pneumococcal disease in Kilifi is supported by the GAVI Alliance. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Introduction

The introduction of pneumococcal conjugate vaccine (PCV) into immunisation programmes in well-resourced countries led to a significant reduction in pneumococcal morbidity and mortality [1,2]. Children in resource-poor countries have a much higher incidence of life-threatening pneumococcal disease [3,4], thus the World Health Organisation recommended that PCV be introduced into developing countries with high childhood mortality, and the GAVI Alliance has provided support for PCV introduction [5]. Pneumococcal disease burden among young children living within the Kilifi District on the coast of Kenya is very high: the annual incidence of clinically-significant pneumococcal bacteraemia among children <5 years of age who presented to the outpatient department of Kilifi District Hospital was estimated at 436 cases per 100,000 [6]; and among children who were <1, <2 and <5 years of age and admitted to hospital the incidence was 241, 213 and 111 cases per 100,000, respectively [7]. Nasopharyngeal carriage prevalence among healthy children is also high: the overall population-based prevalence in a pre-PCV10 study was 66% (79% for children <1 year of age; 51% among children 4.5-5.0 yrs of age) [8]. Therefore, in January 2011, Kenya introduced the 10-valent PCV (PCV10), which contains serotypes 1, 4, 5, 6B, 7F, 9V, 14, 18C, 19F and 23F, into its childhood immunisation programme. PCV10 coverage was estimated to be 42% among carriage isolates and >70% among invasive pneumococci [8,9]. In January and March 2011 the Kenyan Government conducted a two-dose catch-up campaign for all children aged <5 years in Kilifi District to accelerate the conditions of a mature vaccination programme. Population-based surveillance established in 2001 at Kilifi District Hospital will be used to evaluate the impact of PCV10 introduction.

In well-resourced countries, PCV introduction led to a significant overall reduction in the incidence of invasive pneumococcal disease due to vaccine serotypes. PCV also led to a profound reduction in the prevalence of nasopharyngeal carriage of vaccine serotypes among healthy children, but with a compensatory rise in the prevalence of nonvaccine serotypes. This resulted in fundamental changes in the transmission patterns of serotypes, which had an effect on disease: the reduction in transmission of vaccine serotype pneumococci led to a herd-protection effect, benefiting older unvaccinated individuals, and the increased circulation of nonvaccine serotypes led to ‘serotype replacement disease’, which attenuated the net benefit of PCV introduction. However, many nonvaccine serotypes appear to be inherently less invasive than vaccine serotypes, thus the reduction in invasive disease caused by vaccine serotypes has, in most populations, exceeded the increase in serotype replacement disease [10-18]. Changes in the circulating serotypes also resulted in a concomitant change in the circulating genotypes in developed countries, since the serotype and multilocus sequence typing (MLST) genotype are closely associated, with known exceptions [14,19-22]. Much of the focus centred around changes in the prevalence of serotype 19A-associated genotypes, since 19A was the predominant non-vaccine serotype that increased in prevalence after PCV7 introduction, but there were changes in the prevalence of other genotypes as well [13,23,24].

In resource-poor countries, much less is known about which serotypes and genotypes are circulating at a population level, thus the impact of PCV introduction on the pneumococcal population structure and the potential for evolution in response to vaccine selective pressure cannot easily be predicted based on the experience of introducing PCV to well-resourced countries. Therefore, in this study we used MLST to provide the first large-scale characterisation of >1100 invasive and colonising pneumococcal isolates from a single resource-poor country, with the aim of revealing the population structure prior to PCV10 introduction.

MethodsEthical Statement

The project, SSC#1357, entitled “The effect of routine immunization with Pneumococcal Conjugate Vaccine in children on the strain structure of invasive and carriage isolates of S. pneumoniae in children and adults in Kilifi District”, was reviewed by the KEMRI National Ethical Review Committee and they approved the analysis of pneumococcal isolates collected though routine surveillance at Kilifi District Hospital. They also approved the analysis of carriage isolates collected in specific research studies; for these studies individual informed consent was obtained from every parent/guardian. All data in this genotyping study were analysed anonymously.

Selection of Isolates for Genotyping

Invasive isolates were recovered from the blood, cerebrospinal fluid or pleural fluid of ill children <15 years of age from 1994-2008 (total n = 628). All available unique patient isolates from 1994-2002 were included and isolates were systematically selected (e.g. every other isolate in the line listing of isolates) from years 2003-2008 to obtain 25 or 50 unique patient isolates in each year’s sample for genotyping. Three isolate samples from 2003-2008 were mixed or nonviable and were not genotyped.

Carriage isolates (n = 486) were collected during two previous studies conducted in the Kilifi District. The first study was performed in 2004 and sampled the nasopharynx of healthy persons of all age groups [25]; all available unique patient isolates recovered from children <5 yrs of age (n = 170) were included in the current genotyping study. The second study was performed from 2006-2008 and sampled healthy children <5 years of age in a rolling cross-sectional study of nasopharyngeal carriage; 320 isolates were randomly selected for genotyping from the entire collection of 1868 isolates, 4 of which were later removed as freezer stocks were nonviable [8].

In all studies, samples were cultured and pneumococci identified using standard microbiological methods. Isolates were serogrouped by latex agglutination and serotyped using the Quellung method in the Kilifi Laboratory. All invasive isolates were also serotyped by PCR using a published protocol, modified to better suit the Kilifi serotype distribution [26]. Any discrepancies that arose between the Quellung and PCR results were retested until consensus was achieved and the consensus serotypes were used in this study.

MLST, Data Confirmation and Analyses

MLST was performed according to the S. pneumoniae MLST protocol [27] and alleles and sequence types (STs) were assigned using the MLST website [19]. Unusual combinations of serotype and genotype were verified by repeating the genotyping and/or confirming the serotype. Serotype confirmation was done either by repeat testing in Kilifi or by PCR serotyping in Oxford using PCR serotyping primers and protocols adapted from previously published methods [26,28]. We were unable to PCR-amplify the gdh locus for one invasive serotype 23F isolate and the ddl locus for one carriage serotype 19F isolate despite repeated attempts and redesigned primers (presumably due to divergent sequence in the primer binding regions), so these two isolates were removed from genotyping analyses. The data were compiled and analysed using Microsoft Excel and STATA v. 11. STs were clustered into clonal complexes (CCs) using Phyloviz [29] with the following settings: Dataset type, Multi-locus Sequence Typing; Distance, eBURST Distance; and Level, SLV. The Kilifi dataset was combined with the entire MLST database (as of February 2012; total combined n = 16,070 isolates) for the Phyloviz analyses. The entire Kilifi dataset was also submitted to the MLST database.

ResultsSerotype and Sequence Type Distributions and Identification of Clonal Complexes

1114 invasive and carriage isolates were recovered from children in Kilifi from 1994-2008 and genotyped by MLST (Table 1). Isolate representatives of all PCV10 serotypes were found; 8 of the 10 most abundant serotypes in the invasive collection were PCV10 serotypes, as were 5 of the 10 most common serotypes in the carriage collection. Although serotype 7F was a minor serotype in this collection (n = 3 isolates), all other PCV10 serotypes were well represented (29-161 isolates characterised per serotype). Vaccine serotypes 1, 6B, 14 and 19F were each represented by >100 isolates.

10.1371/journal.pone.0081539.t001Serotype distributions among the invasive (1994-2008) and carriage (2004, 2006-2008) collections of pneumococci that were genotyped by MLST.
AllInvasive isolates
Carriage isolates
Serotypeaisolates0-1y2-4y5-14yAll ages0-1y2-4yAll ages
1161333788158033
6B1073713656252651
19F104174425384179
1410351121376151227
6A93287641252752
23F86338445192241
5453051045000
18C39197632167
9V3112221651015
429841123246
19A28125219729
35B26600691120
15B/C23210371320
32081615325
11A20210371017
10A1960396410
23B18200231316
12F1652714112
15A1623057411
131520026713
34113003358
16F92013336
7C90101448
19B80101347
2070000167
2161001145
24F72114033
23A50000235
35A50000235
2942013101
3841102112
10F42103011
17F40000224
10B30000033
33B31001112
7F33003000
9L30022011
820000112
3720000112
18F22002000
28F21001101
33D20000202
35F20000022
211001000
1810000011
3110000011
15F10000011
18B11001000
22A11001000
9N10000101
nontypeable10000101
Total1114338117173628217269486

Serotypes in boldface font are included in PCV10. Invasive isolates were serotyped by both Quellung and PCR; discordant results were generally at the level of serotype (e.g. 6A vs. 6B). Testing was repeated until consensus was achieved and the consensus serotypes were used in this study.

There were 21 major (>15 isolate representatives) CCs identified among the entire collection of pneumococci, 13 of which were of PCV10 serotypes and 8 of which were found among non-PCV10 serotypes (Tables 2 and 3). An initial observation was that the STs and CCs circulating in Kilifi were largely different to those found circulating among most developed countries (as reported to the MLST database [19]). We submitted 40 novel alleles and 212 novel STs to the MLST database; the novel STs described 26.3% of the invasive isolates and 30.9% of the carriage isolates. The clonal complexes and sequence types found among the invasive and carriage collections, stratified by serotype, are described in Tables S1 and S2, respectively.

10.1371/journal.pone.0081539.t002Major clonal complexes<sup>a</sup> identified among isolates of PCV10 serotypes.
CCSTTotal  Invasive  Carriage  Predominant serotypeOther serotype(s)
2171--16115831 (100%)-
21710610421-
613272611-
614262601-
Other2201-
2464--242044 (95.8%)2
853201734-
Other4314 (75.0%)2
2895--444405 (100%)-
245262605-
289141405-
Other4405-
7069V--175129V (100%)-
7067259V-
52836249V-
Other4139V-
138118C--3528718C (100%)-
13813124718C-
Other44018C (75.0%)18B
8546B--3013176B (86.7%)6A
8542813156B (85.7%)6A
Other2026B-
27136B--241776B (91.7%)6A, 23F
271310736B-
53026606B-
Other8446B (75.0%)6A, 23F
6314--61421914 (100%)-
84244271714-
271655014-
Other1210214-
23014,3--4031914 (67.5%)3
2302420414-
70011833-
Other53214 (60.0%)3
84419F--50104019F (92.0%)11A, 14, 23A, 23F
8441961319F (94.7%)14
53391521319F (80.0%) 11A, 23A, 23F
53675--519F-
Other112919F-
271519F--2191219F (95.2%)14
2715114719F-
608864219F (83.3%)14
Other41319F-
271423F--43232023F (100%)-
271439221723F-
Other41323F-
98823F--2516923F (88.0%) 4, 14, 19F, 35B
9881710723F (94.1%)19F
Other86223F (62.5%) 4, 14, 35B

Note: CC = clonal complex; ST = sequence type; ‘-’ indicates that no other serotypes were detected.

CCs with >15 isolate representatives shown here; STs within a CC that have fewer than 5 isolate representatives are grouped as “Other”.

10.1371/journal.pone.0081539.t003Major clonal complexes<sup>a</sup> identified among isolates of non-PCV10 serotypes.
CCSTTotal  Invasive  Carriage  Predominant serotype(s)Other serotypes
70113,15BC--3152613 (48.4%), 15BC (16.1%) 6B, 9V, 23F
7011851313 (61.1%), 6B (16.7%) 9V (11.1%), 15BC (11.1%)
53406--615BC-
Other70713 (57.1%), 15BC (28.6%)23F
590215A--2781915A (68.8%)9V, 11A, 15BC, 18C
8408359V (37.5%)35A (25.0%), 11A, 15BC, 35B
533673415A (71.4%) 15BC, 18C
Other1221015A (50.0%), 35A (25.0%)11A (16.7%), 9V
84719A--2618819A (100%)-
8472013719A-
Other65119A-
4996A--257186A (88.0%)6A/B, 15BC
499216156A (95.2%) 15BC
Other4136A (75.0%), 6B (25.0%)-
114635B--2552035B (84.0%) 17F, 23F, 29
11461641235B (87.5%) 23F, 29
Other91835B (77.8%)17F
85210A--1991010A (89.5%)10B
8521881010A (88.9%)10B
530411010A-
9146A--183156A (88.9%) 6B
5354132116A (84.6%)6B
Other5146A-
98912F--1614212F (100%)-
98999012F-
535264212F-
579711012F-

Note: CC = clonal complex; ST = sequence type; ‘-’ indicates that no other serotypes were detected.

CCs with >15 isolate representatives shown here; STs within a CC that have fewer than 5 isolate representatives are grouped as “Other”.

Population Genetic Structure Among PCV10 Serotypes

161 serotype 1 isolates from Kilifi were genotyped and all were members of CC2171 (note: STserotype; Table 2). Only 3 of the serotype 1 isolates were recovered from carriage, the rest were isolated from patients with invasive disease. It was previously reported that serotype 1 genotypes had a distinctive phylogeography [30]. Ancestral ST2171 and its closely related single locus variants (SLVs) were shown to be the predominant serotype 1 STs in Kenya and Israel, and later were also shown to be common in other African countries, namely Ghana, Burkina Faso and The Gambia [31-33]. Other serotype 1 genotypes appear to predominate elsewhere [19,30]. Based on available data CC2171 still appears to be largely an African CC [19].

Vaccine serotype 4 isolates from Kilifi were predominantly invasive isolates of ST853, which are members of the widely distributed CC2464 (Table 2). CC2464 has been of recent interest since unique recombinants (a serotype 4 to 19A capsular switch) in the United States emerged from this CC [34-36]. Only 1 isolate of CC246 in Kilifi expressed an alternative serotype and this was serotype 2, a member of a cluster of serotype 2 isolates within CC246 that have also previously been reported in Asia and Africa [19].

All of the vaccine serotype 5 isolates (n = 44) were from cases of invasive disease, and all were members of CC2895. ST289 is the widely-distributed Pneumococcal Molecular Epidemiology Network (PMEN) clone represented by the Columbia5-19 strain [37] and although it was a common ST (n = 14 isolates), ST2455, an SLV of ST2895, was more prevalent in Kilifi (n = 26 isolates).

One major serotype 9V CC was identified (CC7069V) and 100% of the CC706 isolates expressed serotype 9V. 39 vaccine serotype 18C isolates were characterised in this study and 35 of these were members of CC138118C. ST138118C shares only the xpt MLST allele with ST11318C (PMEN clone Netherlands18C-36), which is a widely-distributed genotype associated with serotype 18C [19].

Among vaccine serotypes 6B, 14, 19F and 23F there were two major CCs circulating for each serotype and in each case the two CCs defined a large proportion of the isolate representatives for each serotype: 45% of 107 serotype 6B isolates; 85% of 103 serotype 14 isolates; 64% of 104 serotype 19F isolates; and 76% of 86 serotype 23F isolates (Table 2). One notable observation was that all 24 Kilifi isolates of ST230 expressed serotype 14 and 83% of them were causing invasive disease. CC230 is a widely-distributed CC that is predominantly serotypes 14 and 19A, although ST230 (PMEN clone Denmark14-32) isolates have been reported with alternative serotypes [19]. ST23019A increased in prevalence in the United States (US) and Spain after PCV7 implementation [35,36,38]. In Kilifi, a double locus variant (DLV) of ST230 – ST 700 – was also prevalent, isolates of which exclusively expressed serotype 3. ST7003 was the predominant genotype associated with serotype 3 in this Kilifi collection and ST7003 isolates have been found elsewhere in Africa. The common globally-distributed serotype 3 genotype is ST1803 (PMEN clone Netherlands3-31), and it is only distantly related to ST7003 (MLST alleles gki and spi are shared) [19].

Population Genetic Structure Among Non-PCV10 Serotypes

CCs 70113,15BC, 4996A, 114635B, 85210A, 9146A and 590215A were major clones found in Kilifi, predominantly comprised of isolates from carriage rather than invasive disease (Table 3). CC84719A (n = 26 isolates; 100% serotype 19A) is only distantly related (STs share ≤2 MLST alleles) to any of the other major serotype 19A genotypes (e.g. CCs 199, 320/271/236, 695) that circulate in many countries, which have significantly increased post-PCV introduction in the US and elsewhere [35,36,39,40].

All serotype 12F isolates (n = 16) from Kilifi were members of CC98912F and isolates of this CC have also been found elsewhere in Africa [19]. The predicted ancestral ST, ST98912F, shares only the spi MLST allele with the widely-distributed ST21812F (PMEN clone Denmark12F-34). Serotype 12F is generally uncommon among carriage isolates, which is also true of this Kilifi dataset (14 of 16 serotype 12F isolates were invasive), and has historically been shown to cause outbreaks [41,42], but it is not included in any currently available PCV.

Temporal Variation of CCs Containing Invasive Isolates Collected From 1994-2008.

Since the invasive isolate dataset in Kilifi spanned 15 years, we were able to explore temporal variation among circulating STs and CCs over that extended period of time. 12 CCs were comprised of >10 invasive isolate representatives in total over the time period, which represented 68% (424/628) of the entire invasive collection (Figure 1). The prevalence of every CC varied over time, to a greater or lesser extent. Most obvious was the dominance of CC2171 (n = 158), which ranged in prevalence from 8 - 52% of all invasive isolates characterised each year. The next most prevalent CC was CC2895 (n = 44), which fluctuated in prevalence from 0 - 22% over the period of surveillance. Other major CCs also varied in prevalence over the surveillance years. These longitudinal data will be particularly important in analyses that monitor changes in the prevalence of individual CCs, particularly CCs of nonvaccine serotypes such as CC230 (ST7003), CC98912F and CC84719A, which ranged from 0-4%, 0-9% and 0-13% of all invasive isolates characterised each year, respectively, prior to PCV10 introduction. Note that CCs not depicted in Figure 1 made only minor contributions to the overall pre-PCV population structure and contributed 0-3 isolates in any given year (data not shown), but these data will be re-evaluated in future post-PCV10 analyses.

10.1371/journal.pone.0081539.g001Prevalence of major clonal complexes (CCs) associated with invasive isolates collected in Kilifi.

The total number of isolates genotyped each year from 1994-2008 is stated in parentheses in the y-axis labels.

Figures 2 and 3 depict temporal variation in the prevalence of several of the major STs circulating among children with invasive disease in Kilifi over the 15 year surveillance period. These STs were selected for analysis because there were multiple major STs (≥10 isolates) that expressed the same serotype. Serotype 1 isolates were recovered every year of surveillance and the three major STs associated with serotype 1 are all closely related: ST217 is the predicted ancestor and ST613 and ST614 are SLVs of ST217. Either two or three of these STs were represented in 10 of the 15 surveillance years, and in the remaining five years only one ST was detected (Figure 2A). Figure 2B depicts the two major serotype 5 STs (which are SLVs of each other), but despite their genetic similarity at an ST level these STs were more restricted with respect to when they were detected, i.e. ST2895 was predominant in the early years of surveillance, but was replaced by ST2455 in the latter years. In contrast, Figures 3A and 3B depict the prevalence of pairs of unrelated serotype 14 and 23F genotypes, respectively, in each calendar year. The two serotype 14 STs were circulating concomitantly in half of the surveillance years, and both serotype 23F STs were detected in four of the surveillance years (note that only 10 isolates of ST98823F were detected overall).

10.1371/journal.pone.0081539.g002Fluctuation in the prevalence of major sequence types (STs) associated with invasive isolates from Kilifi.

Invasive isolates were collected from 1994-2008 and the major STs associated with serotypes 1 (Figure 2A) and 5 (Figure 2B) are depicted in each panel. The total number of isolates genotyped each year (see Methods) is stated in parentheses in the x-axis labels.

10.1371/journal.pone.0081539.g003Fluctuation in the prevalence of major sequence types (STs) associated with invasive isolates from Kilifi.

Invasive isolates were collected from 1994-2008 and the major STs associated with serotypes 14 (Figure 3A) and 23F (Figure 3B) are depicted in each panel. The total number of isolates genotyped each year (see Methods) is stated in parentheses in the x-axis labels.

Figure 4 completes the picture for the circulating major STs, demonstrating that although only one major genotype predominated for each of these serotypes, none of these major genotypes circulated every year throughout the surveillance period. Although a frequency of ≥10 isolates was used as the criterion for “major” ST and only nine isolates of ST98912F were recovered over the study period, the ST98912F data are shown to demonstrate that there was no serotype 12F epidemic over this surveillance period.

10.1371/journal.pone.0081539.g004Five additional major sequence types (STs) associated with invasive isolates that circulated in Kilifi.

Major STs associated with serotypes 12F, 19A, 6B, 4 and 18C are depicted for each surveillance year from 1994-2008. The total number of isolates genotyped each year (see Methods) is stated in parentheses in the x-axis labels.

Discussion

The introduction of PCV10 is vitally important for Kenya since the morbidity and mortality associated with pneumococcal disease is high. The incidence of invasive disease associated with PCV10 serotypes has decreased since vaccination began [43], with a marked attenuation of morbidity and mortality among vaccinated children, and potentially also among unvaccinated older children and adults if herd protection proves to be sufficient. It is expected that introduction of PCV10 in Kenya will perturb the pneumococcal population genetic structure, but a key question is whether or not these perturbations will diminish the overall benefit of vaccine introduction by resulting in increased nonvaccine serotype disease. By comparison, the US has used PCV for the longest period of time and the most notable nonvaccine serotype increase has been that of serotype 19A disease, which was in part due to the emergence of a novel genotype [34-36]. Importantly, this new recombinant (ST69519A) quickly became the third most common serotype 19A CC causing invasive disease in the US [35], which emphasises the importance of detecting such population-based genetic changes.

Post-vaccine changes in the US were revealed because the pre-vaccine pneumococcal population genetic structure in the US was well defined. Active, population- and laboratory-based invasive disease surveillance has been on-going in the US since 1995 (5 years before PCV was introduced) and two studies specifically characterised large, representative collections of pneumococci to identify the pre-PCV7 baseline set of genotypes circulating in the US [14,20]. The Kilifi pneumococcal invasive disease surveillance programme is also active, population- and laboratory-based, and collected isolates and data for 15 years prior to PCV introduction. Several carefully designed pneumococcal carriage studies have also been performed during that time [8,25]. The data (patient demographics, clinical outcome, serotype) combined with the genotyping work described here provide a comprehensive description of the pre-vaccine pneumococcal population in Kilifi.

MLST has been used to characterise many thousands of pneumococci collected across Europe, North and South America, Australia and Asia and many of the major STs within countries have disseminated across continents, e.g. STs 8123F,19F, 906B, 1569V, 914, 12414, 11318C, 21812F, 1917F, 1803, 19919A [19]. Many such STs have been identified as PMEN clones, in part defined by their widespread nature [19,37]. Our study demonstrated that although the major serotypes found in Kilifi and Africa are similar to those found elsewhere in the world, most of the major STs/CCs are different, at least as far as we can tell from the data in the MLST database. It is important to note that the contribution of data to the MLST database is voluntary – investigators must submit new alleles and new STs to the database for assignment, but only rarely do investigators submit their entire study dataset. Once an ST is added to the database, there is no obligation for another investigator who subsequently detects that ST elsewhere to indicate this to the curators. Thus, the MLST database reliably captures allelic diversity, but we do not have all the information about all isolates genotyped to know for certain whether some of the STs detected in Kilifi are truly African in origin, even though that is what the database suggests. However, with those caveats, the MLST database currently contains over 9,000 STs and more than 21,000 isolates recovered from all over the globe, and we do know that Africa is generally under-represented in the MLST database as compared to other geographical regions. Therefore these data significantly increase our understanding of the pneumococcal population structure in Kilifi, and possibly in Africa more broadly. However, since many of the major STs identified in Kilifi were different to the major STs of the same serotypes that circulate globally, we have little additional knowledge on which to predict their likely increase or evolution and therefore follow-up studies will be essential.

An interesting observation was the temporal fluctuation in STs that expressed the same serotypes, accepting that for some of these serotypes the numbers observed each year were small. There were no major changes in Kilifi District (e.g. changes in laboratory methods or practice, antimicrobial use or stewardship, early vaccine uptake, etc.) over the surveillance period that would have contributed to big fluctuations in the circulating serotypes or genotypes. The prevalence of circulating serotypes is known to vary and thus temporal fluctuations in circulating serotypes and genotypes in Kilifi were not surprising [44-47]; however, our data might suggest that for serotypes 1, 14 and 23F it mattered less which of the different STs associated with these serotypes was circulating. Previously published papers have debated whether it is the serotype or MLST genotype that plays a more important role in the potential for an isolate to cause invasive disease [22,48-50]. The argument favouring a primary role of serotype is supported by the fact that in Kilifi different major genotypes expressing the same serotype were detected concomitantly. All three serotype 1 STs were closely related, but the pairs of serotype 14 and 23F STs were unrelated. Alternatively, the two serotype 5 STs, although closely related at the MLST loci, appeared to be more restricted in their circulation. Perhaps the two serotype 5 STs differ markedly elsewhere in the genome, and as a result one serotype 5 genotype can outcompete the other and make co-circulation less likely. It is known that the immunological response differs for each pneumococcal serotype [51], and thus it may be the serotype-specific immune response that largely determines the circulation of serotypes regardless of which genotypic backbone they maintain. More likely, it is a particular combination of serotype (and corresponding immunity within the human population) and genotype that is the best explanation for the temporal patterns of ST circulation we observed in Kilifi.

Several major CCs were comprised of nonvaccine serotype isolates (10A, 13, 15A, 15BC and 35B) primarily recovered from healthy children, although they were also recovered from children with invasive disease so the potential for serotype replacement disease remains. Serotypes 10A, 15BC and 35B increased in prevalence in the US post-PCV7 vaccine [14,15,23,52], but the predominant STs were different to those found in Kilifi so the predictive power based on genotype is minimal [14,23,24]. Nonvaccine serotypes 3, 12F and 19A were already prevalent prior to vaccine introduction and their ability to cause invasive disease in Kilifi and elsewhere is clear. Serotypes 3 and 19A significantly increased in prevalence post-vaccine implementation in the US, although serotype 12F significantly decreased [14,15]. Serotype 19A has also significantly increased elsewhere [40,53].

In this study we established the baseline set of genotypes in Kilifi prior to PCV10 introduction, which will allow for the detection of changes in prevalence of pre-existing STs and the identification of new nonvaccine STs (as putative imports or new recombinants). It will be essential that any perceived increases or decreases in serotype or genotype prevalence after PCV10 vaccination is established in Kilifi be considered in the context of the pre-vaccine genotypic landscape.

Supporting Information

Clonal complexes and sequence types found in the invasive pneumococcal collection, stratified by serotype.

(DOCX)

Click here for additional data file.

Clonal complexes and sequence types found in the carriage pneumococcal collection, stratified by serotype.

(DOCX)

Click here for additional data file.

The authors acknowledge the use of the pneumococcal MLST database, which is located at Imperial College London and is funded by the Wellcome Trust; in particular, we thank Cynthia Bishop for her effort in entering the many new alleles and sequence types that were identified in this study. This paper is published with the permission of the Director of the Kenya Medical Research Institute. ABB and JAGS are both Wellcome Trust Research Fellows. Co-author Eva Wanjiru tragically passed away in 2012 shortly after delivering her baby. She was an important part of the phenotyping team in this study and with the permission of her husband is included in this paper as a posthumous author.

ReferencesFlascheS, Van HoekAJ, SheasbyE, WaightP, AndrewsNet al. (2011) Effect of pneumococcal conjugate vaccination on serotype-specific carriage and invasive disease in England: a cross-sectional study. PLoS Med8: e1001017 PubMed: 2148371821483718Centers for Disease Control and Prevention (CDC) (2005) Direct and indirect effects of routine vaccination of children with 7-valent pneumococcal conjugate vaccine on incidence of invasive pneumococcal disease--United States, 1998-2003. MMWR Morb Mortal Wkly Rep54: 893-897 PubMed: 1616326216163262O'BrienKL, WolfsonLJ, WattJP, HenkleE, Deloria-KnollMet al. (2009) Burden of disease caused by Streptococcus pneumoniae in children younger than 5 years: global estimates. Lancet374: 893-902.10.1016/S0140-6736(09)61204-6 PubMed: 1974839819748398BlackRE, CousensS, JohnsonHL, LawnJE, RudanIet al. (2010) Global, regional, and national causes of child mortality in 2008: a systematic analysis. Lancet375: 1969-1987.10.1016/S0140-6736(10)60549-1 PubMed: 2046641920466419Centers for Disease Control and Prevention (CDC) (2013) Progress in introduction of pneumococcal conjugate vaccine - worldwide, 2000-2012. MMWR Morb Mortal Wkly Rep62: 308-311 PubMed: 2361567423615674BrentAJ, AhmedI, NdirituM, LewaP, NgetsaCet al. (2006) Incidence of clinically significant bacteraemia in children who present to hospital in Kenya: community-based observational study. Lancet367: 482-488.10.1016/S0140-6736(06)68180-4 PubMed: 1647312516473125BerkleyJA, LoweBS, MwangiI, WilliamsT, BauniEet al. (2005) Bacteremia among children admitted to a rural hospital in Kenya. N Engl J Med352: 39-47.10.1056/NEJMoa040275 PubMed: 1563511115635111AbdullahiO, KaraniA, TigoiCC, MugoD, KunguSet al. (2012) The prevalence and risk factors for pneumococcal colonization of the nasopharynx among children in Kilifi District, Kenya. PLOS ONE7: e30787.10.1371/journal.pone.0030787 PubMed: 2236348922363489MoïsiJC, KaraniA, KosiV, LewaP, MwarumbaSet al. (2012) Secular trends in serotypes causing invasive pneumococcal disease in Kilifi, Kenya, 1994-2010. 8th International Symposium on Pneumococci and Pneumococcal Diseases, Brazil, Poster no. 101De WalsP, RobinE, FortinE, ThibeaultR, OuakkiMet al. (2008) Pneumonia after implementation of the pneumococcal conjugate vaccine program in the province of Quebec, Canada. Pediatr Infect Dis J27: 963-968.10.1097/INF.0b013e31817cf76f PubMed: 1884598218845982RodenburgGD, de GreeffSC, JansenAG, de MelkerHE, SchoulsLMet al. (2010) Effects of pneumococcal conjugate vaccine 2 years after its introduction, The Netherlands. Emerg Infect Dis16: 816-823.10.3201/eid1605.091223 PubMed: 2040937220409372HuangSS, HinrichsenVL, StevensonAE, Rifas-ShimanSL, KleinmanKet al. (2009) Continued impact of pneumococcal conjugate vaccine on carriage in young children. Pediatrics124: e1-11.10.1542/peds.2008-3099 PubMed: 1956425419564254VestrheimDF, HøibyEA, AabergeIS, CaugantDA (2010) Impact of a pneumococcal conjugate vaccination program on carriage among children in Norway. Clin Vaccine Immunol17: 325-334.10.1128/CVI.00435-09 PubMed: 2010700620107006BeallB, McEllistremMC, GertzRE, WedelS, BoxrudDJet 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 Microbiol44: 999-1017.10.1128/JCM.44.3.999-1017.2006 PubMed: 1651788916517889HicksLA, HarrisonLH, FlanneryB, HadlerJL, SchaffnerWet 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 Dis196: 1346-1354.10.1086/521626 PubMed: 1792239917922399WhitneyCG, FarleyMM, HadlerJ, HarrisonLH, BennettNMet al. (2003) Decline in invasive pneumococcal disease after the introduction of protein-polysaccharide conjugate vaccine. N Engl J Med348: 1737-1746.10.1056/NEJMoa022823 PubMed: 1272447912724479(2010) Changing epidemiology of pneumococcal serotypes after introduction of conjugate vaccine: July 2010 report. Wkly Epidemiol Rec85: 434-436 PubMed: 2103871621038716ScottJR, MillarEV, LipsitchM, MoultonLH, WeatherholtzRet al. (2012) Impact of more than a decade of pneumococcal conjugate vaccine use on carriage and invasive potential in Native American communities. J Infect Dis205: 280-288.10.1093/infdis/jir730 PubMed: 2212831522128315Streptococcuspneumoniae MLST website. Available: http://spneumoniae.mlst.net/. Accessed 2013 Oct 26GertzRE, McEllistremMC, BoxrudDJ, LiZ, SakotaVet 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 Microbiol41: 4194-4216.10.1128/JCM.41.9.4194-4216.2003 PubMed: 1295824712958247Active Bacterial Core surveillance (ABCs) website. Available: http://www.cdc.gov/abcs/index.html. Accessed 2013 Oct 26BrueggemannAB, GriffithsDT, MeatsE, PetoT, CrookDWet al. (2003) Clonal relationships between invasive and carriage Streptococcus pneumoniae and serotype- and clone-specific differences in invasive disease potential. J Infect Dis187: 1424-1432.10.1086/374624 PubMed: 1271762412717624GonzalezBE, HultenKG, LamberthL, KaplanSL, MasonEOJr. (2006) Streptococcus pneumoniae serogroups 15 and 33: an increasing cause of pneumococcal infections in children in the United States after the introduction of the pneumococcal 7-valent conjugate vaccine. Pediatr Infect Dis J25: 301-305.10.1097/01.inf.0000207484.52850.38 PubMed: 1656798016567980GertzRE, LiZ, PimentaFC, JacksonD, JuniBAet al. (2010) Increased penicillin nonsusceptibility of nonvaccine-serotype invasive pneumococci other than serotypes 19A and 6A in post-7-valent conjugate vaccine era. J Infect Dis201: 770-775.10.1086/650496 PubMed: 2017813920178139AbdullahiO, NyiroJ, LewaP, SlackM, ScottJA (2008) The descriptive epidemiology of Streptococcus pneumoniae and Haemophilus influenzae nasopharyngeal carriage in children and adults in Kilifi district, Kenya. Pediatr Infect Dis J27: 59-64.10.1097/INF.0b013e31814da70c PubMed: 1816294018162940PaiR, GertzRE, BeallB (2006) Sequential multiplex PCR approach for determining capsular serotypes of Streptococcus pneumoniae isolates. J Clin Microbiol44: 124-131.10.1128/JCM.44.1.124-131.2006 PubMed: 1639095916390959EnrightMC, SprattBG (1998) A multilocus sequence typing scheme for Streptococcus pneumoniae: identification of clones associated with serious invasive disease. Microbiology144: 3049-3060.10.1099/00221287-144-11-3049 PubMed: 98467409846740SahaSK, DarmstadtGL, BaquiAH, HossainB, IslamMet al. (2008) Identification of serotype in culture negative pneumococcal meningitis using sequential multiplex PCR: implication for surveillance and vaccine design. PLOS ONE3: e3576.10.1371/journal.pone.0003576 PubMed: 1897488718974887FranciscoAP, VazC, MonteiroPT, Melo-CristinoJ, RamirezMet al. (2012) PHYLOViZ: Phylogenetic inference and data visualization for sequence based typing methods. BMC Bioinformatics13: 87.10.1186/1471-2105-13-87 PubMed: 2256882122568821BrueggemannAB, SprattBG (2003) Geographic distribution and clonal diversity of Streptococcus pneumoniae serotype 1 isolates. J Clin Microbiol41: 4966-4970.10.1128/JCM.41.11.4966-4970.2003 PubMed: 1460512514605125LeimkugelJ, Adams ForgorA, GagneuxS, PflügerV, FlierlCet al. (2005) An outbreak of serotype 1 Streptococcus pneumoniae meningitis in northern Ghana with features that are characteristic of Neisseria meningitidis meningitis epidemics. J Infect Dis192: 192-199.10.1086/431151 PubMed: 1596221315962213AntonioM, HakeemI, AwineT, SeckaO, SankarehKet al. (2008) Seasonality and outbreak of a predominant Streptococcus pneumoniae serotype 1 clone from The Gambia: expansion of ST217 hypervirulent clonal complex in West Africa. BMC Microbiol8: 198.10.1186/1471-2180-8-198 PubMed: 1901461319014613YaroS, LourdM, TraoréY, Njanpop-LafourcadeBM, SawadogoAet al. (2006) Epidemiological and molecular characteristics of a highly lethal pneumococcal meningitis epidemic in Burkina Faso. Clin Infect Dis43: 693-700.10.1086/506940 PubMed: 1691294116912941BrueggemannAB, PaiR, CrookDW, BeallB (2007) Vaccine escape recombinants emerge after pneumococcal vaccination in the United States. PLoS Pathog3: e168.10.1371/journal.ppat.0030168 PubMed: 1802070218020702BeallBW, GertzRE, HulkowerRL, WhitneyCG, MooreMRet al. (2011) Shifting genetic structure of invasive serotype 19A pneumococci in the United States. J Infect Dis203: 1360-1368.10.1093/infdis/jir052 PubMed: 2139839521398395PaiR, MooreMR, PilishviliT, GertzRE, WhitneyCGet al. (2005) Postvaccine genetic structure of Streptococcus pneumoniae serotype 19A from children in the United States. J Infect Dis192: 1988-1995.10.1086/498043 PubMed: 1626777216267772McGeeL, McDougalL, ZhouJ, SprattBG, TenoverFCet al. (2001) Nomenclature of major antimicrobial-resistant clones of Streptococcus pneumoniae defined by the pneumococcal molecular epidemiology network. J Clin Microbiol39: 2565-2571.10.1128/JCM.39.7.2565-2571.2001 PubMed: 1142756911427569ArdanuyC, TubauF, PallaresR, CalatayudL, DomínguezMAet al. (2009) Epidemiology of invasive pneumococcal disease among adult patients in Barcelona before and after pediatric 7-valent pneumococcal conjugate vaccine introduction, 1997-2007. Clin Infect Dis48: 57-64.10.1086/597108 PubMed: 1903577919035779MooreMR, GertzRE, WoodburyRL, Barkocy-GallagherGA, SchaffnerWet al. (2008) Population snapshot of emergent Streptococcus pneumoniae serotype 19A in the United States, 2005. J Infect Dis197: 1016-1027.10.1086/528996 PubMed: 1841953918419539ReinertR, JacobsMR, KaplanSL (2010) Pneumococcal disease caused by serotype 19A: review of the literature and implications for future vaccine development. Vaccine28: 4249-4259.10.1016/j.vaccine.2010.04.020 PubMed: 2041626620416266HodgesRG, MacLC (1946) Epidemic pneumococcal pneumonia; description of the epidemic. Am J Hyg44: 183-192 PubMed: 2100301521003015RobinsonDA, TurnerJS, FacklamRR, ParkinsonAJ, BreimanRFet al. (1999) Molecular characterization of a globally distributed lineage of serotype 12F Streptococcus pneumoniae causing invasive disease. J Infect Dis179: 414-422.10.1086/314589 PubMed: 98780269878026Trust Kemri-Wellcome The Pneumococcal Conjugate Vaccine Impact Study (PCVIS) website. Available: http://www.kemri-wellcome.org/index.php/en/studies_inner/75. Accessed 2013 October 26FinlandM, BarnesMW (1977) Changes in occurrence of capsular serotypes of Streptococcus pneumoniae at Boston City Hospital during selected years between 1935 and 1974. J Clin Microbiol5: 154-166 PubMed: 14971 14971LodaFA, CollierAM, GlezenWP, StrangertK, ClydeWA Jret al. (1975) Occurrence of Diplococcus pneumoniae in the upper respiratory tract of children. J Pediatr87: 1087-1093.10.1016/S0022-3476(75)80120-X PubMed: 241779241779JefferiesJM, SmithAJ, EdwardsGFS, McMenaminJ, MitchellTJet al. (2010) Temporal analysis of invasive pneumococcal clones from Scotland illustrates fluctuations in diversity of serotype and genotype in the absence of pneumococcal conjugate vaccine. J Clin Microbiol48: 87-96.10.1128/JCM.01485-09 PubMed: 1992348819923488HarboeZB, BenfieldTL, Valentiner-BranthP, HjulerT, LambertsenLet al. (2010) Temporal trends in invasive pneumococcal disease and pneumococcal serotypes over 7 decades. Clin Infect Dis50: 329-337.10.1086/649872 PubMed: 2004747820047478MelinM, TrzcińskiK, MeriS, KäyhtyH, VäkeväinenM (2010) The capsular serotype of Streptococcus pneumoniae is more important than the genetic background for resistance to complement. Infect Immun78: 5262-5270.10.1128/IAI.00740-10 PubMed: 2085551320855513SandgrenA, SjöströmK, LiljequistBO, ChristenssonB, SamuelssonAet al. (2004) Effect of clonal and serotype-specific properties on the invasive capacity of Streptococcus pneumoniae. J Infect Dis189: 785-796.10.1086/381686 PubMed: 1497659414976594HyamsC, TrzcinskiK, CamberleinE, WeinbergerDM, ChimalapatiSet al. (2013) Streptococcus pneumoniae capsular serotype invasiveness correlates with the degree of factor H binding and opsonization with C3b/iC3b. Infect Immun81: 354-363.10.1128/IAI.00862-12 PubMed: 2314703823147038RennelsMB, EdwardsKM, KeyserlingHL, ReisingerKS, HogermanDAet al. (1998) Safety and immunogenicity of heptavalent pneumococcal vaccine conjugated to CRM197 in United States infants. Pediatrics101: 604-611.10.1542/peds.101.4.604 PubMed: 95219419521941PilishviliT, LexauC, FarleyMM, HadlerJ, HarrisonLHet al. (2010) Sustained reductions in invasive pneumococcal disease in the era of conjugate vaccine. J Infect Dis201: 32-41.10.1086/648593 PubMed: 1994788119947881MillerE, AndrewsNJ, WaightPA, SlackMP, GeorgeRC (2011) Herd immunity and serotype replacement 4 years after seven-valent pneumococcal conjugate vaccination in England and Wales: an observational cohort study. Lancet Infect Dis11: 760-768.10.1016/S1473-3099(11)70090-1 PubMed: 2162146621621466