INTRODUCTION
Streptococcus pneumoniae (pneumococcus) is an important human pathogen that may cause a wide range of infections including severe invasive disease such as bacteremia and meningitis. Introduction of the seven-valent pneumococcal conjugate vaccine (PCV7) in the USA in 2000 and, later, in other countries has strongly reduced the incidence of PCV7 vaccine serotype invasive pneumococcal disease (IPD) there [Reference Pilishvili1–Reference van der Linden3]. At the same time, a growing rate of non-PCV7 vaccine serotypes, primarily serotype 19A, has been documented among circulating pneumococci in the post-PCV7 period [Reference Richter4, Reference Hackel5]. In the USA, the incidence of serotype 19A-related IPD has increased from 0·7 to 2·5–2·6 cases per 100 000 population between 1999 and 2008 [Reference Beall6]. According to data from Germany, the proportion of serotype 19A isolates among IPD pneumococci has increased up to 15% in 2010–2011 from a basal range of 1·7–4·2% observed in the period between 1997 and 2006 [Reference van der Linden7]. A similar trend was observed in other European countries after implementation of universal PCV7 vaccination [Reference Aguiar8–Reference Del Grosso10]. These increases were attributed to the PCV-mediated serotype replacement phenomenon, i.e. the proliferation of pre-existing minor non-PCV7 serotype pneumococcal populations under vaccine pressure [Reference Weinberger, Malley and Lipsitch11]. However, in some territories, for instance, in Israel and South Korea, a rise of serotype 19A was observed in the pre-PCV7 period [Reference Choi12–Reference Hsieh14] suggesting that, in addition to serotype, other bacterial factors might be relevant for the expansion.
The proliferation of serotype 19A pneumococci displays clonality. Pneumococcal clones, or clonal complexes (CCs), are readily distinguishable by multilocus sequence typing (MLST) [Reference Enright and Spratt15, 16]. They demonstrate a variable distribution over time and territory as well as diversity in terms of virulence, antimicrobial resistance, and overall fitness [Reference Henriques-Normark17]. Reportedly, only a few clones have contributed to the expansion of serotype 19A including CC199, CC230, CC320, CC695, and CC994 [Reference Beall6–Reference Aguiar8, Reference Del Grosso10, Reference Moore18–Reference Metcalf23]. Most of these clones have reduced antimicrobial susceptibility revealing a multidrug resistant (MDR) phenotype with high level of resistance to penicillin, which is related to alterations of penicillin-binding proteins (PBPs) [Reference Hakenbeck24]. The clone of sequence type (ST) 320, the major, globally emerged serotype 19A clone, and CC199 (ST416) have been shown to carry pili that may provide an additional selective advantage [Reference Del Grosso10, Reference Metcalf23, Reference Zähner, Gudlavalleti and Stephens25].
Data on pneumococcal molecular epidemiology in Russia are limited. A recent report from our group documented the presence of serotype 19A CC320 in Russia before implementation of PCV vaccination in the National immunization program in 2014 [Reference Mayanskiy26]. In the present study, we retrospectively examined a set of non-invasive serotype 19A pneumococci collected in Russia during 2002–2013. We describe clonal diversity changing over this time period in conjunction with antimicrobial susceptibility, pilus prevalence, and PBP sequence alterations.
METHODS
Isolate collection
This retrospective study included serotype 19A pneumococcal isolates obtained from children from different parts of Russia during 2002–2013. One part of this collection was selected from a set of non-invasive isolates recovered by the Institute of Antimicrobial Chemotherapy, Smolensk, Russia, in 2002–2008 [Reference Stratchounski27]. Another part of the collection comprised of serotype 19A pneumococci discovered during a multicenter study conducted in Moscow in 2010–2013 at the Scientific Centre of Children's Health [Reference Mayanskiy28]. The isolates were obtained from respiratory tract specimens, including nasopharyngeal swabs, middle ear fluid, sputum, and bronchoalveolar lavage, i.e. were non-invasive. No invasive serotype 19A isolates were recovered in these studies.
Serotyping and antimicrobial susceptibility testing
Laboratory examination of the isolates was performed at the Scientific Centre of Children's Health. All isolates underwent repeated serotyping by Quellung reaction using Staten Serum Institut antisera (Copenhagen, Denmark). Antimicrobial susceptibility for penicillin (PEN), amoxicillin (AMX), ceftriaxone (CTX), and erythromycin (ERY) was examined by the E-test (bioMerieux, France); the remaining antimicrobials including clindamycin (CLI), chloramphenicol (CHL), trimethoprim/sulfamethoxazole (SXT), and tetracycline (TET) were tested by the disk diffusion method (Bio-Rad, Hercules, CA). Minimum inhibitory concentration (MIC) category and zone diameter breakpoint interpretations were based on updated EUCAST-2013 standards. Pneumococcal isolates with PEN MICs >0·06 mg/l and >2 mg/l were regarded as PEN-nonsusceptible and PEN-resistant, respectively. A PCR was used to detect the erm(B) and mef determinants in ERY-resistant pneumococci, as described previously [Reference Mayanskiy28]. An MDR phenotype was defined as nonsusceptibility to three or more classes of antibiotics.
MLST procedure
MLST analysis was performed according to the S. pneumoniae MLST protocol [Reference Enright and Spratt15, 16]. Alleles and STs were assigned using the software at the pneumococcal web page (http://pubmlst.org/spneumoniae). Analysis of the STs and assignment to CC was performed against all STs found in the online database using the eBURST program. The STs that shared at least five of seven allelic variants composed a CC. Trace files for six new MLST alleles (recP, n = 2; dll, n = 4) and eight new MLST profiles were submitted to the MLST database for designations.
Detection of the pilus islets 1 and 2 (PI1 and PI2)
The presence of PI1 was detected by PCR amplification of an internal fragment of the rlrA gene using primers RLRA-F (5′-TCTGATAGATGAGACGCTGTTG-3′) and RLRA-R (5′-CTCCGCTTCTTTCTACTACAAG-3′), as described elsewhere [Reference Aguiar29]. For detection of PI2, primers for PCR amplification of conserved regions within the PI2 were modified from Bagnoli et al. [Reference Bagnoli30] as follows: PI2-F (5′-CGTGGGTATCAGGTGTCCTATG-3′) and PI2-R (5′-TGCAGTGAATAGCTTTTTAAAGAA-3′).
Amplification, sequencing, and genetic analyses of pbp1a, pbp2b, and pbp2x genes
Fragments of pbp1a, pbp2b, and pbp2x including the coding regions of the penicillin-binding domains and adjacent sequences were amplified from genomic DNA and sequenced using the Sanger method, as described elsewhere [Reference Granger31]. Amplified fragments corresponded to the following amino acid positions: PBP1a, positions 341–582; PBP2b, positions 375–652; PBP2x, positions 337–661. Deduced amino acid sequences were compared with the corresponding sequences of the reference strain S. pneumoniae R6; pneumococcal clones of Taiwan19F-14 (ST236), Spain9V-3 (ST156), and MDR19A (ST320). We used the same GenBank accession numbers, as described earlier [Reference Ip32]. The relatedness of discovered amino acid sequences was depicted using neighbor-joining trees constructed with Molecular Evolutionary Genetics Analysis (MEGA) software [Reference Ip32].
Statistical analysis
Statistical analysis was performed using IBM SPSS Statistics for Windows, version 20·0 (IBM Corp, Armonk, NY). Contingency table analysis for comparing proportions was done by the χ 2 test, or by means of the z-test, where appropriate. The tests were considered statistically significant at P < 0·05.
RESULTS
Description of the serotype 19A pneumococcal collection
This study included 49 serotype 19A pneumococcal non-invasive isolates obtained from children (median age, 2·76 years; IQR (interquartile range) 2·59 years). Eighteen isolates were collected in 2003 in European Russia (cities of Moscow, Saint Petersburg, and Smolensk); nine isolates were obtained in Asian Russia in 2002–2008 (cities of Anadyr, Irkutsk, Novokuznetsk, and Yekaterinburg); 22 isolates were discovered in Moscow in 2010–2013 (Supplementary Table S1). The overall rate of 19A pneumococci was 6·4% (30 of 468; 27 isolates were available for further analyses) and 2·2% (22 of 1002) for the 2002–2008 and 2010–2013 collections, respectively (Z = 3·7, P < 0·001). The isolates were recovered from the nasopharynx (67%), middle ear fluid (18%), and lower respiratory tract specimens (sputum and bronchoalveolar lavage fluid; 15%).
Clonal evolution and antimicrobial resistance
The clonal lineage and ST diversity of serotype 19A pneumococci from time and geographic perspective is presented in Table 1 and Fig. 1. Overall, 25 STs were discovered, including 14 newly described STs, that belonged to 6 major clonal lineages; three STs were singletons. The distribution of clonal lineages over the study period was significantly different (χ 2 = 29·60, P < 0·001). The European-2003 isolates were represented by only two clonal lineages, CC663 (n = 12) and CC156 (n = 6). In contrast, among isolates recovered in 2010–2013, these clones had a minor prevalence being largely replaced by CC230 and CC320, which collectively comprised 73% (16 of 22) of the 2010–2013 isolates. Five of nine strains from Asian Russia were singletons; three isolates belonged to CC230.
* PMEN clones and their identifier STs (in brackets) related to the corresponding CC are indicated in bold typeface. For newly described STs, the relatedness to the corresponding clone identifier ST is shown.
† Significant differences between the proportions of 2003 and 2010–2013 (P < 0·05).
‡ Significant differences between the proportions of 2002–2008 and 2003, 2010–2013 (P < 0·05).
The most abundant clone, the clone of ST663 (n = 14, 28% of the entire collection), was represented by four related STs (Table 1). CC663 predominated in European Russia in 2003 (67%); in 2010–2013, only two CC663 isolates were recovered. The majority of CC663 representatives (79%) had an MDR phenotype with high β-lactam MICs (Tables 2 and 3). All CTX-resistant isolates (n = 4; MIC > 2 mg/l) belonged to this CC (Table 3). Moreover, the isolates of ST10434 (n = 5) had an extremely drug-resistant phenotype being nonsusceptible to five out of six antibiotic groups tested, including PEN, ERY, CLI, CHL, and TET (Supplementary Table S1). Among the examined antibiotics, the ST10434 isolates were only susceptible to SXT.
AMX, amoxicillin; CC, clonal complex; CHL, chloramphenicol; CLI, clindamycin; CTX, ceftriaxone; ERY, erythromycin; MDR, multidrug resistant; PEN, penicillin; SXT, trimethoprim/sulfamethoxazole; TET, tetracycline.
* The percentage indicates proportions of isolates carrying a macrolide resistance determinant among ERY-resistant isolates. R, resistant.
† PI, pilus islet.
‡ PEN-S, a group of six PEN-susceptible isolates.
PBP, penicillin-binding protein; CC, clonal complex; MIC, minimum inhibitory concentration; PEN, penicillin; AMX, amoxicillin; CTX, ceftriaxone; S, susceptible; I, intermediate; R, resistant.
* The proportion of altered amino acids was calculated by comparing amino acid sequences of the conserved motifs of PBP1a (positions 341–582), PBP2b (positions 375–652), and PBP2x (positions 337–661) of the study isolates with those of the R6 pneumococcal reference strain. PEN-S isolates were excluded from calculations of the overall mean.
† In six out of seven isolates, an MIC for AMX was higher than an MIC for PEN.
‡ PEN-S group, a group of six PEN-susceptible isolates; described in Table 2.
The second most prevalent clonal lineage was CC230 (n = 12, 25% of the entire collection) that included eight STs (Table 1). The majority of CC230 isolates were recovered in Moscow during 2010–2013 contributing to 41% for this part of the collection; three isolates was found in Asian Russia. CC230 was characterized by moderate-to-low β-lactam MICs being largely susceptible to CLI, CHL, and TET, but resistant to SXT (Tables 2 and 3). An MDR phenotype was evident in five (42%) of the CC230 isolates that were ERY-resistant (Supplementary Table S1).
CC156 was represented by seven serotype 19A pneumococci of four STs accounting for 14% of the entire collection. Six out of seven isolates were collected in European Russia in 2003 composing a subclone of three closely related new STs including ST10437, ST10438, and ST10514 (Table 1). These 19A pneumococci possessed a remarkable phenotype being fully resistant to AMX with MICs 4–8 mg/l that exceeded PEN MICs (Table 3, Supplementary Table S1). In addition, these isolates were resistant to SXT, but susceptible to ERY, CLI, CHL, and TET.
CC320 (seven isolates, 14% of the entire collection) was evident exclusively in 2010–2013 having a 32% prevalence (Table 1). It was represented by two STs: ST320 (n = 3) and ST9656 (n = 4; a new SLV of ST320 at the recP locus). All CC320 isolates displayed an MDR phenotype having high PEN and AMX MICs (range 1·0–4·0 mg/l), resistance to ERY and SXT (Tables 2 and 3). All but one isolate had higher MICs for AMX than for PEN; all isolates were susceptible to CHL (Supplementary Table S1).
Nine isolates represented various clonal lineages or unrelated STs. All six PEN-susceptible serotype 19A pneumococci belonged to this group of isolates (Tables 2 and 3). None of the PEN-susceptible isolates had an MDR phenotype. One singleton (ST10436) had high PEN and AMX MICs (2 mg/l) (Table 3).
In the European isolates, the overall nonsusceptibility rate to tested β-lactams was similar in 2003 and 2010–2013; however, the proportion of fully PEN-resistant isolates (MIC > 2 mg/l) declined from 28% (5 of 18) to 5% (1 of 22) during this period (P = 0·041). Although the prevalence of ERY-resistant serotype 19A pneumococci remained unchanged (67–68%), the proportion of ERY-resistant isolates carrying both erm(B) and mef increased from 17% (2 of 12) in 2003 to 60% (9 of 15) in 2010–2013 (P = 0·023). The resistance rate to CHL significantly declined over time from 33% (6 of 18) to 5% (1 of 22) (P = 0·017). In contrast, the prevalence of SXT-resistant isolates increased from 44% (8 of 18) to 82% (18 of 22) (P = 0·014) during the study period.
The prevalence of PI carriage in relation to clonal distribution and antimicrobial resistance
In total, 29 isolates (59%) showed the presence of PI1, of which seven isolates carried PI2 in addition; 20 isolates (41%) were PI-negative (Table 2, Supplementary Table S1). The prevalence of PI1 and PI2 had a clear clonal pattern. CC663 and CC156 pneumococci showed the presence of PI1 but not PI2, whereas the CC320 isolates harbored both PIs (Table 2). Among the remaining pneumococci, no isolates carried a PI except for the singleton of ST10436 which was PI1-positive. The rate of MDR phenotypes was significantly higher among PI-positive isolates (Table 4). The vast majority of piliated isolates were nonsusceptible to β-lactam antibiotics, whereas all PI-negative isolates were susceptible to AMX and CTX and had low-to-moderate PEN MICs (Table 4, Supplementary Table S1). In addition, the prevalence of ERY- and CLI-resistant pneumococci was higher among piliated isolates (Table 4).
AMX, amoxicillin; CHL, chloramphenicol; CLI, clindamycin; CTX, ceftriaxone; ERY, erythromycin; MDR, multidrug resistant; PEN, penicillin; SXT, trimethoprim/sulfamethoxazole; TET, tetracycline.
* PI, pilus islet. PI-negative, isolates carrying none of the PI; PI-positive, isolates carrying PI1 or both PI1 and PI2.
† Statistically significant differences in proportions of nonsusceptible isolates among PI-negative and PI-positive pneumococci are in bold typeface.
Sequence analyses of PBPs
Amino acid sequences of the examined PBP fragments were compared with the corresponding sequences from PEN-susceptible S. pneumoniae R6 strain. The number of amino acid substitutions differed between PEN-susceptible and PEN-nonsusceptible isolates, as well as between clonal groups (Table 3, Supplementary Table S2). PEN-susceptible isolates (n = 6) had no alterations in PBP2b and PBP2x; only one amino acid substitution was observed in PBP1a. In PEN-nonsusceptible pneumococci, the average rate of altered amino acids was 13·9%, 7·2%, and 8·7% for PBP1a, PBP2b, and PBP2x, respectively.
The rate of amino acid alteration in PBP1a was similar among different clonal groups of PEN-nonsusceptible serotype 19A isolates (Table 3, Supplementary Table S2). In 33 PEN-nonsusceptible isolates (77%), sequence alterations were found in the conserved motif 370STMK373 (T371A/S) and in flanking amino acids of the 428SRN430 motif (P432T) of PBP1a. In all but one PEN-nonsusceptible isolates, a four-amino acid substitution was present (574TSQF577 replaced by NTGY) (Supplementary Table S2).The amino acid sequence of the conserved motif 557KTG559 of PBP1a was intact.
In PBP2b, sequence alterations were more prevalent among CC156 and CC320 isolates (Table 3). All PEN-nonsusceptible isolates had the substitution of T445A in the conservative motif 442SSNT445 and the mutation T488A/S. All 14 AMX-resistant isolates contained a unique set of 10 amino acid substitutions (A591S, G596P, N605D, L608T, A618G, D624G, Q627E, T629N, S639T, and D640E) in the 591–640 amino acid region of PBP2b that have been associated with reduced susceptibility to AMX [Reference du Plessis, Bingen and Klugman33, Reference Kosowska34]. No mutation was found in, or close to, the conserved motif 385SVVK388 of PBP2b in any study isolate.
Similarly to PBP2b, amino acid alterations in PBP2x were more abundant in CC156 and CC320 (Table 3). The conservative motif 337STMK340 was altered in 31 PEN-nonsusceptible isolate (72%) by the substitution T338A/P. Other prevalent mutations included I371T, R384G/S, Q552E, and N605T (Supplementary Table S2).
Next, we assessed the relatedness of amino acid sequences of PBPs of our study isolates to a number of reference sequences, constructing a neighbor-joining tree (Fig. 2). All PEN-susceptible isolates clustered together with the R6 strain (Fig. 2a–c). The amino acid sequences of PBP1a from the study isolates of CC320 were 100% identical to the reference sequences of MDR19A (ST320) and its precursor, Taiwan19F-14 (ST236) (Fig. 2a). One highly resistant pneumococcus of ST10434, CC663 (Supplementary Table S1, #242; MICs: PEN and AMX, 3 mg/l; CTX, 8 mg/l) clustered to the same isolate branch. CC230 representatives were localized near this group of isolates, too. Amino acid sequences from isolates of CC156 and CC663 had an intermediate position between MDR19A (ST320) and Spain9V-3 (ST156) (Fig. 2a).
For PBP2b, the amino acid sequences of CC320 and CC156 (with the exception of ST143, which was identical to Spain9V-3, ST156) were equal or very similar to the MDR19A (ST320) sequences (Fig. 2b). CC230 and CC663, except for one extremely resistant isolate of ST663 (see below), comprised a big cluster of similar sequences in proximity to R6 and the PEN-susceptible group. Of note, the group of highly resistant serotype 19A pneumococci of CC156 and CC320 possessing the A681G mutation in PBP2b clustered together. One of these isolates, an isolate of ST663 with an extremely high β-lactam resistance (isolate #58, MICs: PEN, 16 mg/ml; AMX, 12 mg/ml; CTX, 32 mg/ml; Supplementary Tables S1 and S2), clustered separately from the remaining CC663 together with MDR19A (ST320) (Fig. 2b).
PBP2x sequences of CC156 were similar to those of Taiwan19F-14 (ST236) and Spain9V-3 (ST156) (Fig. 2c). Study isolates of CC320 and MDR19A (ST320) had fully equivalent PBP2x sequences. This group of isolates clustered separately from CC663, which was more similar to R6.
DISCUSSION
In the present study, we analyzed a set of serotype 19A pneumococci collected over 2002–2013 in several geographic locations in Russia. Four clonal groups including CC156, CC230, CC320, and CC663 composed 80% of the entire collection. These lineages related to major international clones that have different distributions over the world. CC156 has a minor representation among serotype 19A pneumococci, and this has not changed significantly under PCV pressure [Reference Beall6, Reference van der Linden7, Reference Dagan13, Reference Moore18, Reference Setchanova21]. The same was true for the PMEN clone Colombia23F-26, to which our ST663 clonal group was attributed by the eBURST analysis [Reference Beall6, Reference Moore18]. Actually, ST663 has rarely been reported in the literature, although it was associated with high β-lactam resistance and severe IPD [Reference Shouval35]. ST663 pneumococci had noticeable pre-PCV7 prevalence (9·6%) in Bulgaria [Reference Metcalf23]. In the present study, CC663 dominated in 2003 in European Russia. Interestingly, the earliest serotype 19A isolate of ST663 reported to the pubmlst.org/spneumoniae database, was discovered in Moscow in 1987 [16]. Thus, it is possible that ST663 and related genotypes were locally disseminated in the Eastern Europe in the 1990s and early 2000s. Of note, within the study period, we did not recover any isolate of CC199 that had high prevalence in the pre-PCV7 period in the USA and Europe [Reference van der Linden7, Reference Moore18].
In contrast to the above discussed clones, CC230 and CC320 have significantly expanded in the post-PCV7 period. While the proliferation of CC230 was reported only in a few countries, e.g. Germany and Portugal [Reference van der Linden7, Reference Aguiar8], the CC320 lineage is currently among the most prevalent serotype 19A pneumococci in many territories including the USA, Canada, Italy, and Spain [Reference Beall6, Reference Moore18, Reference Gene20, Reference Golden22]. Our results demonstrated that, at present, in European Russia, CC230 and CC320 represent more than 70% of circulating 19A pneumococci; these clones have largely replaced CC156 and CC663. Pneumococci of the ST320 clonal group have been recovered in Russia starting from 2010 [Reference Mayanskiy26]. We found only two CC320-related STs, the identifier ST320 and the new ST9656 that is emerging in Russia. Since our collection has been assembled before the introduction of PCV in the National immunization program in 2014, the genetic structure of serotype 19A pneumococci was evolving in a PCV-independent manner, as it has been reported from some other regions [Reference Choi12–Reference Hsieh14].
The expansion of only a small number of serotype 19A clones implies their competitive advantage at the genotype level. Besides vaccine use, clonal success of CC320 and other serotype 19A clones is attributed to antibiotic selective pressure [Reference Henriques-Normark17]. Indeed, in the present study we observed that all CC320 isolates had an MDR phenotype with high β-lactam MICs, which is in agreement with previous studies [Reference Beall6, Reference Dagan13, Reference Golden22]. However, in European Russia, CC320 has replaced another highly resistant clone, CC663. This suggested that additional bacterial factors may be important for conferring competitive advantages. For instance, serotype 19A pneumococci of ST320 are capable of more effective colonization than their ancestral ST236 strain, at least in a murine model [Reference Hsieh14]. The increased colonization capacity of CC320 may be associated with the presence of PI1 and PI2 [Reference Del Grosso10, Reference Metcalf23, Reference Zähner, Gudlavalleti and Stephens25, Reference Sjöström36], which have been shown to mediate adhesion of pneumococci to eukaryotic cells [Reference Bagnoli30]. In the present study, CC320 was the only clonal lineage carrying both PIs, whereas CC663 isolates were PI1 positive only. Reportedly, the presence of pili is a clonal property [Reference Zähner, Gudlavalleti and Stephens25, Reference Aguiar29], and majority of antimicrobial-nonsusceptible clones does possess these structures [Reference Hjálmarsdóttir37]. This may suggest that piliation and antimicrobial resistance operate in concert conferring competitive benefits.
All isolates of CC156, CC320, and CC663 had PEN MICs ⩾1 mg/l, whereas CC230 representatives had moderate-to-low PEN MICs and were fully susceptible to AMX and CTX. This indicates that β-lactam susceptibility among examined serotype 19A pneumococci has a clonal pattern. In pneumococci, the principal mechanism of β-lactam resistance is mediated through modification of PBPs [Reference Hakenbeck24]. Sequence analyses of penicillin-binding domains of these proteins performed in the present study demonstrated that isolates of related STs and clonal lineages had highly similar amino acid sequences explaining similar β-lactam susceptibility profiles. Clonal lineages had a variable degree of sequence alterations but most of PBP changes found in our collection were typical of those previously described. For instance, the A618G substitution close to the KTG motif of PBP2b has been previously associated with elevated AMX MICs [Reference du Plessis, Bingen and Klugman33, Reference Kosowska34]. In the present study, the A618G mutation conferred an AMX MIC > PEN MIC phenotype in the majority of isolates that carried it. In addition, this mutation was found in two isolates with extremely high β-lactam MICs.
Starting from 2010, the extended-valent PCV13 that covers serotype 19A has replaced PCV7. PCV13 has demonstrated a high effectiveness in reducing the incidence of serotype 19A-related IPD [Reference Moore38, Reference van der Linden39]. Being effective in the repulsion of vaccine serotypes from the circulation, vaccination with PCVs is having much broader impact on the microbial community than has been anticipated. In a recent paper by Galanis et al. [Reference Galanis40], the authors demonstrated that PCV vaccination promoted expansion of new or minor pneumococcal clones and increased clonal diversity not only in vaccinees, but also in non-vaccinated individuals. These findings suggest a possibility of the spread of some serotypes in a PCV-naïve population. Most likely, the changing epidemiology of pneumococcal serotypes, including serotype 19A that we have described in the present study, reflects pressures of several factors such as vaccination and antibiotic use, as well as naturally occurring fluctuations.
SUPPLEMENTARY MATERIAL
The supplementary material for this article can be found at https://doi.org/10.1017/S0950268817000541
Acknowledgements
The authors thank A. Pakhomov for his help in performing MLST.
DECLARATION OF INTEREST
None.