Spatial and temporal invasion dynamics of the 2014-2017 Zika and chikungunya epidemics in Colombia

The effects of seasonal human mobility and Aedes aegypti habitat suitability on Zika virus epidemic severity in Colombia

The Zika virus epidemic of 2015-16, which caused over 1 million confirmed or suspected human cases in the Caribbean and Latin America, was driven by a combination of movement of infected humans and availability of suitable habitat for mosquito species that are key disease vectors. Both human mobility and mosquito vector abundances vary seasonally, and the goal of our research was to analyze the interacting effects of disease vector densities and human movement across metapopulations on disease transmission intensity and the probability of super-spreader events. Our research uses the novel approach of combining geographical modeling of mosquito presence with network modeling of human mobility to offer a comprehensive simulation environment for Zika virus epidemics that considers a substantial number of spatial and temporal factors compared to the literature. Specifically, we tested the hypotheses that 1) regions with the highest probability of mosquito presence will have more super-spreader events during dry months, when mosquitoes are predicted to be more abundant, 2) regions reliant on tourism industries will have more super-spreader events during wet months, when they are more likely to contribute to network-level pathogen spread due to increased travel. We used the case study of Colombia, a country with a population of about 50 million people, with an annual calendar that can be partitioned into overlapping cycles of wet and dry seasons and peak tourism and off tourism seasons that drive distinct cyclical patterns of mosquito abundance and human movement. Our results show that whether the first infected human was introduced to the network during the wet versus dry season and during the tourism versus off tourism season profoundly affects the severity and trajectory of the epidemic. For example, Zika virus was first detected in Colombia in October of 2015. Had it originated in January, a dry season month with high rates of tourism, it likely could have infected up to 60% more individuals and up to 40% more super-spreader events may have occurred. In addition, popular tourism destinations such as Barranquilla and Cartagena have the highest risk of super-spreader events during the winter, whereas densely populated areas such as Medellín and Bogotá are at higher risk of sustained transmission during dry months in the summer. Our research demonstrates that public health planning and response to vector-borne disease outbreaks requires a thorough understanding of how vector and host patterns vary due to seasonality in environmental conditions and human mobility dynamics. This research also has strong implications for tourism policy and the potential response strategies in case of an emergent epidemic.

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.

Spatio-temporal clusters and patterns of spread of dengue, chikungunya, and Zika in Colombia

Background

Colombia has one of the highest burdens of arboviruses in South America. The country was in a state of hyperendemicity between 2014 and 2016, with co-circulation of several Aedes-borne viruses, including a syndemic of dengue, chikungunya, and Zika in 2015.

Methodology/Principal findings

We analyzed the cases of dengue, chikungunya, and Zika notified in Colombia from January 2014 to December 2018 by municipality and week. The trajectory and velocity of spread was studied using trend surface analysis, and spatio-temporal high-risk clusters for each disease in separate and for the three diseases simultaneously (multivariate) were identified using Kulldorff’s scan statistics. During the study period, there were 366,628, 77,345 and 74,793 cases of dengue, chikungunya, and Zika, respectively, in Colombia. The spread patterns for chikungunya and Zika were similar, although Zika’s spread was accelerated. Both chikungunya and Zika mainly spread from the regions on the Atlantic coast and the south-west to the rest of the country. We identified 21, 16, and 13 spatio-temporal clusters of dengue, chikungunya and Zika, respectively, and, from the multivariate analysis, 20 spatio-temporal clusters, among which 7 were simultaneous for the three diseases. For all disease-specific analyses and the multivariate analysis, the most-likely cluster was identified in the south-western region of Colombia, including the Valle del Cauca department.

Conclusions/Significance

The results further our understanding of emerging Aedes-borne diseases in Colombia by providing useful evidence on their potential site of entry and spread trajectory within the country, and identifying spatio-temporal disease-specific and multivariate high-risk clusters of dengue, chikungunya, and Zika, information that can be used to target interventions.

Dengue, chikungunya, and Zika are diseases transmitted to humans by the bite of infected Aedes mosquitoes. Between 2014 and 2016 chikungunya and Zika viruses started causing outbreaks in Colombia, one of the countries historically most affected by dengue. We used case counts of the diseases by municipality and week to study the spread trajectory of chikungunya and Zika within Colombia’s territory, and to identify space-time high-risk clusters, i.e., the areas and time periods that dengue, chikungunya, and Zika were more present. Chikungunya and Zika spread similarly in Colombia, but Zika spread faster. The Atlantic coast, a famous touristic destination in the country, was likely the place of entry of chikungunya and Zika in Colombia. The south-western region was identified as a high-risk cluster for all three diseases in separate and simultaneously. This region has a favorable climate for the Aedes mosquitoes and other characteristics that facilitate the diseases’ transmission, such as social deprivation and high population mobility. Our results provide useful information on the locations that should be prioritized for interventions to prevent the entry of new diseases transmitted by Aedes and to reduce the burden of dengue, chikungunya and Zika where they are established.