Projecting vaccine demand and impact for emerging zoonotic pathogens.

BACKGROUND: Despite large outbreaks in humans seeming improbable for a number of zoonotic pathogens, several pose a concern due to their epidemiological characteristics and evolutionary potential. To enable effective responses to these pathogens in the event that they undergo future emergence, the Coalition for Epidemic Preparedness Innovations is advancing the development of vaccines for several pathogens prioritized by the World Health Organization. A major challenge in this pursuit is anticipating demand for a vaccine stockpile to support outbreak response. METHODS: We developed a modeling framework for outbreak response for emerging zoonoses under three reactive vaccination strategies to assess sustainable vaccine manufacturing needs, vaccine stockpile requirements, and the potential impact of the outbreak response. This framework incorporates geographically variable zoonotic spillover rates, human-to-human transmission, and the implementation of reactive vaccination campaigns in response to disease outbreaks. As proof of concept, we applied the framework to four priority pathogens: Lassa virus, Nipah virus, MERS coronavirus, and Rift Valley virus. RESULTS: Annual vaccine regimen requirements for a population-wide strategy ranged from > 670,000 (95% prediction interval 0-3,630,000) regimens for Lassa virus to 1,190,000 (95% PrI 0-8,480,000) regimens for Rift Valley fever virus, while the regimens required for ring vaccination or targeting healthcare workers (HCWs) were several orders of magnitude lower (between 1/25 and 1/700) than those required by a population-wide strategy. For each pathogen and vaccination strategy, reactive vaccination typically prevented fewer than 10% of cases, because of their presently low R0 values. Targeting HCWs had a higher per-regimen impact than population-wide vaccination. CONCLUSIONS: Our framework provides a flexible methodology for estimating vaccine stockpile needs and the geographic distribution of demand under a range of outbreak response scenarios. Uncertainties in our model estimates highlight several knowledge gaps that need to be addressed to target vulnerable populations more accurately. These include surveillance gaps that mask the true geographic distribution of each pathogen, details of key routes of spillover from animal reservoirs to humans, and the role of human-to-human transmission outside of healthcare settings. In addition, our estimates are based on the current epidemiology of each pathogen, but pathogen evolution could alter vaccine stockpile requirements.

Burden is in the eye of the beholder: Sensitivity of yellow fever disease burden estimates to modeling assumptions.

Estimates of disease burden are important for setting public health priorities. These estimates involve numerous modeling assumptions, whose uncertainties are not always well described. We developed a framework for estimating the burden of yellow fever in Africa and evaluated its sensitivity to modeling assumptions that are often overlooked. We found that alternative interpretations of serological data resulted in a nearly 20-fold difference in burden estimates (range of central estimates, 8.4 × 104 to 1.5 × 106 deaths in 2021­2030). Uncertainty about the vaccination status of serological study participants was the primary driver of this uncertainty. Even so, statistical uncertainty was even greater than uncertainty due to modeling assumptions, accounting for a total of 87% of variance in burden estimates. Combined with estimates that most infections go unreported (range of 95% credible intervals, 99.65 to 99.99%), our results suggest that yellow fever's burden will remain highly uncertain without major improvements in surveillance.

Trade-offs between individual and ensemble forecasts of an emerging infectious disease.

Probabilistic forecasts play an indispensable role in answering questions about the spread of newly emerged pathogens. However, uncertainties about the epidemiology of emerging pathogens can make it difficult to choose among alternative model structures and assumptions. To assess the potential for uncertainties about emerging pathogens to affect forecasts of their spread, we evaluated the performance 16 forecasting models in the context of the 2015-2016 Zika epidemic in Colombia. Each model featured a different combination of assumptions about human mobility, spatiotemporal variation in transmission potential, and the number of virus introductions. We found that which model assumptions had the most ensemble weight changed through time. We additionally identified a trade-off whereby some individual models outperformed ensemble models early in the epidemic, but on average the ensembles outperformed all individual models. Our results suggest that multiple models spanning uncertainty across alternative assumptions are necessary to obtain robust forecasts for emerging infectious diseases.

Leveraging multiple data types to estimate the size of the Zika epidemic in the Americas.

Several hundred thousand Zika cases have been reported across the Americas since 2015. Incidence of infection was likely much higher, however, due to a high frequency of asymptomatic infection and other challenges that surveillance systems faced. Using a hierarchical Bayesian model with empirically-informed priors, we leveraged multiple types of Zika case data from 15 countries to estimate subnational reporting probabilities and infection attack rates (IARs). Zika IAR estimates ranged from 0.084 (95% CrI: 0.067-0.096) in Peru to 0.361 (95% CrI: 0.214-0.514) in Ecuador, with significant subnational variability in every country. Totaling infection estimates across these and 33 other countries and territories, our results suggest that 132.3 million (95% CrI: 111.3-170.2 million) people in the Americas had been infected by the end of 2018. These estimates represent the most extensive attempt to determine the size of the Zika epidemic in the Americas, offering a baseline for assessing the risk of future Zika epidemics in this region.

Modeling human migration across spatial scales in Colombia.

Human mobility, both short and long term, are important considerations in the study of numerous systems. Economic and technological advances have led to a more interconnected global community, further increasing the need for considerations of human mobility. While data on human mobility are better recorded in many developed countries, availability of such data remains limited in many low- and middle-income countries around the world, particularly at the fine temporal and spatial scales required by many applications. In this study, we used 5-year census-based internal migration microdata for 32 departments in Colombia (i.e., Admin-1 level) to develop a novel spatial interaction modeling approach for estimating migration, at a finer spatial scale, among the 1,122 municipalities in the country (i.e., Admin-2 level). Our modeling approach addresses a significant lack of migration data at administrative unit levels finer than those at which migration data are typically recorded. Due to the widespread availability of census-based migration microdata at the Admin-1 level, our modeling approach opens up for the possibilities of modeling migration patterns at Admin-2 and Admin-3 levels across many other countries where such data are currently lacking.

Heterogeneous local dynamics revealed by classification analysis of spatially disaggregated time series data.

Time series data provide a crucial window into infectious disease dynamics, yet their utility is often limited by the spatially aggregated form in which they are presented. When working with time series data, violating the implicit assumption of homogeneous dynamics below the scale of spatial aggregation could bias inferences about underlying processes. We tested this assumption in the context of the 2015-2016 Zika epidemic in Colombia, where time series of weekly case reports were available at national, departmental, and municipal scales. First, we performed a descriptive analysis, which showed that the timing of departmental-level epidemic peaks varied by three months and that departmental-level estimates of the time-varying reproduction number, R(t), showed patterns that were distinct from a national-level estimate. Second, we applied a classification algorithm to six features of proportional cumulative incidence curves, which showed that variability in epidemic duration, the length of the epidemic tail, and consistency with a cumulative normal density curve made the greatest contributions to distinguishing groups. Third, we applied this classification algorithm to data simulated with a stochastic transmission model, which showed that group assignments were consistent with simulated differences in the basic reproduction number, R0. This result, along with associations between spatial drivers of transmission and group assignments based on observed data, suggests that the classification algorithm is capable of detecting differences in temporal patterns that are associated with differences in underlying drivers of incidence patterns. Overall, this diversity of temporal patterns at local scales underscores the value of spatially disaggregated time series data.

Inter-annual variation in seasonal dengue epidemics driven by multiple interacting factors in Guangzhou, China.

Vector-borne diseases display wide inter-annual variation in seasonal epidemic size due to their complex dependence on temporally variable environmental conditions and other factors. In 2014, Guangzhou, China experienced its worst dengue epidemic on record, with incidence exceeding the historical average by two orders of magnitude. To disentangle contributions from multiple factors to inter-annual variation in epidemic size, we fitted a semi-mechanistic model to time series data from 2005-2015 and performed a series of factorial simulation experiments in which seasonal epidemics were simulated under all combinations of year-specific patterns of four time-varying factors: imported cases, mosquito density, temperature, and residual variation in local conditions not explicitly represented in the model. Our results indicate that while epidemics in most years were limited by unfavorable conditions with respect to one or more factors, the epidemic in 2014 was made possible by the combination of favorable conditions for all factors considered in our analysis.

Local and regional dynamics of chikungunya virus transmission in Colombia: the role of mismatched spatial heterogeneity.

BACKGROUND: Mathematical models of transmission dynamics are routinely fitted to epidemiological time series, which must inevitably be aggregated at some spatial scale. Weekly case reports of chikungunya have been made available nationally for numerous countries in the Western Hemisphere since late 2013, and numerous models have made use of this data set for forecasting and inferential purposes. Motivated by an abundance of literature suggesting that the transmission of this mosquito-borne pathogen is localized at scales much finer than nationally, we fitted models at three different spatial scales to weekly case reports from Colombia to explore limitations of analyses of nationally aggregated time series data. METHODS: We adapted the recently developed Disease Transmission Kernel (DTK)-Dengue model for modeling chikungunya virus (CHIKV) transmission, given the numerous similarities of these viruses vectored by a common mosquito vector. We fitted versions of this model specified at different spatial scales to weekly case reports aggregated at different spatial scales: (1) single-patch national model fitted to national data; (2) single-patch departmental models fitted to departmental data; and (3) multi-patch departmental models fitted to departmental data, where the multiple patches refer to municipalities within a department. We compared the consistency of simulations from fitted models with empirical data. RESULTS: We found that model consistency with epidemic dynamics improved with increasing spatial granularity of the model. Specifically, the sum of single-patch departmental model fits better captured national-level temporal patterns than did a single-patch national model. Likewise, multi-patch departmental model fits better captured department-level temporal patterns than did single-patch departmental model fits. Furthermore, inferences about municipal-level incidence based on multi-patch departmental models fitted to department-level data were positively correlated with municipal-level data that were withheld from model fitting. CONCLUSIONS: Our model performed better when posed at finer spatial scales, due to better matching between human populations with locally relevant risk. Confronting spatially aggregated models with spatially aggregated data imposes a serious structural constraint on model behavior by averaging over epidemiologically meaningful spatial variation in drivers of transmission, impairing the ability of models to reproduce empirical patterns.

Spatiotemporal incidence of Zika and associated environmental drivers for the 2015-2016 epidemic in Colombia.

Despite a long history of mosquito-borne virus epidemics in the Americas, the impact of the Zika virus (ZIKV) epidemic of 2015-2016 was unexpected. The need for scientifically informed decision-making is driving research to understand the emergence and spread of ZIKV. To support that research, we assembled a data set of key covariates for modeling ZIKV transmission dynamics in Colombia, where ZIKV transmission was widespread and the government made incidence data publically available. On a weekly basis between January 1, 2014 and October 1, 2016 at three administrative levels, we collated spatiotemporal Zika incidence data, nine environmental variables, and demographic data into a single downloadable database. These new datasets and those we identified, processed, and assembled at comparable spatial and temporal resolutions will save future researchers considerable time and effort in performing these data processing steps, enabling them to focus instead on extracting epidemiological insights from this important data set. Similar approaches could prove useful for filling data gaps to enable epidemiological analyses of future disease emergence events.

Assessing the population at risk of Zika virus in Asia - is the emergency really over?

On November 18, 2016, the WHO ended its designation of the Zika virus (ZIKV) epidemic as a Public Health Emergency of International Concern (PHEIC). At the same time, ZIKV transmission continues in Asia, with the number of Asian countries reporting Zika cases increasing over the last 2 years. Applying a method that combines epidemiological theory with data on epidemic size and drivers of transmission, we characterised the population at risk of ZIKV infection from Aedes aegypti mosquitoes in 15 countries in Asia. Projections made under the assumption of no pre-existing immunity suggest that up to 785 (range: 730-992) million people in Asia would be at risk of ZIKV infection under that scenario. Assuming that 20% of ZIKV infections are symptomatic, this implies an upper limit of 146-198 million for the population at risk of a clinical episode of Zika. Due to limited information about pre-existing immunity to ZIKV in the region, we were unable to make specific numerical projections under a more realistic assumption about pre-existing immunity. Even so, combining numerical projections under an assumption of no pre-existing immunity together with theoretical insights about the extent to which pre-existing immunity may lower epidemic size, our results suggest that the population at risk of ZIKV infection in Asia could be even larger than in the Americas. As a result, we conclude that the WHO's removal of the PHEIC designation should not be interpreted as an indication that the threat posed by ZIKV has subsided.

Temperature modulates dengue virus epidemic growth rates through its effects on reproduction numbers and generation intervals.

Epidemic growth rate, r, provides a more complete description of the potential for epidemics than the more commonly studied basic reproduction number, R0, yet the former has never been described as a function of temperature for dengue virus or other pathogens with temperature-sensitive transmission. The need to understand the drivers of epidemics of these pathogens is acute, with arthropod-borne virus epidemics becoming increasingly problematic. We addressed this need by developing temperature-dependent descriptions of the two components of r-R0 and the generation interval-to obtain a temperature-dependent description of r. Our results show that the generation interval is highly sensitive to temperature, decreasing twofold between 25 and 35°C and suggesting that dengue virus epidemics may accelerate as temperatures increase, not only because of more infections per generation but also because of faster generations. Under the empirical temperature relationships that we considered, we found that r peaked at a temperature threshold that was robust to uncertainty in model parameters that do not depend on temperature. Although the precise value of this temperature threshold could be refined following future studies of empirical temperature relationships, the framework we present for identifying such temperature thresholds offers a new way to classify regions in which dengue virus epidemic intensity could either increase or decrease under future climate change.

Model-based projections of Zika virus infections in childbearing women in the Americas.

Zika virus is a mosquito-borne pathogen that is rapidly spreading across the Americas. Due to associations between Zika virus infection and a range of fetal maladies(1,2), the epidemic trajectory of this viral infection poses a significant concern for the nearly 15 million children born in the Americas each year. Ascertaining the portion of this population that is truly at risk is an important priority. One recent estimate(3) suggested that 5.42 million childbearing women live in areas of the Americas that are suitable for Zika occurrence. To improve on that estimate, which did not take into account the protective effects of herd immunity, we developed a new approach that combines classic results from epidemiological theory with seroprevalence data and highly spatially resolved data about drivers of transmission to make location-specific projections of epidemic attack rates. Our results suggest that 1.65 (1.45-2.06) million childbearing women and 93.4 (81.6-117.1) million people in total could become infected before the first wave of the epidemic concludes. Based on current estimates of rates of adverse fetal outcomes among infected women(2,4,5), these results suggest that tens of thousands of pregnancies could be negatively impacted by the first wave of the epidemic. These projections constitute a revised upper limit of populations at risk in the current Zika epidemic, and our approach offers a new way to make rapid assessments of the threat posed by emerging infectious diseases more generally.

Temperature and population density determine reservoir regions of seasonal persistence in highland malaria.

A better understanding of malaria persistence in highly seasonal environments such as highlands and desert fringes requires identifying the factors behind the spatial reservoir of the pathogen in the low season. In these 'unstable' malaria regions, such reservoirs play a critical role by allowing persistence during the low transmission season and therefore, between seasonal outbreaks. In the highlands of East Africa, the most populated epidemic regions in Africa, temperature is expected to be intimately connected to where in space the disease is able to persist because of pronounced altitudinal gradients. Here, we explore other environmental and demographic factors that may contribute to malaria's highland reservoir. We use an extensive spatio-temporal dataset of confirmed monthly Plasmodium falciparum cases from 1995 to 2005 that finely resolves space in an Ethiopian highland. With a Bayesian approach for parameter estimation and a generalized linear mixed model that includes a spatially structured random effect, we demonstrate that population density is important to disease persistence during the low transmission season. This population effect is not accounted for in typical models for the transmission dynamics of the disease, but is consistent in part with a more complex functional form of the force of infection proposed by theory for vector-borne infections, only during the low season as we discuss. As malaria risk usually decreases in more urban environments with increased human densities, the opposite counterintuitive finding identifies novel control targets during the low transmission season in African highlands.