Comparisons of two co-circulating clusters of WNV-2a in Europe.

(a) The mean time of the most recent ancestor (TMRCA) and 95%HPD interval for each cluster. (b) The mean clock rate (substitutions per site per year, subst./site/yr) and 95% HPD interval for each cluster was estimated using an uncorrelated relaxed molecular clock model. (c) The Mean number of Markov jump between countries and 95%HPD interval for each cluster were estimated using a continuous-time Markov chain (CTMC) model. (d) Estimation of effective population size via time and a 95% HPD interval using the Skygrid coalescent model. The logarithmic effective number of infections (Ne) vs. viral generation time (t), representing effective transmissions is plotted over time. (e) The mean of weighted dispersal velocity averaged over the branches of the entire tree (km/yr) and (f) The weighted dispersal velocity over time (km/yr) with a 95% HPD interval estimated using the continuous phylogenetic diffusion model.</p

Results of drivers on viral dispersal velocity*.

BackgroundWest Nile virus (WNV) outbreaks in birds, humans, and livestock have occurred in multiple areas in Europe and have had a significant impact on animal and human health. The patterns of emergence and spread of WNV in Europe are very different from those in the US and understanding these are important for guiding preparedness activities.MethodsWe mapped the evolution and spread history of WNV in Europe by incorporating viral genome sequences and epidemiological data into phylodynamic models. Spatially explicit phylogeographic models were developed to explore the possible contribution of different drivers to viral dispersal direction and velocity. A “skygrid-GLM” approach was used to identify how changes in environments would predict viral genetic diversity variations over time.FindingsAmong the six lineages found in Europe, WNV-2a (a sub-lineage of WNV-2) has been predominant (accounting for 73% of all sequences obtained in Europe that have been shared in the public domain) and has spread to at least 14 countries. In the past two decades, WNV-2a has evolved into two major co-circulating clusters, both originating from Central Europe, but with distinct dynamic history and transmission patterns. WNV-2a spreads at a high dispersal velocity (88km/yr–215 km/yr) which is correlated to bird movements. Notably, amongst multiple drivers that could affect the spread of WNV, factors related to land use were found to strongly influence the spread of WNV. Specifically, the intensity of agricultural activities (defined by factors related to crops and livestock production, such as coverage of cropland, pasture, cultivated and managed vegetation, livestock density) were positively associated with both spread direction and velocity. In addition, WNV spread direction was associated with high coverage of wetlands and migratory bird flyways.ConclusionOur results suggest that—in addition to ecological conditions favouring bird- and mosquito- presence—agricultural land use may be a significant driver of WNV emergence and spread. Our study also identified significant gaps in data and the need to strengthen virological surveillance in countries of Central Europe from where WNV outbreaks are likely seeded. Enhanced monitoring for early detection of further dispersal could be targeted to areas with high agricultural activities and habitats of migratory birds.</div

Time-scaled MCC trees of WNV mapping with countries.

Trees were reconstructed using sequences of ns3 gene (a) and ns5 gene (b) mapping with countries (n = 14). (DOCX)</p

Distribution of predictors for viral genetic diversity over time.

Data in four months (January, April, July, and October) were shown for each predictor. Definition for each predictor is shown in the (S3 Table). The unit of each predictor is shown after the predictor name above each panel. The European shapefile was created using the R package “rworldmap” (https://cran.r-project.org/web/packages/rworldmap/). (DOCX)</p

Geographic distributions of WNV-2a sequences dataset.

a) ns3 gene, b) ns5 gene, c) full genomes and d) the number of sequences for each dataset. The European shapefile was created using the R package “maps” (https://cran.r-project.org/web/packages/maps/). (DOCX)</p

Predictors of West Nile Virus Dispersal in Europe Assessed for Quality.

Predictors of West Nile Virus Dispersal in Europe Assessed for Quality.</p

Sequences and metadata of WNV lineage 2a used in this study.

Sequences and metadata of WNV lineage 2a used in this study.</p

Quantified transmission network of WNV-2a between European countries and within Greece inferred using discrete trait models.

The shape of colors on the map indicates the number of samples; the edge weight indicates the median number of transmissions between pairs of countries/regions; the arrow on the edge indicates transmission direction; color of the edge from light to dark indicates Bayes Factor (BF) support from low to high only transmissions with BF>3 are shown). The European shapefile was created using the R package “rnaturalearth” (https://cran.r-project.org/web/packages/rnaturalearth/).</p

Predictors for viral genetic diversity over time.

BackgroundWest Nile virus (WNV) outbreaks in birds, humans, and livestock have occurred in multiple areas in Europe and have had a significant impact on animal and human health. The patterns of emergence and spread of WNV in Europe are very different from those in the US and understanding these are important for guiding preparedness activities.MethodsWe mapped the evolution and spread history of WNV in Europe by incorporating viral genome sequences and epidemiological data into phylodynamic models. Spatially explicit phylogeographic models were developed to explore the possible contribution of different drivers to viral dispersal direction and velocity. A “skygrid-GLM” approach was used to identify how changes in environments would predict viral genetic diversity variations over time.FindingsAmong the six lineages found in Europe, WNV-2a (a sub-lineage of WNV-2) has been predominant (accounting for 73% of all sequences obtained in Europe that have been shared in the public domain) and has spread to at least 14 countries. In the past two decades, WNV-2a has evolved into two major co-circulating clusters, both originating from Central Europe, but with distinct dynamic history and transmission patterns. WNV-2a spreads at a high dispersal velocity (88km/yr–215 km/yr) which is correlated to bird movements. Notably, amongst multiple drivers that could affect the spread of WNV, factors related to land use were found to strongly influence the spread of WNV. Specifically, the intensity of agricultural activities (defined by factors related to crops and livestock production, such as coverage of cropland, pasture, cultivated and managed vegetation, livestock density) were positively associated with both spread direction and velocity. In addition, WNV spread direction was associated with high coverage of wetlands and migratory bird flyways.ConclusionOur results suggest that—in addition to ecological conditions favouring bird- and mosquito- presence—agricultural land use may be a significant driver of WNV emergence and spread. Our study also identified significant gaps in data and the need to strengthen virological surveillance in countries of Central Europe from where WNV outbreaks are likely seeded. Enhanced monitoring for early detection of further dispersal could be targeted to areas with high agricultural activities and habitats of migratory birds.</div

Time-scaled phylogeny of WNV-2a genomes in Europe.

(a) Time-scaled MCC (maximum clade credibility) tree of WNV full genome sequences isolated in Europe (n = 192), the two clusters A and B are labelled on the right. A distinct phylogeographic analysis has been based on the NS3 gene by continuous phylogeographic inference based on 1,000 posterior trees. Spatiotemporal diffusion of all WNV- 2a in Europe (b), of Cluster A (c) and Cluster B (d). These MCC trees are superimposed on 80% of the highest posterior density (HPD) interval reflecting phylogeographic uncertainty. Nodes of the trees, as well as HPD regions, are colored by timescale from red (the time to the most recent common ancestor, TMRCA) to green (most recent sampling time), and the oldest nodes (and corresponding HPD regions) are here plotted on top of youngest nodes. The European shapefile was created using the R package “raster” (https://cran.r-project.org/web/packages/raster/).</p