J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 317 http://jad.tums.ac.ir Published Online: January 06, 2016 Original Article Impact of Climate and Environmental Factors on West Nile Virus Circulation in Iran *Farzaneh Ahmadnejad 1, 2, Vahid Otarod 3, Amanollah Fathnia 4, Ali Ahmadabadi 5, Mohammad H. Fallah 6, Alireza Zavareh 1, Nargess Miandehi 1, Benoit Durand 7, Philippe Sabatier 2 1Viral Vaccines Production Department, Pasteur Institute, Tehran, Iran 2TIMC-IMAG Team EPSP, VetAgroSup, Campus Vétérinaire de Lyon, France 3Quarantine and Biosafety Directorate General, Iran Veterinary Organization, Tehran, Iran 4Department of Geography, Razi University, Kermanshah, Iran 5Faculty of Geography, Kharazmi University, Tehran, Iran 6Department of Poultry Viral Diseases, Razi Vaccine and Serum Research Institute, Alborz, Iran 7University Paris Est, Anses, Laboratory of Animal Health, Epidemiology unit, Maisons-Alfort, France (Received 31 May 2014; accepted 17 Dec 2014) Abstract Background: Geographic distribution of West Nile virus (WNV) is heterogeneous in Iran by a high circulation in the southern-western areas. The objective of our study was to determine environmental and climatic factors associ- ated with the risk of WNV equine seropositivity in Iran. Methods: Serological data were obtained from a serosurvey conducted in equine population in 260 districts in Iran. The climate and environmental parameters included in the models were distance to the nearest wetland area, type of stable, Normalized Difference Vegetation Index (NDVI), annual mean temperature, humidity and precipitation. Results: The important risk factors included annual mean temperature, distance to wetlands, local and seasonal NDVI differences. The effect of local NDVI differences in spring was particularly notable. This was a normalized difference of average NDVI between two areas: a 5 km radius area centered on the stable and the 5–10 km sur- rounding area. Conclusion: The model indicated that local NDVI’s contrast during spring is a major risk factor of the transmission of West-Nile virus in Iran. This so-called oasis effect consistent with the seasonal production of vegetation in spring, and is associated to the attractiveness of the local NDVI environment for WNV vectors and hosts. Keywords: West Nile Virus, Environment, Climate, NDVI, Wetland Introduction West Nile virus (WNV) is not restricted by international borders and involves all continents, revealing adaptation of the virus with different ecological zones. Mosquitoes and birds, as vector and reservoir respec- tively, play an important role in the virus circulation. Several studies have defined the potential importance of environmental fac- tors on the WNV-vector-host interaction, the virus replication and on the population/ habitat of the vector/reservoir (Bolling et al. 2005, Kinney et al. 2006). Climate variation, as an environmental factor, has two different effects on ecologi- cal processes: (i) direct effects on virus/ hosts/vectors, for example, by selecting the variants sensitive, or controlling the physiol- ogy of organism such as metabolic processes and reproduction (Gubler et al. 2001, Reiter 2001), and (ii) indirect effect on the habitats, the ecosystems and the relationships among the organisms (Stenseth et al. 2002). Effect of climatic factors on human WNV infection has been studied by analyzing a spectrum of *Corresponding authors: Dr Farzaneh Ahmadnejad, E-mail: F_ahmadnejad@yahoo.com J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 318 http://jad.tums.ac.ir Published Online: January 06, 2016 weather factors for different states in the United States, and demonstrated a relation between these factors and incidence of human WNV cases (Soverow et al. 2009). Climate variations have changed the geo- graphic and seasonal patterns of WNV with shifting toward new geographic areas and occurring earlier in the transmission season. Gibbs et al. demonstrated that temperature, housing density, urban/suburban land use, and physiographic region are important var- iables associated with the geographic dis- tributions of WNV in Georgia. In addition, the models according climate and landscape change were used to predict the future expansion of WNV in North America (Gibbs et al. 2006, Chen et al. 2013). Environmental factors such as elevation range, physiographic region, and percentage of vegetation cover have been evaluated to determine human WNV risk in Chicago (Ruiz et al. 2004). In western Kenya, mosquito larval presence was associated with lower elevations, greater wetness, short distances to water, and land use (Bian et al. 2006). Among climatic factors, temperature is the most important extrinsic factor affecting the population and dynamic of virus/vector/ hosts. Humidity is another factor, which has a direct effect on the survival of mosquitoes, by increasing vector flight and host-seeking behaviors. Survival rate might be reduced when hot weather is accompanied by low humidity. Rain provides the breeding sites for mosquitoes and helps create a humid en- vironment, which prolongs the life of vec- tors. However, response to precipitation de- pends on the geographic location, season and the mosquito species (Reiter 2001, Ruiz et al. 2010). Normalized difference vegetation index (NDVI) is a simple numerical indicator that can be used to assess vegetation density. As plants respond quickly to climate variations, NDVI is related to the climatic variations such as increase or decrease of precipitation, temperature and sunshine level (Stow et al. 2004).The first evidence of WN virus circu- lation in Iran, a country with climatic and environmental diversity, goes back in 1970s when existence of WNV antibodies among healthy human population has been demon- strated (Saidi S. 1976). After a several dec- ades gap, in a study carried out among 1054 animal sera collected from different prov- inces, in 2008–2009, a relatively high circu- lation of the virus has been found in the country. The study showed a highly hetero- geneous prevalence in different provinces, so that a higher prevalence was detected in the southwestern areas while lower preva- lence was in the northern areas (Ahmadnejad et al. 2011). According to the last checklist, Iranian mosquitoes include 64 species and 3 subspecies classified in 7 genera and 12 un- verified species. The most important vector of WNV, Culex spp, includes 19 species (Azari-Hamidian 2007). Although, there are many studies on Iranian mosquito fauna, there is no information about WNV vectors in the country. This study addressed the dynamic and the spatial patterns of WNV transmission by an- alyzing the geographical and ecology pa- rameters affect the WNV circulation in the infected areas in Iran. Materials and Methods Study area Iran is bordered by the Gulf of Oman, the Persian Gulf, and the Caspian Sea. It has arid or semiarid climates mostly character- ized by low rainfall and high potential evap- otranspiration and its location cause it to re- ceive less than a third of the world average precipitation. The complex physical condi- tions of Iran including topography and land- scape have created a diverse climate pattern, so that it led to the formation of different ecological zones with various species of plants and animals (Heshmati 2007). J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 319 http://jad.tums.ac.ir Published Online: January 06, 2016 Furthermore, the climate is influenced by Caspian Sea in north, coastal areas of south of the country, Mediterranean area and Red sea. Alborz and Zagros mountains also play an important role in determining the non- uniform spatial and temporal distribution of precipitation in the whole country. Iran is located at high latitudes rendering vast differences between seasons in the country. Temperature in Iran is intense func- tion of altitude, latitude and moisture content of the atmosphere. The summer is extremely hot with temperatures in the interior; over 55 °C has been recorded in some places. Hu- midity prevents temperature fluctuations in the north and south of the country (Alijani 1995). The rainy period in most of the country is from November to May followed by dry pe- riod between May and October with rare precipitation. The average annual rainfall of the country is about 240 mm with maximum amounts in the Caspian Sea plains. Collection and processing of data Serological data: Equine data and their geographic coordinates were obtained from a previous seroprevalence study. The sera (n= 1054) were collected from 260 stables lo- cated within 27 of 30 provinces in the coun- try. Diagnosis of WNV antibodies in the sera was performed by using plaque reduction neutralization test (PRNT) (Fig. 1). Antibodies to WNV were detected in 249 (23.7%) of the samples (Ahmadnejad et al. 2011). Landscape and Climatic data Climatic data: We analyzed three climatic factors such as rainfall, temperature and rel- ative humidity, which are very important factors to affect the mosquito breeding activ- ities. Mean monthly temperature, precipitation and relative humidity data of 10-year period (1996–2005) from 103 stations were ob- tained from Iran Meteorological Organiza- tion (IRIMO). Numerical layers of sampling places and meteo-stations were inputted into ArcMap 9.3 software. The stations and their nearest sampling places were defined by Thiessen method (Fortes et al. 2005). Mean yearly data was used for analysis between WNV prevalence and climatic parameters. Average temperature and precipitation is given in degrees Celicius and millimeter, respectively. Relative humidity is expressed as a percentage. Normalized Difference Vegetation Index (NDVI) is a simple numerical indicator that can be used to analyze remote sensing meas- urements and assess vegetation density and whether the target being observed contains live green vegetation or not (Brown et al. 2008). NDVI time series data were acquired by Advanced Very High Resolution Radiomete (AVHRR) sensor onboard the National Oce- anic and Atmospheric Administrationr satel- lite (NOAA). The images have been extracted from http://www.class.ngdc.noaa.gov based on AVHRR level 1B and resolution of 1.1 km. As satellite images cannot be directly com- bined with other geographic information system, so geometric correction is necessary to bring the images from ground or slant range geometry into a map reference, and for this purpose, we used ENVI 4.3 software. Composites were used to calculate NDVI values using ERDAS 8.7 software and using near infrared (NIR) and red bands as fol- lows: dNIR dNIR NDVI Re Re    NDVI raster files were imported into a Geographic Information System (GIS, ESRI1, ArcGISTM 9.3). Buffers of 5 and 5– 10 km around each sampled point were cre- ated with the GIS software and mean NDVI was extracted within each circular buffer. J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 320 http://jad.tums.ac.ir Published Online: January 06, 2016 These values were exported into a dbase IV table for further analysis. Three NDVI indicators were used for analysis, 1- Annual mean NDVI, 2- Seasonal NDVI differences, which is normalized differences of average NDVI for successive seasons in a 5km radius area around each stable and, 3- Local NDVI differences, which is nor- malized differences of average NDVI be- tween two areas: a 5 km radius area centered on each stable and the 5–10 km surrounding area (Fig. 1). Dispersal of adult mosquitoes is an im- portant factor in the spatial scale of trans- mission of vector-borne pathogens. The buff- ers of 5 and 5–10 km were chosen based on estimates of dispersal distances of Culex mosquitoes and their blood-sucking activi- ties on which birds (infected) at the first and then other hosts, such as equines. We as- sume that to bite a bird, the maximum esti- mated dispersal of a female mosquito is about 2.5km and similarly in another host (2.5×2=5km). The second buffer area was chosen by the same radius of 5km buffer. Moreover, NDVI contrast between two buff- ers could be an indicator of localized mos- quitoes as well as birds (Venkatesan and Rasgon 2010, Hamer et al. 2014). Distance to the nearest Wetland area Iranian wetlands possess a great diversity and provide a habitat for more than 140 mi- gratory and sedentary bird species. Among 42 kinds of wetlands in the world character- ized by the Ramsar Convention, all, except one, are found in Iran, indicating an im- mense variety of Iranian wetlands. Due to the extension of the country, there are a large number of wetlands with a broad range of size: area between a few to 500000 hectares (Bagherzadeh-Karimi and Rouhani-Rankouhi 2007). The data on wetlands was acquired from Iran Environment protection organization. For this study, we collected the data of 60 important wetlands area, of which 22 are Ramsar sites. The area of these wetlands extends around 2954970.20 hectares (Fig. 2). To define proximity to the wetlands, the shape files of wetlands and sampling sites of animals was imported into ArcMap version 9.3. Using the ‘Euclidean Distance’ func- tion, from the spatial analyst extension, near- est wetland and its distance to each sampling place was calculated. Statistical analysis Univariate analyses were performed in order to test the relationship between the outcome variable and each explanatory vari- able. Finally, logistic regression analyses were performed in order to identify which variables were involved in the stables seropositivity. Explanatory variables used in logistic re- gression modeling included annual average temperature, relative humidity and precipi- tation, distance to the nearest wetland, an- nual average NDVI, Normalized seasonal differences of NDVI, normalized local dif- ferences of NDVI and the type of stable (farm or club) (Table 1). Normalized sea- sonal differences of NDVI for two succes- sive seasons S1 and S2 and for stable S was calculated as: (NDVIS,S2–NDVIS,S1)–(1/N) Σ(NDVI*,S2–NDVI*,S1), where NDVIx,y is the average NDVI value observed in a 5km ra- dius area around stable x during season y, and N is the number of stables. Normalized local differences of NDVI were calculated between two areas; a 5 km radius area centered on the stable and the 5– 10 km surrounding area. For stable x and season y, value was calculated as: (NDVIx,y,0– 5km–NDVIx,y,5–10km)–(1/N) Σ(NDVIx,*,0–5km– NDVIx,*,5–10km), where NDVIx,y,z is the average NDVI value observed around stable x during season y in area z, and N is the number of stables. As the previous study showed a high effect J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 321 http://jad.tums.ac.ir Published Online: January 06, 2016 of the type of stables on the seropositivity (Ahmadnejad et al. 2011), we included this variable in the present analyses. Accuracy (percent of testing sites cor- rectly classified), sensitivity (percent of pos- itive testing sites correctly classified), and specificity (percent of negative testing sites correctly classified) were computed for the model. In addition, the area under the re- ceiver operating characteristic curve (AUC ROC) was calculated as indices of the fit of the model. All the statistical analyses were con- ducted using R 2.10.1, and maps were pro- duced with a Geographic Information Sys- tem (ESRI1, ArcGISTM 9.0). Results The stables with at least one positive animal considered as positive stable and 108 of 260 stables (41.5%) were positive for WNV infection. In the training set, on univariate analysis, all the variables were associated with the presence of WNV antibody (P< 0.05). The odds ratio and 95% confidence in- tervals were calculated for seropositivity and all risk factors. Multivariate logistic regres- sion analyses revealed a strong relation be- tween WNV antibody prevalence and annual mean temperature, distance to the nearest wetland, seasonal NDVI differences and Lo- cal NDVI differences variables (Table 1). The average minimum and maximum an- nual temperature were 5 °C and 27.7 °C, re- spectively. Annual mean temperature for infected and uninfected places was 19.97 °C (95% CI 18.97–20.97) and 15.58 (95% CI 15.04–16.11), respectively. The average minimum and maximum an- nual humidity were 25.6% and 83.8%, re- spectively. Annual mean humidity for in- fected places was 44.57% (95% CI 42.07– 47.06) and for uninfected places was 47.94% (95% CI 45.98–49.89). The average minimum and maximum an- nual precipitation were 50.8 mm and 1759.1 mm, respectively. Annual mean precipitation for infected places was 258.74 mm (95% CI 210.58–306.91) and for uninfected places was 299.18 mm (95% CI 266.83–331.53). The risk of being WNV seropositive in- creased for stables located in the area with higher temperature, as 5 °C increase in mean annual temperature was associated with a statistically significant 100% higher preva- lence of WNV infection. Fig. 3 Fig. 3 shows the aggregation of infected stables within the area with higher temperature. No significant association was observed between type of stable, annual mean humid- ity, precipitation, annual mean NDVI and seropositivity. It is notable that type of stable is different from their location, as the former more is indicating the management of the stable and the latter indicating geographic location. There was a significant negative correla- tion between distance to the nearest wetland and seropositivity of stables. The risk of be- ing WNV seropositive for stables increased (27%) with 10km decrease of distance from the nearest wetland area (Fig. 4). A logistic regression model showed an association between seasonal/ local NDVI and seropositivity. A higher significant dif- ference between NDVI of spring and sum- mer within a 5 km radius area around the stable increased the seropositivity four times higher. The difference of NDVI between two ar- eas, a 5 km radius area centered on the stable and the 5–10 km surrounding area, in spring for infected stables was four times higher that of uninfected stables. Fig. 5 shows annual mean NDVI in spring and distribution of the stables according to the normalized local NDVI differences. J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 322 http://jad.tums.ac.ir Published Online: January 06, 2016 The ROC analysis resulted in an AUC of 0.838, indicating the model had relatively good discrimination ability. The model computes probabilities for a cell to have a high WNV circulation status. Computing sensitivity and specificity of this prediction for various probability cut-offs clearly showed that at a decision threshold of 0.4. The model had a sensitivity of 75 % and a specificity of 83.6 %. Corresponding posi- tive and negative predictive values were 76.4 % and 82.5 %, respectively (Fig. 6). Fig. 1. The studied NDVI indicators for the sampling stables Table 1. Multivariate logistic model of the presence of anti-WNV seropositive animals in Iranian stables, 2008– 2009 Variable Odds-ratio (95% CI) p Type of stable (village of club) NS2 0.07 Annual mean temperature Δ=5°C 2.03 (1.10-3.92) 0.03 Annual mean %humidity1 <38% Reference 38–46% NS 0.12 46–52% NS 0.27 >52% NS 0.34 Annual mean precipitation1 <154 mm Reference 154–247 mm NS 0.84 247–327 mm NS 0.67 >327 mm NS 0.71 Distance to the nearest wetland area Δ=10 km 0.73 (0.54–0.99) 0.05 Annual mean NDVI NS 0.77 Seasonal NDVI differences Δspring-summer>0 2.37 (1.07–5.41) 0.04 Δsummer-autumn>0 NS 0.96 Δautumn-winter>0 NS 0.38 Δwinter-spring>0 NS 0.92 Local NDVI differences spring>0 4.26 (1.07–17.98) 0.04 summer>0 NS 0.11 autumn>0 NS 0.35 winter>0 NS 0.23 1Classes bounds based on distribution quantiles 2NS= Non- significant J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 323 http://jad.tums.ac.ir Published Online: January 06, 2016 Fig. 2. Important studied wetlands and sampling sites Fig. 3. Annual mean temperature and geographical distribution of WNV infected and uninfected stables in Iran Fig. 4. Distribution of positive stables (thick line) and of negative stables (thin line) according to the distance to the nearest wetland Fig. 5. Annual mean NDVI in spring in Iran and the geographic distribution of WNV infected and uninfected stables Fig. 6. Receiver operating curve of the multivariate logistic model of the presence of anti-WNV seropositive animals in Iranian stables, 2008–2009 Discussion The evidence presented here reveals the geographic factors associated with WNV circulation in Iran, which is reported for the first time in the country. Four studied fac- tors, temperature, distance to wetlands, local and regional differences NDVI were corre- lated to the WNV infection in equine. J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 324 http://jad.tums.ac.ir Published Online: January 06, 2016 Positive anomalies of the temperature in some southwestern province appear to have facilitated the mosquito abundance and, con- sequently, WNV infection in equine. We ob- serve a north-south gradient of seropositivity following the gradient of annual average tem- perature. Annual mean temperature above 22 °C was notable in which all, except three, of the stables by this temperature were sero- positive. These results are consistent with previous observations made in Egypt (Taylor et al. 1956). In southern provinces of the coun- try, where the climate is arid, the infection reaches a high level. In Northern provinces, the situation is quite different: warm, wet summers are common, but droughts occa- sionally occur in late summer and early fall; and outbreaks of WNV can occur at that time. Temperature in Iran gets warmer from west to east and from north to south. In- crease of temperature from west to east is due to concentration of the mountains at western part of the country while increase from north to south is because of approach- ing the equator and increasing of the solar radiation angle. Overall, northern and moun- tainous parts of the country have higher annual fluctuating and southern areas have relative stability (Alijani 1995). Temperature influences the mosquito life cycle, its reproduction rates, development of eggs within the mosquito (gonotrophic cy- cle), the vector–host contact rate and conse- quently, virus transmission. Although, Patz and et al. (2000) suggested that higher tem- peratures may increase or reduce survival rate, depending on the vector, its behavior, ecology, and many other factors. Thus, the probability of transmission may or may not be increased by higher temperatures (Taylor et al. 1956, Patz et al. 2000), but it has been shown that even a small increase in temper- atures can have a significant impact on trans- mission of WNV by mosquitoes (Kilpatrick et al. 2008). Effect of high temperature on increasing of mosquito abundance and vec- tor competency, has also been experimental- ly proven (Dohm et al. 2002). WNV is endemic in people living in the semi-arid regions, where the weather most of the year is decidedly hot and dry. Cases has been regularly reported in sub-desertic areas around the Mediterranean basin: (i) along the Nile Valley (Egypt, 1950, 1994) (ii) along the Syrian-African Rift (Israel 1951, 1998– 2000, Jordan 2003), (iii) in the Timimoun oasis in the Central Sahara (Algeria 1994) (Taylor et al. 1956, Murgue et al. 2001). The warmer condition of these area, along with directly effects on vector and host, cause the standing water with high rich organic ma- terials attracting the mosquito and birds, con- sequently increase the interaction and circu- lation of WNV between them (Epstein and Defiippo 2001). Beside the vectors, temperature often plays a significant role in bird life cycles, avail- ability of food and habitat and also in the distribution of migratory birds (Greenberg and Marra 2005). Precipitation and humidity could have a major influence on the distribution of mosquito species. Humidity increases vector flight and host-seeking behaviors and precipitation is necessary for the formation of mosquito breed- ing habitats (Shaman et al. 2002). However, we could not find a significant relation be- tween WNV infection and annual mean hu- midity/rainfall. Evidences on relationship be- tween WNV/mosquito abundance and precip- itation and humidity are variable, some indicate a weak, some strong correlations (Soverow et al. 2009, Wan Norafikah et al. 2009). Seropositivity was negatively correlated to the distance to wetlands, so, near to wet- land can increase infection of the animals to the virus. This is consistent with the previ- ous studies, which showed WNV infection is more strictly linked to the wetlands with abundant bird populations, especially migra- tory birds. Wetlands provide suitable breed- ing habitat for mosquitoes and birds and J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 325 http://jad.tums.ac.ir Published Online: January 06, 2016 consider as a natural foci of WNV infections in Palearctic. Outbreaks in temperate area, especially in Europe, often have occurred in or near wetlands (Hubalek and Halouzka 1999, Jourdain et al. 2007). In 1998, WNV was isolated from horses suffering from neurologic disease and residing in a large wetland area in Italy. In Camargue (France), wetland variations identified as important risk factors in WNV spillover in horses. (Autorino et al. 2002, Pradier et al. 2014). Moreover, proximity to wetland increase mosquito abundance, as the first three stages of their life is aquatic dependent. Iranian wetlands serve as wintering or staging area for migratory birds coming from WNV en- demic area and using the West Siberian- Caspian-East African and Central Siberian- Indus-South Asian flyways (UNDP 2004). The strongest association was detected between seropositivity and NDVI, at seasonal and, especially, local differences. While this relation could not be found with annual mean value of NDVI. The local NDVI’s contrast during spring is a major risk factor of the transmission of West-Nile virus in Iran. During the early drought season, the stable seropositivity is increasing with the contrast between the local and regional NDVI environment of the stable. This so- called oasis effect is associated to the at- tractiveness of the local NDVI environment for WNV’s vectors (mosquitoes and birds). This oasis effect disappeared in autumn, and the local NDVI’s contrast becomes protector (negative impact). This oasis effect is consistent with the seasonal production of vegetation in spring, which is more important compared to the national average (Contrast spring- sum- mer). The higher level of seropositivity of the stables can be related to the attraction of many birds to the resources around these stables (Bock et al. 2008). The extreme cli- matic situation of the southwestern prov- inces of Iran encourages birds to congregate around shrinking water sites, encouraging viral circulation among birds and mosquitoes, while heat accelerates viral maturation. Droughts and heat wave were found to induce WNV amplification by bringing the hosts and the vectors together. In opposition with the dog- ma that increasing precipitation predicts mosquito abundance, some authors hypoth- esize that wild populations should generate outbreaks under drought conditions. Stand- ing water pools become richer in the organic material that Culex vectors needs to thrive and the mosquito predators, such as amphib- ians and dragonflies, are fewer in number (Epstein 2004). These findings are consistent with the hypothesis that spring= nesting (in an "oasis", the birds will nest closer to the horses) and autumn = dispersal of juveniles (juveniles are going further if they leave from an "oasis"). These findings are also consistent with the fact that viral circulation is related to the presence of chicks, seden- tary and still bare skin in altricial species (Marra et al. 2004). High values of NDVI are related to higher photosynthetic activity and improved ecological condition. This suitable condition increases the main food supply, such as ar- thropods, for most migratory birds. Several studies have used NDVI for assessments of WNV infection foci and its vector habitat. Seasonal difference in NDVI was one the best predictive factors for determining mos- quito distribution abundance in the model studied by Jacob et al (Jacob et al. 2009). Ward et al. investigated association between NDVI and cases of equine WNV encepha- lomyelitis. They found that the mean NDVI in biweekly periods with reported cases was significantly higher than the mean NDVI in periods without cases (Ward 2009). Climatic condition of Iran varies through- out the country, resulting in a complex hy- drological and vegetative landscape, with arid and semiarid areas in the center, south- east and southwest and temperate rain forest in the north. Although areas with high WNV J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 326 http://jad.tums.ac.ir Published Online: January 06, 2016 infection rates, such as Khuzestan, have a dry climate, but their water sources are mainly from rivers originated from mountains and precipitations, which play an important role in producing oasis areas in these areas. Although these results provide a new un- derstanding of some ecological parameters effects on WNV circulation in Iran, addi- tional work is needed. Arboviruses circula- tion is multifactorial and studies are needed to evaluate various mosquito species for their potential to transmit WNV in Iran. There are many studies on mosquito fauna in Iran, however there is no specific infor- mation about WNV vectors. WNV has been identified in numerous mosquito species, including members of the genera Culex, Ae- des and Ochlerotatus, but mosquitoes be- longing to the Culex species are the main vectors for WNV. Culex pipiens, the main vector of WNV in US, has spread in differ- ent ecological zones of Iran and its distribu- tion is very similar to its climatic distribution in North, South America and Africa. The virus has been isolated from Cx. vishnui complex in Pakistan. There are some doubt- ful records of this mosquito in Iran; how- ever, based on the record of this species in Pakistan, it seems that Cx. vishnui may occur in southeastern Iran (Zaim et al. 1985, Hubalek and Halouzka 1999, Vatandoost et al. 2004, Dehghan et al. 2010). If we consider Iranian WNV vector within Culex spp, it is speculated that other than several species may be involved in WNV circulation in south and southwestern parts of Iran, where there is a high circulation of the virus. However, vector competency studies should be conducted on various species of mosquitoes in order to determine the ability of the mosquito species to transmit WNV and detection of the virus from mosquitoes. Another study was conducted to deter- mine WNV antibody and genome in differ- ent species of wild water birds captured from wetlands in Iran. Samples were collected from 26 different species, 15% of which were serologically positive, while no WNV viral RNA-positive samples were found in this study. The majority of positive birds were common coot (Fulica atra) (Fereidouni et al. 2011). They did not find any positive samples from Khuzestan province where the most important region for the virus circula- tion is. Moreover, they had a few samples (n=4) from this province, all of which were birds other than common coot, maybe due to lower abundance of common coot (Fulica atra) in this area. The area identified with high risk for WNV circulation may be used for entomological or epidemiological surveillance. Conclusion The model indicated that local NDVI’s contrast during spring is a major risk factor of the transmission of West-Nile virus in Iran. This so-called oasis effect consistent with the seasonal production of vegetation in spring, and is associated to the attractiveness of the local NDVI environment for WNV vectors and hosts. Further studies are needed to understand better the ecology of WNV in Iran. The reservoir hosts and vectors are the keys of circulation of the virus and will be the focus of the future work. Acknowledgements We thank all the personnel of IR of Iran Meteorological Organization (IRIMO), espe- cially Dr Ghaemi. We also thank the person- nel of Iran Environment Protection Organi- zation, especially Ms Elahi-Rad and Mr Amini. The authors declare that there is no conflict of interests. References Ahmadnejad F, Otarod V, Fallah MH, Lowenski S, Sedighi-Moghaddam R, J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 327 http://jad.tums.ac.ir Published Online: January 06, 2016 Zavareh A, Durand B, Lecollinet S, Sabatier P (2011) Spread of West Nile virus in Iran: a cross-sectional serosurvey in equines, 2008–2009. Epidemiol In- fect. 139(10): 1587–1593. Alijani B (1995) Climate of Iran. Piame Noor University, Tehran. Autorino GL, Battisti A, Deubel V, Ferrari G, Forletta R, Giovannini A, Lelli R, Murri S, Scicluna MT (2002) West Nile virus epidemic in horses, Tuscany region, Italy. Emerg Infect Dis. 8(12): 1372–1378. Azari-Hamidian S (2007) Checklist of Iranian mosquitoes (Diptera: Culicidae). J Vec- tor Ecol. 32(2): 235–242. Bagherzadeh-Karimi M, Rouhani-Rankouhi M (2007) A guide to the Iranian wetlands registered in Ramsar Convention. Ruze No, Tehran. Bian L, Li L, Yan G (2006) Combining global and local estimates for spatial distribution of mosquito larval habitats. GIsci Remote Sens. 43(2): 128–141. Bock CE, Jones ZF, Bock JH (2008) The oasis effect: response of birds to ex- urban development in a southwestern savanna. Ecol Appl. 18(5): 1093–1106. Bolling BG, Kennedy JH, Zimmerman EG (2005) Seasonal dynamics of four potential West Nile vector species in north-central Texas. J Vector Ecol. 30 (2): 186–194. Brown H, Duik-Wasser M, Andreadis T, Fish D (2008) Remotely-sensed vege- tation indices identify mosquito clus- ters of West Nile virus vectors in an urban landscape in the northeastern United States. Vector Borne Zoonotic Dis. 8(2): 197–206. Chen CC, Jenkins E, Epp T, Waldner C, Curry PS, Soos C (2013) Climate change and West Nile virus in a highly endemic region of North America. Int J Environ Res Public Health. 10(7): 3052–3071. Dehghan H, Sadraei J, Moosa-Kazemi S (2010) The morphological variations of Culex pipiens larvae (Diptera: Culicidae) in Yazd Province, central Iran. Iran J Arthropod Borne Dis. 4(2): 42–49. Dohm DJ, O'Guinn ML, Turell MJ (2002) Effect of environmental temperature on the ability of Culex pipiens (Diptera: Culicidae) to transmit West Nile virus. J Med Entomol. 39(1): 221–225. Epstein PR (2004) Climate change and pub- lic health: emerging infectious diseases. Encyclopedia of Energy. 1: 381–392. Epstein PR, Defiippo C (2001) West Nile virus and drought. Global Change and Human health. 2(2): 105–107. Fereidouni SR, Ziegler U, Linke S, Niedrig M, Modirrousta H, Hoffmann B, Groschup MH (2011) West Nile virus monitoring in migrating and resident water birds in Iran: are common coots the main reservoirs of the virus in wetlands? Vector Borne Zoonotic Dis. 11(10): 1377–1381. Fortes PS, Platonov AE, Pereira LS (2005) GISAREG- A GIS based irrigation scheduling simulation model to suport improved water use. Agr Water Ma- nage. 77(1): 159–179. Gibbs SE, Wimberly MC, Madden M, Masour J, Yabsley MJ, Stallknecht DE (2006) Factors affecting the geog- raphic distribution of West Nile virus in Georgia, USA: 2002–2004. Vector Borne Zoonotic Dis. 6(1): 73–82. Greenberg R, Marra PP (2005) Birds of two worlds: the ecology and evolution of migration. JHU Press, Baltimore. Gubler DJ, Reiter P, Ebi KL, Yap W, Nasci R, Patz JA (2001) Climate variability and change in the United States: po- tential impacts on vector- and rodent- borne diseases. Environ Health Pers- pect. 109 (Suppl 2): 223–233. Hamer GL, Anderson TK, Donovan DJ, Brawn J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 328 http://jad.tums.ac.ir Published Online: January 06, 2016 JD, Krebs BL, Gardner AM, Ruiz MO, Brown WM, Kitron UD, Newman CM, Goldberg TL, Walker ED (2014) Dispersal of adult Culex mosquitoes in an urban west nile virus hotspot: a mark-capture study incor-porating stable isotope enrichment of natural larval habitats. PLoS Negl Trop Dis. 8(3): e2768. Heshmati GA (2007) Vegetation characteristics of four ecological zones of Iran. Int J Plant Prod. 2: 215–224. Hubalek Z, Halouzka J (1999) West Nile fever- a reemerging mosquito-borne viral disease in Europe. Emerg Infect Dis. 5(5): 643–650. Jacob BJ, Gu W, Caamano EX, Novak RJ (2009) Developing operational algo- rithms using linear and non-linear squares estimation in Python for the identification of Culex pipiens and Culex restuans in a mosquito abatement district (Cook County, Illinois, USA). Geospat Health. 3(2): 157–176. Jourdain E, Gauthier-Clerc M, Bicout DJ, Sabatier P (2007) Bird migration routes and risk for pathogen dispersion into western Mediterranean wetlands. Emerg Infect Dis. 13(3): 365–372. Kilpatrick AM, Meola MA, Moudy RM, Kramer LD (2008) Temperature, viral genetics, and the transmission of West Nile virus by Culex pipiens mos- quitoes. PLoS Pathog. 4(6): e1000092. Kinney RM, Huang CY, Whiteman MC, Bowen RA, Langevin SA, Miller BR, Brault AC (2006) Avian virulence and thermostable replication of the North American strain of West Nile virus. J Gen Virol. 87(Pt 12): 3611–3622. Marra PP, Griffing S, Caffrey C, Kilpatrick AM, McLean R, Brand C, Saito E, Dupuis AP, Kramer L, Novak R (2004) West Nile virus and wildlife. Bio Science. 54(5): 393–402. Murgue B, Murri S, Triki H, Deubel V, Zeller HG (2001) West Nile in the Medite- rranean basin: 1950–2000. Ann N Y Acad Sci. 951: 117–126. Patz JA, McGeehin MA, Bernard SM, Ebi KL, Epstein PR, Grambsch A, Gubler DJ, Reither P, Romieu I, Rose JB, Samet JM, Trtanj J (2000) The poten- tial health impacts of climate varia- bility and change for the United States: executive summary of the report of the health sector of the U.S. National As- sessment. Environ Health Perspect. 108 (4): 367–376. Pradier S, Sandoz A, Paul MC, Lefebvre G, Tran A, Maingault J, Lecollinet S, Leblond A (2014) Importance of wetlands management for West Nile Virus circulation risk, Camargue, Southern France. Int J Environ Res Public Health. 11(8): 7740–7754. Reiter P (2001) Climate change and mosquito- borne disease. Environ Health Perspect. 109(1): 141–161. Ruiz MO, Chaves LF, Hamer GL, Sun T, Brown WM, Walker ED, Haramis L, Goldberg TL, Kitron UD (2010) Local impact of temperature and precipi- tation on West Nile virus infection in Culex species mosquitoes in northeast Illinois, USA. Parasit Vectors. 3(1): 19. Ruiz MO, Tedesco C, McTighe TJ, Austin C, Kitron U (2004) Environmental and social determinants of human risk during a West Nile virus outbreak in the greater Chicago area, 2002. Int J Health Geogr. 3(1): 8. Saidi S. TR, Javadian E, Nadim A (1976) The prevalence of human infection with West Nile virus in Iran. Iran J Publ Hlth. 5(1): 8–13. Shaman J, Day JF, Stieglitz M (2002) Drought-induced amplification of Saint Louis encephalitis virus, Florida. Emerg Infect Dis. 8(6): 575–580. Soverow JE, Wellenius GA, Fisman DN, Mittleman MA (2009) Infectious disease J Arthropod-Borne Dis, September 2016, 10(3): 317–329 F Ahmadnejad et al.: Impact of Climate … 329 http://jad.tums.ac.ir Published Online: January 06, 2016 in a warming world: how weather in- fluenced West Nile virus in the United States (2001–2005). Environ Health Perspect. 117(7): 1049–1052. Stenseth NC, Mysterud A, Ottersen G, Hurrell JW, Chan KS, Lima M (2002) Ecological effects of climate fluctuations. Science. 297(5585): 1292–1296. Stow DA, Hope A, McGuire D, Verbyla D, Gamon J, Huemmrich F, Houston S, Racine C, al. E (2004) Remote sensing of vegetation and land-cover change in Arctic Tundra Ecosystems Remote Sens Environ. 89(3): 281–308. Taylor RM, Work TH, Hurlbut HS, Rizk F (1956) A study of the ecology of West Nile virus in Egypt. Am J Trop Med Hyg. 5(4): 579–620. United Nations Development Programme (UNDP) (2004) Conservation of Iranian Wetlands. Available at: http://www. thegef.org/gef/sites/thegef.org/files/rep ository/Iran_Conservation_of_Iranian_ Wetlands.pdf. Vatandoost H, Ezeddinloo L, Mahvi A, Abai M, Kia E, Mobedi I (2004) Enhanced tolerance of house mosquito to dif- ferent insecticides due to agricultural and household pesticides in sewage system of Tehran, Iran. Iran J Environ Health Sci Eng. 1(1): 46–50. Venkatesan M, Rasgon JL (2010) Population genetic data suggest a role for mos- quito-mediated dispersal of West Nile virus across the western United States. Mol Ecol. 19(8): 1573–1584. Wan Norafikah O, Chen CD, Soh HN, Lee HL, Nazni WA, Sofian-Azirun M (2009) Surveillance of Aedes mosquitoes in a university campus in Kuala Lumpur, Malaysia. Trop Biomed. 26(2): 206–215. Ward MP (2009) Equine West Nile virus disease occurrence and the Normalized Difference Vegetation Index. Prev Vet Med. 88(3): 205–212. Zaim M, Manouchehri A, Yaghoobi-Ershadi M (1985) Mosquito fauna of Iran 2- Culex. Iran J Publ Health. 14: 1–12.