J Gen Virol Email Content Delivery
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Gen Virol 87 (2006), 895-907; DOI 10.1099/vir.0.81236-0

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary tables
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 CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by von Lindern, J. J.
Right arrow Articles by Barrett, A. D. T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by von Lindern, J. J.
Right arrow Articles by Barrett, A. D. T.
Agricola
Right arrow Articles by von Lindern, J. J.
Right arrow Articles by Barrett, A. D. T.
© 2006 Society for General Microbiology

Genome analysis and phylogenetic relationships between east, central and west African isolates of Yellow fever virus

Jana J. von Lindern1, Sarah Aroner2, Nicholas D. Barrett3, Jason A. Wicker4, C. Todd Davis1 and Alan D. T. Barrett1,4,5

1 Department of Pathology, University of Texas Medical Branch at Galveston, 301 University Blvd, Galveston, TX 77555-0609, USA
2 University of California at Berkeley, Berkeley, CA 94720, USA
3 School of Natural Sciences, University of Texas, Austin, TX 78712, USA
4 Department of Microbiology and Immunology, University of Texas Medical Branch at Galveston, 301 University Blvd, Galveston, TX 77555-0609, USA
5 Sealy Center for Vaccine Development and Center for Biodefense and Emerging Infectious Diseases, University of Texas Medical Branch at Galveston, 301 University Blvd, Galveston, TX 77555-0609, USA

Correspondence
Alan D. T. Barrett
abarrett{at}utmb.edu


   ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Yellow fever virus (YFV), a reemerging disease agent in Africa and South America, is the prototype member of the genus Flavivirus. Based on examination of the prM/M, E and 3' non-coding regions of the YFV genome, previous studies have identified seven genotypes of YFV, including the Angolan, east/central African and east African genotypes, which are highly divergent from the prototype strain Asibi. In this study, full genome analysis was used to expand upon these genetic relationships as well as on the very limited full genome database for YFV. This study was the first to investigate genomic sequences of YFV strains from east and central Africa (Angola71, Uganda48a and Ethiopia61b). All three viruses had genomes of 10 823 nt in length. Compared with the prototype strain Asibi (from west Africa) they were approximately 25 % divergent in nucleotide sequence and 7 % divergent in amino acid sequence. Comparison of multiple flaviviruses in the N-terminal region of NS4B showed that amino acid sequences were variable and that west African strains of YFV had an amino acid deletion at residue 21. Additionally, N-linked glycosylation sites were conserved between viral genotypes, while codon usage varied between strains.

Supplementary tables are available in JGV Online.


   INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Yellow fever virus (YFV) is a reemerging mosquito-borne virus of major public health importance in Africa and South America, whose incidence has been increasing over the past 25 years (Robertson et al., 1996Down; Monath, 1999Down). Infection of humans with YFV may result in a haemorrhagic fever with high mortality (Tomori, 1999Down). Transmission cycles involve both a jungle/sylvatic cycle and an urban cycle, with an additional intermediate ‘savannah’ cycle in Africa.

YFV is the prototype member of the family Flaviviridae, the genus Flavivirus, consisting of a single-stranded, positive-sense RNA genome, almost 11 kb in length. A type I 5'-terminal cap (m7GpppAmp) precedes the genome. The genome is flanked by 5' and 3' non-coding regions (NCR) and translated as a single open reading frame (ORF) that is co- and post-translationally cleaved into 10 proteins, including three structural proteins (C, prM/M and E) and seven non-structural proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B and NS5) (Barrett & Monath, 2003Down).

Understanding genetic relationships of this medically important viral pathogen is valuable for monitoring of future YFV outbreaks. Based on the examination of prM/M, E and 3'NCR of the YFV genome, previous studies have identified seven genotypes of YFV (Mutebi et al., 2001Down; Wang et al., 1996Down). Five African genotypes have been proposed, which include west African genotype I (strains from Nigeria and Ivory Coast), west African genotype II (strains from Senegal, plus additional Nigerian and Ivory Coast strains), east African genotype (strains from Kenya and Uganda), east/central African genotype [strains from Zaire (Democratic Republic of Congo), Uganda, Ethiopia and central African Republic] and Angolan genotype (Mutebi et al., 2001Down). Two South American genotypes have also been proposed: South American genotype I (Brazil, Columbia, Panama and Trinidad) and South American genotype II (predominantly Peru, plus a few isolates from Brazil and Trinidad) (Wang et al., 1996Down; Bryant & Barrett, 2003Down; Barrett & Monath, 2003Down).

In the present study, we sought to expand upon these genetic relationships using full genome analysis of east African (Uganda48a), east/central African (Ethiopia61b) and Angolan (Angola71) genotypes. Strain Angola71 was isolated from a human in early 1971, during a YF outbreak that occurred in this region during this time period (Pinto & Filipe, 1973Down). Prior to this outbreak, there had been no reported outbreaks in this country since 1872. Factors that led to the emergence of a YF outbreak in a region that had previously been free of the disease for nearly a century are unknown; however, genomic analysis of this particular virus strain may reveal potential contributing molecular characteristics. In addition, YFV strain Ethiopia61b was isolated during a human epidemic that occurred in this region in 1961 and 1962 (Andral et al., 1968Down; Serie et al., 1968Down), while Uganda48a was isolated in Uganda in 1948 from an unknown source, likely from a sylvan outbreak in this region at the time (Smithburn et al., 1949Down). Using these strains, we sought to expand upon the very limited full genome database for YFV, of which this report is the first to present genomic sequences of strains from east and central Africa.


   METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Viruses.
Origins, year of collection and references for viruses are shown in Table 1Down. Three virus genomes were sequenced in this study: the east African strain Uganda48a, the east/central African strain Ethiopia61b and the Angolan strain Angola71. These strains were obtained as lyophilized stocks from the World Reference Center for Emerging Viruses and Arboviruses at the University of Texas Medical Branch at Galveston. Viruses were reconstituted and passaged once in Vero cells to generate seed stocks, and were then given one additional passage in Vero cells to produce working stocks.


View this table:
[in this window]
[in a new window]
 
Table 1. Viruses

 
RNA extraction, RT-PCR and sequencing.
Methods for growth, extraction, amplification and sequencing of viral RNA have been described previously (Bryant & Barrett, 2003Down). Briefly, viral RNA was extracted using the Qiagen viral RNA extraction kit according to the manufacturer's instructions, followed by Superscriptase II or Titan RT-PCR protocols. Amplified fragments were visualized by ethidium bromide agarose gel staining, extracted from the gels and either sequenced directly if sufficient quantities of the DNA were recovered or cloned into the pGEM T-easy vector (Promega). For cloned PCR products, three separate clones were sequenced to ensure that the sequence was representative of the consensus sequence. YFV primers used in this study are listed in Supplementary Table S1 (available in JGV Online). Oligonucleotide primers were obtained from Sigma-Genosys.

Phylogenetic and sequence analyses.
Sequence analysis was performed using Contig Express. Translations and alignments were performed using Vector NTI. Percentage similarities and differences were calculated using MegAlign (DNASTAR). NetNGlyc (Gupta et al., 2004; http://www.cbs.dtu.dk/services/NetNGlyc/) was utilized for glycosylation site prediction. PAUP (Swofford, 2003Down) was utilized for the construction of phylogenetic trees by the neighbour-joining method (HKY85 parameter). In addition, GCUA version 1.1 (McInerney, 1998Down) was used for codon usage analysis and percentage GC calculations. For genomic codon usage analysis, the first 10 000 nt of the ORF was utilized (GCUA 1.1 is limited to 10 000 nt).


   RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Nucleotide sequence of the genomes of strains Uganda48a, Ethiopia61b and Angola71
To date, seven genotypes of YFV have been identified (Wang et al., 1996Down; Mutebi et al., 2001Down; Bryant et al., 2003Down), including five in Africa: west Africa I, west Africa II, Angola, east Africa and east/central Africa genotypes. Genomic sequences are available for the west Africa I genotype (Ivory Coast strain 85-82H; Pisano et al., 1997Down) and the west Africa II genotype: Asibi (Hahn et al., 1987Down), Gambia01 (Colebunders et al., 2002Down) and French viscerotropic virus (FVV) (Wang et al., 1995Down), while there are no genomic sequences available for the other three African genotypes. We determined the nucleotide sequence of the genomes of Uganda48a (east Africa genotype), Ethiopia61b (east/central Africa genotype) and Angola71 (the only known member of the Angolan genotype). All three viruses had genomes of 10 823 nt in length, with two copies of the repeat nucleotide sequence elements (RYF) in the 3'NCR that has been shown previously to characterize strains of this lineage (Mutebi et al., 2004Down).

Nucleotide sequence variation among diverse YFV strains
Uganda48a, Ethiopia61b and Angola71 were compared to previously published African YFV genomic sequences. In the neighbour-joining phylogenetic tree shown in Fig. 1Down, Uganda48a and Ethiopia61b showed a high level of identity with each other but were both distinct from Angola71. Together, these newly sequenced strains formed a separate clade from the west African strains. In addition, phylogenetic relationships similar to those shown in Fig. 1Down were observed when phylogenetic trees were constructed for each structural and non-structural protein gene independently, as well as the 3'NCR (data not shown).


Figure 1
View larger version (16K):
[in this window]
[in a new window]
 
Fig. 1. YFV genome phylogenetic tree. The tree was constructed by the neighbour-joining method using PAUP, with branch lengths drawn to scale. Yokose virus was used as an outgroup. Numbers at nodes are the scaled branch lengths, with bootstrap values of 1000 pseudoreplicates shown in parentheses.

 
Genomic sequence variation between the various strains ranged from 0·2 to 24·6 % compared with the prototype west African strain Asibi (Table 2Down). All three newly sequenced strains were approximately 25 % divergent in nucleotide sequence and 7 % divergent in amino acid sequence from the prototype Asibi strain (Table 2Down). Consistent with their classification as separate genotypes, Uganda48a and Angola71 were 17·1 % divergent in nucleotide sequence and 3·3 % divergent in amino acid sequence from each another. Interestingly, sequence divergence between Uganda48a and Ethiopia61b was much less pronounced, with these strains differing by 7·3 % at the nucleotide level and 1·5 % in amino acid sequence. The west Africa I genotype was 9 % divergent in nucleotide composition from the west Africa II genotype, 24·6 % divergent from the Angolan and east/central African genotypes and 24·2 % divergent from the east African genotype.


View this table:
[in this window]
[in a new window]
 
Table 2. Percent divergence between YFV genotypes

Strains: 1, Angola71 (Angola genotpe); 2, Asibi (west Africa II); 3, Ethiopia61b (east/central Africa); 4, FVV (west Africa II); 5, Gambia01 (west Africa II); 6, Ivory Coast (west Africa I); 7, Uganda48a (east Africa). Percentage divergence (nucleotide) is shown at the top of the table and percentage divergence (amino acid) is shown at the bottom of the table.

 
YFV genotypes were defined as viruses that differed by more than 7 % in genomic nucleotide composition, and the five African genotypes, as defined by Mutebi et al. (2001)Down, were supported by this study. In addition, various nucleotides differentiated the genotypes identified by Mutebi et al. (2001)Down. Lineage and genotype-specific nucleotides for the NS2A gene, a highly variable area of the genome, are shown in Supplementary Table S2 (available in JGV Online). Compared to the Asibi prototype, Angola71, Uganda48a and Ethiopia61b strains possessed 59 common nucleotide substitutions in NS2A. In addition, Uganda48a possessed 15 genotype-specific nucleotide substitutions in this region, while Ethiopia61b had 19 genotype-specific nucleotide substitutions and Angola71 had 35 genotype-specific nucleotide substitutions. West Africa I genotype Ivory Coast possessed 24 genotype-specific nucleotides in NS2A that differentiated it from west Africa II genotype strains.

Amino acid sequence variation between various YFV genotypes
Differences in amino acid composition ranged from 0·2 to 7·6 % for the various genotype-representative YFV strains compared with the prototype Asibi strain (Table 2Up). Alignment of the ORF of the various strains is shown in Fig. 2Down. Several amino acid substitutions characterized strains from west Africa, east Africa and Angola. Angola71, Uganda48a and Ethiopia61b shared many common amino acid substitutions and together comprised the east and central African lineage. The Angolan genotype possessed 52 additional genotype-specific amino acid substitutions (six in the C protein, one in prM/M, six in NS1, seven in NS2A, two in NS2B, seven in NS3, three in NS4A, two in NS4B and 18 in NS5), while the east African genotype-representative strain Uganda48a had 21 genotype-specific substitutions (four in the C protein, one in prM/M, one in E, four in NS1, two in NS2A, two in NS4B and seven in NS5) and Ethiopia61b had 23 (two in the C protein, three in prM/M, one in E, two in NS1, four in NS2A, four in NS3, one in NS4A and six in NS5). The west Africa I-representative strain, Ivory Coast, had 37 aa substitutions that differentiated it from west Africa II genotype strains (three in C, six in E, two in NS1, two in NS2A, one in NS2B, four in NS3, one in NS4A, four in NS4B and 14 in NS5; Supplementary Table S3 available in JGV Online).


Figure 2
Figure 2
Figure 2
Figure 2
View larger version (197K):
[in this window]
[in a new window]
 
Fig. 2. Alignment of YFV ORFs. Dots indicate consensus. Dashes indicate the absence of an amino acid at that position. First amino acids of individual viral proteins are indicated by arrows.

 
The N terminus of NS4B was highly variable between strains from eastern or central Africa and those from western Africa. Of particular note was a deletion of Thr at residue 21 of NS4B in west African strains of YFV, which was revealed upon alignment of the first 35 aa of NS4B of west African strains with Uganda48a, Ethiopia61b, Angola71 and the previously published NS4B sequences of strains from the central African Republic (CAR77-883 and CAR77-900) (Wang et al., 1996Down) (Fig. 3Downa). Since this region of NS4B is highly variable between YFV strains, we also compared these strains to other available flavivirus NS4B sequences in the alignment shown in Fig. 3(b)Down. As Yokose virus (a flavivirus isolated from a bat in Japan) is a close relative to YFV (Fig. 4Down), it is notable that this virus also possesses a cluster of amino acids in NS4B similar to the east and central African strains of YFV, but not west African strains (Fig. 3bDown). In addition to this variable region in NS4B, there was also a highly variable cluster of amino acids in the C protein of Uganda48a, Ethiopia61b and Angola71 (residues 103–110) and several regions of variability in NS5 of these strains (residues 287–295, 524–527, 556–562 and 632–660) compared with west African strains.


Figure 3
View larger version (60K):
[in this window]
[in a new window]
 
Fig. 3. (a) Alignment of the first 35 aa of NS4B from various YFV strains. Strains from east Africa and Angola (Uganda48a, CAR77-883, CAR77-900 and Angola71) are compared to strains from west Africa (Asibi, Gambia01 and Ivory Coast). Dashes indicate the absence of an amino acid at that position, while dots indicate identity with the consensus sequence. (b)Alignment of NS4B of various flaviviruses. Alignment was performed using MegAlign. IC, Ivory Coast strain of YFV; YOK, Yokose; WN, West Nile; KUN, Kunjin; SLE, St. Louis encephalitis; MVE, Murray valley encephalitis; USU, Usutu; JE, Japanese encephalitis; DEN, Dengue; LGT, Langat; OHF, Omsk haemorrhagic fever; TBE, Tick-borne encephalitis; ALK, Alkhurma; KAR, Karshi; Apoi, Apoi virus; LI, Louping Ill; MLL, Montana myotis leukoencephalitis; MOD, Modoc; KRV, Kamiti river virus; CFA, cell fusing agent.

 

Figure 4
View larger version (10K):
[in this window]
[in a new window]
 
Fig. 4. Nc plot comparing individual genes and ORFs of representative west African (Asibi) and central African (Ethiopia61b) YFV. The solid curve indicates expected Nc values if codon bias is due solely to base compositional constraints. See Table 4Up for the values used to construct the plot, included for purposes of distinguishing individual points.

 
YFV genotypes exhibit differences in codon usage
Codon usage of different genotypes was evaluated to assess relative synonymous codon usage (RSCU) values, defined as the number of times a particular codon is used, divided by the expected frequency of codon usage in the absence of a bias (i.e. RSCU=1). Upon examination of individual viral genes, RSCU values were determined and found to differ for the various genotypes. Values for the representative viral gene NS2A are shown in Table 3Down. NS2A possessed high levels of intergenotype variation in codon usage, with 19 aa exhibiting differences in codon usage. In addition to differences observed when individual viral genes were assessed, examination of genomic ORF (nt 1–10 000) revealed more similar RSCU values between the various genotypes, as well as RSCU values closer to 1 for several of the amino acids when compared to values for single genes (e.g. Asp, Asn, Gln, Phe, Lys and Cys) (Table 3Down).


View this table:
[in this window]
[in a new window]
 
Table 3. Codon usage

IC, Ivory Coast.

 
In addition, percentage G+C content was near 50 % for all genes examined (data not shown). In order to assess the potential relationship between codon usage bias and base composition (Fig. 4Up), the percentage GC at third base codon positions (GC3s) for each gene and ORF were plotted against Nc values for the representative west African strain Asibi and east/central African strain Ethiopia61b. Nc is a measure of codon usage bias that is independent of amino acid composition and gene length, and ranges from 20 to 61 when there is extreme bias and no bias, respectively (Wright, 1990Down). Consistent with previous studies that investigated codon usage bias in flaviviruses (Jenkins et al., 2001Down; Jenkins & Holmes, 2003Down), in the plot of Nc against GC3s for YFV genomic sequences (Fig. 4Up; Table 4Down), some points nearly overlap the curve, while others lie farther away, suggesting that base composition as well as other factors may contribute to the observed codon usage bias.


View this table:
[in this window]
[in a new window]
 
Table 4. Nc and GC3 values for representative west Africa II (Asibi) and east/central Africa (Ethiopia61b) genotypes

Values used for construction of Nc plot (see Fig. 4Up). Values for individual YFV proteins and ORF are listed; C (C protein), prM/M (pre-membrane/membrane protein), E (E protein), NS1 through NS5 (non-structural proteins 1 through 5).

 
Predicted glycosylation sites are similar between YFV genotypes
Glycosylation has been shown to be an important post-translational modification for the function of many different proteins. YFV is known to possess N-linked glycosylation sites in the prM/M, E and NS1 proteins (Chambers et al., 1990Down; Ballinger-Crabtree & Miller, 1990Down; Muylaert et al., 1996Down). Here, we used YFV genomic sequences to examine potential N-linked glycosylation sites by sequence motif prediction (Asn–X–Ser/Thr). As expected, all virus strains examined possessed two N-linked glycosylation motifs in prM at amino acid positions 13 and 29, plus an additional Asn–X–Ser/Thr sequence at position 51 that may not be glycosylated in vivo, due to the presence of a proline immediately following the motif. Consistent with previous results, a glycosylation motif was also present at position 55 of M; however, this site is also unlikely to be glycosylated due to the hydrophobicity of this region (Ballinger-Crabtree & Miller, 1990Down).

In the E protein, all YFV strains examined had one site at position 309 and one site at position 469 which, like position 55 of M, may not be glycosylated due to its location in a hydrophobic region of the protein. N-linked glycosylation sites in NS1 were at Asn 130 and Asn 208, consistent with previous findings (Muylaert et al., 1996Down). All strains also possessed two possible sites in NS4B (positions 63/64 and 89/90) and two possible sites in NS5 (positions 226/227 and 232/233).


   DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Using full-length YFV genomes, this study confirms the genotypic classification of YFV in previous studies that were based on the examination of a 688 nt region of the prM/M–E junction (Mutebi et al., 2001Down) and extends upon previous findings that suggested high levels of genetic variation among YFV isolates (Lepiniec et al., 1994Down; Chang et al., 1995Down; Wang et al., 1996Down; Mutebi et al., 2001Down; Bryant et al., 2003Down). Nucleotide and amino acid variation between genomes (Table 2Up) supports previous observations that proposed at least five genotypes of YFV in Africa. Genotype-specific amino acid substitutions (Supplementary Table S3 available in JGV Online) may aid in identification of amino acids important for the virulence of a particular strain, once phenotypic traits of these viruses have been compared.

In particular, a deletion of Thr (or His when compared to CAR77-900) identified in the N terminus of NS4B of strains from west Africa may be of potential significance (Fig. 3Up). While little is known about the structure and function of YFV NS4B, for other members of the family Flaviviridae there is evidence that the protein may play a role in vector competence (Hanley et al., 2003Down), may serve as an interferon antagonist (Munoz-Jordan et al., 2003Down; Jones et al., 2005Down) and may facilitate virus replication through association with intracellular membranes (Elazar et al., 2004Down; Gao et al., 2004Down). Furthermore, the deletion is located in an area of the protein that is highly variable for east and central African strains compared with the published west African YFV strains (Figs 2 and 3UpUp) and is predicted to be outside the membrane in hydrophobicity plots (data not shown). Comparison of this region of YFV NS4B with the same region of other flaviviruses (Fig. 3bUp) revealed that Yokose virus, a close relative of YFV, also possesses an amino acid at the position of the additional amino acid in the east and central African strains, and has an amino acid signature that is more similar to these strains than west African strains of YFV. In addition, regions of the ORF that were highly variable between eastern and central African strains, and strains from western Africa, were also observed in C and NS5 proteins. The significance of this variability is unknown; however, C protein substitutions were located in the C-terminal hydrophobic sequence that is cleaved from the protein in the mature virion (Chambers et al., 1990Down), and NS5 substitutions were located outside the polymerase motifs (Lai et al., 1999Down), which were conserved for all strains examined.

Results of codon usage analysis revealed more similar codon usage when genomic sequences were examined than when individual genes were studied (Table 3Up). The cause of codon usage bias in YFV may be a combination of factors, including mutational bias (due to genomic compositional constraints) and/or translational selection on areas of the YFV genome. GC3s were therefore plotted against Nc values and compared to the expected Nc if base composition is the only factor involved (solid curve) (Fig. 4Up). Results suggested that while base compositional constraints may account for most of the codon usage bias in some genes (e.g. NS2B of Ethiopia61b that overlaps the curve), it is unlikely to be the single determining factor, as some points are located farther from the expected Nc curve. This is consistent with previous studies suggesting that while base composition may account for the majority of codon usage bias between different flavivirus groups, other factors may also be involved (Jenkins et al., 2001Down; Jenkins & Holmes, 2003Down). It is therefore possible that translational selection, involving increased translational efficiency due to preferential usage of a codon for which there is an abundant cognate tRNA, may also contribute to the observed codon usage bias. For example, vector and vertebrate host species vary according to geographical location for YFV. Since these different cell types may exhibit very different codon/translational preferences, this could contribute to the differences in codon usage bias for various genotypes. This possibility may be unlikely, however, since previous studies did not find a correlation between codon usage bias and arthropod association (tick vs mosquito vs no vector) (Jenkins et al., 2001Down). In order to fully explore the potential importance of translational selection in the codon usage bias of YFV, experiments must be undertaken to determine the cellular levels of tRNA availability in response to viral gene expression.

Since glycosylation of YFV M, E and NS1 proteins may be important for the function of these proteins (Ruiz-Linares et al., 1989Down; Ballinger-Crabtree & Miller, 1990Down; Muylaert et al., 1996Down), we also analysed sequences of the various YFV genomes for potential N-linked glycosylation sites. Positions of N-linked glycosylation motifs in E, prM/M and NS1 were conserved for all virus strains. In addition, conserved N-linked glycosylation motifs were also detected in the NS4B and NS5 proteins. Whether or not these additional sites are glycosylated in vivo remains to be determined. While topology predictions of NS4B indicate three major hydrophobic segments (which are therefore unlikely to be glycosylated), predicted N-linked glycosylation sites are located outside these regions. YFV NS5 serves as the viral RNA polymerase and is therefore unlikely to be glycosylated.

In conclusion, this study supports the existence of previously identified genotypes of YFV, on the basis of full-genome phylogenetic analysis. High degrees of genetic diversity are found within the YFV species, while potential glycosylation sites are generally conserved. Phenotypic evaluation of differences among these viruses remains to be assessed, in particular, the potential significance of highly variable regions in C and NS5, as well as the Thr/His deletion in NS4B of west African viruses. Identification of genetic relationships and potential corresponding phenotypic diversity of YFV is important for monitoring of future outbreaks.


   ACKNOWLEDGEMENTS
 
We wish to thank Juliet Bryant and Amber Engel for helpful discussions, Robert Tesh, Robert Shope and the World Reference Center for Emerging Viruses and Arboviruses (NIH contract N01 A1 30027) for provision of virus strains and the UTMB sequence core lab (S. Smith and E. Surriga).


   REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Andral, L., Bres, P., Serie, C., Casals, J. & Panthier, R. (1968). Studies on yellow fever in Ethiopia. 3. Serological and virological studies of the woodland fauna. Bull W H O 38, 855–861.[Medline]

Ballinger-Crabtree, M. E. & Miller, B. R. (1990). Partial nucleotide sequence of South American yellow fever virus strain 1899/81: structural proteins and NS1. J Gen Virol 71, 2115–2121.[Abstract/Free Full Text]

Barrett, A. D. & Monath, T. P. (2003). Epidemiology and ecology of yellow fever virus. Adv Virus Res 61, 291–315.[CrossRef][Medline]

Bryant, J. E. & Barrett, A. D. (2003). Comparative phylogenies of yellow fever isolates from Peru and Brazil. FEMS Immunol Med Microbiol 39, 103–118.[CrossRef][Medline]

Bryant, J., Wang, H., Cabezas, C., Ramirez, G., Watts, D., Russell, K. & Barrett, A. (2003). Enzootic transmission of yellow fever virus in Peru. Emerg Infect Dis 9, 926–933.[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]

Chang, G. J., Cropp, B. C., Kinney, R. M., Trent, D. W. & Gubler, D. J. (1995). Nucleotide sequence variation of the envelope protein gene identifies two distinct genotypes of yellow fever virus. J Virol 69, 5773–5780.[Abstract]

Colebunders, R., Mariage, J. L., Coche, J. C. & 9 other authors (2002). A Belgian traveler who acquired yellow fever in the Gambia. Clin Infect Dis 35, e113–e116.[CrossRef][Medline]

Elazar, M., Liu, P., Rice, C. M. & Glenn, J. S. (2004). An N-terminal amphipathic helix in hepatitis C virus (HCV) NS4B mediates membrane association, correct localization of replication complex proteins, and HCV RNA replication. J Virol 78, 11393–11400.[Abstract/Free Full Text]

Gao, L., Aizaki, H., He, J. W. & Lai, M. M. (2004). Interactions between viral nonstructural proteins and host protein hVAP-33 mediate the formation of hepatitis C virus RNA replication complex on lipid raft. J Virol 78, 3480–3488.[Abstract/Free Full Text]

Hahn, C. S., Dalrymple, J. M., Strauss, J. H. & Rice, C. M. (1987). Comparison of the virulent Asibi strain of yellow fever virus with the 17D vaccine strain derived from it. Proc Natl Acad Sci U S A 84, 2019–2023.[Abstract/Free Full Text]

Hanley, K. A., Manlucu, L. R., Gilmore, L. E., Blaney, J. E., Jr, Hanson, C. T., Murphy, B. R. & Whitehead, S. S. (2003). A trade-off in replication in mosquito versus mammalian systems conferred by a point mutation in the NS4B protein of dengue virus type 4. Virology 312, 222–232.[CrossRef][Medline]

Jenkins, G. M. & Holmes, E. C. (2003). The extent of codon usage bias in human RNA viruses and its evolutionary origin. Virus Res 92, 1–7.[CrossRef][Medline]

Jenkins, G. M., Pagel, M., Gould, E. A., de A Zanotto, P. M. & Holmes, E. C. (2001). Evolution of base composition and codon usage bias in the genus Flavivirus. J Mol Evol 52, 383–390.[Medline]

Jones, M., Davidson, A., Hibbert, L., Gruenwald, P., Schlaak, J., Ball, S., Foster, G. R. & Jacobs, M. (2005). Dengue virus inhibits alpha interferon signaling by reducing STAT2 expression. J Virol 79, 5414–5420.[Abstract/Free Full Text]

Lai, V. C., Kao, C. C., Ferrari, E., Park, J., Uss, A. S., Wright-Minogue, J., Hong, Z. & Lau, J. Y. (1999). Mutational analysis of bovine viral diarrhea virus RNA-dependent RNA polymerase. J Virol 73, 10129–10136.[Abstract/Free Full Text]

Lepiniec, L., Dalgarno, L., Huong, V. T., Monath, T. P., Digoutte, J. P. & Deubel, V. (1994). Geographic distribution and evolution of yellow fever viruses based on direct sequencing of genomic cDNA fragments. J Gen Virol 75, 417–423.[Abstract/Free Full Text]

McInerney, J. O. (1998). GCUA: general codon usage analysis. Bioinformatics 14, 372–373.[Abstract/Free Full Text]

Monath, T. P. (1999). Facing up to re-emergence of urban yellow fever. Lancet 353, 1541.[Medline]

Munoz-Jordan, J. L., Sanchez-Burgos, G. G., Laurent-Rolle, M. & Garcia-Sastre, A. (2003). Inhibition of interferon signaling by dengue virus. Proc Natl Acad Sci U S A 100, 14333–14338.[Abstract/Free Full Text]

Mutebi, J. P., Wang, H., Li, L., Bryant, J. E. & Barrett, A. D. (2001). Phylogenetic and evolutionary relationships among yellow fever virus isolates in Africa. J Virol 75, 6999–7008.[Abstract/Free Full Text]

Mutebi, J. P., Rijnbrand, R. C. A., Wang, H., Ryman, K. D., Wang, E., Fulop, L. D., Titball, R. & Barrett, A. D. T. (2004). Genetic relationships and evolution of genotypes of yellow fever virus and other members of the yellow fever virus group within the Flavivirus genus based on the 3' noncoding region. J Virol 78, 9652–9665.[Abstract/Free Full Text]

Muylaert, I. R., Chambers, T. J., Galler, R. & Rice, C. M. (1996). Mutagenesis of the N-linked glycosylation sites of the yellow fever virus NS1 protein: effects on virus replication and mouse neurovirulence. Virology 222, 159–168.[CrossRef][Medline]

Pinto, M. R. & Filipe, A. R. (1973). Arbovirus studies in Luanda, Angola. 1. Virological and serological studies during a yellow fever epidemic. Bull W H O 49, 31–35.[Medline]

Pisano, M. R., Nicoli, J. & Tolou, H. (1997). Homogeneity of yellow fever virus strains isolated during an epidemic and a post-epidemic period in West Africa. Virus Genes 14, 225–234.[CrossRef][Medline]

Robertson, S. E., Hull, B. P., Tomori, O., Bele, O., LeDuc, J. W. & Esteves, K. (1996). Yellow fever: a decade of reemergence. JAMA 276, 1157–1162.[Abstract]

Ruiz-Linares, A., Cahour, A., Despres, P., Girard, M. & Bouloy, M. (1989). Processing of yellow fever virus polyprotein: role of cellular proteases in maturation of the structural proteins. J Virol 63, 4199–4209.[Abstract/Free Full Text]

Serie, C., Casals, J., Panthier, R., Bres, P. & Williams, M. C. (1968). Studies on yellow fever in Ethiopia. 2. Serological study of the human population. Bull W H O 38, 843–854.[Medline]

Smithburn, K. C., Haddow, A. J. & Lumsden, W. H. R. (1949). An outbreak of sylvan yellow fever in Uganda with Aedes (Stegomyia) africanus Theobald as principal vector and insect host of the virus. Am J Trop Med 43, 74–89.

Swofford, D. L. (2003). PAUP*: Phylogenetic analysis using parsimony (* and other methods), version 4.0b10. Sunderland, MA: Sinauer Associates.

Tomori, O. (1999). Impact of yellow fever on the developing world. Adv Virus Res 53, 5–34.[Medline]

Wang, E., Ryman, K. D., Jennings, A. D., Wood, D. J., Taffs, F., Minor, P. D., Sanders, P. G. & Barrett, A. D. (1995). Comparison of the genomes of the wild-type French viscerotropic strain of yellow fever virus with its vaccine derivative French neurotropic vaccine. J Gen Virol 76, 2749–2755.[Abstract/Free Full Text]

Wang, E., Weaver, S. C., Shope, R. E., Tesh, R. B., Watts, D. M. & Barrett, A. D. (1996). Genetic variation in yellow fever virus: duplication in the 3' noncoding region of strains from Africa. Virology 225, 274–281.[CrossRef][Medline]

Wright, F. (1990). The ‘effective number of codons' used in a gene. Gene 87, 23–29.[CrossRef][Medline]

Received 3 June 2005; accepted 1 December 2005.



This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary tables
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 CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by von Lindern, J. J.
Right arrow Articles by Barrett, A. D. T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by von Lindern, J. J.
Right arrow Articles by Barrett, A. D. T.
Agricola
Right arrow Articles by von Lindern, J. J.
Right arrow Articles by Barrett, A. D. T.


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