Identifying the geographic distribution pattern of venomous snakes and regions of high snakebite risk in Iran

We describe species richness patterns of venomous snakes in Iran in order to produce snakebite risk prediction maps and identify gaps in regional health care centers capable of managing snakebites. We digitized distribution maps from the literature, Global Biodiversity Information Facility (GBIF), and the results of our own field studies of 24 terrestrial venomous snake species (including 4 endemic to Iran). Species richness patterns were associated with eight environmental factors. The variables have been extracted from the WorldClim dataset (bio12 = annual precipitation, bio15 = precipitation seasonality, bio17 = precipitation of the driest quarter, bio2 = mean diurnal range, bio3 = isothermality (bio2/bio7), bio4 = temperature seasonality, bio9 = mean temperature of the driest quarter and slope). Based on spatial analyses, species richness in Iran is highly affected by three environmental variables (bio12, 15, and 17) associated with precipitation. The relationship patterns among these predictors and species richness were strong and linear. The hotspot regions for venomous snakes species are concentrated on the western to southwestern and north to northeastern regions of Iran, which is partially consistent with the known Irano-Anatolian biodiversity hotspot. Because of the high number of endemic species and climatic conditions on the Iranian Plateau, the venoms of snakes distributed in those areas may contain novel properties and components.