Next Article in Journal
Escherichia coli Heat-Labile Enterotoxin B Subunit Combined with Ginsenoside Rg1 as an Intranasal Adjuvant Triggers Type I Interferon Signaling Pathway and Enhances Adaptive Immune Responses to an Inactivated PRRSV Vaccine in ICR Mice
Next Article in Special Issue
Etiology of Clinical Community-Acquired Pneumonia in Swedish Children Aged 1–59 Months with High Pneumococcal Vaccine Coverage—The TREND Study
Previous Article in Journal
No Such Thing as a Free-Rider? Understanding Drivers of Childhood and Adult Vaccination through a Multicountry Discrete Choice Experiment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pneumococcal Competition Modulates Antibiotic Resistance in the Pre-Vaccination Era: A Modelling Study

1
Department of Zoology, University of Oxford, Oxford OX1 3SZ, UK
2
School of Public Health, Faculty of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
3
Porter School of the Environment and Earth Sciences, Faculty of Exact Sciences, Tel Aviv University, Tel Aviv 69978, Israel
4
NIHR Global Health Research Unit on Mucosal Pathogens, Division of Infection and Immunity, University College London, London WC1E 6BT, UK
*
Authors to whom correspondence should be addressed.
Vaccines 2021, 9(3), 265; https://doi.org/10.3390/vaccines9030265
Submission received: 23 February 2021 / Revised: 6 March 2021 / Accepted: 9 March 2021 / Published: 16 March 2021

Abstract

:
The ongoing emergence of antibiotic resistant strains and high frequencies of antibiotic resistance of Streptococcus pneumoniae poses a major public health challenge. How and which ecological and evolutionary mechanisms maintain the coexistence of antibiotic resistant and susceptible strains remains largely an open question. We developed an individual-based, stochastic model expanding on a previous pneumococci modelling framework. We explore how between- and within-host mechanisms of competition can sustain observed levels of resistance to antibiotics in the pre-vaccination era. Our framework considers that within-host competition for co-colonization between resistant and susceptible strains can arise via pre-existing immunity (immunological competition) or intrinsic fitness differences due to resistance costs (ecological competition). We find that beyond stochasticity, population structure or movement, competition at the within-host level can explain observed resistance frequencies. We compare our simulation results to pneumococcal antibiotic resistance data in the European region using approximate Bayesian computation. Our results demonstrate that ecological competition for co-colonization can explain the variation in co-existence of resistant and susceptible pneumococci observed in the pre-vaccination era. Furthermore, we show that within-host pneumococcal competition can facilitate the maintenance of resistance in the pre-vaccination era. Accounting for these competition-related components of pneumococcal dynamics can improve our understanding of drivers for the emergence and maintenance of antibiotic resistance in pneumococci.

1. Introduction

The Gram-positive bacteria Streptococcus pneumoniae (or pneumococci) remain the most common cause of pneumonia in children and the elderly. Pneumococcal invasive infections are a major cause of death in children, and are likely responsible for over 500,000 annual deaths [1]. Most carriers of pneumococci are asymptomatic, and only a small proportion develops an invasive disease such as meningitis or bacteraemia [2]. Vaccination remains the most successful public health strategy to date of reducing pneumococcal carriage and invasive disease [3,4].
Pneumococci are characterized by a polysaccharide capsule comprising over 100 distinct antigenic types, termed serotypes [5,6]. Pneumococcal conjugate vaccines target seven (PCV7), ten (PCV10) or thirteen (PCV13) serotypes. Since vaccine introduction, a rapid decrease in pneumococcal carriage and disease of the targeted serotypes has been observed [7,8,9,10,11]. Concomitantly, increases in the frequencies of non-targeted serotypes have been recorded both in carriage and disease [12,13,14,15,16,17,18,19]. This phenomenon is known as vaccine-induced serotype replacement (VISR), whereby the removal of targeted serotypes allows non-targeted serotypes to fill their immunological and ecological niches due to the removal of cross-immunity or direct resource competition [20,21,22,23,24]. Changes in the genetic composition of pneumococci have also been observed in the post-vaccination era [25,26]. These genetic changes can reflect perturbations in competition for resources among pneumococci sharing metabolic niches, also termed vaccine-induced metabolic shifts (VIMS) [27,28,29].
A second major route of combating pneumococcal carriage indirectly is by antibiotic treatment of invasive disease [30]. Similarly to vaccination, the effectiveness of antibiotics is often impaired by emergent antibiotic-resistant strains. Generally, vaccination has had a beneficial impact on the frequencies of antibiotic resistance among pneumococci, since vaccine types (serotypes targeted by vaccination, VT) often have increased antibiotic resistance [31,32,33,34]. However, post-vaccination increases in resistance among non-VTs (NVT) have been reported [19,35,36].
Previously, we have put forward hypotheses regarding the possible drivers of post-vaccination increases in resistance among NVTs. VIMS, for example, posits that the genetic profiles associated with resistance in VTs can shift to NVTs enjoying relaxed metabolic competition pressures in the post-vaccination era [28]. Furthermore, in a study by Obolski et al. [37], we have shown that vaccination can reduce within-host competition between co-colonizing serotypes. This reduction in competition may benefit antibiotic resistant strains disproportionately by relaxing the fitness cost of resistance. Vaccination can thus ameliorate the antibiotic resistance cost for NVTs and facilitate the increase in frequency of antibiotic-resistant NVTs, as has been observed [19,35,36].
In this study, we study the co-existence of pneumococcal susceptible and resistant strains in the pre-vaccination era. We do so by generalizing the original framework into an individual-based model. Such a generalization allows us to explore the theoretical effects of stochasticity, host population structure and mobility. We perform a series of sensitivity analyses based on the equilibrium states of stochastic simulations of the pre-vaccination era dynamics. We explore how two within-host competition mechanisms determine the levels of susceptible and resistant strains. Our model reproduces observed variation in the frequencies of resistance to penicillins and macrolides in Europe. In this region, ecological competition between susceptible and resistant strains explains more of the observed differences between countries than does immunological competition. Our results contribute to a better understanding of the processes driving the maintenance of antibiotic resistance in the absence of vaccination.

2. Materials and Methods

2.1. Epidemiological Framework

We model the transmission dynamics of pneumococcal strains according to the epidemiological framework by Obolski et al. [37]. We define pneumococcal strains by the combination of their serotype and antibiotic resistance profile to two distinct antibiotics. Similarly to the nomenclature used in previous studies, each strain’s genotype is conceptualised by a tuple { i , j } , where i determines the serotype and j the antibiotic resistance profile. We model two serotypes only, such that i { a , b } . For the entirety of our study, we define serotype a as a vaccine type (VT) and serotype b as a non-vaccine type (NVT). We also model two antibiotics, x and y , and consider only resistance to each independently, with j { 00 , 01 , 10 } , where j   =   00 refers to sensitive strains, j   =   01 to strains resistant only to antibiotic y , and j   =   10 to strains resistant only to antibiotic x . For instance, a strain defined by i   =   a ,   j   =   00 is of VT serotype a and is sensitive to both antibiotics; while a strain i   =   b ,   j   =   01 is of NVT serotype b , sensitive to antibiotic x and resistant to antibiotic y . For simplicity, we do not consider double resistance, i.e., j   =   11 in our model. This is consistent with relatively low estimates of double resistance to penicillin and macrolides identified in Europe in pneumococci before the pneumococcal vaccine introduction [38,39].
We keep the terminology of the original model formulation regarding the epidemiological parameters: γ is serotype-specific immunity; ψ is the probability that a carried susceptible strain ( j   =   00 ) will suppress host co-colonization by a resistant strain ( j   =   01 or j   =   10 ) due to the fitness cost of antibiotic resistance (details below); 1 / μ is the host life-span; 1 / σ is the host carriage duration; R0 is the basic reproduction number and β is the transmission rate.
Individuals are born naive to all strains, with the size of the host population kept constant (deaths being replaced by newborns). Colonization results in carriage for an average duration of 1 / σ days, from which recovery (clearance) may lead to complete serotype immunity if γ   =   1 or partial life-long immunity if γ < 1 . Thus, if γ < 1 , hosts can be re-colonized by the same serotype throughout their lifetime, with probability 1 γ . Co-colonization of up to two strains is allowed, unless γ   =   1 and strains belong to the same serotype. Within-host strain competition interferes with co-colonization if ψ > 0 , which captures the degree to which intrinsic fitness differences (e.g., growth rates) between resistant and sensitive strains may allow a currently carried sensitive strain ( j   =   00 ) to suppress co-colonization by a resistant strain ( j   =   01 , j   =   10 ). We emphasize that ψ models a form of competition between bacterial strains that is not mediated by immunity—from now onwards referred to ecological competition. The modelled processes related to an individual’s colonization and epidemiological states are summarised in Figure 1A. The full model description is provided in the Supplementary Information.
In this study we use a stochastic, individual-based model (IBM) implementation [40] of the original epidemiological framework [37]. We implement model events by sampling from a Poisson distribution with mean equal to the deterministic expectation of state transitions: e.g., at each time step, the number of deaths is obtained from a Poisson distribution with mean μ × N , where N is the total population size, after which individuals are randomly selected for a death event. The time unit is a day.
The original framework is expanded by including host-population structure and host movement within a meta-population formulation as in our previously published IBMs [41,42,43]. Although we model a closed population, we allow for random introductions of any strain under a fixed frequency 1 / η ; i.e., with probability η per time step a random host and strain are selected and colonization is attempted (which depends on the host’s epidemiological state as in Figure 1A,B). In summary, population structure is included by subdividing individual hosts into a spatially organised set of L × L communities of population size N c , forming a squared lattice of L columns and L rows. Individuals are assumed to mix homogeneously within each community. Local movement is considered, such that the force of infection of a strain in a particular community is dependent on its carriage levels within that community and in the neighbouring communities. Neighbouring communities are defined as the ones directly up, down, left, and right. Long distance host mobility is also modelled and assumed to take place randomly across the meta-population. Host movement is modelled by allowing individuals to transmit to individuals at a different community with probability ω by selecting the source and target communities and individuals randomly; that is, every time step a Poisson distributed number of individuals, with mean ω × y i , j , will attempt to transmit strain i , j to a random individual in the meta-population. The modelled processes related to population structure are summarised in Figure 1C. The model parametrization is given in the Model parameterization section of the Supplementary Information.

2.2. Antibiotic Resistance and Consumption Data

Data source used was the European Centre for Disease Prevention and Control (ECDC) from the surveillance atlas (antimicrobial resistance of invasive isolates reported to the EARS-NET) for the year 2005, vastly representative of a pre-PCV era across EU/EEA European countries (available at [44]). The year 2005 is the earliest made available by the ECDC. The data variables used were (i) the percent of samples sensitive to penicillins and macrolides (independently), and (ii) consumption defined by daily doses per 1000 people per day. We considered solely countries which had a number of samples larger than 50, which included 19 countries. Resistance level was aggregated across countries (EU) using a logistic regression model with random-effect intercepts for each country. Raw data on resistance per country and aggregated EU is presented in Supplementary Figures S1 and S2; for consumption in Supplementary Figure S3.

2.3. Approximate Bayesian Computation

To compare the model output to real-world estimates on the prevalence of antibiotic resistance, we employed approximate Bayesian computation (ABC). This method allowed us to approximate the posterior density of simulation parameters related to competition ( ψ , γ ), given priors on the parameters and data, when there is no explicit likelihood function [45], as is the case for our IBM framework. The priors used for the parameters ( ψ , γ ) were both uniformly distributed in [0.1, 0.9]. These priors were chosen for two reasons: First, we are not aware of any specific prior knowledge of these values - except that they are probably not close to zero or one. Second, the model run time is non-negligible, and is substantially increased when either of the competition parameters is close to zero or one. Hence, constrained by the relatively long computation time the model required, we sampled the prior a total of 6741 times (each with a unique combination of ψ and γ ) . We used three data points as “targets” for the ABC: (i) the observed proportion of resistance to penicillins or macrolides aggregated at the EU level; (ii) the observed relative proportion of NVT in some EU-related epidemiological contexts (Supplementary Figure S4); (iii) the proportion of carriers with co-infection. The ABC distance function employed to our model was the Euclidean distance from the simulation results to targets (i–iii). The threshold ϵ for acceptance under the ABC framework was set to either 0.025 or 0.05, so that the closest 2.5 or 5% of simulations to the targets (i–iii) would not be rejected. Further details of the implementation of the ABC approach are given in the section Approximate Bayesian computation of the Supplementary Information.

3. Results

3.1. Pneumococcal Strain Dynamics

We first assessed average model behaviour with regards to the levels of ecological and immunological competition, for which the parameters ψ and γ were marginally varied across the ranges [0.1, 0.9]. We focused on average model behaviour under the default parameter set (see Model parameterization in the Supplementary Information), but alternative sets including examples of temporal dynamics are presented in Supplementary Figures S5–S8.
When immunological competition γ was high and above γ 0.6 , independently of the level of ecological competition ψ , model dynamics were generally incompatible with observed pneumococcal strain dynamics. In particular, carriage levels of VT and NVT serotypes tended to be similar (Figure 2A), there was no host co-infection (Figure 2B), and strains resistant to antibiotics dominated the population (Figure 2C). Such epidemiological dynamics are incompatible with observed pneumococcal dynamics, which ubiquitously include high co-infection [46,47], reinfection [48,49,50], and co-circulation of NVT and VT serotypes with dominance of VT [10,47,48,51,52]. This dynamic behaviour is commonly observed in multi-locus, multi-strain pathogen modelling frameworks such as the one used here, in which high immunological competition has been shown to universally lead to the emergence of discrete strain structures (e.g., [28]).
In contrast, when immunological competition γ was low to intermediate, the model better replicated known properties of pneumococcal strain dynamics, resulting in dominance of VT serotypes and strains sensitive to antibiotics, as well as allowing for host coinfection (Figure 2A–C). Generally, for low-intermediate γ (i.e., 0.35 0.55 ), the carriage level of NVT serotypes was extremely low (Figure 2A); there was no co-infection (Figure 2B); while strains resistant to antibiotics did not dominate the population, their prevalence generally decreased for higher ecological competition ψ (Figure 2C). In contrast, when γ was lower (i.e., < 0.35 ), carriage levels of NVT serotypes increased with ecological competition (Figure 2A); host coinfection was possible with varying levels of ecological competition (Figure 2B); both coinfection and sensitivity to antibiotics decreased with increasing ecological competition (Figure 2B,C). These general behaviours across the competition parameter space ( ψ , γ ) did not change substantially in sensitivity analyses with different host-population structure and mobility values (Supplementary Figures S7 and S8). This suggests that variation of demographic parameters, which often exists between different regions, does not play a significant role in changing the levels of co-existence between susceptible and resistant strains in our modelling framework.
Resorting to data representative of the pre-vaccination era in Europe (from ECDC, see the Antibiotic resistance and consumption data section), we explored under what ecological and immunological competition levels the model better approximated observed levels and variation of resistance to macrolides and penicillins (independently). According to ECDC data, aggregated resistance across EU countries was approximately 18.1% (95% CI 12.19–26.51%) for macrolides and 8.1% (95% CI 5.02–12.95%) for penicillins (Supplementary Figures S1 and S2). While a wide range of γ and ψ values appeared compatible with these levels of resistance (Figure 2C), some of those combinations would nonetheless result in widely different levels of the NVT versus VT prevalence (Figure 2A) and levels of co-infection (Figure 2B). Using an ABC algorithm applied to a large number of numerical simulations exploring the model’s parameter space, we next set out to identify which combinations of γ and ψ would reproduce the observed levels of resistance while respecting prior knowledge on the relative proportions of NVT versus VT and host co-infection (see Materials and Methods section for details).

3.2. Strain Competition as Determinants of Antibiotic Resistance

As illustrated in Figure 3A,B, reproduction of pre-vaccination resistance levels, while respecting priors related to observed NVT versus VT prevalence and strain co-infection, identified particular competition parameter regions. The latter are characterized by a combination of low immunological ( γ < 0.25 ) and low to intermediate ecological competition ( ψ > 0.2 ,   ψ < 0.7 ), which are compatible with observations of recurrent pneumococcal infection [48,49,50], co-circulation of NVT and VT serotypes [10,47,48,51,52] as well as frequent co-infection [46,47].
For penicillins, a restricted range for immunological competition γ (0.12–0.18) and wide range of ecological competition ψ (0.3–0.7) were compatible with observed resistance levels (Figure 3C). In contrast, resistance to macrolides was compatible with the model output under a wider range of γ (0.1–0.25) and stricter and lower range of ψ (0.22–0.38). The intersection of the two compatible parameter regions for macroIides and penicillins was for γ 0.175 and ψ 0.375 . Doubling the percentage of accepted simulations under the ABC framework ( ϵ   =   5 % ) or the target co-infection level (from 0.2 to 0.3) had no significant impact on the identified parameter regions (Supplementary Figures S9 and S10). In contrast, accepting higher margins of error (e.g., ϵ   =   10 % ) resulted in the ABC identifying substantially wider regions of the parameter space as compatible, thus making it difficult to infer relationships between the model parameterization and observed resistance levels.
We further examined the associations between the competition parameters and resistance to penicillins and macrolides, performing an ABC which utilized each EU country’s antibiotic resistance frequencies to each of the two mentioned antibiotic classes (see Supplementary Information for method details). For either of the antibiotic classes, ecological competition ψ explained a higher percent of the variance of resistance than immunological competition γ in a simple linear model (penicillins, R 2 of 0.9 versus 0.001; macrolides, R 2 of 0.71 versus 0.52, respectively). Furthermore, for either of the antibiotic classes, a likelihood ratio test concluded that a bivariate model including both γ and ψ did not have a significantly better goodness of fit over a model containing ψ alone. Hence, in our model, variations in ecological competition could better explain differences of antibiotic resistance frequencies in the pre-vaccination era.

4. Conclusions

How do ecological and evolutionary forces drive the coexistence of pathogen strains susceptible and resistant to treatments or prophylaxis? This question remains a conundrum of high public health concern. Examples of such coexistence include some of the most relevant human pathogens: antiviral resistance of the immunodeficiency virus (HIV) [53,54], antimalarial resistance to malaria parasites [55], drug and vaccine escape by Mycobacterium tuberculosis [56], and antimicrobial resistance in Streptococcus pneumoniae [19,33,57,58].
Specifically, the understanding of what allows resistant bacteria to persist in human populations with different epidemiological, demographic, antibiotic consumption and epidemiological control conditions remains elusive. Antibiotic consumption is considered a major driver of resistance at the population level for both commensal and pathogenic bacteria [57,59,60]. However, while antibiotic consumption undoubtedly selects for antibiotic resistance, a major fraction of the observed spatio-temporal variation in resistance prevalence remains to be explained. Other population characteristics, e.g., gross domestic product (GDP) per capita or public health-care spending, have also been observed to affect variation in resistance levels [61]. Indeed, assuming no substantial time lag between consumption and prevalence of resistance [62], the ECDC data here explored shows that country-specific resistance is only mildly correlated with antibiotic consumption across countries, with consumption explaining 0.14 ( R 2 ) of the variation for macrolide resistance ( p   =   0.06 ) and 0.38 of penicillin resistance ( p   =   0.002 ) (Supplementary Figure S3).
Mathematical models have the potential to explore possible mechanistic drivers of resistance by explicitly considering factors modulating selective pressures. In this context, most pneumococcal modelling frameworks commonly overlook within-host interactions between phenotypically divergent strains (although for exceptions see [21,63]). The concept of “no coexistence for free” is also key for interpretation of modelling output. That is, candidate models cannot include intrinsic mechanisms allowing for stable coexistence of indistinguishable strains. Only under the assumption of such neutral models is it possible to directly discern the effect of the explicitly modelled factors in defining the existence of a stable state characterised by coexistence [64]. Indeed, the model presented here does not offer coexistence for free, since modelling indistinguishable strains ( γ   =   0 , ψ   =   0 ) dictates dynamics subjected only to drift (stochasticity). In this study, we expand a modelling framework previously used to explore within-host competition between pneumococcal strains as a driver of post-vaccination increases in antibiotic resistance [37]. Informed by available data for European countries, we expand the framework to a pre-vaccination era. We explore how the combination of within-host ecological and immunological competition may drive observed population-level coexistence between resistant and susceptible pneumococcal strains. Other possible drivers of coexistence (e.g., age-dependent assortative mixing [60]) were not explored in this study, although we have shown that high-level heterogeneities such as population structure and host mobility do not significantly affect model output (Supplementary Figures S6 and S7). Furthermore, our study did not consider double resistance to penicillins and macrolides. Unfortunately, the data we had access to did not contain the rates of double resistance to these antibiotics in the studied period. As such, we only model the marginal distributions of resistance to penicillins and macrolides. Although other data suggest that double resistance to macrolides and penicillins in this study’s setting are not high in pneumococci, they might not be negligible [38,39]. Future endeavors, informed by higher resolution epidemiological and molecular data on the interactions between the two resistance types, should also consider double resistance.
Our main results demonstrate the complex relationship between the frequency of NVT and VT pneumococci, co-infection, antibiotic resistance, and ecological and immunological competition factors. Based on data from European countries, our model indicates that pre-vaccination era antibiotic resistance levels are negatively associated with ecological competition. In our model, high ecological competition decreases antibiotic resistance levels, as this parameter governs the co-colonization competence of antibiotic resistant strains against susceptible strains, due to possible fitness costs incurred by acquisition of resistance [65]. The variation of ecological competition levels between susceptible and resistant strains across countries can be due to differences in both the bacterial and human populations. Different strains, with different genetic loads and, hence, resistance fitness costs, can expand across countries due to chance and founder effects [66,67]. Additionally, differences in human population density and demography, environmental conditions, per capita hospital beds, and other factors can substantially affect ecological competition and its influence on co-existence of different strains [68,69,70].
In contrast to the identified role of ecological competition in our model, we do not find that immunological competition significantly contributes to explaining the variation in observed resistance levels across European countries. In reality, immunological competition between strains is determined by a complex, largely unknown network of antigenic relationships between pneumococcal strains; with capsular polysaccharides believed to be the major determinant of strain-specific and cross-strain immunological responses [71]. As such, immunological competition levels would have to substantially differ between countries if they were a factor affecting rates of co-existence between resistant and susceptible strains. However, there is a lack of support for such substantial local differences in the antigenic determinants and overall prevalence levels of circulating antigenic groups [72,73]. Hence, the results of our model are compatible with the notion that population-level antigenic relationships and diversity are not dissimilar enough across countries to justify immunological competition as a major determinant of the co-existence of resistant and susceptible strains in the pre-vaccination era.
To conclude, we present a general yet explicit framework to model pneumococci, considering mechanisms of both between-host transmission and within-host competition as determinants of the pathogen’s population dynamics. In a series of simulated scenarios, we demonstrate that ecological competition during co-colonization between resistant and susceptible pneumococcal strains can explain much of the variation of resistance frequencies observed at the country level in the pre-vaccination era. Our results demonstrate the framework’s general potential to explore the drives of pneumococcal resistance dynamics. Elucidating the role of underlying competition-related components of pneumococcal dynamics improves the understanding of the mechanistic drivers for the emergence and maintenance of antibiotic resistance. Consequently, frameworks such as this can be used to project the influence of changes in vaccination or antibiotic treatment strategies on antibiotic resistance.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-393X/9/3/265/s1, Figure S1: ECDC resistance data for macrolides per country and aggregated estimation for EU, Figure S2: ECDC resistance data for penicillins per country and aggregated estimation for EU; Figure S3: ECDC resistance versus consumption data per country; Figure S4: Cross-country relative NVT frequency in the population; Figure S5: Examples of model output with differences in R0 between VT and NVT; Figure S6: Examples of model output with no differences in R0 between VT and NVT; Figure S7: Model strain dynamics under variations to ecological and immunological competition under variations of mobility; Figure S8: Model strain dynamics under variations to ecological and immunological competition under variations of population structure; Figure S9: Competition parameter space compatible with observed, pre-vaccination resistance levels when varying ABC sensitivity for co-infection target 0.2 and NVT ratio 0.15; Figure S10: Competition parameter space compatible with observed, pre-vaccination resistance levels when varying ABC sensitivity for co-infection target 0.3 and NVT ratio 0.15;

Author Contributions

J.L. and U.O. designed and implemented the model experiments. Y.D., A.G., J.L. and U.O. analysed model output. All authors wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by a lectureship by the Department of Zoology, University of Oxford (JL); The Raymond and Beverly Sackler Post-Doctoral Scholarship and a post-doctoral fellowship from the Tel Aviv University Center for Combating Pandemics (YD); NIHR Mucosal Pathogens Research Unit and the Bill and Melinda Gates Foundation (INC-004301—AG).

Data Availability Statement

All data are in the public domain.

Disclaimer

The views and opinions of the authors expressed herein do not necessarily state or reflect those of ECDC. The accuracy of the authors’ statistical analysis and the findings they report are not the responsibility of ECDC. ECDC is not responsible for conclusions or opinions drawn from the data provided. ECDC is not responsible for the correctness of the data and for data management, data merging and data collation after provision of the data. ECDC shall not be held liable for improper or incorrect use of the data.

Conflicts of Interest

Authors declare no conflict of interest.

References

  1. O’Brien, K.L.; Wolfson, L.J.; Watt, J.P.; Henkle, E.; Deloria-Knoll, M.; McCall, N.; Lee, E.; Mulholland, K.; Levine, O.S.; Cherian, T. Burden of disease caused by Streptococcus pneumoniae in children younger than 5 years: Global estimates. Lancet 2009, 374, 893–902. [Google Scholar] [CrossRef]
  2. Simell, B.; Auranen, K.; Käyhty, H.; Goldblatt, D.; Dagan, R.; O’Brien, K.L.; for the Pneumococcal Carriage Group. The fundamental link between pneumococcal car-riage and disease. Expert Rev. Vaccines 2012, 11, 841–855. [Google Scholar] [CrossRef] [Green Version]
  3. Tsaban, G.; Ben-Shimol, S. Indirect (herd) protection, following pneumococcal conjugated vaccines introduction: A systematic review of the literature. Vaccine 2017, 35, 2882–2891. [Google Scholar] [CrossRef] [PubMed]
  4. De Oliveira, L.H.; Camacho, L.A.B.; Coutinho, E.S.F.; Martinez-Silveira, M.S.; Carvalho, A.F.; Ruiz-Matus, C.; Toscano, C.M. Impact and Effec-tiveness of 10 and 13-Valent Pneumococcal Conjugate Vaccines on Hospitalization and Mortality in Children Aged Less than 5 Years in Latin American Countries: A Systematic Review. PLoS ONE 2016, 11, e0166736. [Google Scholar] [CrossRef] [PubMed]
  5. Hausdorff, W.P.; Feikin, D.R.; Klugman, K.P. Epidemiological differences among pneumococcal serotypes. Lancet Infect. Dis. 2005, 5, 83–93. [Google Scholar] [CrossRef]
  6. Bentley, S.D.; Aanensen, D.M.; Mavroidi, A.; Saunders, D.; Rabbinowitsch, E.; Collins, M.; Donohoe, K.; Harris, D.; Murphy, L.; Quail, M.A.; et al. Genetic Analysis of the Capsular Biosynthetic Locus from All 90 Pneumococcal Serotypes. PLoS Genet. 2006, 2, e31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Flasche, S.; Van Hoek, A.J.; Sheasby, E.; Waight, P.; Andrews, N.; Sheppard, C.; George, R.; Miller, E. Effect of Pneumococcal Conjugate Vaccination on Serotype-Specific Carriage and Invasive Disease in England: A Cross-Sectional Study. PLoS Med. 2011, 8, e1001017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Moore, M.R.; Link-Gelles, R.; Schaffner, W.; Lynfield, R.; Lexau, C.; Bennett, N.M.; Petit, S.; Zansky, S.M.; Harrison, L.H.; Reingold, A.; et al. Effect of use of 13-valent pneumococcal con-jugate vaccine in children on invasive pneumococcal disease in children and adults in the USA: Analysis of multisite, popula-tion-based surveillance. Lancet Infect. Dis. 2015, 15, 301–309. [Google Scholar] [CrossRef] [Green Version]
  9. Gidding, H.F.; Sheridan, S.; Fathima, P.; Moore, H.C.; Liu, B.; McIntyre, P.B.; Palmu, A.A. Impact of Childhood Pneumococcal Conjugate Vaccine on Nonnotified Clinically Suspected Invasive Pneumococcal Disease in Australia. Pediatr. Infect. Dis. J. 2019, 38, 860–865. [Google Scholar] [CrossRef] [Green Version]
  10. Lourenço, J.; Obolski, U.; Swarthout, T.D.; Gori, A.; Bar-Zeev, N.; Everett, D.; Kamng’Ona, A.W.; Mwalukomo, T.S.; Mataya, A.A.; Mwansambo, C.; et al. Determinants of high residual post-PCV13 pneumococcal vaccine-type carriage in Blantyre, Malawi: A modelling study. BMC Med. 2019, 17, 219. [Google Scholar] [CrossRef] [Green Version]
  11. Swarthout, T.D.; Fronterre, C.; Lourenço, J.; Obolski, U.; Gori, A.; Bar-Zeev, N.; Everett, D.; Kamng’Ona, A.W.; Mwalukomo, T.S.; Mataya, A.A.; et al. High residual carriage of vaccine-serotype Streptococcus pneumoniae after introduction of pneumococcal conjugate vaccine in Malawi. Nat. Commun. 2020, 11, 2222. [Google Scholar] [CrossRef]
  12. Moore, M.R.; Gertz, R.E., Jr.; Woodbury, R.L.; Barkocy-Gallagher, G.A.; Schaffner, W.; Lexau, C.; Gershman, K.; Reingold, A.; Farley, M.; Harrison, L.H.; et al. Population Snapshot of EmergentStreptococcus pneumoniaeSerotype 19A in the United States, 2005. J. Infect. Dis. 2008, 197, 1016–1027. [Google Scholar] [CrossRef] [Green Version]
  13. Hanage, W.P.; Huang, S.S.; Lipsitch, M.; Bishop, C.J.; Godoy, D.; Pelton, S.I.; Goldstein, R.; Huot, H.; Finkelstein, J.A. Diversity and antibiotic resistance among nonvac-cine serotypes of Streptococcus pneumoniae carriage isolates in the post-heptavalent conjugate vaccine era. J. Infect. Dis. 2007, 195, 347–352. [Google Scholar] [CrossRef] [Green Version]
  14. Kyaw, M.H.; Lynfield, R.; Schaffner, W.; Craig, A.S.; Hadler, J.; Reingold, A.; Thomas, A.R.; Harrison, L.H.; Bennett, N.M.; Farley, M.M.; et al. Effect of introduction of the pneumococcal conju-gate vaccine on drug-resistant Streptococcus pneumoniae. N. Engl. J. Med. 2006, 354, 1455–1463. [Google Scholar] [CrossRef] [PubMed]
  15. Park, S.Y.; Moore, M.R.; Bruden, D.L.; Hyde, T.B.; Reasonover, A.L.; Harker-Jones, M.; Rudolph, K.M.; Hurlburt, D.A.; Parks, D.J.; Parkinson, A.J. Impact of conjugate vaccine on transmis-sion of antimicrobial-resistant Streptococcus pneumoniae among Alaskan children. Pediatr. Infect. Dis. J. 2008, 27, 335–340. [Google Scholar] [CrossRef] [PubMed]
  16. Huang, S.S.; Platt, R.; Rifas-Shiman, S.L.; Pelton, S.I.; Goldmann, D.; Finkelstein, J.A. Post-PCV7 changes in colonizing pneumococcal serotypes in 16 Massachusetts communities, 2001 and 2004. Pediatrics 2005, 116, e408–e413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Gherardi, G.; D’Ambrosio, F.; Visaggio, D.; Dicuonzo, G.; Del Grosso, M.; Pantosti, A. Serotype and Clonal Evolution of Penicillin-Nonsusceptible Invasive Streptococcus pneumoniae in the 7-Valent Pneumococcal Conjugate Vaccine Era in Italy. Antimicrob. Agents Chemother. 2012, 56, 4965–4968. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Farrell, D.J.; Klugman, K.P.; Pichichero, M. Increased Antimicrobial Resistance among Nonvaccine Serotypes of Streptococcus pneumoniae in the Pediatric Population After the Introduction of 7-Valent Pneumococcal Vaccine in the United States. Pediatr. Infect. Dis. J. 2007, 26, 123–128. [Google Scholar] [CrossRef] [PubMed]
  19. Pillai, D.R.; Shahinas, D.; Buzina, A.; A Pollock, R.; Lau, R.; Khairnar, K.; Wong, A.; Farrell, D.J.; Green, K.; McGeer, A.; et al. Genome-wide dissection of globally emergent multi-drug resistant serotype 19A Streptococcus pneumoniae. BMC Genom. 2009, 10, 642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Lipsitch, M. Vaccination against colonizing bacteria with multiple serotypes. Proc. Natl. Acad. Sci. USA 1997, 94, 6571–6576. [Google Scholar] [CrossRef] [Green Version]
  21. Cobey, S.; Lipsitch, M. Niche and Neutral Effects of Acquired Immunity Permit Coexistence of Pneumococcal Serotypes. Science 2012, 335, 1376–1380. [Google Scholar] [CrossRef] [Green Version]
  22. Bottomley, C.; Roca, A.; Hill, P.C.; Greenwood, B.; Isham, V. A mathematical model of serotype replacement in pneumococcal car-riage following vaccination. J. R. Soc. Interface 2013, 10, 20130786. [Google Scholar] [CrossRef]
  23. Flasche, S.; Edmunds, W.J.; Miller, E.; Goldblatt, D.; Robertson, C.; Choi, Y.H. The impact of specific and non-specific immunity on the ecology of Streptococcus pneumoniae and the implications for vaccination. Proc. R. Soc. B Boil. Sci. 2013, 280, 20131939. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Melegaro, A.; Choi, Y.H.; George, R.; Edmunds, W.J.; Miller, E.; Gay, N.J. Dynamic models of pneumococcal carriage and the impact of the Heptavalent Pneumococcal Conjugate Vaccine on invasive pneumococcal disease. BMC Infect. Dis. 2010, 10, 90. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Beall, B.W.; Gertz, R.E.; Hulkower, R.L.; Whitney, C.G.; Moore, M.R.; Brueggemann, A.B. Shifting Genetic Structure of Invasive Serotype 19A Pneumococci in the United States. J. Infect. Dis. 2011, 203, 1360–1368. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Croucher, N.J.; Finkelstein, J.A.; Pelton, S.I.; Mitchell, P.K.; Lee, G.M.; Parkhill, J.; Bentley, S.D.; Hanage, W.P.; Lipsitch, M. Population genomics of post-vaccine changes in pneumococcal epidemiology. Nat. Genet. 2013, 45, 656–663. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Gupta, S.; Ferguson, N.M.; Anderson, R.M. Vaccination and the population structure of antigenically diverse pathogens that ex-change genetic material. Proc. Biol. Sci. 1997, 264, 1435–1443. [Google Scholar] [CrossRef] [Green Version]
  28. Watkins, E.R.; Penman, B.S.; Lourenço, J.; Buckee, C.O.; Maiden, M.C.J.; Gupta, S. Vaccination Drives Changes in Metabolic and Viru-lence Profiles of Streptococcus pneumoniae. PLoS Pathog. 2015, 11, e1005034. [Google Scholar] [CrossRef] [Green Version]
  29. Lourenço, J.; Watkins, E.R.; Obolski, U.; Peacock, S.J.; Morris, C.; Maiden, M.C.J.; Gupta, S. Lineage structure of Streptococcus pneumoniae may be driven by immune selection on the groEL heat-shock protein. Sci. Rep. 2017, 7, 9023. [Google Scholar] [CrossRef] [Green Version]
  30. Parks, T.; Barrett, L.; Jones, N. Invasive streptococcal disease: A review for clinicians. Br. Med Bull. 2015, 115, 77–89. [Google Scholar] [CrossRef] [Green Version]
  31. Dagan, R.; Givon-Lavi, N.; Zamir, O.; Fraser, D. Effect of a nonavalent conjugate vaccine on carriage of antibiotic-resistant Strep-tococcus pneumoniae in day-care centers. Pediatr. Infect. Dis. J. 2003, 22, 532–540. [Google Scholar] [CrossRef]
  32. Mbelle, N.; Huebner, R.E.; Wasas, A.D.; Kimura, A.; Chang, I.; Klugman, K.P. Immunogenicity and impact on nasopharyngeal car-riage of a nonavalent pneumococcal conjugate vaccine. J. Infect. Dis. 1999, 180, 1171–1176. [Google Scholar] [CrossRef] [Green Version]
  33. Mishra, R.P.N.; Oviedo-Orta, E.; Prachi, P.; Rappuoli, R.; Bagnoli, F. Vaccines and antibiotic resistance. Curr. Opin. Microbiol. 2012, 15, 596–602. [Google Scholar] [CrossRef] [PubMed]
  34. Chiba, N.; Morozumi, M.; Shouji, M.; Wajima, T.; Iwata, S.; Sunakawa, K.; Ubukata, A.; Invasive Pneumococcal Diseases Surveillance Study Group. Rapid Decrease of 7-Valent Conjugate Vaccine Coverage for Invasive Pneumococcal Diseases in Pediatric Patients in Japan. Microb. Drug Resist. 2013, 19, 308–315. [Google Scholar] [CrossRef]
  35. Dagan, R.; Klugman, K.P. Impact of conjugate pneumococcal vaccines on antibiotic resistance. Lancet Infect. Dis. 2008, 8, 785–795. [Google Scholar] [CrossRef]
  36. Danino, D.; Givon-Lavi, N.; Ben-Shimol, S.; Greenberg, D.; Dagan, R. Understanding the Evolution of Antibiotic-nonsusceptible Pneumococcal Nasopharyngeal Colonization Following Pneumococcal Conjugate Vaccine Implementation in Young Children. Clin. Infect. Dis. 2019, 69, 648–656. [Google Scholar] [CrossRef]
  37. Obolski, U.; Lourenço, J.; Thompson, C.; Thompson, R.; Gori, A.; Gupta, S. Vaccination can drive an increase in frequencies of anti-biotic resistance among nonvaccine serotypes of. Proc. Natl. Acad. Sci. USA 2018, 115, 3102–3107. [Google Scholar] [CrossRef] [Green Version]
  38. Reinert, R. The antimicrobial resistance profile of Streptococcus pneumoniae. Clin. Microbiol. Infect. 2009, 15, 7–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Annual Epidemiological Report 2009 [2007 Data]. European Centre for Disease Control [Internet]. 12 October 2009. Available online: https://www.ecdc.europa.eu/en/publications-data/annual-epidemiological-report-2009-2007-data (accessed on 21 September 2020).
  40. Willem, L.; Verelst, F.; Bilcke, J.; Hens, N.; Beutels, P. Lessons from a decade of individual-based models for infectious disease transmission: A systematic review (2006–2015). BMC Infect. Dis. 2017, 17, 612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Lourenço, J.; Recker, M. Natural, Persistent Oscillations in a Spatial Multi-Strain Disease System with Application to Dengue. PLoS Comput. Biol. 2013, 9, e1003308. [Google Scholar] [CrossRef] [Green Version]
  42. Lourenço, J.; Recker, M. Dengue serotype immune-interactions and their consequences for vaccine impact predictions. Epidemics 2016, 16, 40–48. [Google Scholar] [CrossRef] [Green Version]
  43. Flasche, S.; Jit, M.; Rodríguez-Barraquer, I.; Coudeville, L.; Recker, M.; Koelle, K.; Milne, G.; Hladish, T.J.; Perkins, T.A.; Cummings, D.A.T.; et al. The Long-Term Safety, Public Health Impact, and Cost-Effectiveness of Routine Vaccination with a Recombinant, Live-Attenuated Dengue Vaccine (Dengvaxia): A Model Comparison Study. PLoS Med. 2016, 13, e1002181. [Google Scholar] [CrossRef] [Green Version]
  44. European Centre for Disease Control (ECDC) Antibiotic Resistance Data. ECDC [Internet]. Available online: https://www.ecdc.europa.eu/en/antimicrobial-resistance (accessed on 22 July 2020).
  45. Minter, A.; Retkute, R. Approximate Bayesian Computation for infectious disease modelling. Epidemics 2019, 29, 100368. [Google Scholar] [CrossRef] [PubMed]
  46. Tabatabaei, S.; Karimi, A.; Fallah, F.; Shiva, F.; Shamshiri, M.; Gooya, M.; Zahraei, M. Rate of co-colonization with serotypes of strep pneumonia isolated from nasopharyngeal swab. Int. J. Infect. Dis. 2014, 21, 418. [Google Scholar] [CrossRef] [Green Version]
  47. Kamng’Ona, A.W.; Hinds, J.; Bar-Zeev, N.; Gould, K.A.; Chaguza, C.; Msefula, C.; Cornick, J.E.; Kulohoma, B.W.; Gray, K.; Bentley, S.D.; et al. High multiple carriage and emergence of Streptococcus pneumoniae vaccine serotype variants in Malawian children. BMC Infect. Dis. 2015, 15, 234. [Google Scholar] [CrossRef] [Green Version]
  48. Melegaro, A.; Choi, Y.; Pebody, R.; Gay, N. Pneumococcal Carriage in United Kingdom Families: Estimating Serotype-specific Transmission Parameters from Longitudinal Data. Am. J. Epidemiol. 2007, 166, 228–235. [Google Scholar] [CrossRef] [PubMed]
  49. Melegaro, A.; Gay, N.J.; Medley, G.F. Estimating the transmission parameters of pneumococcal carriage in households. Epidemiol. Infect. 2004, 132, 433–441. [Google Scholar] [CrossRef]
  50. Lipsitch, M.; Abdullahi, O.; DʼAmour, A.; Xie, W.; Weinberger, D.M.; Tchetgen Tchetgen, E.; Scott, J.A.G. Estimating rates of carriage acqui-sition and clearance and competitive ability for pneumococcal serotypes in Kenya with a Markov transition model. Epidemiology 2012, 23, 510–519. [Google Scholar] [CrossRef] [Green Version]
  51. Adetifa, I.M.O.; Antonio, M.; Okoromah, C.A.N.; Ebruke, C.; Inem, V.; Nsekpong, D.; Bojang, A.; Adegbola, R.A. Pre-Vaccination Nasopharyngeal Pneumococcal Carriage in a Nigerian Population: Epidemiology and Population Biology. PLoS ONE 2012, 7, e30548. [Google Scholar] [CrossRef]
  52. Heinsbroek, E.; Tafatatha, T.; Phiri, A.; Swarthout, T.D.; Alaerts, M.; Crampin, A.C.; Chisambo, C.; Mwiba, O.; Read, J.M.; French, N. Pneumococcal carriage in households in Karonga District, Malawi, before and after introduction of 13-valent pneumococcal conjugate vaccination. Vaccine 2018, 36, 7369–7376. [Google Scholar] [CrossRef] [PubMed]
  53. Gupta, R.K.; Gregson, J.; Parkin, N.; Haile-Selassie, H.; Tanuri, A.; Forero, L.A.; Kaleebu, P.; Watera, C.; Aghokeng, A.; Mutenda, N.; et al. HIV-1 drug resistance before initiation or re-initiation of first-line antiretroviral therapy in low-income and middle-income countries: A systematic review and meta-regression analysis. Lancet Infect. Dis. 2018, 18, 346–355. [Google Scholar] [CrossRef] [Green Version]
  54. Puertas, M.C.; Ploumidis, G.; Ploumidis, M.; Fumero, E.; Clotet, B.; Walworth, C.M.; Petropoulos, C.J.; Martinez-Picado, J. Pan-resistant HIV-1 emergence in the era of integrase strand-transfer inhibitors: A case report. Lancet Microbe 2020, 1, e130–e135. [Google Scholar] [CrossRef]
  55. Gao, B.; Saralamba, S.; Lubell, Y.; White, L.J.; Dondorp, A.M.; Aguas, R. Determinants of MDA impact and designing MDAs towards malaria elimination. eLife 2020, 9. [Google Scholar] [CrossRef] [Green Version]
  56. Kremer, K.; Van Der Werf, M.J.; Au, B.K.Y.; Anh, D.D.; Kam, K.M.; Van Doorn, H.R.; Borgdorff, M.W.; Van Soolingen, D. Vaccine-induced Immunity Circumvented by Typical Mycobacterium tuberculosis Beijing Strains. Emerg. Infect. Dis. 2009, 15, 335–339. [Google Scholar] [CrossRef]
  57. Austin, D.J.; Kristinsson, K.G.; Anderson, R.M. The relationship between the volume of antimicrobial consumption in human communities and the frequency of resistance. Proc. Natl. Acad. Sci. USA 1999, 96, 1152–1156. [Google Scholar] [CrossRef] [Green Version]
  58. Lehtinen, S.; Blanquart, F.; Croucher, N.J.; Turner, P.; Lipsitch, M.; Fraser, C. Evolution of antibiotic resistance is linked to any genet-ic mechanism affecting bacterial duration of carriage. Proc. Natl. Acad. Sci. USA 2017, 114, 1075–1080. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Colijn, C.; Cohen, T.; Fraser, C.; Hanage, W.; Goldstein, E.; Givon-Lavi, N.; Dagan, R.; Lipsitch, M. What is the mechanism for persistent coexistence of drug-susceptible and drug-resistant strains of Streptococcus pneumoniae? J. R. Soc. Interface 2009, 7, 905–919. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Cobey, S.; Baskerville, E.B.; Colijn, C.; Hanage, W.; Fraser, C.; Lipsitch, M. Host population structure and treatment frequency maintain balancing selection on drug resistance. J. R. Soc. Interface 2017, 14, 20170295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Collignon, P.; Beggs, J.J.; Walsh, T.R.; Gandra, S.; Laxminarayan, R. Anthropological and socioeconomic factors contributing to global antimicrobial resistance: A univariate and multivariable analysis. Lancet Planet. Health 2018, 2, e398–e405. [Google Scholar] [CrossRef]
  62. Blanquart, F.; Lehtinen, S.; Fraser, C. An evolutionary model to predict the frequency of antibiotic resistance under seasonal antibiotic use, and an application to Streptococcus pneumoniae. Proc. R. Soc. B Boil. Sci. 2017, 284, 20170679. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Davies, N.G.; Flasche, S.; Jit, M.; Atkins, K.E. Within-host dynamics shape antibiotic resistance in commensal bacteria. Nat. Ecol. Evol. 2019, 3, 440–449. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Lipsitch, M.; Colijn, C.; Cohen, T.; Hanage, W.P.; Fraser, C. No coexistence for free: Neutral null models for multistrain pathogens. Epidemics 2009, 1, 2–13. [Google Scholar] [CrossRef]
  65. Melnyk, A.H.; Wong, A.; Kassen, R. The fitness costs of antibiotic resistance mutations. Evol. Appl. 2015, 8, 273–283. [Google Scholar] [CrossRef] [PubMed]
  66. Janoir, C.; Cohen, R.; Levy, C.; Bingen, E.; Lepoutre, A.; Gutmann, L.; Varon, E.; the Observatoires Régionaux du Pneumocoque (ORP) Network. Clonal Expansion of the Macrolide Resistant ST386 within Pneumococcal Serotype 6C in France. PLoS ONE 2014, 9, e90935. [Google Scholar] [CrossRef] [PubMed]
  67. Keenan, J.D.; Klugman, K.P.; McGee, L.; Vidal, J.E.; Chochua, S.; Hawkins, P.; Cevallos, V.; Gebre, T.; Tadesse, Z.; Emerson, P.M.; et al. Evidence for clonal expansion after antibiotic se-lection pressure: Pneumococcal multilocus sequence types before and after mass azithromycin treatments. J. Infect. Dis. 2015, 211, 988–994. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. MacFadden, D.R.; McGough, S.F.; Fisman, D.; Santillana, M.; Brownstein, J.S. Antibiotic Resistance Increases with Local Tempera-ture. Nat. Clim. Chang. 2018, 8, 510–514. [Google Scholar] [CrossRef]
  69. Bruinsma, N.; Hutchinson, J.M.; van den Bogaard, A.E.; Giamarellou, H.; Degener, J.; Stobberingh, E.E. Influence of population densi-ty on antibiotic resistance. J. Antimicrob. Chemother. 2003, 51, 385–390. [Google Scholar] [CrossRef] [Green Version]
  70. Blommaert, A.; Marais, C.; Hens, N.; Coenen, S.; Muller, A.; Goossens, H.; Beutels, P. Determinants of between-country differences in ambulatory antibiotic use and antibiotic resistance in Europe: A longitudinal observational study. J. Antimicrob. Chemother. 2013, 69, 535–547. [Google Scholar] [CrossRef] [Green Version]
  71. Brown, J.; Hammerschmidt, S.; Orihuela, C. Streptococcus Pneumoniae: Molecular Mechanisms of Host-Pathogen Interactions; Academic Press: Cambridge, MA, USA, 2015. [Google Scholar]
  72. Jefferson, T.; Ferroni, E.; Curtale, F.; Rossi, P.G.; Borgia, P. Streptococcus pneumoniae in western Europe: Serotype distribution and incidence in children less than 2 years old. Lancet Infect. Dis. 2006, 6, 405–410. [Google Scholar] [CrossRef]
  73. Feikin, D.R.; Klugman, K.P. Historical Changes in Pneumococcal Serogroup Distribution: Implications for the Era of Pneumo-coccal Conjugate Vaccines. Clin. Infect. Dis. 2002, 35, 547–555. [Google Scholar] [CrossRef]
Figure 1. Model transmission assumptions and host-population structure. (a) Serotype-specific immunity as modelled by the parameter γ . Immunity can hinder the acquisition of a strain of a particular serotype (e.g., i   =   a ) if (top) the receiving host is currently carrying a strain of that serotype ( i   =   a ), or (bottom) the receiving host is clear of infection (carriage) but has previously carried a strain of that serotype ( i   =   a ). (b) Resistance-related, within-host competition as modelled by the parameter ψ . Competition from an already established sensitive strain ( j   =   00 ) can hinder the acquisition of a strain carrying resistance to antimicrobials ( j   =   01 or j   =   10 ). (a,b) The effects of immunity and competition are not mutually exclusive, and can take place at the same time if conditions shown in subplot a (top) and subplot b are met. (c) Schematic of the meta-population used, with L × L communities of size N c hosts, with frequent mixing between adjacent communities, and sporadic events between distant communities with probability modelled by the parameter ω .
Figure 1. Model transmission assumptions and host-population structure. (a) Serotype-specific immunity as modelled by the parameter γ . Immunity can hinder the acquisition of a strain of a particular serotype (e.g., i   =   a ) if (top) the receiving host is currently carrying a strain of that serotype ( i   =   a ), or (bottom) the receiving host is clear of infection (carriage) but has previously carried a strain of that serotype ( i   =   a ). (b) Resistance-related, within-host competition as modelled by the parameter ψ . Competition from an already established sensitive strain ( j   =   00 ) can hinder the acquisition of a strain carrying resistance to antimicrobials ( j   =   01 or j   =   10 ). (a,b) The effects of immunity and competition are not mutually exclusive, and can take place at the same time if conditions shown in subplot a (top) and subplot b are met. (c) Schematic of the meta-population used, with L × L communities of size N c hosts, with frequent mixing between adjacent communities, and sporadic events between distant communities with probability modelled by the parameter ω .
Vaccines 09 00265 g001
Figure 2. Model strain dynamics under variations to ecological and immunological competition. (A) Relative frequency (ratio) of total number of individuals carrying NVT versus carrying any type (NVT + VT). (B) Proportion of infected hosts carrying more than one strain (co-infection). (C) Relative frequency (ratio) of total number of individuals carrying resistant strains to one antibiotic (VT j = 10, VT j = 01, NVT j = 10, NVT j = 01) and carrying any strain. (A–C) All model parameters as in default parameter set except ψ , γ , varied in the y and x axis, respectively. Results presented are the mean over the last 5 years of a simulation, for particular combinations of ψ , γ .
Figure 2. Model strain dynamics under variations to ecological and immunological competition. (A) Relative frequency (ratio) of total number of individuals carrying NVT versus carrying any type (NVT + VT). (B) Proportion of infected hosts carrying more than one strain (co-infection). (C) Relative frequency (ratio) of total number of individuals carrying resistant strains to one antibiotic (VT j = 10, VT j = 01, NVT j = 10, NVT j = 01) and carrying any strain. (A–C) All model parameters as in default parameter set except ψ , γ , varied in the y and x axis, respectively. Results presented are the mean over the last 5 years of a simulation, for particular combinations of ψ , γ .
Vaccines 09 00265 g002
Figure 3. Competition parameter space compatible with observed, pre-vaccination resistance levels. Approximate Bayesian Computation (ABC) output when varying immunological ( γ ) and ecological ( ψ ) competition, attempting to reproduce observed levels of resistance to penicillins (A) and macrolides (B) in the European region (see Data section for details). The colour scale is the density in ABC posterior estimate when applied to the aggregated European region resistance levels. Ellipses mark the 50 (dotted) and 95 (full) percentiles in ABC output. ABC priors and targets as detailed in the Methods section, the number of simulations was 6741, ϵ was set to 2.5%. (C) Ellipses marking the percentiles in panels A (green) and B (purple).
Figure 3. Competition parameter space compatible with observed, pre-vaccination resistance levels. Approximate Bayesian Computation (ABC) output when varying immunological ( γ ) and ecological ( ψ ) competition, attempting to reproduce observed levels of resistance to penicillins (A) and macrolides (B) in the European region (see Data section for details). The colour scale is the density in ABC posterior estimate when applied to the aggregated European region resistance levels. Ellipses mark the 50 (dotted) and 95 (full) percentiles in ABC output. ABC priors and targets as detailed in the Methods section, the number of simulations was 6741, ϵ was set to 2.5%. (C) Ellipses marking the percentiles in panels A (green) and B (purple).
Vaccines 09 00265 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lourenço, J.; Daon, Y.; Gori, A.; Obolski, U. Pneumococcal Competition Modulates Antibiotic Resistance in the Pre-Vaccination Era: A Modelling Study. Vaccines 2021, 9, 265. https://doi.org/10.3390/vaccines9030265

AMA Style

Lourenço J, Daon Y, Gori A, Obolski U. Pneumococcal Competition Modulates Antibiotic Resistance in the Pre-Vaccination Era: A Modelling Study. Vaccines. 2021; 9(3):265. https://doi.org/10.3390/vaccines9030265

Chicago/Turabian Style

Lourenço, José, Yair Daon, Andrea Gori, and Uri Obolski. 2021. "Pneumococcal Competition Modulates Antibiotic Resistance in the Pre-Vaccination Era: A Modelling Study" Vaccines 9, no. 3: 265. https://doi.org/10.3390/vaccines9030265

APA Style

Lourenço, J., Daon, Y., Gori, A., & Obolski, U. (2021). Pneumococcal Competition Modulates Antibiotic Resistance in the Pre-Vaccination Era: A Modelling Study. Vaccines, 9(3), 265. https://doi.org/10.3390/vaccines9030265

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop