developing snakebite risk model using venomous snake ... · 06.04.2020 · 22 snakebite envenoming...
TRANSCRIPT
1
1 Title
2 Developing snakebite risk model using venomous snake habitat
3 suitability as an indicating factor: An application of species distribution
4 models in public health research
5 Authors:
6 Masoud Yousefi 1 ⃰, Anooshe Kafash 1, Ali Khani2, Nima Nabati3
7 Affiliations:
8 1 Department of Environmental Science, Faculty of Natural Resources, University of
9 Tehran, Iran
10 2 Khorasan-e-Razavi Provincial Office of the Department of the Environment, Mashhad,
11 Iran.
12 3Shahid Sadoughi Hospital, Shahid Sadoughi University of Medical Sciences, Yazd, Iran.
13 *Corresponding author:
14 Masoud Yousefi Email address: [email protected]
15
16
17
18
19
20
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
2
21 Abstract
22 Snakebite envenoming is an important public health problem in Iran, despite its risk not
23 being quantified. This study aims to use venomous snakes’ habitat suitability as an
24 indicator of snakebite risk, to identify high-priority areas for snakebite management
25 across the country. Thus, an ensemble approach using five distribution modeling
26 methods: Generalized Boosted Models, Generalized Additive Models, Maximum Entropy
27 Modeling Generalized Linear Models, and Random Forest was applied to produce a
28 spatial snakebite risk model for Iran. To achieve this, four venomous snakes’ habitat
29 suitability (Macrovipera lebetina, Echis carinatus, Pseudocerastes persicus and Naja
30 oxiana) were modeled and then multiplied. These medically important snakes are
31 responsible for the most snakebite incidents in Iran. Multiplying habitat suitability
32 models of the four snakes showed that the northeast of Iran (west of Khorasan-e-Razavi
33 province) has the highest snakebite risk in the country. In addition, villages that were at
34 risk of envenoming from the four snakes were identified. Results revealed that 51,112
35 villages are at risk of envenoming from M. lebetina, 30,339 from E. carinatus, 51,657
36 from P. persicus and 12,124 from N. oxiana. This paper demonstrates application of
37 species distribution modeling in public health research and identified potential snakebite
38 risk areas in Iran by using venomous snakes’ habitat suitability models as an indicating
39 factor. Results of this study can be used in snakebite and human–snake conflict
40 management in Iran. We recommend increasing public awareness of snakebite
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
3
41 envenoming and education of local people in areas which identified with the highest
42 snakebite risk.
43 Introduction
44 Snakebite envenoming is known as an important public health problem and the cause of
45 medical emergencies around the globe [1−14]. On Earth, between 421,000 and 1.2
46 million people are envenomed by venomous snakes annually and around 125,000 deaths
47 per year are attributable to snakebite envenoming [1, 6, 14, 15]. Snakebite envenoming is
48 mostly described as a neglected public health issue in the tropics, [2, 6, 16, 17] however,
49 it is also an important challenge for public health in temperate areas like Iran [18].
50 Iran is home to 81 snake species [19, 14] of which 25 are venomous (Nine species are sea
51 snakes and 16 species are terrestrial snakes). Macrovipera lebetina, Echis carinatus,
52 Pseudocerastes persicus and Naja oxiana are widespread in Iran and are responsible for
53 the most snakebite incidents in the country [18, 21]. A study reported 53,787 cases of
54 snake bites between 2002 and 2011 in Iran [18]. Despite considerable research into the
55 phylogeny, taxonomy, morphology and ecology of venomous snakes in Iran [22, 31]
56 snakebite envenoming has received less attention [18, 32]. In fact, snakebite is an
57 important uninvestigated public health problem and conservation challenge in Iran
58 [33−35]. Thus, more effort should be made to identify areas with high snakebite risk and
59 reduce envenoming risk from snakes.
60 Species distribution models (SDMs) have found an important application in
61 biodiversity research [36−38]. They are employed in studying habitat suitability [39−41],
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
4
62 identifying environmental drivers of species distribution [42−45] and predicting impacts
63 of climate change on biodiversity [46−51]. SDMs are successfully used to identify
64 suitable habitats of species even in areas with no distribution records [52−54]. Thus,
65 these models can be used to identify suitable habitats of venomous snakes as proxies of
66 snakebite risk [12, 55−57] in data poor regions like Iran.
67 The main goal of this paper was to apply SDMs and produce a spatial risk model
68 for snakebite in Iran. To do this, firstly, five distribution modeling methods [38]:
69 Generalized Boosted Models, Generalized Additive Models, Maximum Entropy
70 Modeling Generalized Linear Models, and Random Forest and distribution data of M.
71 lebetina, E. carinatus, P. persicus and N. oxiana were to produce their habitat suitability
72 models. Secondly, the five habitat suitability models [38−58] of each species were
73 combined by ensemble approach and finally the four species ensemble models were
74 multiplied to identify potential snakebite risk. We also in addition, number of villages
75 that are at risk of envenoming by these four snakes determined in Iran.
76
77 Materials and methods
78 Occurrence data
79 Distribution records of the M. lebetina, E. carinatus, P. persicus and N. oxiana (Fig 1)
80 were obtained from published books and papers [31, 59−66], and online databases (GBIF
81 and VertNet). These four snakes were selected because they are responsible for the most
82 snakebite incidents in Iran [18] and have the widest distribution range across the country
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
5
83 [19−20]. By combining presence records from the two sources 89 distribution records
84 were obtained for M. lebetina, 68 records for E. carinatus, 54 records for P. persicus and
85 37 records for N. oxiana (Fig 2).
86
87 Fig 1. Photos of the four medically important venomous snakes modeled in this
88 research to map snakebite in Iran. Saw-scaled Viper, Echis carinatus (a), Central
89 Asian Cobra, Naja oxiana (b), Levantine Viper, Macrovipera lebetina (c), and Persian
90 Horned Viper, Pseudocerastes persicus (d). Photos by Masoud Yousefi.
91
92 Fig 2. Distribution records of Macrovipera lebetina, Echis carinatus, Naja oxiana,
93 and Pseudocerastes persicus in a topographic overview of Iran. Provinces numbers
94 should read as follows; (1) West Azerbaijan, (2) East Azerbaijan, (3) Ardabil, (4)
95 Kurdistan, (5) Zanjan, (6) Gilan, (7) Kermanshah, (8) Hamedan, (9) Qazvin, (10) Alborz,
96 (11) Mazandaran, (12) Ilam, (13) Lorestan, (14) Markazi, (15) Qom, (16) Tehran, (17)
97 Khuzestan, (18) Chahar Mahaal and Bakhtiari, (19) Isfahan, (20) Semnan, (21) Golestan,
98 (22) Khorasan-e-Shomali (Northern Khorasan), (23) Khorasan-e-Razavi, (24) Yazd, (25)
99 Khorasan-e-Jonobi (Southern Khorasan), (26) Kohgiluyeh and Boyer-Ahmad, (27) Fars,
100 (28) Kerman, (29) Sistan and Baluchestan, (30) Bushehr, (31) Hormozgan.
101
102 Environmental data
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
6
103 Seven uncorrelated environmental variables (Table 1) related to climate, topography,
104 vegetation, and human footprint were used to develop the four snakes’ habitat suitably
105 models [12, 17, 28, 55, 67−70]. Climatic variables were downloaded from the WorldClim
106 database at 30-seconds spatial resolution [71]. Normalized Difference Vegetation Index
107 (NDVI) was considered as an indicator of resource availability [72]. Snakebite risk is
108 associated with human population density and activities [55], thus, human footprint index
109 was included in models [69]. Human footprint index was produced by combining data on
110 the extent of built environments, population density, electric infrastructure, crop lands,
111 pasture lands, roads, railways, and navigable waterways [70]. Topographic heterogeneity
112 was used as topography variable by measuring the standard deviation of elevation values
113 in area grid cells of 1km from a 90 m resolution in the Raster package [73] in the R
114 environment v.3.4.3 [74]. Elevation layer was obtained from the Shuttle Radar
115 Topography Mission (SRTM) elevation model [75].
116 Table 1. List of Environmental variables. Name, source and the variance inflation
117 factor (VIF) for the seven environmental variables which were used for developing
118 habitat suitability of venomous snakes in Iran.
Variable Description (abbreviation) References VIF
Climate Temperature seasonality (bio4) [71] 1.438
Annual precipitation (bio12) [71] 2.871
Precipitation seasonality (bio15) [71] 2.389
Precipitation of driest quarter (bio17) [71] 3.049
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
7
Topography Topographic heterogeneity (SD) [75] 1.315
Vegetation Normalized Difference Vegetation
Index (NDVI)
[72] 3.33
Anthropogenic Human Footprint (HFP) [69-70] 1.466
119
120 Snakebite risk modeling
121 An ensemble approach [38-58] was applied to model habitat suitability of the M.
122 lebetina, E. carinatus, P. persicus and N. oxiana, using five methods: Generalized
123 Boosted Models (GBM; [76]), Generalized Additive Models (GAM; [77]), Maximum
124 Entropy modeling (Maxent; [78]), Generalized Linear Models (GLM; [79]), Random
125 Forest (RF; [80]). Since these methods need background data points we generated a
126 randomly drawn sample of 10,000 background points (e.g., pseudo-absence points) from
127 the extent of the study area using the PresenceAbsence package [81]. Then areas
128 associated with high snakebite risk in Iran were identified by multiplying habitat
129 suitability models of the four species.
130 To quantify snakebite risk in Iran in more detail number of villages that are at risk
131 of envenoming from the four snakes determined and area of snakebite risk calculated in
132 the Raster package [73]. For this continuous habitat suitability models were converted to
133 suitable/unsuitable maps using maximum test sensitivity with a specificity threshold
134 [82]), and then overlaying 185,000 (population in these villages range from less than 50
135 to 5000 individuals) villages with each snake model.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
8
136 Model performance
137 Assessing model performance is an important step in distribution modeling studies [37].
138 Several metrics were introduced for estimating the predictive powers of models [37, 78,
139 83, 84]. Three metrics were considered to assess the produced habitat suitability models’
140 performances. The true skills statistic (TSS), area under the receiver operating
141 characteristic curve (AUC), and the Boyce index [37, 83, 85].
142 Analyses were carried out using the packages GISTools (https://rdrr.io/cran/
143 GISTools/), dismo (https://rdrr.io/cran/dismo/), biomod2 (https://cran.rproject.org/web/
144 packages/ biom od2/index.html), maptools (http://r-forge.r-project.org/projects/maptools/)
145 , SDMTools (htt ps://cran.r-project.org/web/packages/SDMTools/index.html), and ecospat
146 (https://cran.rproject.org/web/packages/ecospat/ecospat.pdf), in the R environment
147 (v.3.4.3).
148 Results
149 Model performance
150 All models developed in this study performed well based on the three model performance
151 evaluation metrics, AUC, TSS and Boyce index (See figs 3-6).
152 Habitat suitability
153 Echis carinatus
154 Based on ensemble model, southern part of Iran, north of Persian Gulf and vast areas in
155 central parts of the country are identified to have highest suitability for E. carinatus (Fig
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
9
156 3). Precipitation seasonality and NDVI and were the most important determinant of
157 habitat suitability of the species across the Iran.
158
159 Fig 3. Variables importance (a), models performance (b) and habitat suitability
160 model (c) of Echis carinatus in Iran.
161 Macrovipera lebetina
162 The most suitable habitats of M. lebetina are located in Zagros Mountains, Alborz
163 Mountains, Kopat-Dagh Mountains as well as in some isolated mountains in central Iran
164 (Fig 4). Precipitation seasonality, precipitation of the driest quarter, and human footprint
165 were the most important predictors of suitable habitats for the species.
166
167 Fig 4. Variables importance (a), models performance (b) and habitat suitability (c)
168 of Macrovipera lebetina in Iran.
169 Pseudocerastes persicus
170 Central, southwest and northeast of Iran have highest suitability for the P. persicus.
171 While, northern parts of the country are not suitable for this species (Fig 5). Results
172 showed that precipitation seasonality and NDVI and were the most important determinant
173 of habitat suitability of the species across the Iran.
174
175 Fig 5. Variables importance (a), models performance (b) and habitat suitability (c)
176 of Pseudocerastes persicus in Iran.
177 Naja oxiana
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
10
178 Naja oxiana’s most suitable habitats are located in north east of Iran around Kopat-Dagh
179 Mountains as well as some isolated patches in western parts of the country (Fig 6).
180 Precipitation of the driest quarter was the most important predictor of suitable habitats of
181 the species.
182
183 Fig 6. Variables importance (a), models performance (b) and habitat suitability (c)
184 of Naja oxiana in Iran.
185
186 Snakebite risk model
187 The four venomous snakes’ ensemble habitat suitability models were combined to
188 develop a snakebite risk model for Iran (Fig 7). Results showed that Khorasan-e-Razavi,
189 east of Semnan, north of Khorasan-e-Jonobi and south of Khorasan-e-Shomali provinces
190 have highest snakebite risk in Iran. West of Khorasan-e-Razavi province has high
191 suitability for the four venomous snakes.
192
193 Fig 7. Snakebite envenoming risk model in Iran. The snakebite risk model was
194 developed based on combined habitat suitability models’ of Macrovipera lebetina, Echis
195 carinatus, Pseudocerastes persicus and Naja oxiana.
196 Villages at risk of envenoming
197 Number of villages that are at the risk of envenoming by each of the four snakes (Table 2,
198 Fig 8) were determined. Results revealed that 51,112 villages are at risk of envenoming
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
11
199 from M. lebetina, 30,339 from E. carinatus, 51,657 from P. persicus and 12,124 from N.
200 oxiana. Area of envenoming risk by each species was estimated (Table 2), M. lebetina
201 and N. oxiana are identified with largest (362,558 km2) and smallest (121,803 km2) area,
202 respectively.
203
204 Fig 8. Villages at risk of envenoming from Macrovipera lebetina, Echis carinatus,
205 Pseudocerastes persicus and Naja oxiana. Villages at high risk are shown with red
206 circle and villages at moderate risk are shown with yellow circle.
207
208 Table 2. Area (km2) and number of villages that are at the snakebite risk from
209 Macrovipera lebetina, Echis carinatus, Pseudocerastes persicus and Naja oxiana in
210 Iran. Results are based on the ensemble models.
Species Villages at high risk Villages at moderate risk Area
Macrovipera lebetina 9,824 41,288 362,558
Echis carinatus 5,917 24,422 321,089
Pseudocerastes persicus 2,150 49,507 655,415
Naja oxiana 3,091 9,033 121,803
211
212 Discussion
213 With this research the first snakebite envenoming risk model was produced at fine
214 resolution (~1 km2) in Iran by modeling and multiplying habitat suitability of four
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
12
215 medically important venomous snakes which are responsible for the most snakebite
216 incidents in the country [18]. Northeastern parts of Iran were identified to have highest
217 snakebite risk in the country. Results showed that thousands of villages are located in
218 suitable ranges of the four venomous snakes, for example 51,657 villages are in the
219 suitable range of P. persicus. Thus these villages and villagers are at risk of envenoming
220 from these snakes. The snakebite risk model shows which parts of Iran are at risk of
221 envenoming from two, three or even four of snakes. All provinces of Iran, except those in
222 northwest of the country, are at risk of envenoming from at least two venomous snakes.
223 This highlights importance of snakebite envenoming as public health problem in Iran.
224 Under climate change some venomous species may expand their distribution
225 ranges [56-86], thus, envenoming risk will likely vary [12, 17, 87]. For instance, Nori et
226 al. [12] modeled distribution of five venomous snakes in Argentina for 2030 and 2080.
227 They found that the snakes’ suitable climate spaces will increase in human populated
228 areas of the country. In another study, Zacarias & Loyola [57] modeled current and future
229 distributions of 13 snakes in Mozambique and showed that venomous snake distribution
230 will change under climate change. They concluded patterns of snakebite risk may change
231 due to climatic changes [57]. Results of current research revealed that climatic variables
232 play an important role in shaping the distribution of four venomous snakes in Iran, thus
233 their distribution may alter with changing climate. It is predicted that suitable habitats of
234 Echis carinatus will increase in Iran [88]. This species is an important source of
235 snakebite in the country and its distributional range will likely increase under climate
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
13
236 change putting more populations and human settlements at risk of envenoming by this
237 snake until 2070.
238 Venomous snake populations are declining and many of them are listed by the
239 IUCN Red List as Vulnerable, Endangered or Critically Endangered [89]. Conservation
240 of snakes especially venomous snakes is a big challenge [89] as it is not easy to convince
241 people to conserve venomous snakes which are a significant cause of human mortality
242 and morbidity [4, 6, 14, 15]. It is necessary to identify areas with high risk of snakebite
243 envenoming and prioritize those areas for snakebite risk management in each country. In
244 this study potential risk areas were identified by using habitat suitability models as an
245 indicator of snakebite risk. Future snakebite monitoring in areas with high risk can show
246 the model performance and its power in predicting envenoming events across the
247 country.
248 Results of this study can be used to reduce snakebite and venomous snake
249 conflicts with local people, farmers and shepherds in the country. Snakebite risk can be
250 reduced through community education [90-92]. We encourage education of local people
251 about snakebite prevention measures in areas with highest snakebite risk. There are
252 simple solutions to prevent snakebite envenoming in areas with high risk, like villages in
253 west of Khorasan-e-Razavi province. For example, using bed nets and protecting feet,
254 ankles and lower legs by wearing boots can significantly reduce snakebite envenoming
255 [91-93]. We also suggest that areas where snakebite envenoming risk is high should be
256 monitored to determine envenoming events and villagers in these areas must always have
257 access to antivenom supplies [14].
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
14
258 Species distribution models are becoming important tools in public health research
259 [12, 94-103]. We encourage public hearth researchers to apply species distribution
260 models in developing snakebite risk map using venomous snakes’ habitat suitability as an
261 indicator, especially in data poor regions of the world [27, 52-54]. Our approach has the
262 potential for practical application in other countries whit high snakebite risk.
263 Acknowledgments
264 We thank Ollie Thomas for improving the English of the manuscript.
265 References
266 1. Chippaux JP. Snake-bites: Appraisal of the global situation. Bulletin of the World
267 Health Organization. 1998; 76: 515.9868843.
268 2. White J. Bites and stings from venomous animals: a global overview. Ther Drug
269 Monit. 2000; 22: 65–68.
270 3. Juckett G. Hancox JG. Venomous snakebites in the United States: Management
271 review and update. Am Fam Physician. 2002; 65: 1367–1378. PMID: 11996419.
272 4. Kasturiratne A, Pathmeswaran A, Fonseka MMD, Lalloo DG, Brooker S, De Silva
273 HJ. Estimates of disease burden due to land-snake bite in Sri Lankan hospitals.
274 Southeast. Asian J Trop Med Public Health. 2005; 36: 733. PMID: 16124448.
275 5. Gutiérrez JM. Theakston RDG, Warrell DA. Confronting the neglected problem of
276 snake bite envenoming: The need for a global partnership. PLoS Med. 2006; 3:
277 e150. https://doi.org/10.1371/journal.pmed.0030150.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
15
278 6. Kasturiratne A, Wickremasinghe AR, De Silva N, Gunawardena NK,
279 Pathmeswaran A, et al. The global burden of snakebite: A literature analysis and
280 modeling based on regional estimates of envenoming and deaths. PLoS Med. 2008;
281 5: 1591–1604. https://doi.org/10.1371/journal.pmed.0050218.
282 7. Cruz LS, Vargas R, Lopes AA. Snakebite envenomation and death in the developing
283 world. Ethn Dis. 2009; 19: 42. PMID: 19484874
284 8. Warrell D. Snake bite. The Lancet. 2010; 375: 77–88.
285 https://doi.org/10.1016/S0140-6736(09)61754-2.
286 9. Chippaux JP. Estimate of the burden of snakebites in sub-Saharan Africa: A meta-
287 analytic approach. Toxicon. 2011; 57: 586–599.
288 https://doi.org/10.1016/j.toxicon.2010.12.022.
289 10. Gutiérrez JM. Snakebite envenoming: a public health perspective. In Public health-
290 methodology, environmental and systems issues. InTech, 2012.
291 11. Hansson E, Sasa M, Mattisson K, Robles A, Gutiérrez JM. Using geographical
292 information systems to identify populations in need of improved accessibility to
293 antivenom treatment for snakebite envenoming in Costa Rica. PLoS Negl Trop Dis.
294 2013; 7: e2009. https://doi.org/10.1371/journal.pntd.0002009.
295 12. Nori J, Carrasco PA, Leynaud GC. Venomous snakes and climate change: ophidism
296 as a dynamic problem. Clim Change. 2014; 122: 67–80.
297 https://doi.org/10.1007/s10584-013-1019-6.
298 13. Chippaux JP. Incidence and mortality due to snakebite in the Americas. PLoS Negl
299 Trop Dis. 2017; 11: e0005662–39. https://doi.org/10.1371/journal.pntd.0005662.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
16
300 14. Longbottom J, Shearer FM, Devine M, Alcoba G, Chappuis F, Weiss DJ, et al.
301 Vulnerability to snakebite envenoming: a global mapping of hotspots. The Lancet
302 2018; 392: 673–684. https://doi.org/10.1016/S0140-6736(18)31224-8.
303 15. Gutiérrez JM, Calvete JJ, Habib AG, Harrison RA, Williams DJ, Warrell DA.
304 Snakebite envenoming. Nature. 2017; 3: 17063.
305 https://doi.org/10.1038/nrdp.2017.79.
306 16. Alirol E, Sharma SK, Bawaskar HS, Kuch U, Chappuis F. Snake Bite in South Asia:
307 A Review. PLoS Negl Trop Dis. 2010; 4: e603.
308 https://doi.org/10.1371/journal.pntd.0000603.
309 17. Yañez-Arenas C, Díaz-Gamboa L, Patrón-Rivero C, López-Reyes K, Chiappa-
310 Carrara X, Estimating geographic patterns of ophidism risk in Ecuador, Neotropical.
311 Biodiversity. 2018; 4: 55–61. https://doi.org/10.1080/23766808.2018.1454762.
312 18. Dehghani R, Fathi B, Panjeh Shahi M, Jazayeri M. Ten years of snakebites in Iran.
313 Toxicon. 2014; 90: 291–298. https://doi.org/10.1016/j.toxicon.2014.08.063.
314 19. Latifi M. The snakes of Iran. Tehran: Department of Environment; 2000.
315 20. Rajabizadeh M. Snakes of Iran. Tehran: IranShenasi Publishing; 2017.
316 21. Monzavi SM, Dadpour B, Afshari R. Snakebite management in Iran: Devising a
317 protocol. J Res Med Sci. 2014;19: 153–163. PMID: 24778670.
318 22. Rastegar-Pouyani E, Oraie H, Khosravani A, Kaboli M, Mobaraki A, Yousefi M, et
319 al. A re-evaluation of taxonomic status of Montivipera (Squamata: Viperidae) from
320 Iran using a DNA barcoding approach. Biochem Syst Ecol. 2014; 57: 350–356.
321 https://doi.org/10.1016/j.bse.2014.09.017.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
17
322 23. Moradi N, Rastegra-Pouyani N, Rastegra-Pouyani E. Geographic variation in the
323 morphology of Macrovipera lebetina (Linnaeus, 1758) (Ophidia: Viperidae) in
324 Iran. Acta Herpetol. 2014; 9: 187–202. https://doi.org/10.13128/Acta_Herpetol-
325 14384.
326 24. Fathinia B, Rastegar-Pouyani N, Rastegar-Pouyani E, Toodeh-Dehghan F,
327 Rajabizadeh M. Molecular systematics of the genus Pseudocerastes (Ophidia:
328 Viperidae) based on the mitochondrial cytochrome b gene. Turk J Zool. 2014; 38:
329 575–581. https://doi.org/10.3906/zoo-1308-25
330 25. Stümpel N, Rajabizadeh M, Avcı A, Wüster W, Joger U. Phylogeny and
331 diversification of mountain vipers (Montivipera, Nilson et al., 2001) triggered by
332 multiple Plio–Pleistocene refugia and high-mountain topography in the Near and
333 Middle East. Mol Phylogenet Evol. 2016; 101: 336–51.
334 https://doi.org/10.1016/j.ympev.2016.04.025.
335 26. Behrooz R, Kaboli M, Nourani E, Ahmadi M, Shabani AA, Yousefi M, et al. Habitat
336 modeling and conservation of the endemic Latifi's viper (Montivipera latifii) in Lar
337 National Park, Northern Iran. Herpetol Conserv Biol. 2015; 10: 572–582.
338 27. Yousefi M, Ahmadi M, Nourani E, Behrooz R, Rajabizadeh M, Geniez P, et al.
339 Upward altitudinal shifts in habitat suitability of mountain vipers since the last
340 glacial maximum. PloS One. 2015; 10: e0138087.
341 https://doi.org/10.1371/journal.pone.0138087.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
18
342 28. Khani Sh, Kami HG, Rajabizadeh M. Geographic variation of Gloydius halys
343 caucasicus (Serpentes: Viperidae) in Iran. Zool Middle East. 2017; 63: 303–310.
344 https://doi.org/10.1080/09397140.2017.1388491.
345 29. Fathinia B, Rastegar‐Pouyani N, Rastegar‐Pouyani E. Molecular phylogeny and
346 historical biogeography of genera Eristicophis and Pseudocerastes (Ophidia,
347 Viperidae). Zool Scr. 2018; 47: 673–685. https://doi.org/10.1111/zsc.12311.
348 30. Rastegar-Pouyani E, Oraie H, Khosravani A, Akbari A. Phylogenetic position of
349 Iranian pitvipers (Viperidae, Crotalinae, Gloydius) inferred from mitochondrial
350 cytochrome b sequences. Trop. Zool. 2018; 31: 55–67.
351 https://doi.org/10.1080/03946975.2018.1442288.
352 31. Oraie H, Rastegar-Pouyani E, Khosravani A, Moradi N, Akbari A, Sehhatisabet
353 ME. Molecular and morphological analyses have revealed a new species of blunt-
354 nosed viper of the genus Macrovipera in Iran. Salamandra. 2018; 54: 233–238.
355 http://doi.org/10.5281/zenodo.2586857.
356 32. Eslamian L, Mobaiyen H, Bayat-Makoo Z, Piri R, Benisi R, Naghavi-Behzad M.
357 Snake bite in Northwest Iran: A retrospective study. J Anal Res Clin Med. 2016; 4:
358 133–8.
359 33. Dehghani R, Rabani D, Panjeh-Shahi, M, Jazayeri M, Sabahi Bidgoli M. Incidence
360 of snake bite in Kashan, IR Iran during an eight-year period (2004-2011). Archives
361 of Trauma Research. 2012; 1: 67–71. http://doi.org/10.5812/ATR.6445.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
19
362 34. Farzaneh E, Fouladi N, Shafaee Y, Mirzamohammadi Z, Naslseraji F, Mehrpour O.
363 Epidemiological study of snakebites in Ardabil Province (Iran). Electron. Physician.
364 2017; 9: 3986–90. http://doi.org/10.19082/3986.
365 35. Ebrahimi V, Hamdami E, Khademian MH, Moemenbellah-Fard MD,
366 Vazirianzadeh B. Epidemiologic prediction of snake bites in tropical south Iran:
367 Using seasonal time series methods. Clin Epidemiol Glob. Health. 2018; 6: 208–
368 215. http://doi.org/10.15171/jarcm.2016.022.
369 36. Guisan A, Thuiller W. Predicting species distribution: offering more than simple
370 habitat models. Ecol Lett. 2005; 8: 993–1009. https://doi.org/10.1111/j.1461-
371 0248.2005.00792.x.
372 37. Guisan A, Thuiller W, Zimmermann NE. Habitat suitability and distribution
373 models: With Applications in R. Cambridge; Cambridge University Press: 2017.
374 38. Araújo MB, Anderson RP, Márcia Barbosa A, Beale CM, Dorman CF, Early R, et
375 al. Standards for distribution models in biodiversity assessments. Sci Adv. 2019; 5:
376 eaat4858. https://doi.org/10.1126/sciadv.aat4858
377 39. Milanesi P, Holderegger R, Caniglia R, Fabbri E, Randi E. Different habitat
378 suitability models yield different least‐cost path distances for landscape genetic
379 analysis. Basic Appl Ecol. 2015; 17: 61–
380 71. https://doi.org/10.1016/j.baae.2015.08.008.
381 40. Moradi S, Sheykhi Ilanloo S, Kafash A, Yousefi Y. Identifying high-priority
382 conservation areas for avian biodiversity using species distribution modeling. Ecol
383 Indic. 2019; 97: 159–164. https://doi.org/10.1016/j.ecolind.2018.10.003.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
20
384 41. Farrell A, Wang G, Rush SA, Martin JA, Belant JL, Butler AB, et al. Machine
385 learning of large-scale spatial distributions of wild turkeys with high-dimensional
386 environmental data. Ecol Evol. 2019; 9: 5938–5949.
387 https://doi.org/10.1002/ece3.5177.
388 42. Kalle R, Ramesh T, Qureshi Q, Sankar K. Predicting the distribution pattern of small
389 carnivores in response to environmental factors in the Western Ghats. PloS
390 One 2013; 8: e79295. https://doi.org/10.1371/journal.pone.0079295.
391 43. Bradie J, Leung B. A quantitative synthesis of the importance of variables used in
392 MaxEnt species distribution models. J Biogeogr. 2017; 44: 1344–1361.
393 https://doi.org/10.1111/jbi.12894.
394 44. Ashoori A, Kafash A, Varasteh Moradi H, Yousefi M, Kamyab H, Behdarvand N,
395 et al. Habitat modeling of the Common Pheasant Phasianus colchicus (Galliformes:
396 Phasianidae) in a highly modified landscape: application of species distribution
397 models in the study of a poorly documented bird in Iran. Europ Zool J. 2018; 85:
398 373–381. https://doi.org/10.1080/24750263.2018.1510994.
399 45. Sutton LJ, Puschendorf R. Climatic niche of the Saker Falcon Falco cherrug:
400 predicted new areas to direct population surveys in Central Asia. Ibis. 2020; 162:
401 27–41. https://doi.org/10.1111/ibi.12700.
402 46. Qin A, Liu B, Guo Q, Bussmann RW, Ma F, Jian Z, et al. Maxent modeling for
403 predicting impacts of climate change on the potential distribution of Thuja
404 sutchuenensis Franch., an extremely endangered conifer from southwestern China.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
21
405 Glob Ecol Conserv. 2017; 10: 139–146.
406 https://doi.org/10.1016/j.gecco.2017.02.004.
407 47. Wang R, Li Q, He S, Liu Y, Wang M, Jiang G. Modeling and mapping the current
408 and future distribution of Pseudomonas syringae pv. actinidiae under climate
409 change in China. PLoS ONE. 2018; 13: e0192153.
410 https://doi.org/10.1371/journal.pone.0192153.
411 48. Kafash A, Ashrafi S, Ohler A, Yousefi M, Malakoutikhah Sh, Koehler G, et al.
412 Climate change produces winners and losers: Differential responses of amphibians
413 in mountain forests of the Near East. Glob Ecol Conserv. 2018; 16: e00471.
414 https://doi.org/10.1016/j.gecco.2018.e00471.
415 49. Schivo F, Bauni V, Krug P, Quintana RD. Distribution and richness of amphibians
416 under different climate change scenarios in a subtropical region of South America.
417 Appl Geogr. 2019; 103: 70–89. https://doi.org/10.1016/j.apgeog.2019.01.003.
418 50. Ashrafzadeh MR, Naghipour AA, Haidarian M, Kusza S, Pilliod DS. Effects of
419 climate change on habitat and connectivity for populations of a vulnerable, endemic
420 salamander in Iran. Glob Ecol Conserv. 2019; 19: e00637.
421 https://doi.org/10.1016/j.gecco.2019.e00637.
422 51. Yousefi M, Jouladeh-Rodbar A, Kafash A. Using endemic freshwater fishes as
423 proxies of their ecosystems to identify high priority rivers for conservation under
424 climate change. Ecol Indic. 2020; 112: 106137.
425 https://doi.org/10.1016/j.ecolind.2020.106137.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
22
426 52. Segura S, Coppens d'Eeckenbrugge G, López L, Grum M, Guarino L. Mapping the
427 potential distribution of five species of Passiflora in Andean countries. Genet
428 Resour Crop Evol. 2003; 50: 555–566. https://doi.org/10.1023/A:1024424012134.
429 53. Rebelo H, Jones G. Ground validation of presence‐only modelling with rare species:
430 a case study on barbastelles Barbastella barbastellus (Chiroptera:
431 Vespertilionidae). J Appl Ecol. 2010; 47: 410–420. https://doi.org/10.1111/j.1365-
432 2664.2009.01765.x.
433 54. Rhoden CM, Peterman WE, Taylor CA. Maxent-directed field surveys identify new
434 populations of narrowly endemic habitat specialists. PeerJ. 2017; 5: e3632.
435 55. Yañez-Arenas C, Peterson AT, Mokondoko P, Rojas-Soto O, Martínez-Meyer E.
436 The use of ecological niche modeling to infer potential risk areas of snakebite in the
437 Mexican State of Veracruz. PLoS One. 2014; 9: e100957.
438 https://doi.org/10.1371/journal.pone.0100957.
439 56. Yañez-Arenas C, Peterson AT, Rodriguez-Medina K, Barve, N. Mapping current
440 and future potential snakebite risk in the new world. Clim Change. 2016; 134: 697–
441 711. https://doi.org/10.1007/s10584-015-1544-6.
442 57. Zacarias D, Loyola R. Climate change impacts on the distribution of venomous
443 snakes and snakebite risk in Mozambique. Clim Change. 2019; 152: 195.
444 https://doi.org/10.1007/s10584-018-2338-4.
445 58. Araújo MB, New M. Ensemble forecasting of species distributions. Trends Ecol
446 Evol. 2007; 22: 42–47. https://doi.org/10.1016/j.tree.2006.09.010.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
23
447 59. Moradi N, Shafiei S, Sehhatisabet ME. The snake fauna of Khabr National Park,
448 southeast of Iran. IJAB. 2013; 9: 41–55.
449 60. Hosseinian Yousefkhani SS, Yousefi M, Khani A, Rastegar Pouyani E. Snake fauna
450 of Shirahmad wildlife refuge and Parvand protected area, Khorasan Razavi
451 province, Iran. Herpetol Notes. 2014; 7: 75–82.
452 61. Sabbaghzadeh A, Mashayekhi M. Survey of reptile fauna of Nazmabad of Arak,
453 Markazi Province, Iran. AJBAS. 2015; 7: 101–108.
454 62. Yousefi M, Khani A, Shaykhi Ilanloo S, Kafash A. Rastegar Pouyani E. Reptile
455 fauna of the Khajeh protected area, with assessing its similarities with
456 physiogeographical area of the Iranian Lizards. TBJ. 2015; 22: 13–22.
457 63. Mozafari O, Kamali, K, Fahimi H. Atlas of reptiles of Iran. Tehran: Department of
458 the Environment; 2016.
459 64. Kami HG, Babaei Savasari R. new records of the Turkmenian Fat-Tailed Gecko,
460 Eublepharis turcmenicus Darevsky, 1978, from Khorasan Razavi Province, Iran
461 (Squamata: Eublepharidae). IJAB. 2017; 13: 263–273.
462 65. Fathinia B, Rastegar-Pouyani E, Rastegar-Pouyani N, Darvishnia H. A new species
463 of the genus Rhynchocalamus Günther, 1864 (Reptilia: Squamata: Colubridae) from
464 Ilam Province in western Iran. Zootaxa. 2017; 4282: 473–486.
465 http://dx.doi.org/10.11646/zootaxa.4282.3.3
466 66. Salemi A, Heydari N, Jahan Mahin M. A new distribution record for the rare
467 Maynard’s Longnose Sand Snake, Lytorhynchus maynardi Alcock and Finn, 1896
468 from Nikshahr, southeastern Iran. Herpetol Notes. 2018; 11: 617–619.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
24
469 67. Pettorelli N, Vik JO, Mysterud A, Gaillard JM, Tucker CJ, Stenseth NC. Using the
470 satellite-derived NDVI to assess ecological responses to environmental change.
471 Trends Ecol Evol. 2005; 20: 503–510. https://doi.org/10.1016/j.tree.2005.05.011.
472 68. Terribile LC, Olalla-Tárraga MÁ, Morales-Castilla I. Rueda M, Vidanes RM,
473 Rodríguez MÁ, et al. Global richness patterns of venomous snakes reveal
474 contrasting influences of ecology and history in two different clades. Oecologia.
475 2009; 159: 617. https://doi.org/10.1007/s00442-008-1244-2.
476 69. Venter O, Sanderson EW, Magrach A, Allan JR, Beher J, Jones KR, et al. Sixteen
477 years of change in the global terrestrial human footprint and implications for
478 biodiversity conservation. Nat Commun. 2016; 7: 12558.
479 https://doi.org/10.1038/ncomms12558.
480 70. Venter O, Sanderson EW, Magrach A, Allan JR, Beher J, Jones KR, et al. Global
481 terrestrial Human Footprint maps for 1993 and 2009. Sci Data. 2016; 3: 160067.
482 https://doi.org/10.1038/sdata.2016.67.
483 71. Fick SE, Hijmans RJ. WorldClim 2: new 1-km spatial resolution climate surfaces
484 for global land areas. Int J Climatol. 2017; 37: 4302–4315.
485 https://doi.org/10.1002/joc.5086.
486 72. Broxton PD, Zeng X, Schefic W, Troch PA. A MODIS-Based global 1-km
487 maximum green vegetation fraction dataset. J Appl Meteorol Clim. 2014; 53:
488 1996–2004. https://doi.org/10.1175/JAMC-D-13-0356.1.
489 73. Hijmans RJ. Raster: geographic data analysis and modeling. R package, 2015.
490 Available from: https://cran.r-project.org/web/packages/raster/index.html.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
25
491 74. R Core Team R: A language and environment for statistical computing. Vienna,
492 Austria: R Foundation for Statistical Computing; 2014.
493 75. Jarvis A, Reuter HI, Nelson A, Guevara E. Hole-filled SRTM for the globe
494 Version 4; 2008. Available from: the CGIAR-CSI SRTM 90m Database
495 (http://srtm.csi. cgiar.org).
496 76. Ridgeway G. The state of boosting. Comput. Stat. 1999; 31: 172–181.
497 77. Hastie TJ, Tibshirani R. Generalized Additive Models. London: Chapman and
498 Hall; 1990.
499 78. Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species
500 geographic distributions. Ecol Model. 2006; 190: 231–259.
501 https://doi.org/10.1016/j.ecolmodel.2005.03.026.
502 79. McCullagh P, Nelder JA. Generalized Linear Models. London: Chapman and
503 Hall; 1989.
504 80. Breiman L. Random forests. Mach Learn. 2001; 45: 5–32.
505 81. Freeman EA, Moisen G. PresenceAbsence: an r package for presence-absence
506 model analysis. J Stat Softw. 2008; 23:1–31. https://doi.org/10.18637/jss.v023.i11
507 82. Liu C, White M, Newell G. Selecting thresholds for the prediction of species
508 occurrence with presence-only data. J Biogeogr. 2013; 40: 778–789.
509 https://doi.org/10.1111/jbi.12058.
510 83. Boyce MS, Vernier PR, Nielsen SE, Schmiegelow FKA. Evaluating resource
511 selection functions. Ecol Modell. 2002; 157: 281–300.
512 https://doi.org/10.1016/S0304-3800(02)00200-4.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
26
513 84. Hirzel, A.H., Le Lay, G., Helfer, V., Randin, C., Guisan, A., 2006. Evaluating the
514 ability of habitat suitability models to predict species presences. Ecol. Modell. 199,
515 142–152. https://doi.org/10.1016/j.ecolmodel.2006.05.017.
516 85. Fielding AH. Bell JF. A review of methods for the assessment of prediction errors
517 in conservation presence–absence models. Environ Conserv. 1997; 24: 38–49.
518 https://doi.org/10.1017/S0376892997000088.
519 86. Yousefi M, Kafash A, Valizadegan N, Sheykhi Ilanloo S, Rajabizadeh M,
520 Malekoutikhah S, et al. Climate change is a major problem for biodiversity
521 conservation: A systematic review of recent studies in Iran. Contemp Probl Ecol.
522 2019; 12: 394–403. https://doi.org/10.1134/S1995425519040127.
523 87. Hannah L. Climate change biology (Second Edition). Cambridge: Academic Press;
524 2015.
525 88. Yousefi M, Kafash A. Venomous snakes of Iran under climate change. SCCS
526 Europe– 4th Hungarian Student Conference on Conservation Science, 4–8
527 September, 2018.
528 89. Maritz B, Penner J, Martins M, Crnobrnja-Isailović J, Spear S, Alencar LRV, et al.
529 Identifying global priorities for the conservation of vipers. Biol Conserv. 2016; 204:
530 94–102. https://doi.org/10.1016/j.biocon.2016.05.004.
531 90. Sharma SK, Bovier P, Jha N, Alirol E, Loutan L, Chappuis F. Effectiveness of rapid
532 transport of victims and community health education on snake bite fatalities in rural
533 Nepal. Am J Trop Med Hyg. 2013; 89: 145–150. https://doi.org/10.4269/ajtmh.12-
534 0750.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
27
535 91. WHO. Guidelines for the management of snakebites. Geneva: World Health
536 Organization; 2016.
537 92. Trogridou A, Graham S, Santra V, Owens JB, Bharti O, Malhotra A. Prevention is
538 better than cure: snakebite education in India. Episthmes Agogis; 2018.
539 93. Chappuis F, Sharma SK, Jha N, Loutan L, Bovier PA. Protection against snake bites
540 by sleeping under a bed net in southeastern Nepal. Am J Trop Med Hyg. 2007; 77:
541 197–9. DOI: https://doi.org/10.4269/ajtmh.2007.77.197.
542 94. Khatchikian C, Sangermano F, Kendell D, Livdahl T. Evaluation of species
543 distribution model algorithms for fine-scale container-breeding mosquito risk
544 prediction. Med Vet Entomol. 2011; 25: 268–275. 10.1111/j.1365-
545 2915.2010.00935.x.
546 95. Escobar LE, Peterson AT, Favi M, Yung V, Pons DJ, Medina-Vogel G. Ecology
547 and geography of transmission of two bat-borne rabies lineages in Chile. PLoS Negl
548 Trop Dis. 2013; 7: e2577. https://doi.org/10.1371/journal.pntd.0002577.
549 96. Du Z, Wang Z, Liu Y, Wang H, Xue F, Liu Y. Ecological niche modeling for
550 predicting the potential risk areas of severe fever with thrombocytopenia syndrome.
551 Int J Infect Dis. 2014; 26: 1–8. https://doi.org/10.1016/j.ijid.2014.04.006.
552 97. Peterson T., Mapping disease transmission risk: Enriching models using
553 biogeography and ecology. Baltimore: Johns Hopkins University Press; 2014.
554 98. Escobar LE, Craft ME. Advances and limitations of disease biogeography using
555 ecological niche modeling. Front Microbiol. 2016; 7: 1174.
556 https://doi.org/10.3389/fmicb.2016.01174.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
28
557 99. Escobar LE, Peterson AT, Papeş M, Favi M, Yung V, Restif O, et al. Ecological
558 approaches in veterinary epidemiology: mapping the risk of bat-borne rabies using
559 vegetation indices and night-time light satellite imagery. Vet. Res. 2015; 46: 92.
560 https://doi.org/10.1186/s13567-015-0235-7.
561 100. Nyakarahuka L, Ayebare S, Mosomtai G, Kankya C, Lutwama J, Mwiine FN, et al.
562 Ecological niche modeling for filoviruses: A risk map for Ebola and Marburg virus
563 disease outbreaks in Uganda. PLOS Currents Outbreaks. 2017 Sep 5. Edition 1
564 2017.https://doi.org/10.1371/currents.outbreaks.07992a87522e1f229c7cb023270a
565 2af1.
566 101. Soucy JPR, Slatculescu AM, Nyiraneza C, Ogden NH, Leighton PA, Kerr
567 JT, et al. High-resolution ecological niche modeling of Ixodes scapularis ticks
568 based on passive surveillance data at the Northern Frontier of Lyme disease
569 emergence in North America. Vector Borne Zoonotic Dis. 2018; 18: 235–42.
570 https://doi.org/10.1089/vbz.2017.2234
571 102. Little EAH, Anderson JF, Stafford KC, Eisen L, Eise RJ, Molaei G. Predicting
572 spatiotemporal patterns of Lyme disease incidence from passively collected
573 surveillance data for Borrelia burgdorferi sensu lato-infected Ixodes scapularis
574 ticks. Ticks. Ticks Tick Borne Dis. 2019; 10: 970–980.
575 https://doi.org/10.1016/j.ttbdis.2019.04.010.
576 103. Johnson EE, Escobar LE, Zambrana-Torrelio C. An Ecological Framework
577 for Modeling the Geography of Disease Transmission. Trends Ecol Evol. 2019; 34:
578 655–668. https://doi.org/10.1016/j.tree.2019.03.004.
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint
.CC-BY 4.0 International licenseavailable under awas not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made
The copyright holder for this preprint (whichthis version posted April 6, 2020. ; https://doi.org/10.1101/2020.04.06.027342doi: bioRxiv preprint