- Short report
- Open Access
West Nile virus lineage 2 in Romania, 2015–2016: co-circulation and strain replacement
Parasites & Vectorsvolume 11, Article number: 562 (2018)
West Nile virus (WNV) is endemic in southeastern Romania and, after the unprecedented urban epidemic in Bucharest in 1996 caused by lineage 1 WNV, cases of West Nile fever have been recorded every year. Furthermore, a new outbreak occurred in 2010, this time produced by a lineage 2 WNV belonging to the Eastern European clade (Volgograd 2007-like strain), which was detected in humans and mosquitoes in the following years.
We report here, for the first time, the emergence, in 2015, of lineage 2 WNV belonging to the monophyletic Central/Southern European group of strains which replaced in 2016, the previously endemized lineage 2 WNV Volgograd 2007-like strain in mosquito populations. The emerged WNV strain harbors H249P (NS3 protein) and I159T (E glycoprotein) substitutions, which have been previously associated in other studies with neurovirulence and efficient vector transmission.
In 2016, both early amplification of the emerged WNV and complete replacement in mosquito populations of the previously endemized WNV occurred in southeastern Romania. These events were associated with a significant outbreak of severe West Nile neuroinvasive disease in humans.
West Nile virus (WNV) is a mosquito-borne flavivirus and has been circulating in Romania since 1996 when an outbreak with 352 human cases of West Nile neuroinvasive disease (WNND) occurred in southeastern Romania . The outbreak was caused by a WNV strain belonging to genetic lineage 1, and the epidemic vector was identified as Culex pipiens . A second significant outbreak (54 cases of WNND in humans), occurred in 2010 and was caused by a WNV of the genetic lineage 2, similar to the strain causing an outbreak in 2007 in Volgograd, Russia . Molecular evidence for the maintenance of this Eastern European WNV for another three transmission seasons (2011–2013) was obtained from mosquito vectors and human patients, suggesting its endemization in southeastern Romania . WNV vertical transmission in Culex pipiens was documented for the first time in Europe , and the same virus strain was detected in a Hyalomma marginatum tick collected from a song thrush , suggesting virus maintenance and introduction mechanisms. In 2016, a third significant WNND outbreak in humans was recorded in Romania, with 93 neurological cases .
Here, we report, for the first time, detection in Romania, in 2015, of a lineage 2 WNV, belonging to the Central/Southern European monophyletic group of strains [7, 8], co-circulating with previously endemized lineage 2 virus (Volgograd 2007-like strain), and the replacement of the latter by the emerged WNV, in 2016, when a significant outbreak occurred.
Mosquitoes were collected and processed for WNV detection as previously described , from May to September 2015–2016, and from two different habitats in southeastern Romania: urban Bucharest, and a natural lagoon habitat, Tulcea county. Two human samples were included in the study and analyzed as previously described : a urine sample, from a WNND case occurring in central Romania, Sibiu county, in 2015, and a plasma sample from an asymptomatic blood donor, detected in 2016, during the outbreak in Bucharest (Table 1 and Fig. 1). Partial NS5 and E glycoprotein sequences were obtained and used for genotyping and phylogeny [4, 10], and partial E and NS3 sequences for identification of non-synonymous substitutions [10, 11]. The other sequences used for inferring phylogeny were retrieved from GenBank.
In the transmission seasons 2015–2016, we collected and analyzed 19,390 mosquitoes distributed across 746 pools. Most of them (80%) were captured in Tulcea and only 20% in Bucharest. In Bucharest and Tulcea, respectively, the species composition was: Culex pipiens (s.l.) (86.97% and 62.52%); Culex modestus (7.57% and 3.59%); and Coquillettidia richiardii (5.46% and 33.89%).
In 2015, the TaqMan assay used to detect WNV genome was positive in four out of 126 mosquito pools from Bucharest (three Culex pipiens and one Culex modestus positive pools, 2864 mosquitoes analyzed, minimum infection rate/1000 (MIR) = 1.40), and in 15 out of 352 mosquito pools from Tulcea (ten Culex pipiens and five Coquillettidia richiardii positive pools, 9385 mosquitoes analyzed, MIR = 1.60). In 2016, nine out of 41 mosquito pools were positive in Bucharest (seven Culex pipiens and two Culex modestus positive pools, 1072 mosquitoes analyzed, MIR = 8.40), and 25 out of 227 in Tulcea (25 Culex pipiens positive pools, 6069 mosquitoes analyzed, MIR = 4.12). Seventeen out of these positive mosquito pools yielded amplification products for partial sequences of NS5, twelve pools for NS3, and three pools for E region (Table 1). In 2015, WNV infected mosquitoes were first detected on July 9 and 17 in Bucharest and Tulcea, respectively, while in 2016 the first infected mosquitoes were detected on June 21 in Bucharest, and June 6 in Tulcea, revealing early virus amplification, which led to higher mosquito infection rates in the outbreak year.
In 2015, in the urine sample of a WNND patient, a resident of Sibiu county, central Romania, a WNV belonging to the Central/Southern European clade of genetic lineage 2 was detected. The strain was highly similar to the group of strains having the Nea Santa-Greece-2010 WNV lineage 2 strain as its prototype which has been circulating in the Balkans since 2010 (Table 1, Figs. 1 and 2).
In southeastern Romania, in 2015, the endemically established Volgograd 2007-like strain was found co-circulating in mosquito populations with a WNV strain, also closely related to the Nea Santa 2010 group of strains. However, in 2016, only this newly emerged WNV strain was detected in mosquitoes in the study sites, and a very similar isolate was detected in a blood donor sample collected during the outbreak peak in Bucharest (Table 1, Figs. 1 and 2). All the obtained NS3 sequences of these newly detected isolates displayed the H249P substitution, as seen in the Nea Santa strain, but the E fragment differed from the latter by I159T substitution (Table 1). Despite the limited resolution of NS5 fragment used for inferring phylogeny, the analysis showed that all the isolates of the newly detected virus in southeastern Romania form a cluster, which does not include the WNV derived from the patient in Sibiu, central Romania (Table 1, Fig. 2).
Lineage 2 WNV was first detected in Romania as Eastern European clade  and caused the human outbreak of WNND in 2010 . Molecular evidence for its maintenance between 2011–2015, in mosquitoes and humans, in southeastern Romania, was provided by Dinu et al.  and the present study. Our results showed that in this region, the emerged WNV strain belonging to the Central/Southern European genetic lineage 2 was first found in 2015 in mosquitoes, co-circulating with the Volgograd 2007-like strain, while in 2016 it completely replaced the latter in mosquito populations. The emerged WNV strain was also detected in a human sample collected in Bucharest at the peak of the outbreak, suggesting its involvement in the 2016 epidemic. The records of the European Centre for Disease Prevention and Control  showed that between July 1 - October 10, 2016, 93 human WNND cases were reported in Romania, most of them (68 cases) distributed in the southeast of Romania, including Bucharest city (18 cases); re-emergence of severe WNND cases with high mortality was also observed in this area . Only four WNND cases were recorded beyond the Carpathian Mountains, in Transylvania, in 2016 , where a WNV of the Central/Southern European lineage 2 was detected in 2015 and described in the present study. However, this isolate is distinct from the isolates of 2015–2016 southeastern Romania WNV cluster. Co-circulation of different strains of the Central/Southern European clade was also reported in Serbia . Multiple introduction events and local evolution may lead to a high diversity of WNV, as described in Italy, where co-circulation of distinct lineages and different clades of lineage 2 was documented [8, 13].
Molecular evolution with impact on viral fitness and neurovirulence in lineage 2 WNV strains mirrors that of lineage 1 . In both lineage 2 strains circulating in Romania, the E glycoprotein presents an N-linked glycosylation site (position 154), which was previously related to neuroinvasiveness , and efficient virus transmission in vectors . In the strains emerged in 2015–2016, in southeastern and central Romania, at amino acid position 159 (E glycoprotein), isoleucine was substituted by threonine. In lineage 1, V159A substitution was related to the emergence of the North American genotype WN02 which displaced the introduced NY 1999 strain and was linked to efficient vector transmission and dissemination . Threonine in this position was found in an Italian lineage 2 cluster and was predicted to be subject of positive selection . H249P substitution in NS3 was associated with neurovirulence and increased and sustained viremia in birds, leading to high epizootic potential .
Lineage 2 WNV strains belonging to the monophyletic Central/Southern European group of strains and harboring H249P (NS3 protein), and I159T (E glycoprotein) substitutions, were for the first time detected in Romania in 2015. The emerged WNV in southeastern Romania formed a distinct cluster from the isolate detected in central Romania in 2015. In 2016, early amplification of the emerged WNV, and complete replacement in mosquito populations of the previously endemized virus in southeastern Romania occurred. These events were associated with a human outbreak of severe, high fatality WNND. Although limited, the data obtained in this study contribute to the whole picture of the adaptive evolution in Europe of the WNV clade derived from the Hungarian prototype strain.
Minimum infection rate/1000
Nonstructural protein 3
Nonstructural protein 5
West Nile neuroinvasive disease
West Nile virus
Tsai TF, Popovici F, Cernescu C, Campbell GL, Nedelcu NI. West Nile encephalitis epidemic in southeastern Romania. Lancet. 1998;352:767–71.
Savage HM, Ceianu C, Nicolescu G, Karabatsos N, Lanciotti R, Vladimirescu A, et al. Entomologic and avian investigations of an epidemic of West Nile fever in Romania in 1996, with serologic and molecular characterization of a virus isolate from mosquitoes. Am J Trop Med Hyg. 1999;61:600–11.
Sirbu A, Ceianu CS, Panculescu-Gatej RI, Vazquez A, Tenorio A, Rebreanu R, et al. Outbreak of West Nile virus infection in humans, Romania, July to October 2010. Euro Surveill. 2011;16(2).
Dinu S, Cotar AI, Panculescu-Gatej IR, Falcuta E, Prioteasa FL, Sirbu A, et al. West Nile virus circulation in south-eastern Romania, 2011 to 2013. Euro Surveill. 2015;20(20).
Kolodziejek J, Marinov M, Kiss BJ, Alexe V, Nowotny N. The complete sequence of a West Nile virus lineage 2 strain detected in a Hyalomma marginatum marginatum tick collected from a song thrush (Turdus philomelos) in eastern Romania in 2013 revealed closest genetic relationship to strain Volgograd 2007. PloS One. 2014;9:e109905.
European Centre for Disease Prevention and Control. Transmission of West Nile fever, May to November 2016 - Table of cases, 2016. https://ecdc.europa.eu/en/publications-data/table-transmission-west-nile-fever-may-november-2016-table-cases-2016. Accessed 8 Aug 2018.
Bakonyi T, Ivanics E, Erdelyi K, Ursu K, Ferenczi E, Weissenbock H, et al. Lineage 1 and 2 strains of encephalitic West Nile virus, central Europe. Emerg Infect Dis. 2006;12:618–23.
Barzon L, Papa A, Lavezzo E, Franchin E, Pacenti M, Sinigaglia A, et al. Phylogenetic characterization of Central/Southern European lineage 2 West Nile virus: analysis of human outbreaks in Italy and Greece, 2013–2014. Clin Microbiol Infect. 2015;21:1122.e1–10.
Cotar AI, Falcuta E, Prioteasa LF, Dinu S, Ceianu CS, Paz S. Transmission dynamics of the West Nile Virus in mosquito vector populations under the influence of weather factors in the Danube Delta. Romania. EcoHealth. 2016;13:796–807.
Berthet FX, Zeller HG, Drouet MT, Rauzier J, Digoutte JP, Deubel V. Extensive nucleotide changes and deletions within the envelope glycoprotein gene of Euro-African West Nile viruses. J Gen Virol. 1997;78:2293–7.
Papa A, Politis C, Tsoukala A, Eglezou A, Bakaloudi V, Hatzitaki M, et al. West Nile virus lineage 2 from blood donor. Greece. Emerg Infect Dis. 2012;18:688–9.
Popescu CP, Florescu SA, Cotar AI, Badescu D, Ceianu CS, Zaharia M, et al. Re-emergence of severe West Nile virus neuroinvasive disease in humans in Romania, 2012 to 2017 - implications for travel medicine. Travel Med Infect Dis. 2018;22:30–5.
Ravagnan S, Montarsi F, Cazzin S, Porcellato E, Russo F, Palei M, et al. First report outside eastern Europe of West Nile virus lineage 2 related to the Volgograd 2007 strain, northeastern Italy, 2014. Parasit Vectors. 2015;8:418.
McMullen AR, Albayrak H, May FJ, Davis CT, Beasley DW, Barrett AD. Molecular evolution of lineage 2 West Nile virus. J Gen Virol. 2013;94:318–25.
Shirato K, Miyoshi H, Goto A, Ako Y, Ueki T, Kariwa H, et al. Viral envelope protein glycosylation is a molecular determinant of the neuroinvasiveness of the New York strain of West Nile virus. J Gen Virol. 2004;85:3637–45.
Moudy RM, Zhang B, Shi PY, Kramer LD. West Nile virus envelope protein glycosylation is required for efficient viral transmission by Culex vectors. Virology. 2009;387:222–8.
Mann BR, McMullen AR, Swetnam DM, Barrett AD. Molecular epidemiology and evolution of West Nile virus in North America. Int J Environ Res Public Health. 2013;10:5111–29.
Langevin SA, Bowen RA, Reisen WK, Andrade CC, Ramey WN, Maharaj PD, et al. Host competence and helicase activity differences exhibited by West Nile viral variants expressing NS3-249 amino acid polymorphisms. PloS One. 2014;9:e100802.
This work was supported by the following grants of the Romanian Ministry of Research and Innovation: PN-II-PT-PCCA-2013-4-0718 Grant 126/2014, Core Program PN 16 39 01 01, and by CCCDI-UEFISCDI project number PN-III-P1.2-PCCDI-2017-0005/2018, within PNCD III.
Availability of data and materials
NS5, NS3 and E partial sequences described in this study were submitted to the European Nucleotide Archive (ENA) as part of project PRJEB21861 (acc. no. LT903700-16, LT996415-27; LT996429-32, LT996434, LT999995 and LS398545-6).
Ethics approval and consent to participate
The study was approved by the Ethics Committee of Cantacuzino National Medico-Military Institute for Research and Development.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.