In the three-year study, 244,036 mosquito females of 13 species and five genera were caught.
The mosquito species diversity (H) in the agro ecosystem (H = 0.57 ± 0.28, J = 0.77 ± 0.35) was significantly higher (Tab. S1) than in the urban ecosystem (H = 0.41 ± 0.30, J = 0.68 ± 0.53), and no differences in evenness (J) was found in any environment. The dominant species in all the environments was Cx. pipiens, representing 83–98% of total mosquito collections. In addition to Cx. pipiens, Aedes albopictus, a species typical of metropolitan areas in Italy and representing 5–12% of total mosquito collections in this study, was collected in urban and peri-urban areas. In the rural areas, especially near wetlands, Aedes caspius and Aedes vexans, representing 12–13% of total mosquitoes collected, were caught especially after the 2022 flood.
In Bologna municipality, no significant differences in diversity and evenness indices between the monitoring environments were recorded (Tab. S2). However, a greater number of species were found in the peri-urban monitoring sites, where the Shannon indices were higher (Tab. S3).
The pie charts of the faunistic composition of Culicidae collected in different environments are shown in Fig. 3.
The microclimate was significantly different in the two ecosystems: a higher average yearly temperature (+ 1.4°C) and lower relative humidity (RH) (-7.9%) were observed in the urban areas compared to the rural areas, while the differences in rainfall were very low (on average + 9% in the urban areas) (Tab. S4).
Data from www.ornitho.it indicate 88 breeding bird species for the Bologna conurbation in 2010–2023, some of which are considered confirmed, probable, or possibly nesting species. This is a significant increase compared to the 65 species found in 1980–1996 (Tinarelli and Boldreghini 1998). The increase occurred all over the province but particularly in the urban areas around the historic centre (+ 65%; 28 new and 4 disappeared species) due to the presence of small wetlands, settlement of birds of prey, and strengthening and maturation of areas with trees and shrubs. The hooded crow (Corvus cornix), which was present in a few pairs in 1980–1996, and the collared parakeet (Psittacula krameri) and the wood pigeon (Columba palumbus), which were previously absent, have colonised the entire territory of the conurbation in recent years with a significant increase in their populations. The dominant species in Bologna urban areas nowadays are the Columbus, the collared dove, the wood pigeon, the blackbird, the starling, and the hooded crow (S5).
Table 2
Number of Cx. pipiens females (Cx. p.) caught by the traps and pools analysed with RT-PCR flavivirus analysis
Year | Environment monitored | No. of pools of Cx. p. | No. of Cx. p. analysed |
2021 | Urban ecological corridor | 23 | 561 |
| Peri-urban near rural area | 24 | 1080 |
| Peri-urban near corvid roost dormitory | 8 | 137 |
| Urban park | 35 | 513 |
| Rural area near small country town | 54 | 4142 |
| Rural area near wetland | 150 | 18596 |
| Total 2021 | 294 | 25029 |
2022 | Urban ecological corridor | 39 | 1031 |
| Peri-urban near rural area | 22 | 946 |
| Peri-urban near corvid roost dormitory | 22 | 247 |
| Urban park | 43 | 1395 |
| Rural area near small country town | 53 | 3949 |
| Rural area near wetland | 141 | 14817 |
| Total 2022 | 320 | 22385 |
2023 | Urban ecological corridor | 93 | 11595 |
| Peri-urban near rural area | 56 | 7884 |
| Peri-urban near corvid roost dormitory | 29 | 3049 |
| Urban park | 58 | 4073 |
| Rural area near small country town | 140 | 22268 |
| Rural area near wetland | 395 | 64817 |
| Total 2023 | 771 | 113686 |
Grand Total | 1385 | 161,100 |
The RT-PCR flavivirus analysis of 161,100 Cx. pipiens females sorted and pooled in 1385 samples (made up of a maximum of 200 females each) was conducted in the three years of the study (Table 2). Thirty-nine pools tested positive for WNV for lineage 2 (2.82%) and only one for lineage 1 (0.07%).
The Kruskal–Wallis analysis of the number of Cx. pipiens females, infection rates, and vector indices showed significant differences (Table 3e 4) between the six monitored environments. All three parameters were higher in stations activated in rural areas (especially in wet environments) than in those activated in urban areas (mostly those activated in urban parks and ecological corridors). Compared to the rural areas, no statistically significant differences in WNV circulation were observed in the urban ecosystem, but the lowest vector index (VI) and a minor density of Cx. pipiens females were found in urban parks and peri-urban hilly corvid canopy roost areas.
Table 3
Kruskal–Wallis analysis of the number of females collected, infection rates (IR) and vector indices (VI)
Parameters | χ² | Gdl | p |
Females | 45.52 | 5 | < .001 |
IR | 30.12 | 5 | < .001 |
VI | 30.51 | 5 | < .001 |
Table 4
Average number of females collected, infection rates (IR) and vector indices (VI) and Dwass–Steel–Critchlow–Fligner (DSCF) pairwise comparisons (*different letters indicate statistically significant differences, p < 0.05)
Parameters | Area | N | Mean | SD | DSCF* |
Females | Urban ecological corridor | 24 | 97.16 | 305.32 | ac |
| Peri-urban | 24 | 164.25 | 478.52 | ac |
| Peri-Urban hilly roost | 21 | 61.40 | 166.84 | b |
| Urban park | 24 | 39.10 | 82.12 | ab |
| Rural | 24 | 157.90 | 346.49 | cd |
| Rural near wetland | 24 | 190.30 | 269.64 | d |
IR | Urban ecological corridor | 24 | 0.0006 | 0.0023 | a |
| Peri-urban | 24 | 0.0024 | 0.0101 | a |
| Peri-urban hilly roost | 21 | 0.0001 | 0.0003 | a |
| Urban park | 24 | 0.0000 | 0.0000 | a |
| Rural | 24 | 0.0015 | 0.0030 | ab |
| Rural near wetland | 24 | 0.0012 | 0.0018 | b |
VI | Urban ecological corridor | 24 | 0.0215 | 0.0727 | a |
| Peri-urban | 24 | 0.0388 | 0.1317 | a |
| Urban hilly roost | 21 | 0.0235 | 0.1079 | a |
| Urban park | 24 | 0.0000 | 0.0000 | a |
| Rural | 24 | 0.0818 | 0.1569 | ab |
| Rural near wetland | 24 | 0.1484 | 0.2899 | b |
In the three years of the study, 32 confirmed human infections were officially reported, including 24 WNND cases, 2 West Nile fever (WNF) cases, and 6 asymptomatic infections among blood donors (WNI). To guarantee uniformity, only WNND cases were considered in the analysis (Table 5).
Table 5
Human WNV (WNI, WNF and WNND) cases recorded in Bologna province in the period 2021–2023
Year | WNV | Age | Sex | Municipality | Date of initial symptoms |
2021 | WNND | 47 | M | Anzola dell'Emilia | 31 Jul |
| WNND | 80 | M | Medicina | 01 Aug |
2022 | WNND | 73 | M | Imola | 20 Jul |
| WNND | 58 | F | Valsamoggia | 29Jul |
| WNND | 47 | F | Minerbio | 01 Aug |
| WNI | 47 | M | Castel Guelfo di Bologna | 3 Aug |
| WNND | 77 | M | Bologna | 4 Aug |
| WNI | 59 | M | Imola | 16 Aug |
| WNND | 72 | M | Mordano | 20 Aug |
| WNND | 56 | M | Castel Guelfo di Bologna | 28 Aug |
| WNND | 70 | F | Bologna | 02 Sep |
| WNI | 50 | M | Budrio | 06 Sep |
| WNF | 78 | F | Valsamoggia | 15 Sep |
| WNND | 64 | F | Galliera | 08 Oct |
| WNF | 72 | M | Minerbio | 09 Oct |
2023 | WNND | 87 | M | Bologna | 23 Jul |
| WNND | 44 | F | Medicina | 30 Jul |
| WNND | 87 | F | Argelato | 01 Aug |
| WNND | 53 | M | Bologna | 01 Aug |
| WNND | 56 | M | Medicina | 05 Aug |
| WNND | 90 | M | Bologna | 07 Aug |
| WNI | 45 | M | Crevalcore | 07 Aug |
| WNI | 47 | M | Granarolo dell'Emilia | 20 Aug |
| WNND | 59 | F | Budrio | 23 Aug |
| WNND | 65 | M | Medicina | 24 Aug |
| WNND | 76 | F | Minerbio | 28 Aug |
| WNND | 77 | F | Imola | 29 Aug |
| WNND | 68 | F | Sala Bolognese | 02 Sep |
| WNND | 71 | M | Imola | 02 Sep |
| WNND | 86 | F | Sala Bolognese | 10 Sep |
| WNND | 84 | M | Castenaso | 11 Sep |
| WNI | 50 | M | Granarolo dell'Emilia | 13 Oct |
In 2021, a high WNV circulation was observed in mosquitoes collected in early August in rural environments near wetlands, while no and low circulation were observed in urban and rural areas respectively (only two WNND cases were observed in the rural area on 31 July and 1 August, i.e. in week 31).
In 2022, WNV circulation in mosquitoes was recorded only in peri-urban and rural and wetlands in weeks 29 and 30, while the WNND cases were observed from 20 July (week 30) in rural areas and from 1 August (week 32) in urban areas.
In 2023, mosquito pools were infected in all environments monitored, except in urban parks. WNV circulation in mosquitoes started in week 26 up to week 27 in rural areas and in week 28 up to week 29 in urban areas. In urban and rural areas, the first human circulation was observed on 23 July (week 30) and 30 July (week 31) respectively (Fig. 4).
Table 6
Surveillance data recorded in urban and rural areas of Bologna province from 2021 to 2023.
Years | Area | WNND | WNV circulation |
No. cases | Incidence per 100,000 inh. | Cx. pipiens |
Pool No. | Pool + | %pool + | VI Average 25–41 w | VI Max |
2021 | Urban | 0 | 0.00 | 90 | 0 | 0.00% | 0.00 | 0.00 |
Rural | 2 | 0.20 | 205 | 11 | 5.37% | 0.12 | 1.37 |
2022 | Urban | 2 | 0.51 | 126 | 1 | 0.79% | 0.01 | 0.44 |
Rural | 6 | 0.59 | 194 | 6 | 3.09% | 0.06 | 0.25 |
2023 | Urban | 3 | 0.77 | 236 | 4 | 1.69% | 0.05 | 0.49 |
Rural | 11 | 1.08 | 536 | 18 | 3.36% | 0.17 | 0.57 |
The WNND incidence was higher in rural areas than in urban areas (on average 1.01 ± 0.59SD and 0.43 ± 0.32SD cases per 100,000 inh. respectively), and differences were observed in the VI values (VI = 0.12 ± 0.04 SD and VI = 0.02 ± 0.02SD respectively), showing a higher risk in rural areas (Table 6).