Phylogenetic analysis of West Nile virus, Nuevo Leon State, Mexico.

West Nile virus RNA was detected in brain tissue from a horse that died in June 2003 in Nuevo Leon State, Mexico. Nueleotide sequencing and phylogenetic analysis of the premembrane and envelope genes showed that the virus was most closely related to West Nile virus isolates collected in Texas in 2002.

**********

West Nile virus (WNV), a mosquitoborne virus in the genus Flavivirus (family Flaviviridae), was first recognized in the Western Hemisphere during an outbreak in New York in 1999 (1). WNV rapidly disseminated across North America, and its geographic range now encompasses 47 of the 48 contiguous United States (2), 7 Canadian provinces (3), and several Mexican states (4-7).

A phylogenetic analysis of the prototype New York strain (WN-NY99-flamingo382-99), isolated from a dead flamingo from the Bronx Zoo in 1999, showed that this strain was most similar to an Israeli isolate from 1998 (8). WNV isolates collected in the northeastern United States in 2000 were similar to isolates collected in 1999 (9-12). However, studies performed with WNV isolates collected after 2000 suggest that genetically distinct populations have emerged in the United States (13,14). For example, up to 12 nucleotide substitutions (0.60% divergence) were present in the premembrane and envelope protein (prM-E) genes of isolates collected from inland and southeast coastal areas of Texas in 2002 (13).

More recently, Estrada-Franco et al. (5) reported the first isolation of WNV from Mexico. The isolate (TM171-03) was from a corvid that died on May 5, 2003, in Tabasco State, southern Mexico. We identified WNV RNA in the brain of a dead horse from Nuevo Leon State, northern Mexico. Nucleotide sequencing and phylogenetic analysis of the prM-E genes showed that this WNV from Mexico was most similar to isolates collected from noncoastal areas of Texas in 2002.

The Study

Cerebellar tissue was taken from a dead 12-year-old stallion from a privately owned ranch in the municipality of Juarez in Nuevo Leon State, Mexico, approximately 240 km south of the Texas border. The horse was first observed with neurologic symptoms on June 20, 2003, and it was euthanized 7 days later. The horse had never been outside the state of Nuevo Leon and had not been vaccinated against WNV. The tissue sample was immediately placed on dry ice and transported to the biosafety-level-3 facilities at Colorado State University for processing. Although we were unable to isolate virus from the sample by passing brain homogenate in Vero cells, we successfully amplified viral RNA.

Total RNA was extracted from approximately 100 [micro]g of cerebellar tissue with Trizol reagent (Invitrogen, Carlsbad, CA), according to the manufacturer's instructions. The prM-E genes were amplified as two fragments by reverse transcription--polymerase chain reaction (RT-PCR) by using primers designed from the nucleotide sequence of the prototype WN-NY99 strain (GenBank accession no. AF 196835). PCR amplifications were performed by using Ex Taq DNA polymerase (Takara Biomcdicals, Shiga, Japan), which has 3' [right arrow] 5' exonuclease activity. Amplification products were separated by agarose gel electrophoresis, visualized with crystal violet, and extracted by using the rapid gel extraction system (Invitrogen, Carlsbad, CA). The resulting DNA fragments were reamplified by PCR because of the low RNA copy number in the original material and purified by using the QIAquick PCR purification kit (Qiagen, Valencia, CA). Purified DNAs were sequenced on both strands with an ABI 377 DNA sequencer (Davis Sequencing, Davis, CA) and eight pairs of WNV-specific primers.

The nucleotide sequence of the prM-E genes of the WNV from Nuevo Leon State, Mexico (designated MexNL-03) was submitted to GenBank (GenBank accession no. AY426741). This region comprises 2004 nucleotides and corresponds to nucleotides 466 to 2469 of the genomic RNA of the WN-NY99 strain (8). Alignment of the MexNL-03 sequence with other known sequences in the GenBank database showed that it was most closely related to the homologous regions of three WNV isolates collected in Harris County, Texas, in June 2002 (strains 119, 123, and V1151; GenBank accession numbers AY185908, AY185909, AY185911 respectively). The MexNL-03 sequence differed from the Harris County isolates in three nucleotide positions (0.15% divergence; Table). in all cases, one change was in the prM gene at position 549, and two changes were in the E gene at positions 1179 and 1356. All substitutions were in the third codon position, and none resulted in an amino acid change.

The nucleotide sequence of MexNL-03 differed from that of the WN-NY99 strain (GenBank accession no. AF196835) in six positions (0.30% divergence; Table). Two mutations were in the prM gene (positions 549 and 660), and four mutations were in the E gene (positions 1179, 1356, 1442, and 2466). The U to C substitution at 1442 resulted in an amino acid change (Val [right arrow] Ala); all other substitutions were silent. The U to C substitution at 549 and A to G substitution at 1179 have not been reported in any WNV isolates from the United States. However, an isolate collected in Illinois in 2002 (GenBank accession no. AY428521) has a U to A substitution at position 549. Similarly, an isolate from Randall County, Texas (GenBank accession no. AY428519), has an A to C substitution at position 1179. Several more divergent strains of WNV, such as a Kunjin virus isolated in Australia in 1960 (GenBank accession no. D00246), have a G at position 1179. The prM-E genes of MexNL-03 differed from TM171-03 (GenBank accession no. AY371271) in 13 nucleotide positions (0.65% divergence; Table). Five mutations were in the prM gene, and eight mutations were in the E gene. Three mutations resulted in amino acid changes.

A phylogenetic tree was constructed by Bayesian analysis using the complete prM-E gene sequences of 49 WNV strains, including this WN V from Nuevo Leon State, Mexico (Figure). Phylogenetic trees were also generated using neighbor-joining (NJ), maximum parsimony (MP), and maximum likelihood (ML) analyses (data not shown). In the Bayesian tree, the WNV isolates from North America formed a monophyletic group consisting of two sister clades (denoted as clade 1 and 2). MexNL-03 shared a close phylogenetic relationship with isolates from inland Texas, consistent with our nucleotide sequence alignments. These viruses, along with isolates from Colorado, Illinois, Alabama, and Tabasco State, belonged to a nested clade (denoted as 1A) within clade 1. The statistical support for clade 1A by parsimony and distance bootstrap analyses were 64% and 66%, respectively. The other WNV isolates that clustered in clade 1 were from the eastern United States. Clade 2 contained WNV isolates from the southeast coastal area of Texas and the northeastern United States. The WNV isolates from coastal Texas formed a nested clade (denoted as 2A), confirming previous results (13,14).

[FIGURE OMITTED]

The parsimony-informative sites, corresponding to nucleotide positions 660, 1442, and 2466, were important in defining the topologic features of the tree (data not shown; see online figure, available from http://www.cdc.gov/ncidod/EID/vol10no7/03-0959G2.htm). All 21 WNV strains in clade 1A contained a unique C to U substitution at 2466. All, except for the five most basal isolates, contained a unique C to T substitution at 660 and all, except for TM171-03, contained a unique T to C substitution at 1442. TM171-03 was basal to the other WNV strains in clade 1A, and this topologic arrangement was due to the single nucleotide difference at 1442. However, we consider it unlikely that TM171-03 was the predecessor of the other viruses in clade IA. In addition, the parsimony and distance bootstrap analyses did not significantly support this topologic arrangement (bootstrap values = 50% and 51%, respectively). Moreover, WNV isolates from the northeastern United States occupied the basal position of clade 1. Similarly, the basal position of clade 2 was occupied by WNV isolates from the northeastern United States. Thus, our findings suggest that the WNV isolates circulating in time United States and Mexico diverged from a common ancestor from the northeastern United States.

The trees generated by NJ, MP, and ML analyses showed the same overall topologic features to the Bayesian tree, except that all of the coastal Texas isolates were basal to the Israeli isolate (data not shown). The bootstrap support for this topologic arrangement ranged from 58% to 73%. Similar findings were also reported by Estrada-Franco et al. (5). Furthermore, the WNV isolates from Maryland and New Jersey, as well as most of the isolates from New York, occupy the basal positions of clade 1 of the NJ/MP/ML trees. As a result, we have not shown the bootstrap values for clades 1 and 2 of the Bayesian tree because their composition does not match exactly to the corresponding clades of the NJ/MP/ML trees. However, the Bayesian analysis provides a more robust and efficient phylogenetic tool compared to more conventional phylogenetic techniques (15). Additional sequencing and phylogenetic analyses will be necessary to clarify these issues.

Conclusions

The data presented here indicate that WNV was introduced into Nuevo Leon State, Mexico, from inland Texas. A likely mode of introduction was by infected birds traveling for short distances (16). Earlier studies have provided serologic evidence of WNV infection in horses or birds in the nearby Mexican states of Coahuila, Tamaulipas, and Chihuahua (4,5,17). Taken together, our sequence data and the findings from the serosurveys indicate that WNV activity is now widespread in northern Mexico, as well as in other regions in Mexico (5,6,18). The geographic distribution of WNV in the Americas will likely continue to expand; thus, enhanced WNV surveillance in Mexico is warranted.

Table. Nucleotide and deduced amino acid differences in the premembrane
and envelope genes of the West Nile virus from Nuevo Leon State
compared with various other West Nile viruses

                               Nucleotide no. (a)

Strain        Geographic
              origin       483   549   660   858      887

WN-NY99 (b)   New York,     C     U    C      C        U
              NY                                   (Ile) (c)
MexNL-03 (d)  Nuevo Leon          C     U
              State,
              northern
              Mexico
TM171-03 (e)  Tabasco       U                 U       C
              State,                                 (Thr)
              southern
              Mexico
119 (f)       Harris Co.,               U
              inland
              Texas
                    Nucleotide no. (a)

Strain        1137  1179  1356  1432   1442

WN-NY99 (b)    C     A     C      U      U
                                (Ser)  (Val)
MexNL-03 (d)         G     U             C
                                       (Ala)
TM171-03 (e)   U                  C
                                (Pro)
119 (f)                                  C
                                       (Ala)

Strain        1626  2328  2388  2466   Ref.

WN-NY99 (b)    C     C     C      C      8
MexNL-03 (d)                      U
TM171-03 (e)   U     U     U      U      5
119 (f)                           U     13

(a) Nucleotide numbers correspond to those of the prototype New York
strain (WN-NY99).

(b) Isolated from a Chilean flamingo (Phoenicopterus chilensis)
(collection date: 06/01/99). GenBank accession no.: AF196835.

(c) Amino acid changes are shown in parentheses. Gaps indicate no
change in the nucleotide sequence.

(d) Isolated from a horse (collection date: 06/27/03). GenBank
accession no.: AY426741.

(e) Isolated from Common Raven (collection date: 05/05/03). GenBank
accession no.: AY371271.

(f) Isolated from a bluejay (Cyanocitta cristata) and passaged once in
Vero cells (collection date: 06/14/02). GenBank accession no. AY185908.

Acknowledgments

We thank Alan Barrett and Scott Weaver for providing unpublished WNV sequences that were used in the phylogenetic analyses and Charles Calisher for critically reading the manuscript.

This study was supported by grant U50 CCU820510 from the Centers for Disease Control and Prevention and in part by grant A145430 and contract N01-A125489 from the National Institutes of Health.

References

(1.) Nash D, Mostashari F, Fine A, Miller J, O'Leary D, Murray K, et al. The outbreak of West Nile virus infection in the New York City area in 1999. N Engl J Med. 2001;344:1807-14.

(2.) Centers for Disease Control and Prevention. West Nile virus--statistics, surveillance and control, September 30, 2003. [cited 2004 May 14]. Available from: http://www.cdc.gov/ncidod/dvbid/westnile/ surv&control.htm

(3.) Health Canada. Population and Public Health Branch WNV surveillance updates, September 30, 2003. [cited 2004 May 14]. Available from: http://www.hc-sc.gc.ca/pphb-dgspsp/wnv-vwn/mon_e.html #sitrep.

(4.) Blitvich BJ, Fernandez-Salas I, Contreras-Cordero JF, Marlenee NL, Gonzalez-Rojas JI, Komar N, et al. Serologic evidence of West Nile virus infection in horses, Coahuila State, Mexico. Emerg Infect Dis. 2003;9:853-6.

(5.) Estrada-Franco JG, Navarro-Lopez R, Beasley DWC, Coffey L, Carrara A-S, Travassos da Rosa A, et al. West Nile virus in Mexico: serologic evidence of widespread circulation since July 2002. Emerg Infect Dis. 2003;9:1604-7.

(6.) Lorono-Pino MA, Blitvich BJ, Farfan-Ale JA, Puerto FI, Blanco JM, Marlenee NL, et al. Serologic evidence for West Nile virus infection in horses, Yucatan State, Mexico. Emerg Infect Dis. 2003;9:857-9.

(7.) Ulloa A, Langevin SA, Mendez-Sanchez JD, Arredondo-Jimenez JI, Raetz JL, Powers AM, et al. Serologic survey of domestic animals for zoonotic arbovirus infections in the Lacandon Forest region of Chiapas, Mexico. Vector Borne Zoonotic Dis. 2003;3:3-9.

(8.) Lanciotti RS, Roehrig JT, Deubel V, Smith J, Parker M, Steele K, et al. Origin of the West Nile virus responsible for an outbreak of encephalitis in the northeastern United States. Science. 1999;286:2333-7.

(9.) Lanciotti RS, Ebel GD, Deubel V, Kerst AJ, Murri S, Meyer R, et al. Complete genome sequences and phylogenetic analysis of West Nile virus strains isolated from the United States, Europe, and the Middle East. Virology. 2002;298:96-105.

(10.) Briese T, Glass WG, Lipkin WI. Detection of West Nile virus sequences in cerebrospinal fluid. Lancet. 2000;355:1614-5.

(11.) Anderson JK Vossbrinck CR, Andreadis TG, Iton A, Beckwith WH. Mayo DR. A phylogenetic approach to following West Nile virus in Connecticut. Proe Natl Acad Sci U S A. 2001;98:12885-9.

(12.) Ebel GD, Dupuis AP, Ngo K, Nicholas D, Kauffman E, Jones SA, et al. Partial genetic characterization of West Nile virus strains, New York State, 2000. Emerg Infect Dis. 2001;7:650-3.

(13.) Beasley DW, Davis CT, Guzman H, Vanlandingham DL, Travassos da Rosa AP, Parsons RE, et al. Limited evolution of West Nile virus has occurred during its southwesterly spread in the United States. Virology. 2003;309:190-5.

(14.) Davis TC, Beasley DWC, Guzman II Raj P, D'Anton M, Novak RJ, et al. Genetic variation among temporally and geographically distinct West Nile virus isolates collected in the United States during 2001 and 2002. Emerg Infect Dis. 2003;9:1423-9.

(15.) Huelsenbeck JP, Ronquist FR. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17:754-5.

(16.) Rappole JH, Hubalek Z. Migratory birds and West Nile virus. J Appl Microbiol. 2003;94(Suppl):47S-58S.

(17.) Fernandez-Salas I, Contreras-Cordero JF, Blitvich BJ, Gonzalez-Rojas JI, Cavazos-Alvarez A, Marlenee NL, et al. Serologic evidence of West Nile virus infection in birds, Tamaulipas State, Mexico. Vector Borne Zoonotic Dis. 2003;3:209-13.

(18.) Farfan-Ale JA, Blitvich B J, Lorono-Pino MA, Marlenee NL, Rosado-Paredes EP, Garcia-Rejon JE, et al. Longitudinal studies of West Nile virus infection in avians, Yucatan State, Mexico. Vector Borne Zoonotic Dis. 2004;4:3-14.

Address for correspondence: Barry J. Beaty, Arthropod-borne and Infectious Diseases Laboratory, Department of Microbiology, Immunology and Pathology, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, Colorado 80523, USA; fax: [phone]; email: [e-mail]

Bradley J. Blitvich,* Ildefonso Fernandez-Salas, ([dagger]) Juan F. Contreras-Cordero, ([dagger]) Maria A. Lorono-Pino, * Nicole L. Marlenee, * Francisco J. Diaz, * Jose I. Gonzalez-Rojas ([dagger]) Nelson Obregon-Martinez, ([double dagger]) Jorge A. Chiu-Garcia, ([double dagger]) William C. Black IV, * and Barry J. Beaty *

* Colorado State University, Fort Collins, Colorado. USA; ([dagger]) Universidad Autonoma de Nuevo Leon, Apartado, San Nicolas de los Garza, Nuevo Leon, Mexico; and ([double dagger]) Secretaria de Salud de Nuevo Leon, Nuevo Leon, Mexico

Dr. Blitvich is a postdoctoral scientist in the Department of Microbiology, Immunology and Pathology at Colorado State University, Fort Collins, Colorado. His research interests include the mechanisms of vector and host interactions in arbovirus transmission cycles.