ORIGINAL RESEARCH article
Carriage Dynamics of Pneumococcal Serotypes in Naturally Colonized Infants in a Rural African Setting During the First Year of Life
- 1Parasites and Microbes Programme, Wellcome Sanger Institute, Cambridge, United Kingdom
- 2Darwin College, University of Cambridge, Cambridge, United Kingdom
- 3Medical Research Council (MRC) Unit The Gambia at the London School of Hygiene and Tropical Medicine, Fajara, Gambia
- 4Respiratory Diseases Branch, Centers for Disease Control and Prevention, Atlanta, GA, United States
- 5Hubert Department of Global Health, Rollins School of Public Health, Emory University, Atlanta, GA, United States
- 6Emory Global Health Institute, Emory University, Atlanta, GA, United States
- 7Department of Infection, Immunity and Inflammation, University of Leicester, Leicester, United Kingdom
- 8RAMBICON Immunisation & Global Health Consulting, Lekki, Nigeria
- 9Warwick Medical School, University of Warwick, Coventry, United Kingdom
- 10Department of Pathology, University of Cambridge, Cambridge, United Kingdom
- 11NIHR Global Health Research Unit on Mucosal Pathogens, Division of Infection and Immunity, University College London, London, United Kingdom
Streptococcus pneumoniae (the pneumococcus) carriage precedes invasive disease and influences population-wide strain dynamics, but limited data exist on temporal carriage patterns of serotypes due to the prohibitive costs of longitudinal studies. Here, we report carriage prevalence, clearance and acquisition rates of pneumococcal serotypes sampled from newborn infants bi-weekly from weeks 1 to 27, and then bi-monthly from weeks 35 to 52 in the Gambia. We used sweep latex agglutination and whole genome sequencing to serotype the isolates. We show rapid pneumococcal acquisition with nearly 31% of the infants colonized by the end of first week after birth and quickly exceeding 95% after 2 months. Co-colonization with multiple serotypes was consistently observed in over 40% of the infants at each sampling point during the first year of life. Overall, the mean acquisition time and carriage duration regardless of serotype was 38 and 24 days, respectively, but varied considerably between serotypes comparable to observations from other regions. Our data will inform disease prevention and control measures including providing baseline data for parameterising infectious disease mathematical models including those assessing the impact of clinical interventions such as pneumococcal conjugate vaccines.
The pneumococcus continues to kill >320,000 children under 5 years old every year globally despite the use of highly effective serotype-specific pneumococcal conjugate vaccines (PCV) (1). Nasopharyngeal pneumococcal carriage is an essential precursor to invasive pneumococcal disease (2, 3). Furthermore, important biological and ecological events such as evolution and transmission of the pneumococcus occur during carriage (4–6). These fundamental processes influence the population-level characteristics of the strain types and serotypes including the distribution, interactions, and phenotypes such as virulence, vaccine escape, antimicrobial resistance, and immune evasion. Therefore, understanding patterns of pneumococcal carriage remains an essential component of pneumococcal epidemiology for inferences and prediction of the strain dynamics.
Mathematical models are increasingly becoming useful for investigating disease dynamics. This includes understanding transmission patterns (7), antibiotic resistance (5, 8), and impact of infant vaccination programs (9). The accuracy of these models depends on the availability of reliable data to initialize the parameters; therefore, it is crucial to conduct studies to generate such required information. Such data include carriage dynamics specifically time to and rates of acquisition/reacquisition, and carriage duration and clearance rates of different strains. For some pathogens, such as S. pneumoniae, the distribution of strain types or serotypes varies geographically (10). Therefore, studies to describe carriage dynamics should be conducted in different countries to account for the geographical heterogeneity, which may result in accurate inferences from the models.
Assessment of pneumococcal carriage dynamics requires conducting both cross-sectional and longitudinal studies. Cross-sectional carriage surveys are generally easier and cheaper to conduct; therefore, are widely used to provide a snapshot of pneumococcal populations. Conversely, longitudinal studies are usually prohibitively expensive to conduct since they require following up the participants for a period of time. Few high-quality longitudinal studies on pneumococcal carriage have been conducted to date in infants from The Gambia (11, 12), Kenya (13), South Africa (14), Thailand (15), and Papua New Guinea (16). However, due to the geographical differences in serotype distribution between settings, additional studies are required to capture the characteristics of ≥100 known pneumococcal serotypes and to better understand their geographical stability, and heterogeneity between different settings (17, 18). Furthermore, in some settings such as The Gambia, the majority of the pneumococcal carriage studies have focused on describing the carriage dynamics of only a few of the most common serotypes (11, 12).
In this study, we used serotyping data from longitudinally sampled pneumococcal isolates to assess acquisition, reacquisition, carriage duration and clearance of serotypes during the first year of life in newborn infants from a rural setting in the Gambia, West Africa. Previous studies have shown high carriage rates reaching 100% among infants in the Gambia, which makes this setting ideal for this study (19).
Study Design and Participants
The Sibanor Nasopharyngeal Microbiome (SNM) study collected nasopharyngeal swabs (n = 1,595) from 102 newborn infants from 21 villages in the Gambia between September 2008 and April 2009 as previously described (20). The infants were from communities with and without exposure to the seven-valent PCV (PCV7). Parents and guardians of the infants gave informed consent before enrolling them into the study. Nasopharyngeal swabs were then taken from the infants starting from the first week after birth followed by swabs every 2 weeks until 6 months (weeks 1, 3, 5 to 27) and then every two months until their first birthday (weeks 35, 43, and 52). Two-thirds of the infants received PCV7 at 2, 3, and 4 months, and the remaining third received at least 1 dose of PCV during the PCV7 implementation nationwide catch-up campaign. Nasopharyngeal (NPS) specimens were stored in skim milk–tryptone-glucose glycerol medium (STGG) and stored at −80°C within 8 hours of collection.
Sample Processing, Serotyping, and Sequencing
For the isolation of S. pneumoniae, broth enrichment of nasopharyngeal swab samples (NPS) using 5 mL of Todd-Hewitt broth (Oxoid, Basingstoke, UK) containing 5% yeast extract with 1 mL rabbit serum (TCS Biosciences Ltd, Botolph Claydon, UK) was performed as described elsewhere (19). We identified pneumococci by their morphology and optochin sensitivity. Sterile saline suspensions of pneumococcal plate sweeps (on gentamicin blood agar) were then used for serotyping by latex agglutination which can detect multiple serotypes (15). Latex agglutination was performed by capsular and factor-typing sera (Statens Serum Institut, Copenhagen, Denmark) as described previously (21). We also utilized genotypically determined serotype information inferred from whole genome sequencing data for some of the NPS samples as previously described (22) through a collaboration with the Global Pneumococcal Sequencing (GPS) project (www.pneumogen.net) (23).
We defined a colonization episode to be either the first acquisition or re-acquisition of a serotype after clearance; whereby clearance is referred to as the detection of negative cultures for any serotype at two consecutive sampling points as previously described (15). For episodes which commenced at the first week after birth and terminating before week 27, acquisition and clearance of the serotype was assumed to occur at the mid-point between consecutive sampling points. From week 35 to week 52 where samples were collected bi-monthly, acquisition and clearance of serotypes was assumed to occur 1 week before and after detection to prevent overestimation of the carriage duration.
We then fitted a time-to-event parametric model with constant hazard rate (exponential distribution) to the data on time until clearance of serotypes, to estimate the clearance rate as the hazard rate of each serotype. The carriage duration was calculated as the inverse of the hazard rate. Similarly, the acquisition rate and time to acquisition for the first and second carriage episodes with a serotype were estimated as the hazard rate and its inverse after fitting the parametric model to the time until first acquisition data for each serotype. The differences between time-to-event estimates were assessed using log-rank test using non-parametric time-to-event models. We used the survival version 2.42.3 to fit the parametric time-to-event model (24, 25).
We also obtained point estimates for the carriage duration of serotypes in other countries where similar studies were conducted, namely Kenya (13) and South Africa (14). Correlation between the estimates were assessed using the Pearson's correlation coefficient (ρ). The median differences between carriage durations were assessed using Wilcoxon test for paired samples. Graphs were plotted using ggplot2 version 3.3.2 (26) while network visualization was done using Cytoscape version 3.8.1 (27). All the statistical analyses were conducted in R version 3.5.3 (28).
Overall, 1,553 samples from 98 out of 102 newborn infants recruited in the SNM study were analyzed (Figure 1 and Table 1). We identified 80 pneumococcal serotypes naturally colonizing the infants, and pneumococcal-culture positivity rate was 79.3% (1,232/1,553). The mean number of serotypes detected per infant was ≈9 (range: 3–15) while an average of ≈9 carriage episodes (range: 2–15) were detected per infant. By counting each serotype once per episode, serotypes associated with the most carriage episodes were 19A (11.4%), 6A (8.74%), non-typeable strains (NT) (5.71%), 15B/C (4.90%), 19F (3.85%), 23B (4.31%), 34 (3.61%), 21 (3.26%), 35B (3.26%), and 11A (3.03%) (Supplementary Figure 1). All the infants experienced at least 1 carriage episode during the first year of life. Approximately 31% of the infants were colonized within the first week of life but the prevalence increased rapidly to 95% by 2 months and nearly 100% by 12 months (Figure 2A). We then assessed the temporal prevalence of colonization with multiple serotypes during the first year of life. The multiple colonization rate was >40% by 12 months (Figures 2B,C). The proportion of infants carrying multiple serotypes started to plateau approximately 10 weeks after birth mirroring the pattern seen for the overall prevalence of serotypes colonizing the infants. Some serotypes appeared to co-colonize the same infant, for example 6A and 19A, NT and 38, 23B and 9L, 6B and NT (Figure 3).
Figure 1. Flow diagram showing the number of infants and pneumococcal samples included in the analysis. The infants in community 1 were from PCV-unexposed (unvaccinated) villages and received PCV7 after 6 months (control group). The infants in community 2 were also from PCV-unexposed (unvaccinated) and they received PCV7 at 2, 3, and 4 months (direct impact of vaccination). The infants from community 3 came from PCV-exposed (vaccinated) villages and they received PCV7 at 2, 3, and 4 months (direct impact of vaccination and herd immunity).
Figure 2. Temporal prevalence of pneumococcal serotypes during the first year of life. (A) Proportion of infants colonized with any serotype at each sampling point. (B) Proportion of infants co-colonized with multiple serotypes at each sampling point. (C) Proportion of infants colonized by a different number of serotypes at each sampling point. The estimates are shown at bi-weekly sampling points from weeks 1 to 27, and then bi-monthly sampling intervals from weeks 35 to 52.
Figure 3. Network graph showing co-colonization of pneumococcal strains among Gambian infants. The width of the edges connecting a pair of nodes (or serotypes) in the network is proportional to the number of co-occurrences of the pair of serotypes at the same sampling point in the same infant. The shades of color from light green to indigo for each connecting line in the network correspond to the number of detected serotype co-occurrences ranging from 0 to 20.
The time to acquisition of any serotype was 24.1 days [95% confidence interval (CI): 19.7–29.3], and the overall carriage duration was 38.0 days (95% CI: 35.5–40.6) (Figures 4A–C). These implied overall pneumococcal acquisition and clearance rates of 0.042 (95% CI: 0.034–0.051) and 0.026 (95% CI: 0.025–0.028) episodes/day, respectively. We then assessed clearance, acquisition and the reacquisition of individual serotypes (Figures 4A–C, Supplementary Table 1 and Supplementary Figures 2–8). The mean carriage duration of different serotypes ranged from 13.9 to 69.3 days. Serotypes 6B, 4, 13, 23F, and 19F showed longest the carriage durations, and therefore, lowest clearance rates while serotypes with the shortest carriage durations included 1, 5, and 47F (Table 2). The quickest time to first acquisition was associated with serotypes 18, 48, and 47 in contrast delayed acquisition was associated with serotypes 10A, 19C, and 16A.
Figure 4. The Kaplan-Meier survival curves and exponential fit to the longitudinal pneumococcal carriage data. The plots show (A) duration from acquisition to clearance of all serotypes, (B) time until first acquisition of any serotype and (C) time until second acquisition or reacquisition of any serotype. The black curve represents the Kaplan-Meier estimates while the yellow line is the fitted survival curve. The rate parameter in the exponential model is shown at the top of the plots (A–C) represents the mean clearance, acquisition and reacquisition rate, respectively. The inverse of the rate in plots A, B and C equates to mean the carriage duration, time until first acquisition and reacquisition of any serotype, respectively.
Table 2. Carriage duration, initial acquisition and reacquisition times of pneumococcal serotypes in the among Gambian infants.
We then assessed whether time to acquisition and reacquisition of serotypes from birth were similar (Figure 5 and Table 2, Supplementary Table 1 and Supplementary Figure 9). There were no differences in the time to acquisition and reacquisitions of the same serotype from birth were similar except for serotypes 6A, 23F, and NT, which are known pediatric serotypes (29, 30) (Figure 5C and Supplementary Figure 10C). Similarity of the carriage duration after initial colonization was also investigated. No differences were observed for all serotypes except for serotype 11A, which showed longer carriage durations after reacquisition as compared to the initial (P = 0.0026) (Figure 5D and Supplementary Figure 10D). Furthermore, the carriage duration and time to first acquisition of serotypes between PCV-exposed and unexposed communities were also similar (Supplementary Figures 9, 10).
Figure 5. Carriage dynamics of pneumococcal serotypes in Gambian infants. Graph showing (A) carriage duration, (B) time to first acquisition, (C) time to acquisition from birth for the first and second colonization episode with the serotype, (D) carriage duration for the first and second colonization episode with the serotype.
We then assessed whether serotypes acquired earliest had the longest carriage duration (Figures 6A,B). There was no association between carriage duration and time to first acquisition (ρ = −0.127, P = 0.381). In contrast, serotypes which were more prevalent were likely to be carried for longer durations than those carried for short durations (ρ = 0.404, P = 0.003), but there was no association between carriage prevalence and time to initial acquisition of serotypes (ρ = −0.047, P = 0.745) (Figures 6C,D).
Figure 6. Relationship of carriage estimates of pneumococcal serotypes in Gambian infants. (A) Scatter plot showing carriage duration and time to first acquisition. (B) Scatter plot showing carriage duration, second acquisition (reacquisition). (C) Scatter plot showing carriage duration and frequency of serotype/episode, (D) time to first acquisition and frequency of serotypes. (E) Scatter plot showing carriage duration of serotypes in the Gambia and Kenya. (F) Scatter plot showing time to acquisition of serotypes in the Gambia and South Africa.
We compared the points estimates for the carriage duration of serotypes in the present study and similar studies in infants from Kenya (13) and South Africa (14) (Figures 6E,F). The carriage duration estimates from Thailand were available for few serotypes only as such were excluded (15). A weak correlation was observed between carriage duration of serotypes between the Gambia and South Africa (ρ = 0.302, P = 0.118), and the Gambia and Kenya (ρ = 0.357, P = 0.094).
We have described the carriage duration and acquisition patterns of pneumococcal serotypes during the first year of life in the Gambia. We show high and rapid colonization of newborn infants, one third of which carried multiple serotypes. Colonization rates exceeded 20% one week after birth mostly due to serotypes 6A, 34, and NTs. Such early acquisition of certain serotypes may reflect rapid loss of maternally derived immunity (31, 32). Carriage duration of serotypes varied greatly but serotypes carried for longer duration such as 6B were not necessarily acquired earlier in life than serotypes carriage for short durations such as serotype 1. We also found that time from birth to first acquisition and to first reacquisition from birth were similar, which suggests continuous exposure leading to reacquisition of serotypes, possibly among household contacts. The absence of differences for the time to acquisition and carriage duration of serotypes among the infants from vaccinated, vaccine-exposed and unvaccinated communities could be indicative of relatively low herd immunity from PCV7 roll-out at the time of the study. Our findings provide a detailed description of the pneumococcal carriage dynamics in The Gambia, which complements previous pneumococcal carriage studies in the same setting (11, 12).
Since the prevalence of serotypes in carriage varies between countries (10), this may affect how the strains interact and compete in different settings, which may influence carriage dynamics within hosts. By comparing durations between the Gambia (this study), Kenya (13) and South Africa (14), we show that carriage duration is stable geographically for the majority of serotypes although substantial heterogeneity was observed for some serotypes, potentially reflecting variability in local selective pressures. Our study also provides carriage duration estimates for serotypes whose carriage duration were previously unknown in our setting including serotypes 5, 7F, 39, 9A, 40, 28F, and 12F (13–15). Availability of this information will lead to reliable inferences for mathematical models for assessing the population effects of clinical interventions such as PCVs and antibiotics in our geographical setting.
Our study has some limitations. Although our serotyping approach detected colonization with multiple serotypes, it did not quantify their abundance and may have missed less abundant serotypes. To generate a better picture of co-colonization of serotypes, future studies should utilize other high resolution techniques including microarray and deep sequencing of plate sweeps enriched for the pneumococcus (33, 34). Due to the small sample size, our study was not powered to detect potential synergistic or antagonistic interactions between co-carried serotypes on carriage duration of each serotype. In addition, larger sample sizes are needed to obtain robust estimates for rare serotypes. Lastly, the observed carriage duration of serotypes sampled at non-consecutive time points such as serotypes 1, 5, 16A, 41, 19C, 9N, and 47F may have been slightly overestimated. This may be due to the assumption that colonization episodes start and terminate at the midpoint between consecutive sampling periods i.e., acquisition at 1 week before detection and 1 week after last detection. Therefore, caution is required when interpreting the estimated carriage duration for these serotypes.
The implementation of higher-valency PCVs has led to remarkable changes in the pneumococcal populations in Sub Saharan Africa (35–38) and globally (39–42). Our study provides the first comprehensive analysis of temporal carriage dynamics of pneumococcal serotypes in West Africa, which will enhance our understanding of the epidemiology in infants, an age group associated with the greatest risk of pneumococcal diseases and mortality, to inform future prevention and control strategies. Furthermore, our findings will provide baseline data for parameterising mathematical models for infectious diseases to assess and forecast population-level effects of PCVs and antibiotic resistance in the region.
Data Availability Statement
The datasets generated for this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
The SNM study was approved by the Medical Research Council Unit The Gambia and the Gambian Government Joint Ethics Committee (approval number: SCC1108). The patients/participants provided their written informed consent to participate in this study.
BK-A, MA, and RA conducted field and sample collection activities in the Sibanor Nasopharyngeal Microbiome (SNM) study. CC, MS, and BK-A conceived and planned the analysis. The Global Pneumococcal Sequencing (GPS) project was led by KK, RBr, LM, and SB. RG and SL performed genome-based serotyping and quality control of sequenced isolates. P-ET, EF-N, RBa, FC, and CO performed microbiology work. CC carried out the analysis. CC and BK-A drafted the initial version of the manuscript. CC, MS, EB, SL, CE, RG, P-ET, RBa, AW, EF-N, FC, CO, LM, KK, RBr, MB, RA, MA, SB, and BK-A reviewed and contributed to discussions of the manuscript. All authors contributed to the article and approved the submitted version.
The study was funded by the Medical Research Council (MRC) Unit The Gambia at the London School of Hygiene and Tropical Medicine, and the Bill and Melinda Gates Foundation (award no. OPP1034556 to KK, RBr, LM, and SB). CC and SB were funded by the Joint Programme Initiative for Antimicrobial Resistance (JPIAMR).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to thank the study participants and their guardians. We would also like to acknowledge support from the Research Molecular Microbiology Team at the Medical Research Council (MRC) Unit The Gambia at the London School of Hygiene and Tropical Medicine, and the Sequencing and Pathogen Informatics, and Genomics of Pneumonia and Meningitis (and Neonatal Sepsis) teams at the Wellcome Sanger Institute. We would also like to acknowledge the reviewers of our manuscript for their excellent suggestions.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2020.587730/full#supplementary-material
1. Wahl B, O'Brien KL, Greenbaum A, Majumder A, Liu L, Chu Y, et al. Burden of Streptococcus pneumoniae and Haemophilus influenzae type b disease in children in the era of conjugate vaccines: global, regional, and national estimates for 2000-15. Lancet Glob Heal. (2018) 6:e744–57. doi: 10.1016/S2214-109X(18)30247-X
4. Chewapreecha C, Harris SR, Croucher NJ, Turner C, Marttinen P, Cheng L, et al. Dense genomic sampling identifies highways of pneumococcal recombination. Nat Genet. (2014) 46:305–9. doi: 10.1038/ng.2895
5. Lehtinen S, Blanquart F, Croucher NJ, Turner P, Lipsitch M, Fraser C. Evolution of antibiotic resistance is linked to any genetic mechanism affecting bacterial duration of carriage. Proc Natl Acad Sci USA. (2017) 114:1075–80. doi: 10.1073/pnas.1617849114
6. Mosser JF, Grant LR, Millar E V., Weatherholtz RC, Jackson DM, Beall B, et al. Nasopharyngeal carriage and transmission of Streptococcus pneumoniae in American Indian households after a decade of pneumococcal conjugate vaccine use. PLoS ONE. (2014) 9:e79578. doi: 10.1371/journal.pone.0079578
7. Hoti F, Erästö P, Leino T, Auranen K. Outbreaks of Streptococcus pneumoniae carriage in day care cohorts in Finland - implications for elimination of transmission. BMC Infect Dis. (2009) 9:102. doi: 10.1186/1471-2334-9-102
9. Flasche S, Ojal J, Le Polain de Waroux O, Otiende M, O'Brien KL, Kiti M, et al. Assessing the efficiency of catch-up campaigns for the introduction of pneumococcal conjugate vaccine: a modelling study based on data from PCV10 introduction in Kilifi, Kenya. BMC Med. (2017) 15:113. doi: 10.1186/s12916-017-0882-9
10. Johnson HL, Deloria-Knoll M, Levine OS, Stoszek SK, Freimanis Hance L, Reithinger R, et al. Systematic evaluation of serotypes causing invasive pneumococcal disease among children under five: the pneumococcal global serotype project. PLoS Med. (2010) 7:e1000348. doi: 10.1371/journal.pmed.1000348
11. Hill PC, Townend J, Antonio M, Akisanya B, Ebruke C, Lahai G, et al. Transmission of Streptococcus pneumoniae in rural gambian villages: a longitudinal study. Clin Infect Dis. (2010) 50:1468–76. doi: 10.1086/652443
12. Darboe MK, Fulford AJC, Secka O, Prentice AM. The dynamics of nasopharyngeal streptococcus pneumoniae carriage among rural Gambian mother-infant pairs. BMC Infect Dis. (2010) 10:195. doi: 10.1186/1471-2334-10-195
13. Abdullahi O, Karani A, Tigoi CC, Mugo D, Kungu S, Wanjiru E, et al. Rates of acquisition and clearance of pneumococcal serotypes in the nasopharynges of children in Kilifi District, Kenya. J Infect Dis. (2012) 206:1020–9. doi: 10.1093/infdis/jis447
14. Dube FS, Ramjith J, Gardner-Lubbe S, Nduru P, Robberts FJL, Wolter N, et al. Longitudinal characterization of nasopharyngeal colonization with Streptococcus pneumoniae in a South African birth cohort post 13-valent pneumococcal conjugate vaccine implementation. Sci Rep. (2018) 8:12497. doi: 10.1038/s41598-018-30345-5
15. Turner P, Turner C, Jankhot A, Helen N, Lee SJ, Day NP, et al. A longitudinal study of streptococcus pneumoniae carriage in a cohort of infants and their mothers on the Thailand-Myanmar border. PLoS ONE. (2012) 7:e38271. doi: 10.1371/journal.pone.0038271
16. Smith T, Lehmann D, Montgomery J, Gratten M, Riley ID, Alpers MP. Acquisition and invasiveness of different serotypes of Streptococcus pneumoniae in young children. Epidemiol Infect. (1993) 111:27–39. doi: 10.1017/S0950268800056648
17. Bentley SD, Aanensen DM, Mavroidi A, Saunders D, Rabbinowitsch E, Collins M, et al. Genetic analysis of the capsular biosynthetic locus from all 90 pneumococcal serotypes. PLoS Genet. (2006) 2:e31. doi: 10.1371/journal.pgen.0020031
18. Ganaie F, Saad JS, McGee L, van Tonder AJ, Bentley SD, Lo SW, et al. A new pneumococcal capsule type, 10D, is the 100th serotype and has a large cps fragment from an oral Streptococcus. MBio. (2020) 11:e00937–20. doi: 10.1128/mBio.00937-20
19. Hill PC, Akisanya A, Sankareh K, Cheung YB, Saaka M, Lahai G, et al. Nasopharyngeal carriage of Streptococcus pneumoniae in Gambian Villagers. Clin Infect Dis. (2006) 43:673–9. doi: 10.1086/506941
20. Kwambana-Adams B, Hanson B, Worwui A, Agbla S, Foster-Nyarko E, Ceesay F, et al. Rapid replacement by non-vaccine pneumococcal serotypes may mitigate the impact of the pneumococcal conjugate vaccine on nasopharyngeal bacterial ecology. Sci Rep. (2017) 7:8127. doi: 10.1038/s41598-017-08717-0
21. Turner P, Turner C, Jankhot A, Phakaudom K, Nosten F, Goldblatt D. Field evaluation of culture plus latex sweep serotyping for detection of multiple pneumococcal serotype colonisation in infants and young children. PLoS ONE. (2013) 8:e67933. doi: 10.1371/journal.pone.0067933
22. Chaguza C, Senghore M, Bojang E, Gladstone RA, Lo SW, Tientcheu PE, et al. Within-host microevolution of Streptococcus pneumoniae is rapid and adaptive during natural colonisation. Nat Commun. (2020) 11:3442. doi: 10.1038/s41467-020-17327-w
23. Gladstone RA, Lo SW, Lees JA, Croucher NJ, van Tonder AJ, Corander J, et al. International genomic definition of pneumococcal lineages, to contextualise disease, antibiotic resistance and vaccine impact. EBioMedicine. (2019) 43:338–46. doi: 10.1016/j.ebiom.2019.04.021
24. Therneau TM. A Package for Survival Analysis in R. (2020) Available online at: https://cran.r-project.org/package=survival (accessed November 16, 2020).
27. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303
28. R Core Team. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. (2020). Available online at: https://www.R-project.org/
29. Zangwill KM, Vadheim CM, Vannier AM, Hememvay LS, Greenberg DP, Ward JI. Epidemiology of invasive pneumococcal disease in Southern California: implications for the design and conduct of a pneumococcal conjugate vaccine efficacy trial. J Infect Dis. (1996) 174:752–9. doi: 10.1093/infdis/174.4.752
30. Nielsen SV, Henrichsen J. Incidence of invasive pneumococcal disease and distribution of capsular types of pneumococci in Denmark, 1989-94. Epidemiol Infect. (1996) 117:411–6. doi: 10.1017/S0950268800059057
31. Voysey M, Pollard AJ, Sadarangani M, Fanshawe TR. Prevalence and decay of maternal pneumococcal and meningococcal antibodies: a meta-analysis of type-specific decay rates. Vaccine. (2017) 35:5850–7. doi: 10.1016/j.vaccine.2017.09.002
32. Weinberger DM, Dagan R, Givon-Lavi N, Regev-Yochay G, Malley R, Lipsitch M. Epidemiologic evidence for serotype-specific acquired immunity to pneumococcal carriage. J Infect Dis. (2008) 197:1511–8. doi: 10.1086/587941
33. Kandasamy R, Gurung M, Thapa A, Ndimah S, Adhikari N, Murdoch DR, et al. Multi-serotype pneumococcal nasopharyngeal carriage prevalence in vaccine naive nepalese children, assessed using molecular serotyping. PLoS ONE. (2015) 10:e0114286. doi: 10.1371/journal.pone.0114286
34. Knight JR, Dunne EM, Mulholland EK, Saha S, Satzke C, Tothpal A, et al. Determining the serotype composition of mixed samples of pneumococcus using whole genome sequencing. bioRxiv. (2019) 741603. doi: 10.1101/741603
35. Mackenzie GA, Hill PC, Jeffries DJ, Hossain I, Uchendu U, Ameh D, et al. Effect of the introduction of pneumococcal conjugate vaccination on invasive pneumococcal disease in The Gambia: a population-based surveillance study. Lancet Infect Dis. (2016) 16:703–11. doi: 10.1016/S1473-3099(16)00054-2
36. Hammitt LL, Etyang AO, Morpeth SC, Ojal J, Mutuku A, Mturi N, et al. Effect of ten-valent pneumococcal conjugate vaccine on invasive pneumococcal disease and nasopharyngeal carriage in Kenya: a longitudinal surveillance study. Lancet. (2019) 393:2146–54. doi: 10.1016/S0140-6736(18)33005-8
37. Swarthout TD, Fronterre C, Lourenço J, Obolski U, Gori A, Bar-Zeev N, et al. High residual carriage of vaccine-serotype Streptococcus pneumoniae after introduction of pneumococcal conjugate vaccine in Malawi. Nat Commun. (2020) 11:2222. doi: 10.1038/s41467-020-15786-9
38. von Gottberg A. Effectiveness of the 13-valent pneumococcal conjugate vaccine against invasive pneumococcal disease in South African children: a case-control study. Lancet Glob Health. (2017) 5:e359–69. doi: 10.1016/S2214-109X(17)30043-8
39. Dunne EM, Satzke C, Ratu FT, Neal EFG, Boelsen LK, Matanitobua S, et al. Effect of ten-valent pneumococcal conjugate vaccine introduction on pneumococcal carriage in Fiji: results from four annual cross-sectional carriage surveys. Lancet Glob Heal. (2018) 6:E1375–85. doi: 10.1016/S2214-109X(18)30383-8
40. Lo SW, Gladstone RA, van Tonder AJ, Lees JA, du Plessis M, Benisty R, et al. Pneumococcal lineages associated with serotype replacement and antibiotic resistance in childhood invasive pneumococcal disease in the post-PCV13 era: an international whole-genome sequencing study. Lancet Infect Dis. (2019) 19. doi: 10.1016/S1473-3099(19)30297-X
41. Moore MR, Link-Gelles R, Schaffner W, Lynfield R, Lexau C, Bennett NM, et al. Effect of use of 13-valent pneumococcal conjugate vaccine in children on invasive pneumococcal disease in children and adults in the USA: analysis of multisite, population-based surveillance. Lancet Infect Dis. (2015) 15:301–9. doi: 10.1016/S1473-3099(14)71081-3
42. Agudelo CI, Castañeda-Orjuela C, Brandileone MC de C, Echániz-Aviles G, Almeida SCG, Carnalla-Barajas MN, et al. The direct effect of pneumococcal conjugate vaccines on invasive pneumococcal disease in children in the Latin American and Caribbean region (SIREVA 2006-17): a multicentre, retrospective observational study. Lancet Infect Dis. (2020). doi: 10.1016/S1473-3099(20)30489-8. [Epub ahead of print].
Keywords: pneumococcus, serotype, carriage duration, acquisition, Africa, infants
Citation: Chaguza C, Senghore M, Bojang E, Lo SW, Ebruke C, Gladstone RA, Tientcheu P-E, Bancroft RE, Worwui A, Foster-Nyarko E, Ceesay F, Okoi C, McGee L, Klugman KP, Breiman RF, Barer MR, Adegbola RA, Antonio M, Bentley SD and Kwambana-Adams BA (2021) Carriage Dynamics of Pneumococcal Serotypes in Naturally Colonized Infants in a Rural African Setting During the First Year of Life. Front. Pediatr. 8:587730. doi: 10.3389/fped.2020.587730
Received: 07 August 2020; Accepted: 30 November 2020;
Published: 08 January 2021.
Edited by:Robert Cohen, Independent Researcher, France
Reviewed by:Philipp Agyeman, University of Bern, Switzerland
Reinhard Berner, Technische Universität Dresden, Germany
Copyright © 2021 Chaguza, Senghore, Bojang, Lo, Ebruke, Gladstone, Tientcheu, Bancroft, Worwui, Foster-Nyarko, Ceesay, Okoi, McGee, Klugman, Breiman, Barer, Adegbola, Antonio, Bentley and Kwambana-Adams. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors have contributed equally to this work