Emerg Infect DisEmerging Infect. DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention21529408332178110-156910.3201/eid1705.101569DispatchPhylogenetic Analysis of West Nile Virus Isolates, Italy, 2008–2009Running head: West Nile Virus Isolates, Italy, 2008–2009RossiniGiada1CarlettiFabrizio1BordiLiciaCavriniFrancescaGaibaniPaoloLandiniMaria P.PierroAnnaCapobianchiMaria R.Di CaroAntoninoSambriVittorioAuthor affiliations: University of Bologna, Bologna, Italy (G. Rossini, F. Cavrini, P. Gaibani, M.P. Landini, A. Pierro, V. Sambri);National Institute for Infectious Diseases “L. Spallanzani,” Rome, Italy (F. Carletti, L. Bordi, M.R. Capobianchi, A. Di Caro)Address for correspondence: Vittorio Sambri, Alma Mater Studiorum, Università di Bologna, Department of Hematology and Oncology, Section Microbiology, Via Massarenti 9, Bologna 40138, Italy; email: vittorio.sambri@unibo.it52011175903906

To determine the lineage of West Nile virus that caused outbreaks in Italy in 2008 and 2009, several West Nile virus strains were isolated from human specimens and sequenced. On the basis of phylogenetic analyses, the strains isolated constitute a distinct group within the western Mediterranean cluster.

Keywords: VirusesWest Nile virusflavivirusvector-borne infectionsphylogeneticsItalydispatch

West Nile virus (WNV) is an arthropod-borne flavivirus affecting a wide range of vertebrates, including birds and mammals. The natural cycle of infection involves birds and mosquitoes (1); many species of wild birds act as amplifying hosts, whereas humans and horses are considered dead-end hosts (2). First identified in tropical Africa, WNV infection has been found in northern Africa, Israel, India, and Australia (3) and has progressively spread in the Americas since 1999 (4,5). WNV has been the cause of outbreaks and sporadic cases in central, eastern, and Mediterranean Europe for >45 years (6). Most strains responsible for the European and Mediterranean basin outbreaks are in lineage 1, with most grouped in the so-called European Mediterranean/Kenyan cluster (7,8).

In August 2008, an outbreak involving wild birds, horses, and humans affected 8 provinces in 3 regions in Italy (Emilia-Romagna, Veneto, Lombardy) (9,10). In 2009, a new epidemic was reported in the same region and in other neighboring regions in Italy, with up to 17 confirmed human cases of West Nile neuroinvasive disease (11).

The Study

We obtained nucleotide sequences of the envelope (E) and nonstructural (NS)3/NS5 protein regions of 6 WNV human strains isolated during the 2008–2009 outbreaks in Italy and analyzed their phylogenetic relationships with other WNV strains isolated in Europe and the Mediterranean basin. Virus isolation on Vero E6 cells was performed from human serum samples obtained from patients with neuroinvasive WNV infection: 5 samples were collected during the 2009 outbreak from patients residing in the provinces of Ferrara (2 patients), Modena (1), and Mantua (2); 1 sample was collected during the 2008 outbreak from a patient living in Ferrara. Cytopathic effects were observed 6 days after the serum was injected into cell cultures; WNV infection was identified by immunofluorescence staining and confirmed by real-time reverse transcription–PCR (9).

Virus isolates were designed as follows: TOS-09 (Ferrara), NAL-09 (Ferrara), CHI-09 (Mantua), BAL-09 (Mantua), FAN-09 (Modena), and MAN-08 (Ferrara). After 1 passage, supernatants of infected tissue cultures were processed for viral RNA extraction and PCR amplification. Nucleotide sequences were obtained for the E gene (1,503 nt) and for a region of 5,759 nt encompassing the NS3, NS4a, NS4b, and NS5 genes, which are considered to be the most phylogenetic informative regions (12,13). Overlapping amplicons were amplified and bidirectionally sequenced. Primer sequences are available upon request.

Multiple alignments of the 6 new WNV sequences of E and NS3/NS5 regions, along with the corresponding sequences from 30 other WNV strains and isolates available in GenBank and Usutu virus as outgroup virus, were generated by using ClustalW software (www2.ebi.ac.uk/clustalw). Phylogenetic analyses of the E (Figure 1) and NS3/NS5 (Figure 2) sequences generated highly comparable topologies and confirmed that strains from Italy and all western Mediterranean strains belonging to lineage 1, clade 1a, were highly related. In addition, on the basis of the phylogenetic trees of E and NS3/NS5 sequences, we found that the 6 new sequences clustered with other strains isolated in Italy in 2008 and 2009, constituting a distinct and well-supported group within the western Mediterranean cluster (100% bootstrap on E and NS3/NS5 regions) (Figure1, Figure 2).

Phylogenetic trees of West Nile virus strains isolated during outbreaks in Italy, 2008–2009, based on nucleotide sequences of the complete envelope gene. Phylogenetic tree and distance matrices were constructed by using nucleotide alignment, the Kimura 2-parameter algorithm, and the neighbor-joining method implemented in MEGA version 4.1 (www.megasoftware.net/mega4/mega41.html). The tree was rooted by using Usutu virus as the outgroup virus. The robustness of branching patterns was tested by 1,000 bootstrap pseudoreplications. The percentage of successful bootstrap replicates is indicated at nodes, showing only values >80. The scale bar indicates nucleotide substitutions per site. Strains sequenced in this study are TOS-09, FAN-09, BAL-09, CHI-09, NAL-09, and MAN-08. Envelope nucleotide sequence GenBank accession nos: NAL-09, HM991272; TOS-09, HM991273; BAL-09, HM991274; CHI-09, HM991275; MAN-09, HM991276; FAN-09, HM991277. Boldface indicates the 6 newly sequenced strains.

Genomic region encompassing nonstructural (NS) 3, NS4a, NS4b, and NS5 genes. Phylogenetic tree and distance matrices were constructed by using nucleotide alignment, the Kimura 2-parameter algorithm, and the neighbor-joining method implemented in MEGA version 4.1 (www.megasoftware.net/mega4/mega41.html). The tree was rooted by using Usutu virus as the outgroup virus. The robustness of branching patterns was tested by 1,000 bootstrap pseudoreplications. The percentage of successful bootstrap replicates is indicated at nodes, showing only values >80. The scale bar indicates nucleotide substitutions per site. Strains sequenced in this study are TOS-09, FAN-09, BAL-09, CHI-09, NAL-09 and 191 MAN-08. NS3/NS5 region nucleotide sequence GenBank accession nos.: NAL-09, HM641230; TOS-09, HM641225; BAL-09, HM641226; CHI-09, HM641227; MAN-09, HM641229; FAN-09, HM641228. Boldface indicates the 6 newly sequenced strains. *Strains from Italy carrying the T249P mutation.

For the E gene, the mean genetic distance within the 6 new isolates from Italy was 0.08% at the nucleotide level, whereas the derived amino acid sequences were totally conserved. For the NS3/NS5 region, the mean genetic distance within the 6 new isolates from Italy was 0.10% at the nucleotide and amino acid levels.

Within the 2008–2009 cluster from Italy, the 6 newly sequenced human strains exhibited a 0.14% mean nucleotide distance with the 2 strains isolated in the Emilia Romagna region in 2008 from magpies (Italy 15803/08 and Italy 15217/08) and 0.15% with the first human WNV strain isolated in Italy in 2009 (Italy/09) in the Veneto region (14), in the E and NS3/NS5 genomic regions. On the basis of NS3/NS5 sequences, all strains from Italy isolated during the 2008–2009 outbreaks showed the highest relatedness to Morocco/03, France/04, and Portugal/04 with 1.08%, 1.28%, and 1.45% mean nucleotide distance, respectively, whereas they show a mean nucleotide distance of 1.56% with Italy/98. The close relatedness to Morocco/03, France/04, and Portugal/04 was also confirmed on E gene sequences (1.44%, 1.51%, and 1.17% mean nucleotide distance, respectively), compared with a mean nucleotide distance of 2.06% with Italy/98.

Comparison of the amino acid sequences of the E and NS3/NS5 region showed that all 2008–2009 strains from Italy exhibit a set of amino acid substitutions in the NS4a and NS5 proteins not shared with other western Mediterranean strains that could be considered as molecular signatures, e.g., I2209 and S2224 in the NS4a protein, A2786 and K2950 in the NS5 protein (Table A1). No relevant amino acid changes were observed in the E region.

Notably, the NS3 Thr249Pro substitution, suggested to be associated with increased virulence for American crows (15), was not uniformly present in all the 2008–2009 isolates in Italy. It was detected in all 2008 isolates, but heterogeneously distributed among the 2009 isolates. In particular, among the newly sequenced strains, P at NS3–249 was present in the BAL-09, FAN-09, and MAN-08 isolates, whereas T was present in the TOS-09, NAL-09, and CHI-09 isolates. Thus, this substitution could not be considered a hallmark of the well-defined 2008–2009 cluster in Italy.

Conclusions

Although in the past decades several outbreaks of WNV infection have occurred in Europe and in the Mediterranean Basin (6), few viral isolates and nucleotide sequences are available from this geographic area, thus preventing an in-depth analysis of the molecular epidemiology of West Nile infection in the region. We provide 6 new nucleotide sequences of the E and NS3/NS5 regions from WNV strains isolated in 2008–2009 from patients in Italy with West Nile neuroinvasive disease and describe the isolates’ phylogenetic relationships with other WNV strains in the Western Mediterranean cluster.

The analysis of the E and NS3/NS5 sequences enabled a new phylogenetic reconstruction of the Western Mediterranean cluster, suggesting that the 2008–2009 strains from Italy constitute a distinct and well-supported group within the European–Mediterranean clade, characterized by specific amino acid changes in NS4a and NS5 proteins. The structure of phylogenetic trees, together with the presence of amino acid signatures, is consistent with a common origin of the strains circulating 2008–2009, supporting the concept of endemic circulation of WNV in Italy (6).

Furthermore, the T249P change in WNV-NS3 helicase has been suggested as increasing the virulence of WNV in birds in North America (15), although recent data do not support this hypothesis (8). Our results indicate that, at least in 2009 in Italy, there was a co-circulation of WNV strains carrying either a P or a T in this position. On the basis of our study, it is not possible to establish whether this was because of different virus introductions in Italy or reversion of some originally mutated strains. However, our findings do not support the hypothesis that the recent emergence of this mutation might have contributed to rapid spread of WNV infection with the occurrence of human cases in the Mediterranean Basin. More accurate studies are required to define the role of the T249P mutation in the pathogenesis and evolutionary history of WNV in recent outbreaks.

Suggested citation for this article: Rossini G, Carletti F, Bordi L, Cavrini F, Gaibani P, Landini MP, et al. Phylogenetic analysis of West Nile virus isolates, Italy, 2008–2009. Emerg Infect Dis [serial on the Internet]. 2011 May [date cited]. http://dx.doi.org/10.3201/eid1705.101569

These authors contributed equally to this article.

This study was supported by “Fondi Finalizzati Centro di Riferimento Regionale per le Emergenze Microbiologiche,” Regione Emilia Romagna, and by grants to the National Institute for Infectious Diseases (Ricerca Corrente, Ricerca Finalizzata/Progetto ordinario 07.109 [9AAF/P2]).

Comparison of amino acid substitutions of West Nile virus strains isolated in Italy, 2008–2009, to different European and Mediterranean West Nile virus strains and to reference strain EG101*
AA position in WNV polyproteinAA position in WNV proteinTOS‐09NAL‐09CHI‐09BAL‐09FAN‐09MAN‐08Italy 15803/08 (FJ483548)Italy 15217/08 (FJ483549)Italy‐09 (GU011992)Italy/98 (AF404757)Morocco/03 (AY701413)France 405/04 (DQ786572)Spain GE2o/V (FJ766332)Portugal/04 (AJ965626)EG101 (AF260968)
341E‐51AAAAAAAAAAAATAA
378E‐88PPPPPPPPPPPPSPP
443E‐153GGGGGGGGGGGRGGG
602E‐312LLLLLLLLLLFLLLL
1754NS3‐249TTTPPPPPPTTTPPP
2209NS4a‐85IIIIIIIIIVVVVVV
2224NS4a‐100SSSSSSSSSPPPPPP
2265NS4a‐141MMMLMMMMMMMMMML
2359NS4b‐86FIFFFFFFFFFFFFF
2516NS4b‐243IIITTTTTTTTTTTT
2786NS5‐258AAAAAAAAAVVVVVV
2920NS5‐292PPPPSSSSSSSSSSR
2950NS5‐422KKKKKKKKKRRRRRR
3169NS5‐641TTTTAATTTTTTTTT
3407NS5‐879IIIHIIIIIIIIIII

*The amino acid position numbers refer to the EG101 strain. GenBank accession numbers are shown in parentheses. WNV, West Nile virus; E, envelope; NS, nonstructural.

Dr Rossini is a postdoctoral researcher at the University of Bologna, Unit of Clinical Microbiology. Her research interests include emerging arboviruses transmitted by mosquitoes.

ReferencesKulasekera VL, Kramer L, Nasci RS, Mostashari F, Cherry B, Trock SC, West Nile virus infection in mosquitoes, birds, horses, and humans, Staten Island, New York, 2000. Emerg Infect Dis. 2001;7:7225 10.3201/eid0704.01042111589172Komar N. West Nile viral encephalitis. Rev Sci Tech. 2000;19:1667611189714Zeller HG, Schuffenecker I. West Nile virus: an overview of its spread in Europe and the Mediterranean basin in contrast to its spread in the Americas. Eur J Clin Microbiol Infect Dis. 2004;23:14756 10.1007/s10096-003-1085-114986160Calisher CH. West Nile virus in the New World: appearance, persistence, and adaptation to a new econiche–an opportunity taken. Viral Immunol. 2000;13:4114 10.1089/vim.2000.13.41111192287Ebel GD, Carricaburu J, Young D, Bernard KA, Kramer LD. Genetic and phenotypic variation of West Nile virus in New York, 2000–2003. Am J Trop Med Hyg. 2004;71:49350015516648Calistri P, Giovannini A, Hubalek Z, Ionescu A, Monaco F, Gavini G, Epidemiology of West Nile in Europe and in the Mediterranean basin. Open Virol J. 2010;4:293720517490Schuffenecker I, Peyrefitte CN, el Harrak M. Murri Si, Leblond A, Zeller HG. West Nile virus in Morocco, 2003. Emerg Infect Dis. 2005;11:306915752452Sotelo E, Fernandez-Pinero J, Llorente F, Aguero M, Hoefle U, Blanco JM, Characterization of West Nile virus isolates from Spain: new insights into the distinct West Nile virus eco-epidemiology in the Western Mediterranean. Virology. 2009;395:28997 10.1016/j.virol.2009.09.01319833373Rossini G, Cavrini F, Pierro A, Macini P, Finarelli A, Po C, First human case of West Nile virus neuroinvasive infection in Italy, September 2008— case report. Euro Surveill. 2008;13:pii: 19002.Calistri P, Giovannini A, Savini G, Monaco F, Bonfanti L, Ceolin C, West Nile virus transmission in 2008 in north-eastern Italy. Zoonoses Public Health. 2010;57:2119 10.1111/j.1863-2378.2009.01303.x20042066Rizzo C, Vescio F, Declich S, Finarelli AC, Macini P, Mattivi A, West Nile virus transmission with human cases in Italy, August–September 2009. Euro Surveill. 2009;14:pii: 19353.Charrel RN, Brault AC, Gallian P, Lemasson JJ, Murgue B, Murri S, Evolutionary relationship between Old World West Nile virus strains. Evidence for viral gene flow between Africa, the Middle East, and Europe. Virology. 2003;315:3818 10.1016/S0042-6822(03)00536-114585341Gray RR, Veras NM, Santos LA, Salemi M. Evolutionary characterization of the West Nile Virus complete genome. Mol Phylogenet Evol. 2010;56:195200 10.1016/j.ympev.2010.01.01920102743Barzon L, Franchin E, Squarzon L, Lavezzo E, Toppo S, Martello T, Genome sequence analysis of the first human West Nile virus isolated in Italy in 2009. Euro Surveill. 2009;14:pii: 19384.Brault AC, Huang CY, Langevin SA, Kinney RM, Bowen RA, Ramey WN, A single positively selected West Nile viral mutation confers increased virogenesis in American crows. Nat Genet. 2007;39:11626 10.1038/ng209717694056