Incubation periods impact the spatial predictability of cholera and Ebola outbreaks in Sierra Leone

EpiFusion: Joint inference of the effective reproduction number by integrating phylodynamic and epidemiological modelling with particle filtering

Accurately estimating the effective reproduction number (Rt) of a circulating pathogen is a fundamental challenge in the study of infectious disease. The fields of epidemiology and pathogen phylodynamics both share this goal, but to date, methodologies and data employed by each remain largely distinct. Here we present EpiFusion: a joint approach that can be used to harness the complementary strengths of each field to improve estimation of outbreak dynamics for large and poorly sampled epidemics, such as arboviral or respiratory virus outbreaks, and validate it for retrospective analysis. We propose a model of Rt that estimates outbreak trajectories conditional upon both phylodynamic (time-scaled trees estimated from genetic sequences) and epidemiological (case incidence) data. We simulate stochastic outbreak trajectories that are weighted according to epidemiological and phylodynamic observation models and fit using particle Markov Chain Monte Carlo. To assess performance, we test EpiFusion on simulated outbreaks in which transmission and/or surveillance rapidly changes and find that using EpiFusion to combine epidemiological and phylodynamic data maintains accuracy and increases certainty in trajectory and Rt estimates, compared to when each data type is used alone. We benchmark EpiFusion's performance against existing methods to estimate Rt and demonstrate advances in speed and accuracy. Importantly, our approach scales efficiently with dataset size. Finally, we apply our model to estimate Rt during the 2014 Ebola outbreak in Sierra Leone. EpiFusion is designed to accommodate future extensions that will improve its utility, such as explicitly modelling population structure, accommodations for phylogenetic uncertainty, and the ability to weight the contributions of genomic or case incidence to the inference.

Prevention and control of Ebola virus transmission: mathematical modelling and data fitting

The Ebola virus disease (EVD) has been endemic since 1976, and the case fatality rate is extremely high. EVD is spread by infected animals, symptomatic individuals, dead bodies, and contaminated environment. In this paper, we formulate an EVD model with four transmission modes and a time delay describing the incubation period. Through dynamical analysis, we verify the importance of blocking the infection source of infected animals. We get the basic reproduction number without considering the infection source of infected animals. And, it is proven that the model has a globally attractive disease-free equilibrium when the basic reproduction number is less than unity; the disease eventually becomes endemic when the basic reproduction number is greater than unity. Taking the EVD epidemic in Sierra Leone in 2014-2016 as an example, we complete the data fitting by combining the effect of the media to obtain the unknown parameters, the basic reproduction number and its time-varying reproduction number. It is shown by parameter sensitivity analysis that the contact rate and the removal rate of infected group have the greatest influence on the prevalence of the disease. And, the disease-controlling thresholds of these two parameters are obtained. In addition, according to the existing vaccination strategy, only the inoculation ratio in high-risk areas is greater than 0.4, the effective reproduction number can be less than unity. And, the earlier the vaccination time, the greater the inoculation ratio, and the faster the disease can be controlled.

Key epidemiological indicators and spatial autocorrelation patterns across five waves of COVID-19 in Catalonia

This research studies the evolution of COVID-19 crude incident rates, effective reproduction number R(t) and their relationship with incidence spatial autocorrelation patterns in the 19 months following the disease outbreak in Catalonia (Spain). A cross-sectional ecological panel design based on n = 371 health-care geographical units is used. Five general outbreaks are described, systematically preceded by generalized values of R(t) > 1 in the two previous weeks. No clear regularities concerning possible initial focus appear when comparing waves. As for autocorrelation, we identify a wave's baseline pattern in which global Moran's I increases rapidly in the first weeks of the outbreak to descend later. However, some waves significantly depart from the baseline. In the simulations, both baseline pattern and departures can be reproduced when measures aimed at reducing mobility and virus transmissibility are introduced. Spatial autocorrelation is inherently contingent on the outbreak phase and is also substantially modified by external interventions affecting human behavior.

Cholera past and future in Nigeria: Are the Global Task Force on Cholera Control's 2030 targets achievable?

BACKGROUND

Understanding and continually assessing the achievability of global health targets is key to reducing disease burden and mortality. The Global Task Force on Cholera Control (GTFCC) Roadmap aims to reduce cholera deaths by 90% and eliminate the disease in twenty countries by 2030. The Roadmap has three axes focusing on reporting, response and coordination. Here, we assess the achievability of the GTFCC targets in Nigeria and identify where the three axes could be strengthened to reach and exceed these goals.

METHODOLOGY/PRINCIPAL FINDINGS

Using cholera surveillance data from Nigeria, cholera incidence was calculated and used to model time-varying reproduction number (R). A best fit random forest model was identified using R as the outcome variable and several environmental and social covariates were considered in the model, using random forest variable importance and correlation clustering. Future scenarios were created (based on varying degrees of socioeconomic development and emission reductions) and used to project future cholera transmission, nationally and sub-nationally to 2070. The projections suggest that significant reductions in cholera cases could be achieved by 2030, particularly in the more developed southern states, but increases in cases remain a possibility. Meeting the 2030 target, nationally, currently looks unlikely and we propose a new 2050 target focusing on reducing regional inequities, while still advocating for cholera elimination being achieved as soon as possible.

CONCLUSION/SIGNIFICANCE

The 2030 targets could potentially be reached by 2030 in some parts of Nigeria, but more effort is needed to reach these targets at a national level, particularly through access and incentives to cholera testing, sanitation expansion, poverty alleviation and urban planning. The results highlight the importance of and how modelling studies can be used to inform cholera policy and the potential for this to be applied in other contexts.

Tobler's law and wavefront patterns in the spatial spread of COVID-19 across Europe during the Delta and Omicron waves

AIMS

Epidemic wavefront models predict the spread of medieval pandemics such as the plague well. Our aim was to explore whether they contribute to understanding the spread of COVID-19, the first truly global pandemic of the 21st century with its fast and frequent international travel links.

METHODS

We analysed the spatial spread of reaching a threshold of very high incidence of new daily infections of the virus across European countries in the autumn of 2021 in which the Delta variant was dominant, as well as an even higher threshold of incidence in the subsequent spread of infections across the same set of countries during the winter of 2021/2022 when the Omicron variant of the virus became dominant.

RESULTS

We found patterns that are consistent with wavefront models for both periods of the pandemic in Europe.

CONCLUSIONS

Modern means of transportation strongly accelerated the spread of the virus and typically generated diffusion patterns along bidirectional constrained mobility networks in addition to stochastic diffusion processes. However, since the majority of mobility, including mobility across international borders, is over short distances, wavefront patterns in the spread of a pandemic are still to be expected.

Mathematical Models for Cholera Dynamics-A Review

Cholera remains a significant public health burden in many countries and regions of the world, highlighting the need for a deeper understanding of the mechanisms associated with its transmission, spread, and control. Mathematical modeling offers a valuable research tool to investigate cholera dynamics and explore effective intervention strategies. In this article, we provide a review of the current state in the modeling studies of cholera. Starting from an introduction of basic cholera transmission models and their applications, we survey model extensions in several directions that include spatial and temporal heterogeneities, effects of disease control, impacts of human behavior, and multi-scale infection dynamics. We discuss some challenges and opportunities for future modeling efforts on cholera dynamics, and emphasize the importance of collaborations between different modeling groups and different disciplines in advancing this research area.

Connecting the dots: understanding how human mobility shapes TB epidemics

Tuberculosis (TB) remains a leading infectious cause of death worldwide. Reducing TB infections and TB-related deaths rests ultimately on stopping forward transmission from infectious to susceptible individuals. Critical to this effort is understanding how human host mobility shapes the transmission and dispersal of new or existing strains of Mycobacterium tuberculosis (Mtb). Important questions remain unanswered. What kinds of mobility, over what temporal and spatial scales, facilitate TB transmission? How do human mobility patterns influence the dispersal of novel Mtb strains, including emergent drug-resistant strains? This review summarizes the current state of knowledge on mobility and TB epidemic dynamics, using examples from three topic areas, including inference of genetic and spatial clustering of infections, delineating source-sink dynamics, and mapping the dispersal of novel TB strains, to examine scientific questions and methodological issues within this topic. We also review new data sources for measuring human mobility, including mobile phone-associated movement data, and discuss important limitations on their use in TB epidemiology.

Fundamental limits on inferring epidemic resurgence in real time using effective reproduction numbers

We find that epidemic resurgence, defined as an upswing in the effective reproduction number (R) of the contagion from subcritical to supercritical values, is fundamentally difficult to detect in real time. Inherent latencies in pathogen transmission, coupled with smaller and intrinsically noisier case incidence across periods of subcritical spread, mean that resurgence cannot be reliably detected without significant delays of the order of the generation time of the disease, even when case reporting is perfect. In contrast, epidemic suppression (where R falls from supercritical to subcritical values) may be ascertained 5-10 times faster due to the naturally larger incidence at which control actions are generally applied. We prove that these innate limits on detecting resurgence only worsen when spatial or demographic heterogeneities are incorporated. Consequently, we argue that resurgence is more effectively handled proactively, potentially at the expense of false alarms. Timely responses to recrudescent infections or emerging variants of concern are more likely to be possible when policy is informed by a greater quality and diversity of surveillance data than by further optimisation of the statistical models used to process routine outbreak data.

Optimal strategies to protect a sub-population at risk due to an established epidemic

Epidemics can particularly threaten certain sub-populations. For example, for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the elderly are often preferentially protected. For diseases of plants and animals, certain sub-populations can drive mitigation because they are intrinsically more valuable for ecological, economic, socio-cultural or political reasons. Here, we use optimal control theory to identify strategies to optimally protect a ‘high-value’ sub-population when there is a limited budget and epidemiological uncertainty. We use protection of the Redwood National Park in California in the face of the large ongoing state-wide epidemic of sudden oak death (caused by Phytophthora ramorum) as a case study. We concentrate on whether control should be focused entirely within the National Park itself, or whether treatment of the growing epidemic in the surrounding ‘buffer region’ can instead be more profitable. We find that, depending on rates of infection and the size of the ongoing epidemic, focusing control on the high-value region is often optimal. However, priority should sometimes switch from the buffer region to the high-value region only as the local outbreak grows. We characterize how the timing of any switch depends on epidemiological and logistic parameters, and test robustness to systematic misspecification of these factors due to imperfect prior knowledge.

The impact of social and environmental extremes on cholera time varying reproduction number in Nigeria

Nigeria currently reports the second highest number of cholera cases in Africa, with numerous socioeconomic and environmental risk factors. Less investigated are the role of extreme events, despite recent work showing their potential importance. To address this gap, we used a machine learning approach to understand the risks and thresholds for cholera outbreaks and extreme events, taking into consideration pre-existing vulnerabilities. We estimated time varying reproductive number (R) from cholera incidence in Nigeria and used a machine learning approach to evaluate its association with extreme events (conflict, flood, drought) and pre-existing vulnerabilities (poverty, sanitation, healthcare). We then created a traffic-light system for cholera outbreak risk, using three hypothetical traffic-light scenarios (Red, Amber and Green) and used this to predict R. The system highlighted potential extreme events and socioeconomic thresholds for outbreaks to occur. We found that reducing poverty and increasing access to sanitation lessened vulnerability to increased cholera risk caused by extreme events (monthly conflicts and the Palmers Drought Severity Index). The main limitation is the underreporting of cholera globally and the potential number of cholera cases missed in the data used here. Increasing access to sanitation and decreasing poverty reduced the impact of extreme events in terms of cholera outbreak risk. The results here therefore add further evidence of the need for sustainable development for disaster prevention and mitigation and to improve health and quality of life.

Trip duration drives shift in travel network structure with implications for the predictability of spatial disease spread

Human travel is one of the primary drivers of infectious disease spread. Models of travel are often used that assume the amount of travel to a specific destination decreases as cost of travel increases with higher travel volumes to more populated destinations. Trip duration, the length of time spent in a destination, can also impact travel patterns. We investigated the spatial patterns of travel conditioned on trip duration and find distinct differences between short and long duration trips. In short-trip duration travel networks, trips are skewed towards urban destinations, compared with long-trip duration networks where travel is more evenly spread among locations. Using gravity models to inform connectivity patterns in simulations of disease transmission, we show that pathogens with shorter generation times exhibit initial patterns of spatial propagation that are more predictable among urban locations. Further, pathogens with a longer generation time have more diffusive patterns of spatial spread reflecting more unpredictable disease dynamics.

During an epidemic of an infectious pathogen, cases of disease can be imported to new locations when people travel. The amount of time that an infected person spends in a destination (trip duration) determines how likely they are to infect others while travelling. In this study, we analyzed travel data and found specific spatial patterns in trip duration, where short-duration trips are more common between urban destinations and long-duration trips are evenly spread out among locations. To show how this spatial pattern impacts the spread of infectious diseases, we used data-driven models and simulations to show that pathogens with shorter generation times have patterns of spatial spread that are more predictable among urban locations. However, pathogens with longer generation times tend to spread along the long-duration travel networks that are more evenly distributed among locations giving them more unpredictable disease dynamics.

Discrete simulation analysis of COVID-19 and prediction of isolation bed numbers

Background

The outbreak of COVID-19 has been defined by the World Health Organization as a pandemic, and containment depends on traditional public health measures. However, the explosive growth of the number of infected cases in a short period of time has caused tremendous pressure on medical systems. Adequate isolation facilities are essential to control outbreaks, so this study aims to quickly estimate the demand and number of isolation beds.

Methods

We established a discrete simulation model for epidemiology. By adjusting or fitting necessary epidemic parameters, the effects of the following indicators on the development of the epidemic and the occupation of medical resources were explained: (1) incubation period, (2) response speed and detection capacity of the hospital, (3) disease healing time, and (4) population mobility. Finally, a method for predicting the number of isolation beds was summarized through multiple linear regression. This is a city level model that simulates the epidemic situation from the perspective of population mobility.

Results

Through simulation, we show that the incubation period, response speed and detection capacity of the hospital, disease healing time, degree of population mobility, and infectivity of cured patients have different effects on the infectivity, scale, and duration of the epidemic. Among them, (1) incubation period, (2) response speed and detection capacity of the hospital, (3) disease healing time, and (4) population mobility have a significant impact on the demand and number of isolation beds (P <0.05), which agrees with the following regression equation: N = P × (−0.273 + 0.009I + 0.234M + 0.012T1 + 0.015T2) × (1 + V).

How lessons learned from the 2015 Middle East respiratory syndrome outbreak affected the response to coronavirus disease 2019 in the Republic of Korea

The Republic of Korea (ROK) experienced a public health crisis due to Middle East respiratory syndrome (MERS) in 2015 and is currently going through the coronavirus disease 2019 (COVID-19) pandemic. Lessons learned from the disastrous MERS outbreak were ref lected in the preparedness system, and the readiness capabilities that were subsequently developed enabled the country to successfully flatten the epidemic curve of COVID-19 in late February and March 2020. In this review, we summarize and compare the epidemiology and response of the ROK to the 2015 MERS outbreak and the COVID-19 epidemic in early 2020. We emphasize that, because further COVID-19 waves seem inevitable, it is urgent to develop comprehensive preparedness and response plans for the worst-case scenarios of the COVID-19 pandemic. Simultaneously strengthening healthcare capacity to endure the peak demand and implementing smart strategies to sustain social distancing and public hygiene are necessary until safe and effective therapeutics and vaccines against COVID-19 are available.

COVID-19 and the 1918 influenza pandemics: a concise overview and lessons from the past

The coronavirus disease 2019 (COVID-19), which is caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), was first reported in December, 2019, in Wuhan, China. Even the public health sector experts could not anticipate that the virus would spread rapidly to create the worst worldwide crisis in more than a century. The World Health Organization (WHO) declared COVID-19 a public health emergency on January 30, 2020, but it was not until March 11, 2020 that the WHO declared it a global pandemic. The epidemiology of SARS-CoV-2 is different from the SARS coronavirus outbreak in 2002 and the Middle East Respiratory Syndrome (MERS) in 2012; therefore, neither SARS nor MERS could be used as a suitable model for foreseeing the future of the current pandemic. The influenza pandemic of 1918 could be referred to in order to understand and control the COVID-19 pandemic. Although influenza and the SARS-CoV-2 are from different families of viruses, they are similar in that both silently attacked the world and the societal and political responses to both pandemics have been very much alike. Previously, the 1918 influenza pandemic and unpredictability of the second wave caused distress among people as the first wave of that outbreak (so-called Spanish flu) proved to be relatively mild compared to a much worse second wave, followed by smaller waves. As of April, 2021, the second wave of COVID-19 has occurred around the globe, and future waves may also be expected, if the total population of the world is not vaccinated. This article aims to highlight the key similarities and differences in both pandemics. Similarly, lessons from the previous pan-demics and various possibilities for the future course of COVID-19 are also highlighted.

Risk of Secondary Infection Waves of COVID-19 in an Insular Region: The Case of the Balearic Islands, Spain

The Spanish government declared the lockdown on March 14th, 2020 to tackle the fast-spreading of COVID-19. As a consequence, the Balearic Islands remained almost fully isolated due to the closing of airports and ports, these isolation measures and the home-based confinement have led to a low prevalence of COVID-19 in this region. We propose a compartmental model for the spread of COVID-19 including five compartments (Susceptible, Exposed, Presymptomatic Infective, Diseased, and Recovered), and the mobility between municipalities. The model parameters are calibrated with the temporal series of confirmed cases provided by the Spanish Ministry of Health. After calibration, the proposed model captures the trend of the official confirmed cases before and after the lockdown. We show that the estimated number of cases depends strongly on the initial dates of the local outbreak onset and the number of imported cases before the lockdown. Our estimations indicate that the population has not reached the level of herd immunization necessary to prevent future outbreaks. While the low prevalence, in comparison to mainland Spain, has prevented the saturation of the health system, this low prevalence translates into low immunization rates, therefore facilitating the propagation of new outbreaks that could lead to secondary waves of COVID-19 in the region. These findings warn about scenarios regarding after-lockdown-policies and the risk of second outbreaks, emphasize the need for widespread testing, and could potentially be extrapolated to other insular and continental regions.

Early detection of cholera epidemics to support control in fragile states: estimation of delays and potential epidemic sizes

BACKGROUND

Cholera epidemics continue to challenge disease control, particularly in fragile and conflict-affected states. Rapid detection and response to small cholera clusters is key for efficient control before an epidemic propagates. To understand the capacity for early response in fragile states, we investigated delays in outbreak detection, investigation, response, and laboratory confirmation, and we estimated epidemic sizes. We assessed predictors of delays, and annual changes in response time.

METHODS

We compiled a list of cholera outbreaks in fragile and conflict-affected states from 2008 to 2019. We searched for peer-reviewed articles and epidemiological reports. We evaluated delays from the dates of symptom onset of the primary case, and the earliest dates of outbreak detection, investigation, response, and confirmation. Information on how the outbreak was alerted was summarized. A branching process model was used to estimate epidemic size at each delay. Regression models were used to investigate the association between predictors and delays to response.

RESULTS

Seventy-six outbreaks from 34 countries were included. Median delays spanned 1-2 weeks: from symptom onset of the primary case to presentation at the health facility (5 days, IQR 5-5), detection (5 days, IQR 5-6), investigation (7 days, IQR 5.8-13.3), response (10 days, IQR 7-18), and confirmation (11 days, IQR 7-16). In the model simulation, the median delay to response (10 days) with 3 seed cases led to a median epidemic size of 12 cases (upper range, 47) and 8% of outbreaks ≥ 20 cases (increasing to 32% with a 30-day delay to response). Increased outbreak size at detection (10 seed cases) and a 10-day median delay to response resulted in an epidemic size of 34 cases (upper range 67 cases) and < 1% of outbreaks < 20 cases. We estimated an annual global decrease in delay to response of 5.2% (95% CI 0.5-9.6, p = 0.03). Outbreaks signaled by immediate alerts were associated with a reduction in delay to response of 39.3% (95% CI 5.7-61.0, p = 0.03).

CONCLUSIONS

From 2008 to 2019, median delays from symptom onset of the primary case to case presentation and to response were 5 days and 10 days, respectively. Our model simulations suggest that depending on the outbreak size (3 versus 10 seed cases), in 8 to 99% of scenarios, a 10-day delay to response would result in large clusters that would be difficult to contain. Improving the delay to response involves rethinking the integration at local levels of event-based detection, rapid diagnostic testing for cluster validation, and integrated alert, investigation, and response.