Quantifying the population-level impact of expanded antibiotic treatment for cholera outbreak management

Sharia M. Ahmed, Cormac R. LaPrete, Iza Ciglenecki, Andrew S. Azman, Daniel T. Leung, Lindsay T. Keegan
doi: https://doi.org/10.1101/2024.11.04.24316579

Abstract

Background Since 2021, there has been a resurgence in the number of cholera cases, countries affected, and the case fatality risk. Partly due to concerns about antibiotic resistance, current cholera treatment guidelines reserve antibiotics for severely symptomatic cases, recommending only supportive care (e.g., oral rehydration) for non-severe cases. However, it has been suggested that the reduction in transmissibility from antibiotic treatment may result in circumstances under which treating mild or moderate cases with antibiotics may have population-level benefits. We developed a compartmental model of cholera transmission in a non-endemic setting to quantify the potential impact of expanded antibiotic treatment on disease burden and antibiotic use. Through simulations, we evaluated different outbreak scenarios, by varying the reproductive number, care-seeking behavior, and proportion of non-severe cases receiving antibiotics. We found that expanding antibiotic treatment could significantly reduce the final outbreak size under certain outbreak characteristics. Under these different transmission scenarios, treating non-severe symptomatic infections with antibiotics decreased cholera transmission and, in some cases, the total number of antibiotic doses used. In high transmission settings, the benefits of expanded treatment are less pronounced and the strategy may lead to increased antibiotic use, potentially increasing the risk of antibiotic resistance. We show that the effectiveness of expanded antibiotic treatment is highly dependent on achieving high care-seeking rates among non-severely symptomatic infections and tailoring the approach to specific outbreak conditions. While expanding antibiotic eligibility could enhance outbreak control in some settings, careful consideration of antibiotic resistance risks is necessary in high-transmission contexts.

Introduction

Cholera, caused by the toxigenic bacterium Vibrio cholerae O1/O139, remains a significant public health threat [1]. It is characterized by severe acute watery diarrhea that can cause death within hours if left untreated [1,2]. Despite its severe clinical manifestations, the true burden of cholera is not well characterized, as only approximately 10% of cases experience severe symptoms and it has been estimated that the vast majority of infections are asymptomatic or unreported [3–6]. A 2015 global burden analysis estimated 2.86 million cases (uncertainty range 1.3–4 million cases) and 95,000 deaths (uncertainty range 21,000–143,000 deaths) in endemic countries annually [7]. Although the global burden of cholera has been declining in recent decades, a notable resurgence in cases, countries affected, and the case fatality rate has been described since 2021 [8].

Patterns of cholera transmission vary between endemic regions and non-endemic regions, in part due to protective immunity [9]. Importantly, the estimated durability of natural immunity varies from several months to 10 years [10]. Subclinical cholera infections may confer lower protection than clinical infections [11]. As a result, in endemic regions where local transmission has been detected over the past 3 years, incidence of cholera is highest among young children, as young children are the least likely to have previous exposure and therefore immunity. Whereas in non-endemic (outbreak) settings where cholera does not regularly occur, attack rates among children and adults are similar, as the entire population has a similarly low or non-existent level of pre-existing immunity [12,13].

Cholera is typically an easily treatable disease, with rehydration serving as the cornerstone of treatment. Non-severe cholera cases can be successfully managed with oral rehydration solution and severe cases with intravenous rehydration. When administered promptly, rehydration therapy can reduce the case fatality risk to below 1% [1,6,14,15].

Antibiotics, though available for cholera treatment, are generally reserved for severely dehydrated patients or those with high-risk conditions such as pregnancy or severe acute malnutrition, in part due to concerns about antibiotic resistance [14–16]. While antibiotics reduce the duration of symptoms and stool volume, their use is not typically recommended for non-severely dehydrated cases as the benefit to the patient is modest relative to the risk of development of antibiotic resistance at the population level [15,17,18]. Antibiotics do however confer benefits related to transmissibility: infections untreated with antibiotics shed V. cholerae for up to 10 days after symptom resolution, contributing to community transmission [1]. In contrast, antibiotic treatment reduces the duration of bacterial shedding by up to 90% and likely the concentration of infectious bacteria in the stool, thereby considerably reducing cases’ transmission potential [15,16].

Currently, treatment guidelines reserve antibiotics for the most severe infections in part due to the modest benefit conferred to non-severely symptomatic infections and partly due to concerns about antibiotic resistance [15]. When indicated, the recommended treatment for cholera is a single dose of doxycycline, though azithromycin and ciprofloxacin can be used in certain situations [15]. Antibiotics are often under-regulated in areas where cholera is most likely to occur, and people frequently self-medicate with antibiotics before presenting to formal clinical care [19–21]. These self-medicated antibiotics are often inappropriate to treat cholera and further contribute to the development of antibiotic resistance [19]. Recently, it has been proposed that the reduction in transmissibility associated with antibiotic treatment could offer public health benefits by curtailing outbreak transmission [22]. Expanding antibiotic treatment guidelines to include mild and moderately symptomatic infections may result in fewer antibiotic doses used over the course of an outbreak by reducing transmission, resulting in fewer infections compared to current treatment guidelines. In this study, we develop and apply a compartmental model of cholera transmission to quantify the impact of expanding antibiotic treatment to include mild or moderately symptomatic infections. Our objective was to assess how treatment of those with moderate/some dehydration could impact the overall burden of cholera in an outbreak as well as the risk of developing antibiotic resistance.

Methods

Model

To assess the impact of expanded antibiotic treatment guidelines on cholera outbreaks, we developed and analyzed a compartmental model of cholera transmission in a non- endemic setting (Figure 1). Susceptible (S) individuals, with no immunity from prior infection or vaccination become exposed (E) upon successful transmission and then progress to one of the infectious compartments (I). We differentiate within the infectious compartments based on symptomatology: asymptomatic (IA), non-severely symptomatic infections, including both mild and moderately symptomatic (IM), (hereafter referred to collectively as “moderate”), and severely symptomatic (IS) as defined by dehydration. Additionally, we differentiate within the non-severely symptomatic infections based on care seeking behavior: seeking care (IU) or not seeking care (I%).

Figure 1:
Figure 1: Compartmental diagram of the cholera transmission model.

All individuals start as susceptible (S) and become exposed (E) at a rate λ. Exposed individuals transition to the infected (I) compartment and we differentiate by symptoms (IA, IM, IS, asymptomatic, non-severely, or severely symptomatic, respectively) and by care seeking behavior (IH,IT,, not care seeking and care seeking, respectively). All severely symptomatic infections are treated with antibiotics whereas not all non-severely symptomatic infections who seek treatment receive antibiotics. The proportion of healthcare seeking non-severely symptomatic infections who receive antibiotics is governed by q. Untreated infections, both non-severe and severe, continue to shed for a longer duration following the resolution of symptoms (IMsh, ISsh), occurring at rate αM, αS, respectively. Non-severely symptomatic infections who are treated with antibiotics continue to shed for a shorter duration following treatment (IMabx), occurring at rate δθ, whereas severely symptomatic infections who are treated with antibiotics remain in a treatment facility and do not contribute to transmission (Isabx). Compartments with a dark grey outline indicate they contribute to transmission. Infectious individuals either recover (R, at rates γA, γM, γMabx, γS, γSabx) or die (D, at rates μM, μS) and we differentiate between individuals who recover without antibiotic treatment (Run) and those who recover with antibiotics (Rabx) to compare the number of doses used under different treatment scenarios.

We do not attempt to parameterize care seeking behavior, rather, we vary the proportion of non-severely symptomatic infections who seek care between simulations. Due to the nature of their symptoms, we assume that all severely symptomatic infections who seek care (ISH) receive antibiotic treatment (ISabx) and that severely symptomatic infections are much more likely to seek treatment. We model the effect of varying the proportion of non- severely symptomatic infections that seek (IMH) and receive antibiotic treatment (IMabx).

We assume that asymptomatic infections never seek treatment and thus do not further divide this compartment by treatment.

Since infected individuals untreated with antibiotics continue to shed V. cholerae for up to 10 days post-symptom resolution [23], we assume that they continue to contribute to transmission after the conclusion of symptoms, but at a reduced rate (IMsh, ISsh, for non- severely and severely symptomatic infections, respectively). For non-severely symptomatic infections treated with antibiotics (IMabx), we assume that they contribute to transmission after receiving antibiotic treatment but at a reduced rate and for a reduced duration. In contrast, we assume that severely symptomatic infections that seek care (ISabx) are admitted to a treatment center and therefore no longer shed into the community, and thus do not contribute to transmission after receiving antibiotic treatment. In general, we parameterize our model to ensure key rates for non-severely symptomatic infections, such as treatment seeking (θ) and proportion in each symptom class (ρA, ρM, ρS) are proportional to those for severe infections. For full model equations, see supplemental material.

Scenarios

We evaluate the impact of treating mild and moderately (non-severe) symptomatic infections with antibiotics over different outbreak scenarios. These scenarios are characterized by varying the effective reproductive number, the proportion care seeking, and proportion treated with appropriate antibiotics.

Because cholera outbreaks have high variability in the reported effective reproductive number by outbreak setting, we designed three scenarios based on these different transmission characteristics at the start of the outbreak: low Re = 1.3 − 1.5 and intermediate Re = 1.6 − 2.0 based off outbreaks from Africa [26–30] and high Re = 2.3 −2.8 based off outbreaks from the Americas [24,25,31,32].

Our primary interest lies in understanding the impact of expanding antibiotic treatment eligibility guidelines to include non-severely symptomatic infections. However, the proportion of non-severely symptomatic infections that receive appropriate antibiotics is governed both by the proportion that seek care (εMH) as well as the proportion of care- seeking non-severely symptomatic infections who are treated (q). Care-seeking behavior among non-severely symptomatic infections is not well characterized and is likely influenced by the current treatment guidelines that limit antibiotics to severely symptomatic infections. As such, we explore five scenarios for care seeking: 5%, 25%, 50%, 75%, and 100% of non-severely symptomatic infections seeking care (five values of εMH).

Within the scenarios described above, we also vary the proportion of non-severely symptomatic infections who have sought care that receive antibiotic treatment (“proportion treated with antibiotics”, derived from q, see supplement) from no one (current treatment guidelines) to everyone (all care-seeking non-severely symptomatic infections receive antibiotics).

Simulations

We simulate our model using R statistical software [45] over the range of parameter values presented in Table 1 assuming a total population of 4 million, based on average cholera outbreak sizes[46,47]. Since we are modeling a non-endemic outbreak, we initialize our model with 10 exposed individuals and all others susceptible. We simulate each scenario 1000 times until stochastic extinction.

View this table:
Table 1: Summary of the model parameters (and available sources) used in model.

Many key parameters are not well characterized for cholera. As such, we use Latin hypercube sampling (LHS) to sample across the full parameter uncertainty range [48,49]. LHS ensures comprehensive exploration of the parameter space by sampling each interval only once, allowing for the distinction between variability in outcomes driven by the stochastic nature of the model from that arising due to parameter uncertainty [49]. In other words, if our range of parameter value encompasses the true parameter value, the outcome metric should also encompass the truth. We constrain all parameters to ensure that no draws from the LHS yield epidemiologically implausible values.

Epidemiologic Outcomes

We evaluate the scenarios across several metrics, including the total number of infections (final size), the total number of antibiotic doses used over the course of the outbreak, the number of infections under expanded treatment minus the number of infections under current treatment guidelines (infections averted), as well as the number of doses used under expanded treatment minus the number of doses used under current treatment guidelines (additional doses used).

Results

Through simulation, we show that expanding antibiotic treatment guidelines to include non-severely symptomatic infections can substantially reduce the burden of cholera in low and intermediate transmission settings, especially when rates of care-seeking behavior is high. In high transmission settings, the impact of expanded antibiotic treatment guidelines is less pronounced but still offers marginal reductions in the final outbreak size.

We find that treating non-severely symptomatic infections with antibiotics decreases the final size of the outbreak across all care seeking scenarios and transmission scenarios except in high Re outbreaks with very low proportion of care seeking (5%) among non- severely symptomatic infections (Figure 2). While treating non-severely symptomatic infections with antibiotics almost always decreases the final outbreak size, the effect is largest when more non-severely symptomatic infections seek care and receive appropriate antibiotics as well as when the effective reproductive number is low.

Figure 2:
Figure 2:

Plot of the final outbreak size by the proportion of care-seeking non-severely symptomatic infections treated with antibiotics. Each plot shows the final proportion of the population infected by the proportion of care-seeking non-severely symptomatic infections who receive treatment for low (yellow), intermediate (pink), and high (purple) effective reproductive numbers, for a different percent of non-severely symptomatic infections who seek treatment (5%, 25%, 50%, 75%, 100%). The solid line indicates the mean estimate under expanded treatment guidelines, the shaded region represents the 25% and 75% quantiles, and the dashed line shows the final size of the outbreak under current antibiotic treatment guidelines (treating no non-severely symptomatic infections).

While treating non-severely symptomatic infections with antibiotics can substantially reduce cholera transmission, we show that expanded antibiotic eligibility alone is not sufficient to significantly reduce transmission unless it is coupled with high care-seeking rates.

The development of antibiotic resistance remains a significant concern that informs current cholera treatment guidelines. Although we did not explicitly model the evolution of resistance, we use the number of antibiotic doses administered as a proxy for selective pressure. As such, we evaluated our scenarios based on the number of antibiotic doses used and categorized them by the scale of impact: isolated, expanding, and broad. Our simulations reveal three distinct regions within the range of tested parameters (Figure 3). In one region, extending antibiotic treatment guidelines to include non-severely symptomatic infections confer only isolated benefits, as using additional doses does not reduce transmission sufficiently to avert more than one infection over the course of the outbreak (yellow points, Figure 3). The next region is such that extending antibiotic treatment criteria provides an expanding benefit, whereby treating non-severely symptomatic infections with antibiotics reduces transmission sufficiently such that each additional dose averts more than one infection (green points, Figure 3). In the final region extending antibiotic treatment criteria results in broad benefits where fewer antibiotic doses are used over the course of an outbreak compared to current treatment practices (purple points, Figure 3). Which region a simulation falls into is governed by both Re and the proportion of non-severely symptomatic infections treated with antibiotics. As Re increases, the benefits of treating non-severely symptomatic infections with antibiotics shifts from population-level benefits and antibiotic dose reduction to primarily individual-level benefits. Likewise, as the proportion of non-severely symptomatic infections treated with antibiotics increases, the population-level benefit of expanded antibiotic use becomes more pronounced.

Figure 3:
Figure 3:

Plot of the population-level impact of expanded antibiotic treatment guidelines. Each plot compares the relative infections averted per population to the relative antibiotic usage per population. The relative reduction in infections is calculated by the ratio of the number of infections in each expanded antibiotic treatment scenario presented compared to the simulation of current antibiotic treatment guidelines (treating no non-severely symptomatic infections with antibiotics) using the same LHS sampled parameters, normalized by population size. Similarly, the relative antibiotic usage is calculated by the ratio of antibiotic doses used in each scenario presented compared to the simulation of current antibiotic treatment guidelines (treating no non-severely symptomatic infections with antibiotics) using the same LHS sampled parameters, normalized by population size. Each plot represents a different proportion of non-severely symptomatic infections seeking care who receive antibiotic treatment (5%, 25%, 50%, 75%, 100%) and a different Re scenario (low (Re = 1.3 − 1.5), intermediate (Re = 1.6 − 2.0), high (Re = 2.3 − 2.8)). The outcomes are split into three regions by the impact of expanded eligibility criteria: between the dashed line along the x-axis (expanding criteria averts no infections) and the diagonal dashed line (each addition dose of antibiotics deployed averts one infection), yellow points represent simulations in which expanded eligibility only has isolated benefits; between the diagonal dashed line and the vertical dashed line (no additional doses are used to avert infections), green points represent simulations in which expanded eligibility results in each dose preventing more than one additional infection; and in the region left of the vertical dashed line, purple points represent simulations in which expanded eligibility results in fewer doses used over the course of the outbreak than compared to current antibiotic treatment guidelines.

To further explore the impact of expanded treatment on the number of antibiotic doses used over the course of the outbreak, we evaluate how the proportion of care-seeking non- severe individuals who receive appropriate antibiotics impacts the proportion of the total population who receive appropriate antibiotics (Figure 4). As with infections averted, when the transmission rate is lower and when more non-severely symptomatic infections are treated, treating non-severely symptomatic infections with antibiotics has greater public health benefits. Indeed, when Re is low, increasing the proportion of non-severely symptomatic infections treated with antibiotics results in fewer antibiotic doses used. For intermediate and high transmission settings, increasing the proportion of non-severely symptomatic infections treated with antibiotics only results in an increase in the proportion of the total population treated with antibiotics.

Figure 4:
Figure 4:

Plot of population-level antibiotic treatment. Each plot shows the proportion of the total population who receive appropriate antibiotics, by the proportion of care- seeking non-severely symptomatic infections that receive appropriate antibiotics, separated by Rescenario (low (Re = 1.3 − 1.5), intermediate (Re = 1.6 − 2.0), high (Re = 2.3 − 2.8)). For each Re scenario, we show five values for the proportion of non-severely symptomatic infections that seek care, 5% (yellow), 25% (green), 50% (teal), 75% (blue), 100% (navy), and vary the proportion of care-seeking non-severely symptomatic infections who receive antibiotics from 0 – 100%.

Discussion

This study builds on prior research suggesting that, under certain circumstances, expanding antibiotic treatment of cholera cases can provide population-level benefits by reducing transmission and limiting an outbreak size such that fewer antibiotic doses are used over the course of the outbreak. While previous research proposed a theoretical mechanism, our work explores the specific conditions where these benefits can be achieved for cholera. Our findings show how treating non-severe cholera cases with antibiotics can help control outbreaks and highlight the critical role of linking these cases to healthcare in achieving impact. We detail how outbreak specific quantities such as the effective reproductive number and the proportion of non-severely symptomatic infections that seek care impact the effectiveness of this strategy.

The effectiveness of expanding antibiotic treatment criteria improves as the proportion of non-severely symptomatic infections seeking care increases. In the majority of outbreaks with low to intermediate transmission rates (measured by the reproductive number), such as those seen in most cholera outbreaks [21–25], reducing transmission can substantially lower the proportion of the population infected, supporting the theoretical findings that expanded antibiotic access can enhance outbreak control.

A key distinction emerges between low and high Re settings. In lower Resettings, treating a high proportion of non-severely symptomatic infections reduces both cholera burden and number of antibiotic doses used over the course of the outbreak. However, as Re increases, these benefits disappear, though an intermediate range exists where each additional dose of antibiotics used averts more than one infection, still providing a public health benefit. In high Resettings, treating non-severely symptomatic infections may reduce the outbreak size, but it requires treating a larger portion of the population. These results suggest that treating non-severe cholera cases may be effective in low Re settings but may carry a risk of antibiotic resistance in high-transmission contexts, underscoring the importance of tailoring treatment strategies to the specific outbreak context. In this paper, the only intervention we consider is expanding antibiotic treatment guidelines, however this could be combined with other interventions as a component in a large outbreak containment strategy.

Across all simulations, the most important factor governing the effectiveness of treating non-severely symptomatic infections with antibiotics is achieving high antibiotic treatment rates. Often the focus on clinical care during outbreaks is at cholera treatment centers, typically located at secondary and tertiary health structures. While there are gaps in our understanding of care seeking behavior for non-severely symptomatic infections, experience in endemic settings suggest that most people with non-severe disease will seek care at pharmacies or lower-level health facilities, frequently obtaining inappropriate or non-indicated antibiotics [20,21,50–52]. Developing strategies to expand access to appropriate antibiotics could have a critical impact on achieving sufficiently high coverage of antibiotics among non-severe cases to reduce cholera transmission while also improving antibiotic stewardship. Expanded access could be coupled with the use of cholera rapid diagnostic tests, which are typically not used for clinical decision making, to improve the specificity of targeting people with diarrhea caused by Vibrio cholerae O1.

Addressing barriers to care seeking behavior is necessary to achieve sufficiently high treatment rates. Cholera case-area targeted interventions are an active and promising area of research and may offer a strategy to achieve these desired antibiotic treatment rates [53–55].

Some of the limitations of our study include that many key model parameters are not well characterized for cholera. Accurate estimates of, in particular, Re and the proportion of non-severely symptomatic infections seeking treatment are crucial. Our model only considers non-endemic settings, assuming no prior immunity. Further analysis is needed to assess how prior immunity or vaccination affects our results. While we do not model immunity directly, in a supplementary analysis, we evaluate the impact of varying the ratio of asymptomatic to symptomatic infections, one of the expected main differences between an outbreak setting to an endemic setting. We find that our final size results are robust to variations in the proportion asymptomatic, but for outbreaks with a very high proportion of infections that are asymptomatic (consistent with highly endemic settings), the broad impacts from expanding antibiotic eligibility disappear (Figure SI1). Conversely, for outbreaks with a very low proportion of infections that are asymptomatic (consistent with totally naïve settings), not only are broad and expanded impacts more pronounced (Figure SI2), in low Re outbreak settings, expanding antibiotic treatment guidelines has such a substantial benefit that it halts transmission so effectively that it can prevent an outbreak from ever taking off (Figure SI3). Thus, we expect expanded treatment guidelines may have a smaller impact in endemic settings than outbreak settings.

Additionally, we only consider three symptom categories – severely symptomatic, non- severely symptomatic, and asymptomatic infections. A more granular distinction within the non-severely symptomatic compartment may refine treatment estimates. Since infectiousness correlates with symptom severity [5], it may suffice to target the most symptomatic non-severely symptomatic infections. Achieving a high enough treatment coverage may be a barrier to utilizing this strategy as an outbreak containment strategy, particularly in resource-limited settings where both antibiotics and healthcare infrastructure may be constrained.

Our study builds on phenomenological findings which suggest that treating non-severely symptomatic infections can reduce the number of antibiotic doses needed and may therefore lower the selective pressure for cholera to develop resistance. However, uncertainties in parameter estimates raise concerns about balancing these benefits with the risks of promoting antimicrobial resistance. Future studies aimed at improving characterization of key parameters or a clinical trial aimed at validating these model results and assessing the feasibility of implementing expanded treatment guidelines are a critical next step towards implementing these findings. Careful monitoring of key parameters and resistance patterns will be crucial in ensuring that the benefits of expanded antibiotic treatment criteria do not increase the risk of developing antibiotic resistant cholera.

Data Availability

This manuscript only uses publicly available data.

Funding

SMA, LTK acknowledge funding from the CDC Centers for Forecasting and Outbreak Analytics (CFA) (Grant #1NU38FT000009-01-00). CRL, LTK acknowledge funding from the Centers for Disease Control and Prevention (CDC) TRANSMIT: Training Research Acumen in Students Modeling Infectious Threats (Grant #1U01CK000675). ASA acknowledges funding from the National Institutes of Health (Grant #1R01AI179917-01A1) and funding from the Bill and Melinda Gates Foundation (Grant # INV-044856).

Supplemental Information

Methods

Model Equations

Proportion of care-seeking non-severely symptomatic infections that receive antibiotics

As derived from LaPrete et al. [22], the proportion of care-seeking non-severely symptomatic infections that receive antibiotics, Mabx, is given by:

Embedded Image
View this table:
Table S1:

Summary of the transmission parameters (and available sources) used in model.

Assessment of the ratio of asymptomatic to symptomatic infections: Although our study aims to explore non-endemic settings, since the ratio of asymptomatic to symptomatic infections varies by location, we conducted an additional sensitivity analysis to quantify the impact of the ratio of asymptomatic to symptomatic infections on disease burden and number of antibiotic doses used under expanded antibiotic treatment guidelines. Our results in the main text show asymptomatic infections representing 65%-85% of all infections. Here we explore increasing this ratio to 90%-99%, as was reported in Bangladesh[56] and lowering it to 17%-33%, as was found in Haiti [57].

Results

In a supplementary analysis, we vary the ratio of asymptomatic to symptomatic infections from 65%-85% up to 90%-99%, as well as down to 17%-33%. Overall, we find that our final size results are robust to changes in the ratio of asymptomatic to symptomatic infections and that treating non-severely symptomatic infections reduces the final size regardless of this ratio. However, we find that the broad impacts from expanded eligibility disappear for very high proportion asymptomatic (90%-99%) (Figure SI1). Simultaneously, for very high proportion of asymptomatic, we find that the total proportion of the population receiving antibiotics is very small, since most infections are asymptomatic and therefore do not seek treatment. Conversely, when asymptomatic infections represent a smaller proportion of the population (17%-33%), the proportion of simulations with expanded or broad impacts from expanded eligibility become more pronounced (Figure SI2). Indeed, for some simulations, the benefit is so substantial that it can halt the outbreak before it can take off (Figure SI3).

Figure SI1:
Figure SI1: Plot of the population-level impact of expanded antibiotic treatment guidelines with 90%-99% asymptomatic infections.

Each plot compares the relative infections averted per population to the relative antibiotic usage per population. The relative reduction in infections is calculated by the ratio of the number of infections in each expanded antibiotic treatment scenario presented compared to the simulation of current antibiotic treatment guidelines (treating no non-severely symptomatic infections with antibiotics) using the same LHS sampled parameters, normalized by population size. Similarly, the relative antibiotic usage is calculated by the ratio of antibiotic doses used in each scenario presented compared to the simulation of current antibiotic treatment guidelines (treating no non-severely symptomatic infections with antibiotics) using the same LHS sampled parameters, normalized by population size. Each plot represents a different proportion of non-severely symptomatic infections seeking care who receive antibiotic treatment (5%, 25%, 50%, 75%, 100%) and a different Re scenario (low (Re = 1.3 − 1.5), intermediate (Re = 1.6 − 2.0), high (Re = 2.3 − 2.8)). The outcomes are split into three regions by the impact of expanded eligibility criteria: between the dashed line along the x-axis (expanding criteria averts no infections) and the diagonal dashed line (each addition dose of antibiotics deployed averts one infection), yellow points represent simulations in which expanded eligibility only has isolated benefits; between the diagonal dashed line and the vertical dashed line (no additional doses are used to avert infections), green points represent simulations in which expanded eligibility results in each dose preventing more than one additional infection; and in the region left of the vertical dashed line, purple points represent simulations in which expanded eligibility results in fewer doses used over the course of the outbreak than compared to current antibiotic treatment guidelines.

Figure SI2:
Figure SI2: Plot of the population-level impact of expanded antibiotic treatment guidelines with 17%-33% asymptomatic infections.

Figure 3: Plot of the population- level impact of expanded antibiotic treatment guidelines. Each plot compares the relative infections averted per population to the relative antibiotic usage per population. The relative reduction in infections is calculated by the ratio of the number of infections in each expanded antibiotic treatment scenario presented compared to the simulation of current antibiotic treatment guidelines (treating no non-severely symptomatic infections with antibiotics) using the same LHS sampled parameters, normalized by population size. Similarly, the relative antibiotic usage is calculated by the ratio of antibiotic doses used in each scenario presented compared to the simulation of current antibiotic treatment guidelines (treating no non-severely symptomatic infections with antibiotics) using the same LHS sampled parameters, normalized by population size. Each plot represents a different proportion of non-severely symptomatic infections seeking care who receive antibiotic treatment (5%, 25%, 50%, 75%, 100%) and a different Rescenario (low (Re =1.3 − 1.5), intermediate (Re = 1.6 − 2.0), high (Re = 2.3 − 2.8)). The outcomes are split into three regions by the impact of expanded eligibility criteria: between the dashed line along the x-axis (expanding criteria averts no infections) and the diagonal dashed line (each addition dose of antibiotics deployed averts one infection), yellow points represent simulations in which expanded eligibility only has isolated benefits; between the diagonal dashed line and the vertical dashed line (no additional doses are used to avert infections), green points represent simulations in which expanded eligibility results in each dose preventing more than one additional infection; and in the region left of the vertical dashed line, purple points represent simulations in which expanded eligibility results in fewer doses used over the course of the outbreak than compared to current antibiotic treatment guidelines.

Figure SI3:
Figure SI3:

Plot of the final outbreak size by the proportion of care-seeking non- severely symptomatic infections treated with appropriate antibiotics for small proportion asymptomatic (17%-33%). Each plot shows the final proportion of the population infected by the proportion of care-seeking non-severely symptomatic infections who receive appropriate antibiotic treatment for low (yellow), intermediate (pink), and high (purple) effective reproductive numbers, for a different percent of non- severely symptomatic infections who seek care (5%, 25%, 50%, 75%, 100%). The solid line indicates the mean estimate under expanded antibiotic treatment guidelines, the shaded region represents the 25% and 75% quantiles, and the dashed line shows the final size of the outbreak under current antibiotic treatment guidelines (treating no non-severely symptomatic infections).

Acknowledgements

The authors thank Kate Alberti for helpful discussions regarding framing the research project.

References

  1. 1.↵
    World Health Organization. Cholera Fact Sheet. [cited 23 Sep 2024]. Available: https://www.who.int/news-room/fact-sheets/detail/cholera
  2. 2.↵
    Kanungo S, Azman AS, Ramamurthy T, Deen J, Dutta S. Cholera. The Lancet. 2022;399: 1429–1440. doi:10.1016/S0140-6736(22)00330-0/ASSET/68E02E63-1885-4E8D-9ADB-4E6BFE4C3044/MAIN.ASSETS/GR2.JPG
  3. 3.↵
    Azman AS, Lauer SA, Bhuiyan TR, Luquero FJ, Leung DT, Hegde ST, et al. Vibrio cholerae O1 transmission in Bangladesh: insights from a nationally representative serosurvey. Lancet Microbe. 2020;1: e336–e343. doi:10.1016/S2666-5247(20)30141-5
    OpenUrlCrossRef
  4. 4.
    Mosley WH, Benenson AS, Barui R. A serological survey for cholera antibodies in rural East Pakistan. Bull World Health Organ. 1968;38: 335. Available:/pmc/articles/PMC2554488/?report=abstract
  5. 5.↵
    Weil AA, Begum Y, Chowdhury F, Khan AI, Leung DT, LaRocque RC, et al. Bacterial Shedding in Household Contacts of Cholera Patients in Dhaka, Bangladesh. Am J Trop Med Hyg. 2014;91: 738. doi:10.4269/AJTMH.14-0095
    OpenUrlCrossRef
  6. 6.↵
  7. 7.↵
    Ali M, Nelson AR, Lopez AL, Sack DA. Updated Global Burden of Cholera in Endemic Countries. PLoS Negl Trop Dis. 2015;9: e0003832. doi:10.1371/JOURNAL.PNTD.0003832
    OpenUrlCrossRef
  8. 8.↵
    Cholera – Global situation. [cited 23 Sep 2024]. Available: https://www.who.int/emergencies/disease-outbreak-news/item/2022-DON426
  9. 9.↵
    Cholera is endangering children globally | UNICEF. [cited 23 Sep 2024]. Available: https://www.unicef.org/stories/cholera-is-endangering-children-globally
  10. 10.↵
    Hoque M, Deok M, Kim R, Ahmmed F, Islam T, Im J, et al. Sustained Protection Conferred by Natural Cholera against Recurrent Cholera in an Endemic Setting. Lancet Microbe.
  11. 11.↵
    Ali M, Emch M, Park JK, Yunus M, Clemens J. Natural Cholera Infection–Derived Immunity in an Endemic Setting. J Infect Dis. 2011;204: 912–918. doi:10.1093/INFDIS/JIR416
    OpenUrlCrossRefPubMed
  12. 12.↵
    Deen JL, von Seidlein L, Sur D, Agtini M, Lucas MES, Lopez L, et al. The High Burden of Cholera in Children: Comparison of Incidence from Endemic Areas in Asia and Africa. PLoS Negl Trop Dis. 2008;2. doi:10.1371/JOURNAL.PNTD.0000173
    OpenUrlCrossRef
  13. 13.↵
    Final Report of the Independent Panel of Experts on the Cholera Outbreak in Haiti - Haiti | ReliefWeb. [cited 23 Sep 2024]. Available: https://reliefweb.int/report/haiti/final-report-independent-panel-experts-cholera-outbreak-haiti
  14. 14.↵
    Global Task Force on Cholera Control. Cholera Outbreak Response Field Manual – Global Task Force on Cholera Control. Available: https://www.gtfcc.org/resources/gtfcc-cholera-outbreak-response-field-manual/
  15. 15.↵
    Global Task Force on Cholera Control. Interim Technical Note Use of antibiotics for the treatment and control of cholera. 2022 Oct. Available: https://www.gtfcc.org/wp-content/uploads/2019/10/gtfcc-technical-note-on-use-of-antibiotics-for-the-treatment-of-cholera-1.pdf
  16. 16.↵
    Leibovici-Weissman Y, Neuberger A, Bitterman R, Sinclair D, Salam MA, Paul M. Antimicrobial drugs for treating cholera. Cochrane Database Syst Rev. 2014;2014. doi:10.1002/14651858.CD008625.PUB2
    OpenUrlCrossRef
  17. 17.↵
    Magill SS, Edwards JR, Beldavs ZG, Dumyati G, Janelle SJ, Kainer MA, et al. Prevalence of antimicrobial use in US acute care hospitals, May-September 2011. JAMA. 2014;312: 1438–1446. doi:10.1001/JAMA.2014.12923
    OpenUrlCrossRef
  18. 18.↵
    Mo Y, Oonsivilai M, Lim C, Niehus R, Cooper BS. Implications of reducing antibiotic treatment duration for antimicrobial resistance in hospital settings: A modelling study and meta-analysis. PLoS Med. 2023;20: e1004013. doi:10.1371/JOURNAL.PMED.1004013
    OpenUrlCrossRefPubMed
  19. 19.↵
    Alexandrova L, Haque F, Rodriguez P, Marrazzo AC, Grembi JA, Ramachandran V, et al. Identification of Widespread Antibiotic Exposure in Patients With Cholera Correlates With Clinically Relevant Microbiota Changes. J Infect Dis. 2019;220: 1655–1666. doi:10.1093/INFDIS/JIZ299
    OpenUrlCrossRefPubMed
  20. 20.↵
    Islam MT, Hegde ST, Khan AI, Bhuiyan MTR, Khan ZH, Ahmmed F, et al. National Hospital-Based Sentinel Surveillance for Cholera in Bangladesh: Epidemiological Results from 2014 to 2021. Am J Trop Med Hyg. 2023;109: 575–583. doi:10.4269/AJTMH.23-0074
    OpenUrlCrossRefPubMed
  21. 21.↵
    Mannan A, Chakma K, Dewan G, Saha A, Chy NUHA, Mehedi HMH, et al. Prevalence and determinants of antibiotics self-medication among indigenous people of Bangladesh: a cross-sectional study. BMJ Open. 2024;14. doi:10.1136/BMJOPEN-2022-071504
    OpenUrlCrossRef
  22. 22.↵
    Laprete CR, Ahmed SM, Toth DJA, Reimer JR, Vaughn VM, Adler FR, et al. A theoretical framework to quantify the tradeoff between individual and population benefits of expanded antibiotic use. 2024. doi:10.1101/2024.08.28.24312731
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    World Health Organization. Cholera Fact Sheet. 2024 Sep. Available: https://www.who.int/news-room/fact-sheets/detail/cholera
  24. 24.↵
    Smirnova A, Sterrett N, Mujica OJ, Munayco C, Suárez L, Viboud C, et al. Spatial dynamics and the basic reproduction number of the 1991–1997 Cholera epidemic in Peru. PLoS Negl Trop Dis. 2020;14: 1–22. doi:10.1371/JOURNAL.PNTD.0008045
    OpenUrlCrossRef
  25. 25.↵
    Chao DL, Halloran ME, Longini IM. Vaccination strategies for epidemic cholera in Haiti with implications for the developing world. Proc Natl Acad Sci U S A. 2011;108: 7081– 7085. doi:10.1073/PNAS.1102149108/-/DCSUPPLEMENTAL
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    Che EN, Kang Y, Yakubu AA. Risk structured model of cholera infections in Cameroon. Math Biosci. 2020;320. doi:10.1016/J.MBS.2019.108303
    OpenUrlCrossRef
  27. 27.
    Mukandavire Z, Liao S, Wang J, Gaff H, Smith DL, Morris JG. Estimating the reproductive numbers for the 2008-2009 cholera outbreaks in Zimbabwe. Proc Natl Acad Sci U S A. 2011;108: 8767–8772. doi:10.1073/PNAS.1019712108/-/DCSUPPLEMENTAL/PNAS.201019712SI.PDF
    OpenUrlAbstract/FREE Full Text
  28. 28.
    Peak CM, Reilly AL, Azman AS, Buckee CO. Prolonging herd immunity to cholera via vaccination: Accounting for human mobility and waning vaccine effects. PLoS Negl Trop Dis. 2018;12: e0006257. doi:10.1371/JOURNAL.PNTD.0006257
    OpenUrlCrossRefPubMed
  29. 29.
    Yang C, Wang J. A cholera transmission model incorporating the impact of medical resources. Math Biosci Eng. 2019;16: 5226–5246. doi:10.3934/MBE.2019261
    OpenUrlCrossRefPubMed
  30. 30.↵
    Mukandavire Z, Manangazira P, Nyabadza F, Cuadros DF, Musuka G, Morris JG. Stemming cholera tides in Zimbabwe through mass vaccination. Int J Infect Dis. 2020;96: 222–227. doi:10.1016/J.IJID.2020.03.077
    OpenUrlCrossRefPubMed
  31. 31.↵
    Mukandavire Z, Smith DL, Morris JG. Cholera in Haiti: Reproductive numbers and vaccination coverage estimates. Scientific Reports 2013 3:1. 2013;3: 1–8. doi:10.1038/srep00997
    OpenUrlCrossRefPubMed
  32. 32.↵
    Tuite AR, Tien J, Eisenberg M, Earn DJD, Ma J, Fisman DN. Cholera epidemic in Haiti, 2010: using a transmission model to explain spatial spread of disease and identify optimal control interventions. Ann Intern Med. 2011;154: 593–601. doi:10.7326/0003-4819-154-9-201105030-00334
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.
    Azman AS, Rudolph KE, Cummings DAT, Lessler J. The incubation period of cholera: a systematic review. J Infect. 2013;66: 432–438. doi:10.1016/J.JINF.2012.11.013
    OpenUrlCrossRefPubMed
  34. 34.
    Weil AA, Khan AI, Chowdhury F, LaRocque RC, Faruque ASG, Ryan ET, et al. Clinical outcomes in household contacts of patients with cholera in Bangladesh. Clin Infect Dis. 2009;49: 1473–1479. doi:10.1086/644779
    OpenUrlCrossRefPubMed
  35. 35.
    Hartley DM, Morris JG, Smith DL. Hyperinfectivity: A Critical Element in the Ability of V. cholerae to Cause Epidemics? PLoS Med. 2006;3: 63–69. doi:10.1371/JOURNAL.PMED.0030007
    OpenUrlCrossRefWeb of Science
  36. 36.
    Response to Haiti’s Earthquake in 2010 - PAHO/WHO | Pan American Health Organization. [cited 23 Sep 2024]. Available: https://www.paho.org/en/response-haitis-earthquake-2010
  37. 37.
    Bertuzzo E, Mari L, Righetto L, Gatto M, Casagrandi R, Blokesch M, et al. Prediction of the spatial evolution and effects of control measures for the unfolding Haiti cholera outbreak. Geophys Res Lett. 2011;38. doi:10.1029/2011GL046823
    OpenUrlCrossRef
  38. 38.
    Meszaros VA, Miller-Dickson MD, Junior FBA, Almagro-Moreno S, Brandon Ogbunugafor C. Direct transmission via households informs models of disease and intervention dynamics in cholera. PLoS One. 2020;15. doi:10.1371/JOURNAL.PONE.0229837
    OpenUrlCrossRef
  39. 39.
    Greenough WB, Rosenberg IS, Gordon RS, Davies BI, Benenson AS. TETRACYCLINE IN THE TREATMENT OF CHOLERA. Lancet. 1964;1: 355–357. doi:10.1016/S0140-6736(64)92099-9
    OpenUrlCrossRefPubMedWeb of Science
  40. 40.
    Lindenbaum J, Greenough WB, Islam MR. Antibiotic therapy of cholera in children. Bull World Health Organ. 1967;37: 529. Available: /pmc/articles/PMC2554369/?report=abstract
  41. 41.
    Wallace CK, Anderson PN, Brown TC, Khanra SR, Lewis GW, Pierce NF, et al. Optimal antibiotic therapy in cholera. Bull World Health Organ. 1968;39: 239. Available: /pmc/articles/PMC2554550/?report=abstract
  42. 42.
    Lee EC, Kelly MR, Ochocki BM, Akinwumi SM, Hamre KES, Tien JH, et al. Model distinguishability and inference robustness in mechanisms of cholera transmission and loss of immunity. J Theor Biol. 2017;420: 68–81. doi:10.1016/J.JTBI.2017.01.032
    OpenUrlCrossRefPubMed
  43. 43.
    Dangbé E, Irépran D, Perasso A, Békollé D. Mathematical modelling and numerical simulations of the influence of hygiene and seasons on the spread of cholera. Math Biosci. 2018;296: 60–70. doi:10.1016/J.MBS.2017.12.004
    OpenUrlCrossRefPubMed
  44. 44.
    Meszaros VA, Miller-Dickson MD, Junior FBA, Almagro-Moreno S, Brandon Ogbunugafor C. Direct transmission via households informs models of disease and intervention dynamics in cholera. PLoS One. 2020;15: e0229837. doi:10.1371/JOURNAL.PONE.0229837
    OpenUrlCrossRef
  45. 45.↵
    R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria; 2021. Available: https://www.R-project.org/
  46. 46.↵
    Barkat A, Majid M, Mahiyuddin G, Poddar A, Tanvir K, Faisal H, et al. REPORT ON SURVEY FOR PROJECT AREAS IN URBAN SLUMS OF MIRPUR AND KARAIL, DHAKA United Nations Children Fund (UNICEF) Dhaka, Bangladesh. [cited 23 Sep 2024]. Available: www.hdrc-bd.com
  47. 47.↵
    Hasan MZ, Hasan AMR, Rabbani MG, Selim MA, Mahmood SS. Knowledge, attitude, and practice of Bangladeshi urban slum dwellers towards COVID-19 transmission- prevention: A cross-sectional study. PLOS Global Public Health. 2022;2. doi:10.1371/JOURNAL.PGPH.0001017
    OpenUrlCrossRef
  48. 48.↵
    Drake JM, Rohani P. Sensitivity analysis of deterministic models through Latin hypercube sampling: A model for transmission of HIV among homosexual men * 1 A model for the transmission of HIV among homosexual men.
  49. 49.↵
    Wu J, Dhingra R, Gambhir M, Remais J V. Sensitivity analysis of infectious disease models: methods, advances and their application. J R Soc Interface. 2013;10. doi:10.1098/RSIF.2012.1018
    OpenUrlCrossRef
  50. 50.↵
    Biswas M, Roy MN, Manik MIN, Hossain MS, Tapu STA, Moniruzzaman M, et al. Self medicated antibiotics in Bangladesh: a cross-sectional health survey conducted in the Rajshahi City. BMC Public Health. 2014;14. doi:10.1186/1471-2458-14-847
    OpenUrlCrossRefPubMed
  51. 51.
    Jamil S, Shariff S, Kantawala B, Satapathy P, Nazli Khatib M, Gaidhane S, et al. Self- medication in Bangladesh: significance of drug control and distribution by pharmacists. International Journal of Surgery: Global Health. 2024;7. doi:10.1097/GH9.0000000000000420
    OpenUrlCrossRef
  52. 52.↵
    Alexandrova L, Haque F, Rodriguez P, Marrazzo AC, Grembi JA, Ramachandran V, et al. Identification of Widespread Antibiotic Exposure in Patients With Cholera Correlates With Clinically Relevant Microbiota Changes. J Infect Dis. 2019;220: 1655–1666. doi:10.1093/INFDIS/JIZ299
    OpenUrlCrossRefPubMed
  53. 53.↵
    Bwire G, Waniaye JB, Otim JS, Matseketse D, Kagirita A, Orach CG. Cholera risk in cities in Uganda: understanding cases and contacts centered strategy (3CS) for rapid cholera outbreak control. Pan Afr Med J. 2021;39. doi:10.11604/PAMJ.2021.39.193.27794
    OpenUrlCrossRef
  54. 54.
    Ouamba JP, Mbarga NF, Ciglenecki I, Ratnayake R, Tchiasso D, Finger F, et al. Implementation of targeted cholera response activities, Cameroon. Bull World Health Organ. 2023;101: 170–178. doi:10.2471/BLT.22.288885
    OpenUrlCrossRefPubMed
  55. 55.↵
    Ratnayake R, Peyraud N, Ciglenecki I, Gignoux E, Lightowler M, Azman AS, et al. Effectiveness of case-area targeted interventions including vaccination on the control of epidemic cholera: protocol for a prospective observational study. BMJ Open. 2022;12. doi:10.1136/BMJOPEN-2022-061206
    OpenUrlCrossRef
  56. 56.↵
    Hegde ST, Khan AI, Perez-Saez J, Khan II, Hulse JD, Islam MT, et al. Clinical surveillance systems obscure the true cholera infection burden in an endemic region. Nature Medicine 2024 30:3. 2024;30: 888–895. doi:10.1038/s41591-024-02810-4
    OpenUrlCrossRef
  57. 57.↵
    Finger F, Lemaitre J, Juin S, Jackson B, Funk S, Lessler J, et al. Inferring the proportion of undetected cholera infections from serological and clinical surveillance in an immunologically naive population. medRxiv. 2023 [cited 20 Oct 2024]. doi:10.1101/2023.11.01.23297461
    OpenUrlAbstract/FREE Full Text
Back to top
Posted November 05, 2024.
Download PDF
Data/Code
Email
Quantifying the population-level impact of expanded antibiotic treatment for cholera outbreak management
Sharia M. Ahmed, Cormac R. LaPrete, Iza Ciglenecki, Andrew S. Azman, Daniel T. Leung, Lindsay T. Keegan
medRxiv 2024.11.04.24316579; doi: https://doi.org/10.1101/2024.11.04.24316579
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area