Research Article

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

Institut Pasteur, Université Paris Cité, Epidemiology and Modelling of Antibiotic Evasion (EMAE) unit, France
Université Paris-Saclay, Université de Versailles Saint-Quentin-en-Yvelines, Inserm U1018, CESP, Anti-infective evasion and pharmacoepidemiology team, France
Modélisation, épidémiologie et surveillance des risques sanitaires (MESuRS), Conservatoire national des arts et métiers, France
Health Economics Research Centre, Nuffield Department of Health, University of Oxford, United Kingdom
PACRI unit, Institut Pasteur, Conservatoire national des arts et métiers, France
Centre National de Référence des Pneumocoques, Centre Hospitalier Intercommunal de Créteil, France
Institut Mondor de Recherche Biomédicale-Groupe de Recherche Clinique Groupe d’Etude des Maladies Infectieuses Néonatales et Infantiles (IMRB-GRC GEMINI), Université Paris Est, 94000, France
Groupe de Pathologie Infectieuse Pédiatrique (GPIP), 06200, France
Unité Court Séjour, Petits Nourrissons, Service de Néonatologie, Centre Hospitalier, Intercommunal de Créteil, France
Association Clinique et Thérapeutique Infantile du Val-de-Marne (ACTIV), 94000, France
Association Française de Pédiatrie Ambulatoire (AFPA), 45000, France

Mar 7, 2024

Open access
Copyright information

Abstract
Editor's evaluation
Introduction
Results
Discussion
Methods
Appendix 1
Appendix 2
Data availability
References
Article and author information
Metrics

Abstract

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 R₀ 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.

Editor's evaluation

The mathematical modeling approach disentangles the impacts of the COVID-19 pandemic on antibiotic resistance in Streptococcus pneumoniae in the community setting and identifies the three most plausible driving mechanisms responsible for the observed trends in invasive isolates and pneumococcal carriage. The significance of the findings is important for our understanding of the changed epidemiology of invasive pneumococcal infections during the COVID19 pandemic. The strength of the evidence is graded as solid, derived from a well-described mathematical model evaluating multiple scenarios.

https://doi.org/10.7554/eLife.85701.sa0

Introduction

In the early 2020, international responses to the coronavirus disease 2019 (COVID-19) pandemic led to unprecedented worldwide change in population mixing, healthcare-seeking behavior, and infection prevention and control practices. This modified the ecology and epidemiology of many infectious diseases at a global scale. Strong impacts of COVID-19 on infectious disease dynamics have been reported for common viral and bacterial respiratory infections, sexually transmitted pathogens like HIV, vector-borne diseases like dengue, and even non-communicable diseases (Braunstein et al., 2020; Brueggemann et al., 2021; Chen et al., 2022; Palmer et al., 2020). Antimicrobial resistance (AMR), however, remains one of the leading threats to global health. In 2019, estimates showed that AMR in clinically relevant bacteria was associated with 4.95 million deaths, of which 1.27 million were described as directly attributable to resistance (Murray et al., 2022). Impacts of the COVID-19 pandemic on AMR dynamics remain relatively poorly understood.

A joint report from the WHO and European Centre for Disease Prevention and Control (ECDC) has reported 2020 AMR trends across 29 European countries for eight antibiotic-resistant bacterial pathogens of concern, including S. pneumoniae (European Centre for Disease Prevention and Control and World Health Organization, 2022). While the situation varies widely across bacterial species, antimicrobial groups, and regions, most European countries, including France, documented an increase in pneumococcal resistance to both penicillin and macrolides between 2019 and 2020. The resistance rates rose from 12.2% in 2019 to 15.6% in 2020 for penicillin and from 14.5% in 2019 to 16.9% in 2020 for macrolides, as reported in the EU/EEA (European Centre for Disease Prevention and Control and World Health Organization, 2022). However, increased pneumococcal resistance was accompanied by a sharp worldwide decline in invasive pneumococcal disease (IPD) incidence (Brueggemann et al., 2021; Shaw et al., 2023).

Similar declines in bacterial disease during early waves of COVID-19 have been observed in the context of sentinel community-acquired infections in New Zealand (Duffy et al., 2021), IPDs in Taiwan (Chien et al., 2021) and Hong Kong (Teng et al., 2022), and lower respiratory tract infections in China (Chen et al., 2021). Yet, surprisingly, a growing number of studies have reported mostly stable pneumococcal carriage throughout the COVID-19 pandemic containment, including among infants in Belgium (Willen et al., 2021), children in Vietnam (Nation et al., 2023), Serbia (Petrović et al., 2022), France (Rybak et al., 2022), South Africa (Olwagen et al., 2024), and Israel (Dagan et al., 2023), adults in Connecticut (Wyllie et al., 2023), and households in Seattle (Bennett et al., 2023). In contrast, a study conducted in Denmark reported a decrease in pneumococcal carriage among older adults during the COVID-19 lockdown (Tinggaard et al., 2023).

Understanding the cause of these trends is not straightforward, as many responses to the COVID-19 pandemic, such as the implementation of non-pharmaceutical interventions (NPIs), changes in healthcare-seeking behavior, and alterations in antibiotic prescribing, occurred over the period (Knight et al., 2021). To gain a comprehensive understanding of the changes in AMR epidemiology during the COVID-19 pandemic, it is essential to simultaneously consider a range of scales and indicators. These include the rates of incidence of invasive bacterial diseases (IBDs), the proportion of antibiotic-resistant isolates among total invasive bacterial isolates, and the prevalence of asymptomatic bacterial carriage in healthy individuals.

Several mechanisms may underlie the explanation of these contrasting observations. First, NPIs implemented to block SARS-CoV-2 transmission, such as lockdowns and mask mandates, may have led to reduced bacterial transmission. Containment measures also massively reduced circulation of common respiratory viruses, which are known to be associated with IBD (Domenech de Cellès et al., 2019; Smith and Opatowski, 2021). Second, the lockdown was associated with reductions in primary care consultations (Homeniuk and Collins, 2021; Read et al., 2023; Zhang et al., 2021) leading to a global decrease of antibiotic prescriptions (Högberg et al., 2021). In contrast, frequent antibiotic prescribing to COVID-19 outpatients may have exacerbated AMR (Clancy et al., 2020; Knight et al., 2021). Differences in the duration of pneumococcal carriage may have also played a role (Lehtinen et al., 2017). Finally, potential within-host interactions between SARS-CoV-2 and S. pneumoniae could also have an impact on infection risk (Amin-Chowdhury et al., 2021), although strong evidence for such interactions remains limited (Wong et al., 2023).

Mathematical models incorporating the co-transmission of multiple pathogens within the same host population provide a framework for investigating different hypotheses that underlie the observed patterns in antibiotic resistance and incidence of IPD in S. pneumoniae and help to enhance our understanding of the mechanisms involved. Co-circulation models have been used previously to disentangle the public health consequences of interactions between pathogens such as influenza and S. pneumoniae (Arduin et al., 2017; Domenech de Cellès et al., 2019; Shrestha et al., 2013) and could similarly be used to understand impacts of the COVID-19 pandemic on pathogens coinciding with SARS-CoV-2. However, in a systematic PubMed search conducted on 4 December 2023, we identified no epidemiological models describing the simultaneous transmission of SARS-CoV-2 and antibiotic-resistant bacteria specific to the community setting (Appendix 1).

Here, to disentangle how the COVID-19 pandemic has impacted the epidemiological dynamics of antibiotic resistance in S. pneumoniae, we propose a mathematical model that formalizes the transmission of SARS-CoV-2 and both antibiotic-sensitive and -resistant strains of S. pneumoniae in the community setting, and which includes mechanistic impacts of COVID-19 burden on epidemiological parameters. Through simulation, we assess all possible combinations of these mechanisms to evaluate their overall impact on IPD incidence, antibiotic resistance, and the prevalence of pneumococcal carriage. Furthermore, we assess the changes in the incidence of antibiotic-resistant IPD as we vary the basic reproduction number (R₀) of SARS-CoV-2 during the first COVID-19 outbreak in Europe. We also consider assumed within-host pathogen interactions between SARS-CoV-2 and S. pneumoniae.

Results

Antibiotic resistance trends and incidence of invasive pneumococcal disease in 2020

In routine surveillance data reported to the European Antimicrobial Resistance Surveillance Network (EARS-Net), most European countries reported an increase in antibiotic resistance in S. pneumoniae from 2019 to 2020, as indicated by increases in the proportion of invasive isolates with phenotypic resistance to both penicillin and macrolides (Figure 1A). On the contrary, the total number of reported invasive isolates in the EU/EEA decreased by 44.3% from 2019 to 2020 (European Centre for Disease Prevention and Control and World Health Organization, 2022) suggesting a decrease in incidence of IPD (Appendix 2—table 1).

Figure 1 with 1 supplement see all

Download asset Open asset

Antibiotic resistance trends in invasive *Streptococcus pneumoniae* isolates for the years 2019 and 2020.

(A) The proportion of invasive *S. pneumoniae* isolates resistant to penicillin and macrolides (azithromycin/ clarithromycin/ erythromycin) reported to EARS-Net (European Antimicrobial Resistance Surveillance Network) for 24 European countries. Error bars show 95% confidence intervals. (B) The proportion of invasive *S. pneumoniae* isolates resistant to penicillin (MIC >0.064 mg/L) and macrolides (erythromycin) according to age. Error bars show 95% confidence intervals. The total number of invasive pneumococcal isolates reported in France decreased by 45.1% from 2019 to 2020 (from 1119–614). Data are provided by the French National Reference Center for Pneumococci.

Figure 1—source data 1 European trends in antibiotic resistance in Streptococcus pneumoniae invasive isolates (2019-2020).: https://cdn.elifesciences.org/articles/85701/elife-85701-fig1-data1-v1.xlsx
Download elife-85701-fig1-data1-v1.xlsx

Invasive pneumococcal isolate data for France provided by the French National Reference Center for Pneumococci (CNRP) revealed similar trends. In France, the total number of reported invasive pneumococcal isolates decreased by 45.1% from 2019 to 2020 (from 1119–614), while antibiotic resistance in S. pneumoniae isolates to penicillin and macrolides showed an increasing trend from 26.2% in 2019 to 35.5% in 2020 for penicillin, and from 20.9% in 2019 to 23.0% in 2020 for macrolides (Figure 1B). General decreasing trend in antibiotic resistance from 2017 to 2019 in S. pneumoniae was interrupted in 2020 (Figure 1—figure supplement 1). These variations in antibiotic resistance manifested differently across age, with some age groups showing an increase in antibiotic resistance in 2020 compared to 2019, while others showed no significant change (Figure 1B).

Coinfection model of SARS-CoV-2 and Streptococcus pneumoniae

As mentioned above, several mechanisms may underlie the explanation of these contrasting observations (Figure 2A). COVID-19 NPIs may have led to reduced person-to-person bacterial transmission, potentially contributing to reduced rates of IPD incidence. These containment measures also massively reduced circulation of common respiratory viruses and the incidence of influenza-like-illnesses (ILIs). Respiratory viruses are known triggers and risk factors for developing an IBD from otherwise asymptomatic carriage; in that context, their reduction may have led to reduced infection risk (Domenech de Cellès et al., 2019; Smith and Opatowski, 2021). Due to reductions in primary care consultations in 2020, 26 European countries reported an estimated average decrease of 18.3% in overall antibiotic consumption, aligning with the global trend of reduced antibiotic prescriptions compared to 2019 (Högberg et al., 2021). On the other hand, frequent prescribing of azithromycin, a macrolide antibiotic initially hypothesized to be effective in COVID-19 treatment, has raised concerns for pandemic-associated antimicrobial overuse or misuse and may have exacerbated AMR during and following the first wave of the pandemic (Clancy et al., 2020; Knight et al., 2021; Kournoutou and Dinos, 2022; Langford et al., 2021; PRINCIPLE Trial Collaborative Group, 2021; Rusic et al., 2021). There are still uncertainties about pneumococcal ecology and the evolutionary processes that enable the robust coexistence of strains sensitive and resistant to antibiotics. The role of carriage duration, along with the impact of antibiotic consumption, is also not fully understood in this context. Longer carriage duration of antibiotic-resistant pneumococcal strains is a proposed explanation for this coexistence (Lehtinen et al., 2017). If so, antibiotic-resistant pneumococcal strains may have had an advantage during the lockdown period due to smaller clearance rates, ultimately leading to an increase in antibiotic resistance. Finally, among individuals with COVID-19, potential within-host interactions between SARS-CoV-2 and S. pneumoniae could also have had an impact on bacterial colonization and infection dynamics (Amin-Chowdhury et al., 2021).

Figure 2

Download asset Open asset

A modelling framework describing the transmission of SARS-CoV-2 and *Streptococcus pneumoniae* in the community setting, in the context of both general antibiotic prescribing and azithromycin prescribing for COVID-19 infected individuals.

(A) Non-pharmaceutical interventions (NPIs) implemented to control SARS-CoV-2 transmission (lockdown, face mask use, improved hygiene practices, travel restrictions, quarantine, telemedicine, and physical distancing) may also modify transmission of other pathogens, in addition to impacting antibiotic prescribing due to altered inter-individual contact and health-care seeking behavior. (B) SEIR (Susceptible-Exposed-Infected-Recovered) model with antibiotic treatment compartments depicts interaction between SARS-CoV-2 infection and antibiotic prescribing, including both general community prescribing and azithromycin prescribing among individuals infected with SARS-CoV-2. (C) Diagram depicting how pneumococcal colonization and the community antibiotic prescribing are affected by the COVID-19 pandemic impacts. Initiation of antibiotic treatment is assumed independent of bacterial carriage, reflecting widespread bystander selection for commensal bacteria like *S. pneumoniae*. For a complete modeling framework, see section S2 in Supporting Information.

To test mechanistic impacts of responses to the COVID-19 pandemic on pneumococcal epidemiology, we developed a compartmental, deterministic transmission model describing infection with SARS-CoV-2 being introduced on 1 Jan 2020 (Figure 2B) after colonization with S. pneumoniae reached an equilibrium in a large, well-mixed human population (Figure 2C). Two lockdowns were implemented in the model in agreement with the two lockdowns implemented in France in 2020. The model was parameterized to S. pneumoniae and five mechanisms were built in into the model: (1) a population-wide reduction of antibiotic prescriptions in the community by 18% due to the reduced healthcare-seeking behavior, (2) lockdown reducing pneumococcal transmission by 25%, (3) a reduced risk of developing an IPD from asymptomatic carriage due to the absence of common respiratory viruses during the first lockdown (reduced by a factor IPD_risk = 0.2), which continues after the first lockdown, albeit at a diminished level (IPD_risk = 0.4), (4) community azithromycin use in 10% of COVID-19 infected individuals, (5) and a longer carriage duration of antibiotic-resistant pneumococcal strains giving them a fitness advantage over antibiotic-resistant strains (40 vs 30 days).

Exploring the mechanisms and identifying the optimal scenario for explaining reported trends

We conducted assessments on five distinct hypotheses, each characterized by a precise underlying mechanism, and explored these hypotheses in combination within 31 pandemic scenarios, along with two pre-pandemic (baseline) scenarios, which assume no SARS-CoV-2 circulation in the population and allow for the same 30-day carriage duration (pre-pandemic 1) of both antibiotic-sensitive and -resistant strains (d_s=d_R) or a longer, 40-day carriage duration (pre-pandemic 2) of -resistant strains (d_s> d_R) (Table 1).

Table 1

Five mechanisms implemented in 31 pandemic scenarios proposed to explain the reported trends of IPD incidence, antibiotic resistance, and pneumococcal carriage in S. pneumoniae.

Scenarios explore all possible combinations of mechanisms proposed to test hypotheses that can explain the reported trends of annual invasive pneumococcal disease incidence (annual no. of cases per 100,000 inhabitants), antibiotic resistance (% of annual antibiotic-resistant IPD cases among total IPD cases), and % change in the pneumococcal carriage prevalence at the end of the first 60-day lockdown compared the prevalence before the lockdown. Model simulations were initiated assuming the initial 20% antibiotic resistance. Two pre-pandemic scenarios assume no SARS-CoV-2 circulation in the population and allow for the same 30-day carriage duration of both antibiotic-sensitive and -resistant strains (d_S = d_R) or a longer, 40-day carriage duration of -resistant strains (d_R > d_S) . When implemented, these five mechanisms assume 18% reduction in community antibiotic prescribing, a reduced risk of developing an IPD during (0.2) and after the first lockdown (0.4), a 25% reduction in transmission of pneumococcal carriage during the first lockdown, a 10% of azithromycin use among COVID-19 infected individuals, and a longer 40-day carriage duration of -resistant strains. For a full list of parameters see Appendix 2—table 2. Reported trends in European countries showed a decrease in annual IPD incidence by 44.3% on average, an increase in antibiotic resistance, and generally stable asymptomatic pneumococcal carriage in healthy individuals during the first lockdown period. Only scenarios S19 and S29 fulfill all three reported trends during the COVID-19 pandemic in 2020 simultaneously while accounting for the reported reduction in community antibiotic prescribing ( $d_{S}$ = carriage duration of antibiotic-sensitive pneumococcal strains; $d_{R}$ = carriage duration of antibiotic-resistant pneumococcal strains; PENI = penicillin; ERY = erythromycin).


Scenarios	Mechanisms					IPD inc.	AR (%)	Sp. (%)
Scenarios	1	2	3	4	5	IPD inc.	AR (%)	Sp. (%)
Pre-pandemic 1: (d_S = d_R)						10.8	20.0	NA
Pre-pandemic 2: (d_R > d_S)					x	11.3	20.0	NA
Pandemic: S1	x					10.9	19.2	+1.3
S2		x				8.9	20.1	–36.1
S3			x			5.9	20.0	0
S4				x		9.9	23.7	–9.1
S5					x	11.3	20.0	0
S6	x	x				9.1	19.4	–35.2
S7	x		x			6.0	19.4	+1.3
S8	x			x		10.1	22.9	–8.0
S9	x				x	11.5	19.3	+1.3
S10		x	x			5.2	20.0	–36.1
S11		x		x		8.9	20.1	–36.1
S12		x			x	9.4	20.9	–34.3
S13			x	x		5.6	22.5	–9.1
S14			x		x	6.2	20.0	0
S15				x	x	10.4	23.4	–9.1
S16	x	x	x			5.3	19.6	–35.2
S17	x	x		x		8.3	22.4	–41.3
S18	x	x			x	9.6	20.3	–33.5
S19	x		x	x		5.7	22.0	–8.0
S20	x		x		x	6.3	19.5	+1.3
S21	x			x	x	10.6	22.7	–7.9
S22		x	x	x		5.0	22.0	–42.1
S23		x	x		x	5.5	20.6	–34.3
S24		x		x	x	8.7	23.9	–40.2
S25			x	x	x	5.9	22.3	–9.1
S26	x	x	x	x		5.0	21.6	–41.3
S27	x	x	x		x	5.6	20.1	–33.5
S28	x	x		x	x	8.8	23.2	–39.4
S29	x		x	x	x	5.9	21.8	–7.9
S30		x	x	x	x	5.2	22.5	–40.2
S31	x	x	x	x	x	5.3	22.1	–39.4

REPORTED TRENDS:	IPD inc.				AR (%)			Sp. (%)
Pre-pandemic (FR, 2019)	10.5 [10.3–10.7]				26.2 (PENI) and 20.9 (ERY)			NA
Pandemic (FR, 2020)	5.8 [5.7–5.9]				35.5 (PENI) and 23.0 (ERY)			Stable
*Pandemic (EU/EEA, 2020)* *General trends*	*Decrease by 44.3% on avg.*				*Majority of EU countries report an increase*			*Generally stable*

Mechanisms: (1) Reduced community antibiotic prescribing; (2) Lockdown effect on reducing transmission of S. pneumoniae; (3) Reduced risk of developing an IPD; (4) Community azithromycin use in COVID-19 infected individuals; (5) Longer carriage duration of antibiotic-resistant pneumococcal strains.

We assessed how different combinations of mechanisms may impact: (i) a change in the annual IPD incidence as compared to the pre-pandemic (baseline) period, (ii) antibiotic resistance rate in IPDs, defined as the annual number of antibiotic-resistant IPD cases over the total number of IPD cases, and (iii) daily prevalences of antibiotic-resistant and total pneumococcal carriage in a simulated population of 100,000 individuals (see Appendix 2—table 2 for parameter values). To identify scenarios most compatible with the reported trends, results from model simulations were compared to reported data trends from France in 2020 and more broadly to general EU/EEA reported trends that followed similar patterns. Surprisingly only two scenarios were compatible with reported trends. Scenarios S19 and S29 univocally reproduced increased antibiotic resistance in the general population (AR%) accompanied by a reduction in the annual IPD incidence by almost 50% (IPD inc) with generally stable pneumococcal carriage prevalence in healthy individuals during lockdown (Sp.). In contrast, model simulations revealed that a reduction in the community antibiotic consumption alone (–18%) could not explain the reported trends and generally led to a reduction of antibiotic resistance (Table 1, see scenario S1). Assuming a longer duration of antibiotic-resistant pneumococcal carriage alone did not explain either the rise in antibiotic resistance (Table 1, see scenario S5). Hypothesizing that lockdown reduced the transmission of pneumococcal carriage (by 25%) in addition to a reduced community antibiotic prescribing did not seem probable since, in simulations, this yielded a major reduction in pneumococcal carriage during containment measures in all scenarios where this mechanism was implemented. On the other hand, considering an indirect impact of lockdown on pneumococcal carriage where we implemented a reduction factor for the risk of developing and IPD from otherwise asymptomatic carriage due to the absence of viral respiratory infections during (IPD_risk = 0.2) and after lockdown (IPD_risk = 0.4) reproduced the reported reduction in the annual IPD incidence while maintaining a stable prevalence of pneumococcal carriage during lockdown (Table 1, see scenario S3). By itself however, this scenario did not allow to observe an increase in antibiotic resistance.

When we combined reduced antibiotic prescribing and a reduced risk of developing an IPD with community azithromycin use in a proportion of COVID-19 infected individuals, which remains in the body for an additional 15.5 days after the last dose, in a single scenario, this scenario satisfied the observed trends in AMR (Table 1, see scenario S19). A similar outcome was observed in scenario S29 when adding a longer carriage duration of antibiotic resistant strains on top of this, however, in the absence of community azithromycin use in COVID-19 infected (Table 1, see scenario S1, S5, S3, S19, and S20) trends of increasing antibiotic resistance cannot be reproduced. Therefore, our best model scenario for describing the observed trends combined: (1) a reduction in the overall community antibiotic consumption; (2) the assumption that lockdown effectively reduced SARS-CoV-2 transmission including transmission of other respiratory viruses, but not pneumococcal carriage transmission, indirectly reducing the risk of developing an IPD; (3) either identical or longer carriage durations of antibiotic-resistant strains compared to antibiotic-sensitive strains, and (4) the community azithromycin use in a proportion of COVID-19 infected individuals.

Effect of age

Next, we used the pandemic scenario S19 that best explains the reported trends to test the model using different parameter combinations to mimic different subpopulations (children and the elderly) considering that SARS-CoV-2 infection risk, pneumococcal disease risk, disease severity, bacterial carriage prevalence, and antibiotic prescribing are all highly heterogeneous across age groups. Using scenario S19, we initialized the model with lower and higher baseline carriage prevalence (10%, 20%, and 30%) (Cohen et al., 2023; Rose et al., 2021; Rybak et al., 2022; Tinggaard et al., 2023; Wang et al., 2017), we varied durations of pneumococcal carriage (20, 30, and 45 days), pneumococcal invasion rate, and considered reductions of antibiotic consumption at various levels (–13%, –18%, and –39%) consistent with the French data along with a range of community azithromycin use in COVID-19 infected (0–20%). For a full list of parameters see Appendix 2—table 2. Simulations showed that annual IPD incidence decreased between 43% and 51% compared to the pre-pandemic (baseline) scenario for children, the elderly, and the general population (Figure 3, grey bars). Although the overall antibiotic prescribing in the community was reduced (between 13% and 39%), antibiotic resistance is expected to increase (from 20.1% up to 23.6% in the elderly and from 32.8% up to 36.0% in children) compared to the pre-pandemic period in all age groups and in all scenarios where azithromycin was used in COVID-19 infected individuals (Figure 3, red bars). Daily prevalence of total pneumococcal carriage remained relatively stable, exhibiting higher levels of decrease with increased azithromycin use, while the prevalence of antibiotic-resistant pneumococcal carriage is expected to increase since clearance of antibiotic-susceptible strains due to azithromycin use shifts the competitive balance in favor of the existing resistant strains (Figure 3, third panel).

Figure 3

Download asset Open asset

Annual incidence of invasive pneumococcal disease (IPD), antibiotic resistance (AR%), and pneumococcal carriage prevalence for three different subpopulations.

(A) The elderly (≥65 years-old) (B) general population (all ages), and (C) children (<5 years-old). Using pandemic scenario S19, which includes a combination of three different mechanisms: reduced community antibiotic prescribing, a reduced risk of developing an IPD, and community azithromycin use in COVID-19 infected individuals, we ran model simulations for three different subpopulations. For a full list of parameter values see Appendix 2—table 2. Annual IPD incidence (grey bars) decreased between 43% and 51% relative to the pre-pandemic (baseline) period with magnitude of a decrease depending on an age group and the level of azithromycin use in COVID-19 infected individuals. Antibiotic resistance (red bars) increased compared to the pre-pandemic (baseline) period in all age groups whenever azithromycin was used in COVID-19 infected. Black arrows indicate model outcomes that approximate the reported trends in antibiotic resistance in France for different age groups. Daily prevalence of total pneumococcal carriage remained relatively stable (solid-colored lines), exhibiting higher levels of decrease with increased azithromycin use. The prevalence of antibiotic-resistant pneumococcal carriage increased (dashed colored lines) over time in relation to SARS-CoV-2 outbreak (black dashed line) and higher azithromycin use. Highlighted time intervals (days 75–135 and 305–365) represent two lockdown periods.

General trends produced in model simulations using scenario S19 remained unchanged across different age groups. The extent of the impact depended on the combined magnitude of a decrease in the general antibiotic use in the community and a degree of azithromycin use in COVID-19 infected individuals belonging to a particular age group or a subpopulation. In the elderly (≥65 years-old) and the general population, antibiotic resistance is expected to increase due to azithromycin use in COVID-19 infected. Black arrows indicate model outcomes that approximate the reported trends in antibiotic resistance in France for different age groups including general population (Figure 3). Only in instances when there was no azithromycin use in COVID-19 infected individuals, we observed a decrease in antibiotic resistance relative to the pre-pandemic period (e.g. children <5 years-old). When combining the largest decrease in overall antibiotic use with no or minimal azithromycin use in COVID-19 infected individuals, we expect to see the largest decrease or no change in antibiotic resistance relative to the pre-pandemic period.

Effect of SARS-CoV-2 basic reproduction number (R₀) and within-host pathogen interactions on AMR

Considering that model simulations reproduced an absolute increase in antibiotic resistance comparable to that of 2% reported for macrolides in France but did not reproduce the reported larger increase in penicillin resistance, which was more than a 9% rise (35.5% relative increase) in France, we explored additional factors that may have amplified this increase. Using model scenario S19, we show that an association between higher values of SARS-CoV-2 R₀ and a greater percentage of COVID-19 infected individuals taking azithromycin leads to increased cumulative incidence of antibiotic-resistant IPDs and elevated antibiotic resistance (Figure 4A). For example, if pre-lockdown R₀ of SARS-CoV-2 was 3.8 instead of 3.2, model simulations predict an increase of 3.5% (23.5%) in antibiotic resistance from the pre-pandemic levels instead of 2%. As the R₀ value increases, the impact of azithromycin use becomes more pronounced.

Figure 4

Download asset Open asset

The impact of varying SARS-CoV-2 R₀ and percentage of COVID-19 infected individuals taking azithromycin in scenario S19 on antibiotic resistance (%) and the annual incidence of antibiotic-resistant invasive pneumococcal disease (IPD).

Hypothetical within-host interactions contribute to an excess incidence of antibiotic-resistant IPDs. (A) Cumulative incidence of antibiotic-resistant IPDs and antibiotic resistance increase with greater values of SARS-CoV-2 R₀ and higher percentage of the COVID-19 infected individuals taking azithromycin. The reproduction number for SARS-CoV-2 (R₀) in the community corresponds to the most common estimates of R₀ in France and other European countries ranging from R₀=2–4 (Allieta et al., 2022; D’Arienzo and Coniglio, 2020; Di Domenico et al., 2020; Flaxman et al., 2020; Liu et al., 2020; Roux et al., 2020; Salje et al., 2020). (B) Annual excess in cumulative antibiotic-resistant IPD incidence in scenario S19 due to synergistic within-host ecological interactions compared to the same scenario with no within-host interactions and no azithromycin use (1.17 resistant IPD cases/100,000 inhabitants). A rate of disease progression increased by a factor $ψ_{c} = 1$ (no within-host interaction) and $ψ_{c} = 40$ in scenario S19 applied to the general population assuming azithromycin use in 10% of the infected individuals resulted in approximately 0.06 and 0.75 additional cases of antibiotic-resistant disease per 100,000 inhabitants over the course of 1 year, respectively, compared to the scenario S19 assuming no within-host interaction and no azithromycin use (indicated by the black arrow). For more details, see Appendix 2—figure 1.

Assuming within-host interactions where SARS-CoV-2 infection favors progression from pneumococcal colonization to disease $(ψ_{c} > 1)$ , we found that surges in COVID-19 cases accompanied by increasing levels of azithromycin use lead to excess number of cases caused by antibiotic-resistant strains. Indeed, a rate of disease progression increased by a factor $ψ_{c} = 40$ in scenario S19 with 10% of infected using azithromycin applied to the general population results in approximately 0.75 additional cases of antibiotic-resistant disease per 100,000 inhabitants over the course of 1 year compared to 0.06 additional cases if there are no within-host interactions (Figure 4B). This represents 5% rise in resistance from the pre-pandemic levels (25% relative increase).

Discussion

We propose a novel co-circulation model describing the spread of SARS-CoV-2 and antibiotic-resistant bacteria in a community setting to show how human behavioral responses to the COVID-19 pandemic can differentially impact antibiotic resistance. Our model simulations assessed different hypotheses proposed to explain the observed trends of antibiotic resistance, IPD incidence, and pneumococcal carriage. We identified the most plausible mechanisms underlying the observed patterns of resistance and disease incidence, showing how lockdowns indirectly substantially reduce the incidence of IPD, while surges in COVID-19 cases accompanied by antibiotic prescribing in COVID-19 infected individuals increase antibiotic resistance.

Many studies have reported trends on the incidence of community-acquired bacterial infections since the onset of the pandemic (Brueggemann et al., 2021; Shaw et al., 2023). There was a significant reduction in the risk of invasive disease caused by S. pneumoniae (risk ratio 0·47; 95% CI 0·40–0·55; Shaw et al., 2023). Initially, this observation seemed to support the hypothesis that NPIs implemented to control SARS-CoV-2 transmission may have simultaneously reduced the incidence of bacterial infections by preventing bacterial transmission and acquisition (Brueggemann et al., 2021; Kadambari et al., 2022). Indeed, the scenario of lockdown impact on pneumococcal transmission reproduced such trends. However, incorporating a mechanism of reduced risk for developing an IPD due to the absence of circulation of common respiratory viruses led to similar estimates of the relative reduction in IPD incidence as reported in the EU/EEA for 2020 (Brueggemann et al., 2021; European Centre for Disease Prevention and Control and World Health Organization, 2022). This finding, coupled with the outcome of other studies that found a generally stable pneumococcal carriage prevalence in healthy individuals, both children and adults, during COVID-19 containment measures (Nation et al., 2023; Petrović et al., 2022; Rybak et al., 2022; Willen et al., 2021; Wyllie et al., 2023), supports the alternative hypothesis. This explanation accounts for the decreased incidence of IPD, rather than attributing it to reduced pneumococcal transmission, which resulted in a significant reduction in carriage according to the simulations (Smith and Opatowski, 2021). Furthermore, a study in Vietnam found that reductions in IPD associated with NPIs may be due to reductions in overall pneumococcal carriage density rather than carriage prevalence, driven by reductions in capsular pneumococcal carriage density frequently implicated in IPD (Nation et al., 2023). Considering that common respiratory viruses such as influenza increase pneumococcal carriage density, which contributes to transmission and disease, this hypothesis seems plausible (Alpkvist et al., 2015; Diavatopoulos et al., 2010; McCullers et al., 2010; Short et al., 2012; Wolter et al., 2014).

Globally, community antibiotic consumption dropped during the first year of the COVID-19 pandemic compared to the pre-pandemic period. Decreasing temporal trends were observed in England (Hussain et al., 2021), Canada (Mamun et al., 2021), the United States (Buehrle et al., 2021), China (Zhang et al., 2021), South Korea (Ryu et al., 2021), New Zealand (Duffy et al., 2021), and across European countries (Högberg et al., 2021). In France in particular, the number of antibiotic prescriptions decreased by 18.2% in the general population; however, this reduction ranged from 13% to 39% for the oldest and youngest age groups, respectively (Bara et al., 2022). These trends in antibiotic prescribing may largely be explained by reduced incidence of seasonal respiratory tract infections and reduced primary care consultations (Andrews et al., 2022; Homeniuk and Collins, 2021). On the other hand, the advent of telemedicine, pandemic-induced patient stress, and increased antibiotic demand may have partly offset prescription reductions due to decreased consultations and healthcare-seeking behavior (Hsu, 2020; Read et al., 2023). In a global analysis of antimicrobial sales, Khouja et al. found that antibiotic consumption initially increased by approximately 7% in March 2020, prior to subsequent declines through to August 2020 (Khouja et al., 2022). While overall antibiotic prescribing may have decreased, prescription of specific antibiotics has increased, particularly those associated with COVID-19 patient management. Across continents, a rise of 10% in monthly COVID-19 cases exhibited a correlative trend with elevated macrolide sales of 0.8%, 1.3%, and 1.5% in Europe, North America, and Africa, respectively (Nandi et al., 2023).

Community consumption of azithromycin, a macrolide, increased during the first year of the pandemic in multiple countries with significant variation across geographic locations and with greatest prescribing among older patients (Bara et al., 2022; Bednarčuk et al., 2023; Bogdanić et al., 2022; Crisafulli et al., 2021; Parveen et al., 2020; Weill et al., 2021). In an outpatient setting in southern Italy between February 2020 and January 2021, azithromycin represented 42.1% of all drug prescriptions to individuals diagnosed with COVID-19, while all other antibiotics combined represented just 20.9% (Crisafulli et al., 2021). A study in northwest London across two epidemic waves between January and August 2020 found that, among COVID-19 patients prescribed an antibiotic by a general practitioner during the study period, 31.5% received their prescription within 14 days of a positive SARS-CoV-2 test (Zhu et al., 2021). Two large USA-based studies have also described early pandemic antibiotic prescribing among COVID-19 patients. From April 2020 to April 2021, approximately 30% of outpatient COVID-19–related visits among Medicare beneficiaries (≥65 years-old) have resulted in a filled antibiotic prescription, 50.7% of which were for azithromycin (Tsay et al., 2022). For 0-to-5 year-olds and 45-to-64 year-olds, 4% and 16% of outpatient COVID-19–related visits have resulted in a filled antibiotic prescription, respectively (Wittman et al., 2023). In the Alsace region in France, there was a clear peak azithromycin prescribing during the first wave of the COVID-19 (Danion et al., 2023). During the first lockdown in France, community azithromycin consumption increased by 25.9%, with the increase varying from 13.4% to 47.3% depending on the week (Weill et al., 2021), while the overall number of azithromycin prescriptions across France in 2020 increased by 10.1% relative to 2019 (Bara et al., 2022). Azithromycin treatment usually lasts 3–5 days depending on the disease, but the drug stays in the system for about 15.5 days after the last dose due to the long half-life of more than 60 hr (Foulds et al., 1990; Girard et al., 2005). On the other hand, penicillin has an elimination half-life of approximately 1.4 hr and leaves the body in 7.7 hr after the last dose. This suggests that if azithromycin consumption increased during the first year of the pandemic, antibiotic exposure time also increased as a result, although the overall number of antibiotic prescriptions decreased. Moreover, the use of azithromycin has been associated with selection of both macrolide and non-macrolide resistance (Doan et al., 2020). In a study investigating the direct effect of antibiotic exposure on resistance in the oral streptococcal flora of healthy volunteers, use of azithromycin (500 mg once daily for 3 days) significantly increased the proportion of macrolide-resistant streptococci in healthy individuals (Malhotra-Kumar et al., 2007). Resistance peaked at day four in the azithromycin group and this increase remained significantly higher in the azithromycin group than in the placebo group until day 180 (Malhotra-Kumar et al., 2007). A clinical trial of mass azithromycin distributions for treating trachoma in Ethiopia resulted in an increase in resistant S. pneumoniae isolates among children under the age of 10 (Keenan et al., 2018; Keenan et al., 2015).

Our model simulations show that antibiotic resistance increases with surges in SARS-CoV-2 infections when there is a corresponding increase in azithromycin use, but that lockdowns can moderate this increasing trend by effectively limiting transmission of SARS-CoV-2 (Salje et al., 2020). Conversely, surges in azithromycin prescribing during SARS-CoV-2 outbreaks in the absence of effective measures to prevent transmission, as reported in certain regions and pandemic periods, may cause substantial increases in antibiotic resistance. Our model successfully captured the main trends of antibiotic resistance and IPD incidence observed in Europe in 2020 for S. pneumoniae. However, not all European countries reported an increase in antibiotic resistance. This inter-country heterogeneity may not be due only to heterogeneity of antibiotic use as shown in our model but may be attributed to other pandemic factors not directly implemented or assumed in the model scenario, such as different adherence to COVID-19 control measures across countries and different age groups, including impacts on disease surveillance and data reporting during the pandemic. Real-life scenarios are significantly more complicated and involve multiple alterations of many pandemic factors at different points in time and heterogeneity across populations (e.g. antibiotic prescribing increases in some demographic groups and decreases in others, multiple lockdowns, curfews, or telework).

In our model simulations, we used SARS-CoV-2 parameter value R₀=3.2 (Liu et al., 2020) in the absence of population immunity, best reflecting epidemiological dynamics from early in the pandemic. The most common estimates of SARS-CoV-2 R₀ in France and other European countries ranged from R₀=2–4 (Flaxman et al., 2020; Liu et al., 2020). Modeling results suggest that higher SARS-CoV-2 R₀ estimates combined with higher proportion of COVID-19 infected individuals using azithromycin exacerbated impacts of COVID-19 on antibiotic resistance (Figure 4A). However, the overall impacts of COVID-19 on AMR are difficult to predict, likely vary over the short, medium, and long term, and depend on the organism, setting, and subpopulation considered.

SARS-CoV-2 bacterial coinfection has been reported relatively rarely over the course of the pandemic, suggesting that most COVID-19 patients probably do not require antibiotic therapy (Garcia-Vidal et al., 2021; Karami et al., 2021; Langford et al., 2020), although extensive prophylactic antibiotic use may have limited observed co-infection incidence. The inflammatory immune response resulting from COVID-19 likely predisposes patients to subsequent progression to an IBD to some extent (Sender et al., 2021), but antibiotic use may also favor progression to IBD for patients colonized with drug-resistant strains (Baggs et al., 2018). We do not explicitly model the dynamics of interaction since strong evidence for such interactions remains limited (Wong et al., 2023). The results presented in Figure 4B suggest that such within-host interactions could have important consequences for the resistant IPD incidence during COVID-19 waves, especially in the elderly and high-risk groups. The model’s structure allows for easy integration of mechanistic interactions as more information becomes available on this phenomenon.

Our study focused on the general community, but COVID-19 distinctly influenced AMR in hospitals and long-term care facilities. Extensive antibiotic use in COVID-19 patients and disruptions to antibiotic stewardship programs may have increased antibiotic-resistant carriage in these settings. A meta-analysis conducted on studies published until June 2020 found that 68–81% of hospitalized COVID-19 patients and 74–94% in intensive care received antibiotics (Monnet and Harbarth, 2020). The disorganization in hospitals during the COVID-19 pandemic might have reduced antibiotic resistance surveillance, allowing resistant organisms to spread. However, the early implementation of antibiotic stewardship programs in March 2020, patient isolation, and widespread use of personal protective equipment (PPE) have mitigated this increase to some degree (Henig et al., 2021; Monnet and Harbarth, 2020; Seaton et al., 2020; Van Laethem et al., 2021). Models analyzing these impacts in hospitals contribute to understanding COVID-19’s specific role in the antibiotic resistance burden in different settings (Smith et al., 2023).

A limitation of our model is the lack of age structure and contact patterns between age groups, as SARS-CoV-2 infection risk, pneumococcal disease risk, disease severity, bacterial carriage prevalence and antibiotic prescribing are all highly heterogeneous across age groups. While this choice was made to keep the model as simple as possible, we tested the model using different parameter combinations to mimic different subpopulations (children and ≥65 years-old). This included varying durations of pneumococcal carriage, initializing the model with lower and higher baseline carriage prevalence, considering reductions of general antibiotic consumption at various levels, and varying a percentage of COVID-19 infected individuals using azithromycin. Simulations of the different age groups individually interestingly reproduced realistic trends by age.

In conclusion, we introduce the first epidemiological model outlining the impact of the COVID-19 pandemic on the dynamics of AMR in the community. Our work demonstrates the utility of mathematical modeling approach in unraveling the complex effects of the COVID-19 pandemic responses AMR dynamics. While our model was structured and parameterized based upon S. pneumoniae, its adaptability allows for application to various bacteria and epidemiological scenarios in the community (e.g. impacts of SARS-CoV-2-bacteria interactions in the context of seasonal outbreaks of endemic pathogens). Future research would benefit from fitting the model to real-world data for different bacterial species to enhance our understanding of AMR trends.

Methods

Streptococcus pneumoniae surveillance data

Antibiotic resistance trends reported in 2019 and 2020, provided by EARS-Net (European Antimicrobial Resistance Surveillance Network) were acquired from a joint 2022 report on AMR during 2020 by WHO and ECDC (European Centre for Disease Prevention and Control and World Health Organization, 2022). The annual incidence of S. pneumoniae invasive isolates for 2019 and 2020 was measured as the number of invasive isolates from blood or cerebrospinal fluid. The proportion of resistant isolates represents the proportion of isolates with phenotypic resistance to penicillin and macrolides using standardized bacterial culture methods and EUCAST breakpoints. Out of 28 European countries that reported antibiotic resistance data for S. pneumoniae, 24 countries had enough samples to establish 2019–2020 resistance trends for penicillin and macrolides. The resistance data for France, which were subsequently analyzed, were provided by the CNRP (The French National Reference Center for Pneumococci).

Model structure

We developed a pathogen co-circulation model (Appendix 2—figure 2) written using systems of ordinary differential equations (ODEs) (Appendix 2-Equations; code available on GitHub, copy archived at Kovacevic, 2024). The model simultaneously describes potential infection with SARS-CoV-2 and colonization with antibiotic-sensitive and/or -resistant strains of S. pneumoniae in a well-mixed community population. SARS-CoV-2 infection is modeled by a Susceptible-Exposed-Infectious-Recovered (SEIR) process where individuals become exposed to SARS-CoV-2 at rate β_C upon contact with other infected individuals. Infection begins with a non-infectious exposed period lasting 1/εdays and is followed by an infectious period lasting $1 / γ^{C}$ days, eventually leading to recovery and immunization against future re-infection. Waning immunity and competitive multi-strain SARS-CoV-2 dynamics are not considered.

Individuals in S, E, I, and R compartments can be uncolonized with S. pneumoniae (U), colonized with either a drug-sensitive (C^S) or a drug-resistant strain (C^R), or co-colonized with two strains (C^SS, C^RR, C^SR). Colonization with each respective strain is acquired at rates $β_{S}$ and $β_{s} f$ upon contact with other colonized individuals (Appendix 2—table 2). We assume a metabolic cost of resistance, whereby the drug-resistant strain has a reduced intrinsic transmission rate relative to the drug-sensitive strain due to reduced fitness, f. Bacterial carriage is cleared naturally after an average duration of $\frac{1}{γ^{S}} = \frac{1}{γ^{R}} = \frac{1}{γ^{S R}} = \frac{1}{γ^{S S}} = \frac{1}{γ^{R R}}$ days, which we assume to be the same for all types of carriers in our baseline scenario (in the scenarios assuming longer carriage duration of antibiotic-resistant strains, $\frac{1}{γ^{S}} = \frac{1}{γ^{S S}}$ and $\frac{1}{γ^{R}} = \frac{1}{γ^{S R}} = \frac{1}{γ^{R R}}$ ). We further assume that some share of the population is exposed to antibiotics at any given time, independent of bacterial carriage, with individuals initiating antibiotic therapy at rate $τ$ , which lasts for an average duration of $\frac{1}{r}$ days. Another model assumption is that a proportion $p_{a z}$ of those infected with COVID-19 in the community (between 0% and 20% of individuals in an I compartment) receive azithromycin prescription from general practitioner reflecting azithromycin prescriptions in the early pandemic, while the rest of the infectious individuals ( $1 - p_{a z}$ ) are exposed to the baseline antibiotic therapy. We assume baseline treatment duration of 7 days, on average, regardless of the antibiotic prescribed and without any assumed persistence of the antibiotic in the system after the last dose $\frac{1}{r}$ . In case of antibiotic treatment with azithromycin for COVID-19 infected individuals we assume the treatment lasts three days with antibiotics remaining in the system for additional 15.5 days after the last dose for a total of 18.5 days of antibiotic exposure where COVID-19 recovered individuals $R_{a z}$ treated with azithromycin retain azithromycin in their system for an additional 11.5 days $\frac{1}{r_{a z}}$ after COVID-19 recovery. Individuals treated with antibiotics are unable to acquire the sensitive strain. Antibiotics are assumed to clear colonization with sensitive strains at a rate $ω$ while having no direct impact on colonization with resistant strains. This bacterial colonization process results in antibiotic selection for resistance via competition for limited hosts, facilitates epidemiological coexistence between strains and is adapted from previous models of S. pneumoniae (Colijn et al., 2010; Lipsitch et al., 2009; Mulberry et al., 2020). For a full list of parameter values see Appendix 2—table 2.

Simulation in an early COVID-19 pandemic context

ODEs were integrated numerically using the R package deSolve to simulate and quantify epidemiological dynamics (Soetaert et al., 2010). First, bacterial dynamics were simulated until endemic equilibrium was achieved, under the assumption that S. pneumoniae was at endemic equilibrium upon the emergence of COVID-19. Second, using equilibrium states as initial conditions and re-initializing simulation time to t=0, a single SARS-CoV-2 infected individual was introduced into the population and ODEs were again integrated numerically to t=365 days. Parameter values used for simulation were taken from prior studies prioritizing French data and are provided in Appendix 2—table 2.

These simulations were conducted in the context of an ‘early pandemic scenario’ coinciding with the implementation of population-wide NPIs to slow SARS-CoV-2 transmission. This was conceived as the implementation of two 60 day lockdown periods starting on day 75 and on day 305 in response to the simulated surge in COVID-19 cases. Lockdowns were assumed to have three major potential impacts on population behavior and, in turn, the transmission dynamics of SARS-CoV-2 and S. pneumoniae. These impacts were incorporated into simulations by modifying epidemiological parameters in the model coincident with lockdowns. Three such modifications were considered and switched on and off, considering all possible combinations. First, lockdown led to reduced SARS-CoV-2 transmissibility by a factor $θ_{c}$ . Second, lockdown led to a population-wide change in antibiotic initiation rate by a factor a (representing modified healthcare-seeking behavior leading to a reduction in the number of antibiotic prescriptions). Finally, lockdowns changed the pneumococcal disease risk by a factor ${I P D}_{r i s k}$ (representing a reduced risk of developing an IPD due to the absence of other respiratory viruses).

Effect of SARS-CoV-2 basic reproduction number (R₀) on AMR

Impacts of SARS-CoV-2 on antibiotic-resistant IPD incidence may also depend on the characteristics of locally circulating SARS-CoV-2 R₀. To account for potential impacts of SARS-CoV-2 transmissibility and azithromycin use in the community, in simulations we varied (i) values of R₀ (basic reproduction number) for SARS-CoV-2 in France ( $2 \leq R_{0} \leq 4$ ) and (ii) the proportion of the COVID-19 infected individuals using azithromycin at simulation outset (from 0% to 20%).

Effect of within-host interactions on AMR

SARS-CoV-2 infection may impact progression from bacterial colonization to IBD at the within-host level. To incorporate this mechanism in our model, we included a within-host interaction term in scenario S19: the ecological interaction term $(ψ_{c})$ increases the rate of progression to invasive disease among colonized individuals who are also infected with SARS-CoV-2. The equations for calculating daily IPD incidence assuming within-host interactions due to SARS-CoV-2 co-infection with accompanying details can be found in Appendix 2.

Appendix 1

Literature search

We searched PubMed on December 4, 2023, for published mathematical modeling studies describing coinfection with SARS-CoV-2 and antibiotic-resistant bacteria and including epidemiological transmission of both pathogens. We used the following search query:

("model*"[Title/Abstract] OR "models, theoretical"[MeSH Terms] OR "simulat*"[Title/Abstract])

AND

("mathematic*"[MeSH Terms] OR "compartmental"[Title/Abstract] OR "agent-based"[Title/Abstract] OR "individual-based"[Title/Abstract] OR "stochastic*"[Title/Abstract] OR "deterministic"[Title/Abstract] OR "computer-generated"[Title/Abstract] OR "computer simulation"[MeSH Terms] OR "computer"[Title/Abstract] OR "computational"[Title/Abstract] OR "mathematical*"[Title/Abstract] OR "theoretical"[Title/Abstract] OR "prediction"[Title/Abstract] OR "simulat*"[Title/Abstract])

AND

("outbreak*"[Title/Abstract] OR "dynamic*"[Title/Abstract] OR "disseminat*"[Title/Abstract] OR "propagat*"[Title/Abstract] OR "transmi*"[Title/Abstract] OR "cross infection/transmission"[MeSH Terms] OR "spread*"[Title/Abstract] OR "diffus*"[Title/Abstract] OR "circulati*"[Title/Abstract] OR "disease activity"[Title/Abstract] OR "epidemiology"[Title/Abstract] OR "co-infection"[Title/Abstract])

("COVID-19"[Title/Abstract] OR "SARS-CoV-2"[Title/Abstract])

("bacteria*"[Title/Abstract] OR "antibiotic resist*"[Title/Abstract] OR "antibiotic-resist*"[Title/Abstract] OR "antimicrobial resist*"[Title/Abstract] OR "antimicrobial-resist*"[Title/Abstract]),

which returned 153 papers. One published study developed a mathematical model including the transmission of both SARS-CoV-2 and antibiotic-resistant bacteria but is specific to the hospital setting (Smith et al., 2023). Two bacterial transmission models implicitly accounted for COVID-19 by simulating pandemic-associated social distancing interventions/lockdowns, and corresponding reductions in (i) close-proximity contacts, in the context of transmission of antibiotic-resistance genes among human microbiota (Rebelo et al., 2021), and (ii) sexual contacts, in the context of transmission of HIV, gonorrhea and chlamydia (Jenness et al., 2021). A third model of Neisseria meningitidis transmission also evaluated colonization and infection dynamics over the pandemic period (Cascante-Vega et al., 2023). However, none of these models included SARS-CoV-2 infection and transmission, and only the first considered antibiotic resistance. One further transmission modelling study included coinfection with SARS-CoV-2 and influenza-like illness but did not consider bacteria (Bhowmick et al., 2023). Our search also returned some statistical models used to analyze epidemiological trends in bacterial coinfection subsequent to SARS-CoV-2 infection (Zhou et al., 2023) and changes in the incidence of bacterial disease during the pandemic (Shaw et al., 2023), but these did not use transmission modelling approaches to understand or quantify the underlying epidemiological mechanisms driving these trends. Several other modelling studies have described within-host molecular dynamics or airborne fluid particle dynamics, but do not describe pathogen transmission at the host population level. No other relevant articles were identified.

Appendix 2

Model equations:

Compartment notation:

X_{j}^{i}

X=SARS-CoV-2 status (S=susceptible, E=exposed, I=infected, and R=recovered)

i=bacterial carriage (no superscript =no carriage, S=antibiotic-sensitive carriage, R=antibiotic-resistant carriage, SR =dual carriage, SS =dual antibiotic-sensitive carriage, RR =dual antibiotic-resistant carriage)

j=antibiotic exposure (no subscript =no antibiotic exposure, a=exposed to baseline antibiotic treatment, az =exposed to azithromycin treatment)

Total population size:

\begin{array}{ll} N = & S + S^{S} + S^{R} + S^{S R} + S^{S S} + S^{R R} + S_{a} + S_{a}^{S} + S_{a}^{R} + S_{a}^{S R} + S_{a}^{S S} + S_{a}^{R R} + E + E^{S} + E^{R} \\ + E^{S R} + E^{S S} + E^{R R} + E_{a} + E_{a}^{S} + E_{a}^{R} + E_{a}^{S R} + E_{a}^{S S} + E_{a}^{R R} + I + I^{S} + I^{R} + I^{S R} \\ + I^{S S} + I^{R R} + I_{a} + I_{a}^{S} + I_{a}^{R} + I_{a}^{S R} + I_{a}^{S S} + I_{a}^{R R} + I_{a z} + I_{a z}^{S} + I_{a z}^{R} + I_{a z}^{S R} + I_{a z}^{S S} \\ + I_{a z}^{R R} + R + R^{S} + R^{R} + R^{S R} + R^{S S} + R^{R R} + R_{a} + R_{a}^{S} + R_{a}^{R} + R_{a}^{S R} + R_{a}^{S S} + R_{a}^{R R} \\ + I_{a}^{R R} + R_{a z} + R_{a z}^{S} + R_{a z}^{R} + R_{a z}^{S R} + R_{a z}^{S S} + R_{a z}^{R R} \end{array}

Forces of infection:

S. pneumoniae

\begin{array}{ll} λ_{S} = & β_{S} θ_{β} (S^{S} + 2 q p_{s} S^{S R} + 2 q S^{S S} + S_{a}^{S} + 2 q p_{s} S_{a}^{S R} + 2 q S_{a}^{S S} + E^{S} + 2 q p_{s} E^{S R} + 2 q E^{S S} + E_{a}^{S} \\ + 2 q p_{s} E_{a}^{S R} + 2 q E_{a}^{S S} + I^{S} + 2 q p_{s} I^{S R} + 2 q I^{S S} + I_{a}^{S} + 2 q p_{s} I_{a}^{S R} + 2 q I_{a}^{S S} + I_{a z}^{S} \\ + 2 q p_{s} I_{a z}^{S R} + 2 q I_{a z}^{S S} + R^{S} + 2 q p_{s} R^{S R} + 2 q R^{S S} + R_{a}^{S} + 2 q p_{s} R_{a}^{S R} + 2 q R_{a}^{S S} + R_{a z}^{S} \\ + 2 q p_{s} R_{a z}^{S R} + 2 q R_{a z}^{S S}) / N \end{array}

\begin{array}{ll} λ_{R} = & β_{S} f θ_{β} (S^{R} + 2 q (1 - p_{s}) S^{S R} + 2 q S^{R R} + S_{a}^{R} + 2 q (1 - p_{s}) S_{a}^{S R} + 2 q S_{a}^{R R} + E^{R} \\ + 2 q (1 - p_{s}) E^{S R} + 2 q E^{R R} + E_{a}^{R} + 2 q (1 - p_{s}) E_{a}^{S R} + 2 q E_{a}^{R R} + I^{R} + 2 q (1 - p_{s}) I^{S R} \\ + 2 q I^{R R} + I_{a}^{R} + 2 q (1 - p_{s}) I_{a}^{S R} + 2 q I_{a}^{R R} + I_{a z}^{R} + 2 q (1 - p_{s}) I_{a z}^{S R} + 2 q I_{a z}^{R R} + R^{R} \\ + 2 q (1 - p_{s}) R^{S R} + 2 q R^{R R} + R_{a}^{R} + 2 q (1 - p_{s}) R_{a}^{S R} + 2 q R_{a}^{R R} + R_{a z}^{R} \\ + 2 q (1 - p_{s}) R_{a z}^{S R} + 2 q R_{a z}^{R R}) / N \end{array}

\begin{array}{ll} λ_{X S} = & β_{S} θ_{β} (S^{S} + p_{s i n g l e} (2 q p_{s} S^{S R} + 2 q S^{S S}) + S_{a}^{S} + p_{s i n g l e} (2 q p_{s} S_{a}^{S R} + 2 q S_{a}^{S S}) + E^{S} \\ + p_{s i n g l e} (2 q p_{s} E^{S R} + 2 q E^{S S}) + E_{a}^{S} + p_{s i n g l e} (2 q p_{s} E_{a}^{S R} + 2 q E_{a}^{S S}) + I^{S} \\ + p_{s i n g l e} (2 q p_{s} I^{S R} + 2 q I^{S S}) + I_{a}^{S} + p_{s i n g l e} (2 q p_{s} I_{a}^{S R} + 2 q I_{a}^{S S}) + I_{a z}^{S} \\ + p_{s i n g l e} (2 q p_{s} I_{a z}^{S R} + 2 q I_{a z}^{S S}) + R^{S} + p_{s i n g l e} (2 q p_{s} R^{S R} + 2 q R^{S S}) + R_{a}^{S} \\ + p_{s i n g l e} (2 q p_{s} R_{a}^{S R} + 2 q R_{a}^{S S}) + R_{a z}^{S} + p_{s i n g l e} (2 q p_{s} R_{a z}^{S R} + 2 q R_{a z}^{S S})) / N \end{array}

\begin{array}{ll} λ_{X R} = & β_{S} f θ_{β} (S^{R} + p_{s i n g l e} (2 q (1 - p_{s}) S^{S R} + 2 q S^{R R}) + S_{a}^{R} \\ + p_{s i n g l e} (2 q (1 - p_{s}) S_{a}^{S R} + 2 q S_{a}^{R R}) + E^{R} + p_{s i n g l e} (2 q (1 - p_{s}) E^{S R} + 2 q E^{R R}) \\ + E_{a}^{R} + p_{s i n g l e} (2 q (1 - p_{s}) E_{a}^{S R} + 2 q E_{a}^{R R}) + I^{R} \\ + p_{s i n g l e} (2 q (1 - p_{s}) I^{S R} + 2 q I^{R R}) + I_{a}^{R} + p_{s i n g l e} (2 q (1 - p_{s}) I_{a}^{S R} + 2 q I_{a}^{R R}) \\ + I_{a z}^{R} + p_{s i n g l e} (2 q (1 - p_{s}) I_{a z}^{S R} + 2 q I_{a z}^{R R}) + R^{R} \\ + p_{s i n g l e} (2 q (1 - p_{s}) R^{S R} + 2 q R^{R R}) + R_{a}^{R} + p_{s i n g l e} (2 q (1 - p_{s}) R_{a}^{S R} + 2 q R_{a}^{R R}) \\ + R_{a z}^{R} + p_{s i n g l e} (2 q (1 - p_{s}) R_{a z}^{S R} + 2 q R_{a z}^{R R})) / N \end{array}

\begin{array}{cc} λ_{X S R} = ((1 - p_{s i n g l e}) (β_{S} θ_{β} p_{s} + β_{S} θ_{β} f (1 - p_{s})) 2 q (S^{S R} + S_{a}^{S R} + E^{S R} + E_{a}^{S R} + I^{S R} + I_{a}^{S R} + I_{a z}^{S R} + R^{S R} + R_{a}^{S R} + R_{a z}^{S R})) / N \end{array}

\begin{array}{cc} λ_{X S S} = ((1 - p_{s i n g l e}) β_{S} θ_{β} 2 q (S^{S S} + S_{a}^{S S} + E^{S S} + E_{a}^{S S} + I^{S S} + I_{a}^{S S} \\ + I_{a z}^{S S} + R^{S S} + R_{a}^{S S} + R_{a z}^{S S})) / N \end{array}

\begin{array}{cc} λ_{X R R} = ((1 - p_{s i n g l e}) β_{S} θ_{β} f 2 q (S^{R R} + S_{a}^{R R} + E^{R R} + E_{a}^{R R} + I^{R R} + I_{a}^{R R} \\ + I_{a z}^{R R} + R^{R R} + R_{a}^{R R} + R_{a z}^{R R})) / N \end{array}

SARS-CoV-2:

\begin{array}{l} _{C} = β_{C} θ_{c} ((I + I^{S} + I^{R} + I^{S R} + I^{S S} + I^{R R} + I_{a} + I_{a}^{S} + I_{a}^{R} + I_{a}^{S R} + I_{a}^{S S} + I_{a}^{R R} \\ + I_{a z} + I_{a z}^{S} + I_{a z}^{R} + I_{a z}^{S R} + I_{a z}^{S S} + I_{a z}^{R R})) / N \end{array}

Clearance rates of pneumococcal carriage:

γ^{S} = γ^{R} = γ^{S S} = γ^{R R} = γ^{S R}

The equations for SARS-CoV-2 susceptible individuals (no antibiotic treatment):

\begin{array}{ll} \frac{d S}{d t} = & - (λ_{X S} + λ_{X R} + λ_{X S R} + λ_{X S S} + λ_{X R R}) S + γ^{S} (S^{S} + S^{S S}) + γ^{R} (S^{R} + S^{S R} + S^{R R} + S_{a}^{R} + S_{a}^{R R}) \\ + (γ^{S} + ω) (S_{a}^{S} + S_{a}^{S S}) - τ a S + r S_{a} - λ_{C} S \end{array}

\frac{d S^{S}}{d t} = λ_{X S} S - γ^{S} S^{S} - k λ_{S} S^{S} - k λ_{R} S^{S} - τ a S^{S} - λ_{C} S^{S}

\frac{d S^{R}}{d t} = λ_{X R} S - γ^{R} S^{R} - k λ_{S} S^{R} - k λ_{R} S^{R} - τ a S^{R} + r S_{a}^{R} + (γ^{S} + ω) S_{a}^{S R} - λ_{C} S^{R}

\begin{array}{ll} \frac{d S^{S R}}{d t} = & λ_{X S R} S - γ^{R} S^{S R} + k λ_{R} S^{S} + k λ_{S} S^{R} - k c λ_{S} S^{S R} - k c λ_{R} S^{S R} + 2 k c \\ λ_{R} S^{S S} + 2 k c λ_{S} S^{R R} - τ a S^{S R} - λ_{C} S^{S R} \end{array}

\frac{d S^{S S}}{d t} = λ_{X S S} S - γ^{S} S^{S S} + k λ_{S} S^{S} + k c λ_{S} S^{S R} - 2 k c λ_{R} S^{S S} - τ a S^{S S} - λ_{C} S^{S S}

\frac{d S^{R R}}{d t} = λ_{X R R} S - γ^{R} S^{R R} + k λ_{R} S^{R} + k c λ_{R} S^{S R} - 2 k c λ_{S} S^{R R} - τ a S^{R R} + r S_{a}^{R R} - λ_{C} S^{R R}

The equations for SARS-CoV-2 susceptible individuals (with antibiotic treatment):

\frac{d S_{a}}{d t} = - (λ_{X R} + λ_{X R R}) S_{a} + τ a S - r S_{a} - λ_{C} S_{a}

\frac{d S_{a}^{S}}{d t} = - (γ^{S} + ω) S_{a}^{S} + τ a S^{S} - k λ_{R} S_{a}^{S} - λ_{C} S_{a}^{S}

\frac{d S_{a}^{R}}{d t} = - γ^{R} S_{a}^{R} + λ_{X R} S_{a} + τ a S^{R} - r S_{a}^{R} - λ_{C} S_{a}^{R}

\frac{d S_{a}^{S R}}{d t} = - (γ^{S} + ω) S_{a}^{S R} + k λ_{R} S_{a}^{S} + 2 k c λ_{R} S_{a}^{S S} + τ a S^{S R} - λ_{C} S_{a}^{S R}

\frac{d S_{a}^{S S}}{d t} = - (γ^{S} + ω) S_{a}^{S S} - 2 k c λ_{R} S_{a}^{S S} + τ a S^{S S} - λ_{C} S_{a}^{S S}

\frac{d S_{a}^{R R}}{d t} = λ_{X R R} S_{a} + τ a S^{R R} - r S_{a}^{R R} - γ^{R} S_{a}^{R R} - λ_{C} S_{a}^{R R}

The equations for SARS-CoV-2 exposed individuals (no antibiotic treatment):

\begin{array}{ll} \frac{d E}{d t} = & - (λ_{X S} + λ_{X R} + λ_{X S R} + λ_{X S S} + λ_{X R R}) E + γ^{S} (E^{S} + E^{S S}) + γ^{R} (E^{R} + E^{S R} + E^{R R} + E_{a}^{R} + E_{a}^{R R}) \\ + (γ^{S} + ω) (E_{a}^{S} + E_{a}^{S S}) - τ a E + r E_{a} + λ_{C} S - ε E \end{array}

\frac{d E^{S}}{d t} = λ_{X S} E - γ^{S} E^{S} - k λ_{S} E^{S} - k λ_{R} E^{S} - τ a E^{S} + λ_{C} S^{S} - ε E^{S}

\frac{d E^{R}}{d t} = λ_{X R} E - γ^{R} E^{R} - k λ_{S} E^{R} - k λ_{R} E^{R} - τ a E^{R} + r E_{a}^{R} + λ_{C} S^{R} - ε E^{R} + (γ^{S} + ω) E_{a}^{S R}

\begin{array}{ll} \frac{d E^{S R}}{d t} = & λ_{X S R} E - γ^{R} E^{S R} + k λ_{R} E^{S} + k λ_{S} E^{R} - k c λ_{S} E^{S R} - k c λ_{R} E^{S R} + 2 k c λ_{R} E^{S S} + 2 k c λ_{S} E^{R R} \\ - τ a E^{S R} + λ_{C} S^{S R} - ε E^{S R} \end{array}

\frac{d E^{S S}}{d t} = λ_{X S S} E - γ^{S} E^{S S} + k λ_{S} E^{S} + k c λ_{S} E^{S R} - 2 k c λ_{R} E^{S S} - τ a E^{S S} + λ_{C} S^{S S} - ε E^{S S}

\frac{d E^{R R}}{d t} = λ_{X R R} E - γ^{R} E^{R R} + k λ_{R} E^{R} + k c λ_{R} E^{S R} - 2 k c λ_{S} E^{R R} - τ a E^{R R} + r E_{a}^{R R} + λ_{C} S^{R R} - ε E^{R R}

The equations for SARS-CoV-2 exposed individuals (with antibiotic treatment):

\frac{d E_{a}}{d t} = - (λ_{X R} + λ_{X R R}) E_{a} + τ a E - r E_{a} + λ_{C} S_{a} - ε E_{a}

\frac{d E_{a}^{S}}{d t} = - (γ^{S} + ω) E_{a}^{S} + τ a E^{S} - k λ_{R} E_{a}^{S} + λ_{C} S_{a}^{S} - ε E_{a}^{S}

\frac{d E_{a}^{R}}{d t} = - γ^{R} E_{a}^{R} + λ_{X R} E_{a} + τ a E^{R} - r E_{a}^{R} + λ_{C} E_{a}^{R} - ε E_{a}^{R}

\frac{d E_{a}^{S R}}{d t} = - (γ^{S} + ω) E_{a}^{S R} + k λ_{R} E_{a}^{S} + 2 k c λ_{R} E_{a}^{S S} + τ a E^{S R} + λ_{C} S_{a}^{S R} - ε E_{a}^{S R}

\frac{d E_{a}^{S S}}{d t} = - (γ^{S} + ω) E_{a}^{S S} - 2 k c λ_{R} E_{a}^{S S} + τ a E^{S S} + λ_{C} S_{a}^{S S} - ε E_{a}^{S S}

\frac{d E_{a}^{R R}}{d t} = λ_{X R R} E_{a} + τ a E^{R R} - r E_{a}^{R R} - γ^{R} E_{a}^{R R} + λ_{C} S_{a}^{R R} - ε E_{a}^{R R}

The equations for SARS-CoV-2 infected individuals (no antibiotic treatment):

\begin{array}{ll} \frac{d I}{d t} = & - (λ_{X S} + λ_{X R} + λ_{X S R} + λ_{X S S} + λ_{X R R}) I + γ^{S} (I^{S} + I^{S S}) + γ^{R} (I^{R} + I^{S R} + I^{R R} + I_{a}^{R} + I_{a}^{R R}) \\ + (γ^{S} + ω) (I_{a}^{S} + I_{a}^{S S}) - (1 - p_{a z}) τ a I + r I_{a} + ε E - γ^{C} I - p_{a z} I + γ^{R} (I_{a z}^{R} + I_{a z}^{R R}) \end{array}

\frac{d I^{S}}{d t} = λ_{X S} I - γ^{S} I^{S} - k λ_{S} I^{S} - k λ_{R} I^{S} - (1 - p_{a z}) τ a I^{S} + ε E^{S} - γ^{C} I^{S} - p_{a z} I^{S}

\frac{d I^{R}}{d t} = λ_{X R} I - γ^{R} I^{R} - k λ_{S} I^{R} - k λ_{R} I^{R} + (γ^{S} + ω) I_{a}^{S R} - (1 - p_{a z}) τ a I^{R} + r I_{a}^{R} + ε E^{R} - γ^{C} I^{R} - p_{a z} I^{R}

\begin{array}{ll} \frac{d I^{S R}}{d t} = & λ_{X S R} I - γ^{R} I^{S R} + k λ_{R} I^{S} + k λ_{S} I^{R} - k c λ_{S} I^{S R} - k c λ_{R} I^{S R} + 2 k c λ_{R} I^{S S} + 2 k c λ_{S} I^{R R} \\ - (1 - p_{a z}) τ a I^{S R} + ε E^{S R} - γ^{C} I^{S R} - p_{a z} I^{S R} \end{array}

\frac{d I^{S S}}{d t} = λ_{X S S} I - γ^{S} I^{S S} + k λ_{S} I^{S} + k c λ_{S} I^{S R} - 2 k c λ_{R} I^{S S} - (1 - p_{a z}) τ a I^{S S} + ε E^{S S} - γ^{C} I^{S S} - p_{a z} I^{S S}

\frac{d I^{R R}}{d t} = λ_{X R R} I - γ^{R} I^{R R} + k λ_{R} I^{R} + k c λ_{R} I^{S R} - 2 k c λ_{S} I^{R R} - (1 - p_{a z}) τ a I^{R R} + r I_{a}^{R R} + ε E^{R R} - γ^{C} I^{R R} - p_{a z} I^{R R}

The equations for SARS-CoV-2 infected individuals (with antibiotic treatment):

\frac{d I_{a}}{d t} = - (λ_{X R} + λ_{X R R}) I_{a} + (1 - p_{a z}) τ a I - r I_{a} + ε E_{a} - γ^{C} I_{a}

\frac{d I_{a}^{S}}{d t} = - (γ^{S} + ω) I_{a}^{S} + (1 - p_{a z}) τ a I^{S} - k λ_{R} I_{a}^{S} + ε E_{a}^{S} - γ^{C} I_{a}^{S}

\frac{d I_{a}^{R}}{d t} = - γ^{R} I_{a}^{R} + λ_{X R} I_{a} + (1 - p_{a z}) τ a I^{R} - r I_{a}^{R} + ε E_{a}^{R} - γ^{C} I_{a}^{R}

\frac{d I_{a}^{S R}}{d t} = - (γ^{S} + ω) I_{a}^{S R} + k λ_{R} I_{a}^{S} + 2 k c λ_{R} I_{a}^{S S} + (1 - p_{a z}) τ a I^{S R} + ε E_{a}^{S R} - γ^{C} I_{a}^{S R}

\frac{d I_{a}^{S S}}{d t} = - (γ^{S} + ω) I_{a}^{S S} - 2 k c λ_{R} I_{a}^{S S} + (1 - p_{a z}) τ a I^{S S} + ε E_{a}^{S S} - γ^{C} I_{a}^{S S}

\frac{d I_{a}^{R R}}{d t} = λ_{X R R} I_{a} + (1 - p_{a z}) τ a I^{R R} - r I_{a}^{R R} - γ^{R} I_{a}^{R R} + ε E_{a}^{R R} - γ^{C} I_{a}^{R R}

The equations for SARS-CoV-2 infected individuals (with azithromycin treatment):

\frac{d I_{a z}}{d t} = p_{a z} I - γ^{C} I_{a z} - (λ_{X R} + λ_{X R R}) I_{a z} + (γ^{S} + ω) (I_{a z}^{S} + I_{a z}^{S S})

\frac{d I_{a z}^{S}}{d t} = p_{a z} I^{S} - γ^{C} I_{a z}^{S} - (γ^{S} + ω) I_{a z}^{S} - k λ_{R} I_{a z}^{S}

\frac{d I_{a z}^{R}}{d t} = p_{a z} I^{R} - γ^{C} I_{a z}^{R} + λ_{X R} I_{a z} - γ^{R} I_{a z}^{R} + (γ^{S} + ω) I_{a z}^{S R}

\frac{d I_{a z}^{S R}}{d t} = p_{a z} I^{S R} - γ^{C} I_{a z}^{S R} - (γ^{S} + ω) I_{a z}^{S R} + k λ_{R} I_{a z}^{S} + 2 k c λ_{R} I_{a z}^{S S}

\frac{d I_{a z}^{S S}}{d t} = p_{a z} I^{S S} - γ^{C} I_{a z}^{S S} - (γ^{S} + ω) I_{a z}^{S S} - 2 k c λ_{R} I_{a z}^{S S}

\frac{d I_{a z}^{R R}}{d t} = p_{a z} I^{R R} - γ^{C} I_{a z}^{R R} + λ_{X R R} I_{a z} - γ^{R} I_{a z}^{R R}

The equations for SARS-CoV-2 recovered individuals (no antibiotic treatment):

\begin{array}{ll} \frac{d R}{d t} = & - (λ_{X S} + λ_{X R} + λ_{X S R} + λ_{X S S} + λ_{X R R}) R + γ^{S} (R^{S} + R^{S S}) \\ + γ^{R} (R^{R} + R^{S R} + R^{R R} + R_{a}^{R} + R_{a}^{R R}) + (γ^{S} + ω) (R_{a}^{S} + R_{a}^{S S}) - τ a R + r R_{a} \\ + γ^{C} I + r_{a z} R_{a z} + γ^{R} (R_{a z}^{R} + R_{a z}^{R R}) \end{array}

\frac{d R^{S}}{d t} = λ_{X S} R - γ^{S} R^{S} - k λ_{S} R^{S} - k λ_{R} R^{S} - τ a R^{S} + γ^{C} I^{S} + r_{a z} R_{a z}^{S}

\frac{d R^{R}}{d t} = λ_{X R} R - γ^{R} R^{R} - k λ_{S} R^{R} - k λ_{R} R^{R} + (γ^{S} + ω) R_{a}^{S R} - τ a R^{R} + r R_{a}^{R} + γ^{C} I^{R} + r_{a z} R_{a z}^{R}

\begin{array}{ll} \frac{d R^{S R}}{d t} = & λ_{X S R} R - γ^{R} R^{S R} + k λ_{R} R^{S} + k λ_{S} R^{R} - k c λ_{S} R^{S R} - k c λ_{R} R^{S R} + 2 k c λ_{R} R^{S S} + 2 k c λ_{S} R^{R R} \\ - τ a R^{S R} + γ^{C} I^{S R} + r_{a z} R_{a z}^{S R} \end{array}

\frac{d R^{S S}}{d t} = λ_{X S S} R - γ^{S} R^{S S} + k λ_{S} R^{S} + k c λ_{S} R^{S R} - 2 k c λ_{R} R^{S S} - τ a R^{S S} + γ^{C} I^{S S} + r_{a z} R_{a z}^{S S}

\frac{d R^{R R}}{d t} = λ_{X R R} R - γ^{R} R^{R R} + k λ_{R} R^{R} + k c λ_{R} R^{S R} - 2 k c λ_{S} R^{R R} - τ a R^{R R} + r R_{a}^{R R} + γ^{C} I^{R R} + r_{a z} R_{a z}^{R R}

The equations for SARS-CoV-2 recovered individuals (with antibiotic treatment):

\frac{d R_{a}}{d t} = γ^{C} I_{a} + τ a R - r R_{a} - λ_{X R} R_{a} - λ_{X R R} R_{a}

\frac{d R_{a}^{S}}{d t} = γ^{C} I_{a}^{S} + τ a R^{S} - (γ^{S} + ω) R_{a}^{S} - k λ_{R} R_{a}^{S}

\frac{d R_{a}^{R}}{d t} = γ^{C} I_{a}^{R} + τ a R^{R} - r R_{a}^{R} + λ_{X R} R_{a} - γ^{R} R_{a}^{R}

\frac{d R_{a}^{S R}}{d t} = γ^{C} I_{a}^{S R} + τ a R^{S R} - (γ^{S} + ω) R_{a}^{S R} + k λ_{R} R_{a}^{S} + 2 k c λ_{R} R_{a}^{S S}

\frac{d R_{a}^{S S}}{d t} = γ^{C} I_{a}^{S S} + τ a R^{S S} - (γ^{S} + ω) R_{a}^{S S} - 2 k c λ_{R} R_{a}^{S S}

\frac{d R_{a}^{R R}}{d t} = γ^{C} I_{a}^{R R} + τ a R^{R R} - r R_{a}^{R R} + λ_{X R R} R_{a} - γ^{R} R_{a}^{R R}

The equations for SARS-CoV-2 recovered individuals (with azithromycin treatment):

\frac{d R_{a z}}{d t} = γ^{C} I_{a z} - r_{a z} R_{a z} - λ_{X R} R_{a z} - λ_{X R R} R_{a z} + (γ^{S} + ω) (R_{a z}^{S} + R_{a z}^{S S})

\frac{d R_{a z}^{S}}{d t} = γ^{C} I_{a z}^{S} - r_{a z} R_{a z}^{S} - (γ^{S} + ω) R_{a z}^{S} - k λ_{R} R_{a z}^{S}

\frac{d R_{a z}^{R}}{d t} = γ^{C} I_{a z}^{R} - r_{a z} R_{a z}^{R} + λ_{X R} R_{a z} - γ^{R} R_{a z}^{R} + (γ^{S} + ω) R_{a z}^{S R}

\frac{d R_{a z}^{S R}}{d t} = γ^{C} I_{a z}^{S R} - r_{a z} R_{a z}^{S R} - (γ^{S} + ω) R_{a z}^{S R} + k λ_{R} R_{a z}^{S} + 2 k c λ_{R} R_{a z}^{S S}

\frac{d R_{a z}^{S S}}{d t} = γ^{C} I_{a z}^{S S} - r_{a z} R_{a z}^{S S} - (γ^{S} + ω) R_{a z}^{S S} - 2 k c λ_{R} R_{a z}^{S S}

\frac{d R_{a z}^{R R}}{d t} = γ^{C} I_{a z}^{R R} - r_{a z} R_{a z}^{R R} + τ a R^{R R} - r R_{a}^{R R} + λ_{X R R} R_{a} - γ^{R} R_{a}^{R R}

The equations for invasive pneumococcal disease (IPD) incidence caused by sensitive and resistant strains:

\frac{d I P D^{S}}{d t} = p_{i n v} I P D_{r i s k} (S^{S} + S^{S S} + \frac{1}{2} S^{S R} + E^{S} + E^{S S} + \frac{1}{2} E^{S R} + I^{S} + I^{S S} + \frac{1}{2} I^{S R} + R^{S} + R^{S S} + \frac{1}{2} R^{S R})

\begin{array}{cc} \frac{d I P D^{R}}{d t} = p_{i n v} I P D_{r i s k} (S^{R} + S^{R R} + \frac{1}{2} S^{S R} + E^{R} + E^{R R} + \frac{1}{2} E^{S R} + I^{R} + I^{R R} + \frac{1}{2} I^{S R} + R^{R} + R^{R R} \\ + \frac{1}{2} R^{S R} + S_{a}^{R} + S_{a}^{R R} + S_{a}^{S R} + E_{a}^{R} + E_{a}^{R R} + E_{a}^{S R} + I_{a}^{R} + I_{a}^{R R} + I_{a}^{S R} \\ + R_{a}^{R} + R_{a}^{R R} + R_{a}^{S R} + I_{a z}^{R} + I_{a z}^{R R} + I_{a z}^{S R} + R_{a z}^{R} + R_{a z}^{R R} + R_{a z}^{S R}) \end{array}

Equations for the prevalence of pneumococcal carriage in the population:

Prevalence of antibiotic-sensitive pneumococcal carriage is calculated as:

\begin{array}{ll} S_{p r o p} = & (S^{S} + \frac{1}{2} S^{S R} + S^{S S} + S_{a}^{S} + \frac{1}{2} S_{a}^{S R} + S_{a}^{S S} + E^{S} + \\ \frac{1}{2} E^{S R} + E^{S S} + E_{a}^{S} + \frac{1}{2} E_{a}^{S R} + E_{a}^{S S} + I^{S} + \frac{1}{2} I^{S R} + I^{S S} + I_{a}^{S} + \\ \frac{1}{2} I_{a}^{S R} + I_{a}^{S S} + I_{a z}^{S} + \frac{1}{2} I_{a z}^{S R} + I_{a z}^{S S} + R^{S} + \frac{1}{2} R^{S R} + R^{S S} + R_{a}^{S} + \\ \frac{1}{2} R_{a}^{S R} + R_{a}^{S S} + R_{a z}^{S} + \frac{1}{2} R_{a z}^{S R} + R_{a z}^{S S}) / T o t a l c a r r i a g e \end{array}

Prevalence of antibiotic-resistant pneumococcal carriage is calculated as:

\begin{array}{ll} R_{p r o p} = & (S^{R} + \frac{1}{2} S^{S R} + S^{R R} + S_{a}^{R} + \frac{1}{2} S_{a}^{S R} + S_{a}^{R R} + E^{R} + \\ \frac{1}{2} E^{S R} + E^{R R} + E_{a}^{R} + \frac{1}{2} E_{a}^{S R} + E_{a}^{R R} + I^{R} + \frac{1}{2} I^{S R} + I^{R R} + I_{a}^{R} + \\ \frac{1}{2} I_{a}^{S R} + I_{a}^{R R} + I_{a z}^{R} + \frac{1}{2} I_{a z}^{S R} + I_{a z}^{R R} + R^{R} + \frac{1}{2} R^{S R} + R^{R R} + R_{a}^{R} + \\ \frac{1}{2} R_{a}^{S R} + R_{a}^{R R} + R_{a z}^{R} + \frac{1}{2} R_{a z}^{S R} + R_{a z}^{R R}) / T o t a l c a r r i a g e \end{array}

Total pneumococcal carriage is calculated as:

T o t a l c a r r i a g e = S_{p r o p} + R_{p r o p}

Prevalence of total pneumococcal carriage is calculated as:

T o t a l_{p r o p} = T o t a l c a r r i a g e / N

Equations for within-host interactions:

Invasive pneumococcal disease (IPD) incidence implementing within-host interactions is only accounted for in the equations for IPD incidence outside of the model equations, and hence has no impact on model dynamics.

Daily IPD incidence, for the antibiotic-resistant strain is calculated as:

\begin{array}{ll} \frac{d I P D^{R}}{d t} = & p_{i n v} I P D_{r i s k} (S^{R} + S^{R R} + \frac{1}{2} S^{S R} + E^{R} + E^{R R} + \frac{1}{2} E^{S R} + \\ ψ_{c} (I^{R} + I^{R R} + \frac{1}{2} I^{S R}) + R^{R} + R^{R R} + \frac{1}{2} R^{S R} + S_{a}^{R} + \\ S_{a}^{R R} + S_{a}^{S R} + E_{a}^{R} + E_{a}^{R R} + E_{a}^{S R} + \\ ψ_{c} (I_{a}^{R} + I_{a}^{R R} + I_{a}^{S R}) + R_{a}^{R} + R_{a}^{R R} + R_{a}^{S R} + \\ ψ_{c} (I_{a z}^{R} + I_{a z}^{R R} + I_{a z}^{S R}) + R_{a z}^{R} + R_{a z}^{R R} + R_{a z}^{S R}) \end{array}

where $p_{i n v}$ represents the bacterial invasion rate or progression of carriage to invasive disease and $ψ_{c}$ represents the co-infection interaction term that increases the rate of progression to invasive disease among colonized individuals who are also infected with SARS-CoV-2.

Daily IPD incidence for the antibiotic-sensitive strain is calculated as:

\begin{array}{ll} \frac{d I P D^{S}}{d t} = & p_{i n v} I P D_{r i s k} (S^{R} + S^{R R} + \frac{1}{2} S^{S R} + E^{R} + E^{R R} + \frac{1}{2} E^{S R} \\ + ψ_{c} (I^{R} + I^{R R} + \frac{1}{2} I^{S R}) + R^{R} + R^{R R} + \frac{1}{2} R^{S R} + S_{a}^{R} \\ + S_{a}^{R R} + S_{a}^{S R} + E_{a}^{R} + E_{a}^{R R} + E_{a}^{S R} \\ + ψ_{c} (I_{a}^{R} + I_{a}^{R R} + I_{a}^{S R}) + R_{a}^{R} + R_{a}^{R R} + R_{a}^{S R} \\ + ψ_{c} (I_{a z}^{R} + I_{a z}^{R R} + I_{a z}^{S R}) + R_{a z}^{R} + R_{a z}^{R R} + R_{a z}^{S R}) \end{array}

Due to large uncertainty in the values of $ψ_{c}$ we considered broad range in our analysis estimated from prior modelling studies (Domenech de Cellès et al., 2019; Opatowski et al., 2013). In the context of ecological interactions between S. pneumoniae and influenza, estimations show that influenza co-infection increases the rate of progression from S. pneumonia colonization to invasive pneumococcal disease between 49 and 146-fold, on average, depending on the age group (Domenech de Cellès et al., 2019; Opatowski et al., 2013). Therefore, we implemented the range $1 < ψ_{c} < 100$ for SARS-CoV-2.

Appendix 2—figure 1

Download asset Open asset

Hypothetical within-host interactions in scenario S19 promote the incidence of antibiotic-resistant invasive pneumococcal disease (IPD).

(A) Impacts of ecological interactions between pathogens. When SARS-CoV-2 infection leads to faster progression from pneumococcal colonization to disease , surges in COVID-19 cases accompanied by increasing levels of azithromycin use lead to substantial increases in the daily incidence of antibiotic-resistant IPD. (B) Annual excess in cumulative IPD incidence due to synergistic within-host ecological interactions. A rate of disease progression increased by a factor (no within-host interaction) and in scenario S19 applied to the general population resulted in approximately 0.06 and 0.75 additional cases of antibiotic-resistant disease per 100,000 inhabitants over the course of one year, respectively, compared to the scenario S19 assuming no within-host interaction and no azithromycin use indicated by the black arrow (1.17 cases/100,000 inhabitants).

Appendix 2—figure 2

Download asset Open asset

Model schematic describing co-circulation of SARS-CoV-2 infection and pneumococcal carriage transmission.

(A) Individuals in each of the SARS-CoV-2 transmission model compartments (S, E, I, R) can be either (B) uncolonized (U) or colonized by the bacteria (C^S, C^R, C^SS, C^RR, C^SR) with or without concomitant antibiotic treatment (subscript _atb denotes compartments exposed to antibiotics – both standard baseline antibiotic exposure and azithromycin exposure). Full list of model parameters can be found in Appendix 2—table 2.

Appendix 2—table 1

Total number of invasive S. pneumoniae isolates from blood or cerebrospinal fluid tested in 2019 and 2020 and reported to EARS-Net (European Antimicrobial Resistance Surveillance Network).

	No. of invasive Streptococcus pneumoniae isolates in the EU/EEA
Country	2019	2020	% decrease in the total number of reported invasive isolates from 2019–2020
Austria	550	301	45.3%
Belgium	1548	858	44.6%
Bulgaria	46	28	39.1%
Croatia	156	55	64.7%
Czechia	387	204	47.3%
Denmark	601	351	41.6%
Estonia	161	80	50.3%
Finland	678	293	56.8%
France	1264	668	47.2%
Germany	2035	1314	35.4%
Hungary	222	124	44.1%
Iceland	44	20	54.5%
Ireland	348	136	60.9%
Italy	1351	685	49.3%
Latvia	79	42	46.8%
Lithuania	120	96	20.0%
Luxembourg	38	24	36.8%
Malta	27	16	40.7%
Netherlands	1552	997	35.8%
Norway	507	243	52.1%
Poland	364	165	54.7%
Portugal	983	588	40.2%
Romania	107	42	60.7%
Slovakia	40	15	62.5%
Slovenia	283	172	39.2%
Spain	1038	611	41.1%
Sweden	1071	551	48.6%
United Kingdom	3468	1412	59.3%

Appendix 2—table 2

Summary of the model parameters.

Symbol	Interpretation	Value(s)	References
SARS-CoV-2 infection parameters:
$β_{C}$	transmission rate of SARS-CoV-2	0.46 days^–1	Liu et al., 2020
$γ^{c}$	SARS-CoV-2 recovery rate (mild cases)	1/7 days^–1	Lauer et al., 2020; Rhee et al., 2020
$ϵ$	SARS-CoV-2 incubation rate (latent period from exposed to infectious state)	1/5 days^–1	Elias et al., 2021
$θ_{c}$	Relative risk of SARS-CoV-2 transmission due to lockdown implementation	0.23	Salje et al., 2020
Pneumococcal colonization and invasion parameters:
$β_{s}$	transmission rate of antibiotic-sensitive strain	0.056 days^–10.046 days^–10.034 days^–1	Davies et al., 2019; Olesen et al., 2020
$f$	fitness of antibiotic-resistant strain (assuming there is a fitness cost on transmissibility)	0.9652 0.949 0.926	Dagan et al., 2008; Melnyk et al., 2015
$β_{s} f$	transmission rate of antibiotic-resistant strain		calculated
$θ_{β}$	Relative risk of pneumococcal transmission due to lockdown implementation	1 or 0.75	assumed
$\begin{array}{cc} γ^{S} = γ^{R} = γ^{S S} \\ = γ^{R R} = γ^{S R} \end{array}$	rate of natural bacterial clearance (assumed to be the same for antibiotic-sensitive and -resistant strains)	1/20 days^–11/30 days^–11/45 days^–1	Abdullahi et al., 2012; Davies et al., 2019; Ekdahl et al., 1997; Högberg et al., 2021 Melegaro et al., 2004
$q$	relative infectiousness with each strain for dually colonized	0.5	assumedColijn et al., 2010
$c$	fraction of dually colonized returning to single-colonized upon reinfection	0.5	assumed (Colijn et al., 2010)
$r$	probability of acquiring secondary bacterial carriage	0.5	assumed (Colijn et al., 2010)
$p_{s}$	probability of transmitting antibiotic-sensitive strain	0.5	assumed (Colijn et al., 2010)
$p_{s i n g l e}$	probability of a single infection	0.5	assumed (Colijn et al., 2010)
$p_{i n v}$	pneumococcal invasion rate (summer and winter)	[3x10^-6 day^–1,9x10^-6 day^–1] in the elderly and general population, and [1x10^-6 day^–1, 2.5x10^-6 day^–1] in <5 years-old	Domenech de Cellès et al., 2019; Opatowski et al., 2013
$S^{S} + S^{R}$	initial states – initial prevalence of the total pneumococcal carriage (antibiotic-sensitive and -resistant) in different populations	10% 20% 30%	Cohen et al., 2023; Rose et al., 2021; Rybak et al., 2022; Tinggaard et al., 2023; Wang et al., 2017
Antibiotic exposure parameters:
$ω$	rate of antibiotic-induced pneumococcal clearance for sensitive strains(1/time before antibiotic action)	1/3 days^–1	Kuitunen et al., 2023
$r$	rate of return to antibiotic unexposed compartment (1/duration of antibiotic treatment)	1/7 days^–1	Grant and Saux, 2021; Kuitunen et al., 2023
$r_{a z}$	rate of return to antibiotic unexposed compartment (1/the remainder of how long azithromycin stays in the body)	1/11.5 days^–1	calculated (Foulds et al., 1990; Girard et al., 2005)
$τ$	baseline rate of antibiotic exposure in the community (France)	0.0014 average daily ppc (prescriptions per capita)	Bara et al., 2022
$a$	A reduction factor for antibiotic exposure in the community resulting from changes in healthcare-seeking behavior in response to the COVID-19 pandemic	[0.51, 0.77, 0.84]to represent annual 13%, 18%, and 39% decrease observed in France	Bara et al., 2022
$p_{a z}$	A proportion of COVID-19 infected individuals in the community receiving azithromycin	[0–0.20]testing between 0% and 20%	Tsay et al., 2022; Wittman et al., 2023
Pathogenicity (invasive pneumococcal disease risk):
$I P D_{r i s k}$	A reduction factor for the risk of developing an invasive pneumococcal disease (IPD) due to the absence of influenza-like-illnesses (ILIs) after lockdown implementation	1 (pre-lockdown) 0.2 (lockdown) 0.4 (post-lockdown) for an average of 0.5 in 2020	Shaw et al., 2023

Data availability

Source Data 1 contains the numerical data used to generate Figure 1.

References

1. Abdullahi O
2. Karani A
3. Tigoi CC
4. Mugo D
5. Kungu S
6. Wanjiru E
7. Jomo J
8. Musyimi R
9. Lipsitch M
10. Scott JAG
(2012) Rates of acquisition and clearance of pneumococcal serotypes in the nasopharynges of children in Kilifi District, Kenya
The Journal of Infectious Diseases 206:1020–1029.

https://doi.org/10.1093/infdis/jis447
- PubMed
- Google Scholar
(2022) COVID-19 outbreak in Italy: estimation of reproduction numbers over 2 months prior to phase 2
Zeitschrift Fur Gesundheitswissenschaften = Journal of Public Health 30:2719–2727.

https://doi.org/10.1007/s10389-021-01567-1
- PubMed
- Google Scholar
1. Alpkvist H
2. Athlin S
3. Nauclér P
4. Herrmann B
5. Abdeldaim G
6. Slotved HC
7. Hedlund J
8. Strålin K
(2015) Clinical and microbiological factors associated with high nasopharyngeal pneumococcal density in patients with pneumococcal pneumonia
PLOS ONE 10:e0140112.

https://doi.org/10.1371/journal.pone.0140112
- PubMed
- Google Scholar
1. Amin-Chowdhury Z
2. Aiano F
3. Mensah A
4. Sheppard CL
5. Litt D
6. Fry NK
7. Andrews N
8. Ramsay ME
9. Ladhani SN
(2021) impact of the coronavirus disease 2019 (COVID-19) Pandemic on Invasive Pneumococcal Disease and Risk of Pneumococcal Coinfection With Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2): Prospective National Cohort Study, England
Clinical Infectious Diseases 72:e65–e75.

https://doi.org/10.1093/cid/ciaa1728
- PubMed
- Google Scholar
1. Andrews A
2. Bou-Antoun S
3. Guy R
4. Brown CS
5. Hopkins S
6. Gerver S
(2022) Respiratory antibacterial prescribing in primary care and the COVID-19 pandemic in England, winter season 2020-21
The Journal of Antimicrobial Chemotherapy 77:799–802.

https://doi.org/10.1093/jac/dkab443
- PubMed
- Google Scholar
(2017) An agent-based model simulation of influenza interactions at the host level: insight into the influenza-related burden of pneumococcal infections
BMC Infectious Diseases 17:382.

https://doi.org/10.1186/s12879-017-2464-z
- PubMed
- Google Scholar
(2018) Risk of subsequent sepsis within 90 days after a hospital stay by type of antibiotic exposure
Clinical Infectious Diseases 66:1004–1012.

https://doi.org/10.1093/cid/cix947
- PubMed
- Google Scholar
(2022) Outpatient antibiotic prescriptions in france: patients and providers characteristics and impact of the COVID-19 Pandemic
Antibiotics 11:643.

https://doi.org/10.3390/antibiotics11050643
- PubMed
- Google Scholar
(2023) Antibiotic Utilization during COVID-19: are we over-prescribing?
Antibiotics 12:308.

https://doi.org/10.3390/antibiotics12020308
- Google Scholar
1. Bennett JC
2. Emanuels A
3. Heimonen J
4. O’Hanlon J
5. Hughes JP
6. Han PD
7. Chow EJ
8. Ogokeh CE
9. Rolfes MA
10. Lockwood CM
11. Pfau B
12. Uyeki TM
13. Shendure J
14. Hoag S
15. Fay K
16. Lee J
17. Sibley TR
18. Rogers JH
19. Starita LM
20. Englund JA
21. Chu HY
(2023) Streptococcus pneumoniae nasal carriage patterns with and without common respiratory virus detections in households in Seattle, WA, USA before and during the COVID-19 pandemic
Frontiers in Pediatrics 11:1198278.

https://doi.org/10.3389/fped.2023.1198278
- PubMed
- Google Scholar
(2023) Decoding the double trouble: A mathematical modelling of co-infection dynamics of SARS-CoV-2 and influenza-like illness
Bio Systems 224:104827.

https://doi.org/10.1016/j.biosystems.2023.104827
- PubMed
- Google Scholar
(2022) Azithromycin consumption during the COVID-19 pandemic in Croatia, 2020
PLOS ONE 17:e0263437.

https://doi.org/10.1371/journal.pone.0263437
- PubMed
- Google Scholar
(2020) COVID-19 infection among people with HIV in New York City: a population-level analysis of matched surveillance data
SSRN Electronic Journal 01:3675416.

https://doi.org/10.2139/ssrn.3675416
- Google Scholar
1. Brueggemann AB
2. Jansen van Rensburg MJ
3. Shaw D
4. McCarthy ND
5. Jolley KA
6. Maiden MCJ
7. van der Linden MPG
8. Amin-Chowdhury Z
9. Bennett DE
10. Borrow R
11. Brandileone MCC
12. Broughton K
13. Campbell R
14. Cao B
15. Casanova C
16. Choi EH
17. Chu YW
18. Clark SA
19. Claus H
20. Coelho J
21. Corcoran M
22. Cottrell S
23. Cunney RJ
24. Dalby T
25. Davies H
26. de Gouveia L
27. Deghmane AE
28. Demczuk W
29. Desmet S
30. Drew RJ
31. du Plessis M
32. Erlendsdottir H
33. Fry NK
34. Fuursted K
35. Gray SJ
36. Henriques-Normark B
37. Hale T
38. Hilty M
39. Hoffmann S
40. Humphreys H
41. Ip M
42. Jacobsson S
43. Johnston J
44. Kozakova J
45. Kristinsson KG
46. Krizova P
47. Kuch A
48. Ladhani SN
49. Lâm TT
50. Lebedova V
51. Lindholm L
52. Litt DJ
53. Martin I
54. Martiny D
55. Mattheus W
56. McElligott M
57. Meehan M
58. Meiring S
59. Mölling P
60. Morfeldt E
61. Morgan J
62. Mulhall RM
63. Muñoz-Almagro C
64. Murdoch DR
65. Murphy J
66. Musilek M
67. Mzabi A
68. Perez-Argüello A
69. Perrin M
70. Perry M
71. Redin A
72. Roberts R
73. Roberts M
74. Rokney A
75. Ron M
76. Scott KJ
77. Sheppard CL
78. Siira L
79. Skoczyńska A
80. Sloan M
81. Slotved HC
82. Smith AJ
83. Song JY
84. Taha MK
85. Toropainen M
86. Tsang D
87. Vainio A
88. van Sorge NM
89. Varon E
90. Vlach J
91. Vogel U
92. Vohrnova S
93. von Gottberg A
94. Zanella RC
95. Zhou F
(2021) Changes in the incidence of invasive disease due to Streptococcus pneumoniae, Haemophilus influenzae, and Neisseria meningitidis during the COVID-19 pandemic in 26 countries and territories in the Invasive Respiratory Infection Surveillance Initiative: a prospective analysis of surveillance data
The Lancet. Digital Health 3:e360–e370.

https://doi.org/10.1016/S2589-7500(21)00077-7
- PubMed
- Google Scholar
(2021) Trends in outpatient antibiotic prescriptions in the united states during the covid-19 pandemic in 2020
JAMA Network Open 4:e2126114.

https://doi.org/10.1001/jamanetworkopen.2021.26114
- PubMed
- Google Scholar
1. Cascante-Vega J
2. Galanti M
3. Schley K
4. Pei S
5. Shaman J
(2023) Inference of transmission dynamics and retrospective forecast of invasive meningococcal disease
PLOS Computational Biology 19:e1011564.

https://doi.org/10.1371/journal.pcbi.1011564
- PubMed
- Google Scholar
1. Chen C
2. Zhu P
3. Zhang Y
4. Liu B
(2021) Effect of the “Normalized Epidemic Prevention and Control Requirements” on hospital-acquired and community-acquired infections in China
BMC Infectious Diseases 21:1178.

https://doi.org/10.1186/s12879-021-06886-y
- PubMed
- Google Scholar
1. Chen Y
2. Li N
3. Lourenço J
4. Wang L
5. Cazelles B
6. Dong L
7. Li B
8. Liu Y
9. Jit M
10. Bosse NI
11. Abbott S
12. Velayudhan R
13. Wilder-Smith A
14. Tian H
15. Brady OJ
16. Procter SR
17. Wong KL
18. Hellewell J
19. Davies NG
20. Jarvis CI
21. McCarthy CV
22. Medley G
23. Meakin SR
24. Rosello A
25. Finch E
26. Lowe R
27. Pearson CAB
28. Clifford S
29. Quilty BJ
30. Flasche S
31. Gibbs HP
32. Chapman LAC
33. Atkins KE
34. Hodgson D
35. Barnard RC
36. Russell TW
37. Klepac P
38. Jafari Y
39. Eggo RM
40. Mee P
41. Quaife M
42. Endo A
43. Funk S
44. Hué S
45. Kucharski AJ
46. Edmunds WJ
47. O’Reilly K
48. Pung R
49. Villabona-Arenas CJ
50. Gimma A
51. Abbas K
52. Prem K
53. Knight GM
54. Sun FY
55. Waites W
56. Munday JD
57. Koltai M
58. Sandmann FG
59. Tully DC
(2022) Measuring the effects of COVID-19-related disruption on dengue transmission in southeast Asia and Latin America: a statistical modelling study
The Lancet. Infectious Diseases 22:657–667.

https://doi.org/10.1016/S1473-3099(22)00025-1
- PubMed
- Google Scholar
1. Chien Y-C
2. Lee Y-L
3. Liu P-Y
4. Lu M-C
5. Shao P-L
6. Lu P-L
7. Cheng S-H
8. Lin C-Y
9. Wu T-S
10. Yen M-Y
11. Wang L-S
12. Liu C-P
13. Lee W-S
14. Shi Z-Y
15. Chen Y-S
16. Wang F-D
17. Tseng S-H
18. Chen Y-H
19. Sheng W-H
20. Lee C-M
21. Chen Y-H
22. Ko W-C
23. Hsueh P-R
(2021) National surveillance of antimicrobial susceptibilities to dalbavancin, telavancin, tedizolid, eravacycline, omadacycline and other comparator antibiotics and serotype distribution of invasive Streptococcus pneumoniae isolates in adults: results from the Surveillance of Multicenter Antimicrobial Resistance in Taiwan (SMART) programme in 2017–2020
Journal of Global Antimicrobial Resistance 26:308–316.

https://doi.org/10.1016/j.jgar.2021.07.005
- Google Scholar
(2020) PRO: The COVID-19 pandemic will result in increased antimicrobial resistance rates
JAC-Antimicrobial Resistance 2:dlaa049.

https://doi.org/10.1093/jacamr/dlaa049
- PubMed
- Google Scholar
1. Cohen R
2. Bidet P
3. Varon E
4. Béchet S
5. Cohen JF
6. Bonacorsi S
7. Levy C
(2023) Unprecedentedly high rates of Group A Streptococcus nasopharyngeal carriage in infants and toddlers in France, 2022-2023
Infectious Diseases Now 53:104720.

https://doi.org/10.1016/j.idnow.2023.104720
- PubMed
- Google Scholar
1. Colijn C
2. Cohen T
3. Fraser C
4. Hanage W
5. Goldstein E
6. Givon-Lavi N
7. Dagan R
8. Lipsitch M
(2010) What is the mechanism for persistent coexistence of drug-susceptible and drug-resistant strains of Streptococcus pneumoniae?
Journal of the Royal Society, Interface 7:905–919.

https://doi.org/10.1098/rsif.2009.0400
- PubMed
- Google Scholar
1. Crisafulli S
2. Ientile V
3. L’Abbate L
4. Fontana A
5. Linguiti C
6. Manna S
7. Mercaldo M
8. Pagliaro C
9. Vezzaro M
10. Santacà K
11. Lora R
12. Moretti U
13. Reno C
14. Fantini MP
15. Corrao S
16. Barbato D
17. Tari M
18. Trifirò G
(2021) COVID-19 Patient Management in Outpatient Setting: A Population-Based Study from Southern Italy
Journal of Clinical Medicine 11:51.

https://doi.org/10.3390/jcm11010051
- PubMed
- Google Scholar
1. Dagan R
2. Barkai G
3. Givon-Lavi N
4. Sharf AZ
5. Vardy D
6. Cohen T
7. Lipsitch M
8. Greenberg D
(2008) Seasonality of antibiotic-resistant streptococcus pneumoniae that causes acute otitis media: a clue for an antibiotic-restriction policy?
The Journal of Infectious Diseases 197:1094–1102.

https://doi.org/10.1086/528995
- PubMed
- Google Scholar
(2023) The COVID-19 pandemic as an opportunity for unravelling the causative association between respiratory viruses and pneumococcus-associated disease in young children: a prospective study
EBioMedicine 90:104493.

https://doi.org/10.1016/j.ebiom.2023.104493
- PubMed
- Google Scholar
1. Danion F
2. Margue M
3. Ruch Y
4. Séverac F
5. Hansmann Y
(2023) Seasonal variation in azithromycin prescription
The Lancet. Infectious Diseases 23:277–278.

https://doi.org/10.1016/S1473-3099(23)00009-9
- PubMed
- Google Scholar
1. D’Arienzo M
2. Coniglio A
(2020) Assessment of the SARS-CoV-2 basic reproduction number, R0, based on the early phase of COVID-19 outbreak in Italy
Biosafety and Health 2:57–59.

https://doi.org/10.1016/j.bsheal.2020.03.004
- Google Scholar
1. Davies NG
2. Flasche S
3. Jit M
4. Atkins KE
(2019) Within-host dynamics shape antibiotic resistance in commensal bacteria
Nature Ecology & Evolution 3:440–449.

https://doi.org/10.1038/s41559-018-0786-x
- PubMed
- Google Scholar
(2010) Influenza A virus facilitates Streptococcus pneumoniae transmission and disease
FASEB Journal 24:1789–1798.

https://doi.org/10.1096/fj.09-146779
- PubMed
- Google Scholar
(2020) Impact of lockdown on COVID-19 epidemic in Île-de-France and possible exit strategies
BMC Medicine 18:240.

https://doi.org/10.1186/s12916-020-01698-4
- PubMed
- Google Scholar
1. Doan T
2. Worden L
3. Hinterwirth A
4. Arzika AM
5. Maliki R
6. Abdou A
7. Zhong L
8. Chen C
9. Cook C
10. Lebas E
11. O’Brien KS
12. Oldenburg CE
13. Chow ED
14. Porco TC
15. Lipsitch M
16. Keenan JD
17. Lietman TM
(2020) Macrolide and nonmacrolide resistance with mass azithromycin distribution
The New England Journal of Medicine 383:1941–1950.

https://doi.org/10.1056/NEJMoa2002606
- PubMed
- Google Scholar
(2019) Unraveling the seasonal epidemiology of pneumococcus
PNAS 116:1802–1807.

https://doi.org/10.1073/pnas.1812388116
- PubMed
- Google Scholar
1. Duffy E
2. Thomas M
3. Hills T
4. Ritchie S
(2021) The impacts of New Zealand’s COVID-19 epidemic response on community antibiotic use and hospitalisation for pneumonia, peritonsillar abscess and rheumatic fever
The Lancet Regional Health 12:100162.

https://doi.org/10.1016/j.lanwpc.2021.100162
- PubMed
- Google Scholar
(1997) Duration of nasopharyngeal carriage of penicillin-resistant Streptococcus pneumoniae: experiences from the South Swedish Pneumococcal Intervention Project
Clinical Infectious Diseases 25:1113–1117.

https://doi.org/10.1086/516103
- PubMed
- Google Scholar
1. Elias C
2. Sekri A
3. Leblanc P
4. Cucherat M
5. Vanhems P
(2021) The incubation period of COVID-19: A meta-analysis
International Journal of Infectious Diseases 104:708–710.

https://doi.org/10.1016/j.ijid.2021.01.069
- Google Scholar
Book
1. European Centre for Disease Prevention and Control
2. World Health Organization
(2022)
Antimicrobial Resistance Surveillance in Europe: 2022: 2020 Data

LU: Publications Office.
- Google Scholar
1. Flaxman S
2. Mishra S
3. Gandy A
4. Unwin HJT
5. Mellan TA
6. Coupland H
7. Whittaker C
8. Zhu H
9. Berah T
10. Eaton JW
11. Monod M
12. Ghani AC
13. Donnelly CA
14. Riley S
15. Vollmer MAC
16. Ferguson NM
17. Okell LC
18. Bhatt S
19. Imperial College COVID-19 Response Team
(2020) Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe
Nature 584:257–261.

https://doi.org/10.1038/s41586-020-2405-7
- PubMed
- Google Scholar
(1990) The pharmacokinetics of azithromycin in human serum and tissues
The Journal of Antimicrobial Chemotherapy 25 Suppl A:73–82.

https://doi.org/10.1093/jac/25.suppl_a.73
- PubMed
- Google Scholar
1. Garcia-Vidal C
2. Sanjuan G
3. Moreno-García E
4. Puerta-Alcalde P
5. Garcia-Pouton N
6. Chumbita M
7. Fernandez-Pittol M
8. Pitart C
9. Inciarte A
10. Bodro M
11. Morata L
12. Ambrosioni J
13. Grafia I
14. Meira F
15. Macaya I
16. Cardozo C
17. Casals C
18. Tellez A
19. Castro P
20. Marco F
21. García F
22. Mensa J
23. Martínez JA
24. Soriano A
25. COVID-19 Researchers Group
(2021) Incidence of co-infections and superinfections in hospitalized patients with COVID-19: a retrospective cohort study
Clinical Microbiology and Infection 27:83–88.

https://doi.org/10.1016/j.cmi.2020.07.041
- PubMed
- Google Scholar
1. Girard D
2. Finegan SM
3. Dunne MW
4. Lame ME
(2005) Enhanced efficacy of single-dose versus multi-dose azithromycin regimens in preclinical infection models
The Journal of Antimicrobial Chemotherapy 56:365–371.

https://doi.org/10.1093/jac/dki241
- PubMed
- Google Scholar
1. Grant J
2. Saux NL
(2021) Duration of antibiotic therapy for common infections
Journal of the Association of Medical Microbiology and Infectious Disease Canada = Journal Officiel de l’Association Pour La Microbiologie Medicale et l’infectiologie Canada 6:181–197.

https://doi.org/10.3138/jammi-2021-04-29
- PubMed
- Google Scholar
1. Henig O
2. Kehat O
3. Meijer SE
4. Chikly A
5. Weiss-Meilik A
6. Egoz E
7. Ben-Ami R
8. Paran Y
(2021) Antibiotic Use during the COVID-19 pandemic in a tertiary hospital with an ongoing antibiotic stewardship program
Antibiotics 10:1056.

https://doi.org/10.3390/antibiotics10091056
- PubMed
- Google Scholar
(2021) Decrease in community antibiotic consumption during the COVID-19 pandemic, EU/EEA, 2020
Euro Surveillance 26:2101020.

https://doi.org/10.2807/1560-7917.ES.2021.26.46.2101020
- PubMed
- Google Scholar
1. Homeniuk R
2. Collins C
(2021) How COVID-19 has affected general practice consultations and income: general practitioner cross-sectional population survey evidence from Ireland
BMJ Open 11:e044685.

https://doi.org/10.1136/bmjopen-2020-044685
- Google Scholar
1. Hsu J
(2020) How covid-19 is accelerating the threat of antimicrobial resistance
BMJ 369:m1983.

https://doi.org/10.1136/bmj.m1983
- PubMed
- Google Scholar
(2021) Impact of the COVID-19 Pandemic on the Prescribing Patterns of First-Line Antibiotics in English Primary Care: A Longitudinal Analysis of National Prescribing Dataset
Antibiotics 10:591.

https://doi.org/10.3390/antibiotics10050591
- PubMed
- Google Scholar
1. Jenness SM
2. Le Guillou A
3. Chandra C
4. Mann LM
5. Sanchez T
6. Westreich D
7. Marcus JL
(2021) Projected HIV and bacterial sexually transmitted infection incidence following COVID-19-related sexual distancing and clinical service interruption
The Journal of Infectious Diseases 223:1019–1028.

https://doi.org/10.1093/infdis/jiab051
- PubMed
- Google Scholar
(2022) Indirect effects of the covid-19 pandemic on childhood infection in England: population based observational study
BMJ 376:e067519.

https://doi.org/10.1136/bmj-2021-067519
- PubMed
- Google Scholar
(2021) Few bacterial co-infections but frequent empiric antibiotic use in the early phase of hospitalized patients with COVID-19: results from a multicentre retrospective cohort study in The Netherlands
Infectious Diseases 53:102–110.

https://doi.org/10.1080/23744235.2020.1839672
- Google Scholar
1. Keenan JD
2. Klugman KP
3. McGee L
4. Vidal JE
5. Chochua S
6. Hawkins P
7. Cevallos V
8. Gebre T
9. Tadesse Z
10. Emerson PM
11. Jorgensen JH
12. Gaynor BD
13. Lietman TM
(2015) Evidence for clonal expansion after antibiotic selection pressure: pneumococcal multilocus sequence types before and after mass azithromycin treatments
The Journal of Infectious Diseases 211:988–994.

https://doi.org/10.1093/infdis/jiu552
- PubMed
- Google Scholar
1. Keenan JD
2. Bailey RL
3. West SK
4. Arzika AM
5. Hart J
6. Weaver J
7. Kalua K
8. Mrango Z
9. Ray KJ
10. Cook C
11. Lebas E
12. O’Brien KS
13. Emerson PM
14. Porco TC
15. Leitman TM
(2018) Mass azithromycin distribution for reducing childhood mortality in sub-Saharan Africa
The New England Journal of Medicine 378:1583–1592.

https://doi.org/10.1056/NEJMoa1715474
- Google Scholar
(2022) Global consumption of antimicrobials: impact of the WHO Global Action Plan on Antimicrobial Resistance and 2019 coronavirus pandemic (COVID-19)
Journal of Antimicrobial Chemotherapy 77:1491–1499.

https://doi.org/10.1093/jac/dkac028
- Google Scholar
1. Knight GM
2. Glover RE
3. McQuaid CF
4. Olaru ID
5. Gallandat K
6. Leclerc QJ
7. Fuller NM
8. Willcocks SJ
9. Hasan R
10. van Kleef E
11. Chandler CI
(2021) Antimicrobial resistance and COVID-19: Intersections and implications
eLife 10:e64139.

https://doi.org/10.7554/eLife.64139
- PubMed
- Google Scholar
1. Kournoutou GG
2. Dinos G
(2022) Azithromycin through the Lens of the COVID-19 Treatment
Antibiotics 11:1063.

https://doi.org/10.3390/antibiotics11081063
- PubMed
- Google Scholar
Software
1. Kovacevic A
(2024) antibiotic_resistance, version swh:1:rev:2e1a778298d23145d46d8d3f702a21ec16e57996
Software Heritage.

https://archive.softwareheritage.org/swh:1:dir:db75008f81b8b1494e623f024b566ddc0e59e4a0;origin=https://github.com/alekskovacevic/antibiotic_resistance;visit=swh:1:snp:3399afe72a09fca6904c379bf83475ddb356c27e;anchor=swh:1:rev:2e1a778298d23145d46d8d3f702a21ec16e57996
(2023) Antibiotic treatment duration for community-acquired pneumonia in outpatient children in high-income countries-a systematic review and meta-analysis
Clinical Infectious Diseases 76:e1123–e1128.

https://doi.org/10.1093/cid/ciac374
- PubMed
- Google Scholar
1. Langford BJ
2. So M
3. Raybardhan S
4. Leung V
5. Westwood D
6. MacFadden DR
7. Soucy JPR
8. Daneman N
(2020) Bacterial co-infection and secondary infection in patients with COVID-19: a living rapid review and meta-analysis
Clinical Microbiology and Infection 26:1622–1629.

https://doi.org/10.1016/j.cmi.2020.07.016
- PubMed
- Google Scholar
1. Langford BJ
2. So M
3. Raybardhan S
4. Leung V
5. Soucy JPR
6. Westwood D
7. Daneman N
8. MacFadden DR
(2021) Antibiotic prescribing in patients with COVID-19: rapid review and meta-analysis
Clinical Microbiology and Infection 27:520–531.

https://doi.org/10.1016/j.cmi.2020.12.018
- PubMed
- Google Scholar
1. Lauer SA
2. Grantz KH
3. Bi Q
4. Jones FK
5. Zheng Q
6. Meredith HR
7. Azman AS
8. Reich NG
9. Lessler J
(2020) The Incubation Period of Coronavirus Disease 2019 (COVID-19) from publicly reported confirmed cases: estimation and application
Annals of Internal Medicine 172:577–582.

https://doi.org/10.7326/M20-0504
- PubMed
- Google Scholar
(2017) Evolution of antibiotic resistance is linked to any genetic mechanism affecting bacterial duration of carriage
PNAS 114:1075–1080.

https://doi.org/10.1073/pnas.1617849114
- PubMed
- Google Scholar
1. Lipsitch M
2. Colijn C
3. Cohen T
4. Hanage WP
5. Fraser C
(2009) No coexistence for free: neutral null models for multistrain pathogens
Epidemics 1:2–13.

https://doi.org/10.1016/j.epidem.2008.07.001
- PubMed
- Google Scholar
(2020) The reproductive number of COVID-19 is higher compared to SARS coronavirus
Journal of Travel Medicine 27:021.

https://doi.org/10.1093/jtm/taaa021
- Google Scholar
(2007) Effect of azithromycin and clarithromycin therapy on pharyngeal carriage of macrolide-resistant streptococci in healthy volunteers: a randomised, double-blind, placebo-controlled study
The Lancet 369:482–490.

https://doi.org/10.1016/S0140-6736(07)60235-9
- Google Scholar
1. Mamun AA
2. Saatchi A
3. Xie M
4. Lishman H
5. Blondel-Hill E
6. Marra F
7. Patrick DM
(2021) Community antibiotic use at the population level during the SARS-CoV-2 Pandemic in British Columbia, Canada
Open Forum Infectious Diseases 8:ofab185.

https://doi.org/10.1093/ofid/ofab185
- PubMed
- Google Scholar
(2010) Influenza enhances susceptibility to natural acquisition of and disease due to Streptococcus pneumoniae in ferrets
The Journal of Infectious Diseases 202:1287–1295.

https://doi.org/10.1086/656333
- PubMed
- Google Scholar
(2004) Estimating the transmission parameters of pneumococcal carriage in households
Epidemiology and Infection 132:433–441.

https://doi.org/10.1017/s0950268804001980
- PubMed
- Google Scholar
(2015) The fitness costs of antibiotic resistance mutations
Evolutionary Applications 8:273–283.

https://doi.org/10.1111/eva.12196
- PubMed
- Google Scholar
1. Monnet DL
2. Harbarth S
(2020) Will coronavirus disease (COVID-19) have an impact on antimicrobial resistance?
Euro Surveillance 25:2001886.

https://doi.org/10.2807/1560-7917.ES.2020.25.45.2001886
- PubMed
- Google Scholar
(2020) Systematic comparison of coexistence in models of drug-sensitive and drug-resistant pathogen strains
Theoretical Population Biology 133:150–158.

https://doi.org/10.1016/j.tpb.2019.12.001
- PubMed
- Google Scholar
1. Murray CJL
2. Ikuta KS
3. Sharara F
4. Swetschinski L
5. Robles Aguilar G
6. Gray A
7. Han C
8. Bisignano C
9. Rao P
10. Wool E
11. Johnson SC
12. Browne AJ
13. Chipeta MG
14. Fell F
15. Hackett S
16. Haines-Woodhouse G
17. Kashef Hamadani BH
18. Kumaran EAP
19. McManigal B
20. Achalapong S
21. Agarwal R
22. Akech S
23. Albertson S
24. Amuasi J
25. Andrews J
26. Aravkin A
27. Ashley E
28. Babin F-X
29. Bailey F
30. Baker S
31. Basnyat B
32. Bekker A
33. Bender R
34. Berkley JA
35. Bethou A
36. Bielicki J
37. Boonkasidecha S
38. Bukosia J
39. Carvalheiro C
40. Castañeda-Orjuela C
41. Chansamouth V
42. Chaurasia S
43. Chiurchiù S
44. Chowdhury F
45. Clotaire Donatien R
46. Cook AJ
47. Cooper B
48. Cressey TR
49. Criollo-Mora E
50. Cunningham M
51. Darboe S
52. Day NPJ
53. De Luca M
54. Dokova K
55. Dramowski A
56. Dunachie SJ
57. Duong Bich T
58. Eckmanns T
59. Eibach D
60. Emami A
61. Feasey N
62. Fisher-Pearson N
63. Forrest K
64. Garcia C
65. Garrett D
66. Gastmeier P
67. Giref AZ
68. Greer RC
69. Gupta V
70. Haller S
71. Haselbeck A
72. Hay SI
73. Holm M
74. Hopkins S
75. Hsia Y
76. Iregbu KC
77. Jacobs J
78. Jarovsky D
79. Javanmardi F
80. Jenney AWJ
81. Khorana M
82. Khusuwan S
83. Kissoon N
84. Kobeissi E
85. Kostyanev T
86. Krapp F
87. Krumkamp R
88. Kumar A
89. Kyu HH
90. Lim C
91. Lim K
92. Limmathurotsakul D
93. Loftus MJ
94. Lunn M
95. Ma J
96. Manoharan A
97. Marks F
98. May J
99. Mayxay M
100. Mturi N
101. Munera-Huertas T
102. Musicha P
103. Musila LA
104. Mussi-Pinhata MM
105. Naidu RN
106. Nakamura T
107. Nanavati R
108. Nangia S
109. Newton P
110. Ngoun C
111. Novotney A
112. Nwakanma D
113. Obiero CW
114. Ochoa TJ
115. Olivas-Martinez A
116. Olliaro P
117. Ooko E
118. Ortiz-Brizuela E
119. Ounchanum P
120. Pak GD
121. Paredes JL
122. Peleg AY
123. Perrone C
124. Phe T
125. Phommasone K
126. Plakkal N
127. Ponce-de-Leon A
128. Raad M
129. Ramdin T
130. Rattanavong S
131. Riddell A
132. Roberts T
133. Robotham JV
134. Roca A
135. Rosenthal VD
136. Rudd KE
137. Russell N
138. Sader HS
139. Saengchan W
140. Schnall J
141. Scott JAG
142. Seekaew S
143. Sharland M
144. Shivamallappa M
145. Sifuentes-Osornio J
146. Simpson AJ
147. Steenkeste N
148. Stewardson AJ
149. Stoeva T
150. Tasak N
151. Thaiprakong A
152. Thwaites G
153. Tigoi C
154. Turner C
155. Turner P
156. van Doorn HR
157. Velaphi S
158. Vongpradith A
159. Vongsouvath M
160. Vu H
161. Walsh T
162. Walson JL
163. Waner S
164. Wangrangsimakul T
165. Wannapinij P
166. Wozniak T
167. Young Sharma TEMW
168. Yu KC
169. Zheng P
170. Sartorius B
171. Lopez AD
172. Stergachis A
173. Moore C
174. Dolecek C
175. Naghavi M
(2022) Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis
The Lancet 399:629–655.

https://doi.org/10.1016/S0140-6736(21)02724-0
- Google Scholar
(2023) Global antibiotic use during the COVID-19 pandemic: analysis of pharmaceutical sales data from 71 countries, 2020–2022
eClinicalMedicine 57:101848.

https://doi.org/10.1016/j.eclinm.2023.101848
- Google Scholar
1. Nation ML
2. Manna S
3. Tran HP
4. Nguyen CD
5. Vy LTT
6. Uyen DY
7. Phuong TL
8. Dai VTT
9. Ortika BD
10. Wee-Hee AC
11. Beissbarth J
12. Hinds J
13. Bright K
14. Smith-Vaughan H
15. Nguyen TV
16. Mulholland K
17. Temple B
18. Satzke C
(2023) Impact of COVID-19 Nonpharmaceutical Interventions on Pneumococcal Carriage Prevalence and Density in Vietnam
Microbiology Spectrum 11:e0361522.

https://doi.org/10.1128/spectrum.03615-22
- PubMed
- Google Scholar
(2020) The role of “spillover” in antibiotic resistance
PNAS 117:29063–29068.

https://doi.org/10.1073/pnas.2013694117
- PubMed
- Google Scholar
1. Olwagen CP
2. Downs SL
3. Izu A
4. Tharasimbi L
5. Van Der Merwe L
6. Nunes MC
7. Madhi SA
(2024) Bacterial nasopharyngeal colonisation in children in South Africa before and during the COVID-19 pandemic: an observational study
The Lancet. Microbe 5:e34–e42.

https://doi.org/10.1016/S2666-5247(23)00260-4
- PubMed
- Google Scholar
1. Opatowski L
2. Varon E
3. Dupont C
4. Temime L
5. van der Werf S
6. Gutmann L
7. Boëlle P-Y
8. Watier L
9. Guillemot D
(2013) Assessing pneumococcal meningitis association with viral respiratory infections and antibiotics: insights from statistical and mathematical models
Proceedings of the Royal Society B 280:20130519.

https://doi.org/10.1098/rspb.2013.0519
- Google Scholar
1. Palmer K
2. Monaco A
3. Kivipelto M
4. Onder G
5. Maggi S
6. Michel JP
7. Prieto R
8. Sykara G
9. Donde S
(2020) The potential long-term impact of the COVID-19 outbreak on patients with non-communicable diseases in Europe: consequences for healthy ageing
Aging Clinical and Experimental Research 32:1189–1194.

https://doi.org/10.1007/s40520-020-01601-4
- Google Scholar
1. Parveen M
2. Molla MMA
3. Yeasmin M
4. Nafisa T
5. Barna AA
6. Ghosh AK
(2020) Evidences on Irrational anti-microbial prescribing and consumption among COVID-19 positive patients and possible mitigation strategies: a descriptive cross sectional study
Bangladesh Journal of Infectious Diseases 7:S3–S7.

https://doi.org/10.3329/bjid.v7i00.50155
- Google Scholar
(2022) Pneumococcal nasopharyngeal carriage in children under 5 years of age at an outpatient healthcare facility in Novi Sad, Serbia during the COVID-19 pandemic
IJID Regions 4:88–96.

https://doi.org/10.1016/j.ijregi.2022.07.001
- PubMed
- Google Scholar
1. PRINCIPLE Trial Collaborative Group
(2021) Azithromycin for community treatment of suspected COVID-19 in people at increased risk of an adverse clinical course in the UK (PRINCIPLE): a randomised, controlled, open-label, adaptive platform trial
Lancet 397:1063–1074.

https://doi.org/10.1016/S0140-6736(21)00461-X
- PubMed
- Google Scholar
(2023) Changes in public health-seeking behaviours for self-limiting respiratory tract infections across England during the COVID-19 pandemic
European Journal of Public Health 33:987–993.

https://doi.org/10.1093/eurpub/ckad136
- PubMed
- Google Scholar
(2021) COVID-19 lockdowns may reduce resistance genes diversity in the human microbiome and the need for antibiotics
International Journal of Molecular Sciences 22:6891.

https://doi.org/10.3390/ijms22136891
- PubMed
- Google Scholar
(2020) Incidence of Nosocomial COVID-19 in patients hospitalized at a large us academic medical center
JAMA Network Open 3:e2020498.

https://doi.org/10.1001/jamanetworkopen.2020.20498
- PubMed
- Google Scholar
(2021) nasopharyngeal carriage in children after the introduction of generalized infant pneumococcal conjugate vaccine immunization in Germany
Frontiers in Medicine 8:719481.

https://doi.org/10.3389/fmed.2021.719481
- Google Scholar
Preprint
(2020) COVID-19: One-Month Impact of the French Lockdown on the Epidemic Burden
bioRxiv.

https://doi.org/10.1101/2020.04.22.20075705
- Google Scholar
1. Rusic D
2. Vilovic M
3. Bukic J
4. Leskur D
5. Seselja Perisin A
6. Kumric M
7. Martinovic D
8. Petric A
9. Modun D
10. Bozic J
(2021) Implications of COVID-19 Pandemic on the Emergence of Antimicrobial Resistance: Adjusting the Response to Future Outbreaks
Life 11:220.

https://doi.org/10.3390/life11030220
- PubMed
- Google Scholar
1. Rybak A
2. Levy C
3. Angoulvant F
4. Auvrignon A
5. Gembara P
6. Danis K
7. Vaux S
8. Levy-Bruhl D
9. van der Werf S
10. Béchet S
11. Bonacorsi S
12. Assad Z
13. Lazzati A
14. Michel M
15. Kaguelidou F
16. Faye A
17. Cohen R
18. Varon E
19. Ouldali N
(2022) association of nonpharmaceutical interventions during the COVID-19 pandemic with invasive pneumococcal disease, pneumococcal carriage, and respiratory viral Infections among children in france
JAMA Network Open 5:e2218959.

https://doi.org/10.1001/jamanetworkopen.2022.18959
- PubMed
- Google Scholar
1. Ryu S
2. Hwang Y
3. Ali ST
4. Kim DS
5. Klein EY
6. Lau EHY
7. Cowling BJ
(2021) Decreased use of broad-spectrum antibiotics during the coronavirus disease 2019 epidemic in South Korea
The Journal of Infectious Diseases 224:949–955.

https://doi.org/10.1093/infdis/jiab208
- PubMed
- Google Scholar
1. Salje H
2. Tran Kiem C
3. Lefrancq N
4. Courtejoie N
5. Bosetti P
6. Paireau J
7. Andronico A
8. Hozé N
9. Richet J
10. Dubost CL
11. Le Strat Y
12. Lessler J
13. Levy-Bruhl D
14. Fontanet A
15. Opatowski L
16. Boelle PY
17. Cauchemez S
(2020) Estimating the burden of SARS-CoV-2 in France
Science 369:208–211.

https://doi.org/10.1126/science.abc3517
- PubMed
- Google Scholar
1. Seaton RA
2. Gibbons CL
3. Cooper L
4. Malcolm W
5. McKinney R
6. Dundas S
7. Griffith D
8. Jeffreys D
9. Hamilton K
10. Choo-Kang B
11. Brittain S
12. Guthrie D
13. Sneddon J
(2020) Survey of antibiotic and antifungal prescribing in patients with suspected and confirmed COVID-19 in Scottish hospitals
The Journal of Infection 81:952–960.

https://doi.org/10.1016/j.jinf.2020.09.024
- PubMed
- Google Scholar
(2021) virus-induced changes of the respiratory tract environment promote secondary infections with Streptococcus pneumoniae
Frontiers in Cellular and Infection Microbiology 11:643326.

https://doi.org/10.3389/fcimb.2021.643326
- PubMed
- Google Scholar
1. Shaw D
2. Abad R
3. Amin-Chowdhury Z
4. Bautista A
5. Bennett D
6. Broughton K
7. Cao B
8. Casanova C
9. Choi EH
10. Chu Y-W
11. Claus H
12. Coelho J
13. Corcoran M
14. Cottrell S
15. Cunney R
16. Cuypers L
17. Dalby T
18. Davies H
19. de Gouveia L
20. Deghmane A-E
21. Demczuk W
22. Desmet S
23. Domenech M
24. Drew R
25. du Plessis M
26. Duarte C
27. Erlendsdóttir H
28. Fry NK
29. Fuursted K
30. Hale T
31. Henares D
32. Henriques-Normark B
33. Hilty M
34. Hoffmann S
35. Humphreys H
36. Ip M
37. Jacobsson S
38. Johnson C
39. Johnston J
40. Jolley KA
41. Kawabata A
42. Kozakova J
43. Kristinsson KG
44. Krizova P
45. Kuch A
46. Ladhani S
47. Lâm T-T
48. León ME
49. Lindholm L
50. Litt D
51. Maiden MCJ
52. Martin I
53. Martiny D
54. Mattheus W
55. McCarthy ND
56. Meehan M
57. Meiring S
58. Mölling P
59. Morfeldt E
60. Morgan J
61. Mulhall R
62. Muñoz-Almagro C
63. Murdoch D
64. Murphy J
65. Musilek M
66. Mzabi A
67. Novakova L
68. Oftadeh S
69. Perez-Argüello A
70. Pérez-Vázquez M
71. Perrin M
72. Perry M
73. Prevost B
74. Roberts M
75. Rokney A
76. Ron M
77. Sanabria OM
78. Scott KJ
79. Sheppard C
80. Siira L
81. Sintchenko V
82. Skoczyńska A
83. Sloan M
84. Slotved H-C
85. Smith AJ
86. Steens A
87. Taha M-K
88. Toropainen M
89. Tzanakaki G
90. Vainio A
91. van der Linden MPG
92. van Sorge NM
93. Varon E
94. Vohrnova S
95. von Gottberg A
96. Yuste J
97. Zanella R
98. Zhou F
99. Brueggemann AB
(2023) Trends in invasive bacterial diseases during the first 2 years of the COVID-19 pandemic: analyses of prospective surveillance data from 30 countries and territories in the IRIS Consortium
The Lancet Digital Health 5:e582–e593.

https://doi.org/10.1016/S2589-7500(23)00108-5
- Google Scholar
(2012) Increased nasopharyngeal bacterial titers and local inflammation facilitate transmission of Streptococcus pneumoniae
mBio 3:e00255-12.

https://doi.org/10.1128/mBio.00255-12
- PubMed
- Google Scholar
1. Shrestha S
2. Foxman B
3. Dawid S
4. Aiello AE
5. Davis BM
6. Berus J
7. Rohani P
(2013) Time and dose-dependent risk of pneumococcal pneumonia following influenza: a model for within-host interaction between influenza and Streptococcus pneumoniae
Journal of The Royal Society Interface 10:20130233.

https://doi.org/10.1098/rsif.2013.0233
- Google Scholar
1. Smith DRM
2. Opatowski L
(2021) COVID-19 containment measures and incidence of invasive bacterial disease
The Lancet. Digital Health 3:e331–e332.

https://doi.org/10.1016/S2589-7500(21)00085-6
- PubMed
- Google Scholar
(2023) Collateral impacts of pandemic COVID-19 drive the nosocomial spread of antibiotic resistance: A modelling study
PLOS Medicine 20:e1004240.

https://doi.org/10.1371/journal.pmed.1004240
- PubMed
- Google Scholar
(2010) Solving differential equations in R: package deSolve
Journal of Statistical Software 33:1–25.

https://doi.org/10.18637/jss.v033.i09
- Google Scholar
1. Teng JLL
2. Fok KMN
3. Lin KPK
4. Chan E
5. Ma Y
6. Lau SKP
7. Woo PCY
(2022) Substantial decline in invasive pneumococcal disease during coronavirus disease 2019 pandemic in Hong Kong
Clinical Infectious Diseases 74:335–338.

https://doi.org/10.1093/cid/ciab382
- Google Scholar
(2023) Decreased Pneumococcal carriage among older adults in denmark during the COVID-19 Lockdown
Open Forum Infectious Diseases 10:ofad365.

https://doi.org/10.1093/ofid/ofad365
- PubMed
- Google Scholar
1. Tsay SV
2. Bartoces M
3. Gouin K
4. Kabbani S
5. Hicks LA
(2022) Antibiotic Prescriptions Associated With COVID-19 Outpatient Visits Among Medicare Beneficiaries, April 2020 to April 2021
JAMA 327:2018–2019.

https://doi.org/10.1001/jama.2022.5471
- PubMed
- Google Scholar
1. Van Laethem J
2. Wuyts S
3. Van Laere S
4. Dirkx S
5. Seyler L
6. Mertens R
7. Ilsen B
8. Lacor P
9. Pierard D
10. Allard SD
(2021) Antibiotic Prescriptions Targeting Bacterial Respiratory Infections in Admitted Patients with COVID-19: A Prospective Observational Study
Infectious Diseases and Therapy 10:2575–2591.

https://doi.org/10.1007/s40121-021-00535-2
- PubMed
- Google Scholar
1. Wang L
2. Fu J
3. Liang Z
4. Chen J
(2017) Prevalence and serotype distribution of nasopharyngeal carriage of Streptococcus pneumoniae in China: a meta-analysis
BMC Infectious Diseases 17:765.

https://doi.org/10.1186/s12879-017-2816-8
- PubMed
- Google Scholar
Book
1. Weill A
2. Drouin J
3. Desplas D
4. Cuenot F
5. Dray-Spira R
6. Zureik M
(2021)
Usage Des Médicaments de Ville En France Durant L’Épidémie de La COVID-19 – Point de Situation Jusqu’Au 25 Avril 2021. Étude Pharmaco-Épidémiologique à Partir Des Données de Remboursement Du SNDS

EPI-PHARE (Groupement d’intérêt scientifique ANSM-Cnam.
- Google Scholar
1. Willen L
2. Ekinci E
3. Cuypers L
4. Theeten H
5. Desmet S
(2021) Infant pneumococcal carriage in belgium not affected by COVID-19 containment measures
Frontiers in Cellular and Infection Microbiology 11:825427.

https://doi.org/10.3389/fcimb.2021.825427
- PubMed
- Google Scholar
(2023) Antibiotic Receipt During Outpatient Visits for COVID-19 in the US, From 2020 to 2022
JAMA Health Forum 4:e225429.

https://doi.org/10.1001/jamahealthforum.2022.5429
- PubMed
- Google Scholar
1. Wolter N
2. Tempia S
3. Cohen C
4. Madhi SA
5. Venter M
6. Moyes J
7. Walaza S
8. Malope-Kgokong B
9. Groome M
10. du Plessis M
11. Magomani V
12. Pretorius M
13. Hellferscee O
14. Dawood H
15. Kahn K
16. Variava E
17. Klugman KP
18. von Gottberg A
(2014) High nasopharyngeal pneumococcal density, increased by viral coinfection, is associated with invasive pneumococcal pneumonia
The Journal of Infectious Diseases 210:1649–1657.

https://doi.org/10.1093/infdis/jiu326
- PubMed
- Google Scholar
(2023) The interactions of SARS-CoV-2 with cocirculating pathogens: Epidemiological implications and current knowledge gaps
PLOS Pathogens 19:e1011167.

https://doi.org/10.1371/journal.ppat.1011167
- PubMed
- Google Scholar
1. Wyllie AL
2. Mbodj S
3. Thammavongsa DA
4. Hislop MS
5. Yolda-Carr D
6. Waghela P
7. Nakahata M
8. Stahlfeld AE
9. Vega NJ
10. York A
11. Allicock OM
12. Wilkins G
13. Ouyang A
14. Siqueiros L
15. Strong Y
16. Anastasio K
17. Alexander-Parrish R
18. Arguedas A
19. Gessner BD
20. Weinberger DM
(2023) Persistence of Pneumococcal carriage among older adults in the community despite COVID-19 mitigation measures
Microbiology Spectrum 11:e0487922.

https://doi.org/10.1128/spectrum.04879-22
- PubMed
- Google Scholar
1. Zhang T
2. Shen X
3. Liu R
4. Zhao L
5. Wang D
6. Lambert H
7. Cabral C
(2021) The impact of COVID-19 on primary health care and antibiotic prescribing in rural China: qualitative study
BMC Health Services Research 21:1048.

https://doi.org/10.1186/s12913-021-07082-z
- PubMed
- Google Scholar
1. Zhou C
2. Jiang Y
3. Sun L
4. Li H
5. Liu X
6. Huang L
(2023) Secondary pulmonary infection and co-infection in elderly COVID-19 patients during the pandemics in a tertiary general hospital in Beijing, China
Frontiers in Microbiology 14:1280026.

https://doi.org/10.3389/fmicb.2023.1280026
- PubMed
- Google Scholar
1. Zhu N
2. Aylin P
3. Rawson T
4. Gilchrist M
5. Majeed A
6. Holmes A
(2021) Investigating the impact of COVID-19 on primary care antibiotic prescribing in North West London across two epidemic waves
Clinical Microbiology and Infection 27:762–768.

https://doi.org/10.1016/j.cmi.2021.02.007
- PubMed
- Google Scholar

Article and author information

Author details

Aleksandra Kovacevic
1. Institut Pasteur, Université Paris Cité, Epidemiology and Modelling of Antibiotic Evasion (EMAE) unit, Paris, France
2. Université Paris-Saclay, Université de Versailles Saint-Quentin-en-Yvelines, Inserm U1018, CESP, Anti-infective evasion and pharmacoepidemiology team, Montigny-Le-Bretonneux, France
Contribution
Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing

For correspondence
aleksandra.kovacevic@pasteur.fr

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-9740-6207
David RM Smith
1. Institut Pasteur, Université Paris Cité, Epidemiology and Modelling of Antibiotic Evasion (EMAE) unit, Paris, France
2. Université Paris-Saclay, Université de Versailles Saint-Quentin-en-Yvelines, Inserm U1018, CESP, Anti-infective evasion and pharmacoepidemiology team, Montigny-Le-Bretonneux, France
3. Modélisation, épidémiologie et surveillance des risques sanitaires (MESuRS), Conservatoire national des arts et métiers, Paris, France
4. Health Economics Research Centre, Nuffield Department of Health, University of Oxford, Oxford, United Kingdom
Contribution
Conceptualization, Formal analysis, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing

Competing interests
No competing interests declared
Eve Rahbé
1. Institut Pasteur, Université Paris Cité, Epidemiology and Modelling of Antibiotic Evasion (EMAE) unit, Paris, France
2. Université Paris-Saclay, Université de Versailles Saint-Quentin-en-Yvelines, Inserm U1018, CESP, Anti-infective evasion and pharmacoepidemiology team, Montigny-Le-Bretonneux, France
Contribution
Investigation, Writing – original draft, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-3828-3910
Sophie Novelli

Université Paris-Saclay, Université de Versailles Saint-Quentin-en-Yvelines, Inserm U1018, CESP, Anti-infective evasion and pharmacoepidemiology team, Montigny-Le-Bretonneux, France

Contribution
Investigation, Writing – original draft

Competing interests
No competing interests declared
Paul Henriot
1. Modélisation, épidémiologie et surveillance des risques sanitaires (MESuRS), Conservatoire national des arts et métiers, Paris, France
2. PACRI unit, Institut Pasteur, Conservatoire national des arts et métiers, Paris, France
Contribution
Investigation, Writing – original draft, Writing – review and editing

Competing interests
No competing interests declared
Emmanuelle Varon

Centre National de Référence des Pneumocoques, Centre Hospitalier Intercommunal de Créteil, Créteil, France

Contribution
Writing – review and editing

Competing interests
No competing interests declared
Robert Cohen
1. Institut Mondor de Recherche Biomédicale-Groupe de Recherche Clinique Groupe d’Etude des Maladies Infectieuses Néonatales et Infantiles (IMRB-GRC GEMINI), Université Paris Est, 94000, Créteil, France
2. Groupe de Pathologie Infectieuse Pédiatrique (GPIP), 06200, Nice, France
3. Unité Court Séjour, Petits Nourrissons, Service de Néonatologie, Centre Hospitalier, Intercommunal de Créteil, Créteil, France
4. Association Clinique et Thérapeutique Infantile du Val-de-Marne (ACTIV), 94000, Créteil, France
5. Association Française de Pédiatrie Ambulatoire (AFPA), 45000, Orléans, France
Contribution
Writing – review and editing

Competing interests
received consulting fees from Pfizer, Sanofi, MSD, and GSK, including travel grants from Pfizer and payments from Symposium Pfizer, MSD, GSK, and Sanofi, and participated in advisory/data safety monitoring board at Pfizer, Sanofi, MSD, and GSK
Corinne Levy
1. Institut Mondor de Recherche Biomédicale-Groupe de Recherche Clinique Groupe d’Etude des Maladies Infectieuses Néonatales et Infantiles (IMRB-GRC GEMINI), Université Paris Est, 94000, Créteil, France
2. Groupe de Pathologie Infectieuse Pédiatrique (GPIP), 06200, Nice, France
3. Association Clinique et Thérapeutique Infantile du Val-de-Marne (ACTIV), 94000, Créteil, France
4. Association Française de Pédiatrie Ambulatoire (AFPA), 45000, Orléans, France
Contribution
Writing – review and editing

Competing interests
received travel grants from Pfizer and payments from Symposium Pfizer and MSD
Laura Temime
1. Modélisation, épidémiologie et surveillance des risques sanitaires (MESuRS), Conservatoire national des arts et métiers, Paris, France
2. PACRI unit, Institut Pasteur, Conservatoire national des arts et métiers, Paris, France
Contribution
Conceptualization, Funding acquisition, Methodology, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-8850-5403
Lulla Opatowski
1. Institut Pasteur, Université Paris Cité, Epidemiology and Modelling of Antibiotic Evasion (EMAE) unit, Paris, France
2. Université Paris-Saclay, Université de Versailles Saint-Quentin-en-Yvelines, Inserm U1018, CESP, Anti-infective evasion and pharmacoepidemiology team, Montigny-Le-Bretonneux, France
Contribution
Conceptualization, Supervision, Funding acquisition, Methodology, Writing – review and editing

Competing interests
received a research grant from Pfizer and Sanofi Pasteur on unrelated topics through her institution

Funding

Agence Nationale de la Recherche (ANR-10-LABX-62 IBEID)

Aleksandra Kovacevic
Lulla Opatowski

Université Paris-Saclay (AAP Covid-19 2020)

Lulla Opatowski

Agence Nationale de la Recherche (SPHINX-17-CE36-0008-01)

Laura Temime

Fondation de France (MODCOV Grant 106059)

Lulla Opatowski
Laura Temime

Institut National de la Santé et de la Recherche Médicale

Lulla Opatowski

Institut Pasteur

Lulla Opatowski

Conservatoire National des Arts et Métiers

Laura Temime

Université de Versailles Saint-Quentin-en-Yvelines

Lulla Opatowski

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

797

views
155

downloads
3

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Aleksandra Kovacevic
David RM Smith
Eve Rahbé
Sophie Novelli
Paul Henriot
Emmanuelle Varon
Robert Cohen
Corinne Levy
Laura Temime
Lulla Opatowski

(2024)

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

eLife 13:e85701.

https://doi.org/10.7554/eLife.85701

Share this article

Cite this article

Antibiotic resistance trends in invasive Streptococcus pneumoniae isolates for the years 2019 and 2020.

Figure 1—source data 1

A modelling framework describing the transmission of SARS-CoV-2 and Streptococcus pneumoniae in the community setting, in the context of both general antibiotic prescribing and azithromycin prescribing for COVID-19 infected individuals.

Five mechanisms implemented in 31 pandemic scenarios proposed to explain the reported trends of IPD incidence, antibiotic resistance, and pneumococcal carriage in S. pneumoniae.

Annual incidence of invasive pneumococcal disease (IPD), antibiotic resistance (AR%), and pneumococcal carriage prevalence for three different subpopulations.

The impact of varying SARS-CoV-2 R0 and percentage of COVID-19 infected individuals taking azithromycin in scenario S19 on antibiotic resistance (%) and the annual incidence of antibiotic-resistant invasive pneumococcal disease (IPD).

Hypothetical within-host interactions in scenario S19 promote the incidence of antibiotic-resistant invasive pneumococcal disease (IPD).

Model schematic describing co-circulation of SARS-CoV-2 infection and pneumococcal carriage transmission.

Total number of invasive S. pneumoniae isolates from blood or cerebrospinal fluid tested in 2019 and 2020 and reported to EARS-Net (European Antimicrobial Resistance Surveillance Network).

Summary of the model parameters.

Author details

Aleksandra Kovacevic

Contribution

For correspondence

Competing interests

David RM Smith

Contribution

Competing interests

Eve Rahbé

Contribution

Competing interests

Sophie Novelli

Contribution

Competing interests

Paul Henriot

Contribution

Competing interests

Emmanuelle Varon

Contribution

Competing interests

Robert Cohen

Contribution

Competing interests

Corinne Levy

Contribution

Competing interests

Laura Temime

Contribution

Competing interests

Lulla Opatowski

Contribution

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Further reading

The impact of varying SARS-CoV-2 R₀ and percentage of COVID-19 infected individuals taking azithromycin in scenario S19 on antibiotic resistance (%) and the annual incidence of antibiotic-resistant invasive pneumococcal disease (IPD).