Enhanced West Nile Virus Circulation in the Emilia-Romagna and Lombardy Regions (Northern Italy) in 2018 Detected by Entomological Surveillance

Introduction

The West Nile Virus (WNV) is a flavivirus that circulates in the environment among birds and mosquitoes, but can infect other animals such as equids and humans as dead-end hosts. In humans, virus infection is generally, asymptomatic, but about 20% of infected subjects can develop a febrile illness, which can evolve into a neuro-invasive disease (WNND) in a minority of cases (less than 1% of infected subjects). In addition to symptomatic human cases, the presence of infected, but asymptomatic, blood and organ donors is a primary health issue (1).

After sporadic reports in the '90s in Romania (2) and in Italy (3), as from 2008, the circulation of the virus was recorded at a different extent every season in Europe. In 2018, an alarming recrudescence of the virus circulation was recorded in Europe (4). This recrudescence was also noticed in Northern Italy, including the Emilia-Romagna and Lombardy regions, where an interdisciplinary surveillance plan targeting WNV was set up in 2008 and implemented over the years (5–10).

The plan focuses on the early detection and the pinpointing of the virus circulation area, including those areas where the virus did not circulate the previous season. Over time, the system has shown that it is capable of WNV circulation's early detection at provincial scale (6) and it is economically sustainable (11). Moreover, the plan was able to detect other arboviruses circulating in the surveyed area (12).

The plan has a one-health approach in order to maximize early detection, involving searching for the virus in humans, horses, birds, and mosquitoes. Entomological surveillance is a fundamental part of the plan, characterized by regular scheduling of samplings, and the precise geographic characterization of obtained samples. This permits a fine characterization of the virus circulation on the surveyed area. In this work, we characterize the extraordinary WNV circulation observed in 2018, utilizing data from entomological surveillance obtained in the Emilia-Romagna and Lombardy regions. In particular, we describe the relationship between WNV detection and the abundance of collected mosquitoes, as well as the influence of weather conditions, such as weekly maximum temperatures and number of wet days, on the virus circulation in mosquitoes.

Materials and Methods

Surveyed Area

The surveyed area included the plan area of the Emilia-Romagna and Lombardy regions. This is a large portion of Italy's main plain, the Pianura Padana (or Pianura Padano-Veneta), which includes about 24,000 km² (51.3%) of the total area of 44,700 km². This area encompasses about 10 million inhabitants, with many urban areas, including two of the biggest Italian cities (Milan and Bologna), as well as several important industrial districts. This area has a relevant livestock patrimony and is one of the more important agricultural areas in Italy; distinctive cultivations such as rice fields, vineyards and orchards are highly represented in specific areas.

In general, a strong anthropic modification, with the abundant presence of industrial and urban settlements, characterizes the surveyed area. The rural part of the territory is connoted by an intensive agriculture and animal husbandry, with few hedges, rare scattered trees, and a dense irrigation network. Natural areas are rare, mainly represented by river borders, characterized by riparian vegetation, or re-naturalized areas. The eastern part of the surveyed area ends in the Adriatic Sea, where some of the largest wetlands in Europe are present (Valli di Comacchio and the Po River Delta); pinewood and typical Mediterranean vegetation can be found in this zone.

Trapping sites, with different densities, were selected to cover the entire plain area of the two regions: 95 traps in Emilia-Romagna (surveillance grid of 110 km²) and 39 traps in Lombardy (surveillance grid of 400 km²) (Table 1, Figure 1). These traps worked regularly throughout the surveillance season (seasonal traps). Surveillance was focused on rural areas, semi-natural, rural, or peri-urban locations in Emilia-Romagna, and farms in Lombardy. A supplementary effort was carried out by activating 23 specific traps along the surveillance season in Lombardy, and 9 supplementary traps, which worked once (on August 10th), among selected urban areas in Emilia-Romagna.

TABLE 1

Table 1. Features of surveillance in the two surveyed regions.

FIGURE 1

Figure 1. Map showing the location of traps working throughout the season (circles) and for part of the season (squares), with reference to WNV detections (red), and the reference of the surveyed area on a map of Italy in which the Pianura Padano-Veneta is depicted in gray.

Mosquito Analysis

We collected mosquitoes in fixed geo-referenced stations using attractive traps baited by carbon dioxide, working overnight, roughly from 5:00 p.m. to 9:00 a.m. Every site regularly included in the surveillance was sampled every fortnight. Mosquitoes were identified using morphological characteristics according to the classification key of Severini et al. (13). Due to the used trap model, baited with carbon dioxide, we collected almost exclusively host seeking females, the catching of males was very rare. Identified mosquitoes were counted and pooled according to date, location, and species. The maximum number of specimens per pool was 200 in Emilia-Romagna and 100 in Lombardy. A maximum of 1,000 specimens per species was tested for every sample in Lombardy. Pooled mosquitoes belonging to the genus Culex, and selected pools of other species, were prepared and submitted to biomolecular tests as described elsewhere (14).

In brief, pooled mosquitoes were ground by a vortex mixer with 2-3 4.5-mm-diameter copper-plated round balls (H&N Sport, Münden, Germany) in 2 ml of PBS. After centrifugation, 200 μL aliquots were collected and submitted to biomolecular analysis, for extracting and retrotranscribing RNA. We tested samples by a WNV real-time RT-PCR, according to the method described by Tang et al. (15). This was confirmed by the protocol of Del Amo et al. (16) and by an universal PCR protocol for the identification of flaviviruses (17).

We estimated the number of positive mosquitoes for every site and the day of sampling utilizing the Maximal Likelihood Estimation (MLE) of the PooledInfRate 4.0 excel add-in (18), multiplying the MLE of the infectious rate for the number of collected mosquitoes, and rounding up the obtained result. MLE cannot be estimated when all sample pools test positive. In this case, we approximated the number of infected mosquitoes using the Minimum Infection Rate (MIR), assuming that only one positive mosquito was present in each pool.

We tested the relationship between the fraction of infected mosquitoes collected in each province, with at least two working seasonal traps, and the incidence of WNND cases in the same province through a linear model. Additionally, we checked the hypothesis of normality in the distribution of model residuals through the Shapiro-Wilk test.

Mosquito Abundance and WNV Circulation

We tested the effect of mosquito abundance in a given area (using the catches per trap as a proxy of abundance) on the probability of WNV circulation in the same area in 2018. To guarantee comparable catching efforts, considering the start/end of trapping seasons differed among traps, we used the total number of mosquitoes caught in each trap in July and August (for a total of four sampling sessions) as a measure of mosquito abundance. Traps with less than four sampling sessions in the July-August time window were excluded from the analyses. We built a generalized linear model (GLM) with binomial error distribution and logit link function to estimate the probability of observing a WNV-positive trap in 2018, using the observed number of mosquitoes caught in each trap and the region where the trap was deployed as explanatory variables. In this model, the effect on the GLM of caught mosquitoes (namely, M_OBS) represents the marginal increase in the probability of observing a WNV-positive trap due to an additional mosquito caught (in the logit link).

To test the effect of mosquito abundance on WNV circulation, we built a null model where the probability of each mosquito to test WNV-positive only depends on the region where the mosquito was caught and not on the abundance of mosquitoes in the trap area. These probabilities, namely p_{NULL_ER} and p_{NULL_LO}, were estimated by maximizing the likelihood of obtaining the observed pattern of positive traps in 2018, given the number of mosquitoes caught in each trap in Emilia-Romagna and Lombardy, respectively. We then randomly generated 10,000 synthetic samples where we assigned a WNV-positive/negative status to each mosquito caught in Emilia-Romagna [resp. Lombardy] with probabilities p_{NULL_ER} [resp. p_{NULL_LO}]. Analogously to the observed sample, we estimated for each synthetic sample the effect of the number of mosquitoes caught on the probability of observing a WNV-positive trap (namely, M_NULL) through GLMs. Thus, we generated a distribution of the estimated M_NULL which represents the null hypothesis where mosquito abundance in a given area does not affect the probability of each single mosquito testing WNV-positive. We then tested the observed effect, M_OBS, against the distribution of M_NULL. If M_OBS is significantly higher [resp. smaller] than M_NULL (i.e., it falls outside the 95% range of the M_NULL distribution), then mosquito abundance provides an amplification [resp. dilution] effect on WNV circulation. Analyses were performed with R software.

In order to spatially characterize mosquito abundance, we interpolated the number of mosquitoes sampled in each seasonal trap between July-August (expressed as a natural logarithm) on a map, using the inverse distance weighting method (IDW). The virus circulation's spatial pattern was obtained using positive trap location with the kernel density estimation (KDE) method, with the bandwidth size of 15 km. The area of more intense WNV circulation was estimated through the 50% volume contour of the obtained WNV KDE surface. We used QGIS 3 software (http://www.qgis.org) to perform these analyses.

Meteorological Factors, Mosquito Abundance and Infection Rate

We retrieved daily maximum temperatures and precipitations data from the ECA&D dataset (19) and extracted the data for every surveyed trap site in 2018. We aimed to estimate the association of meteorological parameters with the number of sampled Culex mosquitoes and with the infection rate (namely the proportion of WNV-positive mosquitoes). Since the effect of meteorological parameters on the outcomes of interest might be delayed (lag-effect), we focused on the meteorological parameters recorded up to 4 weeks (lag 1–4) before the night of sampling. Starting from the day before the sampling, we computed the weekly average of maximum temperatures and the weekly number of wet days (number of days in a week with precipitation > 0.5 mm).

Analyses were performed in the framework of the generalized linear mixed models (GLMMs) (20). Specifically, we applied a GLMM including a site-specific and week-specific random effect to account for site and temporal heterogeneity. The abundance of mosquitoes sampled per trap and night was analyzed by applying a linear mixed model with a normal distribution and identity link function. We used the natural logarithm of the count of sampled mosquitoes as the dependent variable and the meteorological parameters recorded up to 4 weeks before the sampling as predictors (lag 1–4). The per trap and night infection rate was analyzed by applying a generalized linear mixed model with Poisson distribution and log link function. We used the estimate of infected mosquitoes as the dependent variable, the number of sampled mosquitoes as offset and the meteorological parameters recorded up to 4 weeks before the sampling as predictors (lag 1–4). The association between temperatures and outcomes was assessed by estimating the coefficient for a unit increase (1°C) in the weekly average of maximum temperatures. The association between precipitations and outcomes was assessed by estimating the coefficient for a unit increase (1 day) in the weekly number of wet days. Analyses were performed with R software (lme4 package).

Results

Mosquito Results

385,293 mosquitoes belonging to 13 species were sampled between June 4 and October 25, (Table 2). Mosquitoes of the Culex genus were grouped into 2,337 pools, of which 232 tested positive for lineage 2 of the WNV. The vast majority of tested pools and all pools which tested positive were composed of Culex (Cx.) pipiens (Table 3). The 134 seasonal traps collected 269,952 Cx. pipiens specimens; of these, 266,854 were tested in 2,254 pools, and 226 were WNV-positive. From among the 134 seasonal traps, 85 collected at least one WNV-positive pool (Figure 1). A total of 1,101 mosquitoes were sampled in the 9 extra traps activated in cities (Bologna, Modena, and Ferrara) in August. All mosquitoes from these samples (823 Cx. pipiens, 171 Aedes (Ae.) albopictus, and 107 Ae. caspius) were tested, and 6 out of the 14 Cx. pipiens pools resulted positive. No positive pools were detected in the other extra traps only activated for part of the season, mainly in northern Lombardy. In order to evaluate the possible involvement of Ae. caspius in virus circulation, we also tested 160 selected pools of this species (for a total of 8,211 specimens) collected between July 17 and August 14. All gave negative results.

TABLE 2

Table 2. Mosquitoes sampled during the 2018 surveillance season.

TABLE 3

Table 3. Tested and WNV-positive mosquito pools for the 2018 surveillance season.

The average of collected mosquitoes per sample at provincial level showed the highest values in the western part of Emilia-Romagna, in Piacenza, Reggio Emilia, and Parma provinces (Table 4). A significant difference in the number of collected mosquitoes was recorded between Lombardy and Emilia-Romagna traps (considering only seasonal traps: μ_{E−R =} 258, CI 233-283, μ_{LOM =} 135, CI 110-159, t = 5.9016, p = 4.6e-9).

TABLE 4

Table 4. Details of entomological surveillance at provincial level during the 2018 surveillance season.

Characterization of Virus Circulation

In 2018, WNV circulation was first detected on June 15 in Emilia-Romagna and on July 17 in Lombardy. We collected the last positive pools on September 6 and October 2 in Lombardy and Emilia-Romagna, respectively (Table 4). The epidemic curve of WNND human cases was postponed with respect to the detection of positive mosquitoes (Figure 2).

FIGURE 2

Figure 2. Average of sampled Cx. pipiens mosquitoes (black line) with CI (dashed lines), positive estimated mosquitoes (white, number above the bar) and number of WNND cases (gray, number above the bar) for Emilia-Romagna (A) and Lombardy (B), during the surveillance.

Only taking seasonal trap data into account, we evaluated the intensity of circulation at the provincial level considering the number of estimated WNV positive mosquitoes on total mosquitoes tested (Figure 3). This highlighted groups of provinces with a different circulation intensity: (i) a group of provinces with a more intense circulation (including central Emilia-Romagna and Pavia); (ii) a group of provinces with an intermediate circulation intensity (in the central and eastern surveillance zones); (iii) a group with a low circulation intensity (Lodi, Brescia, and Piacenza); and (iv) a group in which the virus was undetected (northern Lombardy). Provinces with a more intense circulation often recorded an early and more prolonged virus circulation (Table 4), and a major number of WNND cases.

FIGURE 3

Figure 3. Evaluation of WNV circulation intensity according to number of estimated positive mosquitoes on 1,000 sampled mosquitoes at provincial level (A) and the same on a map (B).

We found the existence of a significant relationship between the rate of infected mosquitoes collected in each province and the incidence of WNND cases in the same province (p = 0.0013). Specifically, we estimated that an increase of 1%0 in the rate of WNV positive pools leads to an additional 1.88 WNND cases per 100,000 inhabitants (Figure 4).

FIGURE 4

Figure 4. Linear model showing the relationship between the fraction of infected mosquitoes at provincial level, and the incidence of WNND cases in the same province.

Relationship Between the Number of Collected Mosquitoes and WNV Incidence

We found a significantly higher probability of finding WNV-positive traps in Emilia-Romagna with respect to Lombardy (p-value = 8.7e-5). We also estimated the marginal increase in the probability of observing a WNV-positive trap due to an additional mosquito caught as M_OBS = 0.00102 (p-value = 0.0021). Figure 5B shows the distribution of the marginal increase in the probability of observing a WNV-positive trap due to an additional mosquito caught obtained with the null model (M_NULL, gray bars) compared to M_OBS (black line). We found that M_OBS is significantly lower than M_NULL (p-value = 0.0232), this means that real number of captured mosquitoes affects less than in the simulations the probability to detect WNV in a trap, suggesting that mosquito abundance in a given site provides a dilution effect on WNV circulation. Figure 5 displays the probabilities of observing a WNV-positive trap as a function of the number of mosquitoes per trap in Lombardy (panel A) and Emilia-Romagna (panel C), estimated through the observed data (black solid lines: best fit, dashed lines: 95% confidence interval) and the null model (white lines: median, gray areas: 95% interval). Panels A and C in Figure 5 show that the probability of finding WNV-positive traps is higher than expected when the abundance of caught mosquitoes is low (i.e., black lines are higher than the white lines), while the opposite is true when the abundance of caught mosquitoes is high (i.e., black lines are lower than the white lines).

FIGURE 5

Figure 5. Probability of observing a WNV-positive trap as a function of the number of mosquitoes per trap in Lombardy (A) and Emilia-Romagna (C) estimated through observed data (black solid lines: best fit, dashed lines: 95% confidence interval) compared with the null model (white lines: median, gray areas: 95% interval). (B) distribution of the marginal increase in the probability to observe a WNV-positive trap due to an additionally caught mosquito obtained with the null model (M_NULL, gray bars) compared to M_OBS (black line).

Despite the probability of finding that WNV increases with the number of mosquitoes (i.e., M_OBS > 0), the probability of a single mosquito being infected decreases as the number of collected mosquitoes increases. Thus, the estimated probability of a single mosquito collected in one site being infected decreased with the increase in the number of collected mosquitoes in the same site.

Fifty percent of the positive traps' KDE revealed a single hot spot in Emilia-Romagna where WNV circulation is more intense (in orange in Figure 6). Overlaying this map onto the map of sampled mosquitoes, directly related to mosquito abundance, highlighted that areas where mosquitoes are most abundant are not necessarily included in this WNV hot spot (Figure 6).

FIGURE 6

Figure 6. Map of the surveyed area with the IDW interpolation of Cx. pipiens females collected between June and August (expressed as natural logarithm) and the hot-spot of WNV circulation (represented by 50% of the KDE of sites with at least one positive pool).

Association Between Weather Data and WNV Circulation

Table 5 and Figure 7 report the β coefficients and 95% CI for the fixed effects of lagged meteorological parameters. Fixed effect coefficients are average estimates from among all sampled traps of the association between a unit increase in the meteorological parameters and the outcomes. β coefficients are estimated under the null hypothesis of absence of association (β = 0), thus a positive estimate of β coefficient suggests that an increase in the meteorological parameter is associated with an increase in the outcome under study (number of sampled mosquitoes or infection rate).

TABLE 5

Table 5. Linear regression coefficient (β) and 95% Confidence Interval (CI) for the logarithm of the number of Culex pipiens Sampled and Lagged Meteorological Parameters.

FIGURE 7

Figure 7. Beta coefficients (black dots) and 95% Confidence intervals (bars) with respect to a 1°C increase in weekly average maximum temperatures in the previous 4 weeks (lag 1–4). Linear regression estimates for the logarithm of the number of sampled Culex pipiens (A); Poisson regression estimates for Infection Rate (B).

As shown in Table 5 and Figure 7, we found an association between the logarithm of the number of sampled mosquitoes and the weekly average of maximum temperatures recorded the week before sampling (β at lag1 : 0.24, 95% CI 0.16–0.31). In contrast, we did not find evidence of association between the logarithm of the count of sampled mosquitoes and the number of wet days at any lag (Table 5).

When analyzing the proportion of WNV positive mosquitoes, we observed a positive association between the infection rate and the weekly average of maximum temperatures recorded two weeks before the sampling (β lag 2: 0.22; 95% CI: 0.02-0.42) (Table 5). No evidence of association was detected between the infection rate and the weekly number of wet days observed in the previous 4 weeks (Table 6).

TABLE 6

Table 6. Poisson regression coefficient for Infection Rate and Lagged Meteorological Parameters.

Discussion

In 2018, in Italy as in other European countries (4), the WNV transmission season began earlier than in previous years. In fact, in Emilia-Romagna, the season started earlier by about one month with respect to previous transmission seasons and the Lombardy region. Circulation intensity and the number of human cases were also the highest observed since the first appearance of WNV in Italy. Starting from the central Provinces of Emilia-Romagna (Modena, Bologna, and Ferrara), WNV circulation progressively extended to both the western and eastern part of the surveyed area covering all the Emilia-Romagna region right up to Lombardy (Pavia, Brescia, Milano, Lodi, and Mantua). However, the virus was undetected in mosquitoes collected in Cremona, Bergamo, Monza Brianza, Varese, and Lecco. The early detection of the virus circulation at provincial level makes it possible to start test blood bags when an area is affected, avoiding infected transfusions.

Reported data confirmed the Cx. pipiens mosquito as the main vector of WNV in the surveyed area. To evaluate the involvement of other species, pools of Ae. caspius and Ae. albopictus were also tested during this survey, without any positive results. The possible involvement of Ae. albopictus as a bridge vector in the WNV cycle was deduced due to WNV positive pools collected in the field in the U.S. (21) and Turkey (22), and the vectorial competence demonstrated in experimental studies (23, 24). The possible vector role of Ae. caspius, since WNV-positive pools of this mosquito were collected in the field in the surveyed area (25), still remains to be demonstrated as we have no clues as to these two mosquitoes' involvement in the WNV cycle in our epidemiological scenario. Independently of intrinsic competence, this was likely due to the host preference of these mosquitoes which, contrary to Cx. pipiens (strictly ornithophilic), are mammophilic feeders. Also, Cx. modestus, a species more competent for WNV than Cx. pipiens in a laboratory study (26), had a marginal role in the monitored area's WNV cycle. This is due to the scarcity of this species, which is strictly linked to the natural environment (13). Results of extra-plan samples in some cities highlighted the intensive viral circulation in 2018, also involving urban areas which are usually less suitable than rural areas for WNV circulation. It is likely that virus circulation in urban areas, also sustained by Cx. pipiens, had an important role in causing the large amount of human cases recorded in 2018.

The surveillance system showed a relevant difference in WNV circulation intensity between the two regions, which was also confirmed by the difference in the number of human cases. Nevertheless, different mosquito sampling efforts between the two regions, and the difference recorded in the GLM, was largely ascribable to the different circulation extent between the two areas. The applied models demonstrate the final number of infected mosquitoes increases when the number of mosquitoes also increases. However, the probability of a single mosquito being infected decreased with mosquito abundance, revealing a non-linear relationship between mosquito abundance and virus circulation. This result agreed with the absence of a direct correlation between mosquito abundance and virus circulation in years with a lower WNV circulation (7). The geographic interpolation of the data utilized to characterize this relationship demonstrated that the rate of infected mosquitoes was not higher where mosquito density was highest. This confirms a less than linear correlation between virus presence and mosquito abundance, temporally and spatially.

The data obtained demonstrated the obvious importance of mosquito abundance but highlighted that this is not the limiting factor for virus circulation in the considered area, which we identified as likely in the number of infected hosts. If we postulate the finite number of infected hosts as a limiting factor, and then the finite capacity to attract mosquitoes, we can hypothesize a dilution effect exerted by the mosquitoes not attracted to them, above a certain number of mosquitoes. In this theoretical scenario, the synchronization of mosquito population peaks and WNV-susceptible birds is fundamental in determining WNV dynamics and the environmental viral load for the season. The observed pattern may also support the hypothesis that increased mosquito density leads to increased avian defensive behaviors, leading to a shift in mosquito feedings toward less defensive hosts, such as mammals which are not WNV-susceptible (27, 28).

A complex interaction between mosquito abundance and susceptible hosts determines the level of virus circulation, which differed in direct relation (namely more mosquitoes, more virus) and is usually retained in vector-borne diseases. We recommend that this result is considered in the construction of epidemiological models and in the cost-efficacy analysis of vector control methods.

In this study, temperature was confirmed as one of the most important factors in determining the virus circulation. We detected an association between temperatures recorded one week and two weeks before mosquito sampling and the mosquito infection rate. Temperature influenced the bionomics of mosquitoes in several ways, for example, by increasing vectorial capacity or shortening the extrinsic incubation period and development times (29). These results are consistent with previous studies that observed a positive association between temperatures and WNV-positive mosquitoes and human cases observed in the following weeks (30–32). This association, already confirmed by several research studies in different ecological settings, is alarming in a global warming scenario, since it implies that increasing temperatures will increase the risk of WNV circulation. The same association can partially explain the recent emergence of the virus, which could be linked to the increase in temperatures recorded in recent years in the surveyed area (33). Interestingly, we observed the difference of one or two weeks in the temperature's ability to influence the number of mosquitoes and the virus infection rate. This observed delay could be explained with the virus' extrinsic incubation period, which consistently lasts about one week at temperatures recorded in the surveyed area (34, 35).

As recorded in another study (31), we were unable to find a consistent link between rainfall (expressed as wet days), mosquito abundance and the mosquito infection rate. An issue related to this process, which can hide the effects of rain, could be the background noise linked to the fall in temperature caused by rainfall. However, water availability under ordinary conditions was not a limiting factor in the study area (due to the widespread presence of irrigating networks, rivers, and water basins).

A human case of WNND is of striking relevance, due to the severity of the symptoms associated with the infection. However, occurrence is rare and not relevant for the environmental persistence of the virus. Despite this, occurrence of WNND was the best available indicator of circulation intensity in people, since there was no organic and standardized system for WNV fever diagnosis. The timing and location of WNND cases are random due to the low rate of WNND cases in infected persons (<1%), the complexity in defining a certain site of patient infection and has human density as one of the main determinants (36), in a particular area. Despite this, we were able to find a clear relationship between the rate of positive mosquitoes and human cases, highlighting the surveillance ability to assess the intensity of virus circulation in a given area at an early stage.

The entomological monitoring described in this work is part of a multidisciplinary surveillance integrating other monitoring systems (wild bird testing and syndromic surveillance of horses). Entomological surveillance often provided the first sign of WNV circulation and, due to the sampling program's standardized regularity, allowed for a fine characterization of the period, area and intensity of circulation. This provides essential data for modulating actions to limit the health risks associated with the circulation of WNV.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

Author Contributions

MCa, PA, RC, DC, MPC, MCh, GG, MT, RB, and DL conceived and designed the surveillance system. MCa, SC, and TT conducted the sampling and mosquito identification. PB and DT performed the molecular analysis. MCa, GM, LB, and AA processed the data and performed statistical analyses. MCa wrote the manuscript. All authors reviewed and approved the final manuscript.

Funding

Data were obtained within the framework of the Regional Plan for Surveillance and Control of Arboviruses−2018 (Emilia-Romagna Region and the Lombardy Region), and the West Nile and Usutu Viruses Integrated National Surveillance and Response plan−2018. Scientific work was partially financed by IZSLER.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

1. Sambri V, Capobianchi M, Charrel R, Fyodorova M, Gaibani P, Gould E, et al. West Nile virus in Europe: emergence, epidemiology, diagnosis, treatment, and prevention. Clin Microbiol Infect. (2013) 19:699–704. doi: 10.1111/1469-0691.12211

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Tsai TF, Popovici F, Cernescu C, Campbell GL, Nedelcu NI. West Nile encephalitis epidemic in southeastern Romania. Lancet. (1998) 352:767–771. doi: 10.1016/S0140-6736(98)03538-7

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Autorino GL, Battisti A, Deubel V, Ferrari G, Forletta R, Giovannini A, et al. West Nile virus epidemic in horses, Tuscany region, Italy. Emerg Infect Dis. (2002) 8:1372–1378. doi: 10.3201/eid0812.020234

PubMed Abstract | CrossRef Full Text | Google Scholar

4. ECDC. Epidemiological update: West Nile virus transmission season in Europe, 2018. Availlable online at: https://ecdc.europa.eu/en/news-events/epidemiological-update-west-nile-virus-transmission-season-europe-2018

Google Scholar

5. Angelini P, Tamba M, Finarelli AC, Bellini R, Albieri A, Bonilauri P, et al. West Nile virus circulation in Emilia-Romagna, Italy: the integrated surveillance system 2009. Euro Surveill. (2010) 15:19547.

PubMed Abstract | Google Scholar

6. Bellini R, Calzolari M, Mattivi A, Tamba M, Angelini P, Bonilauri P, et al. The experience of West Nile virus integrated surveillance system in the Emilia-Romagna region: five years of implementation, Italy, 2009 to 2013. Euro Surveill. (2014) 19:20953. doi: 10.2807/1560-7917.ES2014.19.44.20953

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Calzolari M, Pautasso A, Montarsi F, Albieri A, Bellini R, Bonilauri P, et al. West nile virus surveillance in 2013 via mosquito screening in northern italy and the influence of weather on virus circulation. PLoS ONE. (2015) 10:e0140915. doi: 10.1371/journal.pone.0140915

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Chiari M, Prosperi A, Faccin F, Avisani D, Cerioli M, Zanoni M, et al. West nile virus surveillance in the Lombardy Region, Northern Italy. Transbound Emerg Dis. (2015) 62:343–349 doi: 10.1111/tbed.12375

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Chiari M, Calzolari M, Prosperi A, Perulli S, Faccin F, Avisani D, et al. Surveillance of mosquitoes and selected arthropod-borne viruses in the context of milan EXPO 2015. Int J Environ Res Public Health. (2016) 13:E689. doi: 10.3390/ijerph13070689

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Rovida F, Sarasini A, Campanini G, Percivalle E, Gorini G, Mariani B, et al. West nile virus outbreak in the Lombardy region, northern Italy, summer 2013. Vector Borne Zoonotic Dis. (2015) 15:278–283. doi: 10.1089/vbz.2014.1711

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Paternoster G, Babo Martins S, Mattivi A, Cagarelli R, Angelini P, Bellini R, et al. Economics of one health: costs and benefits of integrated West Nile virus surveillance in Emilia-Romagna. PLoS ONE. (2017) 12:e0188156. doi: 10.1371/journal.pone.0188156

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Calzolari M, Bonilauri P, Bellini R, Albieri A, Defilippo F, Maioli G, et al. Evidence of simultaneous circulation of West Nile and Usutu viruses in mosquitoes sampled in Emilia-Romagna region (Italy) in 2009. PLoS ONE. (2010) 5:e14324. doi: 10.1371/journal.pone.0014324

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Severini F, Toma L, Di Luca M, Romi R. Le zanzare italiane: generalità e identificazione degli adulti (Diptera, Culicidae). Frag Ent. (2009) 41:213–372. doi: 10.4081/fe.2009.92

CrossRef Full Text | Google Scholar

14. Calzolari M, Bonilauri P, Bellini R, Caimi M, Defilippo F, Maioli G, et al. Arboviral survey of mosquitoes in two northern Italian regions in 2007 and 2008. Vector Borne Zoonotic Dis. (2010) 10:875–84. doi: 10.1089/vbz.2009.0176

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Tang Y, Anne Hapip C, Liu B, Fang CT. Highly sensitive TaqMan RT-PCR assay for detection and quantification of both lineages of West Nile virus RNA. J Clin Virol. (2006) 36:177–82. doi: 10.1016/j.jcv.2006.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Del Amo J, Sotelo E, Fernández-Pinero J, Gallardo C, Llorente F, Agüero M, et al. A novel quantitative multiplex real-time RT-PCR for the simultaneous detection and differentiation of West Nile virus lineages 1 and 2, and of Usutu virus. J Virol Methods. (2013) 189:321–7. doi: 10.1016/j.jviromet.2013.02.019

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Scaramozzino N, Crance JM, Jouan A, DeBriel DA, Stoll F, Garin D. Comparison of flavivirus universal primer pairs and development of a rapid, highly sensitive heminested reverse transcription-PCR assay for detection of flaviviruses targeted to a conserved region of the NS5 gene sequences. J Clin Microbiol. (2001) 39:1922–7. doi: 10.1128/JCM.39.5.1922-1927.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Biggerstaff J. PooledInfRate, Version 4.0: a Microsoft® Office Excel© Add-In to compute prevalence estimates from pooled samples. Fort Collins, CO: Centers for Disease Control and Prevention (2009).

Google Scholar

19. Klein Tank AMG, Wijngaard JB, Können GP, Böhm R, Demarée G, Gocheva A, et al. Daily dataset of 20th-century surface air temperature and precipitation series for the European Climate Assessment. Int J of Climatol. (2002) 22:1441–53. doi: 10.1002/joc.773

CrossRef Full Text | Google Scholar

20. Wolfinger R, O'connell M. Generalized linear mixed models a pseudo-likelihood approach. J Stat Comp Simul. (1993) 48:233–43. doi: 10.1080/00949659308811554

CrossRef Full Text | Google Scholar

21. Holick J, Kyle A, Ferraro W, Delaney RR, Iwaseczko M. Discovery of Aedes albopictus infected with west nile virus in southeastern Pennsylvania. J Am Mosq Control Assoc. (2002) 18:131.

PubMed Abstract | Google Scholar

22. Akiner MM, Öztürk M, Başer AB, Günay F, Hacioglu S, Brinkmann A, et al. Arboviral screening of invasive Aedes species in northeastern Turkey: West Nile virus circulation and detection of insect-only viruses. PLoS Negl Trop Dis. (2019) 13:e0007334. doi: 10.1371/journal.pntd.0007334

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Sardelis MR, Turell MJ, O'Guinn ML, Andre RG, Roberts DR. Vector competence of three North American strains of Aedes albopictus for West Nile virus. J Am Mosq Control Assoc. (2002) 18:284–289.

PubMed Abstract | Google Scholar

24. Tiawsirisup S, Platt KB, Evans RB, Rowley WA. Susceptibility of Ochlerotatus trivittatus (Coq.), Aedes albopictus (Skuse), and Culex pipiens (L.) to West Nile virus infection. Vector Borne Zoonotic Dis. (2004) 4:190–7. doi: 10.1089/1530366042162524

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Mancini G, Montarsi F, Calzolari M, Capelli G, Dottori M, Ravagnan S, et al. Mosquito species involved in the circulation of West Nile and Usutu viruses in Italy. Vet Ital. (2017) 53:97–110. doi: 10.12834/VetIt.114.933.4764.2

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Balenghien T, Vazeille M, Grandadam M, Schaffner F, Zeller H, Reiter P, et al. Vector competence of some French Culex and Aedes mosquitoes for West Nile virus. Vector Borne Zoonotic Dis. (2008) 8:589–95. doi: 10.1089/vbz.2007.0266

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Edman JD, Taylor DJ. Culex nigripalpus: seasonal shift in the bird: mammal feeding rates in a mosquito vector of human encephalitis. Science. (1968) 161:67–8. doi: 10.1126/science.161.3836.67

CrossRef Full Text | Google Scholar

28. Thiemann TC, Wheeler SS, Barker CM, Reisen WK. Mosquito host selection varies seasonally with host availability and mosquito density. PLoS Negl Trop Dis. (2011) 5:e1452 doi: 10.1371/journal.pntd.0001452

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Kramer LD. Complexity of virus-vector interactions. Curr Opin Virol. (2016) 21:81–6. doi: 10.1016/j.coviro.2016.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Moirano G, Gasparrini A, Acquaotta F, Fratianni S, Merletti F, Maule M, et al. West Nile Virus infection in Northern Italy: Case-crossover study on the short-term effect of climatic parameters. Environ Res. (2018) 167:544–9. doi: 10.1016/j.envres.2018.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Cotar AI, Falcuta E, Prioteasa LF, Dinu S, Ceianu CS, Paz S, et al. Transmission dynamics of the West Nile virus in mosquito vector populations under the influence of weather factors in the Danube Delta, Romania. EcoHealth. (2016) 13:796–807. doi: 10.1007/s10393-016-1176-y

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Paz S, Malkinson D, Green MS, Tsioni G, Papa A, Danis K, et al. Permissive summer temperatures of the 2010 European West Nile fever upsurge. PLoS ONE. (2013) 8:e56398 doi: 10.1371/journal.pone.0056398

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Antolini G, Pavan V, Tomozeiu R, Marletto V. Atlante Climatico dell'Emilia-Romagna 1961-2015. Edizione 2017. Bologns: Arpae Emilia-Romagna. (2017) pp 32.

Google Scholar

34. Danforth ME, Reisen WK, Barker CM. Extrinsic incubation rate is not accelerated in recent california strains of West Nile Virus in culex tarsalis (Diptera: Culicidae). J Med Entomol. (2015) 52:1083–9. doi: 10.1093/jme/tjv082

CrossRef Full Text | Google Scholar

35. Reisen WK, Fang Y, Martinez VM. Effects of temperature on the transmission of West Nile Virus by Culex tarsalis (Diptera: Culicidae). J Med Entomol. (2006) 43:309–17. doi: 10.1093/jmedent/43.2.309

PubMed Abstract | CrossRef Full Text | Google Scholar

36. 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 Zoonotic Dis. (2019) 19:533–9. doi: 10.1089/vbz.2018.2366

PubMed Abstract | CrossRef Full Text | Google Scholar