Climatic fluctuations and malaria transmission dynamics, prior to elimination, in Guna Yala, República de Panamá

Trypanosoma cruzi infection in dogs along the US-Mexico border: R0 changes with vector species composition

Infection with Trypanosoma cruzi, etiological agent of Chagas disease, is common in US government working dogs along the US-Mexico border. This 3145 km long border comprises four states: Texas (TX), New Mexico (NM), Arizona (AZ) and California (CA) with diverse ecosystems and several triatomine (a.k.a., kissing bug) species, primary vectors of T. cruzi in this region. The kissing bug (Heteroptera: Reduviidae) community ranging from CA to TX includes Triatoma protracta (Uhler), Triatoma recurva (Stål) and Triatoma rubida (Uhler) and becomes dominated by Triatoma gerstaeckeri Stål in TX. Here, we ask if T. cruzi infection dynamics in dogs varies along this border region, potentially reflecting changes in vector species and their vectorial capacity. Using reversible catalytic models of infection, where seropositivity can be lost, we estimated an R0 (Estimate ± S.E.) of 1.192 ± 0.084 for TX and NM. In contrast, seropositivity decayed to zero as dogs aged in AZ and CA. These results suggest that dogs are likely infected by T. cruzi during their training in western TX, with a force of infection large enough for keeping R0 above 1, i.e., the disease endemically established, in TX and NM. In AZ and CA, a lower force of infection, probably associated with different vector species communities and associated vectorial capacity and/or different lineages of T. cruzi, results in dogs decreasing their seropositivity with age.

rspatialdata: a collection of data sources and tutorials on downloading and visualising spatial data using R

Spatial and spatio-temporal data are used in a wide range of fields including environmental, health and social disciplines. Several packages in the statistical software R have been recently developed as clients for various databases to meet the growing demands for easily accessible and reliable spatial data. While documentation on how to use many of these packages exist, there is an increasing need for a one stop repository for tutorials on this information. In this paper, we present rspatialdata a website that provides a collection of data sources and tutorials on downloading and visualising spatial data using R. The website includes a wide range of datasets including administrative boundaries of countries, Open Street Map data, population, temperature, vegetation, air pollution, and malaria data. The goal of the website is to equip researchers and communities with the tools to engage in spatial data analysis and visualisation so that they can address important local issues, such as estimating air pollution, quantifying disease burdens, and evaluating and monitoring the United Nation’s sustainable development goals.

Anopheles drivers of persisting malaria transmission in Guna Yala, Panamá: an operational investigation

Background

Though most of Panamá is free from malaria, localized foci of transmission persist, including in the Guna Yala region. Government-led entomological surveillance using an entomological surveillance planning tool (ESPT) sought to answer programmatically-relevant questions that would enhance the understanding of both local entomological drivers of transmission and gaps in protection that result in persisting malaria transmission to guide local vector control decision-making.

Methods

The ESPT was used to design a sampling plan centered around the collection of minimum essential indicators to investigate the relevance of LLINs and IRS in the communities of Permé and Puerto Obaldía, Guna Yala, as well as to pinpoint any remaining spaces and times where humans are exposed to Anopheles bites (gaps in protection). Adult Anopheles were collected at three time points via human landing catches (HLCs), CDC Light Traps (LT), and pyrethrum spray catches (PSCs) during the rainy and dry seasons. Mosquitoes were identified to species via molecular methods. Insecticide susceptibility testing of the main vector species to fenitrothion was conducted.

Results

In total, 7537 adult Anopheles were collected from both sites. Of the 493 specimens molecularly confirmed to species, two thirds (n = 340) were identified as Nyssorhynchus albimanus, followed by Anopheles aquasalis. Overall Anopheles human biting rates (HBRs) were higher outdoors than indoors, and were higher in Permé than in Puerto Obaldía: nightly outdoor HBR ranged from 2.71 bites per person per night (bpn) (Puerto Obaldía), to 221.00 bpn (Permé), whereas indoor nightly HBR ranged from 0.70 bpn (Puerto Obaldía) to 81.90 bpn (Permé). Generally, peak biting occurred during the early evening. The CDC LT trap yields were significantly lower than that of HLCs and this collection method was dropped after the first collection. Pyrethrum spray catches resulted in only three indoor resting Anopheles collected. Insecticide resistance (IR) of Ny. albimanus to fenitrothion was confirmed, with only 65.5% mortality at the diagnostic time.

Conclusion

The early evening exophagic behaviour of Anopheles vectors, the absence of indoor resting behaviours, and the presence of resistance to the primary intervention insecticide demonstrate limitations of the current malaria strategy, including indoor residual spraying (IRS) and long-lasting insecticidal nets (LLINs), and point to both gaps in protection and to the drivers of persisting malaria transmission in Guna Yala. These findings highlight the need for continued and directed entomological surveillance, based on programmatic questions, that generates entomological evidence to inform an adaptive malaria elimination strategy.

The association between gold mining and malaria in Guyana: a statistical inference and time-series analysis

BACKGROUND

Guyana reported a significant rise in malaria between 2008 and 2014. As there was no evidence of impairment of national malaria control strategies, public health authorities attributed the surge to a temporal increase in gold mining activity in forested regions. However, systematic analysis of this association is lacking because of the difficulties associated with collecting reliable data for both malaria and mining. We aimed to investigate the association between the international gold price and Plasmodium falciparum malaria transmission in Guyana between 2007 and 2019. We also aimed to evaluate the association between P falciparum cases and the El Niño-Southern Oscillation pattern, which has previously been suggested as a major driver of malaria.

METHODS

We used national malaria surveillance data from Guyana to estimate the correlation over time between the international gold price and reported P falciparum infections in individuals who were likely to be involved in mining activities (ie, men and boys aged between 15 and 50 years who were living in mining regions) for each month between 2007 and 2019. We compared the estimates with those obtained from individuals who were unlikely to be directly involved in mining activities (ie, women, children aged 12 years and younger, and adults aged over 70 years) and estimates obtained from individuals living in non-mining regions. We also evaluated the correlation between P falciparum infections and the El Niño-Southern Oscillation pattern in the same subpopulations and time period. Lastly, we evaluated the performance of a statistical model formulated to estimate P falciparum infections in real time using the international gold price as the predictor variable.

FINDINGS

The proportion of P falciparum malaria cases in temporary residents, which was used as a proxy for circulating individuals involved in gold mining, was highest during the years of peak gold price (ie, between 2008 and 2014). Cases of malaria in all demographic groups showed a strong positive correlation with the gold price, but only in regions with mining camps (0·88 [95% CI 0·84-0·89] for boys and men aged between 15 and 50 years and 0·80 [0·73-0·85] for the aggregated population of women, children aged 12 years and younger, and adults older than 70 years). The highest correlation occurred earlier in men and boys aged between 15 and 50 years, the demographic most likely to be miners, suggesting that transmission in mining camps is followed by infections in the community. On the basis of these findings, we were able to reliably forecast P falciparum malaria trends using only the gold price as the predictor variable. A 1% increase in gold price was associated with a 2·13% increase in P falciparum infections after 1 month in the mining populations, and with a 1·63% increase after 2 months in the non-mining populations. Lastly, La Niña climatic events showed an additional, smaller positive correlation with malaria transmission.

INTERPRETATION

Our analysis provides evidence that the P falciparum malaria surge observed in Guyana between 2008 and 2014 was likely to have been driven mainly by an increase in gold mining, while climate factors might have contributed synergistically. We propose that the international gold price over time is a useful indicator of malaria trends. We conclude that the feasibility of malaria elimination in Guyana, and in other areas in the Amazon where malaria and gold mining overlap, should be evaluated against the challenges posed by rapidly rising gold prices.

FUNDING

Ramón Areces Foundation, National Institutes of Health, and National Institute of General Medical Sciences.

Near-term climate change impacts on sub-national malaria transmission

The role of climate change on global malaria is often highlighted in World Health Organisation reports. We modelled a Zambian socio-environmental dataset from 2000 to 2016, against malaria trends and investigated the relationship of near-term environmental change with malaria incidence using Bayesian spatio-temporal, and negative binomial mixed regression models. We introduced the diurnal temperature range (DTR) as an alternative environmental measure to the widely used mean temperature. We found substantial sub-national near-term variations and significant associations with malaria incidence-trends. Significant spatio-temporal shifts in DTR/environmental predictors influenced malaria incidence-rates, even in areas with declining trends. We highlight the impact of seasonally sensitive DTR, especially in the first two quarters of the year and demonstrate how substantial investment in intervention programmes is negatively impacted by near-term climate change, most notably since 2010. We argue for targeted seasonally-sensitive malaria chemoprevention programmes.

Modeling the association between Aedes aegypti ovitrap egg counts, multi-scale remotely sensed environmental data and arboviral cases at Puntarenas, Costa Rica (2017–2018)

Problems with vector surveillance are a major barrier for the effective control of vector-borne disease transmission through Latin America. Here, we present results from a 80-week longitudinal study where Aedes aegypti (L.) (Diptera: Culicidae) ovitraps were monitored weekly at 92 locations in Puntarenas, a coastal city in Costa Rica with syndemic Zika, chikungunya and dengue transmission. We used separate models to investigate the association of either Ae. aegypti-borne arboviral cases or Ae. aegypti egg counts with remotely sensed environmental variables. We also evaluated whether Ae. aegypti-borne arboviral cases were associated with Ae. aegypti egg counts. Using cross-correlation and time series modeling, we found that arboviral cases were not significantly associated with Ae. aegypti egg counts. Through model selection we found that cases had a non-linear response to multi-scale (1-km and 30-m resolution) measurements of temperature standard deviation (SD) with a lag of up to 4 weeks, while simultaneously increasing with finely-grained NDVI (30-m resolution). Meanwhile, median ovitrap Ae. aegypti egg counts increased, and respectively decreased, with temperature SD (1-km resolution) and EVI (30-m resolution) with a lag of 6 weeks. A synchrony analysis showed that egg counts had a travelling wave pattern, with synchrony showing cyclic changes with distance, a pattern not observed in remotely sensed data with 30-m and 10-m resolution. Spatially, using generalized additive models, we found that eggs were more abundant at locations with higher temperatures and where EVI was leptokurtic during the study period. Our results suggest that, in Puntarenas, remotely sensed environmental variables are associated with both Ae. aegypti-borne arbovirus transmission and Ae. aegypti egg counts from ovitraps.

•

We sampled Ae. aegypti eggs using ovitraps for 80 weeks in Puntarenas, Costa Rica.

•

We were able to correlate Ae. aegypti egg-counts and arboviral cases with remotely sensed data.

•

Egg counts and arboviral cases were correlated with temperature and vegetation indices.

•

Arboviral cases were not associated with egg counts.

Plasmodium falciparum Genetic Diversity in Panamá Based on glurp, msp-1 and msp-2 Genes: Implications for Malaria Elimination in Mesoamerica

Panamá, together with all the nations in Mesoamerica, has committed to eliminate malaria from the region by 2020. As these countries approach malaria elimination and local transmission decreases, an active molecular surveillance to identify genotypes circulating along the border areas is particularly needed to accurately infer infection origin, drug resistance and disease propagation patterns in the region. This study evaluated the genetic diversity and allele frequencies of msp-1, msp-2 and glurp genes using different molecular analyses (nested PCR, PCR-restriction fragment length polymorphism (RFLP) and sequencing) from 106 autochthonous and imported P. falciparum isolates collected from different endemic areas in Panamá between 2003 and 2019. We also explored if P. falciparum genotypes assessed with these molecular markers were associated with relevant malaria epidemiological parameters using a multiple correspondence analysis. A strong association of certain local haplotypes with their geographic distribution in endemic areas, but also with parasite load and presence of gametocytes, was evidenced. Few multiclonal infections and low genetic diversity among locally transmitted P. falciparum samples were detected, consequent with the low transmission intensity of this parasite in Panamá, a pattern likely to be extended across Mesoamerica. In addition, several imported cases were genetically dissimilar to local infections and representative of more diverse extra-continental lineages.

Long-term transmission patterns and public health policies leading to malaria elimination in Panamá

Background

The present study provides a countrywide perspective of the malaria situation in Panamá over a long-term framework, with the purpose of identifying historical malaria resurgence events and their potential causes.

Methods

A descriptive-ecological study was conducted by analysing demographic and epidemiological annual malaria time series data in Panamá (1884–2019) using several data sources. Malaria intensity indicators were calculated during the study period. The effects of El Niño Southern Oscillation on malaria transmission were also analysed using a retrospective analysis of malaria cases between 1957 and 2019.

Results

Several factors were identified responsible for malaria resurgence in Panamá, mostly related with Malaria Control Programme weakening. During the past 20 years (2000–2019) malaria has progressively increased in prevalence within indigenous settlements, with a predominance of male cases and a high proportion (15% of total cases) in children less than 5 years old. During this period, a significant and increasing proportion of the Plasmodium falciparum cases were imported. Retrospective analysis (1957–2019) evidenced that ENSO had a significant impact on malaria transmission dynamics in Panamá.

Conclusions

Data analysis confirmed that although authorities have been successful in focalizing malaria transmission in the country, there are still neglected issues to be solved and important intercultural barriers that need to be addressed in order to achieve elimination of the disease by 2022. This information will be useful for targeting strategies by the National Malaria Elimination Programme.

Health policy impacts on malaria transmission in Costa Rica

Costa Rica is near malaria elimination. This achievement has followed shifts in malaria health policy. Here, we evaluate the impacts that different health policies have had on malaria transmission in Costa Rica from 1913 to 2018. We identified regime shifts and used regression models to measure the impact of different health policies on malaria transmission in Costa Rica using annual case records. We found that vector control and prophylactic treatments were associated with a 50% malaria case reduction in 1929-1931 compared with 1913-1928. DDT introduction in 1946 was associated with an increase in annual malaria case reduction from 7.6% (1942-1946) to 26.4% (1947-1952). The 2006 introduction of 7-day supervised chloroquine and primaquine treatments was the most effective health policy between 1957 and 2018, reducing annual malaria cases by 98% (2009-2018) when compared with 1957-1968. We also found that effective malaria reduction policies have been sensitive to natural catastrophes and extreme climatic events, both of which have increased malaria transmission in Costa Rica. Currently, outbreaks follow malaria importation into vulnerable areas of Costa Rica. This highlights the need to timely diagnose and treat malaria, while improving living standards, in the affected areas.

Rainfall Trends and Malaria Occurrences in Limpopo Province, South Africa

This contribution aims to investigate the influence of monthly total rainfall variations on malaria transmission in the Limpopo Province. For this purpose, monthly total rainfall was interpolated from daily rainfall data from weather stations. Annual and seasonal trends, as well as cross-correlation analyses, were performed on time series of monthly total rainfall and monthly malaria cases in five districts of Limpopo Province for the period of 1998 to 2017. The time series analysis indicated that an average of 629.5 mm of rainfall was received over the period of study. The rainfall has an annual variation of about 0.46%. Rainfall amount varied within the five districts, with the northeastern part receiving more rainfall. Spearman's correlation analysis indicated that the total monthly rainfall with one to two months lagged effect is significant in malaria transmission across all the districts. The strongest correlation was noticed in Vhembe (r = 0.54; p-value = <0.001), Mopani (r = 0.53; p-value = <0.001), Waterberg (r = 0.40; p-value =< 0.001), Capricorn (r = 0.37; p-value = <0.001) and lowest in Sekhukhune (r = 0.36; p-value = <0.001). Seasonally, the results indicated that about 68% variation in malaria cases in summer-December, January, and February (DJF)-can be explained by spring-September, October, and November (SON)-rainfall in Vhembe district. Both annual and seasonal analyses indicated that there is variation in the effect of rainfall on malaria across the districts and it is seasonally dependent. Understanding the dynamics of climatic variables annually and seasonally is essential in providing answers to malaria transmission among other factors, particularly with respect to the abrupt spikes of the disease in the province.

Study of the epidemiological behavior of malaria in the Darien Region, Panama. 2015-2017

BACKGROUND

Malaria is endemic in Darién and an assessment of the different factors affecting its epidemiology is crucial for the development of adequate strategies of surveillance, prevention, and disease control. The objective of this study was to determine the main characteristics of the epidemiological behavior of malaria in the Darien region.

METHODS

This research was comprised of a retrospective analysis to determine the incidence and malaria distribution in the Darien region from 2015 to 2017. We evaluated malaria indicators, disease distribution, incidence (by age group and sex), diagnostic methods, treatment, and control measures. In addition, we examined the cross-border migration activity and its possible contribution to the maintenance and distribution of malaria.

RESULTS

During the period of 2015-2017, we examined 41,141 thick blood smear samples, out of which 501 tested positive for malaria. Plasmodium vivax was responsible for 92.2% of those infections. Males comprised 62.7% of the total diagnosed cases. Meanwhile, a similar percentage, 62.7%, of the total cases were registered in economically active ages. The more frequent symptoms included fever (99.4%) and chills (97.4%), with 53.1% of cases registering between 2,000 and 6,000 parasites/μl of blood. The annual parasitic incidence (API) average was 3.0/1,000 inhabitants, while the slide positivity rate (SPR) was 1.2% and the annual blood examination rate (ABER) 22.5%. In Darién there is a constant internal and cross-border migration movement between Panama and Colombia. Malaria control measures consisted of the active and passive search of suspected cases and of the application of vector control measures.

CONCLUSION

This study provides an additional perspective on malaria epidemiology in Darién. Additional efforts are required to intensify malaria surveillance and to achieve an effective control, eventually moving closer to the objective of malaria elimination. At the same time, there is a need for more eco-epidemiological, entomological and migratory studies to determine how these factors contribute to the patterns of maintenance and dissemination of malaria.

The influence of weather and weather variability on mosquito abundance and infection with West Nile virus in Harris County, Texas, USA

Early warning systems for vector-borne diseases (VBDs) prediction are an ecological application where data from the interface of several environmental components can be used to predict future VBD transmission. In general, models for early warning systems only consider average environmental conditions ignoring variation in weather variables, despite the prediction from Schmalhausen's law about the importance of environmental variability for biological systems. We present results from a long-term mosquito surveillance program from Harris County, Texas, USA, where we use time series analysis techniques to study the abundance and West Nile virus (WNV) infection patterns in the local primary vector, Culex quinquefasciatus Say. We found that, as predicted by Schmalhausen's law, mosquito abundance was associated with the standard deviation and kurtosis of environmental variables. By contrast, WNV infection rates were associated with 8-month lagged temperature, suggesting environmental conditions during overwintering might be key for WNV amplification during summer outbreaks. Finally, model validation showed that seasonal autoregressive models successfully predicted mosquito WNV infection rates up to 2 months ahead, but did rather poorly at predicting mosquito abundance, a result that might reflect impacts of vector control for mosquito population reduction, geographic scale, and other artifacts generated by operational constraints of mosquito surveillance systems.

Socioeconomic and environmental factors associated with malaria hotspots in the Nanoro demographic surveillance area, Burkina Faso

BACKGROUND

With limited resources and spatio-temporal heterogeneity of malaria in developing countries, it is still difficult to assess the real impact of socioeconomic and environmental factors in order to set up targeted campaigns against malaria at an accurate scale. Our goal was to detect malaria hotspots in rural area and assess the extent to which household socioeconomic status and meteorological recordings may explain the occurrence and evolution of these hotspots.

METHODS

Data on malaria cases from 2010 to 2014 and on socioeconomic and meteorological factors were acquired from four health facilities within the Nanoro demographic surveillance area. Statistical cross correlation was used to quantify the temporal association between weekly malaria incidence and meteorological factors. Local spatial autocorrelation analysis was performed and restricted to each transmission period using Kulldorff's elliptic spatial scan statistic. Univariate and multivariable analysis were used to assess the principal socioeconomic and meteorological determinants of malaria hotspots using a Generalized Estimating Equation (GEE) approach.

RESULTS

Rainfall and temperature were positively and significantly associated with malaria incidence, with a lag time of 9 and 14 weeks, respectively. Spatial analysis showed a spatial autocorrelation of malaria incidence and significant hotspots which was relatively stable throughout the study period. Furthermore, low socioeconomic status households were strongly associated with malaria hotspots (aOR = 1.21, 95% confidence interval: 1.03-1.40).

CONCLUSION

These fine-scale findings highlight a relatively stable spatio-temporal pattern of malaria risk and indicate that social and environmental factors play an important role in malaria incidence. Integrating data on these factors into existing malaria struggle tools would help in the development of sustainable bottleneck strategies adapted to the local context for malaria control.

Variability in malaria cases and the association with rainfall and rivers water levels in Amazonas State, Brazil

Understanding the relations between rainfall and river water levels and malaria cases can provide important clues on modulation of the disease in the context of local climatic variability. In order to demonstrate how these relations can vary in the same endemic space, a coherence and wavelet phase analysis was performed between environmental and epidemiological variables from 2003 to 2010 for 8 municipalities (counties) in the state of Amazonas, Brazil (Barcelos, Borba, Canutama, Carauari, Coari, Eirunepé, Humaitá, and São Gabriel da Cachoeira). The results suggest significant coherences, mainly on the scale of annual variability, but scales of less than 1 year and of 2 years were also found. The analyses show that malaria cases display a peak at approximately 1 and a half months before or after peak rainfall and on average 1-4 months after peak river water levels in most of the municipalities studied. Each environmental variable displayed distinct local behavior in time and in space, suggesting that other local variables (e.g. topography) may control environmental conditions, favoring different patterns in each municipality. However, when the analyses were performed jointly it was possible to show a non-random order in these relations. Although environmental and climatic factors indicate a certain influence on malaria dynamics, surveillance, prevention, and control issues should not be overlooked, meaning that government public health interventions can mask possible relations with local hydrological and climatic conditions.

Population Dynamics of Anopheles albimanus (Diptera: Culicidae) at Ipetí-Guna, a Village in a Region Targeted for Malaria Elimination in Panamá

Anopheles albimanus Wiedemann is a major malaria vector in Mesoamerica and the Caribbean whose population dynamics, in response to changing environments, has been relatively poorly studied. Here, we present monthly adult and larvae data collected from May 2016 to December 2017 in Ipetí-Guna, a village within an area targeted for malaria elimination in the República de Panamá. During the study period we collected a total of 1678 Anopheles spp. mosquitoes (1602 adults and 76 larvae). Over 95% of the collected Anopheles spp. mosquitoes were An. albimanus. Using time series analysis techniques, we found that population dynamics of larvae and adults were not significantly correlated with each other at any time lag, though correlations were highest at one month lag between larvae and adults and four months lag between adults and larvae. Larvae population dynamics had cycles of three months and were sensitive to changes in temperature with 5 months lag, while adult abundance was correlated with itself (1 month lag) and with the normalized difference vegetation index (NDVI) with three months lag. A key observation from our study is the absence of both larvae and adults of An. albimanus between January and April from environments associated with Guna population’s daily activities, which suggests this time window could be the best time to implement elimination campaigns aimed at clearing Plasmodium spp. parasites from Guna populations using, for example, mass drug administration.

Exploring the Influence of Daily Climate Variables on Malaria Transmission and Abundance of Anopheles arabiensis over Nkomazi Local Municipality, Mpumalanga Province, South Africa

The recent resurgence of malaria incidence across epidemic regions in South Africa has been linked to climatic and environmental factors. An in-depth investigation of the impact of climate variability and mosquito abundance on malaria parasite incidence may therefore offer useful insight towards the control of this life-threatening disease. In this study, we investigate the influence of climatic factors on malaria transmission over Nkomazi Municipality. The variability and interconnectedness between the variables were analyzed using wavelet coherence analysis. Time-series analyses revealed that malaria cases significantly declined after the outbreak in early 2000, but with a slight increase from 2015. Furthermore, the wavelet coherence and time-lagged correlation analyses identified rainfall and abundance of Anopheles arabiensis as the major variables responsible for malaria transmission over the study region. The analysis further highlights a high malaria intensity with the variables from 1998-2002, 2004-2006, and 2010-2013 and a noticeable periodicity value of 256-512 days. Also, malaria transmission shows a time lag between one month and three months with respect to mosquito abundance and the different climatic variables. The findings from this study offer a better understanding of the importance of climatic factors on the transmission of malaria. The study further highlights the significant roles of An. arabiensis on malaria occurrence over Nkomazi. Implementing the mosquito model to predict mosquito abundance could provide more insight into malaria elimination or control in Africa.

Assessing changing weather and the El Niño Southern Oscillation impacts on cattle rabies outbreaks and mortality in Costa Rica (1985-2016)

BACKGROUND

Rabies is a major zoonotic disease affecting humans, domestic and wildlife mammals. Cattle are the most important domestic animals impacted by rabies virus in the New World, leading to thousands of cattle deaths per year and eliciting large economic losses. In the New World, virus transmission in cattle is primarily associated with Desmodus rotundus, the common vampire bat. This study analyses the association of weather fluctuations and the El Niño Southern Oscillation (ENSO), with the occurrence and magnitude, in terms of associated mortality, of cattle rabies outbreaks. Data from the 100 cattle rabies outbreaks recorded between 1985 and 2016 in Costa Rica were analyzed. Periodograms for time series of rabies outbreaks and the El Niño 4 index were estimated. Seasonality was studied using a seasonal boxplot. The association between epidemiological and climatic time series was studied via cross wavelet coherence analysis. Retrospective space-time scan cluster analyses were also performed. Finally, seasonal autoregressive time series models were fitted to study linear associations between monthly number of outbreaks, monthly mortality rates and the El Niño 4 index, temperature, and rainfall.

RESULTS

Large rabies mortality occurred towards the Atlantic basin of the country. Outbreak occurrence and size were not directly associated with ENSO, but were sensitive to weather variables impacted by ENSO. Both, ENSO phases and rabies outbreaks, showed a similar 5 year period in their oscillations. Cattle rabies mortality and outbreak occurrence increased with temperature, whereas outbreak occurrence decreased with rainfall. These results suggest that special weather conditions might favor the occurrence of cattle rabies outbreaks.

CONCLUSIONS

Further efforts are necessary to articulate the mechanisms underpinning the association between weather changes and cattle rabies outbreaks. One hypothesis is that exacerbation of cattle rabies outbreaks might be mediated by impacts of weather conditions on common vampire bat movement and access to food resources on its natural habitats. Further eco-epidemiological field studies could help to understand rabies virus transmission ecology, and to propose sound interventions to control this major veterinary public health problem.