West Nile virus lineage 2 (WNV-2) was mainly confined to sub-Saharan Africa until the early 2000s, when it was identified for the first time in Central Europe causing outbreaks of human and animal infection. The aim of this study was to reconstruct the origin and dispersion of WNV-2 in Central Europe and Italy on a phylodynamic and phylogeographical basis. To this aim, discrete and continuous space phylogeographical models were applied to a total of 33 newly characterised full-length viral genomes obtained from mosquitoes, birds and humans in Northern Italy in the years 2013–2015 aligned with 64 complete sequences isolated mainly in Europe. The European isolates segregated into two highly significant clades: a small one including three sequences and a large clade including the majority of isolates obtained in Central Europe since 2004. Discrete phylogeographical analysis showed that the most probable location of the root of the largest European clade was in Hungary a mean 12.78 years ago. The European clade bifurcated into two highly supported subclades: one including most of the Central/East European isolates and the other encompassing all of the isolates obtained in Greece. The continuous space phylogeographical analysis of the Italian clade showed that WNV-2 entered Italy in about 2008, probably by crossing the Adriatic sea and reaching a central area of the Po Valley. The epidemic then spread simultaneously eastward, to reach the region of the Po delta in 2013, and westward to the border area between Lombardy and Piedmont in 2014; later, the western strain changed direction southward, and reached the central area of the Po valley once again in 2015. Over a period of about seven years, the virus spread all over an area of northern Italy by following the Po river and its main tributaries.

No abstract

Lineage 2 West Nile virus (WNV) caused a vast epidemic in Europe in 2018, with the highest incidence being recorded in Italy. To reconstruct the evolutionary dynamics and epidemiological history of the virus in Italy, 53 envelope gene and 26 complete genome sequences obtained from human and animal samples were characterised by means of next-generation sequencing. Phylogenetic analysis revealed two Italian strains originating between 2010 and 2012: clade A, which apparently became extinct in 2013–2014, and clade B, which was responsible for the 2018 epidemic. The mean genetic distances in clade B increased over time and with the distance between sampling locations. Bayesian birth-death and coalescent skyline plots of the clade B showed that the effective number of infections and the effective reproduction number (Re) increased between 2015 and 2018. Our data suggest that WNV-2 entered Italy in 2011 as a result of one or a few penetration events. Clade B differentiated mainly as a result of genetic drift and purifying selection, leading to the appearance of multiple locally circulating sub-clades for different times. Phylodynamic analysis showed a current expansion of the infection among reservoir birds and/or vectors.

Border disease virus (BDV) affects a wide range of ruminants worldwide, mainly domestic sheep and goat. Since 2001 several outbreaks of disease associated to BDV infection have been described in Pyrenean chamois (Rupicapra pyrenaica pyrenaica) in Spain, France and Andorra. In order to reconstruct the most probable places of origin and pathways of dispersion of BDV among Pyrenean chamois, a phylogenetic analysis of 95 BDV 5’untranslated sequences has been performed on chamois and domestic ungulates, including novel sequences and retrieved from public databases, using a Bayesian Markov Chain Monte Carlo method. Discrete and continuous space phylogeography have been applied on chamois sequences dataset, using centroid positions and latitude and longitude coordinates of the animals, respectively.The estimated mean evolutionary rate of BDV sequences was 2.9×10−3 subs/site/year (95% HPD: 1.5–4.6×10−3). All the Pyrenean chamois isolates clustered in a unique highly significant clade, that originated from BDV-4a ovine clade. The introduction from sheep (dated back to the early 90s) generated a founder effect on the chamois population and the most probable place of origin of Pyrenean chamois BDV was estimated at coordinates 42.42 N and 1.9 E. The pathways of virus dispersion showed two main routes: the first started on the early 90s of the past century with a westward direction and the second arise in Central Pyrenees. The virus spread westward for more than 125 km and southward for about 50km and the estimated epidemic diffusion rate was about 13.1 km/year (95% HPD 5.2–21.4 km/year). The strong spatial structure, with strains from a single locality segregating together in homogeneous groups, and the significant pathways of viral dispersion among the areas, allowed to reconstruct both events of infection in a single area and of migrations, occurring between neighboring areas.