Estimating snakebite incidence from mathematical models: A test in Costa Rica

Snakebite envenoming is a neglected public health challenge that affects mostly economically deprived communities who inhabit tropical regions. In these regions, snakebite incidence data is not always reliable, and access to health care is scare and heterogeneous. Thus, addressing the problem of snakebite effectively requires an understanding of how spatial heterogeneity in snakebite is associated with human demographics and snakes' distribution. Here, we use a mathematical model to address the determinants of spatial heterogeneity in snakebite and we estimate snakebite incidence in a tropical country such as Costa Rica.

We combined a mathematical model that follows the law of mass action, where the incidence is proportional to the exposed human population and the venomous snake population, with a spatiotemporal dataset of snakebite incidence (Data from year 1990 to 2007 for 193 districts) in Costa Rica. This country harbors one of the most dangerous venomous snakes, which is the Terciopelo (Bothrops asper, Garman, 1884). We estimated B. asper distribution using a maximum entropy algorithm, and its abundance was estimated based on field data. Then, the model was adjusted to the data using a lineal regression with the reported incidence. We found a significant positive correlation (R2 = 0.66, p-value < 0.01) between our estimation and the reported incidence, suggesting the model has a good performance in estimating snakebite incidence.

Our model underscores the importance of the synergistic effect of exposed population size and snake abundance on snakebite incidence. By combining information from venomous snakes' natural history with census data from rural populations, we were able to estimate snakebite incidence in Costa Rica. The model was able to fit the incidence data at fine administrative scale (district level), which is fundamental for the implementation and planning of management strategies oriented to reduce snakebite burden.