# Spatio-temporal trends and risk factors affecting west nile virus and related flavivirus exposure in spanish wild ruminants

Ignacio Garcia-Bocanegra, Jorge Paniagua, Ana V. Gutierrez-Guzman, Sylvie Lecollinet, Mariana Boadella, Antonio Arenas-Montes, David Cano-Terriza, Steeve Lowenski, Christian Gortazar, Ursula Hoefle

► To cite this version:

Ignacio Garcia-Bocanegra, Jorge Paniagua, Ana V. Gutierrez-Guzman, Sylvie Lecollinet, Mariana Boadella, et al. Spatio-temporal trends and risk factors affecting west nile virus and related flavivirus exposure in spanish wild ruminants. BMC Veterinary Research, 2016, 12, 10.1186/s12917-016-0876-4. hal-01606863

HAL Id: hal-01606863

https://hal.science/hal-01606863

Submitted on 28 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

Open Access

<figure><img src="image_2.png" /></figure>

CrossMark

Spatio-temporal trends and risk factors
affecting West Nile virus and related
flavivirus exposure in Spanish wild
ruminants

Ignacio García-Bocanegra¹, Jorge Paniagua¹, Ana V. Gutiérrez-Guzmán², Sylvie Lecollinet³, Mariana Boadella⁴,
Antonio Arenas-Montes¹, David Cano-Terriza¹, Steeve Lowenski³, Christian Gortázar² and Ursula Höfle²*

<figure><img src="image_3.png" /></figure>

**Abstract**

**Background:** During the last decade, the spread of many flaviviruses (Genus *Flavivirus*) has been reported, representing an emerging threat for both animal and human health. To further study utility of wild ruminant samples in West Nile virus (WNV) surveillance, we assessed spatio-temporal trends and factors associated with WNV and cross-reacting flaviviruses exposure, particularly Usutu virus (USUV) and Meaban virus (MBV), in wild ruminants in Spain. Serum samples from 4693 wild ruminants, including 3073 free-living red deer (*Cervus elaphus*), 201 fallow deer (*Dama dama*), 125 mouflon (*Ovis aries musimon*), 32 roe deer (*Capreolus capreolus*) and 1262 farmed red deer collected in 2003–2014, were screened for WNV and antigenically-related flavivirus antibodies using a blocking ELISA (bELISA). Positive samples were tested for neutralizing antibodies against WNV, USUV and MBV by virus micro-neutralization tests.

**Results:** Mean flavivirus seroprevalence according to bELISA was 3.4 ± 0.5 % in red deer, 1.0 ± 1.4 % in fallow deer, 2.4 ± 2.7 % in mouflon and 0 % in roe deer. A multivariate logistic regression model revealed as main risk factors for seropositivity in red deer; year (2011), the specific south-coastal bioregion (bioregion 5) and presence of wetlands. Red deer had neutralizing antibodies against WNV, USUV and MBV.

**Conclusions:** The results indicate endemic circulation of WNV, USUV and MBV in Spanish red deer, even in areas without known flavivirus outbreaks. WNV antibodies detected in a free-living red deer yearling sampled in 2010, confirmed circulation this year. Co-circulation of WNV and USUV was detected in bioregions 3 and 5, and of WNV and MBV in bioregion 3. Sampling of hunted and farmed wild ruminants, specifically of red deer yearlings, could be a complementary way to national surveillance programs to monitor the activity of emerging flaviviruses.

**Keywords:** West Nile virus, Usutu virus, Meaban virus, Red deer *Cervus elaphus*, Risk factors, Wild ruminants, Spain

**Background**

The distribution of vector-borne flaviviruses (family *Flaviviridae*) in the world has substantially increased over the last decades. During this period, many flavivirus infections have become a major public health concern due to continuous and growing reporting of outbreaks in humans [1]. Flaviviruses are mainly transmitted within an enzootic cycle involving ornithophilic mosquitoes or ticks as

competent vectors, as well as wild birds as the main ampli-
fying hosts in the wild. Most mammalian species including
humans are considered dead-end or incidental hosts,
because they can get infected but are not thought to be
able to transmit the viruses.

During the last few years, six flaviviruses, including West
Nile virus (WNV), Usutu virus (USUV), tick-borne ence-
phalitis virus (TBEV), Bagaza virus (BAGV), Meaban virus
(MBV) and louping-ill virus (LIV), have been detected in
Europe [2, 3]. Five of them have circulated in Spain in the
last decade. WNV exposure has been documented in
mosquitoes [4], wild birds [5, 6] and different mammalian

<footnote>* Correspondence: ursula.hofle@uclm.es

²Instituto de Investigación en Recursos Cinegéticos IREC, (CSIC-UCLM-JCCM),
Ciudad Real, Spain

Full list of author information is available at the end of the article</footnote>
species, including humans [7], horses [8], dromedary camels (*Camelus dromedarios*) [9], wild boar (*Sus scrofa*) and Iberian pigs and red foxes (*Vulpes vulpes*) [10]. Clinical disease and mortality associated with WNV infection has also been detected in wild birds, horses and humans in this country [11, 12]. USUV has been detected in mosquitoes [4, 13] and in both migratory and resident birds in Spain [14, 15]. Mortality associated with Louping-ill virus (LIV) infection was also detected in sheep and goats in northern Spain [16] and suspected in chamois (*Rupicapra pyrenaica*) [17]. An unusually high mortality due to Bagaza virus (BAGV) infection was also confirmed in free-living game birds in south-western Spain in 2010 [18]. Finally, Meaban virus (MBV) has been found in both yellow-legged gulls (*Larus michaelis*) and ticks (*Ornithodoros maritimus*) in north-eastern Spain [3].

Wild and domestic artiodactyls can be useful sentinel species for monitoring flavivirus activity [19–21]. Antibodies against St. Louis encephalitis virus (SLEV) and WNV have been found in white-tailed deer (*Odocoileus virginianus*) from the United States [22, 23], and against TBEV and WNV in different game species from the Czech Republic [24]. In addition, fatal cases of WN fever were reported in white-tailed deer and reindeer (*Rangifer tarandus*) in North America [25, 26]. Serosurveillance on WNV and related flaviviruses has also been performed in wild and domestic Spanish animals such as wild boar, red deer (*Cervus elaphus*), Iberian pigs, red foxes, and others [8, 10]. Accordingly, 0.2 % yearling red deer from south-western Spain had antibodies against WNV or cross-reacting flaviviruses [27].

As red deer is an important game species in Spain [28] and is also frequently farm raised, it could be an easily accessible, cost-effective species to use as a complementary tool to the national surveillance programs to monitor the activity of mosquito-borne flaviviruses [27]. While the analysis of samples from yearlings is useful for continuous surveillance and early warning systems, the (retrospective) analysis of serum samples of individuals of all age classes is a useful tool to explore temporal and spatial trends of flavivirus activity in a given region [23].

The aim of this study was to monitor seroprevalence of WNV and antigenically-related mosquito and tick-borne flaviviruses, particularly USUV and MBV, in wild ruminants in Spain, and, using red deer as the most distributed, abundant and readily available species, to assess the spatial-temporal trends and risk factors associated with the exposure to these flaviviruses in this species.

# Methods

## Study area

Samples of 4693 wild ruminants, including 3073 free-living and 1262 farmed red deer, 201 fallow deer (*Dama dama*), 125 mouflon (*Ovis aries musimon*) and 32 roe deer (*Capreolus capreolus*) were collected between 2003 and 2014 from

hunting estates across Spain. Red deer were sampled in 130 hunting areas located throughout the five geographical bioregions of Spain previously defined for the national wildlife disease monitoring program [29, 30]: (1) Atlantic coast, characterized by a wet, temperate climate and abundant rainy seasons (32 sampling areas); (2) cereal plains, in which agriculture with cereal crops are dominant (7 sampling areas); (3) Continental Mediterranean ecosystem, with cold winters, hot dry summers, and rainy seasons in autumn and spring (70 sampling areas); (4) Interior mountains, consisting of small mountain chains with a continental Mediterranean climate (9 sampling areas); (5) and south coast, with a humid coastal Mediterranean climate, warm humid winters and hot dry summers (12 sampling areas). Samples of roe deer, fallow deer, and mouflon, were obtained only from bioregions 3 and 5 (south-central Spain).

## Sampling

Blood samples were taken from the thoracic cavity of freshly killed hunted animals or during health inspection in the slaughterhouse. Samples of farmed red deer were obtained by puncture of the jugular vein using a sterile collection system during health inspections. All samples of free-living ruminants were grouped by season and habitat type (sites associated with large permanent wetlands i.e., rivers, lagoons or lakes, and areas without wetlands or large water bodies). Samples were classified by age according to dentition patterns (yearlings, sub-adults, adults) [31]. The animals sampled were classified according to gender and status (farmed or free-living). When kidney samples were available from red deer, the right kidney fat index (KFI) was obtained as an indicator of body condition [32]. Upon arrival at the laboratory, blood samples were centrifuged for 15 min at 1,800X g for serum separation and serum was stored at -20 °C until analysis.

## Laboratory analysis

All ruminants were tested for antibodies against an epitope of the WNV pre-membrane-envelope (prM-E) protein shared with other viruses of the Japanese encephalitis serocomplex. A commercially available blocking enzyme-linked immunosorbent assay (bELISA 10.WNV.K3 INGEZIM West Nile COMPAC®, Ingenasa, Madrid, Spain) was used in accordance with the manufacturer's recommendations. The bELISA was used as a serological screening tool and bELISA-positive and doubtful sera were confirmed by micro virus neutralization test (VNT) for the detection of specific neutralizing antibodies against WNV (Is98 strain, lineage 1) and USUV (It12 strain). Additionally, given the possibility of cross-reaction with antigenically-related flaviviruses not included within the Japanese encephalitis serocomplex such as Meaban virus (MBV), bELISA-positive and doubtful samples were also tested by VNT against this flavivirus (MBV; Brest ART707 strain). VNTs were performed
as previously described [3, 33]. Samples that showed neutralization and absence of cytopathic effect at dilutions ≥ 10 for WNV and USUV and ≥ 20 for MBV were considered positive. Interpretation of results was based on comparison of VNT titers obtained in parallel against the three flaviviruses. The neutralizing immune response observed was considered specific when VNT titers for a given virus was ≥ 4-fold higher than titers obtained for the other viruses. Samples showing VNT titers differences ≤ 2-fold between the viruses examined were considered positive for flaviviruses but not conclusive for any specific virus.

## Statistical analysis

The prevalence of antibodies against antigenically-related flaviviruses was estimated from the ratio of positives to the total number of samples, with the exact binomial confidence intervals of 95 %. Differences between species were analyzed using a Pearson's Chi-square test and a Fisher's exact test (when observations/category were <6). Due to the absence of seropositivity in roe deer and to the limited number of seropositive fallow deer and mouflon, the associated risk factors were only analyzed for red deer, the most widely distributed and abundant wild ruminant species in Spain. A Chi-square or Fisher's exact test were used to test the relevance of the explanatory variables (age class, gender, season, year, status, bioregion and wetland area) to the risk of an red deer being exposed to flaviviruses. Covariates correlated with a P-value < 0.20 in the bivariate analysis were included for further analysis. Biologically plausible confounding factors were assessed using Mantel-Haenszel analysis and confounding was considered to be potentially significant if odd ratios (ORs) shifted appreciably. Variables that altered the coefficients for the independent variables of interest by 30 % or more were removed from the model and were classified as confounding factors. Finally, a multiple logistic regression analysis [34] was performed including risk factors potentially associated with related flavivirus exposure (likelihood-ratio Wald's test, P < 0.05). The goodness of fit was assessed using the Hosmer-Leme-show goodness-of-fit test. SPSS 22.0 software (IBM Corp., Armonk, NY, USA) was used for statistical analyses.

## Results

A total of 153 out of the 4693 (3.3 %, CI<sub>95%</sub>: 2.7–3.8 %) wild ruminants were seropostive against antigenically-related flaviviruses using bELISA. A mean seroprevalence of 3.4 ± 0.5 % (148/4335) was obtained in red deer, 1.0 ± 1.4 % (2/201) in fallow deer, and 2.4 ± 2.7 % (3/125) in mouflon. No seropositivity was detected in the 32 roe deer tested. Statistically significant differences in seroprevalence among species were not observed.

Five (year, bioregion, wetland area, season and status) out of seven variables were selected (P < 0.20) from the bivariate analysis (Table 1). At least one seropositive red deer was

<table><caption>**Table 1** Seroprevalence to WNV and antigenically-related flaviviruses in wild ruminants in Spain</caption><thead><tr><th>Variable</th><th>Categories</th><th>No. examined<sup>a</sup>/positive (%)</th><th>P-value</th></tr></thead><tbody><tr><td rowspan="4">Species</td><td>Red deer</td><td>4335/148 (3.4)</td><td rowspan="4">0.174</td></tr><tr><td>Fallow deer</td><td>201/2 (1.0)</td></tr><tr><td>Mouflon</td><td>125/3 (2.4)</td></tr><tr><td>Roe deer</td><td>32/0 (0.0)</td></tr><tr><td rowspan="3">Age class</td><td>Juveniles</td><td>56/1 (1.8)</td><td rowspan="3">0.628</td></tr><tr><td>Sub-adults</td><td>221/8 (3.6)</td></tr><tr><td>Adults</td><td>607/26 (4.3)</td></tr><tr><td rowspan="2">Sex</td><td>Females</td><td>949/24 (2.5)</td><td rowspan="2">0.217</td></tr><tr><td>Males</td><td>893/29 (3.2)</td></tr><tr><td rowspan="5">Bioregion</td><td>1</td><td>437/1 (0.2)</td><td rowspan="5">&lt;0.001</td></tr><tr><td>2</td><td>506/3 (0.6)</td></tr><tr><td>3</td><td>2164/61 (2.8)</td></tr><tr><td>4</td><td>230/4 (1.7)</td></tr><tr><td>5</td><td>998/79 (7.9)</td></tr><tr><td rowspan="10">Year</td><td>2003</td><td>502/9 (1.8)</td><td rowspan="10">&lt;0.001</td></tr><tr><td>2004</td><td>179/1 (0.6)</td></tr><tr><td>2005</td><td>338/5 (1.5)</td></tr><tr><td>2006</td><td>446/12 (2.7)</td></tr><tr><td>2007</td><td>664/13 (2.0)</td></tr><tr><td>2008</td><td>551/15 (2.7)</td></tr><tr><td>2009</td><td>555/21 (3.8)</td></tr><tr><td>2010</td><td>576/20 (3.5)</td></tr><tr><td>2011</td><td>250/44 (17.6)</td></tr><tr><td>2012</td><td>74/1 (1.4)</td></tr><tr><td>2013</td><td>140/5 (3.6)</td></tr><tr><td>2014</td><td>60/2 (3.3)</td></tr><tr><td rowspan="4">Season</td><td>Autumn</td><td>1754/52 (3.0)</td><td rowspan="4">0.017</td></tr><tr><td>Spring</td><td>377/10 (2.7)</td></tr><tr><td>Summer</td><td>503/29 (5.8)</td></tr><tr><td>Winter</td><td>1701/57 (3.4)</td></tr><tr><td rowspan="2">Wetland area</td><td>Presence</td><td>1634/109 (6.7)</td><td rowspan="2">&lt;0.001</td></tr><tr><td>Absence</td><td>2701/39 (1.4)</td></tr><tr><td rowspan="2">Status</td><td>Free-living</td><td>3073/74 (2.4)</td><td rowspan="2">&lt;0.001</td></tr><tr><td>Farmed</td><td>1262/74 (5.9)</td></tr></tbody></table>

detected in every year analyzed, with a significantly higher seroprevalence detected in 2011 (17.6 %; 44/250; P < 0.001). Forty-one out of the 130 (31.5 %) sampling areas presented at least one seropositive animal. The seropositivity by sampling areas was 3.1 % (1/32) in bioregion 1, 42.9 % (3/7) in bioregion 2, 44.3 % (31/70) in bioregion 3, 11.1 % (1/9) in bioregion 4 and 41.7 % (5/12) in bioregion 5. Even though seropositivity was detected in all bioregions, considering farmed and free-living red deer together, a significantly

<footnote><sup>a</sup>Missing values excluded</footnote>
higher seroprevalence was found in bioregion 5 (7.9%; 79/998) compared to the other bioregions (*P* < 0.001), and in bioregion 3 (2.8%; 61/2164) compared to bioregion 1 (0.2%; 1/437) (*P* < 0.001). Significantly higher seropositivity was also detected in bioregion 5 (3.0%; 8/269) (*P* = 0.003), bioregion 4 (3.0%; 4/134) (*P* = 0.012) and bioregion 3 (2.9%; 59/2017) (*P* < 0.001) compared to bioregion 1 (0.2%; 1/437) when seroprevalence was compared only in free-living red deer. A significantly higher seroprevalence was obtained in summer (5.8%; 29/503) as compared to autumn (3.0%; 52/1754) (*P* = 0.003) and winter (3.4%; 57/1701) (*P* = 0.012). Seroprevalence was significantly higher in farmed 5.9% (74/1262) than in free-living 2.4 ± 0.5% (74/3073) red deer (*P* < 0.001). Seropositivity was confirmed in three out of the five red deer farms tested, with the mean seroprevalence ranging from 0.4% (1/243) to 9.7% (71/729). In the red deer farm with the highest seropositivity, a significantly higher mean seroprevalence was found in autumn (18.5%; 20/108) compared to spring (4.6%; 6/131; *P* = 0.002) and summer (8.4%; 26/309; *P* = 0.012), but not in winter (10.5%; 19/181; *P* = 0.113). The final multivariate logistic regression model showed that the main risk factors associated with exposure to related flaviviruses in red deer were: year (2011), bioregion (5) and wetland areas (presence) (Table 2).

<table><caption>**Table 2** Logistic regression model of potential risk factors associated with seroprevalence to WNV and antigenically-related flaviviruses in red deer in Spain</caption><thead><tr><th>Variable</th><th>Categories</th><th>B</th><th>P-value</th><th>OR</th><th>95 % CI</th></tr></thead><tbody><tr><td rowspan="14">Year</td><td>2003</td><td>a</td><td>a</td><td>a</td><td>a</td></tr><tr><td>2004</td><td>-1.174</td><td>0.271</td><td>0.309</td><td>0.038 2.496</td></tr><tr><td>2005</td><td>-0.028</td><td>0.962</td><td>0.972</td><td>0.313 3.020</td></tr><tr><td>2006</td><td>0.294</td><td>0.525</td><td>1.342</td><td>0.542 3.322</td></tr><tr><td>2007</td><td>0.223</td><td>0.627</td><td>1.250</td><td>0.509 3.072</td></tr><tr><td>2008</td><td>0.430</td><td>0.332</td><td>1.538</td><td>0.645 3.669</td></tr><tr><td>2009</td><td>0.509</td><td>0.219</td><td>1.663</td><td>0.739 3.743</td></tr><tr><td>2010</td><td>0.386</td><td>0.355</td><td>1.471</td><td>0.649 3.335</td></tr><tr><td>2011</td><td>1.991</td><td>&lt;0.001</td><td>7.320</td><td>3.404 15.739</td></tr><tr><td>2012</td><td>-0.741</td><td>0.487</td><td>0.477</td><td>0.059 3.851</td></tr><tr><td>2013</td><td>0.077</td><td>0.892</td><td>1.081</td><td>0.353 3.306</td></tr><tr><td>2014</td><td>0.608</td><td>0.447</td><td>1.837</td><td>0.383 8.822</td></tr><tr><td rowspan="5">Bioregion</td><td>1</td><td>a</td><td>a</td><td>a</td><td>a</td></tr><tr><td>2</td><td>0.602</td><td>0.605</td><td>1.826</td><td>0.186 17.946</td></tr><tr><td>3</td><td>1.913</td><td>0.060</td><td>6.777</td><td>0.920 49.903</td></tr><tr><td>4</td><td>1.784</td><td>0.115</td><td>5.954</td><td>0.648 54.709</td></tr><tr><td>5</td><td>2.192</td><td>0.034</td><td>8.957</td><td>1.180 67.990</td></tr><tr><td rowspan="2">Wetland area</td><td>Absence</td><td>a</td><td>a</td><td>a</td><td>a</td></tr><tr><td>Presence</td><td>1.097</td><td>&lt;0.001</td><td>2.995</td><td>1.888 4.751</td></tr></tbody></table>

Testing by VNT was possible in 140 out of the 153 bELISA seropositive wild ruminants, while 13 samples (11 from red deer and 2 from fallow deer) could not be analyzed by VNT due to serum cytotoxicity. Of these, 25 (22 from red deer and 3 from mouflon) showed negative results against the three flaviviruses tested using VNT. Specific antibodies against WNV, USUV and MBV were confirmed in 103 (69.7%; bioregions 2, 3, 4 and 5), 4 (2.7%; bioregions 1, 3 and 5) and 2 (1.4%; bioregion 3) red deer, respectively (Table 3, Fig. 1). Six red deer (five from bioregion 3 and one from bioregion 5) showed ≤ 2-fold VNT titer differences between WNV and USUV and were therefore considered positive for other related flaviviruses but not conclusive for any of these viruses. Taking into account possible co-infections, the overall antibody prevalence ranged between 2.2 and 2.4% for WNV, 0.1 and 0.2% for USUV and 0.04% for MBV. Specific WNV, USUV and MBV antibodies were confirmed in 31, four and one sampling areas, respectively. WNV and USUV and WNV and MBV co-circulation (different individuals from the same sampling area with specific antibodies against WNV, USUV or MBV) were detected in two (bioregions 3 and 5) and one (bioregion 3) sampling areas, respectively (Fig. 1).

Discussion

The overall prevalence of antibodies against WNV and antigenically-related flaviviruses detected in wild ruminants (3.3%) in the present study was consistent with that recently observed in captive zoo artiodactyls in Spain [15]. Similar seroprevalences were also found in white-tailed deer in the United States [22, 23]. In contrast, a limited number of wild ruminants from south-western Spain tested during surveillance of WNV outbreaks were negative for flavivirus antibodies using the same bELISA [8]. In our study, no antibodies were detected in roe deer, potentially due to the small sample size, while a similar study in the Czech Republic reported flavivirus exposure in roe deer, mouflons, fallow deer and red deer [24]. Remarkably, wild ruminants and specifically red deer showed lower seroprevalence to

<table><caption>**Table 3** Virus neutralization test (VNT) titers obtained in parallel against West Nile virus (WNV), Usutu virus (USUV) and Meaban virus (MBV) from 115 sera positive by bELISA</caption><thead><tr><th>VNT Titers</th><th>WNV</th><th>USUV</th><th>MBV</th><th>WNV and/or USUV</th></tr></thead><tbody><tr><td>10</td><td>28</td><td>3</td><td></td><td>4 (WNV); 3 (USUV)</td></tr><tr><td>20</td><td>31</td><td></td><td>2</td><td>1 (WNV); 2 (USUV)</td></tr><tr><td>40</td><td>18</td><td>1</td><td></td><td></td></tr><tr><td>80</td><td>12</td><td></td><td></td><td></td></tr><tr><td>160</td><td>8</td><td></td><td></td><td></td></tr><tr><td>320</td><td>3</td><td></td><td></td><td>1 (USUV)</td></tr><tr><td>≥640</td><td>3</td><td></td><td></td><td>1 (WNV)</td></tr></tbody></table>

<footnote><sup>a</sup>Reference category</footnote>
<figure><img src="image_4.png" /></figure>

WNV and related flaviviruses than other artiodactyl species such as wild boar from the same regions and periods [10, 27]. In North America, *Culex* family mosquitoes (flavivirus vector) have shown a preference of feeding on deer [35]. Similarly, *Culex* species sampled in Canary Islands (Spain) frequently fed on ruminant species [36], while a study in wetlands of bioregion 5 did show feeding on mammals in addition to birds and reptiles but without identifying the mammal hosts beyond humans and horses [37]. Thus, our results would suggest a reduced susceptibility of red deer to flaviviruses of the encephalitis antigenic complex and other cross reacting flaviviruses.

Multivariate logistic regression identified the year 2011, bioregion 5 and the presence of wetland areas as individual risk factors for antigenically-related flavivirus exposure in red deer. The seroprevalence data obtained indicates widespread but not homogeneous distribution of these flaviviruses in wild ruminants in Spain. Consistently with other studies, seroprevalence of flaviviruses in both farmed and free-living wild ruminants, mainly WNV, was higher in southern Spain, namely in bioregions 3, 4 and 5 [10, 38, 39]. In fact, bioregion 5 includes the provinces with the highest number of WNV outbreaks reported in horses and humans and the area where BAGV infection caused high mortality in wild
game birds in 2010 [18, 40]. The results indicate that southern regions present the highest risk of flavivirus circulation in Spain. Even though a study conducted in 1980 described the presence of antibodies against flaviviruses in rodents in northern Spain [41], the extremely low prevalence (0.2 %) of flavivirus antibodies in north Spain (bioregion 1), which has a colder mean temperature [42], suggests lower activity of flaviviruses in this area.

Animals sampled in wetland areas were three times more likely to be exposed to related flaviviruses, which may be due to the larger populations and greater diversity of competent vectors. In addition, in late summer, the number of wetlands and waterholes in southern Spain becomes greatly reduced, aggregating animals in these hotspots [43], which may also act as breeding areas for mosquito species. Our results are consistent with those previously reported in France [44], and indicate that wetland areas may help to target risk-based surveillance in Spain.

Seropositivity of at least one animal in all sample years supports a more-or-less continuous circulation of antigenically-related flaviviruses in the study period. However, a significantly higher seroprevalence was found in 2011, 1 year after the highest number of WNV and Bagaza outbreaks were reported in Spain [18, 40]. In 2011, seroprevalence was detected in adult and sub-adult deer but not in yearling animals suggesting exposure to related flaviviruses of the wild ruminant populations at the end of 2010 and persistence of antibodies in infected animals. Accordingly, a seropositive yearling free-living red deer sampled in 2010 was positive to WNV by VNT, confirming WNV circulation this year. Similar results were obtained in juvenile red deer from a farm in the 2010 outbreak area [27]. Our results are in agreement with previous studies in which an increase in WNV seroprevalence was observed in red deer and white-tailed deer 2 years after the first outbreak occurred in the area [23, 27].

Intrinsic host factors such as age can influence exposure and susceptibility to flavivirus infection. In the present study, seroprevalence increased with age but no significant differences among age classes were observed. The results may be explained by the limited number of samples from yearlings analyzed, due to the fact that yearlings are not usually hunted. A low seroprevalence in yearling red deer from Spain was also previously reported by [27]. Moreover, lower seroprevalences were detected in yearling while-tailed deer but no association with age was observed [22, 23]. Our findings are also consistent with the increase in WNV and related flavivirus seroprevalence with age found in wild boar [10]. Although flavivirus antibody persistence in artiodactyl species is unknown, a longer time span of possible exposure and antibody persistence could be associated with the higher seroprevalence found in adult individuals. Further research with larger sample sizes for the different age classes is needed to have a better

understanding about the association between age and flavivirus infection.

Significantly higher seroprevalence to antigenically-related flavivirus was found in farmed red deer as compared to free-living ones, although this variable was not retained in the multivariate analysis. In contrast, free-living red deer sampled in open and fenced sampling areas did not differ in seroprevalence. Our results indicate that samples from farmed red deer, which are accessible during all year, may be a valuable surveillance tool. Variations in the seroprevalence among seasons were also detected in the bivariate analysis, with significantly higher seroprevalence in summer. Moreover, significantly higher mean seroprevalence was found in autumn in farmed red deer sampled in bioregion 5. Most of the WNV outbreaks in horses in Spain have been reported in late summer and autumn [40], confirming the higher risk of flavivirus infection during these seasons. In addition, the higher seropositivity observed in autumn in farmed deer is also in agreement with what has been reported in wild birds in Spain [5, 6].

To the authors' knowledge, this is the first report on WNV, USUV and MBV antibodies in red deer, which confirms local circulation of these flaviviruses in Spain. Twenty two sera were positive by bELISA but negative to USUV, WNV and MBV in VNT. The results could be due to the circulation of other related flaviviruses [45]. Mortality due to BAGV and LIV infections have recently been detected in wild game birds and ruminant species respectively in Spain [16, 18]. Further studies are needed to determine the possible impact of these flaviviruses on wild ruminant populations. Unfortunately, VNT against BAGV and LIV could not be carried out in the present study. Six red deer presented positive results to WNV and USUV by VNT with titer differences ≤ 2-fold. This could be due to VNT cross reactions, however co-infections cannot be ruled out. In fact, both WNV and USUV as well as WNV and MBV co-circulation were confirmed (in two and one sampling areas, respectively). These results are of interest because a previous infection with a closely related flavivirus might confer cross-reactive immunity in the animal host, reducing the amplification of new strains with higher virulence introduced into the same region [46]. In fact, co-circulation of different flaviviruses, and the induction of cross-protective antibodies, has been suggested as an explanation for the high rate of seropositivity to WNV in birds and the limited number of cases reported in humans [46]. In Spain, both WNV lineage 1 and the putative new lineage, as well as BAGV, USUV and MBV virus have been confirmed [47]. In contrast, the understanding of flavivirus co-infections both in host and vector species is still very limited.

Specific antibodies against WNV were confirmed in four out of the five bioregions. The higher WNV seroprevalence in bioregion 3 and 5 is in accordance with WNV distribution reported in horses in Spain [40]. Interestingly,
WNV antibodies were found in animals in areas where WNV outbreaks have not been reported to date. Our results constitute the first report of WNV circulation in bioregion 4, where habitat conditions do not favour WNV activity. Yearly re-emergence of WNV in Mediterranean countries is thought to be due to WNV overwintering and becoming endemic in local bird populations [48, 49]. The presence of animals seropositive to WNV throughout the sampling period (2003–2014) suggest an endemic circulation of WNV in southern Spain. In fact, after WNV was first reported in horses from southern Spain in 2010, outbreaks in this species have been reported every year [40].

During the last decade, USUV has rapidly expanded in different European countries, and has been detected in mosquitoes and vertebrate species [50]. Even though human infections are believed to be clinically mild or asymptomatic, neuroinvasive illness associated with USUV infection has been reported in Italy and Croatia [51, 52]. In our study, USUV infection was confirmed in red deer from bioregions 1, 3 and 5 in 2007 and 2008. The seroprevalence found in this species, as well as the USUV infections previously reported in mosquitoes, resident and migratory wild birds [4, 14, 15, 53], indicate circulation of this emerging flavivirus during the last decade in Spain. Further epidemiological and molecular research is required to determine the role of wild ruminants in the epidemiology of USUV and to determine USUV strains and lineages circulating in Spain.

Infection with MBV or a Meaban-like virus has been recently detected in yellow-legged and herring (*Larus argentatus*) gulls from different breeding colonies in France and Spain [3], and in a great frigatebird (*Fregata minor*) and a sooty tern (*Onychoprion fuscatus*) in the Western Indian Ocean [54], although the zoonotic potential of this emerging flavivirus is still unknown. Two adult free-living red deer sampled in Seville province (bioregion 3) in 2010 showed antibodies against this flavivirus, which indicates that MBV is also circulating among wild ruminants in Spain. Interestingly, seropositivity against MBV has been detected precisely in the same province and year in a recent study carried out in waterfowl used as decoys in Spain (unpublished data), which suggests limited spread of this flavivirus in the study area. This result is consistent with a previous report of MBV infection in only two specific locations within the Mediterranean Basin [3].

## Conclusions

The results obtained indicate local circulation of antigenically-related flaviviruses in wild ruminants in Spain. The significance of WNV, USUV and MBV infection in this species is unknown at this time and further research on aspects such as the genetic diversity, infection rates and competent vectors of these emerging flaviviruses is needed. Surveillance of hunted wild ruminants, and

especially on farmed red deer, may be a useful, cost-effective, rapid and complementary tool to the national surveillance programs to monitor the activity of mosquito and tick-borne flaviviruses in Europe.

## Abbreviations

BAGV: Bagaza virus; ELISA: Enzyme-linked immunosorbent assay; KFI: Kidney fat index; LIV: Louping-ill virus; MBV: Meaban virus; ORs: Odd ratios; SLEV: Saint Louis encephalitis virus; TBEV: Tick borne encephalitis virus; USUV: Usutu virus; VNT: Virus neutralization test; WNV: West Nile virus

## Acknowledgements

We acknowledge the help of fellow students with sample collection in hunting drives. We are indebted to Dr. Pelayo Acevedo for his help with Fig. 1.

## Funding

This study has been supported by projects PAC08-0296-7771 (JCCM), AG2008-02504GAN, FAU2008-00019-C03-01 and AGL2013-49159-C2-2-R. Gutierrez-Guzman AV, was a JCCM fellow (PAC08-0296-7771).

## Availability of data and materials

The datasets generated during and/or analysed during the current study are not publicly available due to the fact that they are employed further, in yet unpublished studies and due to their origin from private hunting estates and breeding farms that have agreed on participating in the study on the account of confidential handling of the data, but are available from the corresponding author on reasonable request.

## Authors' contributions

Conceived and designed the experiment: IGB, JP, AVGG, MB, AAM, DCT, CG, UH. Performed the experiment: IGB, JP, AVGG, MB, AAM, DCT, CG, UH. Contributed to/carried out analyses: IGB, JP, AVGG, SL, MB, AAM, DCT, SL, CG, UH. Performed statistical analysis: IGB, UH. Drafted and amended the manuscript: IGB, JP, AVGG, MB, AAM, DCT, CG, UH. All authors read and approved the final manuscript.

## Competing interests

The authors declare that they have no competing interests.

## Consent for publication

Not applicable.

## Ethics approval and consent to participate

This study did not involve purposeful killing of animals. Samples from dead animals were collected from individuals legally hunted during the hunting seasons. These animals were legally hunted under Spanish and EU legislation and all hunters had hunting licenses. Blood from farmed live individuals was obtained from archive samples previously collected during routine handling procedures. None of the blood samples was collected specifically for this study, therefore, no ethical approval was deemed necessary. The collection of all samples was performed following routine procedures before the design of this study in compliance with the Ethical Principles in Animal Research. Protocols, amendments and other resources were done according to the guidelines approved by each Autonomous government following the R.D.1201/2005 of the Ministry of Presidency of Spain (10th October 2005, BOE 21st October 2005) (http://www.umh.es/_web_rwlock/ceie/docs/animales/1201_05%20proteccion%20animales%20experimentacion.pdf).

## Author details

¹Departamento de Sanidad Animal, Facultad de Veterinaria, Universidad de Córdoba-Agrifood Excellence International Campus (ceiA3), Rabanales, 14071 Córdoba, Spain. ²Instituto de Investigación en Recursos Cinegéticos IREC, (CSIC-UCLM-JCCM), Ciudad Real, Spain. ³ANSES, Laboratoire de Santé Animale de Maisons-Alfort, UMR 1161 Virologie, INRA, ANSES, ENVA, Maisons-Alfort F-94703, France. ⁴Sabiotec, Camino de Moledores s.n., Ed. Polivalente UCLM, 13005 Ciudad Real, Spain.

Received: 26 June 2016 Accepted: 29 October 2016

Published online: 09 November 2016
## References

1. Weissenböck H, Hubálek Z, Bakonyi T, Nowotny N. Zoonotic mosquito-borne flavivirus: worldwide presence of agents with proven pathogenicity and potencial candidates of future emerging diseases. Vet Microbiol. 2010;140:271-80.

2. Beck C, Jiménez-Clavero MA, Leblond A, Durand B, Nowotny N, Leparc-Goffart I, Zientara S, Jourdain E, Lecollinet S. Flaviviruses in Europe: complex circulation patterns and their consequences for the diagnosis and control of West Nile disease. Int J Environ Res Public Health. 2013;10:6049-83.

3. Arnal A, Gómez-Díaz E, Cerdà-Cuéllar M, Lecollinet S, Pearce-Duvet J, Busquets N, García-Bocanegra I, Pagés N, Vittecoq M, Hammouda A, Samraoui B, Garnier R, Ramos R, Selmi S, González-Solis J, Jourdain E, Boulinier T. Circulation of a meaban-like virus in yellow-legged gulls and seabird ticks in the Western Mediterranean basin. PLoS ONE. 2014;9, e89601.

4. Vázquez A, Jiménez-Clavero MA, Franco L, Donoso-Mantke O, Sambri V, Niedrig M, Zeller H, Tenorio A. Usutu virus - potential risk of human disease in Europe. Euro Surveill. 2011;16:22-6.

5. Figuerola J, Jiménez-Clavero MA, López G, Rubio C, Soriguer R, Gomez-Tejedor C, Tenorio A. Size matters: West Nile virus neutralizing antibodies in resident and migratory birds in Spain. Vet Microbiol. 2008;132:39-46.

6. García-Bocanegra I, Busquets N, Napp S, Alba A, Zorrilla I, Villalba R, Arenas A. Serosurvey of West Nile virus and other flaviviruses of the Japanese encephalitis antigenic complex in birds from Andalusia, Southern Spain. Vector Borne Zoonotic Dis. 2011;11:1107-13.

7. Bernabeu-Wittel M, Ruiz-Pérez MA, del Toro MD, Aznar J, Muniain A, de Ory F. West Nile virus past infections in the general population of Southern Spain. Enferm Infecc Microbiol Clin. 2007;25:561-5.

8. García-Bocanegra I, Jaén-Téllez JA, Napp S, Arenas-Montes A, Fernández-Morente M, Fernández-Molera V, Arenas A. Monitoring of the West Nile Virus epidemic in Spain between 2010 and 2011. Transbound Emerg Dis. 2012;59:448-55.

9. Mentaberre G, Gutiérrez C, Rodríguez NF, Joseph S, González-Barrio D, Cabezón O, de la Fuente J, Gortazar C, Boadella M. A transversal study on antibodies against selected pathogens in dromedary camels in the Canary Islands, Spain. Vet Microbiol. 2013;167:468-73.

10. Gutiérrez-Guzmán A, Vicente J, Sobrino R, Perez-Ramírez E, Llorente F, Höfle U. Antibodies to West Nile virus and related flaviviruses in wild boar, red foxes and other mesomammals from Spain. Vet Microbiol. 2012;159:291-7.

11. Jiménez-Clavero MA, Sotelo E, Fernandez-Pinero J, Llorente F, Blanco JM, Rodríguez-Ramos J, Perez-Ramírez E, Höfle U. West Nile virus in golden eagles, Spain, 2007. Emerg Infect Dis. 2008;14:1489-91.

12. García-Bocanegra I, Jaén-Téllez JA, Napp S, Arenas-Montes A, Fernández-Morente M, Fernández-Molera V, Arenas A. West Nile outbreak in horses and humans, Spain, 2010. Emerg Infect Dis. 2011;17:2397-9.

13. Busquets N, Alba A, Allepuz A, Aranda C, Núñez JE. Usutu virus sequences in *Culex pipiens* (Diptera: *Culicidae*), Spain. Emerg Infect Dis. 2008;14:861-3.

14. Höfle U, Gamino V, Fernández de Mera IG, Mangold AJ, Ortiz JA, de la Fuente J. Usutu virus in migratory song thrushes, Spain. Emerg Infect Dis. 2013;19:1173-5.

15. Cano-Terriza D, Guerra R, Lecollinet S, Cerdà-Cuéllar M, Cabezón O, Almería S, García-Bocanegra I. Epidemiological survey of zoonotic pathogens in feral pigeons (*Columba livia* var. *domestica*) and sympatric zoo species in Southern Spain. Comp Immunol Microbiol Infect Dis. 2015;43:22-7.

16. Balseiro A, Royo LJ, Pérez-Martínez C, Fernández de Mera IG, Höfle U, Polledo L, Marreros N, Casais R, García-Marín JF. Louping III in goats, Spain, 2011. Emerg Infect Dis. 2012;18:976-8.

17. Ruiz-Fons F, Balseiro A, Willoughby K, Oleaga A, Dagleish MP, Pérez-Ramírez E, Havlíková S, Klempa B, Llorente F, Martín-Hernando MP. Clinical infection of Cantabrian chamois (*Rupicapra pyrenaica parva*) by louping ill virus: New concern for mountain ungulate conservation? Eur J Wildlife Res. 2014;60:691-4.

18. García-Bocanegra I, Zorrilla I, Rodríguez E, Rayas E, Camacho L, Redondo I, Gómez-Guillamón F. Monitoring of the Bagaza virus epidemic in wild bird species in Spain, 2010. Transbound Emerg Dis. 2013;60:120-6.

19. Puerto FI, Lorono-Pino MA, Farfan-Ale JA, García-Rejón JE, Cetina-Trejo RC, Hidalgo-Martinez AC, Ramos C, Rosado-Paredes EP, Flores-Flores LF. Captive animals as sentinels for West Nile virus transmission in zoos from Yucatan and Tabasco states of Mexico. Am J Trop Med Hyg. 2007;77:134.

20. Nitatpattana N, Le Flohic G, Thongchai P, Nakgoi K, Palaboodeewat S, Khin M, Barbazan P, Yoksan S, González JP. Elevated Japanese encephalitis virus activity monitored by domestic sentinel piglets in Thailand. Vector Borne Zoonotic Dis. 2011;11:391-4.

21. Phoutrides E, Josino-Mendez T, Pérez-Medina T, Seda-Lozada R, Garcia-Negron M, Davila-Toro F, Hunsperger E. The utility of animal surveillance in the detection of West Nile Virus activity in Puerto Rico. Vector Borne Zoonotic Dis. 2011;11:447-50.

22. Farajollahi A, Gates R, Crans W, Komar N. Serologic evidence of West Nile virus and St. Louis encephalitis virus infections in while-tailed deer (*Odocoileus virginianus*) from New Jersey, 2001. Vector Borne Zoonotic Dis. 2004;4:379-83.

23. Santaella J, Mclean R, Hall JS, Gill JS, Bowen RA, Hadow HH, Clark L. West Nile virus serosurvailance in Iowa White-tailde deer (199-2003). Am J Trop Med Hyg. 2005;73:1038-42.

24. Juricova Z, Hubalek Z. Serological surveys for arboviruses in the game animals of southern Moravia (Czech Republic). Folia Zool. 1999;48:185-9.

25. Palmer MV, Stoffregen WC, Rogers DG, Hamir AN, Richt JA, Pedersen DD. West Nile virus in reindeer (*Rangifer tarandus*). J Vet Diagn Invest. 2004;16:219-22.

26. Miller DL, Radi ZA, Baldwin C, Ingram D. Fatal West Nile virus infection in a white-tailed deer (*Odocoileus virginianus*). J Wildlife Dis. 2005;41:246-9.

27. Boadella M, Díez-Delgado I, Gutiérrez-Guzmán AV, Höfle U, Gortázar C. Do wild ungulates allow improved monitoring of flavivirus circulation in Spain? Vector Borne Zoonotic Dis. 2012;12:490-5.

28. Acevedo P, Ruiz-Fons F, Vicente J, Reyes-García AR, Alzaga V, Gortázar C. Estimating red deer abundance in a wide range of management situations in Mediterranean habitats. J Zool. 2008;276:37-47.

29. Muñoz P, Boadella M, Arnal M, de Miguel M, Revilla M, Martinez D, Vicente J, Acevedo P, Oleaga A, Ruiz-Fons F, Marín CM, Prieto JM, de la Fuente J, Barral M, Barberán M, Fernández de Luco D, Blasco JM, Gortázar C. Spatial distribution and risk factors of Brucellosis in Iberian wild ungulates. BMC Infect Dis. 2010;10:46.

30. MMARM, (Ministerio de Medio Ambiente y Medio Rural y Marino de España), Plan Nacional de Vigilancia Sanitaria en Fauna Silvestre. 2016. http://rasve.magrama.es/publica/programas/NORMATIVA%20Y%20PROGRAMAS%5CPROGRAMAS%5CFAUNA%20SILVESTRE%5CPLAN%20NACIONAL%20DE%20VIGILANCIA%20EN%20FAUNA%20SILVESTRE.PDF. Accessed 15 Apr 2016.

31. Sáenz de Buruaga M, Lucio-Calero A, Purroy-Iraizoz FJ. Reconocimiento de Sexo y Edad en Especies Cinegéticas. 1st ed. Spain: Edilesa; 2001.

32. Santos JPV, Fernández-de-Mera IG, Acevedo P, Boadella M, Fierro Y, Vicente J, Gortázar C. Optimizing the sampling effort to evaluate body condition in ungulates: a case study on red deer. Ecol Ind. 2013;30:65-71.

33. Chaintoutis SC, Dovas CI, Papanastassopoulou M, Gewehr S, Danis K, Beck C. Evaluation of a West Nile virus surveillance and early warning system in Greece, based on domestic pigeons. Comp Immunol Microbiol Infect Dis. 2014;37:131-41.

34. Hosmer DW, Lemeshow S. Applied Logistic Regression. 2nd ed. New York: Wiley; 2000.

35. Watts SL, Fitzpatrick DM, Maruniak JE. Blood meal identification from Florida mosquitoes (Diptera: *Culicidae*). Florida Entomol. 2009;92:619-22.

36. de la Martínez- Puente J, Martínez J, Ferraguti M, de la Morales- Nuez A, Castro N, Figueroa J. Genetic characterization and molecular identification of the bloodmeal sources of the potential bluetongue vector *Culicoides obsoletus* in the Canary islands, Spain. Parasit Vector. 2012;5:147.

37. Muñoz J, Ruiz S, Soriguer R, Alcaide M, Viana DS, David R, Vázquez A, Figueroa J. Feeding patterns of potential West Nile virus vectors in South-West Spain. PLoS ONE. 2012;7, e39549.

38. Figueroa J, Jiménez-Clavero MA, Rojo G, Gómez C, Soriguer R. Prevalence of West Nile virus neutralizing antibodies in colonial aquatic birds in southern Spain. Avian Pathol. 2007;36:209-12.

39. López G, Jiménez-Clavero MA, Vázquez A, Soriguer R, Gómez-Tejedor C, Tenorio A, Figueroa J. Incidence of West Nile virus in birds arriving in wildlife rehabilitation centers in Southern Spain. Vector Borne Zoonotic Dis. 2010;11:285-90.

40. RASVE (Red de Alerta Sanitaria Veterinaria). 2016. http://rasve.magrama.es/RASVE_2008/Publica/Focos/Consulta.aspx. Accessed 15 May 2016.

41. Chastel C, Launay H, Rogues G, Beaucornu JC. Infections à arbovirus en Espagne: enquête sérologique chez les petits mammifères. Bull Soc Path Exot. 1980;73:384-90.

42. Ninyerola M, Pons X, Roure JM. Atlas climático digital de la Península Ibérica. Metodología y aplicaciones en Bioclimatología y geobotánica. 2005. http://opengis.uab.es/wms/iberia/pdf/acdpi.pdf. Accessed 3 May 2016.

43. Barasona JA, Latham MC, Acevedo P, Armenteros JA, Latham ADM, Gortázar C, Carro F, Soriguer RC, Vicente J. Spatiotemporal interactions between wild boar and cattle: implications for cross-species disease transmission. Vet Res. 2014;45:122.
44. Pradier S, Sandoz A, Paul MC, Lefebvre G, Tran A, Maingault J, Lecollinet S, Leblond A. Importance of wetlands management for West Nile Virus circulation risk, Camargue, Southern France. Int J Environ Res Public Health. 2014;11:7740–54.

45. Blitvich BJ, Fernandez-Salas I, Contreras-Cordero JF, Marlenee NL, González-Rojas JI, Komar N. Serologic evidence of West Nile virus infection in horses, Coahuila State, Mexico. Emerg Infect Dis. 2003;9:853–6.

46. Rizzoli A, Jiménez-Clavero MA, Barzon L, Cordioli P, Figuerola J, Koraka P, Martina B, Moreno A, Nowotny N, Pardigon N, Sanders N, Ulbert S, Tenorio A. The challenge of West Nile virus in Europe: knowledge gaps and research priorities. Euro Surveill. 2015;20(20). Available online: http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId=21135.

47. Vázquez A, Sanchez-Seco MP, Ruiz S, Molero F, Hernandez L, Moreno J, Magallanes A, Tejedor CG, Tenorio A. Putative new lineage of West Nile virus, Spain. Emerg Infect Dis. 2010;16:549–52.

48. Monaco F, Savini G, Calistri P, Polci A, Pinoni C, Bruno R, Lelli R. 2009 West Nile disease epidemic in Italy: first evidence of overwintering in Western Europe? Res Vet Sci. 2011;91:321–6.

49. Di Sabatino D, Bruno R, Sauro F, Danzetta ML, Cito F, Iannetti S, Narcisi V, De Massis F, Calistri P. Epidemiology of West Nile disease in Europe and in the Mediterranean Basin from 2009 to 2013. Biomed Res Int. 2014;2014:907852.

50. Ashraf U, Ye J, Ruan X, Wan S, Zhu B, Cao S. Usutu virus: an emerging flavivirus in Europe. Viruses. 2015;7:219–38.

51. Pecorari M, Longo G, Gennari W, Grottola A, Sabbatini A, Tagliazucchi S, Savini G, Monaco F, Simone M, Lelli R, Rumpianesi F. First human case of Usutu virus neuroinvasive infection, Italy, August-September 2009. Euro Surveill. 2009;14. pii: 19446.

52. Santini M, Vilibic-Cavlek T, Barsic B, Barbic L, Savic V, Stevanovic V, Listes E, Di Gennaro A, Savini G. First cases of human Usutu virus neuroinvasive infection in Croatia, August-September 2013: clinical and laboratory features. J Neurovirol. 2015;21:92–7.

53. Llorente F, Pérez-Ramírez E, Fernández-Pinero J, Soriguer R, Figuerola J, Jiménez-Clavero MA. Flaviviruses in Game Birds, Southern Spain, 2011–2012. Emerg Infect Dis. 2013;19:1023–5.

54. Jaeger A, Lecollinet S, Beck C, Bastien M, Corre MLE, Dellagi K, Pascalis H, Boulinier T, Lebarbenchon C. Serological evidence for the circulation of flavivirus in seabird populations of the western Indian Ocean. Epidemiol Infect. 2015;21:1–9.

Submit your next manuscript to BioMed Central
and we will help you at every step:

* We accept pre-submission inquiries

* Our selector tool helps you to find the most relevant journal

* We provide round the clock customer support

* Convenient online submission

* Thorough peer review

* Inclusion in PubMed and all major indexing services

* Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central