Rosalind M. Eggo

Transmission Dynamics of Zika Virus in Island Populations: A Modelling Analysis of the 2013–14 French Polynesia Outbreak

Between October 2013 and April 2014, more than 30,000 cases of Zika virus (ZIKV) disease were estimated to have attended healthcare facilities in French Polynesia. ZIKV has also been reported in Africa and Asia, and in 2015 the virus spread to South America and the Caribbean. Infection with ZIKV has been associated with neurological complications including Guillain-Barré Syndrome (GBS) and microcephaly, which led the World Health Organization to declare a Public Health Emergency of International Concern in February 2015. To better understand the transmission dynamics of ZIKV, we used a mathematical model to examine the 2013–14 outbreak on the six major archipelagos of French Polynesia. Our median estimates for the basic reproduction number ranged from 2.6–4.8, with an estimated 11.5% (95% CI: 7.32–17.9%) of total infections reported. As a result, we estimated that 94% (95% CI: 91–97%) of the total population of the six archipelagos were infected during the outbreak. Based on the demography of French Polynesia, our results imply that if ZIKV infection provides complete protection against future infection, it would take 12–20 years before there are a sufficient number of susceptible individuals for ZIKV to re-emerge, which is on the same timescale as the circulation of dengue virus serotypes in the region. Our analysis suggests that ZIKV may exhibit similar dynamics to dengue virus in island populations, with transmission characterized by large, sporadic outbreaks with a high proportion of asymptomatic or unreported cases.

Reconstructing disease outbreaks from genetic data: a graph approach

Epidemiology and public health planning will increasingly rely on the analysis of genetic sequence data. In particular, genetic data coupled with dates and locations of sampled isolates can be used to reconstruct the spatiotemporal dynamics of pathogens during outbreaks. Thus far, phylogenetic methods have been used to tackle this issue. Although these approaches have proved useful for informing on the spread of pathogens, they do not aim at directly reconstructing the underlying transmission tree. Instead, phylogenetic models infer most recent common ancestors between pairs of isolates, which can be inadequate for densely sampled recent outbreaks, where the sample includes ancestral and descendent isolates. In this paper, we introduce a novel method based on a graph approach to reconstruct transmission trees directly from genetic data. Using simulated data, we show that our approach can efficiently reconstruct genealogies of isolates in situations where classical phylogenetic approaches fail to do so. We then illustrate our method by analyzing data from the early stages of the swine-origin A/H1N1 influenza pandemic. Using 433 isolates sequenced at both the hemagglutinin and neuraminidase genes, we reconstruct the likely history of the worldwide spread of this new influenza strain. The presented methodology opens new perspectives for the analysis of genetic data in the context of disease outbreaks.

Spatial dynamics of the 1918 influenza pandemic in England, Wales and the United States

There is still limited understanding of key determinants of spatial spread of influenza. The 1918 pandemic provides an opportunity to elucidate spatial determinants of spread on a large scale. To better characterize the spread of the 1918 major wave, we fitted a range of city-to-city transmission models to mortality data collected for 246 population centres in England and Wales and 47 cities in the US. Using a gravity model for city-to-city contacts, we explored the effect of population size and distance on the spread of disease and tested assumptions regarding density dependence in connectivity between cities. We employed Bayesian Markov Chain Monte Carlo methods to estimate parameters of the model for population, infectivity, distance and density dependence. We inferred the most likely transmission trees for both countries. For England and Wales, a model that estimated the degree of density dependence in connectivity between cities was preferable by deviance information criterion comparison. Early in the major wave, long distance infective interactions predominated, with local infection events more likely as the epidemic became widespread. For the US, with fewer more widely dispersed cities, statistical power was lacking to estimate population size dependence or the degree of density dependence, with the preferred model depending on distance only. We find that parameters estimated from the England and Wales dataset can be applied to the US data with no likelihood penalty.

Effectiveness of Ring Vaccination as Control Strategy for Ebola Virus Disease.

Using an Ebola virus disease transmission model, we found that addition of ring vaccination at the outset of the West Africa epidemic might not have led to containment of this disease. However, in later stages of the epidemic or in outbreaks with less intense transmission or more effective control, this strategy could help eliminate the disease.

Comparative analysis of dengue and Zika outbreaks reveals differences by setting and virus

The pacific islands of Micronesia have experienced several outbreaks of mosquito-borne diseases over the past decade. In outbreaks on small islands, the susceptible population is usually well defined, and there is no co-circulation of pathogens. Because of this, analysing such outbreaks can be useful for understanding the transmission dynamics of the pathogens involved, and particularly so for yet understudied pathogens such as Zika virus. Here, we compared three outbreaks of dengue and Zika virus in two different island settings in Micronesia, the Yap Main Islands and Fais, using a mathematical model of transmission dynamics and making full use of commonalities in disease and setting between the outbreaks. We found that the estimated reproduction numbers for Zika and dengue were similar when considered in the same setting, but that, conversely, reproduction number for the same disease can vary considerably by setting. On the Yap Main Islands, we estimated a reproduction number of 8.0–16 (95% Credible Interval (CI)) for the dengue outbreak and 4.8–14 (95% CI) for the Zika outbreak, whereas for the dengue outbreak on Fais our estimate was 28–102 (95% CI). We further found that the proportion of cases of Zika reported was smaller (95% CI 1.4%–1.9%) than that of dengue (95% CI: 47%–61%). We confirmed these results in extensive sensitivity analysis. They suggest that models for dengue transmission can be useful for estimating the predicted dynamics of Zika transmission, but care must be taken when extrapolating findings from one setting to another.

Persistence and clearance of Ebola virus RNA from seminal fluid of Ebola virus disease survivors: a longitudinal analysis and modelling study

BACKGROUND By January, 2016, all known transmission chains of the Ebola virus disease (EVD) outbreak in west Africa had been stopped. However, there is concern about persistence of Ebola virus in the reproductive tract of men who have survived EVD. We aimed to use biostatistical modelling to describe the dynamics of Ebola virus RNA load in seminal fluid, including clearance parameters. METHODS In this longitudinal study, we recruited men who had been discharged from three Ebola treatment units in Guinea between January and July, 2015. Participants provided samples of seminal fluid at follow-up every 3-6 weeks, which we tested for Ebola virus RNA using quantitative real-time RT-PCR. Representative specimens from eight participants were then inoculated into immunodeficient mice to test for infectivity. We used a linear mixed-effect model to analyse the dynamics of virus persistence in seminal fluid over time. FINDINGS We enrolled 26 participants and tested 130 seminal fluid specimens; median follow up was 197 days (IQR 187-209 days) after enrolment, which corresponded to 255 days (228-287) after disease onset. Ebola virus RNA was detected in 86 semen specimens from 19 (73%) participants. Median duration of Ebola virus RNA detection was 158 days after onset (73-181; maximum 407 days at end of follow-up). Mathematical modelling of the quantitative time-series data showed a mean clearance rate of Ebola virus RNA from seminal fluid of -0·58 log units per month, although the clearance kinetic varied greatly between participants. Using our biostatistical model, we predict that 50% and 90% of male survivors clear Ebola virus RNA from seminal fluid at 115 days (90% prediction interval 72-160) and 294 days (212-399) after disease onset, respectively. We also predicted that the number of men positive for Ebola virus RNA in affected countries would decrease from about 50 in January 2016, to fewer than 1 person by July, 2016. Infectious virus was detected in 15 of 26 (58%) specimens tested in mice. INTERPRETATION Time to clearance of Ebola virus RNA from seminal fluid varies greatly between individuals and could be more than 13 months. Our predictions will assist in decision-making about surveillance and preventive measures in EVD outbreaks. FUNDING This study was funded by European Unions Horizon 2020 research and innovation programme, Directorate-General for International Cooperation and Development of the European Commission, Institut national de la santé et de la recherche médicale (INSERM), German Research Foundation (DFG), and Innovative Medicines Initiative 2 Joint Undertaking.

Respiratory virus transmission dynamics determine timing of asthma exacerbation peaks: Evidence from a population-level model.

Significance Asthma exacerbations are frequently triggered by common cold infections. Whereas prior studies have swabbed individual patients to detect viruses, we take a population level approach to investigate common cold circulation and its consequences for asthmatics. We introduce a dynamic model of common cold transmission with different contact patterns for adults and children, which are modified by school vacations. By jointly fitting this and an asthma risk model to daily hospitalization rates in eight large cities in Texas, we found that the common cold is a primary driver of asthma exacerbations, and both are predictably influenced by the school calendar. This provides actionable insight for assessing asthma hospitalization risk and targeting preventive measures. Asthma exacerbations exhibit a consistent annual pattern, closely mirroring the school calendar. Although respiratory viruses—the “common cold” viruses—are implicated as a principal cause, there is little evidence to link viral prevalence to seasonal differences in risk. We jointly fit a common cold transmission model and a model of biological and environmental exacerbation triggers to estimate effects on hospitalization risk. Asthma hospitalization rate, influenza prevalence, and air quality measures are available, but common cold circulation is not; therefore, we generate estimates of viral prevalence using a transmission model. Our deterministic multivirus transmission model includes transmission rates that vary when school is closed. We jointly fit the two models to 7 y of daily asthma hospitalizations in adults and children (66,000 events) in eight metropolitan areas. For children, we find that daily viral prevalence is the strongest predictor of asthma hospitalizations, with transmission reduced by 45% (95% credible interval =41–49%) during school closures. We detect a transient period of nonspecific immunity between infections lasting 19 (17–21) d. For adults, hospitalizations are more variable, with influenza driving wintertime peaks. Neither particulate matter nor ozone was an important predictor, perhaps because of the large geographic area of the populations. The school calendar clearly and predictably drives seasonal variation in common cold prevalence, which results in the “back-to-school” asthma exacerbation pattern seen in children and indirectly contributes to exacerbation risk in adults. This study provides a framework for anticipating the seasonal dynamics of common colds and the associated risks for asthmatics.

Duration of Ebola virus RNA persistence in semen of survivors: population-level estimates and projections.

Ebola virus can persist in semen after recovery, potentially for months, which may impact the duration of enhanced surveillance required after interruption of transmission. We combined recent data on viral RNA persistence with weekly disease incidence to estimate the current number of semen-positive men in affected West African countries. We find the number is low, and since few reported sexual transmission events have occurred, the future risk is also likely low, although sexual health promotion remains critical.

Real-time forecasting of infectious disease dynamics with a stochastic semi-mechanistic model

Highlights • A Bayesian semi-mechanistic model was applied to the Ebola Forecasting Challenge.• Model fits to simulated data were obtained from particle Markov-chain Monte Carlo.• Posterior samples of model parameters were used to generate forecast trajectories.• The forecasts were assessed using subsequently released simulation points.

Cholera epidemic in Yemen, 2016–18: an analysis of surveillance data

Summary Background In war-torn Yemen, reports of confirmed cholera started in late September, 2016. The disease continues to plague Yemen today in what has become the largest documented cholera epidemic of modern times. We aimed to describe the key epidemiological features of this epidemic, including the drivers of cholera transmission during the outbreak. Methods The Yemen Health Authorities set up a national cholera surveillance system to collect information on suspected cholera cases presenting at health facilities. Individual variables included symptom onset date, age, severity of dehydration, and rapid diagnostic test result. Suspected cholera cases were confirmed by culture, and a subset of samples had additional phenotypic and genotypic analysis. We first conducted descriptive analyses at national and governorate levels. We divided the epidemic into three time periods: the first wave (Sept 28, 2016, to April 23, 2017), the increasing phase of the second wave (April 24, 2017, to July 2, 2017), and the decreasing phase of the second wave (July 3, 2017, to March 12, 2018). We reconstructed the changes in cholera transmission over time by estimating the instantaneous reproduction number, Rt. Finally, we estimated the association between rainfall and the daily cholera incidence during the increasing phase of the second epidemic wave by fitting a spatiotemporal regression model. Findings From Sept 28, 2016, to March 12, 2018, 1 103 683 suspected cholera cases (attack rate 3·69%) and 2385 deaths (case fatality risk 0·22%) were reported countrywide. The epidemic consisted of two distinct waves with a surge in transmission in May, 2017, corresponding to a median Rt of more than 2 in 13 of 23 governorates. Microbiological analyses suggested that the same Vibrio cholerae O1 Ogawa strain circulated in both waves. We found a positive, non-linear, association between weekly rainfall and suspected cholera incidence in the following 10 days; the relative risk of cholera after a weekly rainfall of 25 mm was 1·42 (95% CI 1·31–1·55) compared with a week without rain. Interpretation Our analysis suggests that the small first cholera epidemic wave seeded cholera across Yemen during the dry season. When the rains returned in April, 2017, they triggered widespread cholera transmission that led to the large second wave. These results suggest that cholera could resurge during the ongoing 2018 rainy season if transmission remains active. Therefore, health authorities and partners should immediately enhance current control efforts to mitigate the risk of a new cholera epidemic wave in Yemen. Funding Health Authorities of Yemen, WHO, and Médecins Sans Frontières.