West Nile virus spread in Europe: Phylogeographic pattern analysis and key drivers.

BACKGROUND: West 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. METHODS: We 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. FINDINGS: Among 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. CONCLUSION: Our 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.

ICTV Virus Taxonomy Profile: Nairoviridae 2024.

Nairoviridae is a family for negative-sense RNA viruses with genomes of about 17.2-21.1 kb. These viruses are maintained in and/or transmitted by arthropods among birds, reptiles and mammals. Norwaviruses and orthonairoviruses can cause febrile illness in humans. Several orthonairoviruses can infect mammals, causing mild, severe and sometimes, fatal diseases. Nairovirids produce enveloped virions containing two or three single-stranded RNA segments with open reading frames that encode a nucleoprotein (N), sometimes a glycoprotein precursor (GPC), and a large (L) protein containing an RNA-directed RNA polymerase (RdRP) domain. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) report on the family Nairoviridae, which is available at www.ictv.global/report/nairoviridae.

History and classification of Aigai virus (formerly Crimean-Congo haemorrhagic fever virus genotype VI).

Crimean-Congo haemorrhagic fever virus (CCHFV) is the medically most important member of the rapidly expanding bunyaviral family Nairoviridae. Traditionally, CCHFV isolates have been assigned to six distinct genotypes. Here, the International Committee on Taxonomy of Viruses (ICTV) Nairoviridae Study Group outlines the reasons for the recent decision to re-classify genogroup VI (aka Europe-2 or AP-92-like) as a distinct virus, Aigai virus (AIGV).

Early initiation of the respiratory syncytial virus season in 2021-2022, Greece.

Respiratory syncytial virus (RSV) is the most common viral pathogen causing respiratory disease in the pediatric population. An unexpected sudden upsurge of RSV infections among children was observed in September 2021 in Greece. Forty-one rhinopharyngeal samples from children under the age of 2 years with confirmed RSV bronchiolitis were tested to identify the genotype(s) of the RSV strain(s). The children were hospitalized during September-November 2021 in three tertiary hospitals in northern Greece. A one-step RT-PCR which amplifies a fragment of the second hypervariable region of the G protein gene was applied. PCR products were sequenced, and phylogenetic analysis was performed. Most (80.5%) RSV cases were typed as RSV-A, with RSV-B accounting for 19.5% of cases. RSV-A and RSV-B sequences clustered within the ON1 and BA genotypes, respectively. As the same genotypes were detected in cases observed during 2016-2018 in northern Greece, it was suggested that the early upsurge of infections was not related to the emergence of novel strain(s), but it was the result of the absence of immunity among children and their mothers due to the restriction measures taken during the COVID-19 pandemic in the previous RSV season. Awareness is needed to diagnose even the out-of-season RSV infections, while molecular epidemiology plays a key role in monitoring the efficacy of currently available therapeutics and for those under development.

Emergence of ST39 carbapenem-resistant Klebsiella pneumoniae producing VIM-1 and KPC-2.

INTRODUCTION: Carbapenem-resistant Klebsiella pneumoniae (CRKP) causes life-threatening hospital-acquired infections. KPC and VIM carbapenemase production is the main molecular mechanism for carbapenem resistance. The aim of the current study was the genetic characterization of four ST39 CRKP isolates simultaneously producing VIM-1 and KPC-2, obtained in a Greek tertiary hospital. METHODS: Identification and antimicrobial susceptibility testing were performed through VITEK 2. Multiplex PCR, multiplex lateral flow immunoassay, phenotypic tests and next generation sequencing were applied. The sequence reads were de novo assembled and annotated, while antimicrobial resistance genes and plasmids were identified using bioinformatics software. Genomic comparison and core genome single-nucleotide polymorphism-based phylogenetic analysis were also performed. RESULTS: Three isolates were pandrug-resistant, and one was extensively drug-resistant; they all carried blaVIM-1 and blaKPC-2 genes and were assigned to ST39. BlaVIM-1 was integrated in a class 1 integron. They all harboured many antimicrobial resistance genes and various plasmids. The mgrB gene of all isolates was disrupted by an insertion sequence (ISKpn14). Genome comparison and phylogenetic analysis revealed that the isolates were closely related. CONCLUSION: To our knowledge this is the first report on detection of CRKP ST39 isolates simultaneously producing VIM-1 and KPC-2 in addition to colistin resistance. The knowledge of the clonal relatedness of the isolates can lead to the implementation of strict infection control measures absolutely needed to eliminate their spread.

First report of canine Astrovirus and Sapovirus in Greece, hosting both asymptomatic and gastroenteritis symptomatic dogs.

BACKGROUND: Astrovirus, Norovirus and Sapovirus are widely distributed viruses in humans and animals worldwide. They have frequently been associated with disease, mainly of gastroenteric nature. In dogs, these viruses have been detected both in symptomatic and asymptomatic animals, mainly of young age. METHODS: In the present epidemiologic study, we investigated the presence of canine Astrovirus (CAstV), canine Norovirus (canine NoV) and canine Sapovirus (Canine SaV) in saliva and stools of 201 domestic dogs originating from throughout Greece, based on two different molecular methods, i.e. conventional and SYBR-Green Real-time RT-PCR. The samples derived from young and adult asymptomatic and symptomatic animals. CAstV was detected in 15/201 (7.5%) and 29/201 (15%) of the examined dogs using conventional RT-PCR and SYBR-Green Real time RT-PCR, respectively. RESULTS: The prevalence of the virus was higher at healthy dogs, with a slight discrepancy of the two methods on the aspect of age (67% young dogs with the method of conventional RT-PCR, versus 52% adult positive dogs with the method of SYBR-Green Real-time RT-PCR). Canine SaV was detected in 52/201 (23%) of the dogs (mainly young and asymptomatic), with the method of SYBR-Green Real-time RT-PCR only, while canine NoV was not detected in any sample with either of the two methods applied. Sequencing of the CAstV positive samples resulted in the acquisition of one CAstV sequence. Phylogenetic analysis confirmed the results, clustering the CAstV sequence with homologous canine hosting sequences from other countries. CONCLUSIONS: CAstV and Canine SaV were proved to circulate in Greek dogs. SYBR-Green Real time RT-PCR showed greater sensitivity in the detection of these viruses. Additionally, we were able to specify the CAstV strain that circulates in Greece, through phylogenetic analysis. To our knowledge, this is the first epidemiological study of CAstV and canine SaV in dogs in Greece, as well as the first time detected in dogs from Greece.

Changes in molecular epidemiology of carbapenem-resistant Klebsiella pneumoniae in the intensive care units of a Greek hospital, 2018-2021.

The spread of multi-drug resistant (MDR) Gram-negative bacteria, including Klebsiella pneumoniae, constitutes a global threat. The most frequent mechanism of acquired carbapenem resistance is the production of carbapenemases, especially KPC, NDM, VIM, IMP and OXA-48. We analyzed the epidemiological trend of carbapenem resistance genes of carbapenem-resistant K. pneumoniae (CRKP) strains isolated from critically ill patients in a Greek tertiary hospital. The study included 150 CRKP isolates collected from 116 (77.4%) patients hospitalized in the adult ICU and 17 (11.3%) each in the pediatric and the two neonatal ICUs between March 2018 and March 2021. Identification and antimicrobial susceptibility testing were performed using VITEK-2. A multiplex lateral flow immunoassay was used for the detection of carbapenemases, while the detection of bla VIM, bla KPC, bla NDM, bla IMP and bla OXA-48-like genes was achieved by multiplex PCR. The bla NDM was mainly detected in adults (54/116, 46.9%), while in children the most often detected gene was bla KPC (24/34, 70.6%). The predominant carbapenem resistance gene during 2018-2019 was bla KPC alone or in combination with bla VIM, reaching 44.4% in 2019, while during 2020-2021 the detection of bla NDM prevailed significantly, reaching 45.5 and 60.7% for 2020 and 2021, respectively. A shift in the molecular epidemiology of CRKP was seen during 2018-2021, which is probably associated with the recent excessive empiric use of newer antimicrobials. Surveillance studies and proper and strict implementation of infection control measures are highly needed to decrease the spread of MDR bacteria, including CRKP.

Spa diversity and genetic characterization of t127 methicillin-resistant Staphylococcus aureus in a tertiary Greek hospital.

Introduction: Methicillin-resistant Staphylococcus aureus (MRSA) causes severe community and hospital acquired infections. Identification of staphylococcal cassette chromosome mec (SCCmec), multilocus-sequence typing, and sequencing of S. aureus protein A (spa) gene are used for MRSA typing. The aim was to investigate the spa types of MRSA isolates in a tertiary hospital in Greece and analyse the whole genome sequences of two t127 MRSA isolates. Methods: Totally, 39 MRSA isolates collected from July 2019 to June 2020 in "Georgios Gennimatas" General Hospital of Thessaloniki, Greece, were included in the study. Identification and antimicrobial susceptibility testing were performed using VITEK II automated system, and spa typing was performed. A minimum spanning tree was used to display the spa type frequencies and the genetic distances among them. Two t127-MRSA isolates (IM-MRSA and PD-MRSA) were selected for WGS. Results: Six isolates (15.4%) were resistant to mupirocin, 18 (46.2%) to fusidic acid, three (7.7%) to vancomycin and two (5.1%) to teicoplanin. Twenty-two different spa types were detected, with t002, t003, and t422 being the most frequent (5/39, 12.8% each), followed by t1994 (4/39, 10.3%). The isolates presented high genetic diversity and, taking into account the time between hospital admission and sampling, intrahospital spread did not occur. Even the two t127 isolates were assigned to different sequence types, ST9-XII-t127 and ST1-IVa-t127. Plasmids and genes conferring antimicrobial resistance and virulence were also identified. Conclusions: Various spa types were identified and together with the information about the time between hospital admission and sampling supports polyclonal MRSA spread in the hospital excluding a nosocomial infection. WGS provides a more detailed analysis distinguishing even the isolates belonging to the same spa type.

Association Between Upper Respiratory Tract Viral Load, Comorbidities, Disease Severity, and Outcome of Patients With SARS-CoV-2 Infection.

BACKGROUND: There is limited information on the association between upper respiratory tract (URT) viral loads, host factors, and disease severity in SARS-CoV-2-infected patients. METHODS: We studied 1122 patients (mean age, 46 years) diagnosed by polymerase chain reaction (PCR). URT viral load, measured by PCR cycle threshold, was categorized as high, moderate, or low. RESULTS: There were 336 (29.9%) patients with comorbidities; 309 patients (27.5%) had high, 316 (28.2%) moderate, and 497 (44.3%) low viral load. In univariate analyses, compared to patients with moderate or low viral load, patients with high viral load were older, more often had comorbidities, developed Symptomatic disease (COVID-19), were intubated, and died. Patients with high viral load had longer stay in intensive care unit and longer intubation compared to patients with low viral load (P values < .05 for all comparisons). Patients with chronic cardiovascular disease, hypertension, chronic pulmonary disease, immunosuppression, obesity, and chronic neurological disease more often had high viral load (P value < .05 for all comparisons). In multivariate analysis high viral load was associated with COVID-19. Level of viral load was not associated with any other outcome. CONCLUSIONS: URT viral load could be used to identify patients at higher risk for morbidity or severe outcome.

Transmission dynamics of SARS-CoV-2 within families with children in Greece: A study of 23 clusters.

There is limited information on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection clustering within families with children. We aimed to study the transmission dynamics of SARS-CoV-2 within families with children in Greece. We studied 23 family clusters of coronavirus disease 2019 (COVID-19). Infection was diagnosed by reverse-transcriptase polymerase chain reaction in respiratory specimens. The level of viral load was categorized as high, moderate, or low based on the cycle threshold values. There were 109 household members (66 adults and 43 children). The median attack rate per cluster was 60% (range: 33.4%-100%). An adult member with COVID-19 was the first case in 21 (91.3%) clusters. Transmission of infection occurred from an adult to a child in 19 clusters and/or from an adult to another adult in 12 clusters. There was no evidence of child-to-adult or child-to-child transmission. In total 68 household members (62.4%) tested positive. Children were more likely to have an asymptomatic SARS-CoV-2 infection compared to adults (40% vs 10.5%; P = .021). In contrast, adults were more likely to develop a severe clinical course compared with children (8.8% vs 0%; P = .021). In addition, infected children were significantly more likely to have a low viral load while adults were more likely to have a moderate viral load (40.7% and 18.6% vs 13.8% and 51.7%, respectively; P = .016). In conclusion, while children become infected by SARS-CoV-2, they do not appear to transmit infection to others. Furthermore, children more frequently have an asymptomatic or mild course compared to adults. Further studies are needed to elucidate the role of viral load on these findings.

PCR-based next-generation West Nile virus sequencing protocols.

The epidemiology of West Nile virus (WNV) is unpredictable and changing. Availability of whole genome sequences enables the detailed molecular epidemiology studies and the evaluation and design of diagnostic tools. In the present study we provide two PCR-based protocols which can be applied directly on biological samples from hosts infected by WNV strains belonging to lineage 1 or lineage 2. It was shown that the protocols worked successfully even on samples with relatively low viral load.

Genetic characterization of two methicillin-resistant Staphylococcus aureus spa type t127 strains isolated from workers in the dairy production chain in Greece.

Methicillin-resistant Staphylococcus aureus (MRSA) constitutes a constant threat for the public health. Aim of the present study was to analyse the whole genome sequences of two MRSA strains belonging to Staphylococcus protein A (spa) type t127 isolated from humans working in two distantly located dairy production farms in Greece.MRSA strains were isolated from the nasal cavity of a food handler in a milk industry in Epirus, northwestern Greece (E-MRSA), and a person working in a cattle farm in Thrace, northeastern Greece (T-MRSA). Whole genome sequences taken using next generation sequencing were analysed for resistance and virulence genes applying various bioinformatic tools.Both isolates were assigned to ST1-IVa-t127 type, and they were transferring genes conferring resistance to tetracycline, ß-lactams, and aminoglycosides; T-MRSA was carrying additional genes leading to macrolide, lincosamide and streptogramin B (MLSB) resistance. Both isolates were carrying three plasmid replicon types, rep5, rep7 and rep16, while T-MRSA harboured also rep10 and rep15. E-MRSA carried scn and sak genes which were absent from T-MRSA.In conclusion, the genetic characterization of two unrelated ST1-IVa-t127 MRSA strains isolated from humans in close contact with livestock in Greece can be used as basis for further epidemiological and evolutionary studies.

Spread of NDM-producing Klebsiella pneumoniae in a tertiary Greek hospital.

Bacterial carbapenem resistance, especially when mediated by transferable carbapenemases, is of important public health concern. An increased number of metallo-ß-lactamase (MBL)-producing Klebsiella pneumoniae strains isolated in a tertiary hospital in Thessaloniki, Greece, called for further genetic investigation.The study included 29 non-repetitive carbapenem resistant K. pneumoniae isolates phenotypically characterized as MBL-producers collected in a tertiary hospital in Greece. The isolates were screened for the detection of carbapenemase genes (K. pneumoniae carbapenemase (blaKPC), Verona-integron-encoded MBL-1 (blaVIM-1), imipenemase (blaIMP), oxacillinase-48 (blaOXA-48) and New Delhi MBL (blaNDM)). The genetic relationship of the isolates was determined by Random Amplified Polymorphic DNA (RAPD) analysis. The whole genome sequences (WGS) from two NDM-positive K. pneumoniae isolates were further characterized.The presence of New Delhi MBL (blaNDM) gene was confirmed in all K. pneumoniae isolates, while blaKPC and blaVIM-1 genes were co-detected in one and two isolates, respectively. The RAPD analysis showed that the isolates were clustered into two groups. The whole genome sequence analysis of two K. pneumoniae isolates revealed that they belonged to the sequence type 11, they carried the blaNDM-1 gene, and exhibited differences in the number and type of the plasmids and the resistant genes.All MBL-producing K. pneumoniae isolates of the study harbored a blaNDM gene, while WGS analysis revealed genetic diversity in resistance genes. Continuous surveillance is needed to detect the emergence of new clones in a hospital setting, while application of antimicrobial stewardship is the only way to reduce the spread of multi-resistant bacteria.

Delayed Laboratory Response to COVID-19 Caused by Molecular Diagnostic Contamination.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) created an exceptional situation in which numerous laboratories in Europe simultaneously implemented SARS-CoV-2 diagnostics. These laboratories reported in February 2020 that commercial primer and probe batches for SARS-CoV-2 detection were contaminated with synthetic control material, causing delays of regional testing roll-out in various countries.

West Nile neuroinvasive disease. Report of four cases in Northern Greece, 2018.

West Nile virus (WNV) is a mosquito-borne RNA flavivirus which caused several epidemics worldwide. The year 2018 was a WNV record year for Europe, including Greece, with earlier and longer transmission season with higher than the previous number of cases. It has been proposed that some simple biochemical markers may be helpful for the recognition of WNV neuroinvasive disease, its differential from other neurological infectious diseases and prognosis. We describe four cases that suffered from WNV meningitis and/or encephalitis hospitalized in 2018 in a tertiary hospital in Thessaloniki, Greece, and investigate the importance of simple biomarkers for the recognition of WNV etiology.

Laboratory readiness and response for novel coronavirus (2019-nCoV) in expert laboratories in 30 EU/EEA countries, January 2020.

Timely detection of novel coronavirus (2019-nCoV) infection cases is crucial to interrupt the spread of this virus. We assessed the required expertise and capacity for molecular detection of 2019-nCoV in specialised laboratories in 30 European Union/European Economic Area (EU/EEA) countries. Thirty-eight laboratories in 24 EU/EEA countries had diagnostic tests available by 29 January 2020. A coverage of all EU/EEA countries was expected by mid-February. Availability of primers/probes, positive controls and personnel were main implementation barriers.

West Nile virus in humans, Greece, 2018: the largest seasonal number of cases, 9 years after its emergence in the country.

BackgroundHuman cases of West Nile virus (WNV) infection are recorded since 2010 in Greece, with seasonal outbreaks occurring almost annually. Enhanced surveillance has been implemented since 2010, to promptly characterise cases' temporal and geographical distribution and inform authorities for implementation of appropriate measures (mosquito control, health education, blood safety).AimWe describe the epidemiology of WNV human infections in Greece focusing on the 2018 season.MethodsThe National Public Health Organization advised physicians to test all suspect WNV infection cases and refer samples to reference laboratories. Laboratories notified diagnosed cases on a daily basis. Treating physicians, patients, and infected blood donors were interviewed within 48 hours after diagnosis and the probable infection location was identified. Hospitalised cases were followed up until discharge.ResultsA total of 317 autochthonous WNV infection cases were diagnosed in 2018. Among them, 243 cases had neuroinvasive disease (WNND), representing a 23% increase of WNND cases compared with 2010, the previous most intense season. There were 51 deaths. Cases started occurring from week 22, earlier than usual. Both rural and urban areas were affected, with 86 (26% of the total) municipalities belonging to seven (54% of the total) regions recording cases. Two major epicentres were identified in Attica and Central Macedonia regions.ConclusionsThe largest number of human cases of WNV infection ever recorded in Greece occurred in 2018, with a wide geographical distribution, suggesting intense virus circulation. Enhanced surveillance is vital for the early detection of human cases and the prompt implementation of response measures.

Specialist laboratory networks as preparedness and response tool - the Emerging Viral Diseases-Expert Laboratory Network and the Chikungunya outbreak, Thailand, 2019.

We illustrate the potential for specialist laboratory networks to be used as preparedness and response tool through rapid collection and sharing of data. Here, the Emerging Viral Diseases-Expert Laboratory Network (EVD-LabNet) and a laboratory assessment of chikungunya virus (CHIKV) in returning European travellers related to an ongoing outbreak in Thailand was used for this purpose. EVD-LabNet rapidly collected data on laboratory requests, diagnosed CHIKV imported cases and sequences generated, and shared among its members and with the European Centre for Disease Prevention and Control. Data across the network showed an increase in CHIKV imported cases during 1 October 2018-30 April 2019 vs the same period in 2018 (172 vs 50), particularly an increase in cases known to be related to travel to Thailand (72 vs 1). Moreover, EVD-LabNet showed that strains were imported from Thailand that cluster with strains of the ECSA-IOL E1 A226 variant emerging in Pakistan in 2016 and involved in the 2017 outbreaks in Italy. CHIKV diagnostic requests increased by 23.6% between the two periods. The impact of using EVD-LabNet or similar networks as preparedness and response tool could be improved by standardisation of the collection, quality and mining of data in routine laboratory management systems.

Proficiency Testing of Virus Diagnostics Based on Bioinformatics Analysis of Simulated In Silico High-Throughput Sequencing Data Sets.

Quality management and independent assessment of high-throughput sequencing-based virus diagnostics have not yet been established as a mandatory approach for ensuring comparable results. The sensitivity and specificity of viral high-throughput sequence data analysis are highly affected by bioinformatics processing using publicly available and custom tools and databases and thus differ widely between individuals and institutions. Here we present the results of the COMPARE [Collaborative Management Platform for Detection and Analyses of (Re-)emerging and Foodborne Outbreaks in Europe] in silico virus proficiency test. An artificial, simulated in silico data set of Illumina HiSeq sequences was provided to 13 different European institutes for bioinformatics analysis to identify viral pathogens in high-throughput sequence data. Comparison of the participants' analyses shows that the use of different tools, programs, and databases for bioinformatics analyses can impact the correct identification of viral sequences from a simple data set. The identification of slightly mutated and highly divergent virus genomes has been shown to be most challenging. Furthermore, the interpretation of the results, together with a fictitious case report, by the participants showed that in addition to the bioinformatics analysis, the virological evaluation of the results can be important in clinical settings. External quality assessment and proficiency testing should become an important part of validating high-throughput sequencing-based virus diagnostics and could improve the harmonization, comparability, and reproducibility of results. There is a need for the establishment of international proficiency testing, like that established for conventional laboratory tests such as PCR, for bioinformatics pipelines and the interpretation of such results.

Evolutionary dynamics of lineage 2 West Nile virus in Europe, 2004-2018: Phylogeny, selection pressure and phylogeography.

West Nile virus (WNV) is an arbovirus causing neuroinvasive disease to humans and equines. Since 2004, lineage 2 WNV strains have been identified in Europe and have been implicated in severe outbreaks, with that of 2018 exceeding the total number from the previous seven years. The aim of this study was to explore the evolutionary process that shapes the genetic diversity of lineage 2 WNV strains (belonging to the Central European/Hungarian subclade) and reconstruct the origin and transmission routes in Europe, and especially in the Balkans. For this purpose, a high number of whole genome sequences (WGSs) were analyzed, along with newly characterized sequences, including strains from the 2018 WNV transmission season in Greece. Maximum likelihood and Bayesian inference methods were used to perform the phylogenetic and phylodynamic analyses and phylogeographic reconstruction. The majority of the Central European/Hungarian lineage 2 strains are grouped in 2 phylogenetic subgroups (Central/South-West European and Balkan) with bush-like topology. Purifying selection shapes their evolution, however, strong evidence of positive selection was revealed in 7 non-structural protein codons of NS1, NS4B and NS5. Thirty-two amino-acid substitutions were fixed in different phylogenetic subgroups, indicating that random genetic drift is responsible for the majority of evolutionary changes. Virus migration, followed by subsequent local evolution is responsible for continuously evolving strains throughout Europe. In total, 10 virus transitions between discrete geographical locations, involving virus spread from Central Europe to other regions, were highly supported. Three novel, independent introductions from Hungary and Bulgaria were responsible for the 2018 re-emergence of WNV in Northern Greece, indicating that Hungary remains an important ecological niche for the virus and has a central role for the dissemination of novel strains in the Balkans. In Northern Greece, tMRCA estimations indicated that a 1-to 2-year period of silent enzootic transmission precedes spread to dead-end hosts. Reconstruction of WNV population dynamics, from WGS data, revealed epidemic patterns characterized by 3- to 5-year oscillations in Europe. Future studies are necessary to determine the possible driving factors for these fluctuations i.e. avian herd immunity and climatic conditions affecting mosquito and bird populations. Maintaining adequate epidemiological surveillance with emphasis on obtaining WGS data, in areas at risk, is crucial for understanding the epidemiology and transmission patterns of WNV. It can further support integrated programs for risk assessment of virus circulation dynamics, aiming to targeted prevention and response measures for veterinary and public health in Europe.