Skip to main content
Advertisement
  • Loading metrics

Spatial patterns of West Nile virus distribution in the Volgograd region of Russia, a territory with long-existing foci

Abstract

Southern Russia remains affected by West Nile virus (WNV). In the current study, we identified the spatial determinants of WNV distribution in an area with endemic virus transmission, with special reference to the urban settings, by mapping probable points of human infection acquisition and points of virus detection in mosquitoes, ticks, birds, and mammals during 1999–2016. The suitability of thermal conditions for extrinsic virus replication was assessed based on the approach of degree-day summation and their changes were estimated by linear trend analysis. A generalized linear model was used to analyze the year-to-year variation of human cases versus thermal conditions. Environmental suitability was determined by ecological niche modelling using MaxEnt software. Human population density was used as an offset to correct for possible bias. Spatial analysis of virus detection in the environment showed significant contributions from surface temperature, altitude, and distance from water bodies. When indicators of location and mobility of the human population were included, the relative impact of factors changed, with roads becoming most important. When the points of probable human case infection were added, the percentage of leading factors changed only slightly. The urban environment significantly increased the epidemic potential of the territory and created quite favorable conditions for virus circulation. The private building sector with low-storey houses and garden plots located in the suburbs provided a connection between urban and rural transmission cycles.

Author summary

West Nile fever is a zoonotic disease caused by flavivirus infection that dispersed circumglobally by the end of the 20th century. There is the long-existing endemic area in southern Russia (Volgograd region) that remains affected by West Nile virus (WNV) with the maximum number of reported cases in the country. The natural features of the region formed favorable conditions for long-term WNV circulation. At the same time, the urban environment could pose greater hazards for WNV infection in humans. Road density, distance from water bodies, land surface temperature and altitude were identified as the most significant spatial determinants of WNV distribution. The private building sector with low-storey houses and garden plots located in the city suburban area provided connections between urban and rural transmission cycles. Climate change also may favor WNV transmission as increased temperatures affect the sum of effective temperatures necessary for virus extrinsic replication. A statistically significant correlation between WNV human cases and the sum of effective temperatures was found during the last 20 years.

Introduction

West Nile fever (WNF) is a zoonotic disease caused by infection with a flavivirus (Flaviviridae family) that is endemic in Africa, the Middle East, and South Asia, and then dispersed circumglobally by the end of the 20th century [1,2]. Viral transmission originally occurred between mosquitoes and birds, but currently, circulation has involved more hosts and vectors outside the original ecological niche [3]. Transmission via mosquitoes may cause severe infection in humans and horses that are dead-end hosts for the virus.

Since the middle of the 20th century, there has been evidence of a changing geographical distribution of West Nile virus (WNV), as well as a shift in epidemiology and activity of the virus. The emergence of WNF in new regions has been associated with climate change, however transmission, distribution, and incidence of West Nile virus disease are influenced by multiple factors [4]. More information about WNV spatial distribution can be found in worldwide reviews [59] and regional studies in Europe [1013], Africa and Eastern Mediterranean region [1415] and North America [1618].

West Nile virus spatio-temporal distribution in Russia

In Russia, WNV was first isolated from larvae and nymphs of Hyalomma marginatum (syn. H. plumbeum) ticks in the Astrakhan region (hereinafter, the ‘region’ refers to an official territorial administrative unit) in the Volga River Delta in 1963 and was associated with Lineage 1 [19,20]. The first laboratory confirmed WNF cases were detected in the Astrakhan region in 1967 [21]. However, the first official WNF case was not recorded from the same region until 1997. The first documented outbreak occurred in 1999 in the Astrakhan and neighbouring Volgograd regions (475 cases) and again WNV Lineage 1 was detected [20]. The most extensive outbreak occurred in 2010, with 521 cases (Fig 1) [22]. The virus spread to seven regions, with the Volgograd region most affected (413 cases). Additional human cases later were reported from other southern regions of Russia, in addition to the long-existing foci in the Astrakhan and Volgograd regions [23] (Fig 2). However, the Volgograd region has remained the area with the maximum number of reported cases. The second major outbreak in 2010 was associated with Lineage 2 [24]. Sequencing of virus strains recently isolated from European Russia showed they were different from the strains isolated 20–30 years ago [25]. The present strains are similar to the strains from Europe, identificated as Lineage 2 [24].

thumbnail
Fig 1. WNV laboratory confirmed cases among patients with fever or neuroinvasive disease in Russia and the Volgograd region, 1997–2019.

Source: official records of Federal Service for Surveillance on Consumer Rights Protection and Human Wellbeing (Rospotrebnadzor) upon request.

https://doi.org/10.1371/journal.pntd.0010145.g001

thumbnail
Fig 2. Spatial distribution of WNV laboratory confirmed cases among patients with fever or neuroinvasive disease per administrative units in Russia, 1997–2019.

Source: official records of Federal Service for Surveillance on Consumer Rights Protection and Human Wellbeing (Rospotrebnadzor) upon request. Contains information from OpenStreetMap and OpenStreetMap Foundation, which is made available under the Open Database License. The numbers indicate: 1 Adygea; 2 Astrakhan region; 3 Belgorod region; 4 Volgograd region; 5 Voronezh region; 6 Kaluga region; 7 Krasnodar krai; 8 Kursk region; 9 Lipetsk region; 10 Novosibirsk region; 11 Omsk region; 12 Kalmykia; 13 Rostov region; 14 Samara region; 15 Saratov region; 16 Stavropol krai; 17 Tatarstan; 18 Tula region; 19 Ulyanovsk region; 20 Chelyabinsk region.

https://doi.org/10.1371/journal.pntd.0010145.g002

Hosts and vectors

In the natural environment, the transmission cycle involves wild birds as principal hosts and mosquitoes that feed on birds as vectors [26]. Fifty-six bird species have been involved in virus transmission in wetlands and surroundings of the Volga delta [27]. Storks and other Ciconiiformes were most frequently identified as a virus host among wild birds as well as the great cormorant (Phalacrocorax carbo), Eurasian coot (Fulica atra), common moorhen (Gallinula chloropus), great crested grebe (Podiceps cristatus), and seagull and sterna (family Laridae) (based on virus detection among wild birds) [27,28]. Synanthropic birds, primarily rooks, crows, pigeons, and sparrows, have been involved in the transmission cycle in urban areas [29,30].

There are fifteen mosquito species, mainly of the genus Culex and, less commonly, of Aedes and Anopheles genera, that have been considered potential vectors of WNV in Russia [31]. Culex spp. mosquitoes (Culex pipiens L. and Culex modestus Fic.) are most epidemiologically important [3234] and have been confirmed to be competent to transmit WNV [35]. Culex pipiens is the main vector within bird populations, whereas Cx. modestus dominates in reeds near their larval habitat as well as in urban habitats and can be a bridge vector from birds to humans. The autogenous anthropophilic biotype Cx. pipiens molestus has been found in flooded basements and may be involved in urban settings [36].

Ticks from Ixodidae (mainly Hyalomma marginatum) and Argasidae (mainly Ornithodoros coniceps) families have been implicated as having a minor role in transmission [2], with their important role in preserving the virus during winter. The infection rate of all tick metamorphosis phases is much higher compared to mosquitoes [37,38].

Environmental drivers

The transmission of WNV depends on several environmental factors. Ambient temperature is a direct factor that influences viral replication and growth rates of mosquitoes [26]. Virus transmission by vectors is possible provided that a certain amount of heat (effective temperature or degree days) is accumulated. Degree days are calculated by summing the daily excess of the average daily temperature above a certain threshold value [39]. A total of 109-degree days above the threshold of 14.3°C seems necessary for WNV replication [40]. An increase in temperature reduces the length of the mosquito gonotrophic cycle and shortens the extrinsic incubation period of the virus in the vector [4042]. However, some studies showed that increased temperature can reduce adult mosquito longevity [4345]. According to experiments [43], Cx. pipiens f. pipiens longevity decreased with increasing temperature and did not exceed 12 days at 30°C. The maximum mosquito survival was about 130 days at 15°C and the optimal temperatures for the co-occurrence of virus replication and Cx. pipiens can range from 14.3 to 30°C [44,45]. However, shortening of the extrinsic incubation period due to warm temperatures can compensate for the reduced survival of mosquitoes, making transmission more efficient. For example, transmission rates of northern European Cx. pipiens increased from 0% to 33% in the temperature range from 18°C to 28°C [35].

Precipitation can affect mosquitoes and WNV transmission in different ways. On the one hand, precipitation may increase the WNF incidence due to a higher abundance of mosquitoes [26,46] and their ability to survive at a higher humidity. Conversely heavy rainfall can ‘flush out’ larvae from their habitat [47]. Weak correlations and inconsistent patterns between precipitation and Cx. pipiens dynamics were found across Europe [48]. Such a weak correlation may be due to the inconsistent nature of precipitation and also due to the fact that the natural and urban environments react differently to rainfall. Although for wetlands, it is important to consider the accumulation of precipitation in periods preceding adult mosquito sampling, for urban environments this factor does not have much effect as most precipitation is lost from impervious surfaces into drainage systems. In urban environments it may matter in case of light rainfall or landscape irrigation that can improve mosquito production.

The distribution of mosquito habitat is determined by vegetation cover and the presence of water bodies in the natural environment [49,50]. Wetlands are the most important ecosystems for enzootic transmission due to the presence of resident and migratory bird and mosquito populations, particularly Cx. pipiens [51]. Urban areas with a simple avian fauna consisting of high density of competent hosts can provide even more efficient transmission. Indices such as the normalized difference vegetation index (NDVI) and normalized difference water index (NDWI) have been identified as risk indicators for WNV outbreaks in previous studies in Europe [52], North America [53], and Russia [54]. In cities, conditions often favor Cx. pipiens production as they utilize stagnant and polluted water bodies, such as drainage ditches, ponds, abandoned barrels, and other man-made reservoirs, as well as flooded basements of buildings [55,56].

The role of elevation in WNV transmission depends on the ecology of the mosquito vectors involved. Usually, the abundance and mosquito species diversity are the greatest in lowlands [26,57]. Often, the risk of infection decreases with increasing elevation [5760].

WNV transmission can occur in both rural and urban environments. The urban environment can contribute to WNV transmission by involving synanthropic birds and creating additional man-made mosquito breeding sites [31,6163]. Johnson et al. [64] showed that residential areas contained significantly higher proportions of WNV-competent mosquito species and avian host species when compared to adjacent urban wetlands. Urban environment heterogeneity may cause decelerating speed of pathogen spread waves in incipient epidemics [65]. Overall, the urban environment influences the distribution and mobility of the human population and other dimensions of social life.

In this paper, we aimed to identify the spatial determinants of WNV distribution in an endemic area in southern Russia focusing on the urban environment and the assessment of thermal suitability for WNV transmission. We analyzed the contribution of various environmental and urban variables as determinants of WNV distribution according to different scenarios of spatial modelling. Then, we calculated the effective temperatures and length of periods during the years when replication and infection were possible. Finally, we explored trends in weather variables during recent years in association with climate change.

Methods

Study area

The study area included Volgograd city and its suburbs in southern Russia between 48.38° N and 49.07° N and 43.97° E and 45.11° E (Fig 2). This territory consisted of typical multi-storey urban areas as well as traditional villages and groups of summer cottages that belong to city dwellers and are used as summer vacation spots and for small-scale production of fruit and vegetables for personal consumption.

This area has a hot summer continental climate (Dfa) according to the Köppen climate classification [66]. The average temperature for the warmest month (July) is 23.6°C, and the average temperature for the coolest month (February) is -6.6°C (S1 Fig). The average annual precipitation is 348 mm, with the maximum in May and the minimum in February, April, and September [67].

Volgograd is a large industrial center located on the Volga River, with nearly one million inhabitants (Fig 3). The city stretches along the river for approximately 60 km. The immediate vicinity of the city on the right bank of Volga (to the West) is elevated arid steppes. In contrast, the immediate vicinity of the city on the left bank (to the East) is the Volga-Akhtuba floodplain.

thumbnail
Fig 3. Study area–Volgograd city and its suburbs in southern Russia.

Contains information from OpenStreetMap and OpenStreetMap Foundation, which is made available under the Open Database License.

https://doi.org/10.1371/journal.pntd.0010145.g003

Within the city, along with industrial zones and areas of urban development, there are lowlands with wetlands, flooded Volga islands with various rural settlements, and recreational zones. These islands provide a natural buffer between the large industrial city and floodplain ecosystems of the Volga-Akhtuba floodplain [68]. In the city territory, there are several small rivers with hydrophilous vegetation along the banks, many ponds, and a ravine and gully network as well as artificial water bodies of various size. The wetlands in the vicinity of Volgograd are crossed by many ancient river channels, lakes, and bogs and are part of the migratory routes of birds. Thus, in Volgograd, both rural and urban circulation of the virus can be maintained largely due to the diversity of mosquito habitats.

Virus detection data

Data on virus detection in the environment and points of probable human infection were used to determine WNV distribution. The source of the data was the official record of the Federal Service for Surveillance on Consumer Rights Protection and Human Wellbeing (Rospotrebnadzor). Virus detection in the environment was carried out by epidemiologists of Rospotrebnadzor according to the methodological guidelines “MU 3.1.3.2600–10 Activities against West Nile fever in the Russian Federation” and included testing for virus by RT-PCR. Samples for testing included arthropods (ticks and mosquitoes), animal carcasses, and bird feces. Samples were collected from different places in natural and urban environments during 1999–2016. All locations positive for WNV RNA were recorded as virus detection sites.

A total of 29 places with a positive RNA test were found. Additional information about possible points of laboratory confirmed human infection during 2011 was added to strengthen the spatial analysis. Each human disease case was investigated by epidemiologists of Rospotrebnadzor to determine the probable point of infection. 2011 was considered a “typical” year, with an average number of officially reported WNF cases in Volgograd city and its neighbouring areas (Table 1). We assumed that these 55 points of infection corresponded to the most favorable places for WNV presence in the environment.

thumbnail
Table 1. The statistical analysis of human WNV cases in the Volgograd region, 1997–2019.

https://doi.org/10.1371/journal.pntd.0010145.t001

Thus, we used a general dataset on WNV distribution in the environment, whether it was actual detection or potential presence. All 89 sites were geocoded and mapped using ArcGIS software (Fig 4).

thumbnail
Fig 4. WNV detection sites (1996–2016) and possible places of human infection (2011) in the environment in Volgograd and neighbouring areas.

Contains information from OpenStreetMap and OpenStreetMap Foundation, which is made available under the Open Database License.

https://doi.org/10.1371/journal.pntd.0010145.g004

Environmental data

The daily data from the Volgograd weather station were used to analyze the effects of temporal change in temperature, precipitation, and relative humidity on the occurrence of human WNV cases during 1997 to 2019. This station is operated under the standards of World Meteorological Organization (WMO, ID-34561). Raw data were downloaded from the Russian Institute for Hydrometeorological Information—World Data Center (RIHMI-WDC, source: www.meteo.ru).

To focus on the spatial heterogeneity of environmental conditions, we examined the following environmental variables as natural environment predictors of WNV distribution: land surface temperature (LST), normalized difference vegetation index (NDVI), normalized difference water index by McFeeters (NDWI), distance from water bodies, and elevation. LST was calculated through Landsat thermal infrared imagery converting Band 6 spectral radiance to the brightness temperature according to the “Landsat 7 (L7) Data Users Handbook—Version 2.0” (https://www.usgs.gov/media/files/landsat-7-data-users-handbook) and then upgraded by making atmospheric correction using MODTRAN model [69] considering land surface emissivity based on previous methods [6971]. Cloud screening and handling the resulting missing data was carried out in a semi-automatic mode: for each image, a threshold value for the brightness of pixels that correspond to clouds was determined and all pixels with values less than the threshold were excluded from further analysis. LST was used in spatial modelling as a proxy for air temperature, taking into account the difference between air temperature and LST in areas with low vegetation [7274]. NDVI [75] was calculated using red and near-infrared wavelengths. NDWI [76] was determined using green and near-infrared wavelengths. For each of the three variables, the average values for all images were calculated using a raster calculator of ArcGis software.

The urban environment was represented by three variables: building density, road density, and railway density. The inclusion of the railway density was included to take into account frequent summer trips of city dwellers out of the city to personal household plots (gardens) for summer vacations or weekends. Such seasonal migrations are very typical for the urban population of Russia and railways are the main transport mode for these purposes. In this case, ‘city dwellers’ moving to suburban potential places of infection is possible. Variables in the urban environment were based on vector data from OpenStreetMap and were calculated in ArcGis software. The density of linear objects was calculated using the Kernel algorithm, the building density was calculated as the ratio of the area of buildings to the total area on a regular grid cell.

Population density was used as a bias factor to account for the assumption of greater likelihood of WNV detection in places with a higher population density (Table 2). A population density raster (S2 Fig) was created within the study area and incorporated as a bias file into the model. Maxent software factors out the bias by assigning weights to the random background data during modeling [77].

thumbnail
Table 2. Response, explanatory, and offset variables for spatial modelling of heterogeneity of environmental conditions for WNV distribution.

https://doi.org/10.1371/journal.pntd.0010145.t002

Spatial modelling of environmental suitability for WNV transmission

The environmental suitability for WNV was modelled using an ecological niche modelling (ENM) approach using MaxEnt software. Maxent is a machine learning algorithm based on the maximum entropy approach [81]. This method is widely used to determine the habitat suitability for a particular species according to the environmental conditions [82] or to model the suitability of an area for disease occurrence considering limiting and risk factors [83]. This maximum entropy technique is based on the principle of geospatial regression. It describes the known distribution of the study phenomenon using presence-only data [81,84]. The program randomly generates absence points of the phenomenon and takes a set of environmental predictors across a user-defined landscape that is divided into grid cells [85]. As a result, this method produces a map, where each grid cell expresses the suitability as a function of the environmental variables, based on the input presence data.

In the present study 89 sites of WNV detection (29 sites of virus detection and 55 probable human infection sites) in the environment were used as presence data, while environmental data were presented by natural and/or urban environmental variables (Table 2). According to recommendations found elsewhere [86,87], the number of input presence locations was sufficient for accurate model development. According to these studies the minimum limit sample size varied around 15–25 presence observation with the optimum size no less than 30 [86,87]. Despite the presence of unstable results when applying models over a large area, the use of the low sample size in smaller radius around the occurrence points could provide good working models [88]. In contrary to widespread non-parametric tests, such as generalized linear mixed models or generalized additive mixed models, that usually require numerically powerful sampling, MaxEnt can be used specifically for small sample sizes [8689].

Prior to the MaxEnt analysis, the environmental variables were tested for multicollinearity using the ’usdm’ package in the R 4.0.3 software environment [90]. It was revealed that the only pair of correlated variables was the NDVI and NDWI (Pearson correlation coefficient above 0.85), of which the former had a larger variance inflation factor (VIF). According to VIF values NDVI was excluded from the further analysis. However, we checked separately the contribution of NDVI and NDWI in models and results were similar (S1 Table). Three simulation scenarios were developed (Table 3).

thumbnail
Table 3. The set of environmental predictors in different combinations for modelling.

https://doi.org/10.1371/journal.pntd.0010145.t003

For the model parameterization, we chose the default regularization multiplier of 1 with 10,000 background points (S2 Table). Bootstrap validation was chosen with 100 replications in which a maximum of 5,000 iterations were allowed for training to reach a 0.00001 convergence threshold. The mean values and confidence interval limits of the outputs were reported and presented as diagrams. Assessment of the contribution of each variable was made by the ‘jack-knife’ method based on a comparison of model performance with (1) the sequential omitting of each variable and (2) modelling with only one variable. The accuracy of the models was estimated by the area under the curve indicator (AUC), which assesses the ability of the model to correctly distinguish between true presence and pseudoabsence data, with AUC>0.7 considered an acceptable value [91].

Statistical analysis of thermal suitability for WNV transmission

Warming temperature and other climatic parameters connected with climate change should be taken into account when evaluating the risk of WNV disease spread [4]. Since the 1990s, both global and regional trends in climate change that may affect the distribution of vector-borne diseases have been observed in Russia [92]. Therefore, we analyzed the change in meteorological conditions that determine the favorability for WNV transmission during the period of official registration of human cases (1997–2019).

First, the suitability of thermal conditions for WNV transmission was assessed based on the accumulation of 109 degree-days (or sum of the effective temperatures, ETs, 0C) which is necessary for the completion of virus extrinsic replication [40]. The calculation of daily degree days was based on a modified sinusoidal method [93,94]. The duration of the season of effective temperatures (SETs) and season of effective infectiveness of mosquitoes (SEI) in days were calculated in addition to the ET. SET is the period during which the average daily temperatures consistently remained above the threshold of 14.3°C, which corresponds to the warm period of the year (usually from March to October in Volgograd). SEI is the period during which viral replication in mosquitoes is possible and mosquitoes can infect humans, mammals, and birds. The duration of SEI depended on how quickly the effective temperatures were accumulated and how long the period with daily temperatures above 14.3°C lasted. This period starts at the same moment as SET but stops earlier due to a shorter time during which there is favorable temperature for full viral replication at the end of the warm season (the infection of mosquitoes is not effective in this case). Before the end of SET, the necessary amount of heat required for the virus replication must be accumulated. The SEI end date is the date after which this amount of heat does not accumulate and effective mosquito infection is not possible.

Second, to evaluate the contemporary changes in the meteorological conditions, we performed a linear trend analysis of ET, SET, SEI, as well as daily temperature, precipitation, and relative humidity for the period of official registration of the disease from 1997 to 2019.

The Mann-Kendall nonparametric statistical test was used to evaluate the significance level of the trends, and Sen’s slope estimator [95] was used to evaluate the rates of changes. These statistical methods can be applied even if the time series do not conform to a normal distribution [95,96]. Finally, Sen’s slope coefficients (k) and minimum confidence levels (p) at which the trends would be statistically significant according to the Mann-Kendall test were obtained for each analyzed parameter.

Third, we analyzed the correlation between human WNV cases during 1997–2019 in Volgograd region and variables of thermal conditions for WNV transmission. The Kendall rank correlation coefficient was used to identify the correlation and generalized linear negative binomial regression model was used to analyze the year-to-year variation.

Statistical analysis was performed in R software using the external packages “trend”, “modifiedmk”, and “MASS”.

Results

Environmental suitability assessment

As a result of the modelling of WNV distribution, environmentally suitable areas were identified. For each scenario under consideration, a model was obtained that was characterized by sufficiently high quality with the AUC ranging from 0.86 to 0.93 (Table 4).

thumbnail
Table 4. Variable contributions according to various spatial models.

https://doi.org/10.1371/journal.pntd.0010145.t004

The analysis of virus detection sites in the environment (Model 1) showed a significant contribution of surface temperature (39.7%), distance from the water bodies (23.8%), and elevation (22.3%) (Table 4). The contribution of water bodies identified by the NDWI was the smallest (14.1%) in this model.

When variables of urban environment were included (Model 2), the relative impact of factors changed. The density of roads was ranked first with a value of 54%. The distance from the water bodies remained important (12.5%) and it was the leading factor among all natural environmental drivers. The contribution of surface temperature and elevation was reduced significantly to 7.6 and 9.6%, respectively. The characteristics determined by NDWI, building and railway density did not make significant contributions to the model (with values of 6.1, 7.9, and 1.4%, respectively).

When the places of probable human infection sites were added to the sites of virus detection in the environment (Model 3), the road density remained the main factor, although its percentage increased slightly to 58.7%. The contribution of surface temperature increased to 17.6% and it was the second leading factor. The third factor was the distance from the water bodies (7.7%). The contributions of other factors were less than 5%.

Model 1 showed that the entire territory of the Volga-Akhtuba floodplain is suitable for WNV transmission (Fig 5A). The same suitability was observed along the Volga and Don rivers and their tributaries. An increase of suitability was identified in the immediate vicinity of water bodies.

thumbnail
Fig 5. Environmental suitability for WNV distribution.

(A) Model 1 based on the data on virus detection sites in the environment as presence data with natural environmental explanatory variables, (B) Model 2 based on the virus detection sites in the environment with natural and urban environmental explanatory variables, and (C) Model 3 based on the data on possible human infection and virus detection sites in the environment with natural and urban environmental explanatory variables. The colour indicates the degree of suitability.

https://doi.org/10.1371/journal.pntd.0010145.g005

The inclusion of urban environment variables in the modelling generally traced the same trends (Fig 5B). Visually, the suitability of areas near water bodies outside Volgograd seemed important in virus transmission. With the inclusion of probable human infection sites, the general tendency of confinement to water bodies remained the same (Fig 5C). At the same time, local areas with highest suitability were allocated to outside the central and the most populated part of the Volgograd city.

Response curves defined the characteristics of each environmental metric (Figs 68). The highest probability of WNV occurrence was positively correlated with road density (Fig 6A). Increasing road density beyond a certain threshold produced an insignificant increase in transmission risk. In contrast, there was a high probability response to building density up to 40 sq. m. within a cell, with a subsequent drop in the response with a higher building density typical of industrial zones (Fig 6B). Optimum surface temperatures corresponded to approximately 33–40°C, which is typical for areas of residential urban development (Fig 6C). It should be noted that the surface temperature in the city is lower than that in the surrounding steppe areas (S3 Fig). The likelihood of virus detection was the greatest near water bodies at a distance of approximately 100 m and in areas with lower elevation of 50–60 m (Figs 7A and 8). The values of the NDWI showed the optimum in areas with herbaceous vegetation (Figs 7B and S4).

thumbnail
Fig 6.

Response curves reflecting the influence of road density (A), building density (B) and LST (C) on the likelihood of the appearance of WNV distribution. The blue area shows the statistical significance of the response curve.

https://doi.org/10.1371/journal.pntd.0010145.g006

thumbnail
Fig 7.

Response curves reflecting the influence of distance to the water bodies (A), NDWI (B) and railway density (C) on the likelihood of the appearance of WNV distribution. The blue area shows the statistical significance of the response curve.

https://doi.org/10.1371/journal.pntd.0010145.g007

thumbnail
Fig 8. Response curves reflecting the influence of elevation on the likelihood of the appearance of WNV distribution.

The blue area shows the statistical significance of the response curve.

https://doi.org/10.1371/journal.pntd.0010145.g008

Assessment of thermal suitability for WNV transmission

The average values of variables used to assess the suitability of temperature conditions for WNV transmission are presented in Table 5.

thumbnail
Table 5. Descriptive statistics and trend values for variables of thermal suitability and daily meteorological variables for WNV transmission between March and October from 1997 to 2019.

https://doi.org/10.1371/journal.pntd.0010145.t005

Trend analysis of thermal conditions for virus replication showed an increase in ET and slight decreases in SET and SEI. Changes were statistically significant for ET only. Statistically significant weakly expressed tendencies in daily temperature increases and relative humidity decreases were also observed. There were no changes in the precipitation rate (Table 5).

A statistically significant correlation between human cases and variables of thermal conditions for WNV transmission was found only for the sum of ET (tau = 0.36, p-value = 0.0224). The regression model showed the statistical significance of the sum of ET to explain yearly variation of cases at p<0.05 level (Fig 9). The distribution of residuals was nearly normal (p-value>0.1 according to Shapiro-Wilk test). SET and SEI varied widely from year to year and there were no differences between outbreak and non-outbreak years.

thumbnail
Fig 9. Observed vs predicted number of WNF cases as per negative binomial regression results with the sum of ET as predictor.

https://doi.org/10.1371/journal.pntd.0010145.g009

Discussion

Since 1999, WNF cases have been reported in Volgograd almost every year, except in 2017. Outbreaks occurred in 1999 (380 cases), 2007 (63 cases), 2010 (521 cases), 2011 (60 cases), 2012 (210 cases), and 2013 (49 cases) [97,98]. Most likely, the number of human infections was much higher than the number of reported cases. A serological survey of the population after the 1999 outbreak revealed that the number of infected humans was 3–10 times higher than the number of laboratory-confirmed cases [99]. Throughout the entire period of WNF case registration in the Volgograd region, the disease prevailed in the urban older age group (over 60), and males had clinical manifestations more often than females [54], which corresponds to the data in other countries. The elderly, immunocompromised, and young individuals were defined as the most affected groups [100103].

The increase in temperature for the warm period of 1997–2019 was weakly manifested, but affected the growth in the sum of ET required for the replication of the virus (Table 5). Though the 22 year period might be a short time span to consider climate change, our results are in line with other studies that were made for a longer period and indicated a current warming trend in southern Russia [104,105]. At the same time, the durations of the SET and SEI were marginally shorter, although these changes were not statistically significant. These data indicate that Volgograd is a territory that has been consistently favorable for the transmission of the virus.

Platonov et al. [54] demonstrated that outbreaks occurred when mean temperatures in May–July were above 21°С. Recent studies of Shocket et al. [106] observed a peak in incidence at a mean of 24°C in the USA. Our study showed similar results; the 2010 outbreak in Volgograd occurred when the March-October mean temperature was 24°C, which was the maximum observed value for the warm period for 1997–2019, whereas the 1999 outbreak occurred at lower temperatures, with a mean of approximately 20.5°C. The sum of ET was 1188°C during the 1999 outbreak and 1555°C during the 2010 outbreak. The correlation between human cases and the sum of ET allowed us to use ET as an indicator for the possibility of outbreaks. In general, a further increase in temperature may facilitate future outbreaks.

For precipitation, no significant changes in trends were found, while a trend towards a decrease in relative humidity was clearly visible and expected in accordance with the trend towards an increase in temperature. In arid territories, there may be an increase in the incidence after heavy rainfall in some seasons. However, such events may not be reflected in long-term trends. The larvae of the species Cx. modestus are very sensitive to hydroclimatic conditions. If reeds, as the main breeding sites in Volgograd, are drying, this can lead to a decrease in the population of this species. Cx. pipiens are more environmentally flexible due to production in man-made sources [107].

The distribution of WNF cases throughout the year is seasonal, with a maximum in August to September in Volgograd [98]. That is, peak clinical disease occurs earlier than the end of the SEI (Table 5). The period of virus replication begins at the end of March (the beginning of SET and SEI) that happens slightly earlier before the massive arrival of migratory birds, including passerine birds, and might indicate local overwintering of the virus. The period of virus replication lasts approximately 5 months. Outbreaks might be facilitated by mild winters, which may contribute to the vector survival [108]. Up to 100%, mortality was observed in Cx. pipiens populations from Russia, which was associated with average monthly temperatures of -3°C [109]. Date from Southern and Southeastern Europe indicates that temperatures below 2°C and above 6°C have a negative association with WNV infections [110].

Most cases of WNF have been reported in the southern part of Volgograd along major water bodies. In the central part of the city, most cases were concentrated away from the Volga River, such as in green areas along Volga’s small tributaries. Outside the city, WNF cases were aggregated in rural areas, including summer cottage settlements along the Volga River, on the banks of large reservoirs, and near the Volga-Don Canal. Detection of the virus in the environment covered a much larger area, including the Volga-Akhtuba floodplain and islands in the Volga delta.

In addition, the Volga-Akhtuba floodplain can be considered a zone of transfer of virus from the delta areas of the neighbouring Astrakhan region and other territories [111]. The Volga delta is located along one of the main seasonal European-Asiatic avian migration routes (http://www.seen-net.eu), and its wetlands are favorable for nesting [27]. Out of 229 bird species recorded from this area, 227 species migrate to winter localities in Africa, the Indian subcontinent, and the Mediterranean [112] which are WNV-endemic areas. Close contacts between birds during intermediate stops and in wintering places can increase the exposure of birds to vectors, thus facilitating the exchange of pathogens and contribute to the further evolution of the virus.

The leading role of LST and elevation in a model that includes only the ‘natural environment’ and ‘virus detection’ variables confirmed the epidemic significance of the Volga-Akhtuba floodplain. WNV infections were less likely in hot steppe areas (due to the low number of larval mosquito habitats), and more likely at lower elevations of the floodplain. The use of LST as one of the predictors of WNV distribution has been confirmed in other studies [79,113]. A negative correlation with altitude has been found in various areas [41,114,115].

The proximity to water bodies is an important concomitant factor. The likelihood of infection decreases as the distance from water bodies increases. The presence of water bodies and wetlands is especially important in the case of urban-rural virus transmission [116]. However, depending on the type of surface water, the influence can be positive or negative as was found in case of southeast USA [117]. Living in the vicinity of slowly moving water sources was statistically associated with increased risk for human infection, whereas living near moderate moving water systems was associated with decreased odds of human infection. Living near bayous lined with vegetation as opposed to concrete channels also showed increased risk of infection [117]. Similar patterns may be typical for the south of Russia.

Our study showed a gradual movement of the pathogen from the riverine recreational strip to urban areas. A buffer zone might be provided by summer cottage settlements that provide a link between rural and urban transmission. In southern France, it has been demonstrated that WNV foci are associated with heterogeneous agricultural areas [118]. They represent a combination of small agricultural plots separated by natural spaces, which may favor more close contacts of competent vectors and reservoir hosts.

The urban environment could significantly increase the epidemic hazard of the territory by providing favorable conditions for the existence of Cx. pipiens populations and synanthropic avian hosts. Increased risk due to urbanization has also been found in other regions [55,119]. The link between infection and urban environment can be explained by the fact that in the urban WNV cycle, mosquitoes gravitate towards residential buildings, built-up areas, and areas with moderate vegetation [55,56]. For example, urban microclimates can differ substantially from the macroclimatic background in mosquito habitats [120]. A small number of ditches and detention ponds could be consistent "superproducers" of Culex larvae [121]. In southern Quebec and eastern Ontario the importance of the edge of vegetation and mixed or paved areas for the ecology of Cx. pipiens were highlighted [122].

Suitable larval habitat could be found in basements of buildings if they are flooded and humans were available as a blood meal source. Basement-bred autogenous Cx. p. molestus mosquitoes may play a separate role in maintaining vertical transmission of the virus within mosquito populations, which has been confirmed by many studies [123125]. However, vertical passage is not efficient for permanent maintenance over many generations without amplification by horizontal transmission, although this might be an additional driver of transmission. In basements, the optimal microclimate for the year-round development of mosquitoes is present [62], thereby maintaining the persistence of the virus during the interepidemic period. Some studies showed that higher temperatures did not increase the transmission rates of biotype molestus [35]. The fact that a large number of WNV cases was reported in Volgograd from the area of low-rise dwellings with small houses with garden plots suggests that the role of this vector may be negligible.

Synanthropic bird species, primarily crows and pigeons, are epiornitic amplifiers in urban environments [126,127]. Dead bird occurrences might be correlated with human WNV risk as was observed in some US studies [40,128130]. However, massive mortality among crows only was observed during the outbreak in Volgograd in 1999 [25,27]. Mortality among the local species Corvus corone was not observed in the Astrakhan region despite the widespread involvement of this species in the circulation of the virus, perhaps due to long-term interaction with this virus [25].

In our models, urban areas were determined by the density of roads. Roads identify the places where more people live and road system components such as culverts, storm drains and roadside ditches easily become perfect larval habitat. An increased risk of WNF in areas with higher road density has also been found in Mississippi [130]. The specificity of the Russian road network when not all settlements have a connection with major roads [131] should be considered. Moreover, in urbanized biotopes of Volgograd, the abundance of Cx. modestus was higher than that in natural environments [31]. The impact of the road density also may indicate that the detection of human cases occurs in frequently visited locations. The results based on virus detection in the environment also suggests that pathogen amplification is efficient in the urban environment.

Limitation

Unfortunately, the lack of spatial information on avian and mosquito populations in the study area did not allow us to include these variables in the model. Entomological information in Russia is collected extremely unevenly and we do not have qualitative data on the species composition and number of ticks and mosquitoes, as well as vector competence, so we used indirect indicators in this analysis. Another limitation was the lack of environmental data sets with a high spatial and temporal resolution for the study area that could allow for the model to be more accurate. This concerns the improvements in mapping of different types of water bodies and urban landscape heterogeneity using remote sensing data. Limited data availability on the exact locations of infection also complicated the analysis. Moreover, there is lack of full epidemiological investigations for all cases of unexplained fevers.

Conclusion

The natural features of the Volgograd region formed the conditions for stable long-term circulation of the West Nile virus. At the same time, the urban environment could pose greater hazards for WNV in humans. Quite favorable conditions for the circulation of the virus emerged here, due to both the involvement of competent hosts and vectors and the formation of favorable environmental conditions for their existence. The private building sector with low-storey houses and garden plots located in the city suburban area and Volga large islands provided connections between urban and rural transmission cycles. If case detection is insufficient and diagnosis is not timely in such territories, adverse consequences may develop.

Supporting information

S1 Fig. Climatological information for study area–Volgograd city and its suburbs.

https://doi.org/10.1371/journal.pntd.0010145.s001

(TIF)

S2 Fig. Population density map.

Contains information from OpenStreetMap and OpenStreetMap Foundation, which is made available under the Open Database License.

https://doi.org/10.1371/journal.pntd.0010145.s002

(TIF)

S3 Fig. LST distribution in Volgograd and its surroundings.

Contains information from OpenStreetMap and OpenStreetMap Foundation, which is made available under the Open Database License.

https://doi.org/10.1371/journal.pntd.0010145.s003

(TIF)

S4 Fig. NDWI distribution in Volgograd and its surroundings.

Contains information from OpenStreetMap and OpenStreetMap Foundation, which is made available under the Open Database License.

https://doi.org/10.1371/journal.pntd.0010145.s004

(TIF)

S1 Table. Variable contributions according to various spatial models.

https://doi.org/10.1371/journal.pntd.0010145.s005

(DOCX)

S2 Table. Model parameterization in MaxEnt software.

https://doi.org/10.1371/journal.pntd.0010145.s006

(DOCX)

References

  1. 1. Johnston B, Conly J. West Nile virus—where did it come from and where might it go? Can J Infect Dis. 2000;4 11:175–8; pmid:18159286
  2. 2. Hubalek Z, Halouzka J. West Nile fever—a reemerging mosquito-borne viral disease in Europe. Emerging Infectious Diseases. 1999;5 5:643–50; pmid:10511520
  3. 3. Rizzoli A, Jimenez-Clavero MA, Barzon L, Cordioli P, Figuerola J, Koraka P, et al. The challenge of West Nile virus in Europe: knowledge gaps and research priorities. Eurosurveillance. 2015;20 20:28–42; pmid:26027485
  4. 4. Paz S. Climate change impacts on West Nile virus transmission in a global context. Phil. Trans. R. Soc. 2015;B3702013056120130561; pmid:25688020
  5. 5. Campbell GL, Marfin AA, Lanciotti RS, Gubler DJ. West Nile virus. Lancet Infectious Diseases. 2002;2(9):519–529; pmid:12206968
  6. 6. Sejvar JJ. West nile virus: an historical overview. Ochsner J. 2003;5(3):6–10. pmid:21765761
  7. 7. Granwehr BP, Lillibridge KM, Higgs S, Mason PW, Aronson JF, Campbell GA, Barrett AD. West Nile virus: where are we now? Lancet Infectious Diseases. 2004;4(9):547–556; pmid:15336221
  8. 8. Chancey C, Grinev A, Volkova E, Rios M. The Global Ecology and Epidemiology of West Nile Virus. Biomed Research International. 2015; pmid:25866777
  9. 9. David S, Abraham AM. Epidemiological and clinical aspects on West Nile virus, a globally emerging pathogen. Infectious Diseases. 2016;48(8), 571–586; pmid:27207312
  10. 10. Sotelo E, Fernandez-Pinero J, Jimenez-Clavero MA. West Nile fever/encephalitis: re-emergence in Europe and the situation in Spain. Enfermedades Infecciosas Y Microbiologia Clinica. 2012;30(2):75–83; pmid:22074879
  11. 11. Napp S, Petric D, Busquets N. West Nile virus and other mosquito-borne viruses present in Eastern Europe. Pathogens and Global Health. 2018;112(5):233–248; pmid:29979950
  12. 12. Zannoli S, Sambri V. West Nile Virus and Usutu Virus Co-Circulation in Europe: Epidemiology and Implications. Microorganisms. 2019;7(7); pmid:31248051
  13. 13. Camp J.V., Nowotny N. The knowns and unknowns of West Nile virus in Europe: what did we learn from the 2018 outbreak? Expert Review of Anti-Infective Therapy. 2020;18(2):145–154; pmid:31914833
  14. 14. Eybpoosh S, Fazlalipour M, Baniasadi V, Pouriayevali MH, Sadeghi F, Vasmehjani AA, et al. Epidemiology of West Nile Virus in the Eastern Mediterranean region: A systematic review. Plos Neglected Tropical Diseases. 2019;13(1); pmid:30695031
  15. 15. Abdullahi IN, Emeribe AU, Ghamba PE, Omosigho PO, Bello ZM, Oderinde BS, et al. Distribution pattern and prevalence of West Nile virus infection in Nigeria from 1950 to 2020: a systematic review. Epidemiology and Health. 2020;42; pmid:33254358
  16. 16. Gubler DJ. The continuing spread of West Nile virus in the Western Hemisphere. Clinical Infectious Diseases. 2007;45(8):1039–1046; pmid:17879923
  17. 17. Petersen LR, Brault AC, Nasci RS. West Nile Virus: Review of the Literature. Jama-Journal of the American Medical Association. 2013;310(3):308–315; pmid:23860989
  18. 18. Kramer LD, Ciota AT, Kilpatrick AM. Introduction, Spread, and Establishment of West Nile Virus in the Americas. Journal of Medical Entomology. 2019;56(6):1448–145; pmid:31549719
  19. 19. Chumakov M, Belyaeva A, Butenko A, Martyanova L. West Nile virus in the USSR. Communication 1. Isolation of West Nile virus strains from ticks Hyalomma plumbeum plumbeum Panz. Endemic viral Infections: Proceedings of the Institute for Poliomyelitis and Viral Encephalitis. 1968:365–72.
  20. 20. May FJ, Davis CT, Tesh RB, Barrett AD. Phylogeography of West Nile virus: from the cradle of evolution in Africa to Eurasia, Australia, and the Americas. J Virol. 2011 Mar;85(6):2964–74. pmid:21159871
  21. 21. Butenko AC, MP. Stolbov DN. Serological and virological examinations in a natural focus of West Nile fever in the Astrakhan region. Voprosy Medicinskoi Virusologii. 1967;1:208–11.
  22. 22. On the sanitary and epidemiological situation in the Russian Federation in 2010: State report. Moscow: Federal Center for Hygiene and Epidemiology of Rospotrebnadzor; 2011: 431. (In Russian).
  23. 23. Putintseva E.V., Smelyansky V.P., Borodai N.V., Alekseichik I.O., Shakhov L.O., Tkachenko G.A. et al. West Nile Fever Worldwide and in the Territory of the Russian Federation in 2016, and Forecast of Epidemic Situation Development in 2017. Problemy Osobo Opasnykh Infektsii [Problems of Particularly Dangerous Infections].2017;1:29–36;; (In Russian).
  24. 24. Putintseva EV, Lipnitsky AV, Alekseev VV, Smelyansky VP, Antonov VA, Manankov VV, et al. Dissemination of the West Nile Fever in the Russian Federation and in the World in 2010. Problemy Osobo Opasnykh Infektsii [Problems of Particularly Dangerous Infections].2011;107. (In Russian).
  25. 25. L’vov DK, Savchenko ST, Alekseev VV, Lipnitsky AV, Pashanina TP. Epidemiological Situation and Prognostication of the West Nile Fever Morbidity in the Territory of the Russian Federation. Problems of Particularly Dangerous Infections. 2008;1:10–12.
  26. 26. Paz S, Semenza JC. Environmental Drivers of West Nile Fever Epidemiology in Europe and Western Asia-A Review. International Journal of Environmental Research and Public Health. 2013;10 8:3543–62; pmid:23939389
  27. 27. Baturin AA VA. Smelyansky VP., Zhukov KV. Chernobay VF. Kolyakina NN. The Role of Birds as Potential Reservoirs of West Nile Virus in the Territory of the Russian Federation. Problems of Particularly Dangerous Infections. 2012;4 414:18–21;
  28. 28. Viruses that are ecologically related to birds. Omsk. inst. source of origin. inf. Omsk. 1971 (in Russian).
  29. 29. L’vov DK. West Nile Fever. Voprosy Virusologii. 2000;2:4–9.
  30. 30. Lelli R, Calistri P, Bruno R, Monaco F, Savini G, Di Sabatino D, et al. West Nile Transmission in Resident Birds in Italy. Transboundary and Emerging Diseases. 2012;59 5:421–8; pmid:22212727
  31. 31. Fyodorova MV, Lopatina YV, Bezzhonova OV, Platonov AY: A complex of bloodsucking mosquitoes (Diptera, Culicidae) in a West Nile fever focus in the Volgograd Region. I. Species composition, seasonal size trend, habitat distribution. vol. 1: Medical parasitology and parasitic diseases; 2007: 41–6.
  32. 32. Fyodorova MV, Lopatina YV, Khutorecskaya NV, Lazorenko VV, Platonov AE. The study of mosquito fauna (Diptera, Culicidae) in Volgograd city in light of the outbreak of West Nile Fever in Volgograd region, 1999. Parazitologiya. 2004;38 3:209–18.
  33. 33. Medlock JM, Snow KR, Leach S. Potential transmission of West Nile virus in the British Isles: an ecological review of candidate mosquito bridge vectors. Medical and Veterinary Entomology. 2005;19 1:2–21; pmid:15752172
  34. 34. Platonov AE, Fedorova MV, Karan LS, Shopenskaya TA, Platonova OV, Zhuravlev VI. Epidemiology of West Nile infection in Volgograd, Russia, in relation to climate change and mosquito (Diptera: Culicidae) bionomics. Parasitology Research. 2008;103:S45–S53; pmid:19030885
  35. 35. Vogels CBF, Göertz GP, Pijlman GP, Koenraadt CJM. Vector competence of European mosquitoes for West Nile virus, Emerging Microbes & Infections.2017;6(1):1–13; pmid:29116220
  36. 36. Savage HM, Ceianu C, Nicolescu G, Karabatsos N, Lanciotti R, Vladimirescu A, et al. Entomologic and avian investigations of an epidemic of West Nile fever in Romania in 1996, with serologic and molecular characterization of a virus isolate from mosquitoes. Am J Trop Med Hyg. 1999;61(4):600–11. pmid:10548295
  37. 37. Lvov DK, Butenko AM, Gromashevsky VL, Kovtunov AI, Prilipov AG, Kinney R, et al. West Nile virus and other zoonotic viruses in Russia: examples of emerging-reemerging situations. Archives of Virology. 2004:85–96. pmid:15119764
  38. 38. Kolodziejek J, Marinov M, Kiss BJ, Alexe V, Nowotny N. The Complete Sequence of a West Nile Virus Lineage 2 Strain Detected in a Hyalomma marginatum marginatum Tick Collected from a Song Thrush (Turdus philomelos) in Eastern Romania in 2013 Revealed Closest Genetic Relationship to Strain Volgograd 2007. Plos One. 2014;9 10; pmid:25279973
  39. 39. Kilpatrick AM, Meola MA, Moudy RM, Kramer LD. Temperature, Viral Genetics, and the Transmission of West Nile Virus by Culex pipiens Mosquitoes. PLoS Pathog. 2008;4(6): e1000092; pmid:18584026
  40. 40. Reisen WK, Fang Y, Martinez VM. Effects of temperature on the transmission of West Nile virus by Culex tarsalis (Diptera: Culicidae). Journal of Medical Entomology. 2006;43 2:309–17; pmid:16619616
  41. 41. Ruiz MO, Chaves LF, Hamer GL, Sun T, Brown WM, Walker ED, et al. Local impact of temperature and precipitation on West Nile virus infection in Culex species mosquitoes in northeast Illinois, USA. Parasites & Vectors. 2010;3; pmid:20302617
  42. 42. Hartley DM, Barker CM, Le Menach A, Niu T, Gaff HD, Reisen WK. Effects of temperature on emergence and seasonality of West Nile virus in California. Am J Trop Med Hyg. 2012;86(5):884–94. pmid:22556092
  43. 43. Andreadis SS, Dimotsiou OC, Savopoulou-Soultani M. Variation in adult longevity of Culex pipiens f. pipiens, vector of the West Nile Virus. Parasitology Research. 2014;113 11:4315–9; pmid:25284257
  44. 44. Blagrove MSC, Caminade C, Waldmann E, Sutton ER, Wardeh M, Baylis M. Co-occurrence of viruses and mosquitoes at the vectors’ optimal climate range: An underestimated risk to temperate regions? Plos Neglected Tropical Diseases. 2017;11 6; pmid:28617853
  45. 45. Moon SA, Cohnstaedt LW, McVey DS, Scoglio CM. A spatio-temporal individual-based network framework for West Nile virus in the USA: Spreading pattern of West Nile virus. Plos Computational Biology. 2019;15 3; pmid:30865618
  46. 46. Wang JF, Ogden NH, Zhu HP. The Impact of Weather Conditions on Culex pipiens and Culex restuans (Diptera: Culicidae) Abundance: A Case Study in Peel Region. Journal of Medical Entomology. 2011;48 2:468–75; pmid:21485391
  47. 47. Shaman J, Stieglitz M, Stark C, Le Blancq S, Cane M. Using a dynamic hydrology model to predict mosquito abundances in flood and swamp water. Emerging Infectious Diseases. 2002;8 1:6–13. pmid:11749741
  48. 48. Groen TA, L’Ambert G, Bellini R, Chaskopoulou A, Petric D, Zgomba M, et al. Ecology of West Nile virus across four European countries: empirical modelling of the Culex pipiens abundance dynamics as a function of weather. Parasites & Vectors. 2017;10; pmid:29070056
  49. 49. Roiz D, Ruiz S, Soriguer R, Figuerola J. Landscape Effects on the Presence, Abundance and Diversity of Mosquitoes in Mediterranean Wetlands. Plos One. 2015;10 6; pmid:26086804
  50. 50. Parselia E, Kontoes C, Tsouni A, Hadjichristodoulou C, Kioutsioukis I, Magiorkinis G, et al. Satellite Earth Observation Data in Epidemiological Modeling of Malaria, Dengue and West Nile Virus: A Scoping Review. Remote Sensing. 2019;11 16;
  51. 51. Jourdain E, Gauthier-Clerc M, Bicout DJ, Sabatier P. Bird migration routes and risk for pathogen dispersion into western Mediterranean wetlands. Emerging Infectious Diseases. 2007;13 3:365–72; pmid:17552088
  52. 52. Tran A, Sudre B, Paz S, Rossi M, Desbrosse A, Chevalier V, et al. Environmental predictors of West Nile fever risk in Europe. International Journal of Health Geographics. 2014;13; pmid:24986363
  53. 53. Chuang TW, Wimberly MC. Remote Sensing of Climatic Anomalies and West Nile Virus Incidence in the Northern Great Plains of the United States. Plos One. 2012;7 10; pmid:23071656
  54. 54. Platonov AE, Tolpin VA, Gridneva KA, Titkov AV, Platonova OV, Kolyasnikova NM, et al. The Incidence of West Nile Disease in Russia in Relation to Climatic and Environmental Factors. International Journal of Environmental Research and Public Health. 2014;11 2:1211–32; pmid:24464233
  55. 55. Ruiz MO, Walker ED, Foster ES, Haramis LD, Kitron UD. Association of West Nile virus illness and urban landscapes in Chicago and Detroit. International Journal of Health Geographics. 2007;6; pmid:17352825
  56. 56. Yoo EH, Chen D, Diao CY, Russell C. The Effects of Weather and Environmental Factors on West Nile Virus Mosquito Abundance in Greater Toronto Area. Earth Interactions. 2016;20:1–22; pmid:30026656
  57. 57. Chuang TW, Hockett CW, Kightlinger L, Wimberly MC. Landscape-Level Spatial Patterns of West Nile Virus Risk in the Northern Great Plains. American Journal of Tropical Medicine and Hygiene. 2012;86 4:724–31; pmid:22492161
  58. 58. Hess A, Davis JK, Wimberly MC. Identifying Environmental Risk Factors and Mapping the Distribution of West Nile Virus in an Endemic Region of North America. Geohealth. 2018;2 12:395–409; pmid:32159009
  59. 59. Eisen L; Bolling BG, Blair CD, Beaty BJ, Moore CG. Mosquito species richness, composition, and abundance along habitat-climate-elevation gradients in the northern Colorado front range. J. Med. Entomol. 2008;45: 800–811; pmid:18714885
  60. 60. Liu H, Weng Q. Environmental Factors and Risk Areas of West Nile Virus in Southern California, 2007–2009. Environmental Modeling & Assessment. 2012;17(4):441–452;
  61. 61. Pecoraro HL, Day HL, Reineke R, Stevens N, Withey JC, Marzluff JM, Meschke JS. Climatic and landscape correlates for potential West Nile virus mosquito vectors in the Seattle region. Journal of Vector Ecology. 2007;32(1):22–28; pmid:17633422
  62. 62. Porshakov A, Yakovlev S, Zakharov K, Matrosov A, Knyazeva T, Kuznetsov A, et al. Role of mosquitoes, Culex pipiens complex, in West Nile fever virus persistence in urbanized biocoenoses of Saratov. Problems of Particularly Dangerous Infections 2014;2 66–68.
  63. 63. Galimzyanov KM, Mirekina EV, Kuryatnikova GK, Polukhina AL, Frank GN, Bedlinskaya NR, Akmaeva LR, Bikbaeva YN, Nikeshina TV. Modern clinical and epidemiological features of West Nile fever in the Astrakhan region. Astrakhan medical journal. 2018; 4:124–130 (in Russian).
  64. 64. Johnson BJ, Munafo K, Shappell L, Tsipoura N, Robson M, Ehrenfeld J, Sukhdeo MVK. The roles of mosquito and bird communities on the prevalence of West Nile virus in urban wetland and residential habitats. Urban Ecosystems. 2012;15(3):513–531; pmid:25484570
  65. 65. Magori K, Bajwa WI, Bowden S, Drake JM. Decelerating Spread of West Nile Virus by Percolation in a Heterogeneous Urban Landscape. Plos Computational Biology. 2011;7(7); pmid:21829332
  66. 66. Kottek M, Grieser J, Beck C, Rudolf B, Rubel F. World map of the Koppen-Geiger climate classification updated. Meteorologische Zeitschrift. 2006;15 3:259–63;
  67. 67. Climate of Russia. St. Petersburg: Gidrometeoizdat; 2001.
  68. 68. Rulev AS, SS. Kosheleva OYu. Assessment of the influence of the hydrological regime of the Volga River on the dynamics of flooding on Sarpinsky Island. Uchenye Zapiski Kazanskogo Universiteta Seriya Estestvennye Nauki. 2017;159 1:139–51.
  69. 69. Yang H, Zhang L, Zhang X, Fang H, Tong X. (). Algorithm of emissivity spectrum and temperature separation based on TASI data. Journal of Remote Sensing. 2011;15:1242–1254.
  70. 70. Stankevich S, Filipovich V, Lubskyi M, Krylova A, Kritsuk S, Brovkina O, et al. (2015). Intercalibration of methods for the land surface thermodynamic temperature retrieving inside urban area by thermal infrared satellite imaging. Ukrainian journal of remote sensing. 2015:14–23 (in Russian).
  71. 71. Guanter L, Alonso L, Moreno J. A method for the surface reflectance retrieval from PROBA/CHRIS data over land: Application to ESA SPARC campaigns. Ieee Transactions on Geoscience and Remote Sensing. 2005;43(12):2908–17.
  72. 72. Chabot-Couture G, Nigmatulina K, Eckhoff P. An Environmental Data Set for Vector-Borne Disease Modeling and Epidemiology. Plos One. 2014;9 4; pmid:24755954
  73. 73. Vancutsem C, Ceccato P, Dinku T, Connor SJ. Evaluation of MODIS land surface temperature data to estimate air temperature in different ecosystems over Africa. Remote Sensing of Environment. 2010;114 2:449–65;
  74. 74. Hay SI. An overview of remote sensing and geodesy for epidemiology and public health application. Advances in Parasitology, Vol 47: Remote Sensing and Geographical Information Systems in Epidemiology. 2000;47:1–35. pmid:10997203
  75. 75. Rouse JW, Haas RH, Scheel JA, Deering DW. Monitoring Vegetation Systems in the Great Plains with ERTS. Proceedings, 3rd Earth Resource Technology Satellite (ERTS) Symposium. 1974;1:48–62.
  76. 76. McFeeters S.K. The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features, International Journal of Remote Sensing. 1996;17(7):1425–1432.
  77. 77. Dudík M, Steven P, Schapire R. Correcting sample selection bias in maximum entropy density estimation. Advances in neural information processing systems. 2005(18): 323–330
  78. 78. Marcantonio M, Rizzoli A, Metz M, Rosa R, Marini G, Chadwick E, et al. Identifying the Environmental Conditions Favouring West Nile Virus Outbreaks in Europe. Plos One. 2015;10 3; pmid:25803814
  79. 79. Kala AK, Atkinson SF, Tiwari C. Exploring the socio-economic and environmental components of infectious diseases using multivariate geovisualization: West Nile Virus. Peerj. 2020;8; pmid:33194330
  80. 80. Bondarenko M, Nieves JJ, Stevens FR, Gaughan AE, Tatem A, Sorichetta A. wpgpRFPMS: Random Forests population modelling R scripts, version 0.1.0. University of Southampton: Southampton, UK. 2020; https://doi.org/10.5258/SOTON/WP00665
  81. 81. Phillips S, Anderson R, Schapire R. Maximum entropy modeling of species geographic distributions. Ecological Modeling. 2006;190:231–59.
  82. 82. Wiese D, Escalante AA, Murphy H, Henry KA, Gutierrez-Velez VH. Integrating environmental and neighborhood factors in MaxEnt modeling to predict species distributions: A case study of Aedes albopictus in southeastern Pennsylvania. PLoS One. 2019;14(10); pmid:31622396
  83. 83. Donalisio M, Peterson A. Environmental factors affecting transmission risk for hantaviruses in forested portions of southern Brazil. Acta Tropica. 2011;119:125–130; pmid:21605537
  84. 84. Elith J, Phillips SJ, Hastie T, Dudík M, Chee YE, Yates CJ. A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 2011;17:43–57;
  85. 85. Merow C, Smith MJ, Silander JA. A practical guide to MaxEnt for modeling species’ distributions: what it does, and why inputs and settings matter. Ecography. 2013;35(10):1058–1069;
  86. 86. van Proosdij A, Sosef M, Wieringa J, Raes N. Minimum required number of specimen records to develop accurate species distribution models. Ecography. 2016;39:542–52;
  87. 87. Støa B., Halvorsen R, Stokland J, Gusarov V. How much is enough? Influence of number of presence observations on the performance of species distribution models. Sommerfeltia. 2019; 39:1–28;
  88. 88. Deb C, Jamir N, Kikon Z. Distribution Prediction Model of a Rare Orchid Species (Vanda bicolor Griff.) Using Small Sample Size. American Journal of Plant Sciences. 2017;8:1388–139;
  89. 89. Wisz MS, Hijmans R, Li J, Peterson AT, Graham C, Guisan A. NCEAS Predicting Species Distributions Working Group. Effects of sample size on the performance of species distribution models. Diversity and Distributions. 2008;14:763–773;
  90. 90. Naimi B, Hamm N, Groen T, Skidmore A, Toxopeus A. Where is positional uncertainty a problem for species distribution modelling. Ecography. 2014;37:191–203;
  91. 91. Mandrekar J. Receiver operating characteristic curve in diagnostic test assessment. Journal of Thoracic Oncology. 2010;5 9:1315–6; pmid:20736804
  92. 92. Mironova VA, Shartova NV, Beljaev AE. et al. Re-introduction of vivax malaria in a temperate area (Moscow region, Russia): a geographic investigation. Malar J. 2020:19:116; pmid:32188468
  93. 93. Allen JC. Modified sine wave method for calculating degree days. Environmental Entomology. 1976;5 3:388–96;
  94. 94. Zou L, Miller SN, Schmidtmann ET. A GIS tool to estimate west nile virus risk based on a degree-day model. Environmental Monitoring and Assessment. 2007;129 1–3:413–20; pmid:17180426
  95. 95. Helsel, RM. DH. Statistical methods in water resources. vol. 49; 1992.
  96. 96. Meals DW. Spooner JD, SA. Harcum JB. Statistical analysis for monotonic trends. Fairfax; 2011.
  97. 97. Monastyrskiy M, Shestopalov N, Akimkin V, YV D. The epidemiological situation of West Nile fever in Volgoggrad region Zhivye i biokosnye sistemy. 2014;9. http://www.jbks.ru/archive/issue-9/article-16. pmid:24895773
  98. 98. Ioannidi E, Muromtseva A, Bozhko V, Kuvshinova T, Viktorov D, Smelyanskiy V. Perculiarities of West Nile fever mnifestation in the Volgograd region. Journal of VolgSMU. 2019;2 70;
  99. 99. Negodenkoa AO, Luchinin DN, Konovalov PS, Pavlyukova OA, Skrynnikova EA, Prilepskaya DR, et al. A screening for serum markers of arbovirus infections in healthy blood donors from the Volgograd region. Infektsiya I Immunitet. 2019;9 5–6:743–9;
  100. 100. Nedry M, Mahon CR. West Nile virus: an emerging virus in North America. Clinical laboratory science: journal of the American Society for Medical Technology. 2003;16(1):43–49. pmid:12587658
  101. 101. Wilson CA, Bale JF. West Nile Virus Infections in Children. Current Infectious Disease Reports. 2014;16(3); pmid:24562539
  102. 102. Ulbert S, Magnusson SE. Technologies for the development of West Nile virus vaccines. Future Microbiology. 2014;9(10):1221–1232; pmid:25405890
  103. 103. Godaert L, Drame M, Roubaud-Baudron C. Emerging viruses in older population Chikungunya, West Nile fever and Dengue. Aging Clinical and Experimental Research. 2021;33(3):723–727; pmid:31741192
  104. 104. Varentsov M, Shartova N, Grischenko M, Konstantinov P. Spatial Patterns of Human Thermal Comfort Conditions in Russia: Present Climate and Trends. Weather, Climate, and Society. 2020;12:629–42;
  105. 105. A report on climate features on the territory of the Russian Federation in 2020. Russian federal service for hydrometeorology and environmental monitoring (Roshydromet). Moscow. 2021. 104 p. (in Russian) https://www.meteorf.ru/upload/pdf_download/doklad_klimat2020.pdf
  106. 106. Shocket MS, Verwillow AB, Numazu MG, Slamani H, Cohen JM, El Moustaid F, et al. Transmission of West Nile and five other temperate mosquito-borne viruses peaks at temperatures between 23 degrees C and 26 degrees C. Elife. 2020;9; pmid:32930091
  107. 107. Putintseva EV, Udovichenko SK, Boroday NV, Nikitin DN, Baturin AА, Molchanova EV, et al. Peculiarities of Epidemiological Situation on the West Nile Fever in the Russian Federation in 2020 and Forecast for Its Development in 2021. Problemy Osobo Opasnykh Infektsii [Problems of Particularly Dangerous Infections]. 2021; 1:63–72. (In Russian).
  108. 108. Smelyanskiy V, Alekseev V, Lipnitskiy A, Antonov V, Ioannidi E, Bozhko V, et al. Specificities of clinical and epidemiological characteristics of West Nile fever in the season of 2010 in the Volgograd region. Infection diseases. 2011;9 4:75–8.
  109. 109. Vinogradova EB. Culex pipiens pipiens mosquitoes: taxonomy, distribution, ecology, physiology, genetics, applied importance and control. Sofia-Moscow: Pensoft Publishers; 2000.
  110. 110. Watts MJ, Sarto I Monteys V, Mortyn PG, Kotsila P. The rise of West Nile Virus in Southern and Southeastern Europe: A spatial-temporal analysis investigating the combined effects of climate, land use and economic changes. One health. 2021;13:100315–. pmid:34485672
  111. 111. Kovalevskaya A, Vasil’kova O, Agapov B, Kouklev E, Safronov V, Shcherbakova S, et al. Risk-Oriented Characteristics of the Current Epidemiological Situation on West Nile Fever in the Astrakhan Region. Problems of Particularly Dangerous Infections. 2019;2:74–8;
  112. 112. Chernobai V. Birds of the Volgograd Oblast. Volgograd: Peremena; 2004. pmid:15124471
  113. 113. Conte A, Candeloro L, Ippoliti C, Monaco F, De Massis F, Bruno R, et al. Spatio-Temporal Identification of Areas Suitable for West Nile Disease in the Mediterranean Basin and Central Europe. Plos One. 2015;10(12). pmid:26717483
  114. 114. Bisanzio D, Giacobini M, Bertolotti L, Mosca A, Balbo L, Kitron U, et al. Spatio-temporal patterns of distribution of West Nile virus vectors in eastern Piedmont Region, Italy. Parasites & Vectors. 2011;4; pmid:22152822
  115. 115. Ben Hassine T, Conte A, Calistri P, Candeloro L, Ippoliti C, De Massis F, et al. Identification of Suitable Areas for West Nile Virus Circulation in Tunisia. Transboundary and Emerging Diseases. 2017;64(2):449–58. pmid:26032967
  116. 116. Liu H, Weng Q, Gaines D. Spatio-temporal analysis of the relationship between WNV dissemination and environmental variables in Indianapolis, USA. International Journal of Health Geographics. 2008;7; pmid:19094221
  117. 117. Nolan MS, Zangeneh A, Khuwaja SA, Martinez D, Rossmann SN, Cardenas V, Murray KO. Proximity of Residence to Bodies of Water and Risk for West Nile Virus Infection: A Case-Control Study in Houston, Texas. Journal of Biomedicine and Biotechnology. 2012; pmid:22315511
  118. 118. Pradier S, Leblond A, Durand B. Land cover, landscape structure, and West Nile virus circulation in southern France. Vector-Borne and Zoonotic Diseases. 2008;8 2:253–63; pmid:18429693
  119. 119. Talbot B, Ardis M, Kulkarni MA. Influence of Demography, Land Use, and Urban Form on West Nile Virus Risk and Human West Nile Virus Incidence in Ottawa, Canada. Vector-Borne and Zoonotic Diseases. 2019;19 7:533–9; pmid:30615572
  120. 120. Wimberly MC, Davis JK, Evans MV, Hess A, Newberry PM, Solano-Asamoah N, et al. Land cover affects microclimate and temperature suitability for arbovirus transmission in an urban landscape. Plos Neglected Tropical Diseases. 2020;14 9; pmid:32956355
  121. 121. Irwin P, Arcari C, Hausbeck J, Paskewitz S. Urban wet environment as mosquito habitat in the upper midwest. Ecohealth. 2008;5(1):49–57; pmid:18648797
  122. 122. Moua Y, Kotchi SO, Ludwig A, Brazeau S. Mapping the Habitat Suitability of West Nile Virus Vectors in Southern Quebec and Eastern Ontario, Canada, with Species Distribution Modeling and Satellite Earth Observation Data. Remote Sensing. 2021;13(9);
  123. 123. Hubalek Z. European experience with the West Nile virus ecology and epidemiology: Could it be relevant for the new world? Viral Immunology. 2000;13 4:415–26; pmid:11192288
  124. 124. Vogels CBF, Fros JJ, Goertz GP, Pijlman GP, Koenraadt CJM. Vector competence of northern European Culex pipiens biotypes and hybrids for West Nile virus is differentially affected by temperature. Parasites & Vectors. 2016;9; pmid:27388451
  125. 125. Ciota AT, Chin PA, Kramer LD. The effect of hybridization of Culex pipiens complex mosquitoes on transmission of West Nile virus. Parasites & Vectors. 2013;6; pmid:24499581
  126. 126. Chanteuse’s SC, Dovas CI, Papanastassopoulou M, Gewehr S, Danis K, Beck C, et al. Evaluation of a West Nile virus surveillance and early warning system in Greece, based on domestic pigeons. Comparative Immunology Microbiology and Infectious Diseases. 2014;37(2):131–41. pmid:24503179
  127. 127. Sharti M, Amouakbari MJ, Pourjabari K, Hashemzadeh MS, Tat M, Omidifar A, et al. Detection of West Nile virus by real-time PCR in crows in northern provinces of Iran. Journal of Vector Borne Diseases. 2020;57(1):37–9. pmid:33818453
  128. 128. Nielsen CF, Reisen WK. West Nile virus-infected dead corvids increase the risk of infection in Culex mosquitoes (Diptera: Culicidae) in domestic landscapes. Journal of Medical Entomology. 2007;44(6):1067–73. pmid:18047208
  129. 129. LaDeau SL, Allan BF, Leisnham PT, Levy MZ. The ecological foundations of transmission potential and vector-borne disease in urban landscapes. Functional Ecology. 2015;29(7):889–901. pmid:26549921
  130. 130. Cooke Wr, Grala K, Wallis R. Avian GIS models signal human risk for West Nile virus in Mississippi. Int J Health Geogr. 2006;5 36; https://doi.org/10.1186/1476-072X-5-36
  131. 131. Mayburov I, Leontyeva Y. Creating targeted fiscal sources of road construction in Russia. Innovation, Entrepreneurship and Sustainable Value Chain in a Dynamic Environment. 2015:1529–1540.