Temporal variations in international air travel: implications for modelling the spread of infectious diseases

BACKGROUND

The international flight network creates multiple routes by which pathogens can quickly spread across the globe. In the early stages of infectious disease outbreaks, analyses using flight passenger data to identify countries at risk of importing the pathogen are common and can help inform disease control efforts. A challenge faced in this modelling is that the latest aviation statistics (referred to as contemporary data) are typically not immediately available. Therefore, flight patterns from a previous year are often used (referred to as historical data). We explored the suitability of historical data for predicting the spatial spread of emerging epidemics.

METHODS

We analysed monthly flight passenger data from the International Air Transport Association to assess how baseline air travel patterns were affected in outbreaks of MERS, Zika, and SARS-CoV-2 over the past decade. We then used a stochastic discrete time SEIR metapopulation model to simulate global spread of different pathogens, comparing how epidemic dynamics differed in simulations based on historical and contemporary data.

RESULTS

We observed local, short-term disruptions to air travel from South Korea and Brazil for the MERS and Zika outbreaks we studied, whereas global and longer-term flight disruption occurred during the SARS-CoV-2 pandemic.For outbreak events that were accompanied by local, small, and short-term changes in air travel, epidemic models using historical flight data gave similar projections of timing and locations of disease spread as when using contemporary flight data. However, historical data were less reliable to model the spread of an atypical outbreak such as SARS-CoV-2 in which there were durable and extensive levels of global travel disruption.

CONCLUSIONS

The use of historical flight data as a proxy in epidemic models is an acceptable practice except in rare, large epidemics that lead to substantial disruptions to international travel.

Lessons for cross-species viral transmission surveillance from highly pathogenic avian influenza Korean cat shelter outbreaks

In this Comment, the authors describe recent outbreaks of highly pathogenic avian influenza in cat shelters in Seoul, South Korea. They discuss potential routes of transmission and describe implications for surveillance of spillover infections in animals in non-agricultural settings.

Retrospective evaluation of real-time estimates of global COVID-19 transmission trends and mortality forecasts

Since 8th March 2020 up to the time of writing, we have been producing near real-time weekly estimates of SARS-CoV-2 transmissibility and forecasts of deaths due to COVID-19 for all countries with evidence of sustained transmission, shared online. We also developed a novel heuristic to combine weekly estimates of transmissibility to produce forecasts over a 4-week horizon. Here we present a retrospective evaluation of the forecasts produced between 8th March to 29th November 2020 for 81 countries. We evaluated the robustness of the forecasts produced in real-time using relative error, coverage probability, and comparisons with null models. During the 39-week period covered by this study, both the short- and medium-term forecasts captured well the epidemic trajectory across different waves of COVID-19 infections with small relative errors over the forecast horizon. The model was well calibrated with 56.3% and 45.6% of the observations lying in the 50% Credible Interval in 1-week and 4-week ahead forecasts respectively. The retrospective evaluation of our models shows that simple transmission models calibrated using routine disease surveillance data can reliably capture the epidemic trajectory in multiple countries. The medium-term forecasts can be used in conjunction with the short-term forecasts of COVID-19 mortality as a useful planning tool as countries continue to relax public health measures.

From serological surveys to disease burden: a modelling pipeline for Chagas disease

In 2012, the World Health Organization (WHO) set the elimination of Chagas disease intradomiciliary vectorial transmission as a goal by 2020. After a decade, some progress has been made, but the new 2021-2030 WHO roadmap has set even more ambitious targets. Innovative and robust modelling methods are required to monitor progress towards these goals. We present a modelling pipeline using local seroprevalence data to obtain national disease burden estimates by disease stage. Firstly, local seroprevalence information is used to estimate spatio-temporal trends in the Force-of-Infection (FoI). FoI estimates are then used to predict such trends across larger and fine-scale geographical areas. Finally, predicted FoI values are used to estimate disease burden based on a disease progression model. Using Colombia as a case study, we estimated that the number of infected people would reach 506 000 (95% credible interval (CrI) = 395 000-648 000) in 2020 with a 1.0% (95%CrI = 0.8-1.3%) prevalence in the general population and 2400 (95%CrI = 1900-3400) deaths (approx. 0.5% of those infected). The interplay between a decrease in infection exposure (FoI and relative proportion of acute cases) was overcompensated by a large increase in population size and gradual population ageing, leading to an increase in the absolute number of Chagas disease cases over time. This article is part of the theme issue 'Challenges and opportunities in the fight against neglected tropical diseases: a decade from the London Declaration on NTDs'.

Lessons from COVID-19 for rescalable data collection

Novel data and analyses have had an important role in informing the public health response to the COVID-19 pandemic. Existing surveillance systems were scaled up, and in some instances new systems were developed to meet the challenges posed by the magnitude of the pandemic. We describe the routine and novel data that were used to address urgent public health questions during the pandemic, underscore the challenges in sustainability and equity in data generation, and highlight key lessons learnt for designing scalable data collection systems to support decision making during a public health crisis. As countries emerge from the acute phase of the pandemic, COVID-19 surveillance systems are being scaled down. However, SARS-CoV-2 resurgence remains a threat to global health security; therefore, a minimal cost-effective system needs to remain active that can be rapidly scaled up if necessary. We propose that a retrospective evaluation to identify the cost-benefit profile of the various data streams collected during the pandemic should be on the scientific research agenda.

Contrasting seasonality of African swine fever outbreaks and its drivers

The seasonality of African swine fever (ASF) outbreaks in domestic pigs differs between temperate and subtropical/tropical regions. We hypothesise that variations in the importance of wild boar-to-farm and farm-to-farm transmission routes shape these contrasting patterns, and we emphasise the implications for effective ASF control.

Extending EpiEstim to estimate the transmission advantage of pathogen variants in real-time: SARS-CoV-2 as a case-study

The evolution of SARS-CoV-2 has demonstrated that emerging variants can set back the global COVID-19 response. The ability to rapidly assess the threat of new variants is critical for timely optimisation of control strategies. We present a novel method to estimate the effective transmission advantage of a new variant compared to a reference variant combining information across multiple locations and over time. Through an extensive simulation study designed to mimic real-time epidemic contexts, we show that our method performs well across a range of scenarios and provide guidance on its optimal use and interpretation of results. We also provide an open-source software implementation of our method. The computational speed of our tool enables users to rapidly explore spatial and temporal variations in the estimated transmission advantage. We estimate that the SARS-CoV-2 Alpha variant is 1.46 (95% Credible Interval 1.44-1.47) and 1.29 (95% CrI 1.29-1.30) times more transmissible than the wild type, using data from England and France respectively. We further estimate that Delta is 1.77 (95% CrI 1.69-1.85) times more transmissible than Alpha (England data). Our approach can be used as an important first step towards quantifying the threat of emerging or co-circulating variants of infectious pathogens in real-time.

Estimating the epidemic reproduction number from temporally aggregated incidence data: A statistical modelling approach and software tool

The time-varying reproduction number (Rt) is an important measure of epidemic transmissibility that directly informs policy decisions and the optimisation of control measures. EpiEstim is a widely used opensource software tool that uses case incidence and the serial interval (SI, time between symptoms in a case and their infector) to estimate Rt in real-time. The incidence and the SI distribution must be provided at the same temporal resolution, which can limit the applicability of EpiEstim and other similar methods, e.g. for contexts where the time window of incidence reporting is longer than the mean SI. In the EpiEstim R package, we implement an expectation-maximisation algorithm to reconstruct daily incidence from temporally aggregated data, from which Rt can then be estimated. We assess the validity of our method using an extensive simulation study and apply it to COVID-19 and influenza data. For all datasets, the influence of intra-weekly variability in reported data was mitigated by using aggregated weekly data. Rt estimated on weekly sliding windows using incidence reconstructed from weekly data was strongly correlated with estimates from the original daily data. The simulation study revealed that Rt was well estimated in all scenarios and regardless of the temporal aggregation of the data. In the presence of weekend effects, Rt estimates from reconstructed data were more successful at recovering the true value of Rt than those obtained from reported daily data. These results show that this novel method allows Rt to be successfully recovered from aggregated data using a simple approach with very few data requirements. Additionally, by removing administrative noise when daily incidence data are reconstructed, the accuracy of Rt estimates can be improved.

Plastic Entanglement Poses a Potential Hazard to European Hedgehogs Erinaceus europaeus in Great Britain

A questionnaire to gather evidence on the plastic entanglement of the European hedgehog (Erinaceus europaeus) was sent to 160 wildlife rehabilitation centres in Great Britain. Fifty-four responses were received, and 184 individual admissions owing to plastic entanglement were reported. Death was the outcome for 46% (n = 86) of these cases. A high proportion of Britain's hedgehogs enter rehabilitation centres annually (approximately 5% of the national population and potentially 10% of the urban population), providing a robust basis for assessing the minimum impacts at a national level. We estimate that 4000-7000 hedgehog deaths per year are attributable to plastic, with the true rate likely being higher, since many entangled hedgehogs-in contrast to those involved in road traffic accidents-will not be found. Population modelling indicates that this excess mortality is sufficient to cause population declines. Although the scale of the impact is much lower than that attributable to traffic, it is nevertheless an additional pressure on a species that is already in decline and presents a significant welfare issue to a large number of individuals.

Transmission dynamics of lyssavirus in Myotis myotis: mechanistic modelling study based on longitudinal seroprevalence data

We investigated the transmission dynamics of lyssavirus in Myotis myotis and Myotis blythii, using serological, virological, demographic and ecological data collected between 2015 and 2022 from two maternity colonies in northern Italian churches. Despite no lyssavirus detection in 556 bats sampled over 11 events by reverse transcription-polymerase chain reaction (RT-PCR), 36.3% of 837 bats sampled over 27 events showed neutralizing antibodies to European bat lyssavirus 1, with a significant increase in summers. By fitting sets of mechanistic models to seroprevalence data, we investigated factors that influenced lyssavirus transmission within and between years. Five models were selected as a group of final models: in one model, a proportion of exposed bats (median model estimate: 5.8%) became infectious and died while the other exposed bats recovered with immunity without becoming infectious; in the other four models, all exposed bats became infectious and recovered with immunity. The final models supported that the two colonies experienced seasonal outbreaks driven by: (i) immunity loss particularly during hibernation, (ii) density-dependent transmission, and (iii) a high transmission rate after synchronous birthing. These findings highlight the importance of understanding ecological factors, including colony size and synchronous birthing timing, and potential infection heterogeneities to enable more robust assessments of lyssavirus spillover risk.

Drivers of SARS-CoV-2 testing behaviour: a modelling study using nationwide testing data in England

During the COVID-19 pandemic, national testing programmes were conducted worldwide on unprecedented scales. While testing behaviour is generally recognised as dynamic and complex, current literature demonstrating and quantifying such relationships is scarce, despite its importance for infectious disease surveillance and control. Here, we characterise the impacts of SARS-CoV-2 transmission, disease susceptibility/severity, risk perception, and public health measures on SARS-CoV-2 PCR testing behaviour in England over 20 months of the pandemic, by linking testing trends to underlying epidemic trends and contextual meta-data within a systematic conceptual framework. The best-fitting model describing SARS-CoV-2 PCR testing behaviour explained close to 80% of the total deviance in NHS test data. Testing behaviour showed complex associations with factors reflecting transmission level, disease susceptibility/severity (e.g. age, dominant variant, and vaccination), public health measures (e.g. testing strategies and lockdown), and associated changes in risk perception, varying throughout the pandemic and differing between infected and non-infected people.

Gaps in mobility data and implications for modelling epidemic spread: A scoping review and simulation study

Reliable estimates of human mobility are important for understanding the spatial spread of infectious diseases and the effective targeting of control measures. However, when modelling infectious disease dynamics, data on human mobility at an appropriate temporal or spatial resolution are not always available, leading to the common use of model-derived mobility proxies. In this study we reviewed the different data sources and mobility models that have been used to characterise human movement in Africa. We then conducted a simulation study to better understand the implications of using human mobility proxies when predicting the spatial spread and dynamics of infectious diseases. We found major gaps in the availability of empirical measures of human mobility in Africa, leading to mobility proxies being used in place of data. Empirical data on subnational mobility were only available for 17/54 countries, and in most instances, these data characterised long-term movement patterns, which were unsuitable for modelling the spread of pathogens with short generation times (time between infection of a case and their infector). Results from our simulation study demonstrated that using mobility proxies can have a substantial impact on the predicted epidemic dynamics, with complex and non-intuitive biases. In particular, the predicted times and order of epidemic invasion, and the time of epidemic peak in different locations can be underestimated or overestimated, depending on the types of proxies used and the country of interest. Our work underscores the need for regularly updated empirical measures of population movement within and between countries to aid the prevention and control of infectious disease outbreaks. At the same time, there is a need to establish an evidence base to help understand which types of mobility data are most appropriate for describing the spread of emerging infectious diseases in different settings.

Estimating Zika virus attack rates and risk of Zika virus-associated neurological complications in Colombian capital cities with a Bayesian model

Zika virus (ZIKV) is a mosquito-borne pathogen that caused a major epidemic in the Americas in 2015-2017. Although the majority of ZIKV infections are asymptomatic, the virus has been associated with congenital birth defects and neurological complications (NC) in adults. We combined multiple data sources to improve estimates of ZIKV infection attack rates (IARs), reporting rates of Zika virus disease (ZVD) and the risk of ZIKV-associated NC for 28 capital cities in Colombia. ZVD surveillance data were combined with post-epidemic seroprevalence data and a dataset on ZIKV-associated NC in a Bayesian hierarchical model. We found substantial heterogeneity in ZIKV IARs across cities. The overall estimated ZIKV IAR across the 28 cities was 0.38 (95% CrI: 0.17-0.92). The estimated ZVD reporting rate was 0.013 (95% CrI: 0.004-0.024), and 0.51 (95% CrI: 0.17-0.92) cases of ZIKV-associated NC were estimated to be reported per 10 000 ZIKV infections. When we assumed the same ZIKV IAR across sex or age group, we found important spatial heterogeneities in ZVD reporting rates and the risk of being reported as a ZVD case with NC. Our results highlight how additional data sources can be used to overcome biases in surveillance data and estimate key epidemiological parameters.

Real-time estimation of the epidemic reproduction number: Scoping review of the applications and challenges

The time-varying reproduction number (Rt) is an important measure of transmissibility during outbreaks. Estimating whether and how rapidly an outbreak is growing (Rt > 1) or declining (Rt < 1) can inform the design, monitoring and adjustment of control measures in real-time. We use a popular R package for Rt estimation, EpiEstim, as a case study to evaluate the contexts in which Rt estimation methods have been used and identify unmet needs which would enable broader applicability of these methods in real-time. A scoping review, complemented by a small EpiEstim user survey, highlight issues with the current approaches, including the quality of input incidence data, the inability to account for geographical factors, and other methodological issues. We summarise the methods and software developed to tackle the problems identified, but conclude that significant gaps remain which should be addressed to enable easier, more robust and applicable estimation of Rt during epidemics.

Spatiotemporal variations in exposure: Chagas disease in Colombia as a case study

Age-stratified serosurvey data are often used to understand spatiotemporal trends in disease incidence and exposure through estimating the Force-of-Infection (FoI). Typically, median or mean FoI estimates are used as the response variable in predictive models, often overlooking the uncertainty in estimated FoI values when fitting models and evaluating their predictive ability. To assess how this uncertainty impact predictions, we compared three approaches with three levels of uncertainty integration. We propose a performance indicator to assess how predictions reflect initial uncertainty.

In Colombia, 76 serosurveys (1980–2014) conducted at municipality level provided age-stratified Chagas disease prevalence data. The yearly FoI was estimated at the serosurvey level using a time-varying catalytic model. Environmental, demographic and entomological predictors were used to fit and predict the FoI at municipality level from 1980 to 2010 across Colombia.

A stratified bootstrap method was used to fit the models without temporal autocorrelation at the serosurvey level. The predictive ability of each model was evaluated to select the best-fit models within urban, rural and (Amerindian) indigenous settings. Model averaging, with the 10 best-fit models identified, was used to generate predictions.

Our analysis shows a risk of overconfidence in model predictions when median estimates of FoI alone are used to fit and evaluate models, failing to account for uncertainty in FoI estimates. Our proposed methodology fully propagates uncertainty in the estimated FoI onto the generated predictions, providing realistic assessments of both central tendency and current uncertainty surrounding exposure to Chagas disease.

Understanding the incidence and timing of rabies cases in domestic animals and wildlife in south-east Tanzania in the presence of widespread domestic dog vaccination campaigns

The "Zero by 30" strategic plan aims to eliminate human deaths from dog-mediated rabies by 2030 and domestic dog vaccination is a vital component of this strategic plan. In areas where domestic dog vaccination has been implemented, it is important to assess the impact of this intervention. Additionally, understanding temporal and seasonal trends in the incidence of animal rabies cases may assist in optimizing such interventions. Data on the incidence of probable rabies cases in domestic and wild animals were collected between January 2011 and December 2018 in thirteen districts of south-east Tanzania where jackals comprise over 40% of reported rabies cases. Vaccination coverage was estimated over this period, as five domestic dog vaccination campaigns took place in all thirteen districts between 2011 and 2016. Negative binomial generalized linear models were used to explore the impact of domestic dog vaccination on the annual incidence of animal rabies cases, whilst generalized additive models were used to investigate the presence of temporal and/or seasonal trends. Increases in domestic dog vaccination coverage were significantly associated with a decreased incidence of rabies cases in both domestic dogs and jackals. A 35% increase in vaccination coverage was associated with a reduction in the incidence of probable dog rabies cases of between 78.0 and 85.5% (95% confidence intervals ranged from 61.2 to 92.2%) and a reduction in the incidence of probable jackal rabies cases of between 75.3 and 91.2% (95% confidence intervals ranged from 53.0 to 96.1%). A statistically significant common seasonality was identified in the monthly incidence of probable rabies cases in both domestic dogs and jackals with the highest incidence from February to August and lowest incidence from September to January. These results align with evidence supporting the use of domestic dog vaccination as part of control strategies aimed at reducing animal rabies cases in both domestic dogs and jackals in this region. The presence of a common seasonal trend requires further investigation but may have implications for the timing of future vaccination campaigns.

Reservoir dynamics of rabies in south-east Tanzania and the roles of cross-species transmission and domestic dog vaccination

Understanding the role of different species in the transmission of multi-host pathogens, such as rabies virus, is vital for effective control strategies. Across most of sub-Saharan Africa domestic dogs Canis familiaris are considered the reservoir for rabies, but the role of wildlife has been long debated. Here we explore the multi-host transmission dynamics of rabies across south-east Tanzania.Between January 2011 and July 2019, data on probable rabies cases were collected in the regions of Lindi and Mtwara. Hospital records of animal-bite patients presenting to healthcare facilities were used as sentinels for animal contact tracing. The timing, location and species of probable rabid animals were used to reconstruct transmission trees to infer who infected whom and the relative frequencies of within- and between-species transmission.During the study, 688 probable human rabies exposures were identified, resulting in 47 deaths. Of these exposures, 389 were from domestic dogs (56.5%) and 262 from jackals (38.1%). Over the same period, 549 probable animal rabies cases were traced: 303 in domestic dogs (55.2%) and 221 in jackals (40.3%), with the remainder in domestic cats and other wildlife species.Although dog-to-dog transmission was most commonly inferred (40.5% of transmission events), a third of inferred events involved wildlife-to-wildlife transmission (32.6%), and evidence suggested some sustained transmission chains within jackal populations.A steady decline in probable rabies cases in both humans and animals coincided with the implementation of widespread domestic dog vaccination during the first 6 years of the study. Following the lapse of this program, dog rabies cases began to increase in one of the northernmost districts. Synthesis and applications. In south-east Tanzania, despite a relatively high incidence of rabies in wildlife and evidence of wildlife-to-wildlife transmission, domestic dogs remain essential to the reservoir of infection. Continued dog vaccination alongside improved surveillance would allow a fuller understanding of the role of wildlife in maintaining transmission in this area. Nonetheless, dog vaccination clearly suppressed rabies in both domestic dog and wildlife populations, reducing both public health and conservation risks and, if sustained, has potential to eliminate rabies from this region.

Modelling the influence of naturally acquired immunity from subclinical infection on outbreak dynamics and persistence of rabies in domestic dogs

A number of mathematical models have been developed for canine rabies to explore dynamics and inform control strategies. A common assumption of these models is that naturally acquired immunity plays no role in rabies dynamics. However, empirical studies have detected rabies-specific antibodies in healthy, unvaccinated domestic dogs, potentially due to immunizing, non-lethal exposure. We developed a stochastic model for canine rabies, parameterised for Laikipia County, Kenya, to explore the implications of different scenarios for naturally acquired immunity to rabies in domestic dogs. Simulating these scenarios using a non-spatial model indicated that low levels of immunity can act to limit rabies incidence and prevent depletion of the domestic dog population, increasing the probability of disease persistence. However, incorporating spatial structure and human response to high rabies incidence allowed the virus to persist in the absence of immunity. While low levels of immunity therefore had limited influence under a more realistic approximation of rabies dynamics, high rates of exposure leading to immunizing non-lethal exposure were required to produce population-level seroprevalences comparable with those reported in empirical studies. False positives and/or spatial variation may contribute to high empirical seroprevalences. However, if high seroprevalences are related to high exposure rates, these findings support the need for high vaccination coverage to effectively control this disease.

Using syndromic measures of mortality to capture the dynamics of COVID-19 in Java, Indonesia, in the context of vaccination rollout

BACKGROUND

As in many countries, quantifying COVID-19 spread in Indonesia remains challenging due to testing limitations. In Java, non-pharmaceutical interventions (NPIs) were implemented throughout 2020. However, as a vaccination campaign launches, cases and deaths are rising across the island.

METHODS

We used modelling to explore the extent to which data on burials in Jakarta using strict COVID-19 protocols (C19P) provide additional insight into the transmissibility of the disease, epidemic trajectory, and the impact of NPIs. We assess how implementation of NPIs in early 2021 will shape the epidemic during the period of likely vaccine rollout.

RESULTS

C19P burial data in Jakarta suggest a death toll approximately 3.3 times higher than reported. Transmission estimates using these data suggest earlier, larger, and more sustained impact of NPIs. Measures to reduce sub-national spread, particularly during Ramadan, substantially mitigated spread to more vulnerable rural areas. Given current trajectory, daily cases and deaths are likely to increase in most regions as the vaccine is rolled out. Transmission may peak in early 2021 in Jakarta if current levels of control are maintained. However, relaxation of control measures is likely to lead to a subsequent resurgence in the absence of an effective vaccination campaign.

CONCLUSIONS

Syndromic measures of mortality provide a more complete picture of COVID-19 severity upon which to base decision-making. The high potential impact of the vaccine in Java is attributable to reductions in transmission to date and dependent on these being maintained. Increases in control in the relatively short-term will likely yield large, synergistic increases in vaccine impact.

Descriptive analysis of surveillance data for Zika virus disease and Zika virus-associated neurological complications in Colombia, 2015-2017

Zika virus (ZIKV) is a mosquito-borne pathogen that recently caused a major epidemic in the Americas. Although the majority of ZIKV infections are asymptomatic, the virus has been associated with birth defects in fetuses and newborns of infected mothers as well as neurological complications in adults. We performed a descriptive analysis on approximately 106,000 suspected and laboratory-confirmed cases of Zika virus disease (ZVD) that were reported during the 2015-2017 epidemic in Colombia. We also analyzed a dataset containing patients with neurological complications and recent febrile illness compatible with ZVD. Females had higher cumulative incidence of ZVD than males. Compared to the general population, cases were more likely to be reported in young adults (20 to 39 years of age). We estimated the cumulative incidence of ZVD in pregnant females at 3,120 reported cases per 100,000 population (95% CI: 3,077-3,164), which was considerably higher than the incidence in both males and non-pregnant females. ZVD cases were reported in all 32 departments. Four-hundred and eighteen patients suffered from ZIKV-associated neurological complications, of which 85% were diagnosed with Guillain-Barré syndrome. The median age of ZIKV cases with neurological complications was 12 years older than that of ZVD cases. ZIKV-associated neurological complications increased with age, and the highest incidence was reported among individuals aged 75 and older. Even though neurological complications and deaths due to ZIKV were rare in this epidemic, better risk communication is needed for people living in or traveling to ZIKV-affected areas.

Using digital surveillance tools for near real-time mapping of the risk of infectious disease spread

Data from digital disease surveillance tools such as ProMED and HealthMap can complement the field surveillance during ongoing outbreaks. Our aim was to investigate the use of data collected through ProMED and HealthMap in real-time outbreak analysis. We developed a flexible statistical model to quantify spatial heterogeneity in the risk of spread of an outbreak and to forecast short term incidence trends. The model was applied retrospectively to data collected by ProMED and HealthMap during the 2013-2016 West African Ebola epidemic and for comparison, to WHO data. Using ProMED and HealthMap data, the model was able to robustly quantify the risk of disease spread 1-4 weeks in advance and for countries at risk of case importations, quantify where this risk comes from. Our study highlights that ProMED and HealthMap data could be used in real-time to quantify the spatial heterogeneity in risk of spread of an outbreak.

Reduction in mobility and COVID-19 transmission

In response to the COVID-19 pandemic, countries have sought to control SARS-CoV-2 transmission by restricting population movement through social distancing interventions, thus reducing the number of contacts. Mobility data represent an important proxy measure of social distancing, and here, we characterise the relationship between transmission and mobility for 52 countries around the world. Transmission significantly decreased with the initial reduction in mobility in 73% of the countries analysed, but we found evidence of decoupling of transmission and mobility following the relaxation of strict control measures for 80% of countries. For the majority of countries, mobility explained a substantial proportion of the variation in transmissibility (median adjusted R-squared: 48%, interquartile range - IQR - across countries [27-77%]). Where a change in the relationship occurred, predictive ability decreased after the relaxation; from a median adjusted R-squared of 74% (IQR across countries [49-91%]) pre-relaxation, to a median adjusted R-squared of 30% (IQR across countries [12-48%]) post-relaxation. In countries with a clear relationship between mobility and transmission both before and after strict control measures were relaxed, mobility was associated with lower transmission rates after control measures were relaxed indicating that the beneficial effects of ongoing social distancing behaviours were substantial.

Database of epidemic trends and control measures during the first wave of COVID-19 in mainland China

OBJECTIVES

In this data collation study, we aimed to provide a comprehensive database describing the epidemic trends and responses during the first wave of coronavirus disease 2019 (COVID-19) throughout the main provinces in China.

METHODS

From mid-January to March 2020, we extracted publicly available data regarding the spread and control of COVID-19 from 31 provincial health authorities and major media outlets in mainland China. Based on these data, we conducted descriptive analyses of the epidemic in the six most-affected provinces.

RESULTS

School closures, travel restrictions, community-level lockdown, and contact tracing were introduced concurrently around late January but subsequent epidemic trends differed among provinces. Compared with Hubei, the other five most-affected provinces reported a lower crude case fatality ratio and proportion of critical and severe hospitalised cases. From March 2020, as the local transmission of COVID-19 declined, switching the focus of measures to the testing and quarantine of inbound travellers may have helped to sustain the control of the epidemic.

CONCLUSIONS

Aggregated indicators of case notifications and severity distributions are essential for monitoring an epidemic. A publicly available database containing these indicators and information regarding control measures is a useful resource for further research and policy planning in response to the COVID-19 epidemic.

State-level tracking of COVID-19 in the United States

As of 1st June 2020, the US Centres for Disease Control and Prevention reported 104,232 confirmed or probable COVID-19-related deaths in the US. This was more than twice the number of deaths reported in the next most severely impacted country. We jointly model the US epidemic at the state-level, using publicly available death data within a Bayesian hierarchical semi-mechanistic framework. For each state, we estimate the number of individuals that have been infected, the number of individuals that are currently infectious and the time-varying reproduction number (the average number of secondary infections caused by an infected person). We use changes in mobility to capture the impact that non-pharmaceutical interventions and other behaviour changes have on the rate of transmission of SARS-CoV-2. We estimate that R_t was only below one in 23 states on 1st June. We also estimate that 3.7% [3.4%-4.0%] of the total population of the US had been infected, with wide variation between states, and approximately 0.01% of the population was infectious. We demonstrate good 3 week model forecasts of deaths with low error and good coverage of our credible intervals.

State-level tracking of COVID-19 in the United States

As of 1st June 2020, the US Centres for Disease Control and Prevention reported 104,232 confirmed or probable COVID-19-related deaths in the US. This was more than twice the number of deaths reported in the next most severely impacted country. We jointly model the US epidemic at the state-level, using publicly available death data within a Bayesian hierarchical semi-mechanistic framework. For each state, we estimate the number of individuals that have been infected, the number of individuals that are currently infectious and the time-varying reproduction number (the average number of secondary infections caused by an infected person). We use changes in mobility to capture the impact that non-pharmaceutical interventions and other behaviour changes have on the rate of transmission of SARS-CoV-2. We estimate that Rt was only below one in 23 states on 1st June. We also estimate that 3.7% [3.4%-4.0%] of the total population of the US had been infected, with wide variation between states, and approximately 0.01% of the population was infectious. We demonstrate good 3 week model forecasts of deaths with low error and good coverage of our credible intervals.

State-level tracking of COVID-19 in the United States

As of 1st June 2020, the US Centres for Disease Control and Prevention reported 104,232 confirmed or probable COVID-19-related deaths in the US. This was more than twice the number of deaths reported in the next most severely impacted country. We jointly model the US epidemic at the state-level, using publicly available death data within a Bayesian hierarchical semi-mechanistic framework. For each state, we estimate the number of individuals that have been infected, the number of individuals that are currently infectious and the time-varying reproduction number (the average number of secondary infections caused by an infected person). We use changes in mobility to capture the impact that non-pharmaceutical interventions and other behaviour changes have on the rate of transmission of SARS-CoV-2. We estimate that R_t was only below one in 23 states on 1st June. We also estimate that 3.7% [3.4%-4.0%] of the total population of the US had been infected, with wide variation between states, and approximately 0.01% of the population was infectious. We demonstrate good 3 week model forecasts of deaths with low error and good coverage of our credible intervals.

Response to COVID-19 in South Korea and implications for lifting stringent interventions

BACKGROUND

After experiencing a sharp growth in COVID-19 cases early in the pandemic, South Korea rapidly controlled transmission while implementing less stringent national social distancing measures than countries in Europe and the USA. This has led to substantial interest in their "test, trace, isolate" strategy. However, it is important to understand the epidemiological peculiarities of South Korea's outbreak and characterise their response before attempting to emulate these measures elsewhere.

METHODS

We systematically extracted numbers of suspected cases tested, PCR-confirmed cases, deaths, isolated confirmed cases, and numbers of confirmed cases with an identified epidemiological link from publicly available data. We estimated the time-varying reproduction number, R_t, using an established Bayesian framework, and reviewed the package of interventions implemented by South Korea using our extracted data, plus published literature and government sources.

RESULTS

We estimated that after the initial rapid growth in cases, R_t dropped below one in early April before increasing to a maximum of 1.94 (95%CrI, 1.64-2.27) in May following outbreaks in Seoul Metropolitan Region. By mid-June, R_t was back below one where it remained until the end of our study (July 13th). Despite less stringent "lockdown" measures, strong social distancing measures were implemented in high-incidence areas and studies measured a considerable national decrease in movement in late February. Testing the capacity was swiftly increased, and protocols were in place to isolate suspected and confirmed cases quickly; however, we could not estimate the delay to isolation using our data. Accounting for just 10% of cases, individual case-based contact tracing picked up a relatively minor proportion of total cases, with cluster investigations accounting for 66%.

CONCLUSIONS

Whilst early adoption of testing and contact tracing is likely to be important for South Korea's successful outbreak control, other factors including regional implementation of strong social distancing measures likely also contributed. The high volume of testing and the low number of deaths suggest that South Korea experienced a small epidemic relative to other countries. Caution is needed in attempting to replicate the South Korean response in populations with larger more geographically widespread epidemics where finding, testing, and isolating cases that are linked to clusters may be more difficult.

Response to COVID-19 in South Korea and implications for lifting stringent interventions

BACKGROUND

After experiencing a sharp growth in COVID-19 cases early in the pandemic, South Korea rapidly controlled transmission while implementing less stringent national social distancing measures than countries in Europe and the USA. This has led to substantial interest in their "test, trace, isolate" strategy. However, it is important to understand the epidemiological peculiarities of South Korea's outbreak and characterise their response before attempting to emulate these measures elsewhere.

METHODS

We systematically extracted numbers of suspected cases tested, PCR-confirmed cases, deaths, isolated confirmed cases, and numbers of confirmed cases with an identified epidemiological link from publicly available data. We estimated the time-varying reproduction number, Rt, using an established Bayesian framework, and reviewed the package of interventions implemented by South Korea using our extracted data, plus published literature and government sources.

RESULTS

We estimated that after the initial rapid growth in cases, Rt dropped below one in early April before increasing to a maximum of 1.94 (95%CrI, 1.64-2.27) in May following outbreaks in Seoul Metropolitan Region. By mid-June, Rt was back below one where it remained until the end of our study (July 13th). Despite less stringent "lockdown" measures, strong social distancing measures were implemented in high-incidence areas and studies measured a considerable national decrease in movement in late February. Testing the capacity was swiftly increased, and protocols were in place to isolate suspected and confirmed cases quickly; however, we could not estimate the delay to isolation using our data. Accounting for just 10% of cases, individual case-based contact tracing picked up a relatively minor proportion of total cases, with cluster investigations accounting for 66%.

CONCLUSIONS

Whilst early adoption of testing and contact tracing is likely to be important for South Korea's successful outbreak control, other factors including regional implementation of strong social distancing measures likely also contributed. The high volume of testing and the low number of deaths suggest that South Korea experienced a small epidemic relative to other countries. Caution is needed in attempting to replicate the South Korean response in populations with larger more geographically widespread epidemics where finding, testing, and isolating cases that are linked to clusters may be more difficult.

Anonymised and aggregated crowd level mobility data from mobile phones suggests that initial compliance with COVID-19 social distancing interventions was high and geographically consistent across the UK

Background: Since early March 2020, the COVID-19 epidemic across the United Kingdom has led to a range of social distancing policies, which have resulted in reduced mobility across different regions. Crowd level data on mobile phone usage can be used as a proxy for actual population mobility patterns and provide a way of quantifying the impact of social distancing measures on changes in mobility. Methods: Here, we use two mobile phone-based datasets (anonymised and aggregated crowd level data from O2 and from the Facebook app on mobile phones) to assess changes in average mobility, both overall and broken down into high and low population density areas, and changes in the distribution of journey lengths. Results: We show that there was a substantial overall reduction in mobility, with the most rapid decline on the 24th March 2020, the day after the Prime Minister's announcement of an enforced lockdown. The reduction in mobility was highly synchronized across the UK. Although mobility has remained low since 26th March 2020, we detect a gradual increase since that time. We also show that the two different datasets produce similar trends, albeit with some location-specific differences. We see slightly larger reductions in average mobility in high-density areas than in low-density areas, with greater variation in mobility in the high-density areas: some high-density areas eliminated almost all mobility. Conclusions: These analyses form a baseline from which to observe changes in behaviour in the UK as social distancing is eased and inform policy towards the future control of SARS-CoV-2 in the UK.

Case Fatality Ratio Estimates for the 2013-2016 West African Ebola Epidemic: Application of Boosted Regression Trees for Imputation

BACKGROUND

The 2013-2016 West African Ebola epidemic has been the largest to date with >11 000 deaths in the affected countries. The data collected have provided more insight into the case fatality ratio (CFR) and how it varies with age and other characteristics. However, the accuracy and precision of the naive CFR remain limited because 44% of survival outcomes were unreported.

METHODS

Using a boosted regression tree model, we imputed survival outcomes (ie, survival or death) when unreported, corrected for model imperfection to estimate the CFR without imputation, with imputation, and adjusted with imputation. The method allowed us to further identify and explore relevant clinical and demographic predictors of the CFR.

RESULTS

The out-of-sample performance (95% confidence interval [CI]) of our model was good: sensitivity, 69.7% (52.5-75.6%); specificity, 69.8% (54.1-75.6%); percentage correctly classified, 69.9% (53.7-75.5%); and area under the receiver operating characteristic curve, 76.0% (56.8-82.1%). The adjusted CFR estimates (95% CI) for the 2013-2016 West African epidemic were 82.8% (45.6-85.6%) overall and 89.1% (40.8-91.6%), 65.6% (61.3-69.6%), and 79.2% (45.4-84.1%) for Sierra Leone, Guinea, and Liberia, respectively. We found that district, hospitalisation status, age, case classification, and quarter (date of case reporting aggregated at three-month intervals) explained 93.6% of the variance in the naive CFR.

CONCLUSIONS

The adjusted CFR estimates improved the naive CFR estimates obtained without imputation and were more representative. Used in conjunction with other resources, adjusted estimates will inform public health contingency planning for future Ebola epidemics, and help better allocate resources and evaluate the effectiveness of future inventions.

Real-time Epidemic Forecasting: Challenges and Opportunities

Infectious disease outbreaks play an important role in global morbidity and mortality. Real-time epidemic forecasting provides an opportunity to predict geographic disease spread as well as case counts to better inform public health interventions when outbreaks occur. Challenges and recent advances in predictive modeling are discussed here. We identified data needs in the areas of epidemic surveillance, mobility, host and environmental susceptibility, pathogen transmissibility, population density, and healthcare capacity. Constraints in standardized case definitions and timely data sharing can limit the precision of predictive models. Resource-limited settings present particular challenges for accurate epidemic forecasting due to the lack of granular data available. Incorporating novel data streams into modeling efforts is an important consideration for the future as technology penetration continues to improve on a global level. Recent advances in machine-learning, increased collaboration between modelers, the use of stochastic semi-mechanistic models, real-time digital disease surveillance data, and open data sharing provide opportunities for refining forecasts for future epidemics. Epidemic forecasting using predictive modeling is an important tool for outbreak preparedness and response efforts. Despite the presence of some data gaps at present, opportunities and advancements in innovative data streams provide additional support for modeling future epidemics.

Distance is a barrier to recycling - or is it? Surprises from a clean test

The distance of recycling bins from households is often considered important by practitioners, but published evidence for this uses only indirect and self-reported data. This study aims to provide such evidence by obtaining a clean test using measured distances in a walled community with 1200 households with the same building types, local governance, recycling and waste arrangements. The number of deposits each month of food waste for recycling at a designated site are logged via smart-cards allocated per household. The number of days per month that each household deposits showed a highly significant - but small - negative correlation with distance of the bin: fewer householders participate if further away, accounting for 3% of the variation. Surprisingly, there is no variation with distance among those who do participate: their recycling frequency does not vary. This second result is not consistent with the first in terms of cost/benefit concepts assumed by government planners, nor with the static theories of behaviour currently used in waste management research. We recommend that recycling practitioners note the smallness of the contribution of distance to recycling performance, and not overrate it. And we recommend that researchers make better use of non-static models (which model different stages towards behaviour change), which our second result appears to call for.

Outbreak analytics: a developing data science for informing the response to emerging pathogens

Despite continued efforts to improve health systems worldwide, emerging pathogen epidemics remain a major public health concern. Effective response to such outbreaks relies on timely intervention, ideally informed by all available sources of data. The collection, visualization and analysis of outbreak data are becoming increasingly complex, owing to the diversity in types of data, questions and available methods to address them. Recent advances have led to the rise of outbreak analytics, an emerging data science focused on the technological and methodological aspects of the outbreak data pipeline, from collection to analysis, modelling and reporting to inform outbreak response. In this article, we assess the current state of the field. After laying out the context of outbreak response, we critically review the most common analytics components, their inter-dependencies, data requirements and the type of information they can provide to inform operations in real time. We discuss some challenges and opportunities and conclude on the potential role of outbreak analytics for improving our understanding of, and response to outbreaks of emerging pathogens.

This article is part of the theme issue ‘Modelling infectious disease outbreaks in humans, animals and plants: epidemic forecasting and control‘. This theme issue is linked with the earlier issue ‘Modelling infectious disease outbreaks in humans, animals and plants: approaches and important themes’.

Spatiotemporal variability in case fatality ratios for the 2013-2016 Ebola epidemic in West Africa

BACKGROUND

For the 2013-2016 Ebola epidemic in West Africa, the largest Ebola virus disease (EVD) epidemic to date, we aim to analyse the patient mix in detail to characterise key sources of spatiotemporal heterogeneity in the case fatality ratios (CFR).

METHODS

We applied a non-parametric Boosted Regression Trees (BRT) imputation approach for patients with missing survival outcomes and adjusted for model imperfection. Semivariogram analysis and kriging were used to investigate spatiotemporal heterogeneities.

RESULTS

CFR estimates varied significantly between districts and over time over the course of the epidemic. BRT modelling accounted for most of the spatiotemporal variation and interactions in CFR, but moderate spatial autocorrelation remained for distances up to approximately 90 km. Combining district-level CFR estimates and kriged district-level residuals provided the best linear unbiased predicted map of CFR accounting for the both explained and unexplained spatial variation. Temporal autocorrelation was not observed in the district-level residuals from the BRT estimates.

CONCLUSIONS

This study provides new insight into the epidemiology of the 2013-2016 West African Ebola epidemic with a view of informing future public health contingency planning, resource allocation and impact assessment. The analytical framework developed in this analysis, coupled with key domain knowledge, could be deployed in real time to support the response to ongoing and future outbreaks.

Rabies virus-neutralising antibodies in healthy, unvaccinated individuals: What do they mean for rabies epidemiology?

Rabies has been a widely feared disease for thousands of years, with records of rabid dogs as early as ancient Egyptian and Mesopotamian texts. The reputation of rabies as being inevitably fatal, together with its ability to affect all mammalian species, contributes to the fear surrounding this disease. However, the widely held view that exposure to the rabies virus is always fatal has been repeatedly challenged. Although survival following clinical infection in humans has only been recorded on a handful of occasions, a number of studies have reported detection of rabies-specific antibodies in the sera of humans, domestic animals, and wildlife that are apparently healthy and unvaccinated. These 'seropositive' individuals provide possible evidence of exposure to the rabies virus that has not led to fatal disease. However, the variability in methods of detecting these antibodies and the difficulties of interpreting serology tests have contributed to an unclear picture of their importance. In this review, we consider the evidence for rabies-specific antibodies in healthy, unvaccinated individuals as indicators of nonlethal rabies exposure and the potential implications of this for rabies epidemiology. Our findings indicate that whilst there is substantial evidence that nonlethal rabies exposure does occur, serology studies that do not use appropriate controls and cutoffs are unlikely to provide an accurate estimate of the true prevalence of nonlethal rabies exposure.

Complementary Paths to Chagas Disease Elimination: The Impact of Combining Vector Control With Etiological Treatment

Background

The World Health Organization’s 2020 goals for Chagas disease are (1) interrupting vector-borne intradomiciliary transmission and (2) having all infected people under care in endemic countries. Insecticide spraying has proved efficacious for reaching the first goal, but active transmission remains in several regions. For the second, treatment has mostly been restricted to recently infected patients, who comprise only a small proportion of all infected individuals.

Methods

We extended our previous dynamic transmission model to simulate a domestic Chagas disease transmission cycle and examined the effects of both vector control and etiological treatment on achieving the operational criterion proposed by the Pan American Health Organization for intradomiciliary, vectorial transmission interruption (ie, <2% seroprevalence in children <5 years of age).

Results

Depending on endemicity, an antivectorial intervention that decreases vector density by 90% annually would achieve the transmission interruption criterion in 2–3 years (low endemicity) to >30 years (high endemicity). When this strategy is combined with annual etiological treatment in 10% of the infected human population, the seroprevalence criterion would be achieved, respectively, in 1 and 11 years.

Conclusions

Combining highly effective vector control with etiological (trypanocidal) treatment in humans would substantially reduce time to transmission interruption as well as infection incidence and prevalence. However, the success of vector control may depend on prevailing vector species. It will be crucial to improve the coverage of screening programs, the performance of diagnostic tests, the proportion of people treated, and the efficacy of trypanocidal drugs. While screening and access can be incremented as part of strengthening the health systems response, improving diagnostics performance and drug efficacy will require further research.

Sub- or supercritical transmissibilities in a finite disease outbreak: Symmetry in outbreak properties of a disease conditioned on extinction

This work is concerned with the transmissibility of a disease, on observation of an outbreak of limited size. When such an outbreak occurs, an accurate estimate of the transmissibility of the responsible pathogen is essential for an appropriate response to future outbreaks. Transmissibility is usually characterised in terms of the reproduction number, R, which is the mean number of new cases of infection produced by a single infectious individual. A subcritical reproduction number (R < 1) guarantees that an outbreak will eventually die out of its own accord. By contrast, a supercritical reproduction number (R > 1) does not guarantee spread of the disease, since even with appreciable transmissibility, an outbreak may become extinct due to stochastic effects associated with a small number of infected individuals. Once the number of infectious individuals is conditioned on extinction, we show that an exact symmetry of the underlying theory ensures two distinct values of R, one larger than unity, the other smaller than unity, for which all outbreak properties are identical, with no signature of difference. Therefore a disease with a subcritical R, or its supercritical counterpart, when conditioned on extinction, have, for a given outbreak, identical individual likelihoods. In the full likelihood, this symmetry is lost, since the individual likelihood for a subcritical R is weighted by an extinction probability of unity, but the individual likelihood of a supercritical R is weighted by a sub-unity extinction probability. However, the inference can still benefit from the underlying symmetry, since it yields a mapping of all supercritical reproduction numbers onto the subcritical domain (R < 1), thereby speeding up evaluation of the likelihood profile. The symmetry holds in the standard situation, where the distribution of secondary cases is Poisson, as well as where this distribution has a negative binomial form and super-spreading can occur.

Genetic and spatial characterization of the red fox (Vulpes vulpes) population in the area stretching between the Eastern and Dinaric Alps and its relationship with rabies and canine distemper dynamics

Information on the population dynamics of a reservoir species have been increasingly adopted to understand and eventually predict the dispersal patterns of infectious diseases throughout an area. Although potentially relevant, to date there are no studies which have investigated the genetic structure of the red fox population in relation to infectious disease dynamics. Therefore, we genetically and spatially characterised the red fox population in the area stretching between the Eastern and Dinaric Alps, which has been affected by both distemper and rabies at different time intervals. Red foxes collected from north-eastern Italy, Austria, Slovenia and Croatia between 2006–2012, were studied using a set of 21 microsatellite markers. We confirmed a weak genetic differentiation within the fox population using Bayesian clustering analyses, and we were able to differentiate the fox population into geographically segregated groups. Our finding might be due to the presence of geographical barriers that have likely influenced the distribution of the fox population, limiting in turn gene flow and spread of infectious diseases. Focusing on the Italian red fox population, we observed interesting variations in the prevalence of both diseases among distinct fox clusters, with the previously identified Italy 1 and Italy 2 rabies as well as distemper viruses preferentially affecting different sub-groups identified in the study. Knowledge of the regional-scale population structure can improve understanding of the epidemiology and spread of diseases. Our study paves the way for an integrated approach for disease control coupling pathogen, host and environmental data to inform targeted control programs in the future.

A graph-based evidence synthesis approach to detecting outbreak clusters: An application to dog rabies

Early assessment of infectious disease outbreaks is key to implementing timely and effective control measures. In particular, rapidly recognising whether infected individuals stem from a single outbreak sustained by local transmission, or from repeated introductions, is crucial to adopt effective interventions. In this study, we introduce a new framework for combining several data streams, e.g. temporal, spatial and genetic data, to identify clusters of related cases of an infectious disease. Our method explicitly accounts for underreporting, and allows incorporating preexisting information about the disease, such as its serial interval, spatial kernel, and mutation rate. We define, for each data stream, a graph connecting all cases, with edges weighted by the corresponding pairwise distance between cases. Each graph is then pruned by removing distances greater than a given cutoff, defined based on preexisting information on the disease and assumptions on the reporting rate. The pruned graphs corresponding to different data streams are then merged by intersection to combine all data types; connected components define clusters of cases related for all types of data. Estimates of the reproduction number (the average number of secondary cases infected by an infectious individual in a large population), and the rate of importation of the disease into the population, are also derived. We test our approach on simulated data and illustrate it using data on dog rabies in Central African Republic. We show that the outbreak clusters identified using our method are consistent with structures previously identified by more complex, computationally intensive approaches.

Early assessment of infectious disease outbreaks is key to implementing timely and effective control measures. In particular, rapidly recognising whether infected individuals stem from a single outbreak sustained by local transmission, or from repeated introductions, is crucial to adopt effective interventions. In this study, we introduce a new approach which combines different types of data to identify clusters of related cases of an infectious disease. This approach relies on representing each type of data (e.g. temporal, spatial, or genetic) as a graph where nodes are cases, and two nodes are connected if the corresponding cases are closely related for this data. Our method then identifies clusters of cases which likely stem from the same introduction. Furthermore, we can use the size of these clusters to infer transmissibility of the disease and the number of importations of the pathogen into the population. We apply this approach to analyse dog rabies epidemics in Central African Republic. We show that outbreak clusters identified using our method are consistent with structures previously identified by more complex and computationally intensive approaches. Using simulated rabies epidemics, we show that our method has excellent potential for optimally detecting outbreak clusters. We also identify promising areas of research for transforming our method into a routine analysis tool for processing disease surveillance data.

Potential inconsistencies in Zika surveillance data and our understanding of risk during pregnancy

Background

A significant increase in microcephaly incidence was reported in Northeast Brazil at the end of 2015, which has since been attributed to an epidemic of Zika virus (ZIKV) infections earlier that year. Further incidence of congenital Zika syndrome (CZS) was expected following waves of ZIKV infection throughout Latin America; however, only modest increases in microcephaly and CZS incidence have since been observed. The quantitative relationship between ZIKV infection, gestational age and congenital outcome remains poorly understood.

Methodology/Principle findings

We characterised the gestational-age-varying risk of microcephaly given ZIKV infection using publicly available incidence data from multiple locations in Brazil and Colombia. We found that the relative timings and shapes of ZIKV infection and microcephaly incidence curves suggested different gestational risk profiles for different locations, varying in both the duration and magnitude of gestational risk. Data from Northeast Brazil suggested a narrow window of risk during the first trimester, whereas data from Colombia suggested persistent risk throughout pregnancy. We then used the model to estimate which combination of behavioural and reporting changes would have been sufficient to explain the absence of a second microcephaly incidence wave in Bahia, Brazil; a population for which we had two years of data. We found that a 18.9-fold increase in ZIKV infection reporting rate was consistent with observed patterns.

Conclusions

Our study illustrates how surveillance data may be used in principle to answer key questions in the absence of directed epidemiological studies. However, in this case, we suggest that currently available surveillance data are insufficient to accurately estimate the gestational-age-varying risk of microcephaly from ZIKV infection. The methods used here may be of use in future outbreaks and may help to inform improved surveillance and interpretation in countries yet to experience an outbreak of ZIKV infection.

Zika virus (ZIKV) infection is associated with the rise of microcephaly cases observed in Northeast Brazil at the end of 2015. For women in endemic or at-risk areas, understanding how the relationship between time of infection and microcephaly risk varies through pregnancy is important in informing family planning. However, a relatively modest number of congenital Zika syndrome cases have been observed following subsequent waves of ZIKV infection, limiting our understanding of gestational risk. We used a mathematical model to quantify the shape and magnitude of the gestational-age-varying risk to a fetus. Although the risk profile should be conserved regardless of location, we estimated different profiles when using surveillance data from locations in Northeast Brazil and Colombia. Our results suggest that time-dependent reporting changes likely confound the interpretation of currently available surveillance data. Furthermore, we investigated a range of behavioural and reporting rate changes that could explain two waves of ZIKV infection in Bahia, Brazil despite only one wave of microcephaly. Plausible changes in reporting could explain these data whilst remaining consistent with the hypothesis that ZIKV infection carries a significant risk of microcephaly. Further evidence is needed to disentangle the true risk of congenital Zika syndrome from time-varying reporting changes.

Outbreak of Ebola virus disease in the Democratic Republic of the Congo, April-May, 2018: an epidemiological study

BACKGROUND

On May 8, 2018, the Government of the Democratic Republic of the Congo reported an outbreak of Ebola virus disease in Équateur Province in the northwest of the country. The remoteness of most affected communities and the involvement of an urban centre connected to the capital city and neighbouring countries makes this outbreak the most complex and high risk ever experienced by the Democratic Republic of the Congo. We provide early epidemiological information arising from the ongoing investigation of this outbreak.

METHODS

We classified cases as suspected, probable, or confirmed using national case definitions of the Democratic Republic of the Congo Ministère de la Santé Publique. We investigated all cases to obtain demographic characteristics, determine possible exposures, describe signs and symptoms, and identify contacts to be followed up for 21 days. We also estimated the reproduction number and projected number of cases for the 4-week period from May 25, to June 21, 2018.

FINDINGS

As of May 30, 2018, 50 cases (37 confirmed, 13 probable) of Zaire ebolavirus were reported in the Democratic Republic of the Congo. 21 (42%) were reported in Bikoro, 25 (50%) in Iboko, and four (8%) in Wangata health zones. Wangata is part of Mbandaka, the urban capital of Équateur Province, which is connected to major national and international transport routes. By May 30, 2018, 25 deaths from Ebola virus disease had been reported, with a case fatality ratio of 56% (95% CI 39-72) after adjustment for censoring. This case fatality ratio is consistent with estimates for the 2014-16 west African Ebola virus disease epidemic (p=0·427). The median age of people with confirmed or probable infection was 40 years (range 8-80) and 30 (60%) were male. The most commonly reported signs and symptoms in people with confirmed or probable Ebola virus disease were fever (40 [95%] of 42 cases), intense general fatigue (37 [90%] of 41 cases), and loss of appetite (37 [90%] of 41 cases). Gastrointestinal symptoms were frequently reported, and 14 (33%) of 43 people reported haemorrhagic signs. Time from illness onset and hospitalisation to sample testing decreased over time. By May 30, 2018, 1458 contacts had been identified, of which 746 (51%) remained under active follow-up. The estimated reproduction number was 1·03 (95% credible interval 0·83-1·37) and the cumulative case incidence for the outbreak by June 21, 2018, is projected to be 78 confirmed cases (37-281), assuming heterogeneous transmissibility.

INTERPRETATION

The ongoing Ebola virus outbreak in the Democratic Republic of the Congo has similar epidemiological features to previous Ebola virus disease outbreaks. Early detection, rapid patient isolation, contact tracing, and the ongoing vaccination programme should sufficiently control the outbreak. The forecast of the number of cases does not exceed the current capacity to respond if the epidemiological situation does not change. The information presented, although preliminary, has been essential in guiding the ongoing investigation and response to this outbreak.

FUNDING

None.

Heterogeneities in the case fatality ratio in the West African Ebola outbreak 2013-2016

The 2013–2016 Ebola outbreak in West Africa is the largest on record with 28 616 confirmed, probable and suspected cases and 11 310 deaths officially recorded by 10 June 2016, the true burden probably considerably higher. The case fatality ratio (CFR: proportion of cases that are fatal) is a key indicator of disease severity useful for gauging the appropriate public health response and for evaluating treatment benefits, if estimated accurately. We analysed individual-level clinical outcome data from Guinea, Liberia and Sierra Leone officially reported to the World Health Organization. The overall mean CFR was 62.9% (95% CI: 61.9% to 64.0%) among confirmed cases with recorded clinical outcomes. Age was the most important modifier of survival probabilities, but country, stage of the epidemic and whether patients were hospitalized also played roles. We developed a statistical analysis to detect outliers in CFR between districts of residence and treatment centres (TCs), adjusting for known factors influencing survival and identified eight districts and three TCs with a CFR significantly different from the average. From the current dataset, we cannot determine whether the observed variation in CFR seen by district or treatment centre reflects real differences in survival, related to the quality of care or other factors or was caused by differences in reporting practices or case ascertainment.

This article is part of the themed issue ‘The 2013–2016 West African Ebola epidemic: data, decision-making and disease control’.

Modelling historical changes in the force-of-infection of Chagas disease to inform control and elimination programmes: application in Colombia

Background

WHO's 2020 milestones for Chagas disease include having all endemic Latin American countries certified with no intradomiciliary Trypanosoma cruzi transmission, and infected patients under care. Evaluating the variation in historical exposure to infection is crucial for assessing progress and for understanding the priorities to achieve these milestones.

Methods

Focusing on Colombia, all the available age-structured serological surveys (undertaken between 1995 and 2014) were searched and compiled. A total of 109 serosurveys were found, comprising 83 742 individuals from rural (indigenous and non-indigenous) and urban settings in 14 (out of 32) administrative units (departments). Estimates of the force-of-infection (FoI) were obtained by fitting and comparing three catalytic models using Bayesian methods to reconstruct temporal and spatial patterns over the course of three decades (between 1984 and 2014).

Results

Significant downward changes in the FoI were identified over the course of the three decades, and in some specific locations the predicted current seroprevalence in children aged 0-5 years is <1%. However, pronounced heterogeneity exists within departments, especially between indigenous, rural and urban settings, with the former exhibiting the highest FoI (up to 66 new infections/1000 people susceptible/year). The FoI in most of the indigenous settings remain unchanged during the three decades investigated. Current prevalence in adults in these 15 departments varies between 10% and 90% depending on the dynamics of historical exposure.

Conclusions

Assessing progress towards the control of Chagas disease requires quantifying the impact of historical exposure on current age-specific prevalence at subnational level. In Colombia, despite the evident progress, there is a marked heterogeneity indicating that in some areas the vector control interventions have not been effective, hindering the possibility of achieving interruption by 2020. A substantial burden of chronic cases remains even in locations where serological criteria for transmission interruption may have been achieved, therefore still demanding diagnosis and treatment interventions.

Key data for outbreak evaluation: building on the Ebola experience

Following the detection of an infectious disease outbreak, rapid epidemiological assessment is critical for guiding an effective public health response. To understand the transmission dynamics and potential impact of an outbreak, several types of data are necessary. Here we build on experience gained in the West African Ebola epidemic and prior emerging infectious disease outbreaks to set out a checklist of data needed to: (1) quantify severity and transmissibility; (2) characterize heterogeneities in transmission and their determinants; and (3) assess the effectiveness of different interventions. We differentiate data needs into individual-level data (e.g. a detailed list of reported cases), exposure data (e.g. identifying where/how cases may have been infected) and population-level data (e.g. size/demographics of the population(s) affected and when/where interventions were implemented). A remarkable amount of individual-level and exposure data was collected during the West African Ebola epidemic, which allowed the assessment of (1) and (2). However, gaps in population-level data (particularly around which interventions were applied when and where) posed challenges to the assessment of (3). Here we highlight recurrent data issues, give practical suggestions for addressing these issues and discuss priorities for improvements in data collection in future outbreaks.This article is part of the themed issue 'The 2013-2016 West African Ebola epidemic: data, decision-making and disease control'.