What is the mechanism for persistent coexistence of drug-susceptible and drug-resistant strains of Streptococcus pneumoniae?

Recent approaches in computational modelling for controlling pathogen threats

In this review, we assess the status of computational modelling of pathogens. We focus on three disparate but interlinked research areas that produce models with very different spatial and temporal scope. First, we examine antimicrobial resistance (AMR). Many mechanisms of AMR are not well understood. As a result, it is hard to measure the current incidence of AMR, predict the future incidence, and design strategies to preserve existing antibiotic effectiveness. Next, we look at how to choose the finite number of bacterial strains that can be included in a vaccine. To do this, we need to understand what happens to vaccine and non-vaccine strains after vaccination programmes. Finally, we look at within-host modelling of antibody dynamics. The SARS-CoV-2 pandemic produced huge amounts of antibody data, prompting improvements in this area of modelling. We finish by discussing the challenges that persist in understanding these complex biological systems.

When does antimicrobial resistance increase bacterial fitness? Effects of dosing, social interactions, and frequency dependence on the benefits of AmpC β-lactamases in broth, biofilms, and a gut infection model

One of the longstanding puzzles of antimicrobial resistance is why the frequency of resistance persists at intermediate levels. Theoretical explanations for the lack of fixation of resistance include cryptic costs of resistance or negative frequency-dependence but are seldom explored experimentally. β-lactamases, which detoxify penicillin-related antibiotics, have well-characterized frequency-dependent dynamics driven by cheating and cooperation. However, bacterial physiology determines whether β-lactamases are cooperative, and we know little about the sociality or fitness of β-lactamase producers in infections. Moreover, media-based experiments constrain how we measure fitness and ignore important parameters such as infectivity and transmission among hosts. Here, we investigated the fitness effects of broad-spectrum AmpC β-lactamases in Enterobacter cloacae in broth, biofilms, and gut infections in a model insect. We quantified frequency- and dose-dependent fitness using cefotaxime, a third-generation cephalosporin. We predicted that infection dynamics would be similar to those observed in biofilms, with social protection extending over a wide dose range. We found evidence for the sociality of β-lactamases in all contexts with negative frequency-dependent selection, ensuring the persistence of wild-type bacteria, although cooperation was less prevalent in biofilms, contrary to predictions. While competitive fitness in gut infections and broth had similar dynamics, incorporating infectivity into measurements of fitness in infections significantly affected conclusions. Resistant bacteria had reduced infectivity, which limited the fitness benefits of resistance to infections challenged with low antibiotic doses and low initial frequencies of resistance. The fitness of resistant bacteria in more physiologically tolerant states (in biofilms, in infections) could be constrained by the presence of wild-type bacteria, high antibiotic doses, and limited availability of β-lactamases. One conclusion is that increased tolerance of β-lactams does not necessarily increase selection pressure for resistance. Overall, both cryptic fitness costs and frequency dependence curtailed the fitness benefits of resistance in this study.

Exploring factors shaping antibiotic resistance patterns in Streptococcus pneumoniae during the 2020 COVID-19 pandemic

Non-pharmaceutical interventions implemented to block SARS-CoV-2 transmission in early 2020 led to global reductions in the incidence of invasive pneumococcal disease (IPD). By contrast, most European countries reported an increase in antibiotic resistance among invasive Streptococcus pneumoniae isolates from 2019 to 2020, while an increasing number of studies reported stable pneumococcal carriage prevalence over the same period. To disentangle the impacts of the COVID-19 pandemic on pneumococcal epidemiology in the community setting, we propose a mathematical model formalizing simultaneous transmission of SARS-CoV-2 and antibiotic-sensitive and -resistant strains of S. pneumoniae. To test hypotheses underlying these trends five mechanisms were built into the model and examined: (1) a population-wide reduction of antibiotic prescriptions in the community, (2) lockdown effect on pneumococcal transmission, (3) a reduced risk of developing an IPD due to the absence of common respiratory viruses, (4) community azithromycin use in COVID-19 infected individuals, (5) and a longer carriage duration of antibiotic-resistant pneumococcal strains. Among 31 possible pandemic scenarios involving mechanisms individually or in combination, model simulations surprisingly identified only two scenarios that reproduced the reported trends in the general population. They included factors (1), (3), and (4). These scenarios replicated a nearly 50% reduction in annual IPD, and an increase in antibiotic resistance from 20% to 22%, all while maintaining a relatively stable pneumococcal carriage. Exploring further, higher SARS-CoV-2 R0 values and synergistic within-host virus-bacteria interaction mechanisms could have additionally contributed to the observed antibiotic resistance increase. Our work demonstrates the utility of the mathematical modeling approach in unraveling the complex effects of the COVID-19 pandemic responses on AMR dynamics.

Closing the gap in the Janzen-Connell hypothesis: What determines pathogen diversity?

The high tree diversity in tropical forests has long been a puzzle to ecologists. In the 1970s, Janzen and Connell proposed that tree species (hosts) coexist due to the stabilizing actions of specialized enemies. This Janzen-Connell hypothesis was subsequently supported by theoretical studies. Yet, such studies have taken the presence of specialized pathogens for granted, overlooking that pathogen coexistence also requires an explanation. Moreover, stable ecological coexistence does not necessarily imply evolutionary stability. What are the conditions that allow Janzen-Connell effects to evolve? We link theory from community ecology, evolutionary biology and epidemiology to tackle this question, structuring our approach around five theoretical frameworks. Phenomenological Lotka-Volterra competition models provide the most basic framework, which can be restructured to include (single- or multi-)pathogen dynamics. This ecological foundation can be extended to include pathogen evolution. Hosts, of course, may also evolve, and we introduce a coevolutionary model, showing that host-pathogen coevolution can lead to highly diverse systems. Our work unpacks the assumptions underpinning Janzen-Connell and places theoretical bounds on pathogen and host ecology and evolution. The five theoretical frameworks taken together provide a stronger theoretical basis for Janzen-Connell, delivering a wider lens that can yield important insights into the maintenance of diversity in these increasingly threatened systems.

Challenges in Forecasting Antimicrobial Resistance

Antimicrobial resistance is a major threat to human health. Since the 2000s, computational tools for predicting infectious diseases have been greatly advanced; however, efforts to develop real-time forecasting models for antimicrobial-resistant organisms (AMROs) have been absent. In this perspective, we discuss the utility of AMRO forecasting at different scales, highlight the challenges in this field, and suggest future research priorities. We also discuss challenges in scientific understanding, access to high-quality data, model calibration, and implementation and evaluation of forecasting models. We further highlight the need to initiate research on AMRO forecasting using currently available data and resources to galvanize the research community and address initial practical questions.

Interplay between strain fitness and transmission frequency determines prevalence of antimicrobial resistance

The steep rise of infections caused by bacteria that are resistant to antimicrobial agents threatens global health. However, the association between antimicrobial use and the prevalence of resistance is not straightforward. Therefore, it is necessary to quantify the importance of additional factors that affect this relationship. We theoretically explore how the prevalence of resistance is affected by the combination of three factors: antimicrobial use, bacterial transmission, and fitness cost of resistance. We present a model that combines within-host, between-hosts and between-populations dynamics, built upon the competitive Lotka-Volterra equations. We developed the model in a manner that allows future experimental validation of the findings with single isolates in the laboratory. Each host may carry two strains (susceptible and resistant) that represent the host’s commensal microbiome and are not the target of the antimicrobial treatment. The model simulates a population of hosts who are treated periodically with antibiotics and transmit bacteria to each other. We show that bacterial transmission results in strain co-existence. Transmission disseminates resistant bacteria in the population, increasing the levels of resistance. Counterintuitively, when the cost of resistance is low, high transmission frequencies reduce resistance prevalence. Transmission between host populations leads to more similar resistance levels, increasing the susceptibility of the population with higher antimicrobial use. Overall, our results indicate that the interplay between bacterial transmission and strain fitness affects the prevalence of resistance in a non-linear way. We then place our results within the context of ecological theory, particularly on temporal niche partitioning and metapopulation rescue, and we formulate testable experimental predictions for future research.

The potential of genomics for infectious disease forecasting

Genomic technologies have led to tremendous gains in understanding how pathogens function, evolve and interact. Pathogen diversity is now measurable at high precision and resolution, in part because over the past decade, sequencing technologies have increased in speed and capacity, at decreased cost. Alongside this, the use of models that can forecast emergence and size of infectious disease outbreaks has risen, highlighted by the coronavirus disease 2019 pandemic but also due to modelling advances that allow for rapid estimates in emerging outbreaks to inform monitoring, coordination and resource deployment. However, genomics studies have remained largely retrospective. While they contain high-resolution views of pathogen diversification and evolution in the context of selection, they are often not aligned with designing interventions. This is a missed opportunity because pathogen diversification is at the core of the most pressing infectious public health challenges, and interventions need to take the mechanisms of virulence and understanding of pathogen diversification into account. In this Perspective, we assess these converging fields, discuss current challenges facing both surveillance specialists and modellers who want to harness genomic data, and propose next steps for integrating longitudinally sampled genomic data with statistical learning and interpretable modelling to make reliable predictions into the future.

Epidemiological dynamics of bacteriocin competition and antibiotic resistance

Bacteriocins, toxic peptides involved in the competition between bacterial strains, are extremely diverse. Previous work on bacteriocin dynamics has highlighted the role of non-transitive 'rock-paper-scissors' competition in maintaining the coexistence of different bacteriocin profiles. The focus to date has primarily been on bacteriocin interactions at the within-host scale (i.e. within a single bacterial population). Yet in species such as Streptococcus pneumoniae, with relatively short periods of colonization and limited within-host diversity, ecological outcomes are also shaped by processes at the epidemiological (between-host) scale. Here, we first investigate bacteriocin dynamics and diversity in epidemiological models. We find that in these models, bacteriocin diversity is more readily maintained than in within-host models, and with more possible combinations of coexisting bacteriocin profiles. Indeed, maintenance of diversity in epidemiological models does not require rock-paper-scissors dynamics; it can also occur through a competition-colonization trade-off. Second, we investigate the link between bacteriocin diversity and diversity at antibiotic resistance loci. Previous work has proposed that bacterial duration of colonization modulates the fitness of antibiotic resistance. Due to their inhibitory effects, bacteriocins are a plausible candidate for playing a role in the duration of colonization episodes. We extend the epidemiological model of bacteriocin dynamics to incorporate an antibiotic resistance locus and demonstrate that bacteriocin diversity can indeed maintain the coexistence of antibiotic-sensitive and -resistant strains.

Does plasmid-based beta-lactam resistance increase E. coli infections: Modelling addition and replacement mechanisms

Infections caused by antibiotic-resistant bacteria have become more prevalent during past decades. Yet, it is unknown whether such infections occur in addition to infections with antibiotic-susceptible bacteria, thereby increasing the incidence of infections, or whether they replace such infections, leaving the total incidence unaffected. Observational longitudinal studies cannot separate both mechanisms. Using plasmid-based beta-lactam resistant E. coli as example we applied mathematical modelling to investigate whether seven biological mechanisms would lead to replacement or addition of infections. We use a mathematical neutral null model of individuals colonized with susceptible and/or resistant E. coli, with two mechanisms implying a fitness cost, i.e., increased clearance and decreased growth of resistant strains, and five mechanisms benefitting resistance, i.e., 1) increased virulence, 2) increased transmission, 3) decreased clearance of resistant strains, 4) increased rate of horizontal plasmid transfer, and 5) increased clearance of susceptible E. coli due to antibiotics. Each mechanism is modelled separately to estimate addition to or replacement of antibiotic-susceptible infections. Fitness costs cause resistant strains to die out if other strain characteristics are maintained equal. Under the assumptions tested, increased virulence is the only mechanism that increases the total number of infections. Other benefits of resistance lead to replacement of susceptible infections without changing the total number of infections. As there is no biological evidence that plasmid-based beta-lactam resistance increases virulence, these findings suggest that the burden of disease is determined by attributable effects of resistance rather than by an increase in the number of infections.

A stochastic model explains the periodicity phenomenon of influenza on network

Influenza is an infectious disease with obvious periodic changes over time. It is of great practical significance to explore the non-environment-related factors that cause this regularity for influenza control and individual protection. In this paper, based on the randomness of population number and the heterogeneity of population contact, we have established a stochastic infectious disease model about influenza based on the degree of the network, and obtained the power spectral density function by using the van Kampen expansion method of the master equation. The relevant parameters are obtained by fitting the influenza data of sentinel hospitals. The results of the numerical analysis show that: (1) for the infected, the infection period of patients who go to the sentinel hospitals is particularly different from the others who do not; (2) for all the infected, there is an obvious nonlinear relationship between their infection period and the visiting rate of the influenza sentinel hospitals, the infection rate and the degree. Among them, only the infection period of patients who do not go to the sentinel hospitals decreased monotonously with the infection rate (increased monotonously with the visiting rate), while the rest had a non-monotonic relationship.

Modeling the effect of vaccination on selection for antibiotic resistance in Streptococcus pneumonia e

Vaccines against bacterial pathogens can protect recipients from becoming infected with potentially antibiotic-resistant pathogens. However, by altering the selective balance between antibiotic-sensitive and antibiotic-resistant bacterial strains, vaccines may also suppress-or spread-antibiotic resistance among unvaccinated individuals. Predicting the outcome of vaccination requires knowing what drives selection for drug-resistant bacterial pathogens and what maintains the circulation of both antibiotic-sensitive and antibiotic-resistant strains of bacteria. To address this question, we used mathematical modeling and data from 2007 on penicillin consumption and penicillin nonsusceptibility in Streptococcus pneumoniae (pneumococcus) invasive isolates from 27 European countries. We show that the frequency of penicillin resistance in S. pneumoniae can be explained by between-host diversity in antibiotic use, heritable diversity in pneumococcal carriage duration, or frequency-dependent selection brought about by within-host competition between antibiotic-resistant and antibiotic-sensitive S. pneumoniae strains. We used our calibrated models to predict the impact of non-serotype-specific pneumococcal vaccination upon the prevalence of S. pneumoniae carriage, incidence of disease, and frequency of S. pneumoniae antibiotic resistance. We found that the relative strength and directionality of competition between drug-resistant and drug-sensitive pneumococcal strains was the most important determinant of whether vaccination would promote, inhibit, or have little effect upon the evolution of antibiotic resistance. Last, we show that country-specific differences in pathogen transmission substantially altered the predicted impact of vaccination, highlighting that policies for managing antibiotic resistance with vaccines must be tailored to a specific pathogen and setting.

Understanding the evolution of multiple drug resistance in structured populations

The evolution of multidrug resistance (MDR) is a pressing public health concern. Yet many aspects, such as the role played by population structure, remain poorly understood. Here, we argue that studying MDR evolution by focusing upon the dynamical equations for linkage disequilibrium (LD) can greatly simplify the calculations, generate more insight, and provide a unified framework for understanding the role of population structure. We demonstrate how a general epidemiological model of MDR evolution can be recast in terms of the LD equations. These equations reveal how the different forces generating and propagating LD operate in a dynamical setting at both the population and metapopulation levels. We then apply these insights to show how the LD perspective: (i) explains equilibrium patterns of MDR, (ii) provides a simple interpretative framework for transient evolutionary dynamics, and (iii) can be used to assess the consequences of different drug prescription strategies for MDR evolution.

Evolutionary mechanisms that determine which bacterial genes are carried on plasmids

The evolutionary pressures that determine the location (chromosomal or plasmid‐borne) of bacterial genes are not fully understood. We investigate these pressures through mathematical modeling in the context of antibiotic resistance, which is often found on plasmids. Our central finding is that gene location is under positive frequency‐dependent selection: the higher the frequency of one form of resistance compared to the other, the higher its relative fitness. This can keep moderately beneficial genes on plasmids, despite occasional plasmid loss. For these genes, positive frequency dependence leads to a priority effect: whichever form is acquired first—through either mutation or horizontal gene transfer—has time to increase in frequency and thus becomes difficult to displace. Higher rates of horizontal transfer of plasmid‐borne than chromosomal genes therefore predict moderately beneficial genes will be found on plasmids. Gene flow between plasmid and chromosome allows chromosomal forms to arise, but positive frequency‐dependent selection prevents these from establishing. Further modeling shows that this effect is particularly pronounced when genes are shared across a large number of species, suggesting that antibiotic resistance genes are often found on plasmids because they are moderately beneficial across many species. We also revisit previous theoretical work—relating to the role of local adaptation in explaining gene location and to plasmid persistence—in light of our findings.

Using ecological coexistence theory to understand antibiotic resistance and microbial competition

Tackling antibiotic resistance necessitates deep understanding of how resource competition within and between species modulates the fitness of resistant microbes. Recent advances in ecological coexistence theory offer a powerful framework to probe the mechanisms regulating intra- and interspecific competition, but the significance of this body of theory to the problem of antibiotic resistance has been largely overlooked. In this Perspective, we draw on emerging ecological theory to illustrate how changes in resource niche overlap can be equally important as changes in competitive ability for understanding costs of resistance and the persistence of resistant pathogens in microbial communities. We then show how different temporal patterns of resource and antibiotic supply, alongside trade-offs in competitive ability at high and low resource concentrations, can have diametrically opposing consequences for the coexistence and exclusion of resistant and susceptible strains. These insights highlight numerous opportunities for innovative experimental and theoretical research into the ecological dimensions of antibiotic resistance.

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

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.

Predicting the future distribution of antibiotic resistance using time series forecasting and geospatial modelling

Background: Increasing antibiotic resistance in a location may be mitigated by changes in treatment policy, or interventions to limit transmission of resistant bacteria. Therefore, accurate forecasting of the distribution of antibiotic resistance could be advantageous. Two previously published studies addressed this, but neither study compared alternative forecasting algorithms or considered spatial patterns of resistance spread. Methods: We analysed data describing the annual prevalence of antibiotic resistance per country in Europe from 2012 – 2016, and the quarterly prevalence of antibiotic resistance per clinical commissioning group in England from 2015 – 2018. We combined these with data on rates of possible covariates of resistance. These data were used to compare the previously published forecasting models, with other commonly used forecasting models, including one geospatial model. Covariates were incorporated into the geospatial model to assess their relationship with antibiotic resistance. Results: For the European data, which was recorded on a coarse spatiotemporal scale, a naïve forecasting model was consistently the most accurate of any of the forecasting models tested. The geospatial model did not improve on this accuracy. However, it did provide some evidence that antibiotic consumption can partially explain the distribution of resistance. The English data were aggregated at a finer scale, and expected-trend-seasonal (ETS) forecasts were the most accurate. The geospatial model did not significantly improve upon the median accuracy of the ETS model, but it appeared to be less sensitive to noise in the data, and provided evidence that rates of antibiotic prescription and bacteraemia are correlated with resistance. Conclusion: Annual, national-level surveillance data appears to be insufficient for fitting accurate antibiotic resistance forecasting models, but there is evidence that data collected at a finer spatiotemporal scale could be used to improve forecast accuracy. Additionally, incorporating antibiotic prescription or consumption data into the model could improve the predictive accuracy.

Factors favouring the evolution of multidrug resistance in bacteria

The evolution of multidrug antibiotic resistance in commensal bacteria is an important public health concern. Commensal bacteria such as Escherichia coli, Streptococcus pneumoniae or Staphylococcus aureus, are also opportunistic pathogens causing a large fraction of the community-acquired and hospital-acquired bacterial infections. Multidrug resistance (MDR) makes these infections harder to treat with antibiotics and may thus cause substantial additional morbidity and mortality. Here, we develop an evolutionary epidemiology model to identify the factors favouring the evolution of MDR in commensal bacteria. The model describes the evolution of antibiotic resistance in a commensal bacterial species evolving in a host population subjected to multiple antibiotic treatments. We combine statistical analysis of a large number of simulations and mathematical analysis to understand the model behaviour. We find that MDR evolves more readily when it is less costly than expected from the combinations of single resistances (positive epistasis). MDR frequently evolves when bacteria are in contact with multiple drugs prescribed in the host population, even if individual hosts are only treated with a single drug at a time. MDR is favoured when the host population is structured in different classes that vary in their rates of antibiotic treatment. However, under most circumstances, recombination between loci involved in resistance does not meaningfully affect the equilibrium frequency of MDR. Together, these results suggest that MDR is a frequent evolutionary outcome in commensal bacteria that encounter the variety of antibiotics prescribed in the host population. A better characterization of the variability in antibiotic use across the host population (e.g. across age classes or geographical location) would help predict which MDR genotypes will most readily evolve.

Population structure across scales facilitates coexistence and spatial heterogeneity of antibiotic-resistant infections

Antibiotic-resistant infections are a growing threat to human health, but basic features of the eco-evolutionary dynamics remain unexplained. Most prominently, there is no clear mechanism for the long-term coexistence of both drug-sensitive and resistant strains at intermediate levels, a ubiquitous pattern seen in surveillance data. Here we show that accounting for structured or spatially-heterogeneous host populations and variability in antibiotic consumption can lead to persistent coexistence over a wide range of treatment coverages, drug efficacies, costs of resistance, and mixing patterns. Moreover, this mechanism can explain other puzzling spatiotemporal features of drug-resistance epidemiology that have received less attention, such as large differences in the prevalence of resistance between geographical regions with similar antibiotic consumption or that neighbor one another. We find that the same amount of antibiotic use can lead to very different levels of resistance depending on how treatment is distributed in a transmission network. We also identify parameter regimes in which population structure alone cannot support coexistence, suggesting the need for other mechanisms to explain the epidemiology of antibiotic resistance. Our analysis identifies key features of host population structure that can be used to assess resistance risk and highlights the need to include spatial or demographic heterogeneity in models to guide resistance management.

A survey of within-host and between-hosts modelling for antibiotic resistance

Antibiotic resistance is a global public health problem which has the attention of many stakeholders including clinicians, the pharmaceutical industry, researchers and policy makers. Despite the existence of many studies, control of resistance transmission has become a rather daunting task as the mechanisms underlying resistance evolution and development are not fully known. Here, we discuss the mechanisms underlying antibiotic resistance development, explore some treatment strategies used in the fight against antibiotic resistance and consider recent findings on collateral susceptibilities amongst antibiotic classes. Mathematical models have proved valuable for unravelling complex mechanisms in biology and such models have been used in the quest of understanding the development and spread of antibiotic resistance. While assessing the importance of such mathematical models, previous systematic reviews were interested in investigating whether these models follow good modelling practice. We focus on theoretical approaches used for resistance modelling considering both within and between host models as well as some pharmacodynamic and pharmakokinetic approaches and further examine the interaction between drugs and host immune response during treatment with antibiotics. Finally, we provide an outlook for future research aimed at modelling approaches for combating antibiotic resistance.

Systematic comparison of coexistence in models of drug-sensitive and drug-resistant pathogen strains

A number of mathematical models have recently been proposed to explain empirical trends of pathogen diversity. In particular, long-term coexistence of both drug-sensitive and drug-resistant variants of a single pathogen is something of a mystery, given that simple models of pathogens competing for the same ecological niche predict competitive exclusion, and more complex models admitting coexistence require assumptions that may not be justified. Coinfection is among the candidate mechanisms to generate coexistence, as it occurs in many pathogens and provides the opportunity for strains to interact directly. Recently, coinfection and competitive release have been described as creating a form of negative frequency-dependent selection that promotes coexistence, and a range of models containing coinfection have been proposed as having generic stable coexistence of multiple strains. This abundance of new models presents the challenge of comparison and interpretation. To this end, we describe a dimensionless quantity that can be used to compare the amount of coexistence generated by different models. We focus on models that include coinfection, although this framework could be generalized to a larger class of structured models.

On the probability of strain invasion in endemic settings: Accounting for individual heterogeneity and control in multi-strain dynamics

•

Endemic infection can insulate host populations from invasion by mutant variants.

•

The timing of control implementation strongly influences its efficacy.

•

Controls that exacerbate host heterogeneity outperform those that curtail it.

•

Differential control can facilitate strain invasion and eventual replacement.

Pathogen evolution is an imminent threat to global health that has warranted, and duly received, considerable attention within the medical, microbiological and modelling communities. Outbreaks of new pathogens are often ignited by the emergence and transmission of mutant variants descended from wild-type strains circulating in the community. In this work we investigate the stochastic dynamics of the emergence of a novel disease strain, introduced into a population in which it must compete with an existing endemic strain. In analogy with past work on single-strain epidemic outbreaks, we apply a branching process approximation to calculate the probability that the new strain becomes established. As expected, a critical determinant of the survival prospects of any invading strain is the magnitude of its reproduction number relative to that of the background endemic strain. Whilst in most circumstances this ratio must exceed unity in order for invasion to be viable, we show that differential control scenarios can lead to less-fit novel strains invading populations hosting a fitter endemic one. This analysis and the accompanying findings will inform our understanding of the mechanisms that have led to past instances of successful strain invasion, and provide valuable lessons for thwarting future drug-resistant strain incursions.

Mathematical modelling for antibiotic resistance control policy: do we know enough?

BACKGROUND

Antibiotics remain the cornerstone of modern medicine. Yet there exists an inherent dilemma in their use: we are able to prevent harm by administering antibiotic treatment as necessary to both humans and animals, but we must be mindful of limiting the spread of resistance and safeguarding the efficacy of antibiotics for current and future generations. Policies that strike the right balance must be informed by a transparent rationale that relies on a robust evidence base.

MAIN TEXT

One way to generate the evidence base needed to inform policies for managing antibiotic resistance is by using mathematical models. These models can distil the key drivers of the dynamics of resistance transmission from complex infection and evolutionary processes, as well as predict likely responses to policy change in silico. Here, we ask whether we know enough about antibiotic resistance for mathematical modelling to robustly and effectively inform policy. We consider in turn the challenges associated with capturing antibiotic resistance evolution using mathematical models, and with translating mathematical modelling evidence into policy.

CONCLUSIONS

We suggest that in spite of promising advances, we lack a complete understanding of key principles. From this we advocate for priority areas of future empirical and theoretical research.

Negative frequency dependent selection on plasmid carriage and low fitness costs maintain extended spectrum β-lactamases in Escherichia coli

Plasmids may maintain antibiotic resistance genes in bacterial populations through conjugation, in the absence of direct selection pressure. However, the costs and benefits of conjugation for plasmid and bacterial fitness are not well understood. Using invasion and competition experiments with plasmid mutants we explicitly tested how conjugation contributes to the maintenance of a plasmid bearing a single extended-spectrum ß-lactamase (ESBL) gene (bla_CTX-M-14). Surprisingly, conjugation had little impact on overall frequencies, although it imposed a substantial fitness cost. Instead, stability resulted from the plasmid conferring fitness benefits when rare. Frequency dependent fitness did not require a functional bla_CTX-M-14 gene, and was independent of culture media. Fitness benefits when rare are associated with the core plasmid backbone but are able to drive up frequencies of antibiotic resistance because fitness burden of the bla_CTX-M-14 gene is very low. Negative frequency dependent fitness can contribute to maintaining a stable frequency of resistance genes in the absence of selection pressure from antimicrobials. In addition, persistent, low cost resistance has broad implications for antimicrobial stewardship.

Generality of Post-Antimicrobial Treatment Persistence of Borrelia burgdorferi Strains N40 and B31 in Genetically Susceptible and Resistant Mouse Strains

A basic feature of infection caused by Borrelia burgdorferi, the etiological agent of Lyme borreliosis, is that persistent infection is the rule in its many hosts. The ability to persist and evade host immune clearance poses a challenge to effective antimicrobial treatment. A link between therapy failure and the presence of persister cells has started to emerge. There is growing experimental evidence that viable but noncultivable spirochetes persist following treatment with several different antimicrobial agents. The current study utilized the mouse model to evaluate if persistence occurs following antimicrobial treatment in disease-susceptible (C3H/HeJ [C3H]) and disease-resistant (C57BL/6 [B6]) mouse strains infected with B. burgdorferi strains N40 and B31 and to confirm the generality of this phenomenon, as well as to assess the persisters' clinical relevance. The status of infection was evaluated at 12 and 18 months after treatment. The results demonstrated that persistent spirochetes remain viable for up to 18 months following treatment, as well as being noncultivable. The phenomenon of persistence in disease-susceptible C3H mice is equally evident in disease-resistant B6 mice and not unique to any particular B. burgdorferi strain. The results also demonstrate that, following antimicrobial treatment, both strains of B. burgdorferi, N40 and B31, lose one or more plasmids. The study demonstrated that noncultivable spirochetes can persist in a host following antimicrobial treatment for a long time but did not demonstrate their clinical relevance in a mouse model of chronic infection. The clinical relevance of persistent spirochetes beyond 18 months following antimicrobial treatment requires further studies in other animal models.

The evolution of antibiotic resistance in a structured host population

The evolution of antibiotic resistance in opportunistic pathogens such as Streptococcus pneumoniae, Escherichia coli or Staphylococcus aureus is a major public health problem, as infection with resistant strains leads to prolonged hospital stay and increased risk of death. Here, we develop a new model of the evolution of antibiotic resistance in a commensal bacterial population adapting to a heterogeneous host population composed of untreated and treated hosts, and structured in different host classes with different antibiotic use. Examples of host classes include age groups and geographic locations. Explicitly modelling the antibiotic treatment reveals that the emergence of a resistant strain is favoured by more frequent but shorter antibiotic courses, and by higher transmission rates. In addition, in a structured host population, localized transmission in host classes promotes both local adaptation of the bacterial population and the global maintenance of coexistence between sensitive and resistant strains. When transmission rates are heterogeneous across host classes, resistant strains evolve more readily in core groups of transmission. These findings have implications for the better management of antibiotic resistance: reducing the rate at which individuals receive antibiotics is more effective to reduce resistance than reducing the duration of treatment. Reducing the rate of treatment in a targeted class of the host population allows greater reduction in resistance, but determining which class to target is difficult in practice.

Whole Slide Imaging for High-Throughput Sensing Antibiotic Resistance at Single-Bacterium Level and Its Application to Rapid Antibiotic Susceptibility Testing

Since conventional culture-based antibiotic susceptibility testing (AST) methods are too time-consuming (typically 24–72 h), rapid AST is urgently needed for preventing the increasing emergence and spread of antibiotic resistant infections. Although several phenotypic antibiotic resistance sensing modalities are able to reduce the AST time to a few hours or less, concerning the biological heterogeneity, their accuracy or limit of detection are limited by low throughput. Here, we present a rapid AST method based on whole slide imaging (WSI)-enabled high-throughput sensing antibiotic resistance at single-bacterium level. The time for determining the minimum inhibitory concentration (MIC) was theoretically shortest, which ensures that the growth of each individual cell present in a large population is inhibited. As a demonstration, our technique was able to sense the growth of at least several thousand bacteria at single-cell level. Reliable MIC of Enterobacter cloacae against gentamicin was obtained within 1 h, while the gold standard broth dilution method required at least 16 h for the same result. In addition, the application of our method prevails over other imaging-based AST approaches in allowing rapid and accurate determination of antibiotic susceptibility for phenotypically heterogeneous samples, in which the number of antibiotic resistant cells was negligible compared to that of the susceptible cells. Hence, our method shows great promise for both rapid AST determination and point-of-care testing of complex clinical bacteria isolates.

On the evolutionary ecology of multidrug resistance in bacteria

Resistance against different antibiotics appears on the same bacterial strains more often than expected by chance, leading to high frequencies of multidrug resistance. There are multiple explanations for this observation, but these tend to be specific to subsets of antibiotics and/or bacterial species, whereas the trend is pervasive. Here, we consider the question in terms of strain ecology: explaining why resistance to different antibiotics is often seen on the same strain requires an understanding of the competition between strains with different resistance profiles. This work builds on models originally proposed to explain another aspect of strain competition: the stable coexistence of antibiotic sensitivity and resistance observed in a number of bacterial species. We first identify a partial structural similarity in these models: either strain or host population structure stratifies the pathogen population into evolutionarily independent sub-populations and introduces variation in the fitness effect of resistance between these sub-populations, thus creating niches for sensitivity and resistance. We then generalise this unified underlying model to multidrug resistance and show that models with this structure predict high levels of association between resistance to different drugs and high multidrug resistance frequencies. We test predictions from this model in six bacterial datasets and find them to be qualitatively consistent with observed trends. The higher than expected frequencies of multidrug resistance are often interpreted as evidence that these strains are out-competing strains with lower resistance multiplicity. Our work provides an alternative explanation that is compatible with long-term stability in resistance frequencies.

Within-host dynamics shape antibiotic resistance in commensal bacteria

The spread of antibiotic resistance, a major threat to human health, is poorly understood. Simple population-level models of bacterial transmission predict that above a certain rate of antibiotic consumption in a population, resistant bacteria should completely eliminate non-resistant strains, while below this threshold they should be unable to persist at all. This prediction stands at odds with empirical evidence showing that resistant and non-resistant strains coexist stably over a wide range of antibiotic consumption rates. Not knowing what drives this long-term coexistence is a barrier to developing evidence-based strategies for managing the spread of resistance. Here, we argue that competition between resistant and sensitive pathogens within individual hosts gives resistant pathogens a relative fitness benefit when they are rare, promoting coexistence between strains at the population level. To test this hypothesis, we embed mechanistically explicit within-host dynamics in a structurally neutral pathogen transmission model. Doing so allows us to reproduce patterns of resistance observed in the opportunistic pathogens Escherichia coli and Streptococcus pneumoniae across European countries and to identify factors that may shape resistance evolution in bacteria by modulating the intensity and outcomes of within-host competition.

Evolutionary epidemiology models to predict the dynamics of antibiotic resistance

The evolution of resistance to antibiotics is a major public health problem and an example of rapid adaptation under natural selection by antibiotics. The dynamics of antibiotic resistance within and between hosts can be understood in the light of mathematical models that describe the epidemiology and evolution of the bacterial population. "Between-host" models describe the spread of resistance in the host community, and in more specific settings such as hospitalized hosts (treated by antibiotics at a high rate), or farm animals. These models make predictions on the best strategies to limit the spread of resistance, such as reducing transmission or adapting the prescription of several antibiotics. Models can be fitted to epidemiological data in the context of intensive care units or hospitals to predict the impact of interventions on resistance. It has proven harder to explain the dynamics of resistance in the community at large, in particular because models often do not reproduce the observed coexistence of drug-sensitive and drug-resistant strains. "Within-host" models describe the evolution of resistance within the treated host. They show that the risk of resistance emergence is maximal at an intermediate antibiotic dose, and some models successfully explain experimental data. New models that include the complex host population structure, the interaction between resistance-determining loci and other loci, or integrating the within- and between-host levels will allow better interpretation of epidemiological and genomic data from common pathogens and better prediction of the evolution of resistance.

The distribution of antibiotic use and its association with antibiotic resistance

Antibiotic use is a primary driver of antibiotic resistance. However, antibiotic use can be distributed in different ways in a population, and the association between the distribution of use and antibiotic resistance has not been explored. Here, we tested the hypothesis that repeated use of antibiotics has a stronger association with population-wide antibiotic resistance than broadly-distributed, low-intensity use. First, we characterized the distribution of outpatient antibiotic use across US states, finding that antibiotic use is uneven and that repeated use of antibiotics makes up a minority of antibiotic use. Second, we compared antibiotic use with resistance for 72 pathogen-antibiotic combinations across states. Finally, having partitioned total use into extensive and intensive margins, we found that intense use had a weaker association with resistance than extensive use. If the use-resistance relationship is causal, these results suggest that reducing total use and selection intensity will require reducing broadly distributed, low-intensity use.

Managing the emergence of pathogen resistance via spatially targeted antimicrobial use

From agriculture to public health to civil engineering, managing antimicrobial resistance presents a considerable challenge. The dynamics underlying resistance evolution reflect inherently spatial processes. Resistant pathogen strains increase in frequency when a strain that emerges in one locale can spread and replace pathogen subpopulations formerly sensitive to the antimicrobial agent. Moreover, the strength of selection for antimicrobial resistance is in part governed by the extent of antimicrobial use. Thus, altering how antimicrobials are used across a landscape can potentially shift the spatial context governing the dynamics of antimicrobial resistance and provide a potent management tool. Here, we model how the efficacy of adjusting antimicrobial use over space to manage antimicrobial resistance is mediated by competition among pathogen strains and the topology of pathogen metapopulations. For several pathogen migration scenarios, we derive critical thresholds for the spatial extent of antimicrobial use below which resistance cannot emerge, and relate these thresholds to (a) the ability to eradicate antimicrobial-sensitive pathogens locally and (b) the strength of the trade-off between resistance ability and competitive performance where antimicrobial use is absent. We find that in metapopulations where patches differ in connectedness, constraining antimicrobial use across space to mitigate resistance evolution only works if the migration of the resistant pathogen is modest; yet, this situation is reversed if the resistant strain has a high colonization rate, with variably connected metapopulations exhibiting less sensitivity to reducing antimicrobial use across space. Furthermore, when pathogens are alternately exposed to sites with and without the antimicrobial, bottlenecking resistant strains through sites without an antimicrobial is only likely to be effective under a strong competition-resistance trade-off. We therefore identify life-history constraints that are likely to suggest which pathogens can most effectively be controlled by a spatially targeted antimicrobial regime. We discuss implications of our results for managing and thinking about antimicrobial resistance evolution in spatially heterogeneous contexts.

Application of dynamic modelling techniques to the problem of antibacterial use and resistance: a scoping review

Selective pressure exerted by the widespread use of antibacterial drugs is accelerating the development of resistant bacterial populations. The purpose of this scoping review was to summarise the range of studies that use dynamic models to analyse the problem of bacterial resistance in relation to antibacterial use in human and animal populations. A comprehensive search of the peer-reviewed literature was performed and non-duplicate articles (n = 1486) were screened in several stages. Charting questions were used to extract information from the articles included in the final subset (n = 81). Most studies (86%) represent the system of interest with an aggregate model; individual-based models are constructed in only seven articles. There are few examples of inter-host models outside of human healthcare (41%) and community settings (38%). Resistance is modelled for a non-specific bacterial organism and/or antibiotic in 40% and 74% of the included articles, respectively. Interventions with implications for antibacterial use were investigated in 67 articles and included changes to total antibiotic consumption, strategies for drug management and shifts in category/class use. The quality of documentation related to model assumptions and uncertainty varies considerably across this subset of articles. There is substantial room to improve the transparency of reporting in the antibacterial resistance modelling literature as is recommended by best practice guidelines.

Host population structure and treatment frequency maintain balancing selection on drug resistance

It is a truism that antimicrobial drugs select for resistance, but explaining pathogen- and population-specific variation in patterns of resistance remains an open problem. Like other common commensals, Streptococcus pneumoniae has demonstrated persistent coexistence of drug-sensitive and drug-resistant strains. Theoretically, this outcome is unlikely. We modelled the dynamics of competing strains of S. pneumoniae to investigate the impact of transmission dynamics and treatment-induced selective pressures on the probability of stable coexistence. We find that the outcome of competition is extremely sensitive to structure in the host population, although coexistence can arise from age-assortative transmission models with age-varying rates of antibiotic use. Moreover, we find that the selective pressure from antibiotics arises not so much from the rate of antibiotic use per se but from the frequency of treatment: frequent antibiotic therapy disproportionately impacts the fitness of sensitive strains. This same phenomenon explains why serotypes with longer durations of carriage tend to be more resistant. These dynamics may apply to other potentially pathogenic, microbial commensals and highlight how population structure, which is often omitted from models, can have a large impact.

Vaccination can drive an increase in frequencies of antibiotic resistance among nonvaccine serotypes of Streptococcus pneumoniae

The bacterial pathogen Streptococcus pneumoniae is a major public health concern, being responsible for more than 1.5 million deaths annually through pneumonia, meningitis, and septicemia. Available vaccines target only a subset of serotypes, so vaccination is often accompanied by a rise in the frequency of nonvaccine serotypes. Epidemiological studies suggest that such a change in serotype frequencies is often coupled with an increase of antibiotic resistance among nonvaccine serotypes. Building on previous multilocus models for bacterial pathogen population structure, we have developed a theoretical framework incorporating variation of serotype and antibiotic resistance to examine how their associations may be affected by vaccination. Using this framework, we find that vaccination can result in a rapid increase in the frequency of preexisting resistant variants of nonvaccine serotypes due to the removal of competition from vaccine serotypes.

An evolutionary model to predict the frequency of antibiotic resistance under seasonal antibiotic use, and an application to Streptococcus pneumoniae

The frequency of resistance to antibiotics in Streptococcus pneumoniae has been stable over recent decades. For example, penicillin non-susceptibility in Europe has fluctuated between 12% and 16% without any major time trend. In spite of long-term stability, resistance fluctuates over short time scales, presumably in part due to seasonal fluctuations in antibiotic prescriptions. Here, we develop a model that describes the evolution of antibiotic resistance under selection by multiple antibiotics prescribed at seasonally changing rates. This model was inspired by, and fitted to, published data on monthly antibiotics prescriptions and frequency of resistance in two communities in Israel over 5 years. Seasonal fluctuations in antibiotic usage translate into small fluctuations of the frequency of resistance around the average value. We describe these dynamics using a perturbation approach that encapsulates all ecological and evolutionary forces into a generic model, whose parameters quantify a force stabilizing the frequency of resistance around the equilibrium and the sensitivity of the population to antibiotic selection. Fitting the model to the data revealed a strong stabilizing force, typically two to five times stronger than direct selection due to antibiotics. The strong stabilizing force explains that resistance fluctuates in phase with usage, as antibiotic selection alone would result in resistance fluctuating behind usage with a lag of three months when antibiotic use is seasonal. While most antibiotics selected for increased resistance, intriguingly, cephalosporins selected for decreased resistance to penicillins and macrolides, an effect consistent in the two communities. One extra monthly prescription of cephalosporins per 1000 children decreased the frequency of penicillin-resistant strains by 1.7%. This model emerges under minimal assumptions, quantifies the forces acting on resistance and explains up to 43% of the temporal variation in resistance.

Use of mathematical modelling to assess the impact of vaccines on antibiotic resistance

Antibiotic resistance is a major global threat to the provision of safe and effective health care. To control antibiotic resistance, vaccines have been proposed as an essential intervention, complementing improvements in diagnostic testing, antibiotic stewardship, and drug pipelines. The decision to introduce or amend vaccination programmes is routinely based on mathematical modelling. However, few mathematical models address the impact of vaccination on antibiotic resistance. We reviewed the literature using PubMed to identify all studies that used an original mathematical model to quantify the impact of a vaccine on antibiotic resistance transmission within a human population. We reviewed the models from the resulting studies in the context of a new framework to elucidate the pathways through which vaccination might impact antibiotic resistance. We identified eight mathematical modelling studies; the state of the literature highlighted important gaps in our understanding. Notably, studies are limited in the range of pathways represented, their geographical scope, and the vaccine-pathogen combinations assessed. Furthermore, to translate model predictions into public health decision making, more work is needed to understand how model structure and parameterisation affects model predictions and how to embed these predictions within economic frameworks.

An agent-based model simulation of influenza interactions at the host level: insight into the influenza-related burden of pneumococcal infections

BACKGROUND

Host-level influenza virus-respiratory pathogen interactions are often reported. Although the exact biological mechanisms involved remain unelucidated, secondary bacterial infections are known to account for a large part of the influenza-associated burden, during seasonal and pandemic outbreaks. Those interactions probably impact the microorganisms' transmission dynamics and the influenza public health toll. Mathematical models have been widely used to examine influenza epidemics and the public health impact of control measures. However, most influenza models overlooked interaction phenomena and ignored other co-circulating pathogens.

METHODS

Herein, we describe a novel agent-based model (ABM) of influenza transmission during interaction with another respiratory pathogen. The interacting microorganism can persist in the population year round (endemic type, e.g. respiratory bacteria) or cause short-term annual outbreaks (epidemic type, e.g. winter respiratory viruses). The agent-based framework enables precise formalization of the pathogens' natural histories and complex within-host phenomena. As a case study, this ABM is applied to the well-known influenza virus-pneumococcus interaction, for which several biological mechanisms have been proposed. Different mechanistic hypotheses of interaction are simulated and the resulting virus-induced pneumococcal infection (PI) burden is assessed.

RESULTS

This ABM generates realistic data for both pathogens in terms of weekly incidences of PI cases, carriage rates, epidemic size and epidemic timing. Notably, distinct interaction hypotheses resulted in different transmission patterns and led to wide variations of the associated PI burden. Interaction strength was also of paramount importance: when influenza increased pneumococcus acquisition, 4-27% of the PI burden during the influenza season was attributable to influenza depending on the interaction strength.

CONCLUSIONS

This open-source ABM provides new opportunities to investigate influenza interactions from a theoretical point of view and could easily be extended to other pathogens. It provides a unique framework to generate in silico data for different scenarios and thereby test mechanistic hypotheses.

Evolution of antibiotic resistance is linked to any genetic mechanism affecting bacterial duration of carriage

Understanding how changes in antibiotic consumption affect the prevalence of antibiotic resistance in bacterial pathogens is important for public health. In a number of bacterial species, including Streptococcus pneumoniae, the prevalence of resistance has remained relatively stable despite prolonged selection pressure from antibiotics. The evolutionary processes allowing the robust coexistence of antibiotic sensitive and resistant strains are not fully understood. While allelic diversity can be maintained at a locus by direct balancing selection, there is no evidence for such selection acting in the case of resistance. In this work, we propose a mechanism for maintaining coexistence at the resistance locus: linkage to a second locus that is under balancing selection and that modulates the fitness effect of resistance. We show that duration of carriage plays such a role, with long duration of carriage increasing the fitness advantage gained from resistance. We therefore predict that resistance will be more common in strains with a long duration of carriage and that mechanisms maintaining diversity in duration of carriage will also maintain diversity in antibiotic resistance. We test these predictions in S. pneumoniae and find that the duration of carriage of a serotype is indeed positively correlated with the prevalence of resistance in that serotype. These findings suggest heterogeneity in duration of carriage is a partial explanation for the coexistence of sensitive and resistant strains and that factors determining bacterial duration of carriage will also affect the prevalence of resistance.

Interim results of an ecological experiment - Conjugate vaccination against the pneumococcus and serotype replacement

Streptococcus pneumoniae has more than 95 serotypes, each of which presumably can cause sepsis, meningitis, pneumonia, and acute otitis media. Pneumococcal conjugate vaccines (PCV) targeted against a limited number of serotypes have nonetheless revealed an impressive impact on each manifestation of pneumococcal disease. At the same time, growing evidence of significant non-vaccine type (NVT) replacement disease following implementation of infant PCV programs has raised questions about the long-term viability of PCV immunization strategies and how to optimize PCV formulations. We discuss here theoretical and practical considerations regarding serotype replacement, and provide a snapshot of the most important NVT types seen to date after implementation of the 2 higher-valent PCVs.

Harnessing Intra-Host Strain Competition to Limit Antibiotic Resistance: Mathematical Model Results

Antibiotic overuse has promoted the spread of antibiotic resistance. To compound the issue, treating individuals dually infected with antibiotic-resistant and antibiotic-vulnerable strains can make their infections completely resistant through competitive release. We formulate mathematical models of transmission dynamics accounting for dual infections and extensions accounting for lag times between infection and treatment or between cure and ending treatment. Analysis using the Next-Generation Matrix reveals how competition within hosts and the costs of resistance determine whether vulnerable and resistant strains persist, coexist, or drive each other to extinction. Invasion analysis predicts that treatment of dually infected cases will promote resistance. By varying antibiotic strength, the models suggest that physicians have two ways to achieve a particular resistance target: prescribe relatively weak antibiotics to everyone infected with an antibiotic-vulnerable strain or give more potent prescriptions to only those patients singly infected with the vulnerable strain after ruling out the possibility of them being dually infected with resistance. Through selectivity and moderation in antibiotic prescription, resistance might be limited.

How direct competition shapes coexistence and vaccine effects in multi-strain pathogen systems

We describe an integrated modeling framework for understanding strain coexistence in polymorphic pathogen systems. Previous studies have debated the utility of neutral formulations and focused on cross-immunity between strains as a major stabilizing mechanism. Here we convey that direct competition for colonization mediates stable coexistence only when competitive abilities amongst pathogen clones satisfy certain pairwise asymmetries. We illustrate our ideas with nested SIS models of single and dual colonization, applied to polymorphic pneumococcal bacteria. By fitting the models to cross-sectional prevalence data from Portugal (before and after the introduction of a seven-valent pneumococcal conjugate vaccine), we are able to not only statistically compare neutral and non-neutral epidemiological formulations, but also estimate vaccine efficacy, transmission and competition parameters simultaneously. Our study highlights that the response of polymorphic pathogen populations to interventions holds crucial information about strain interactions, which can be extracted by suitable nested modeling.

How competition governs whether moderate or aggressive treatment minimizes antibiotic resistance

Understanding how our use of antimicrobial drugs shapes future levels of drug resistance is crucial. Recently, there has been debate over whether an aggressive (i.e., high dose) or more moderate (i.e., lower dose) treatment of individuals will most limit the emergence and spread of resistant bacteria. In this study, we demonstrate how one can understand and resolve these apparently contradictory conclusions. We show that a key determinant of which treatment strategy will perform best at the individual level is the extent of effective competition between resistant and sensitive pathogens within a host. We extend our analysis to the community level, exploring the spectrum between strict inter-strain competition and strain independence. From this perspective as well, we find that the magnitude of effective competition between resistant and sensitive strains determines whether an aggressive approach or moderate approach minimizes the burden of resistance in the population.

DOI:http://dx.doi.org/10.7554/eLife.10559.001

Antibiotics are chemical compounds used to treat bacterial infections. The discovery of antibiotics, starting with penicillin in 1929, revolutionized medicine, making it possible to cure or prevent life-threatening infections such as tetanus and pneumonia. However, many bacteria have become resistant to one or more antibiotics and so can no longer be killed by these drugs.

The emergence of antibiotic resistance reflects an evolutionary process that occurs during antibiotic treatment. While the antibiotic will kill most bacteria, some bacteria may naturally have a feature or genetic mutation that allows them to survive in the presence of the antibiotic. These bacteria then reproduce and pass on their resistant traits, eventually leading to the emergence of a new antibiotic-resistant strain of bacteria. Once a resistant strain exists it may be able to spread from one person to another.

There is conflicting evidence about how best to prevent antibiotic-resistant bacteria from evolving and spreading. The results of some experiments suggest that treating bacteria with large doses of antibiotics early in an infection is the most effective way to optimize treatment and minimize the risk of an antibiotic-resistant strain developing. However, other studies suggest that exposing bacteria to high levels of antibiotics more efficiently selects for resistance; in this case a more moderate approach should be used when treating bacterial infections.

Here, Colijn and Cohen present a mathematical model that suggests that the natural competition between the antibiotic-resistant and antibiotic-sensitive strains of bacteria influence which treatment strategy should be taken. Strains were modeled both within individual hosts and spreading in a community of individuals. In the models, aggressive antibiotic treatment is most effective (in that it minimizes antibiotic resistance) when the antibiotic-resistant strain either does not experience strong competition from the non-resistant strains of bacteria or is not fit enough to be a good competitor. However, a more moderate treatment is appropriate when the two strains are competing and the antibiotic-resistant strain is a fit competitor. Competition may mean that moderate treatment is best to avoid resistance at the community level, even in situations when aggressive treatment is likely best for individuals.

Two important future challenges are to better understand the diversity of strains in bacterial infections, and to develop tools to measure to what extent strains are effectively competing with each other.

DOI:http://dx.doi.org/10.7554/eLife.10559.002

The community ecology of pathogens: coinfection, coexistence and community composition

Disease and community ecology share conceptual and theoretical lineages, and there has been a resurgence of interest in strengthening links between these fields. Building on recent syntheses focused on the effects of host community composition on single pathogen systems, we examine pathogen (microparasite) communities using a stochastic metacommunity model as a starting point to bridge community and disease ecology perspectives. Such models incorporate the effects of core community processes, such as ecological drift, selection and dispersal, but have not been extended to incorporate host-pathogen interactions, such as immunosuppression or synergistic mortality, that are central to disease ecology. We use a two-pathogen susceptible-infected (SI) model to fill these gaps in the metacommunity approach; however, SI models can be intractable for examining species-diverse, spatially structured systems. By placing disease into a framework developed for community ecology, our synthesis highlights areas ripe for progress, including a theoretical framework that incorporates host dynamics, spatial structuring and evolutionary processes, as well as the data needed to test the predictions of such a model. Our synthesis points the way for this framework and demonstrates that a deeper understanding of pathogen community dynamics will emerge from approaches working at the interface of disease and community ecology.