Genome sequence and attenuating mutations in West Nile virus isolate from Mexico.

The complete genome sequence of a Mexican West Nile virus isolate, TM171-03, included 46 nucleotide (0.42%) and 4 amino acid (0.11%) differences from the NY99 prototype. Mouse virulence differences between plaque-purified variants of TM171-03 with mutations at the E protein glycosylation motif suggest the emergence of an attenuating mutation.

**********

Since its introduction into North America in 1999, West Nile virus (WNV) has spread rapidly across the continent, and evidence for virus circulation has also been detected in the Caribbean and parts of Central America (1). In 2003, WNV was isolated from a dead raven in Villahermosa, in the state of Tabasco, Mexico (2). Nucleotide sequencing of the premembrane (prM) and envelope (E) structural protein genes of this strain, TM171-03, and comparison with sequences from other North American isolates indicated that this virus had accumulated several unique mutations from the New York 1999 strain 382-99 (NY99) prototype sequence. We describe the complete genomic sequence of TM171-03 and its relationship to other North American isolates, as well as the results of virulence phenotype comparisons.

The Study

The isolation and initial characterization of TM171-03 have been described elsewhere (2). For genomic sequencing, RNA was extracted from infected Vero cell culture supernatant (second Vero cell passage from the original brain material, designated V2) using the QiaAmp kit (Qiagen Inc., Valencia, CA), reverse transcribed with AMV reverse transcriptase (RT) (Roche, Indianapolis, IN) and amplified by polymerase chain reaction (PCR) as nine overlapping fragments by using Taq polymerase (Roche). The PCR products were purified from 1.5% TAE/agarose gels by using the QiaQuick kit (Qiagen) and directly sequenced on an AB1 Prism model 3100 DNA sequencer (Applied Biosystems, Foster City, CA) at the University of Texas Medical Branch's Protein Chemistry Core Facility by using the amplifying primers and additional internal primers. The primers used for RT-PCR and sequencing were similar to those used by Lanciotti et al. (3) and Beasley et al. (4) (complete details are available on request). Sequence data were assembled into the complete genome sequence and analyzed as described elsewhere (2,4). In addition, Bayesian analyses were performed by using MRBAYES v3.0 (5) and 100,000 generations. A general time-reversible model was used with empirically estimated base frequencies and either a codon position-specific or a [gamma] distribution of substitution rates.

The genomic sequence for TM171-03 (GenBank accession number AY660002) differed from the NY99 prototype sequence (GenBank AF196835) at 46 nucleotides (nt) (0.42%). As reported previously, sequencing the prM-E genes of TM171-03 (V1 passage) identified nonsynonymous mutations encoding substitutions at prM-141 Ile [right arrow] Thr and E156 Ser [right arrow] Pro (2). The E-156 mutation results in the loss of the E-154-156 "NYS" glycosylation motif. Complete genome sequencing identified only two other encoded amino acid changes from the NY99 sequence at NS4B-245 (Ile [right arrow] Val) and NS5-898 (Thr [right arrow] Ile). However, during the sequencing of the V2 passage material, a reversion from Pro to Ser encoded at E-156 was observed. Analysis of the sequence chromatograms from V1 and V2 passages, and for a PCR product obtained from the original brain material, indicated that this reversion was likely to be the result of a mixed virus population, with overlapping "T" and "C" peaks visible at residue 1432 in the sense or anti-sense sequences for each product. To confirm this finding, PCR products containing the E-156 region from each passage level of TM171-03 were cloned into pGEM-T(Easy) (Promega, Madison, WI), and five clones were sequenced for each. For the original brain tissue, four clones encoded Pro at E-156 and one clone encoded Ser. For products derived from either V1 or V2 passages, two or three clones encoded a Pro at E-156, and the remainder encoded a Ser. in addition, several variants of TM171-03 were purified through two rounds of plaque selection in Veto cells, and nucleotide sequencing of these variants also confirmed a mixed population. Sequences from approximately half of the plaques encoded a Pro at E-156, while the remainder encoded Ser. Western blotting of infected Vero cell lysate antigens for these variants with WNV E protein-specific monoclonal antibody 7H2 (6) showed differences in the electrophoretic mobility of the proteins consistent with the presence or absence of glycosylation (Figure 1).

[FIGURE 1 OMITTED]

Comparison of the TM171-03 nucleotide sequence with other genomic sequences of WNV strains showed that it was most closely related to strain NY00-grouse3282 (GenBank AF404755: Figure 2). The NY00-grouse3282 sequence difl'ered from NY99 at only 13 nt (0.11%), and its relationship to TM171-03 was apparently based on 10 nonstructural protein region nucleotide differences from NY99 that were shared with TM 171-03 (Table 1). None of these mutations encoded amino acid differences, and TM171-03 differed from NY00-grouse3282 at 39 other nucleotides. These were primarily additional mutations that had accumulated in the TM171-03 strain. Genomic sequence data from other East Coast U.S. isolates collected during 2000 and subsequent years are needed to attempt to establish a definitive relationship for TM171-03 with a particular North American isolate.

[FIGURE 2 OMITTED]

To assess the effects of the E-156 Ser [right arrow] Pro mutation on the virulence of TM171-03, serial 10-fold doses from 1,000 to 0.1 PFU of TM171-03 and four plaque-purified (pp) substrains (TM171-03-pp1 and -pp2 encoding Pro at E-156; TM171-03-pp5 and -pp6 encoding Ser) were inoculated intraperitoneally (i.p.) and intracranially (i.e.) into groups of 3- to 4-week-old female NIH Swiss mice to determine mouse neuroinvasiveness and neurovirulence, as described elsewhere (7) and in accordance with guidelines of the University of Texas Medical Branch Institutional Animal Care and Use Committee.

TM171-03 was highly virulent following i.p. and i.c. inoculation, as were the plaque-purified variants, TM171-03-pp5 and -pp6, which encoded the E154-156 NYS glycosylation motif (i.p. and i.c. 50% lethal dose [L[D.sub.50]] values for each [less than or equal to] 2.0 PFU; Table 2). The lethality of these strains was comparable to that of other North American isolates that have been evaluated by using the NIH Swiss mouse model (4,7). In contrast, the nonglycosylated variants TM171-03-pp1 and -pp2 were both attenuated, having i.p. L[D.sub.50] values >1,000 PFU and i.c. L[D.sub.50] values of 32 and 25 PFU, respectively.

To confirm that the mouse virulence differences between the plaque-purified variants could be primarily attributed to the mutation at E-156, regions that encoded the additional consensus amino acid mutations at prM-141, NS4B-245, and NS5-898 were sequenced. All four plaque purified variants encoded the three amino acid mutations that were present in the consensus sequence. No additional mutations encoding amino acid changes were identified in the regions that were sequenced for these strains (equivalent to [approximately equal to] 3,000 nt in total for each). Although the entire genome of each plaque-purified variant was not sequenced, we believe that it is highly unlikely that the mouse virulence differences observed between the variants would be attributable to other amino acid mutations in the unsequenced regions that were present in the two E-156 Pro variants but not the E-156 Ser variants or the parental TM171-03 strain.

Conclusions

These results are somewhat contrary to previously reported data that described attenuated variants with glycosylated E proteins that were derived from a virulent, nonglycosylated Israeli lineage 1 WNV strain (8). However, subsequent studies identified E glycosylated variants of the same strain that retained a virulent phenotype, which suggests that multiple determinants, most probably including mutations in the nonstructural protein genes, were responsible for the observed variations in virulence (9). Other comparisons of wild-type WNV strains suggested that absence of E protein glycosylation might be associated with attenuation of mouse neuroinvasiveness (7). Recently, some of us have shown that mutating the E protein gene of a WNV infectious clone derived from the NY99 prototype strain to prevent glycosylation resulted in a [approximately equal to]200-fold attenuation of neuroinvasiveness, but not neurovirulence, in the NIH Swiss mouse model (D.W.C. Beasley, et al., unpub. data). Given the greater degree of attenuation of neuroinvasiveness and neurovirulence observed for the nonglycosylated TM171-03-pp1 and -pp2 variants described here, we hypothesize that one or more of the other mutations (at prM-141, NS4B-245, or NS5-898) also contributed to the phenotype, but this hypothesis remains to be determined experimentally. All of these mutations in the absence of the E-156 Ser [right arrow] Pro mutation (as occurred in the TM171-03-pp4 and -pp5 variants) did not appear to significantly affect the mouse virulence phenotype.

E protein glycosylation appears to play an important role in flavivirus assembly in mammalian cell culture (10); the mechanism by which this particular mutation would emerge in a wild-type WNV population, as is the case with the TM 171-03 isolate, is not clear. However, the posttranslational processing of glycoproteins differs between mosquito and mammalian cells (11), and adaptation of dengue virus to mosquito cells resulted in loss of the equivalent glycosylation motif (12), which suggests that the presence of carbohydrate on the E protein may be of lesser importance during virus replication in mosquito cells.

Recent data from the Mexican Department of Health indicate that no human cases of encephalitis attributable to local transmission of WNV have occurred during 2004 and, although many WNV-seropositive horses have been identified, few cases of overt clinical disease have been reported (http://www.cenave.gob.mx/vom accessed 26 Aug, 2004). The epidemiology of WNV disease in Mexico is likely to be complicated by preexisting immunity to other flaviviruses, but the emergence of an attenuated WNV strain would be important. WNV nucleotide sequences obtained from infected horses in Mexican states closer to the U.S. border suggest that they are closely related to recent isolates from Texas that do not encode mutations at the E glycosylation motif (13; J. Estrada-Franco, et al., unpub. data).We are unaware of any other WNV isolates from Tabasco or other southern regions of Mexico. Obtaining additional isolates from southern Mexico is important to determine if a nonglycosylated WNV population is emerging and to ascertain what impact this may have on the prevalence of severe WNV disease in Mexico.

Table 1. Summary of nucleotide and amino acid differences
between West Nile virus strains NY99 (AF196835), NY00-grouse3282
(AF404755) and TM171-03 (a)

Nucleotide
(amino acid)        NY99       NY00      TIM 171-03

50                    A          A           G
71                    A          A           G
93                    C          T           C
381                   C          C           T
483                   C          C           T
858                   C          C           T
887 (prM-141)      T (Ile)    T (Ile)     C (Thr)
1137                  C          C           T
1285                  C          T           C
1432 (E-156)b      T (Ser)    T (Ser)     C (Pro)
1626                  C          C           T
2328                  C          C           T
2388                  C          C           T
2466                  C          C           T
2607                  T          T           C
2832                  T          T           C
2865                  C          C           T
3111                  G          G           A
4146                  A          G           G
4212                  T          T           A
4564                  T          C           C
4749                  C          C           T
6120                  C          C           T
6138                  C          T           T
6141                  C          C           T
6279                  G          G           A
6426                  C          T           T
6495                  G          G           A
6771                  A          G           G
7015                  T          C           C
7359                  C          C           T
7648 (NS4B-245)    A (Ile)    A (Ile)     G (Val)
7672                  C          C           T
7938                  T          C           C
8067                  A          G           A
8109                  C          C           T
8676                  A          A           G
8811                  T          C           C
8838                  T          T           C
8994                  T          T           C
9378                  T          T           C
9408                  T          T           C
9453                  C          C           T
10317                 C          C           T
10373 (NS5-898)    C (Thr)    C (Thr)     T (Ile)
10393                 C          T           T
10828                 T          T           G
10851                 A          G           G
10989                 G          G           A

(a) Bold text indicates residues at which NY00-grouse 3282 and
TM 171-03 both differed from NY99.

(b) Strain TM171-03 had a mixed sequence at this residue. Consensus
sequence of early passage material had "C" at nucleotide 1432.

Table 2. Neuroinvasiveness and neurovirulence of TM-171 Mex03
parental isolate and plaque-purified variants after injection
into 3- to 4-week-old female NIH Swiss mice (a)

Virus          E154-156 sequence   i.p. L[D.sub.50] (PFU)

TM171-03             NYP/S                  1.3
TM171-03-pp1          NYP                  >1000
TM171-03-pp2          NYP                  >1000
TM171-03-pp5          NYS                   2.0
TM171-03-pp6          NYS                   2.0

Virus          i.p. AST [+ or -] SD (d) (b)   i.c. L[D.sub.50] (PFU)

TM171-03             7.9 [+ or -] 0.7                  0.8
TM171-03-pp1                NA                        32
TM171-03-pp2                NA                        25
TM171-03-pp5         9.0 [+ or -] 1.4                  2.0
TM171-03-pp6         8.5 [+ or -] 1.7                  1.3

Virus          i.c. AST [+ or -] SD (d) (b)

TM171-03             6.2 [+ or -] 1.8
TM171-03-pp1         6.0 [+ or -] 0.9
TM171-03-pp2         6.7 [+ or -] 1.7
TM171-03-pp5         7.4 [+ or -] 1.2
TM171-03-pp6         7.4 [+ or -] 0.9

(a) NIH, National Institutes of Health; AST, average survival
time; i.p., intraperitoneal; i.c., intracranial; LI[D.sub.50],
50% lethal dose; SD, standard deviation; NA, not applicable.

(b) Average survival time [+ or -] SD was calculated for all animals
that died following inoculation with 1,000 - 0.1 PFU doses of
indicated virus.

Acknowledgments

We thank Igor Romero Sosa and Pedro Paz-Salazar for their contributions to this project.

This study was funded in part by the State of Texas Advanced Research Program and NIH contract NO1-AI25489. Resources for field collections were provided by the Office of Exotic Diseases of the Agricultural Ministry of Mexico. D.B. is the recipient of a fellowship from the James W. McLaughlin Fellowship Fund.

Dr. Beasley is a postdoctoral fellow at the Center for Biodefense and Emerging Infectious Diseases, Department of Pathology, University of Texas Medical Branch, Galveston, Texas. His research interests include the molecular epidemiology of West Nile virus in the Americas and the identification of virulence determinants of West Nile virus and other mosquitoborne flaviviruses.

References

(1.) Gould LH, Fikrig E. West Nile virus: a growing concern? J Clin Invest. 2004:113:1102-7.

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

(3.) Lanciolti 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.

(4.) Beasley DWC, Davis CT, Guzman H, Vanlandingham DL, Travassos da Rosa APA, 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.

(5.) Ronquist F, Huelsenbeck JP. MRBAYES 3: Bayesian phylogenetic reference under mixed models. Bioinformatics. 2003;19:1572-4.

(6.) Beasley DWC, Barrett ADT. Identification of neutralizing epitopes within structural domain III of the West Nile virus envelope protein. J Virol. 2002:76:13097-100.

(7.) Beasley DWC, Li L, Suderman MT, Barrett ADT. Mouse neuroinvasive phenotype of West Nile virus strains varies depending upon virus genotype. Virology. 2002;296:17-23.

(8.) Halevy M, Akov Y, Ben-Nathan D, Kobiler D, Lachmi B, Lustig S. Loss of active neuroinvasiveness in attenuated strains of West Nile virus: pathogenicity in immunocompetent and SCID mice. Arch Virol. 1994;137:355-70.

(9.) Chambers TJ, Halevy M, Nestorowicz A, Rice CM, Lustig S. West Nile virus envelope proteins: nucleotide sequence analysis of strains differing in mouse neuroinvasiveness. J Gen Virol. 1998:79:2375-80.

(10.) Lorenz IC, Kartenbeck J, Mezzacasa A. Allison SL, Heinz FX, Helenius A. Intracellular assembly and secretion of recombinant subviral particles from tick-borne encephalitis virus. J Virol. 2003:77:4370-82.

(11.) Hsieh P, Robbins PW. Regulation of asparagine-linked oligosaccharide processing. Oligosaccharide processing in Aedes albopictus mosquito cells. J Biol Chem. 1984;259:2375-82.

(12.) Lee E, Weir RC, Dalgarno L. Changes in the dengue virus major envelope protein on passaging and their localization on the three-dimensional structure of the protein. Virology. 1997:232:281-90.

(13.) Blitvich BJ, Fernandez-Salas I, Contreras-Cordero JF, Lorono-Pino MA, Marlenee NL, Diaz FJ, et al. Phylogenetic analysis of West Nile virus, Nuevo Leon State, Mexico. Emerg Infect Dis. 2004:10:1314-7.

Address for correspondence: David W. C. Beasley, Deptartment of Pathology, University of Texas Medical Branch, 301 University Blvd., Galveston, TX 77555-0609, USA: fax: [phone]: email: [e-mail]

David W.C. Beasley, * C. Todd Davis, * Jose Estrada-Franco, * Roberto Navarro-Lopez,([dagger]) Arturo Campomanes-Cortes, ([dagger]) Robert B. Tesh, * Scott C. Weaver, * and Alan D.T. Barrett *

* University of Texas Medical Branch, Galveston, Texas, USA; and ([dagger]) Comision Mexico-Estados Unidos para la Prevencion de la Fiebre Aftosa y Otras Enfermedades Exoticas de los Animales, Mexico City, Mexico

All material published in Emerging Infectious Diseases is in the public domain and may be used and reprinted without special permission: proper citation, however, is appreciated.