| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
1 Wadsworth Center, New York State Department of Health, Albany, NY 12201, USA
2 Department of Biomedical Sciences, University at Albany, State University of New York, Albany, NY 12201, USA
Correspondence
Alexander T. Ciota
aciota{at}wadsworth.org
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
). Evolutionary pressures on the virus, therefore, are applied predominately by the mosquito and avian environments. Since its introduction, WNV has spread throughout the USA, Canada, Mexico, the Caribbean and Central and South America (Austin, et al., 2004; Cruz et al., 2005; Dupuis et al., 2003, 2005; Elizondo-Quiroga, 2005; Granwehr et al., 2004; Morales et al., 2006) and has infected a broad range of both mosquito and avian species (Hayes et al., 2005). There have been some recent population level changes in WNV accompanying its expansion in the western hemisphere, but the virus remains highly conserved (Anderson et al., 2001; Beasley et al., 2003; Davis et al., 2005; Ebel et al., 2001, 2004; Lanciotti et al., 2002). This lack of evolutionary change over space and time is common with many arboviruses (Cilnis et al., 1996; Jenkins et al., 2002; Weaver et al., 1992) and is frequently attributed to the adaptive constraints resulting from the virus cycling between vertebrate hosts and invertebrate vectors (Scott et al., 1994; Weaver et al., 1992). This hypothesis is consistent with studies that suggest that WNV and other arboviruses are subject to high levels of purifying selection (Holmes, 2003; Jerzak et al., 2005; Woelk & Holmes, 2002). Despite this genetic stability, WNV has been highly successful in its emergence and establishment in North America (Davis et al., 2005).
WNV, like most RNA viruses studied, has been shown in nature to be a highly genetically diverse population of genotypes around a dominant master sequence (Jerzak et al., 2005). This finding indicates that WNV is most likely perpetuated in nature as a quasispecies. The term quasispecies refers to the diverse mutant spectrum surrounding a master genotype that presumably possesses the highest fitness value (Domingo et al., 1998; Eigen, 1993; Eigen & Biebricher, 1988; Holmes & Moya, 2002). Genetic diversity within virus populations results from the high error rates and lack of proofreading ability of the RNA-dependent RNA polymerase (Holland et al., 1982), together with high rates of virus replication and large population sizes (Drake & Holland, 1999). The quasispecies theory concludes that selection acts on viral populations as a whole, rather than on individual virions (Eigen, 1993; Martinez et al., 1991; Ruiz-Jarabo et al., 2000). Quasispecies have been shown to be important in viral fitness in previous studies (Duarte et al., 1994; Martinez et al., 1991; Ruiz-Jarabo et al., 2002). In studies with vesicular stomatitis virus (VSV), fitness increases have been measured with no changes in consensus sequences (Novella & Ebendick-Corpus, 2004). The existence of genetic diversity within WNV populations therefore suggests that intrapopulation variants may play a large role in viral fitness and replication, and in the ability of WNV to adapt successfully to new and changing environments.
We reported previously that WNV serially passaged in mosquito (C6/36) cell culture displayed large gains in fitness and replicative ability relative to parental virus populations when measured in mosquito cells. However, this did not correspond to similar changes when measured in avian (DF-1) cell culture (Ciota et al., 2007). These cell-specific phenotypic changes were coupled with modest genetic variation in the consensus sequence. Specifically, two non-synonymous changes were identified by full-genome analysis following 39 mosquito cell passages. These included an A
G mutation at position 1712 of the envelope region resulting in a LysArg substitution, and a GT mutation at position 6687 of the NS4A region resulting in a HisGln substitution. The envelope change is fairly conservative in nature and, given its position and orientation, could not be clearly implicated in any obvious function such as receptor binding or membrane fusion. With no known role for NS4A in replication, it is also difficult to implicate the change there in any specific function, although the amino acid substitution (positively charged residue to polar uncharged residue) is potentially significant. Follow-up studies are presented here with the same passage series using: (i) molecular cloning analyses to evaluate the extent of non-consensus variation that correlated with adaptation, and (ii) reverse genetics to evaluate whether replicative advantages in C6/36 cells were a result of consensus mutations alone or whether the mutant spectrum had a significant impact in generating the adaptive phenotype. These are important questions, as the role of genotypic diversity in viral adaptation has not been examined directly. We also investigated whether WNV that had been serially passaged in avian cells displayed similar cell-specific replicative advantages and changes in genetic diversity. Understanding quasispecies dynamics and the resultant phenotypes of viruses in cell culture systems provides a basis to evaluate in vivo results and is the first step towards understanding the importance of this phenomenon in viral emergence and evolution.
|
METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
Virus strains.
WNV NY003356 is a primary isolate from the kidney tissue of an American crow that was collected in 2000 in Staten Island, NY, USA, and was prepared by one round of amplification in Vero cells (Ebel et al., 2001). Biological clones of this strain (3356.1.1.1
[EC]
) were generated by three rounds of plaque purification on Vero cells in order to establish a homogeneous population for all studies. The biological clone was propagated for use in passaging and growth studies by two rounds of Vero cell amplification and was quantified by plaque assay on Vero cells (Payne et al., 2006). A full-length WNV infectious clone (FL-WNV) was constructed as described previously (Shi et al., 2002b) based on the sequence of the same isolate. Multiple aliquots of all viral stocks were stored at –80 °C.
Serial passage of virus.
The biological clone was passaged 40 (C6/36) or 20 (DF-1) times sequentially in cell culture, using an m.o.i. of 0.1 p.f.u. per cell (approx. 5 log10 p.f.u. virus), based on Vero cell titres, for each passage as described previously (Ciota et al., 2007). The C6/36 passage series used in this study was the same series used in the Ciota et al. (2007) study. Briefly, virus was adsorbed to confluent cell monolayers in six-well plates for 1 h at 28 °C (C6/36) or 39 °C (DF-1), with frequent rocking. Following adsorption, the inoculum was removed, cells were washed with MEM, 3 ml maintenance medium was added and plates were incubated at appropriate temperatures. Medium from each well was harvested at 72 h post-inoculation (p.i.) and multiple aliquots were stored at –80 °C. Following each passage, the viral harvest was quantified by plaque assay on Vero cells and diluted to produce an m.o.i. of 0.1 p.f.u. per cell for each subsequent passage. Virus titres varied during passage, but at no time was there greater than a 0.5 log10 difference between C6/36- and DF-1-derived virus. Alternate-passaged WNV was derived from 39 serial passages in C6/36 cells followed by a single passage in DF-1 cells.
Viral growth curve analysis.
Six-well plates containing confluent monolayers of Vero, C6/36 or DF-1 cells were infected with virus, in duplicate, at an m.o.i. of 0.01 p.f.u. per cell, based on Vero cell titres. At the end of a 1 h absorption period at 37 °C (Vero), 28 °C (C6/36) or 39 °C (DF-1), the inoculum was removed, wells were washed and 3 ml of maintenance medium were added to each well. Samples, consisting of 100 µl media, were taken in duplicate at 12, 24, 36, 48, 72, 96 and 120 h p.i., diluted 1 : 10 in BA-1 medium containing 20 % FBS and stored at –80 °C. All samples were titrated in duplicate by plaque assay on Vero cells and growth curves were constructed using the mean titre for each time point.
Plasmid construction.
Standard site-specific mutagenesis by overlap extension was used for the introduction of mutations into a modified full-length WNV cDNA clone (pFLWNV; Shi et al., 2002a). Primers 5'-ACCACACGCCACGAGGCAGTCTGTGATAG-3', 5'-CTATCACAGACTGCCTCGTGGCGTGTGGTTC-3', 5'-CATACCAACAGCCGCTGTCAGTT-3' and 5'-CAAGATTCCGAATACCGCAAG-3' were used to make the single-nucleotide mutant A1712G. The modified pFLWNV vector and the final PCR product were digested with ClaI and StuI. The 1.5 kb fragment was inserted into the digested vector. Primers 5'-TCCTCCTCATGCATCGGAAGGGCATTG-3', 5'-CTCTATGACTTTCGGCATGTAGG-3', 5'-AATGCCCTTCCGATGCATGAGGAG-3' and 5'-GGTTTGGTTTGTGCCTAGTG-3' were used to make the single-nucleotide mutant G6687T. The PCR product was digested with KpnI and the 2.4 kb fragment was inserted into the digested vector. The full-length WNV cDNA vector containing the G6687T mutation was digested with ClaI and StuI and ligated with the 1.5 kb fragment above to make the double mutant. Before ligation, the 1.5 and 2.4 kb fragments were sequenced to confirm the mutations. Sequencing of the plasmids of pFLWNV containing mutations indicated that no other changes were present.
In vitro transcription of RNA and transfection.
Transcription and transfection of RNAs were performed as described previously (Shi et al., 2002a). Briefly, plasmids were purified using a MaxiPrep kit (Qiagen) and digested with XbaI. The linearized plasmids were extracted and resuspended in RNase-free water. An mMessage mMachine T7 kit (Ambion) was used to transcribe full-length WNV RNA in vitro and DNA template was removed by DNase I digestion. RNA was quantified by spectrophotometer and 10 µg of wild-type or mutant RNA was electroporated into 8x106 BHK-21 cells. Transfected cells were seeded into T-150 flasks. Cell supernatants were collected on days 1 to 5 post-transfection and WNV infectious particles were quantified by plaque assay on Vero cells.
Sequencing of mutant strains.
RNA was extracted from WNV using RNeasy spin columns (Qiagen) according to manufacturer’s protocol. Primers for WNV were designed from GenBank accession no. AF260967
[GenBank]
. One-step RT-PCR (Qiagen) was conducted using primers to generate overlapping PCR products. Reverse transcription reactions were carried out at 50 °C for 30 min, followed by inactivation of the transcriptase at 95 °C for 15 min. Amplification was then carried out for 40 cycles of 94 °C for 20 s, 55 °C for 30 s and 72 °C for 2 min, with a final elongation at 72 °C for 10 min. PCR products were visualized on a 1.5 % agarose gel and bands were allowed to run through 1 % Nusieve GTG low melting point agarose (Cambrex BioScience). Sequencing was performed using ABI 3700 automated sequencers (Applied Biosystems) using overlapping primers to give a minimum of twofold redundancy. Sequences were compiled and edited by using the DNASTAR software package.
High-fidelity RT-PCR, cloning and sequencing.
Production and analysis of clones was performed as described previously (Jerzak et al., 2005). RNA was extracted from infected specimens using RNeasy spin columns (Qiagen) and RT-PCR was conducted using primers designed to amplify the 3' 1159 nt of the WNV envelope (E) coding region and the 5' 779 nt of the WNV non-structural protein 1 (NS1) coding region (forward primer WNV1311: 5'-ATGCGCCAAATTTGCCTGCTCTAC-3'; reverse primer WNV3248: 5'-ATGGGCCCTGGTTTTGTGTCTTGT-3'). Reverse transcription of 5 µl RNA was performed with Sensiscript reverse transcriptase (Qiagen) at 45 °C for 40 min. Reverse transcriptase reactions were followed by heat inactivation at 95 °C for 5 min. The resulting cDNA was used as template for PCR amplification. WNV cDNA was then amplified with a ‘high-fidelity’ protocol using PfuUltra (Stratagene), according to the manufacturer's specifications. Amplification was carried out for 40 cycles of 94 °C for 30 s, 50 °C for 30 s and 72 °C for 4 min, with final extension at 72 °C for 10 min. PCR products were visualized on a 1.5 % agarose gel and DNA was recovered using a MinElute gel extraction kit (Qiagen) as specified by the manufacturer. The recovered DNA was ligated into the cloning vector pCR-Script Amp SK(+) and transformed into XL10-Gold ultracompetent cells (Stratagene) according to the manufacturer’s protocol. The blue–white colour-screening method was used to select transformed colonies. White colonies were screened by direct PCR using primers specific for the desired insert. Plasmid DNA was purified by using a QIAprep Spin Miniprep kit (Qiagen) as specified by the manufacturer. Sequencing was carried out using five pairs of overlapping primers and the T3 reverse primer. Sequencing was performed at the Wadsworth Center Molecular Genetics Core using ABI 3700 and 3100 automated sequencers (Applied Biosystems). Between 15 and 20 clones were sequenced per sample.
Sequence analysis of populations.
Sequences were compiled and edited using the SEQMAN module of the DNASTAR software package, and a minimum of twofold redundancy throughout each clone was required for sequence data to be considered complete. Between 15 and 20 clones from each individual sample were aligned using MEGALIGN within DNASTAR. The consensus sequence for each sample was determined and the sequence of each clone was compared with the consensus. Consensus sequences attained by averaging cloning products agreed with those attained by direct PCR sequencing. The percentage of nucleotide mutations (total number of mutations divided by total number of bases sequenced), amino acid mutations (total number of amino acid changes divided by total number of amino acids sequenced) and the sequence diversity (percentage of clones with at least one difference from the consensus) were used as indicators of genetic diversity. Normalized Shannon entropy (Sn) was calculated based on the frequency of genotypes in populations as follows: Shannon entropy (Sn)=
–i Pi lnPi/lnN, where Pi is the frequency of an individual genotype and N is the number of clones sequenced. Sn values range from 0 (completely homogeneous) to 1 (completely heterogeneous). In order to assess selective pressures on populations, dN/dS ratios were calculated using DNASP (Rozas & Rozas, 1999) and Microsoft EXCEL, as described previously (Jerzak et al., 2005). Statistical analyses were performed using both Microsoft EXCEL 2003 and GRAPHPAD PRISM version 4.00.
|
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
). This represented just one mutation in 18 clones analysed (Table 1). A low level of homogeneity was anticipated for this population of clones, as the virus strain was itself a biological clone. This result confirmed the high fidelity of the PfuUltra enzyme used for RT-PCRs and allowed us to have confidence in the low levels of misincorporations in subsequent reactions. Overall, the genetic diversity increased in a reasonably linear fashion with serial passage in C6/36 cell culture (Fig. 1). The level of nucleotide diversity identified for CP19 (0.031 %) was significantly higher than that measured for UNP (Fisher’s exact test, P=0.012) and continued to increase significantly from CP19 to CP39 (P=0.022). Fourteen different genotypes in the cloned region were identified as resulting from C6/36 passage (Table 1). Levels of entropy (Sn) within populations based on the frequency of new genotypes also increased with C6/36 passage, progressing from 0.07 for UNP (highly homogeneous) to 0.70 for CP39 (highly heterogeneous). This represented a total of eight genotypes identified for CP39 with analysis of just 15 clones. As indicated by values for amino acid diversity (Fig. 1b), the majority of mutations were non-synonymous. For the entire region analysed, the dN/dS ratio for CP39 was 0.94, which is very close to the 1.0 predicted by random accumulation of mutations without significant positive or purifying selection (Holmes, 2003). For the E coding region alone, however, 100 % of mutations identified (6/6) were non-synonymous. Furthermore, the two synonymous changes in the NS1 region were only identified in 2/15 clones, whereas the non-synonymous changes found in the E region were identified in 13/15 clones (Table 1). The majority of clones (57 %) differed from the consensus sequence in at least one position by passage 19. Despite this, the single most common genotype (master sequence) remained the consensus sequence (genotype C). By passage 39, 87 % of clones differed from the consensus sequence, and the master sequence (genotype I) contained a non-conservative amino acid change in the envelope protein (F432S) not identified by consensus analysis.
|
|
). Comparisons of overall nucleotide diversity, amino acid diversity, nucleotide sequence diversity and amino acid sequence diversity all indicated statistically significant decreases from CP39 to AP1 (Fisher’s exact test, P<0.05). Nucleotide diversity decreased from 0.08 to 0.02 % and amino acid diversity decreased from 0.22 to 0.03 % (Fig. 1a). The percentage of clones with consensus genotypes increased from 13 % for CP39 to 67 % for AP1 (Fig. 1b) and levels of entropy decreased from 0.70 to 0.36 (Table 1). These values indicated a significant increase in genotypic homogeneity from CP39 after just 48 h of replication in DF-1 cells. Genotypes C (consensus), K and L were identified in both CP39 and AP1. No additional mutations were identified with AP1.
Although measurements of diversity on all levels (nucleotide, amino acid and sequences) for DP20 were slightly lower than those measured for CP19 (Fig. 1), none were statistically significant (Fisher’s exact test). Nine out of 18 clones analysed for DP20 represented the consensus genotype, while 6/15 for CP19 represented the consensus (Table 1). Calculations of entropy for DP20 (0.47) also indicated a slightly lower genotypic heterogeneity compared with CP19 (0.54).
Viral growth patterns of WNV following cell-culture passage
The replicative advantage measured by growth analyses in C6/36 for WNV CP40 was similar to the advantage reported previously following 39 passages (Fig. 2a; Ciota et al., 2007). However, a single passage in DF-1 cells following 39 passages in C6/36 (AP1), significantly decreased the WNV titre compared with CP40 at 24, 48 and 72 h p.i. in C6/36 cells (Fig. 2a; paired t-test, P<0.05), but not in DF-1 cells (Fig. 2b). There was a decrease in the rate of replication of AP1 in C6/36 cells compared with CP40 up to 96 h p.i. However, AP1 titres in C6/36 cells did remain significantly higher than those of UNP up to 120 h p.i (P<0.05).
|
) or DF-1 cell culture (Fig. 3b) compared with DF-1 passage and UNP (paired t-test). This finding differed from the substantial cell-specific increases in peak growth titre observed following passage of WNV in C6/36 cells (Fig. 2a). However, both DP10 and DP20 reached their peak titre a day earlier than UNP in DF-1 cell culture, and virus titres were significantly higher at 48 h p.i. for DP10 and at both 36 and 48 h p.i. for DP20 relative to UNP (paired t-test, P<0.01). Virus titres of DP20 at 72 h p.i. and DP10 at 96 h p.i. were significantly lower than UNP (P<0.01). Although statistically significant, the biological significance of these differences is uncertain. Only one titre difference (DP20 vs UNP at 36 h p.i.) was greater than 1 log10 p.f.u. ml–1. Nevertheless, these data suggest a slightly increased rate of growth in DF-1 cell culture following passage in DF-1.
|
). Although virus titres were slightly higher for WNV A1712G in C6/36 cells, these differences were only statistically significant (paired t-test, P<0.05) at 12 h p.i. Significantly lower titres (paired t-test, P<0.05) for both WNV A1712G and WNV E/NS4 were measured in DF-1 cell culture when compared with wild-type strain at 72 h p.i. to 120 h p.i. (Fig. 4). No differences in growth were measured between WNV G6687T in DF-1 cell culture and the wild-type strain. All trends were verified in replicate growth analyses (data not shown).
|
|
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
; Beasley et al., 2003; Davis et al., 2005; Ebel et al., 2001, 2004; Lanciotti et al., 2002), even as it has expanded its range in the naive environment of North America (Davis et al., 2005). Worldwide, WNV has infected over 75 species of mosquitoes (Higgs et al., 2004) and over 300 species of birds (Marra et al., 2003). Although it is clear that WNV, like all RNA viruses, has the capacity for rapid accumulation of mutations (Drake & Holland, 1999), the impact of non-consensus genetic variants and quasispecies populations on WNV adaptation has not been studied previously.
We reported previously that WNV is capable of cell-specific adaptation with minimal genetic alterations in consensus sequences following passage in mosquito cell culture (Ciota et al., 2007). This cell-specific adaptation results in a significant replicative advantage in mosquito cells (Fig. 2). Here, we have shown that despite minimal consensus change, multiple intrapopulation genotype variants were generated with mosquito cell passage (Table 1) and that the size of the mutant spectrum increased with adaptation to mosquito cells (Fig. 1, Table 1). The level of variation was such that the consensus genotype was no longer the dominant genotype by passage 39. In fact, as the clones sequenced here represent only a portion of the genome (2 kb), it seems likely that a complete genome identical to the consensus sequence would be quite rare. As population diversity increases, it becomes more likely that the consensus sequence may not represent the most common genotype (master sequence). Given that the complete set of sequences and not individual mutations represent the unit of selection, it is very important to address this when studying viral adaptation and evolution. This is particularly true with WNV, where unlike reports for some other flaviviruses, no evidence of recombination has been noted (Holmes et al., 1999; Twiddy & Holmes, 2003). The most common sequence in the C6/36 passage 39 population that we identified was genotype I (Table 1). This genotype includes a TC change at position 1391 of the E coding region, which translates into a potentially significant phenylalanine to serine change in almost 50 % of the population. These data demonstrated that attempting to identify important mutations using consensus change alone can be misleading.
The generation of quasispecies populations in other RNA viruses such as human immunodeficiency virus and Hepatitis C virus has been attributed to the avoidance of host immune responses (Farci et al., 2000, 2002; Frost et al., 2005; Geffin et al., 2003; Weiner et al., 1995). Studies with Foot-and-mouth disease virus (FMDV) have demonstrated that molecular memory is a function of quasispecies structure (Arias et al., 2004; Ruiz-Jarabo et al., 2000). These studies point to the idea that the accumulation of diversity results from variable selective pressures, which, in turn, allow viruses to succeed in dynamic environments. Results presented here clearly indicate that such factors are not a prerequisite for the accumulation of diverse genotypes within populations. This result is consistent with findings from another study with FMDV, which identified genetically diverse populations following sequential cell culture passage (Arias et al., 2001). Diverse genotypes within in vitro-passaged viruses could represent tolerance for neutral mutations or even somewhat deleterious mutations, due to the low levels of selection and lack of bottlenecks that occur with this type of large population passage. Given the well-documented high mutation rate of RNA viruses (approx. 1 mutation per 104 bp replicated), it is easy to see how such high levels of diversity could accumulate in the absence of selection (Drake & Holland, 1999). A previous study with a high-fitness variant of VSV obtained from cell-culture passage showed high levels of phenotypic heterogeneity within the population, with the majority of subclones possessing fitness values lower than that of the parental population (Duarte et al., 1994). In our study, the dN/dS ratio for mutations generated with C6/36 passage (0.94) was not significantly different from what would be predicted by accumulation of random mutations in the absence of high levels of purifying or positive selection (Holmes, 2003). However, for CP39, the majority of the change was found in the E coding region (Table 1) and all of the mutations in this region were non-synonymous. This indicates that positive selection is probably acting on the envelope region, where as purifying selection is acting on the NS1 region, and that selective pressures can vary drastically between regions for WNV. One possible explanation for the apparent contradicting results of high diversity together with selection is that the mutant spectrum somehow provides a replicative advantage and that the unit of selection, as predicted by the concept of quasispecies, is the entire population (Eigen, 1993). In a study by Jerzak et al. (2005), which analysed the same region, no significant bias was identified for mutation in the envelope region in natural isolates. These discrepant results can probably be attributed to the fact that, in the present study, we measured changes in an adapting population. It is not surprising that changes in the E protein are important in creating this adaptive phenotype given this region’s well-defined role in binding and replication (Chambers et al., 1990; Scherret et al., 2001), yet full-genome analysis of large numbers of clones would be necessary to characterize fully the selective pressures on the quasispecies population.
The finding of less cell-specific adaptation (Figs 2a, 3a) with avian cell-passaged WNV relative to mosquito cell-passaged WNV suggests that these disparate cell lines represent different environments for WNV replication and adaptation. Jerzak et al. (2005) demonstrated that mosquito-derived WNV is more heterogeneous than avian-derived WNV in field isolates. The trend we observed towards less variability for WNV in avian cells provides further evidence that the size of the mutant spectrum for WNV may at least be partially dependent on the host cell. Studies with the plant viruses Tobacco mosaic virus and Cucumber mosaic virus have demonstrated host-dependent effects of population diversity (Schneider & Roossinck, 2001). Inherent differences in selective pressures or replicase fidelity between different environments may explain host cell-dependent effects. Further studies with larger sample sizes of WNV populations are necessary for an accurate assessment of this phenomenon and its potential effects on adaptation.
Previous studies have also attempted to measure the level of evolutionary constraint that is imposed by alternating environments through assessment of changes in consensus sequence alone (Chen et al., 2003; Greene et al., 2005; Novella et al., 1999; Weaver et al., 1999). These studies have drawn conclusions using small numbers of mutations without consideration of potentially significant change within the mutant spectra. Here, we have shown that the levels of genetic diversity generated by sequential passage in mosquito cell culture are significantly decreased with a single passage in the ‘bypassed’ alternate cell line (DF-1) (Fig. 1). Furthermore, the phenotypic advantage measured in mosquito cell culture was significantly decreased after a single passage in avian cell culture (Fig. 2). Although the complexity of in vivo systems and differences in replicative strategy between hosts are critical factors to consider when assessing the factors affecting viral evolution, the data presented here support the idea that differential selection in alternating cellular environments constrains the level of WNV mutation maintained in the virus population and, therefore, the rate of genetic change.
Many previous studies have identified consensus mutations associated with alterations in viral fitness, but none has evaluated directly the importance of such changes in generating these phenotypes (Chen et al., 2003; Greene et al., 2005; Novella et al., 1999; Weaver et al., 1999). Our results clearly indicate that the cell-specific replicative advantage acquired with passaging of WNV is not solely the result of consensus change in amino acid sequences (Fig. 4), and, therefore, implicate minority mutations in the adaptive phenotype of the passaged population. Neither the single nor double mutants displayed a significant replicative advantage in any cell line studied. In fact, in avian cells, the A1712G mutation seemed to be less stable in response to the shifting pH seen in the media at later time points (Fig. 4). Given that this difference was not observed with replication of the heterogeneous CP40 population, this result also implicates non-consensus variants in the adaptability of WNV to disparate environments. The phenotypic importance of the mutant spectrum has been noted in previous reports with other viruses (Farci et al., 2002; Sauder et al., 2006). In our study, adaptation must be conferred either exclusively from individual, non-consensus variants within the population possessing high fitness levels or as a result of complementation between genotypes that increases the overall fitness of the population. A recent study with Dengue virus suggests that a defective genome can be maintained within the mutant spectrum by complementation of functional genomes (Aaskov et al., 2006). A separate study with Poliovirus provides evidence for the role of complementation and cooperative interactions in pathogenesis (Vignuzzi et al., 2006). The fitness level conferred by the combination of variants of the adapted WNV population studied here could potentially be greater than that of any individual variant, yet evaluation of the fitness of multiple variants within the population is necessary before the source of adaptation can be characterized accurately. These data provide us with evidence of the importance of the mutant spectrum in generating phenotypes in WNV and point to quasispecies dynamics as an explanation of how this virus has succeeded in changing environments and hosts with little evidence of significant genetic change.
|
ACKNOWLEDGEMENTS |
|---|
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
Anderson, J. F., Vossbrinck, C. R., Andreadis, T. G., Iton, A., Beckwith, W. H. & Mayo, D. R. (2001). Characterization of West Nile virus from five species of mosquitoes, nine species of birds, and one mammal. Ann N Y Acad Sci 951, 328–331.
Arias, A., Lazaro, E., Escarmis, C. & Domingo, E. (2001). Molecular intermediates of fitness gain of an RNA virus: characterization of a mutant spectrum by biological and molecular cloning. J Gen Virol 82, 1049–1060.
Arias, A., Ruiz-Jarabo, C. M., Escarmis, C. & Domingo, E. (2004). Fitness increase of memory genomes in a viral quasispecies. J Mol Biol 339, 405–412.[CrossRef][Medline]
Austin, R. J., Whiting, T. L., Anderson, R. A. & Drebot, M. A. (2004). An outbreak of West Nile virus-associated disease in domestic geese (Anser anser domesticus) upon initial introduction to a geographic region, with evidence of bird to bird transmission. Can Vet J 45, 117–123.[Medline]
Beasley, D. W., Davis, C. T., Guzman, H., Vanlandingham, D. L., Travassos da Rosa, A. P. A., Parsons, R. E., Higgs, S., Tesh, R. B. & Barrett, A. D. T. (2003). Limited evolution of West Nile virus has occurred during its southwesterly spread in the United States. Virology 309, 190–195.[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]
Chen, W.-J., Wu, H.-R. & Chiou, S.-S. (2003). E/NS1 modifications of dengue 2 virus after serial passages in mammalian and/or mosquito cells. Intervirology 46, 289–295.[CrossRef][Medline]
Cilnis, M. J., Kang, W. & Weaver, S. C. (1996). Genetic conservation of Highlands J viruses. Virology 218, 343–351.[CrossRef][Medline]
Ciota, A. T., Lovelace, A. O., Ngo, K. A., Le, A. N., Maffei, J. G., Franke, M. A., Payne, A. F., Jones, S. A., Kauffman, E. B. & Kramer, L. D. (2007). Cell-specific adaptation of two flaviviruses following serial passage in mosquito cell culture. Virology 357, 165–174.[CrossRef][Medline]
Cruz, L., Cardenas, V. M., Abarca, M., Rodriguez, T., Reyna, R. F., Serpas, M. V., Fontaine, R. E., Beasley, D. W. & Da Rosa, A. P. (2005). Short report: serological evidence of West Nile virus activity in El Salvador. Am J Trop Med Hyg 72, 612–615.
Davis, C. T., Ebel, G. D., Lanciotti, R. S.,, Brault, A. C., Guzman, H., Siirin, M., Lambert, A., Parsons, R. E. & Beasley, D. W. (2005). Phylogenetic analysis of North American West Nile virus isolates, 2001–2004: evidence for the emergence of a dominant genotype. Virology 342, 252–265.[CrossRef][Medline]
Domingo, E., Escarmis, C., Sevilla, N. & Baranowski, E. (1998). Population dynamics in the evolution of RNA viruses. Adv Exp Med Biol 440, 721–727.[Medline]
Drake, J. W. & Holland, J. J. (1999). Mutation rates among RNA viruses. Proc Natl Acad Sci U S A 96, 13910–13913.
Duarte, E. A., Novella, I. S., Ledesma, S., Clarke, D. K., Moya, A., Elena, S. F., Domingo, E. & Holland, J. J. (1994). Subclonal components of consensus fitness in an RNA virus clone. J Virol 68, 4295–4301.
Dupuis, A. P., II, Marra, P. P. & Kramer, L. D. (2003). Serologic evidence of West Nile virus transmission, Jamaica, West Indies. Emerg Infect Dis 9, 860–863.[Medline]
Dupuis, A. P., II, Marra, P. P., Reitsma, R., Jones, M. J., Louie, K. L. & Kramer, L. D. (2005). Serologic evidence for West Nile virus transmission in Puerto Rico and Cuba. Am J Trop Med Hyg 73, 474–476.
Ebel, G. D., Dupuis, A. P., II, Ngo, K. A., Nicholas, D., Kauffman, E., Jones, S. A., Young, D., Maffei, J., Shi, P.-Y. & other authors (2001). Partial genetic characterization of West Nile virus strains, New York State. Emerg Infect Dis 7, 650–653.[Medline]
Ebel, G. D., Carricaburu, J., Young, D., Bernard, K. A. & Kramer, L. D. (2004). Genetic and phenotypic variation of West Nile virus in New York, 2000–2003. Am J Trop Med Hyg 71, 493–500.
Eigen, M. (1993). Viral quasispecies. Sci Am 269, 42–49.[Medline]
Eigen, M. & Biebricher, D. K. (1988). Sequence space and quasispecies distribution. In RNA Genetics, pp. 211–245. Edited by E. Domingo, J. J. Holland & P. Ahlquist. Boca Raton, FL: CRC Press.
Elizondo-Quiroga, D., Davis, C. T., Fernandez-Salas, I., Escobar-Lopez, R., Olmos, D. V., Gastalum, L. C. S., Acosta, M. A., Elizondo-Quiroga, A., Gonzalez-Rojas, J. I. & other authors (2005). West Nile virus isolation in human and mosquitoes, Mexico. Emerg Infect Dis 11, 1449–1452.[Medline]
Farci, P., Shimoda, A., Coiana, A., Diaz, G., Peddis, G., Melpolder, J. C., Strazzera, A., Chien, D. Y., Munoz, S. J. & other authors (2000). The outcome of acute hepatitis C predicted by the evolution of the viral quasispecies. Science 288, 339–344.
Farci, P., Strazzera, R., Alter, H. J., Farci, S., Degioannis, D., Coiana, A., Peddis, G., Usai, F., Serra, G. & other authors (2002). Early changes in hepatitis C viral quasispecies during interferon therapy predict the therapeutic outcome. Proc Natl Acad Sci U S A 99, 3081–3086.
Frost, S. D. W., Wrin, T., Smith, D. M., Kosakovsky Pond, S. L., Liu, Y., Paxinos, E., Chappey, C., Galovich, J., Beauchaine, J. & other authors (2005). Neutralizing antibody responses drive the evolution of human immunodeficiency virus type 1 envelope during recent HIV infection. Proc Natl Acad Sci U S A 102, 18514–18519.
Geffin, R., Hutto, C., Andrew, C. & Scott, G. B. (2003). A longitudinal assessment of autologous neutralizing antibodies in children perinatally infected with human immunodeficiency virus type 1. Virology 310, 207–215.[CrossRef][Medline]
Granwehr, B. P., Lillibridge, K. M., Higgs, S., Mason, P. W., Aronson, J. F., Campbell, G. A. & Barrett, A. D. T. (2004). West Nile virus: where are we now?. Lancet Infect Dis 4, 547–556.[CrossRef][Medline]
Greene, I. P., Wang, E., Deardorff, E. R., Milleron, R., Domingo, E. & Weaver, S. C. (2005). Effect of alternating passage on adaptation of Sindbis virus to vertebrate and invertebrate cells. J Virol 79, 14253–14260.
Hayes, E. B., Komar, N., Nasci, R. S., Montgomery, S. P., O'Leary, D. R. & Campbell, G. L. (2005). Epidemiology and transmission dynamics of West Nile virus disease. Emerg Infect Dis 11, 1167–1173.[Medline]
Higgs, S., Snow, K. & Gould, E. A. (2004). The potential for West Nile virus to establish outside of its natural range: a consideration of potential mosquito vectors in the United Kingdom. Trans R Soc Trop Med Hyg 98, 82–87.[CrossRef][Medline]
Holland, J. J., Spindler, K., Horodyski, F., Grabau, E., Nichol, S. & VandePol, S. (1982). Rapid evolution of RNA genomes. Science 215, 1577–1585.
Holmes, E. C. (2003). Patterns of intra- and interhost nonsynonymous variation reveal strong purifying selection in dengue virus. J Virol 77, 11296–11298.
Holmes, E. C. & Moya, A. (2002). Is the quasispecies concept relevant to RNA viruses?. J Virol 76, 460–465.
Holmes, E. C., Worobey, M. & Rambaut, A. (1999). Phylogenetic evidence for recombination in dengue virus. Mol Biol Evol 16, 405–409.[Abstract]
Jenkins, G. M., Rambaut, A., Pybus, O. G. & Holmes, E. C. (2002). Rates of molecular evolution in RNA viruses: a quantitative phylogenetic analysis. J Mol Evol 54, 156–165.[CrossRef][Medline]
Jerzak, G., Bernard, K. A., Kramer, L. D. & Ebel, G. D. (2005). Genetic variation in West Nile virus from naturally infected mosquitoes and birds suggests quasispecies structure and strong purifying selection. J Gen Virol 86, 2175–2183.
Kramer, L. D. & Bernard, K. A. (2001). West Nile virus in the western hemisphere. Curr Opin Infect Dis 14, 519–525.[Medline]
Lanciotti, R. S., Ebel, G. D., Deubel, V., Kerst, A. J., Murri, S., Meyer, R., Bowen, M., McKinney, N., Morrill, W. E. & 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]
Marra, P. P., Griffing, S. M. & McLean, R. G. (2003). West Nile virus and wildlife health. Emerg Infect Dis 9, 898–899.[Medline]
Martinez, M. A., Carrillo, C., Gonzalez-Candelas, F., Moya, A., Domingo, E. & Sobrino, F. (1991). Fitness alteration of foot-and-mouth disease virus mutants: measurement of adaptability of viral quasispecies. J Virol 65, 3954–3957.
Morales, M. A., Barrandeguy, M., Fabbri, C., Garcia, J. B., Vissani, A., Trono, K., Gutierrez, G., Pigretti, S., Menchaca, H. & other authors (2006). West Nile virus isolation from equines in Argentina, 2006. Emerg Infect Dis 12, 1559–1561.[Medline]
Novella, I. S. & Ebendick-Corpus, B. E. (2004). Molecular basis of fitness loss and fitness recovery in vesicular stomatitis virus. J Mol Biol 342, 1423–1430.[CrossRef][Medline]
Novella, I. S., Hershey, C. L., Escarmis, C., Domingo, E. & Holland, J. J. (1999). Lack of evolutionary stasis during alternating replication of an arbovirus in insect and mammalian cells. J Mol Biol 287, 459–465.[CrossRef][Medline]
Payne, A. F., Binduga-Gajewska, I., Kauffman, E. B. & Kramer, L. D. (2006). Quantitation of flaviviruses by fluorescent focus assay. J Virol Methods 134, 183–187.[CrossRef][Medline]
Rozas, J. & Rozas, R. (1999). DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis. Bioinformatics 15, 174–175.
Ruiz-Jarabo, C. M., Arias, A., Baranowski, E., Escarmis, C. & Domingo, E. (2000). Memory in viral quasispecies. J Virol 74, 3543–3547.
Ruiz-Jarabo, C. M., Arias, A., Molina-Paris, C., Briones, C., Baranowski, E., Escarmis, C. & Domingo, E. (2002). Duration and fitness dependence of quasispecies memory. J Mol Biol 315, 285–296.[CrossRef][Medline]
Sauder, C. J., Vandenburgh, K. M., Iskow, R. C., Malik, T., Carbone, K. M. & Rubin, S. A. (2006). Changes in mumps virus neurovirulence phenotype associated with quasispecies heterogeneity. Virology 350, 48–57.[CrossRef][Medline]
Scherret, J. H., MacKenzie, J. S., Khromykh, A. A. & Hall, R. A. (2001). Biological significance of glycosylation of the envelope protein of Kunjin virus. Ann N Y Acad Sci 951, 361–363.
Schneider, W. L. & Roossinck, M. J. (2001). Genetic diversity in RNA virus quasispecies is controlled by host–virus interactions. J Virol 75, 6566–6571.
Scott, T. W., Weaver, S. C. & Mallampalli, V. L. (1994). Evolution of mosquito-borne viruses. In The Evolutionary Biology of Viruses, pp. 293–324. Edited by S. S. Morse. New York: Raven Press Ltd.
Shi, P.--Y., Tilgner, M. & Lo, M. K. (2002). a Construction and characterization of subgenomic replicons of New York strain of West Nile virus. Virology 296, 219–233.[CrossRef][Medline]
Shi, P.-Y., Tilgner, M., Lo, M. K., Kent, K. A. & Bernard, K. A. (2002). b Infectious cDNA clone of the epidemic West Nile virus from New York City. J Virol 76, 5847–5856.
Twiddy, S. S. & Holmes, E. C. (2003). The extent of recombination in members of the genus Flavivirus. J Gen Virol 84, 429–440.
Vignuzzi, M., Stone, J. K., Arnold, J. J., Cameron, C. E. & Andino, R. (2006). Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population. Nature 439, 344–348.[CrossRef][Medline]
Weaver, S. C., Rico-Hesse, R. & Scott, T. W. (1992). Genetic diversity and slow rates of evolution in New World alphaviruses. Curr Top Microbiol Immunol 176, 99–117.[Medline]
Weaver, S. C., Brault, A. C., Kang, W. & Holland, J. J. (1999). Genetic and fitness changes accompanying adaptation of an arbovirus to vertebrate and invertebrate cells. J Virol 73, 4316–4326.
Weiner, A., Erickson, A. L., Kansopon, J., Crawford, K., Muchmore, E., Hughes, A. L., Houghton, M. & Walker, C. M. (1995). Persistent hepatitis C virus infection in a chimpanzee is associated with emergence of a Cytotoxic T Lymphocyte escape variant. Proc Natl Acad Sci U S A 92, 2755–2759.
Woelk, C. H. & Holmes, E. C. (2002). Reduced positive selection in vector-borne RNA viruses. Mol Biol Evol 19, 2333–2336.
Received 29 September 2006;
accepted 24 November 2006.
This article has been cited by other articles:
|
|
 |
|
  A. T. Ciota, A. O. Lovelace, Y. Jia, L. J. Davis, D. S. Young, and L. D. Kramer Characterization of mosquito-adapted West Nile virus J. Gen. Virol., July 1, 2008; 89(7): 1633 - 1642. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 L. L. Coffey, N. Vasilakis, A. C. Brault, A. M. Powers, F. Tripet, and S. C. Weaver Arbovirus evolution in vivo is constrained by host alternation PNAS, May 13, 2008; 105(19): 6970 - 6975. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
A. T. Ciota, A. O. Lovelace, S. A. Jones, A. Payne, and L. D. Kramer Adaptation of two flaviviruses results in differences in genetic heterogeneity and virus adaptability J. Gen. Virol., September 1, 2007; 88(9): 2398 - 2406. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||
), C6/36-passage 40 (
) and alternate-passages (
) WNV in C6/36 and DF-1 cell cultures. AP1 refers to WNV passaged 39 times in C6/36 cells followed by a single passage in DF-1 cells. The m.o.i. for all growth curves was 0.01 p.f.u. per cell. Results are presented as the means±SD of duplicate assays. (a) WNV growth in C6/36 cell culture. (b) WNV growth in DF-1 cell culture.
, DP10; 
, DP20) WNV in C6/36 and DF-1 cell culture. The m.o.i. for all growth curves was 0.01 p.f.u. per cell. Results are presented as the means±SD of duplicate assays. (a) WNV growth in C6/36 cell culture. (b) WNV growth in DF-1 cell culture.