J Gen Virol Tips for Better Browsing
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Gen Virol 87 (2006), 613-622; DOI 10.1099/vir.0.81320-0

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Li, J.
Right arrow Articles by Ng, M.-L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Li, J.
Right arrow Articles by Ng, M.-L.
Agricola
Right arrow Articles by Li, J.
Right arrow Articles by Ng, M.-L.
© 2006 Society for General Microbiology

The glycosylation site in the envelope protein of West Nile virus (Sarafend) plays an important role in replication and maturation processes

J. Li{dagger}, R. Bhuvanakantham{dagger}, J. Howe and M.-L. Ng

Flavivirology Laboratory, Department of Microbiology, 5 Science Drive 2, National University of Singapore, Singapore 117597

Correspondence
M.-L. Ng
micngml{at}nus.edu.sg


   ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
The complete genome of West Nile (Sarafend) virus [WN(S)V] was sequenced. Phylogenetic trees utilizing the complete genomic sequence, capsid gene, envelope gene and NS5 gene/3' untranslated region of WN(S)V classified WN(S)V as a lineage II virus. A full-length infectious clone of WN(S)V with a point mutation in the glycosylation site of the envelope protein (pWNS-S154A) was constructed. Both growth kinetics and the mode of maturation were affected by this mutation. The titre of the pWNS-S154A virus was lower than the wild-type virus. This defect was corrected by the expression of wild-type envelope protein in trans. The pWNS-S154A virus matured intracellularly instead of at the plasma membrane as shown for the parental WN(S)V.

{dagger}These authors contributed equally to this work. Back


   INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
West Nile virus (WNV) belongs to the genus Flavivirus within the family Flaviviridae and has caused increasing global concerns in recent years. Infection can result in fatal encephalitis associated with damage to the central nervous system in humans and animals. Like all other flaviviruses, the RNA genome of WNV is approximately 11 kb in length and encodes a single large polyprotein. The order of these proteins on the polyprotein is NH2-C-prM-E-NS1-NS2a-NS2b-NS3-NS4a-NS4b-NS5-COOH (Westaway et al., 1985Down; Chambers et al., 1990Down).

The complete genome of WN(Sarafend)V [WN(S)V] has been sequenced (GenBank accession no. AY688948). It is believed to be a human isolate from Israel (Scherret et al., 2002Down). Antigenic analysis using monoclonal antibody against WN and Kunjin viruses revealed that WN(S)V belongs to a distinctive antigenic group (Adams et al., 1995Down; Scherret et al., 2001aDown). In addition, the previously reported long deletion in the nucleotide sequence immediately downstream of the stop codon of several lineage II isolates was not found in this WN(S)V isolate. This disparity between WN(S)V and other WN viruses could probably be due to the extensive passage history in the laboratory (Scherret et al., 2002Down). In general, the grouping of WNV into two lineages is based on molecular, antigenic and epidemiology classification. Previous studies could not convincingly explain the differences in pathogenicity of different WNV isolates. Therefore, more studies of the amino acid sequence variations and the resulting phenotype differences of WN(S)V compared with that of other WNV isolates are needed.

The flavivirus envelope (E) protein is proposed to be involved in a number of biological activities including virus assembly, receptor binding and membrane fusion. It contains potential N-linked glycosylation sites (Heinz, 1986Down; Chambers et al., 1990Down). Lad et al. (2000)Down showed that the glycosylation site of the E protein is important for the native conformation of the epitopes and its proper expression. This glycosylation site is also linked to the neuroinvasiveness (Halevy et al., 1994Down; Chambers et al., 1998Down; Beasley et al., 2002Down; Shirato et al., 2004Down) and enhanced replication efficiency (Scherret et al., 2001bDown; Shirato et al., 2004Down) of several strains of WNV. The rationale behind the impaired replication efficiency of the non-glycosylated variants remains to be elucidated. This study shows the importance of the glycosylation site of the WN(S)V E protein in replication and maturation processes of the virus. The multiple roles played by the glycosylation site of the E protein of the WN(S)V are also discussed.


   METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Cells and viruses.
Vero cells were grown in medium 199 (M199) containing 10 % fetal calf serum. Baby hamster kidney (BHK)-21 and C6/36 cells were grown in RPMI 1640 and L15 media, respectively. WN(S)V, a kind gift from E. G. Westaway (Sir Albert Sakzewski Virus Research Centre, Royal Children's Hospital, Brisbane, Queensland 4029, Australia), was used as the parental virus.

Sequence analysis.
The RNA extracted from WN(S)V was used to generate PCR products for DNA sequencing. The 5'- and 3'-ends of the genome were determined by 5'- and 3'-RACE (GeneRacer; Invitrogen). ScanProsite was used to identify the N-linked glycosylation motif (NYST) in the E protein. Alignment of the nucleotide sequence of WN(S)V with other WNV strains was performed with CLUSTALX 1.8 alignment software. The phylogenetic analysis was performed on the aligned nucleotide sequences with a neighbour-joining distance method in the PHYLIP package version 3.6. Bootstrap confidence levels were calculated with 1000 replicates by using the CONSENSUS program. The topologies of the resulting trees were also examined using other algorithms such as DNA parsimony and DNA maximum-likelihood in the PHYLIP package.

Cloning.
The construction of full-length cDNA clone (pWNS) is shown in Fig. 1Down. Modified pBR322 vector was generated as described earlier (Li et al., 2005aDown). Individual fragments corresponding to the four regions of the full-length genome were separately cloned into the modified vector and served as intermediate clones before sequential ligation into the vector to form pWNS. During the cloning, one silent mutation was engineered to knock out the XbaI site in the parental genome (XbaI*, nt 9474, T to A) and served as a recombination marker. The QuikChange Site-Directed Mutagenesis kit (Stratagene) was used to introduce the substitution at Ser154 of the E protein in pWNS to generate pWNS-S154A so that the expressed E protein would not be recognized by glucosyl transferase. The fragment corresponding to the nt 1–3520 was cloned into pGEM-T-Easy vector (Promega) to generate pGEM-C-prM-E for the complementation study.


Figure 1
View larger version (20K):
[in this window]
[in a new window]
 
Fig. 1. Schematic representation of the construction of the full-length cDNA clone of WN(S)V. The genome organization and position of the four synthesized cDNA fragments, by RT-PCR, covering the full-length genome for sequential cloning are indicated. One silent mutation was engineered to knock out XbaI (XbaI*, nt 9474, T to A) in the wild-type virus to form a recombination marker. Two restrictive sites were created at both ends(BsiWI at the 5'-end and XbaI at the 3'-end) of the cDNA sequence for cloning. The full-length cDNA is positioned at the 3'-end of a T7 promoter sequence for in vitro transcription. Plasmid pBR322 was modified by replacing the AatII–SalI fragment in the tetracycline resistance gene (TCr) with a linkersequence 5'-TCGAACGTACGAACCCGATGCATAAGTGGAGTCGACAAGCTTCGCAGGACTAGTTCCGGGATCCCGGACGTCTGCTCTAGATACGT-3' to form the plasmid pBR322-linkerS containing BsiWI, AatII, SalI, NsiI and XbaI (underlined) for restrictive digestions. The resulting pBR322-linkerS is used for individual cDNA fragment cloning and sequencing before sequentially assembled into the full-length cDNA clone, pWNS.

 
In vitro transcription and transfection.
In vitro transcription was performed as described by Li et al. (2005aDown) to generate RNAs from pWNS and pWNS-S154A clones. Twenty micrograms of the purified RNA was transfected into BHK-21 cells by electroporation as described by Shi et al. (2002)Down. RNA quantification was performed using a Nanodrop ND-1000 Spectrophotometer (Peqlab). The regions spanning the S154A mutation and recombination marker were sequenced to confirm the presence of the S154A mutation and the recombination marker in the recovered viruses.

Western blot and complementation analyses.
Western blot analysis was performed as described by Bhuvanakantham & Ng (2005)Down. Briefly, Vero cells were infected with parental, pWNS and pWNS-S154A viruses for 24 h. The cell lysates were then subjected to SDS-PAGE analysis followed by Western blot analysis with anti-E monoclonal antibody (Microbix). Complementation analysis was done by transfecting the RNA synthesized from the plasmid pGEM-C-prM-E into Vero cells by electroporation as described by Brandt et al. (2001)Down. Five hours post-transfection, the cells were infected with the pWNS-S154A virus at an m.o.i. of 5·0. The expression of the glycosylated and non-glycosylated E proteins in the cell lysates was determined by Western blot analysis.

Growth kinetics of the pWNS-S154A virus.
Vero, BHK-21 and C6/36 cells were inoculated with either parental WN(S)V, pWNS or pWNS-S154A virus at an m.o.i. of 5·0. Virus was allowed to be absorbed for 1 h and productive virus in the supernatant was sampled at 6, 12, 18 and 24 h post-infection (p.i.) for plaque assay. Complementation analysis was performed as described above.

Determination of growth kinetics using transcribed RNAs.
Twenty micrograms of RNAs transcribed from the pWNS and pWNS-S154A clones was used to transfect Vero cells. After transfection, the cells were seeded into 24-well plates and the culture supernatant was sampled at 8, 16, 24, 32, 40, 48 and 56 h for plaque assay and real-time RT-PCR. Viral RNAs extracted from the culture supernatant were then reverse transcribed with Supercriptase II (Invitrogen) using a primer corresponding to nt 11000–11020. The reverse transcripts were applied to a real-time PCR assay using iTaq Supermix (Bio-Rad) with forward (5'-GGTTCGCCACATCACTACACTT-3') and reverse (5'-GATCAGTTCCGTGAAATGGTTTGA-3') primers corresponding to nt 9716–9849 of the WNV genome. The kinetics of cDNA amplification were monitored with an ABI Prism 7000 sequence detection system (Applied Biosystems) using a dual labelled probe (5'-TGGAACCTGCTGCCAGTCATACCACC-3') conjugated with 6-carboxyfluorescein at the 5'-end and 6-carboxy-tetramethylrhodamine at the 3'-end. An in vitro transcribed RNA from pWNS clone was used as a standard for quantification.

Transmission electron microscopy (TEM).
Vero cells infected with parental WN(S)V, pWNS or pWNS-S154A virus at an m.o.i. of 5·0 were processed at 14, 18 and 24 h p.i as described by Li et al. (2005aDown) and viewed under a Philips electron microscope (CM 120; BioTwin). Vero cells transfected with wild-type E protein and infected with pWNS-S154A virus (complementation) were also directly processed for TEM analysis.


   RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Sequence analysis
The complete genomic sequencing of WN(S)V showed that the genome was 11057 bp in length, and its open reading frame begins at nt 97 and encodes a 3434 aa polypeptide. The putative anchored capsid (C) protein is located between nt 97 and 465. The prM protein is encoded from nt 466 to 966, with the M protein encoded between nt 742 and 966. The E protein is encoded between nt 967 and 2469, followed by the non-structural proteins: NS1 (nt 2470–3525), NS2a (nt 3526–4218), NS2b (nt 4219–4611), NS3 (nt 4612–6468), NS4a (nt 6469–6846), putative 2K protein (nt 6847–6915), NS4b (nt 6916–7683) and NS5 (nt 7684–10398), respectively.

As an initial step to examine the genetic relationship of WN(S)V to other WNV isolates, multiple nucleotide sequence alignments were made utilizing the corresponding sequences of WNV lineages. The HNY1999 strain was chosen as a representative of lineage Ia viruses and the MRM61C strain represented lineage Ib viruses. The WNFCG strain represented lineage II viruses. The Rabensburg strain was included since it was considered as a new lineage of WNV (Bakonyi et al., 2005Down). When the complete nucleotide sequences of these virus were aligned, WN(S)V exhibited 79 % nucleotide identity to the HNY1999 and MRM61C strains, and 78 % identity to Rabensburg strain and 93 % identity to WNFCG strain. The alignment of C, E and NS5/3' untranslated region (UTR) nucleotide sequences showed 93–97 % identity with WNFCG, 78–85 % with HNY1999, 78–84 % with MRM61C and 78–79 % with Rabensburg. Taken together, these results demonstrate that WN(S)V exhibited the greatest percentage identity with the WNFCG strain.

Phylogenetic trees were constructed from the aligned nucleic acid data by utilizing the entire genome or the selected genes (C, partial E, NS5/3'-UTR) of 27 WNV strains and one Japanese encephalitis virus strain (JEV – used as an outgroup), and the analysis was performed in 1000 bootstrap replicates (Fig. 2Down). Since the complete genomic sequence for most of the lineage II viruses is not available in GenBank, WNFCG was chosen as the representative while constructing the phylogenetic tree using the complete genome, C and NS5/3'-UTR. The tree derived from the partial E gene of 227 bp region (Fig. 2aDown) included 10 additional lineage II strains in addition to WNFCG, since the genomic sequence corresponding to nt 436–663 of the E gene is available for many lineage II strains. The topology of the tree showed two major lineages. Lineage I viruses form a tight cluster with close genetic relationship among the members. MRM61C (Kunjin) virus also stands closely with lineage I strains as reported previously (Lanciotti et al., 2002Down; Bakonyi et al., 2005Down). WN(S)V is grouped together with lineage II strains. In particular, WN(S)V is clustered together with WNFCG. We have also generated the phylogenetic tree from the complete E gene sequence of 27 WNV strains and the resulting tree showed that WN(S)V is closely related to WNFCG (data not shown). This relationship is statistically supported in the complete genome tree (Fig. 2bDown), where WN(S)V clustered with WNFCG (bootstrap value of 100 %). Phylogenetic trees were also generated based on the C, NS5 gene/3'-UTR sequences (Fig. 2c and dDown). In all cases, the tree demonstrated the same overall topology. The use of complete genome and other virus gene sequences eliminated the possibility of limited phylogenetic comparison. Other algorithms such as DNA parsimony and DNA maximum-likelihood were also tested and the resulting trees displayed nearly identical topology to the neighbour-joining analysis (data not shown).


Figure 2
View larger version (46K):
[in this window]
[in a new window]
 
Fig. 2. (a) Phylogenetic tree generated using neighbour-joining algorithm based on 227 nt sequence of the E gene from 37 WNV strains. Japanese encephalitis isolate NC_001437 [GenBank] JEVgp1 was used as the outgroup in the dendrogram. (b) Phylogenetic tree constructed from the nucleotide sequences of complete virus genomes. (c) Phylogenetic tree constructed from the nucleotide sequences of the C gene. (d) Phylogenetic tree constructed from the nucleotide sequences of NS5/3'-UTR. Values near the branches indicate the percentage bootstrap confidence levels. GenBank accession numbers for sequences used to construct the tree are indicated beside the abbreviated name of the virus strain. (e) Alignment of deduced amino acid sequences of the E proteins of HNY1999, WN(S)V and WNFCG (Wengler) strains. HNY1999 was chosen to represent lineage I WNV strains listed in Fig. 2(a)Up. The single letter code for the amino acids is used. Dashed lines indicate identity between sequences. The semi- and non-conserved residues are shown against light and dark grey colours, respectively. The three main domains within the E protein are delineated with reference to Bressanelli et al. (2004)Down. Vertical bars represent domain I, while diagonal bars represent domain II. Brick pattern indicates domain III and diamond pattern indicates the transmembrane domain. Grey shaded region displays the linker region between the domains.

 
To verify the distribution of these mutations across the E protein, we aligned the amino acid sequence of the E proteins of HNY1999, WN(S)V and WNFCG strains. As shown in Fig. 2(e)Up, several substitutions were identified; fourteen of them were conserved. In addition to this, few semi- and non-conserved substitutions were found. The major difference between WN(S)V and WNFCG was the existence of a putative N-linked glycosylation site located at aa 154 of the E protein. There were also two other non-conserved substitutions in WN(S)V compared with WNFCG. These were E to A (E208) and T to K (E332). Twelve positions in the E protein were found to differ in WN(S)V compared with HYN1999. These were E to D (E55), I to T (E126), R to W (E128), T to I (E129), L to Q (E131), V to I (E159), N to S (E199), T to A (E208), V to T (E232), L to A (E312), G to R (E313) and T to K (E332). The majority of the amino acid differences identified among the various strains reside in domain I and II of the E proteins.

Generation of recombinant viruses
To analyse the effect of the glycosylation site of the E protein on the virus phenotype, full-length cDNA clones of pWNS and pWNS-S154A were constructed. After transfection of pWNS and pWNS-S154A RNAs, the progeny viruses were collected from the supernatant. The titre of amplified pWNS virus in Vero cells at 24 h p.i. was 2–6x108 p.f.u. ml–1, which was similar to that of the parental virus. The titre of amplified pWNS-S154A virus was lower but still reached 2–5x107 p.f.u. ml–1 at 48 h p.i. The plaques formed by parental WN(S) and pWNS viruses ranged from 2·1 to 2·6 mm on day 3, while pWNS-S154A virus displayed smaller size plaques (1·5 to 1·7 mm). However, the plaque size of 1·9 to 2·3 mm could be achieved by day 5 for pWNS-S154A virus. The genetic stability of pWNS and pWNS-S154A viruses was monitored by direct sequencing of the full-length genome by RT-PCR. Sequencing of the virus stocks from passages 1 and 3 confirmed that there are no other additional mutations except the intentionally introduced mutations.

Western blot analysis
To examine the potential molecular size change caused by the silencing of E-NYST, the cell lysates from the Vero cells infected with WN(S)V, pWNS or pWNS-S154A viruses were subjected to Western blot analysis. As shown in Fig. 3Down(a), differences in electrophoretic mobility of the E proteins were observed, which indicated the presence (lanes 2 and 3) and absence (lane 4) of glycosylation of the E protein of the above-mentioned viruses. These results are consistent with previous reports showing the differences in the molecular mass of glycosylated and non-glycosyalated forms of the E protein (Johnson et al., 1994Down; Chambers et al., 1998Down; Lorenz et al., 2003Down; Shirato et al., 2004Down). Complementation analysis was performed by transfecting Vero cells with the RNA from pGEM-C-prM-E followed by infection with the pWNS-S154A virus. The expression of both glycosylated and non-glycosylated forms of the E proteins (lane 5) was observed. The non-glycosylated E protein, shown in lower intensity (lane 5), migrated faster than the glycosylated form.


Figure 3
View larger version (26K):
[in this window]
[in a new window]
 
Fig. 3. (a) Western blot analysis showing the glycosylation status of the WNV E proteins. Cell culture supernatants harvested at 24 h p.i. from mock-infected Vero cells (lane 1) and from Vero cells infected with parental WN(S)V (lane 2), pWNS (lane 3) or pWNS-S154A (lane 4) virus, were analysed by SDS-PAGE followed by Western blot analysis. The blot was probed with WNV-specific anti-E monoclonal antibody. The non-glycosylated E protein (lane 4) migrated faster than the glycosylated forms (lanes 2 and 3). Complementation analysis performed by transfecting the Vero cells with the RNA from pGEM-C-prM-E, 5 h prior to infection with the pWNS-S154A virus, showed both the glycosylated and non-glycosylated forms of the E proteins (lane 5). (b) Comparison of the growth kinetics of parental ({blacklozenge}), pWNS ({blacksquare}), pWNS-S154A ({blacktriangleup}) and pWNS-S154A plus complementation (bullet) viruses at an m.o.i. of 5·0 in Vero cells determined at 6, 12, 18 and 24 h after the inoculation. The mean values±SD are representative of three independent experiments. (c and d) Comparison of the growth kinetics for pWNS ({blacksquare}) and pWNS-S154A ({blacktriangleup}) infectious clones derived viruses by transfecting Vero cells with transcribed RNAs. The culture supernatant was sampled at selected time points (8, 16, 24, 32, 40, 48 and 56 h) and were used for titration of infectious virus by plaque assay (c) and for detection of viral RNA copies by real-time RT-PCR (d). The mean values±SD are representative of two independent experiments. The detection limit for real-time assay was >102 copies. (e and f)Comparison of the growth kinetics of parental ({blacklozenge}), pWNS ({blacksquare}) and pWNS-S154A ({blacktriangleup}) viruses at an m.o.i. of 5·0 in BHK-21 and C6/36 cell lines, respectively, determined at 6, 12, 18 and 24 h after the inoculation. The mean values±SD are representative of three independent experiments.

 
Growth kinetics of pWNS-S154A virus
In order to examine whether the efficiency of the virus replication was affected by the S154A mutation, the virus growth kinetics were determined. Fig. 3(b)Up shows the comparison of growth curves between parental WN(S), pWNS and pWNS-S154A viruses. The growth kinetics in Vero cells showed no significant difference between the parental WN(S) and pWNS viruses. However, the titre of the pWNS-S154A virus was approximately 100-fold lower than that of pWNS virus. In addition, the pWNS-S154A replicated relatively slower during the early stages after infection. The impaired replication was reversed by wild-type E protein complementation as shown in Fig. 3(b)Up, which indicated that the expression of the wild-type E protein is able to compensate for the propagation of the pWNS-S154A virus. These results proved that the glycosylation site of the E protein is important for efficient virus replication in mammalian cells.

The lower replication efficiency of the pWNS-S154A virus might be caused by the presence of a higher proportion of defective viruses. To test this possibility, in vitro transcribed RNAs from pWNS and pWNS-S154A plasmids were used for direct transfection. The supernatant from the transfected cell culture was collected for plaque assay. As shown in Fig. 3(c)Up, the pWNS-S154A virus similarly showed impaired propagations, with the titres being 100-fold lower than those of pWNS virus at 56 h post-transfection. Transfection of the infectious RNA also seemed to increase the latent period compared with infection with intact virus. It showed that the pWNS-S154A RNA is less efficient at replication than the RNA from the pWNS clone.

To investigate the possibility that the production of infectious RNA is less adequate for pWNS-S154A virus, the total number of viral RNA copies (positive-strand) was determined using real-time RT-PCR. As shown in Fig. 3(d)Up, the number of viral RNA copies produced from the pWNS-154A clone was less than that from the pWNS clone at 32 and 40 h post-transfection.

To demonstrate that the defect observed is not limited to the Vero cells, we also examined the growth kinetics of parental WN(S), pWNS and pWNS-S154A viruses in BHK-21 and C6/36 cells, which are shown in Fig. 3(e and f)Up, respectively. These results confirmed the impaired replication of pWNS-S154A virus in both mammalian and mosquito cell lines.

Ultrastructural analysis
To assess the possibility that N-linked glycosylation of the E protein also plays a role in the maturation process of the virus, TEM analysis was carried out. As shown in Fig. 4Down(a and b), virus particles egressed by budding at the plasma membrane (arrows) of Vero cells infected by parental WN(S) and pWNS viruses, respectively, which is a distinctive feature of this strain (Ng et al., 1994Down). Fig. 4(c)Down shows a large proportion of the cell cytoplasm (Cy). There was no accumulation of progeny virus particles intracellularly or at the perinuclear region. Virus particles (arrows) were observed to be budding from the plasma membrane similar to Fig. 4(b)Down.


Figure 4
View larger version (179K):
[in this window]
[in a new window]
 
Fig. 4. Localization of parental WN(S) and pWNS viruses in infected Vero cells using TEM. Vero cells were infected with parental and pWNS viruses at an m.o.i. of 5·0 before being processed for TEM (CM 120; BioTwin) at 14 h p.i. Both the parental (a) and pWNS (b and c) virus particles are shown budding out of the plasma membrane (arrows). (c) A larger cytoplasmic (Cy) area is shown to illustrate that there was no noticeable intracellular maturation of the pWNS progeny virus. Nu, Nucleus; RER, rough endoplasmic reticulum.

 
In the pWNS-S154A virus-infected cells, virus particles were formed within the swollen lumen of the endoplasmic reticulum (Fig. 5aDown, arrows). This is similar to the trans-mode of maturation seen in many other flaviviruses. The usual smooth membrane vesicles with thread-like structures (Fig. 5aDown, arrowheads) within swollen endoplasmic reticulum were observed in the infected cells. At higher magnification, mature progeny virus particles (Fig. 5bDown, arrows) were observed scattered among the smooth membrane vesicles (Fig. 5bDown, arrowheads). In Fig. 5(a and b)Down, no virus budding was observed from the plasma membrane (PM) after intensive screening of the cells.


Figure 5
View larger version (201K):
[in this window]
[in a new window]
 
Fig. 5. Localization of pWNS-S154A virus in infected Vero cells using TEM. (a and b) Maturation of pWNS-S154A progeny virus particles is shown in the swollen lumen of the endoplasmic reticulum (arrows) at 24 h p.i. This is similar to the trans-mode of maturation for many other flaviviruses. The virus particles are present among the smooth membrane vesicles with thread-like structures (arrowheads), a typical flavivirus-induced change. No virus particle is observed budding out of the plasma membrane (PM). (c and d) Complementation analysis was performed by transfecting the Vero cells with the wild-type E protein followed by infection with pWNS-S154A virus. The progeny virus particles are shown budding from the plasma membrane (arrows) at 14 h p.i. The earlier time is chosen as the virus in this case replicates faster than pWNS-S154A-infected cells. The inset in (d) shows a budding virus particle. There was minimal intracellular maturation observed (arrowheads). The images shown here are representatives of 25 cells for each virus infection.

 
When complementation analysis was performed by transfecting the Vero cells with the wild-type E protein followed by infection with pWNS-S154A virus, progeny virus particles were observed budding from the plasma membrane (Fig. 5c, dUp and inset, arrows) in the majority of the cells, even though the intracellular maturation (Fig. 5c and dUp, arrowheads) was also observed. This is similar to the parental or the glycosylated recombinant virus [(pWNS) (Fig. 4Up)]. From the observation in this study, it appears that when the virus matures at the plasma membrane, only a small quantity of paracrystals/convoluted membranes and smooth membrane vesicles are observed, unlike in intracellular maturation. TEM was performed at three different timings (14, 18 and 24 h p.i.) for all the viruses (parental, pWNS, pWNS-S154A viruses and complementation to pWNS-S154A viruses) and similar results were observed irrespective of the timings. These observations suggested that the carbohydrate side chain of the E protein is a determinant of the mode of virus maturation in WN(S)V.


   DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Various methods have been employed to address the genomic diversity of different strains of WNV. These mainly include phylogenetic analysis of aligned complete nucleotide sequences and selected gene sequences. Genetic studies have shown that WN viruses can be divided into two lineages. Recently, Rabensburg isolate was suggested to form a new lineage of WNV or a novel member of the JEV group.

Phylogenetic analysis performed previously placed WN(S)V under lineage II (Scherret et al., 2001aDown). However, this analysis was performed using partial nucleotide sequence of WN(S)V. To extend these findings and to determine more precisely the genetic relationships between the WN(S)V and other WNV isolates, the complete genomic sequence of WN(S)V was determined, and its phylogenetic relationship with other WNV strains was analysed utilizing their complete genomes, C, E and NS5/3'-UTR sequences. The outcome of these analyses suggested that WN(S)V was most closely related to WNFCG (Fig. 2Up). This is the first study showing the phylogenetic trees utilizing the complete genome, C, E and NS5/3'-UTR nucleotide sequences of WN(S)V for the analyses.

Several groups (Scherret et al., 2001aDown; Beasley & Barrett, 2002Down; Li et al., 2005bDown; Oliphant et al., 2005Down; Sanchez et al., 2005Down) have shown that amino acids E306 (Ser), E307 (Lys), E330 (Thr) and E332 (Thr) constituted the epitope on the surface of domain III that can be recognized by lineage-differentiating monoclonal antibodies, and the mutation at residue E332 affected neutralization and membrane receptor preparation-binding characteristics of the virus. Li et al. (2005aDown or bDown) demonstrated that T332K and L312A mutations were sufficient to disrupt the epitopes recognized by neutralizing antibodies. In Sarafend, the residues at E312 and E332 are Ala and Lys, respectively. Given the location of these residues in the identified neutralizing epitopes, WN(S)V could probably escape neutralization by the monoclonal antibodies raised against lineage I strain, HNY1999. This awaits further verification.

The mutations identified at residues E126 and E159 were also reported in two other lineage I strains, Italy 1998 and Romania 1996 (Lanciotti et al., 2002Down) and these mutations have been shown to affect the virus phenotypes in terms of virus neutralization or virulence in other flaviviruses (Lobigs et al., 1987Down; McMinn et al., 1995Down; Ryman et al., 1997Down). Berthet et al. (1997)Down and Scherret et al. (2001a)Down defined the signature motifs within the E gene that can differentiate lineage I from lineage II. When we analysed the presence of these signature motifs in WN(S)V, the residues at E128, E129, E131, E172, E199, E205, E208 and E210 concurred with their observations (Fig. 2eUp). The significance of other mutations identified within the E protein remains to be unveiled.

As demonstrated in this study, the glycosylation site of the WN(S)V E protein plays an important role in regulating the rate of virus replication. Abolishing the glycosylation site of the E protein delayed the replication rate at the early phase (Fig. 3bUp). Initially, plaque assay was used to determine the infectious virus particles from the supernatant harvested from Vero cells infected with the pWNS and pWNS-S154A viruses. It has been shown that the loss of E protein glycosylation affected the viral stability and infectivity in Dengue virus and St. Louis encephalitis virus (Pletnev et al., 1993Down; Vorndam et al., 1993Down; Ishak et al., 2001Down). Thus, the growth kinetics were also determined using in vitro transcribed RNAs (Fig. 3cUp). The difference in virus titres between pWNS and pWNS-S154A was observed in both the virus-infected and RNA-transfected samples (Fig. 3b and cUp). In addition, the total amount of viral RNA measured by real-time RT-PCR was higher for the pWNS virus than the pWNS-S154A virus (Fig. 3dUp). To gain an insight into how significant the defect (impaired replication) observed in Vero cells was, we extended our studies to BHK-21 and C6/36 cell lines using the intact viruses (Fig. 3e and fUp). Taken together, these results proved that the glycosylation site of the E protein is important for efficient virus replication in both mammalian and mosquito cell lines.

The exact mechanism behind the impaired virus replication is unknown. Entry of the virus into the cells could possibly be affected by the mutation abolishing the glycosylation site, since the E protein plays a critical role in receptor attachment. The loss of glycosylation could also affect the assembly and secretion of the virus particles (Lorenz et al., 2003Down). Pulse–chase analysis will be performed in the future to demonstrate the defect in the secretion of the non-glycosylated E protein. Another possible mechanism behind reduced virus replication could be decreased virus–E protein interactions mediated by the glycan residues. Rey et al. (1995)Down suggested that the glycan residue may interact with domain II of the other E subunit and stabilize the dimeric structure. Lorenz et al. (2003)Down also linked the function of the carbohydrate moiety to the other glycoprotein, prM, and suggested that this glycan could be necessary in forming the appropriate interface between the E glycoproteins. The glycosylated amino acids were mapped onto the E dimer near domain III (Rey et al., 1995Down; Modis et al., 2003Down; Mukhopadhyay et al., 2003Down; Zhang et al., 2003Down; Zhang et al., 2004Down). Moreover, it has been suggested that the carbohydrate moieties on the virus surface might regulate the specificity of virus binding to the cellular receptors (Navarro-Sanchez et al., 2003Down).

Previously, our study indicated that 5'-UTR and the three structural genes played a role in influencing the mode of maturation for WN(S)V (Li et al., 2005aDown or bDown). This concurred with the current observation that the non-glycosylated virus (pWNS-S154A) showed an altered maturation mode. Parental WN(S) and pWNS viruses were observed budding out of the cell plasma membrane (Fig. 4Up), a characteristic feature of WN(S)V (Ng et al., 1994Down; Li et al., 2005aDown). In contrast, the matured pWNS-S154A viruses were observed in the swollen lumen of the endoplasmic reticulum (Fig. 5a and bUp) similar to the majority of flaviviruses (Hase et al., 1987aDown, bDown; Ng, 1987Down; Hase, 1993Down; Serafino et al., 1997Down; Rahman et al., 1998Down; Lanciotti et al., 2002Down). This result indicated that maturation by budding at the plasma membrane was altered in pWNS-S154A recombinant virus.

A more detailed study will be conducted using immuno-gold labelling to trace the transportation of the synthesized E protein both in the pWNS-S154A-infected cells and in infected cells complemented with wild-type E proteins. These experiments can confirm further whether E proteins synthesized in pWNS-S154A-infected cells are only localized with the virus-induced membranes and not transported to the plasma membrane as in wild-type virus infection.

It was shown previously that the assembly of WN(S)V was affected by brefeldin A (BFA) treatment. Egression at the plasma membrane and filopodia formation, characteristic features of WN(S)V, were not observed in the presence of BFA (Sreenivasan et al., 1993Down). BFA is a fungal metabolite known to affect glycoprotein maturation (Whealy et al., 1991Down). Taken together with the previous reports and the current findings, it is reasonable to conclude that the E protein glycosylation (status) is the major molecular determinant of the budding mode at the plasma membrane for WN(S)V.

To our knowledge, there are no data to show whether the New York (NY) strain of WNV matures intracellularly or at the plasma membrane. It is therefore tempting to postulate that, if WN(NY)V also buds at the plasma membrane like WN(S)V, the same glycosylation site might influence the maturation mode. In addition, the replication efficiency of WN(NY)V was also slowed down by the mutation at the glycosylation site of the E protein (Shirato et al., 2004Down), consistent with our observations in this study.

There have been several studies showing the relationship between N-linked glycosylation in the E protein and neuroinvasiveness of several virus isolates in vivo. The attenuation of neuroinvasiveness after mutation at the glycosylation site of the E protein was proven in NY strains (NY99-6922, BC787 and TM171-03; Beasley et al., 2002Down; Shirato et al., 2004Down). However, the loss of glycosylation site did not affect the neurovirulence. Alternatively, Chambers et al. (1998)Down suggested that glycosylation might not be directly responsible for the attenuation of WNV since neuroinvasiveness was also observed with the non-glycosylated E protein variant of the Israel strain. Wang et al. (2003)Down showed that WN(S)V is neuroinvasive and neurovirulent in murine models and that virus replication in peripheral tissues is required for virus invasion of the central nervous system. This observation is consistent with the study on NY strains of WNV by Shirato et al. (2004)Down. In this study, we have not determined whether this mutation of the E protein might affect the neuroinvasiveness of WN(S)V. This needs further verification in the animal model. However, it is rational to conclude that the impaired replication rate in pWNS-S154A virus resulted from the mutation of the E protein glycosylation site.

In summary, a comprehensive phylogenetic analysis revealed that WN(S)V is closely related to WNFCG, which is a lineage II virus. The presence of the glycosylation site in the E protein of WNV isolates is advantageous for the virus as it is shown here that loss of glycosylation affects the rate of replication. Glycosylation is also shown to play an important role in influencing the maturation mode of the WN(S)V.


   ACKNOWLEDGEMENTS
 
This work was supported by grants from the Biomedical Research Council, Singapore (BMRC/01/1/21/18/003) and the National University of Singapore (R-182-000-055-112). We thank Dr Chan Shih Yen (Singapore General Hospital) and Kunde Ramamoorthy Govindarajan (Genomic Institute of Singapore) for their useful advice on the phylogenetic analyses. We also thank Terence Tan for his support throughout this study.


   REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Adams, S. C., Broom, A. K., Sammels, L. M., Hartnett, A. C., Howard, M. J., Coelen, R. J., Mackenzie, J. S. & Hall, R. A. (1995). Glycosylation and antigenic variation among Kunjin virus isolates. Virology 206, 49–56.[CrossRef][Medline]

Bakonyi, T., Hubalek, Z., Rudolf, I. & Nowotny, N. (2005). Novel flavivirus or new lineage of West Nile virus, central Europe. Emerg Infect Dis 11, 225–231.[Medline]

Beasley, D. W. & Barrett, A. D. (2002). Identification of neutralizing epitopes within structural domain III of the West Nile virus envelope protein. J Virol 76, 13097–13100.[Abstract/Free Full Text]

Beasley, D. W., Li, L., Suderman, M. T. & Barrett, A. D. (2002). Mouse neuroinvasive phenotype of West Nile virus strains varies depending upon virus genotype. Virology 296, 17–23.[CrossRef][Medline]

Berthet, F. X., Zeller, H. G., Drouet, M. T., Rauzier, J., Digoutte, J. P. & Deubel, V. (1997). Extensive nucleotide changes and deletions within the envelope glycoprotein gene of Euro-African West Nile viruses. J Gen Virol 78, 2293–2297.[Abstract]

Bhuvanakantham, R. & Ng, M. L. (2005). Analysis of self-association of West Nile virus capsid protein and the crucial role played by Trp 69 in homodimerization. Biochem Biophys Res Commun 329, 246–255.[Medline]

Brandt, M., Yao, K., Liu, M., Heckert, R. A. & Vakharia, V. N. (2001). Molecular determinants of virulence, cell tropism, and pathogenic phenotype of infectious bursal disease virus. J Virol 75, 11974–11982.[Abstract/Free Full Text]

Bressanelli, S., Stiasny, K., Allison, S. L., Stura, E. A., Duquerroy, S., Lescar, J., Heinz, F. X. & Rey, F. A. (2004). Structure of a flavivirus envelope glycoprotein in its low-pH-induced membrane fusion conformation. EMBO J 23, 728–738.[CrossRef][Medline]

Chambers, T. J., Hahn, C. S., Galler, R. & Rice, C. M. (1990). Flavivirus genome organization, expression, and replication. Annu Rev Microbiol 44, 649–688.[CrossRef][Medline]

Chambers, T. J., Halevy, M., Nestorowicz, A., Rice, C. M. & Lustig, S. (1998). West Nile virus envelope proteins: nucleotide sequence analysis of strains differing in mouse neuroinvasiveness. J Gen Virol 79, 2375–2380.[Abstract]

Halevy, M., Akov, Y., Ben-Nathan, D., Kobiler, D., Lachmi, B. & Lustig, S. (1994). Loss of active neuroinvasiveness in attenuated strains of West Nile virus: pathogenicity in immunocompetent and SCID mice. Arch Virol 137, 355–370.[CrossRef][Medline]

Hase, T. (1993). Virus-neuron interactions in the mouse brain infected with Japanese encephalitis virus. Virchows Arch B Cell Pathol Incl Mol Pathol 64, 161–170.[Medline]

Hase, T., Summers, P. L., Eckels, K. H. & Baze, W. B. (1987a). Maturation process of Japanese encephalitis virus in cultured mosquito cells in vivo and mouse brain cells in vivo. Arch Virol 96, 135–151.[CrossRef][Medline]

Hase, T., Summers, P. L., Eckels, K. H. & Baze, W. B. (1987b). An electron and immunoelectron microscopic study of dengue-2 virus infection of cultured mosquito cells: maturation events. Arch Virol 92, 273–291.[CrossRef][Medline]

Heinz, F. X. (1986). Epitope mapping of flavivirus glycoproteins. Adv Virus Res 31, 103–168.[Medline]

Ishak, H., Takegami, T., Kamimura, K. & Funada, H. (2001). Comparative sequences of two type 1 dengue virus strains possessing different growth characteristics in vivo. Microbiol Immunol 45, 327–331.[Medline]

Johnson, A. J., Guirakhoo, F. & Roehrig, J. T. (1994). The envelope glycoproteins of dengue 1 and dengue 2 viruses grown in mosquito cells differ in their utilization of potential glycosylation sites. Virology 203, 241–249.[CrossRef][Medline]

Lad, V. J., Shende, V. R., Gupta, A. K., Koshy, A. A. & Roy, A. (2000). Effect of tunicamycin on expression of epitopes on Japanese encephalitis virus glycoprotein E in porcine kidney cells. Acta Virol 44, 359–364.[Medline]

Lanciotti, R. S., Ebel, G. D., Deubel, V. & 9 other authors (2002). Complete genome sequences and phylogenetic analysis of West Nile virus strains isolated from the United States, Europe, and the Middle East. Virology 298, 96–105.[CrossRef][Medline]

Li, J., Bhuvanakantham, R., Howe, J. & Ng, M. L. (2005a). Identifying the region influencing the cis-mode of maturation of West Nile (Sarafend) virus using chimeric infectious clones. Biochem Biophys Res Commun 334, 714–720.[CrossRef][Medline]

Li, L., Barrett, A. D. & Beasley, D. W. (2005b). Differential expression of domain III neutralizing epitopes on the envelope proteins of West Nile virus strains. Virology 335, 99–105.[CrossRef][Medline]

Lobigs, M., Dalgarno, L., Schlesinger, J. J. & Weir, R. C. (1987). Location of a neutralization determinant in the E protein of yellow fever virus (17D vaccine strain). Virology 161, 474–478.[CrossRef][Medline]

Lorenz, I. C., Kartenbeck, J., Mezzacasa, A., Allison, S. L., Heinz, F. X. & Helenius, A. (2003). Intracellular assembly and secretion of recombinant subviral particles from tick-borne encephalitis virus. J Virol 77, 4370–4382.[Abstract/Free Full Text]

McMinn, P. C., Lee, E., Hartley, S., Roehrig, J. T., Dalgarno, L. & Weir, R. C. (1995). Murray Valley encephalitis virus envelope protein antigenic variants with altered hemagglutination properties and reduced neuroinvasiveness in mice. Virology 211, 10–20.[CrossRef][Medline]

Modis, Y., Ogata, S., Clements, D. & Harrison, S. C. (2003). A ligand-binding pocket in the dengue virus envelope glycoprotein. Proc Natl Acad Sci U S A 100, 6986–6991.[Abstract/Free Full Text]

Mukhopadhyay, S., Kim, B. S., Chipman, P. R., Rossmann, M. G. & Kuhn, R. J. (2003). Structure of West Nile virus. Science 302, 248.[Free Full Text]

Navarro-Sanchez, E., Altmeyer, R., Amara, A., Schwartz, O., Fieschi, F., Virelizier, J. L., Arenzana-Seisdedos, F. & Despres, P. (2003). Dendritic-cell-specific ICAM3-grabbing non-integrin is essential for the productive infection of human dendritic cells by mosquito-cell-derived dengue viruses. EMBO Rep 4, 723–728.[CrossRef][Medline]

Ng, M. L. (1987). Ultrastructural studies of Kunjin virus-infected Aedes albopictus cells. J Gen Virol 68, 577–582.[Abstract/Free Full Text]

Ng, M. L., Howe, J., Sreenivasan, V. & Mulders, J. J. (1994). Flavivirus West Nile (Sarafend) egress at the plasma membrane. Arch Virol 137, 303–313.[CrossRef][Medline]

Oliphant, T., Engle, M., Nybakken, G. E. & 11 other authors (2005). Development of a humanized monoclonal antibody with therapeutic potential against West Nile virus. Nat Med 11, 522–530.[CrossRef][Medline]

Pletnev, A. G., Bray, M. & Lai, C. J. (1993). Chimeric tick-borne encephalitis and dengue type 4 viruses: effects of mutations on neurovirulence in mice. J Virol 67, 4956–4963.[Abstract/Free Full Text]

Rahman, S., Matsumura, T., Masuda, K., Kanemura, K. & Fukunaga, T. (1998). Maturation site of dengue type 2 virus in cultured mosquito C6/36 cells and Vero cells. Kobe J Med Sci 44, 65–79.[Medline]

Rey, F. A., Heinz, F. X., Mandl, C., Kunz, C. & Harrison, S. C. (1995). The envelope glycoprotein from tick-borne encephalitis virus at 2 Å resolution. Nature 375, 291–298.[CrossRef][Medline]

Ryman, K. D., Ledger, T. N., Weir, R. C., Schlesinger, J. J. & Barrett, A. D. (1997). Yellow fever virus envelope protein has two discrete type-specific neutralizing epitopes. J Gen Virol 78, 1353–1356.[Abstract]

Sanchez, M. D., Pierson, T. C., McAllister, D., Hanna, S. L., Puffer, B. A., Valentine, L. E., Murtadha, M. M., Hoxie, J. A. & Doms, R. W. (2005). Characterization of neutralizing antibodies to West Nile virus. Virology 336, 70–82.[CrossRef][Medline]

Scherret, J. H., Poidinger, M., Mackenzie, J. S., Broom, A. K., Deubel, V., Lipkin, W. I., Briese, T., Gould, E. A. & Hall, R. A. (2001a). The relationships between West Nile and Kunjin viruses. Emerg Infect Dis 7, 697–705.[Medline]

Scherret, J. H., Mackenzie, J. S., Khromykh, A. A. & Hall, R. A. (2001b). Biological significance of glycosylation of the envelope protein of Kunjin virus. Ann N Y Acad Sci 951, 361–363.[Free Full Text]

Scherret, J. H., Mackenzie, J. S., Hall, R. A., Deubel, V. & Gould, E. A. (2002). Phylogeny and molecular epidemiology of West Nile and Kunjin viruses. Curr Top Microbiol Immunol 267, 373–390.[Medline]

Serafino, A., Valli, M. B., Alessandrini, A., Ponzetto, A., Carloni, G. & Bertolini, L. (1997). Ultrastructural observations of viral particles within hepatitis C virus-infected human B lymphoblastoid cell line. Res Virol 148, 153–159.[CrossRef][Medline]

Shi, P. Y., Tilgner, M., Lo, M. K., Kent, K. A. & Bernard, K. A. (2002). Infectious cDNA clone of the epidemic West Nile virus from New York City. J Virol 76, 5847–5856.[Abstract/Free Full Text]

Shirato, K., Miyoshi, H., Goto, A., Ako, Y., Ueki, T., Kariwa, H. & Takashima, I. (2004). Viral envelope protein glycosylation is a molecular determinant of the neuroinvasiveness of the New York strain of West Nile virus. J Gen Virol 85, 3637–3645.[Abstract/Free Full Text]

Sreenivasan, V., Ng, K. L. & Ng, M. L. (1993). Brefeldin A affects West Nile virus replication in Vero cells but not C6/36 cells. J Virol Methods 45, 1–17.[CrossRef][Medline]

Vorndam, V., Mathews, J. H., Barrett, A. D., Roehrig, J. T. & Trent, D. W. (1993). Molecular and biological characterization of a non-glycosylated isolate of St Louis encephalitis virus. J Gen Virol 74, 2653–2660.[Abstract/Free Full Text]

Wang, Y., Lobigs, M., Lee, E. & Mullbacher, A. (2003). CD8+ T cells mediate recovery and immunopathology in West Nile virus encephalitis. J Virol 77, 13323–13334.[Abstract/Free Full Text]

Westaway, E. G., Brinton, M. A., Gaidamovich, S. & 7 other authors (1985). Flaviviridae. Intervirology 24, 183–192.[Medline]

Whealy, M. E., Card, J. P., Meade, R. P., Robbins, A. K. & Enquist, L. W. (1991). Effect of brefeldin A on alphaherpesvirus membrane protein glycosylation and virus egress. J Virol 65, 1066–1081.[Abstract/Free Full Text]

Zhang, W., Chipman, P. R., Corver, J. & 7 other authors (2003). Visualization of membrane protein domains by cryo-electron microscopy of dengue virus. Nat Struct Biol 10, 907–912.[CrossRef][Medline]

Zhang, Y., Zhang, W., Ogata, S., Clements, D., Strauss, J. H., Baker, T. S., Kuhn, R. J. & Rossmann, M. G. (2004). Conformational changes of the flavivirus E glycoprotein. Structure 12, 1607–1618.[Medline]

Received 6 July 2005; accepted 6 November 2005.


This article has been cited by other articles:


Home page
J. Gen. Virol.Home page
D. B. A. Medeiros, M. R. T. Nunes, P. F. C. Vasconcelos, G.-J. J. Chang, and G. Kuno
Complete genome characterization of Rocio virus (Flavivirus: Flaviviridae), a Brazilian flavivirus isolated from a fatal case of encephalitis during an epidemic in Sao Paulo state
J. Gen. Virol., August 1, 2007; 88(8): 2237 - 2246.
[Abstract] [Full Text] [PDF]


Home page
J. Virol.Home page
M. A. Samuel and M. S. Diamond
Pathogenesis of West Nile Virus Infection: a Balance between Virulence, Innate and Adaptive Immunity, and Viral Evasion
J. Virol., October 1, 2006; 80(19): 9349 - 9360.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Li, J.
Right arrow Articles by Ng, M.-L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Li, J.
Right arrow Articles by Ng, M.-L.
Agricola
Right arrow Articles by Li, J.
Right arrow Articles by Ng, M.-L.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS