Previous Article | Next Article ![]()
Journal of Virology, January 2004, p. 1032-1038, Vol. 78, No. 2
0022-538X/04/$08.00+0 DOI: 10.1128/JVI.78.2.1032-1038.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Dennis W. Trent, and Thomas P. Monath*
Acambis, Inc., Cambridge, Massachusetts 02139
Received 4 June 2003/ Accepted 3 October 2003
|
|
|---|
|
|
|---|
10 kb. Mutation rates vary between 10-3 and 10-6 in individual studies, depending on the virus and method used, although none of these methods provide an exact measure of the error rate. In vitro methods that utilize purified viral RNA polymerases may overestimate the error rate due to suboptimal in vitro enzymatic conditions and experimental complexity. In ex vivo methods, such as sequencing of phenotypic mutants generated in cell culture and analysis of monoclonal antibody escape mutants or of rates of reversion of phenotypic markers, the number of rounds of RNA replication in infected cells is not precisely known. Methods that involve a reverse transcription-PCR step with cloning of cDNA in bacteria prior to sequencing can be misleading because of the accumulation of additional mutations introduced by reverse transcription-PCR. The Flavivirus genus of the Flaviviridae family consists of about 70 viruses of which 40 have been associated with human illness (reviewed in reference 2). Flaviviruses are small enveloped plus-strand RNA viruses (reviewed in reference 18). The viral particle contains a nucleocapsid composed of viral RNA and capsid protein C, surrounded by a lipid bilayer in which the envelope protein E and membrane protein M are embedded. The genomic RNA is about 11,000 nt long and contains a single long open reading frame (ORF). Translation of the ORF produces a single polyprotein precursor containing viral proteins in the order: C-prM-E-NS1-NS2A-NS2B-NS3-NS4A-2K-NS4B-NS5, where C through E are the structural components of the virion (prM is a precursor for M). Nonstructural (NS) proteins NS1 through NS5 are required for virus replication in the cytoplasm of infected cells (18). Immature virions are produced by budding into the lumen of the endoplasmic reticulum (ER) and movement of virus to the cell surface through the exocytosis pathway. The structure of the flavivirus particle has been resolved by cryoelectron microscopy and fitting the known structure of the E protein (29) into the electron density map (15).
The ChimeriVax technology has been used to create attenuated flavivirus vaccines using the yellow fever virus (YF) 17D vaccine virus as a vector in which the prM-E envelope protein genes are replaced with those from a heterologous flavivirus (4, 10, 11, 24-26). To produce a chimeric virus, the viral RNA genome is first synthesized by in vitro transcription of an appropriately engineered DNA template with bacteriophage SP6 RNA polymerase. Following transfection of cells with the in vitro RNA transcripts, virus replication begins with genome amplification by the YF 17D RNA-dependent RNA polymerase (YFpol). In this study, yellow fever virus-dengue virus (ChimeriVax-DEN) type 1 to 4 chimeras were prepared by transfecting Vero cells with in vitro RNA transcripts. The resulting viruses were cloned by three consecutive plaque purifications. Both cloned and uncloned viruses were sequentially passaged to evaluate their genetic stability during prolonged growth in cell culture. Each virus was sequenced by the consensus method several times at different passages, including before and after plaque purification of the cloned variants. We observed that YFpol introduced very few mistakes in the genome during plaque purification steps, indicating a high fidelity of the polymerase. Characterization of genomes with YFpol mistakes and mutations that accumulated during multiple passages shed new light on the mechanism of flavivirus evolution, RNA replication, and particle assembly. Data from these studies challenge the widely held concept that RNA genomes accumulate substitution of up to 1 nt for each synthesized genomic RNA molecule.
|
|
|---|
![]() View larger version (18K): [in a new window] |
FIG. 1. Passages of the cloned and uncloned ChimeriVax-DEN22000 viruses. P1 virus was obtained by electroporation of Vero cells with in vitro RNA transcripts synthesized using SP6 RNA polymerase. Uncloned P2 virus was generated by one additional passage. Cloned PMS candidates (clones A and B) were produced by three consecutive plaque purifications starting from P1 virus followed by two additional passages. Genetic stabilities of these cloned and uncloned variants were examined in additional downstream virus passages. ChimeriVax-DEN22001, -DEN3, -DEN4, and -DEN1 viruses were passaged similarly, except that plaque purification passages (P3 to P5) were initiated from the uncloned P2 viruses.
|
5 to 10% cytopathic effect was observed. Virus-containing media were clarified by low-speed centrifugation, diluted with an equal volume of FBS, and stored at -70°C. Up to 10 clones were prepared for each chimera by three sequential plaque purifications. Plaque picks (agarose plugs) were placed in 0.4 ml of medium 199 (Gibco) containing 50% FBS, vortexed for several seconds, and frozen at -70°C. To sequence viruses obtained after the first plaque purification, 100 µl of the indicated plaque pick supernatant (virus eluted from agarose plug) was used to infect cells in T25 flasks, and the amplified viruses were harvested and sequenced. Otherwise, all samples of virus passages shown in Table 1 were sequenced without additional expansion. |
View this table: [in a new window] |
TABLE 1. Nucleotide and amino acid changes in uncloned and cloned ChimeriVax-DEN viruses
|
|
|
|---|
Consensus full-genome sequencing, which does not detect random mutations, of the uncloned P2 viruses revealed no unexpected mutations (Table 1). Mutations accumulated in the uncloned ChimeriVax-DEN1 to -DEN3 variants during their further passages to study genetic stability (see passages after P2 in Table 1).
Since no mutations were detected in the uncloned P2 ChimeriVax-DEN22000 virus, no mutations should be detected in the consensus sequence of P1 virus used to start plaque purifications. Nevertheless, each of the cloned A and B P6 ChimeriVax-DEN22000 viruses had a 1-nt change: G to T at nt 7734, resulting in an R-to-M amino acid change at NS5-29 located in a relatively conserved N-terminal region of NS5 and a silent A-to-T change at nt 6643, respectively (Table 1). Two additional genetic stability passages of clone A virus resulted in accumulation of subpopulations containing changes at nt 515 and 2135 detected at P8. Changes in clone A had potentially significant phenotypic implications. Therefore, further work with this clone was terminated. In clone B virus at P9, an additional G/C heterogeneity was observed at nt 517. The G-to-C nt change at this position resulted in an M-to-I amino acid change at prM-12 in the N-terminal part of prM that is cleaved during virus maturation. The effect of this mutation is similar to that of the nt 515 change in clone A, resulting in a similar mutation (M to L) at the same residue of prM. An M-to-L change at this residue was also seen in uncloned P15 ChimeriVax-DEN22001 virus (Table 1). This heterogeneity became a dominating change in virus at P14. An L-to-F change at amino acid 177 of NS4B protein was selected during virus passage to P14. The P18 of clone B had a consensus sequence identical to that of P14. It was significant that the P2 samples of both clones A and B harvested after the first plaque purification contained the same mutations as the P6 viruses at nt 7734 and 6643, respectively (Table 1).
Starting plaque purification with the uncloned P2 virus, rather than with P1 virus, immediately yielded acceptable cloned vaccine candidates of ChimeriVax-DEN22001 (clone A), -DEN3 (clone A), and -DEN4 (clone B; clone A was not sequenced) (Table 1). Only clone A of the ChimeriVax-DEN3 chimera had one silent nucleotide change at position 6607 detected at both P3 after the first plaque purification and P6 after the third plaque purification. In contrast, 10 clones of ChimeriVax-DEN1 were examined to find one without unacceptable mutations. Only clone J of the ChimeriVax-DEN1 chimera had no mutations at both P3 and P6. The entire genomes of clones A, B, C, and E were sequenced at P6, and the presence of mutations was checked at P3. Only the structural protein coding region of clone D was sequenced at both P3 and P6. In all cases, mutations present at P6 were present in P3 viruses (Table 1). (Fragment I of clones F to I was sequenced only at P3.)
Some clones of the ChimeriVax-DEN1 chimera had identical mutations, such as the E-251 V-to-F change in clones A, B, D, and G; the E-204 K-to-R change in clones E, F, and I; and the E-311 E-to-D and E-351 V-to-L changes in clones C and H. These changes were not detected in the consensus sequence of the starting uncloned P2 virus. Mutations E-204 and E-251 were first detected as heterogeneities in the uncloned P5 virus. The E-204 mutation became dominant in the uncloned P15 virus, and it also appeared in clone J virus by P10 (Table 1). The rate of accumulation of the E-204 mutation was studied in more detail in additional passaging experiments using clone J. These studies indicated this mutation was infrequently introduced by YFpol (C. Penal and F. Mitchell, data not shown). Since the E-204 mutation and all other mutations found in clones A to I at P3 appeared as clear mutations, they must have been present in individual particles of the uncloned P2 virus population before the first plaque purification. The fact that the E-204, E-251, E-311, and E-351 mutations were not detected in the consensus P2 sequence despite the high proportion of these mutations in the plaques picked was most likely due to the low sensitivity of consensus sequencing. Our prior experience with quantitation of mutant subpopulations using the MAPREC method (5) demonstrated that the level of detection of mutations by sequencing is approximately 10 to 15%. During the first plaque purification step, a few slightly larger plaques were picked together with smaller plaques which could account for some of the observed variation. The possibility that any of these mutations originated in plasmid templates is unlikely because of the high fidelity of bacterial DNA synthesis (7, 23). Ten individual plasmid clones from the ChimeriVax-DEN1 DNA template were prepared in Escherichia coli and sequenced. None of the clones contained the mutations detected in clones A to I of the virus or any other mutations (M. Parsons, data not shown).
Clone B of ChimeriVax-DEN22001 was prepared and sequenced at P3 after the first plaque purification, P6 after the third plaque purification, and P15. The P6 virus had two silent nucleotide changes at nt 1030 and 1789 which were not present at P3 (Table 1). These mutations occurred during plaque purification passages, and thus represented mistakes of YFpol. These were the only two mistakes introduced by YFpol during three plaque purifications of all the cloned chimeras described.
The cloned and uncloned ChimeriVax-DEN4 viruses did not accumulate mutations during genetic stability passages. In contrast, both cloned and uncloned variants of the other chimeras accumulated various mutations in genetic stability passages. Most of these mutations resulted in amino acid changes, indicating that virus adaptation is driven by forces operating at the protein function level. Some of the mutations observed among different chimeras coincided (Table 1). For instance, the same E-204 mutation accumulated in uncloned and clone B ChimeriVax-DEN22001 (detected at P15 of both variants), uncloned and clone J of ChimeriVax-DEN1 (at P5 and P10, respectively), and ChimeriVax-DEN3 (the E-202 mutation; see footnote e to Table 1). This suggests convergent mechanisms of adaptation of the ChimeriVax-DEN1 to -DEN3 chimeras to Vero cells. Many diverse mutations accumulated in uncloned and cloned versions of the same chimera or in different clones of a chimera (Table 1). The latter observation indicates that one virus can evolve in different directions under similar growth conditions. It is possible that several different mutations can improve virus replication and that further virus evolution depends on which mutation occurs first.
|
|
|---|
|
View this table: [in a new window] |
TABLE 2. Assigned classification of mutations detected during preparation of ChimeriVax-DEN viruses and the genetic stability passages shown in
|
Some mutations in ChimeriVax-DEN1 clones A to I detected after the first and the third plaque purification (Table 1) did not clearly fit the classification in Table 2. The E-204 and E-251 mutations were detected in more than one clone and also accumulated in genetic stability passages of the uncloned virus. The E-204 mutation also accumulated during continuous passage of ChimeriVax-DEN1 clone J. Thus, these mutations are beneficial for the ChimeriVax-DEN1 chimera replication. Since mutations are introduced by YFpol infrequently, these changes were likely mistakes of SP6 polymerase amplified prior to the first plaque purification.
Of the 17 nucleotide changes considered to be mistakes of SP6 polymerase, 5 were transitions and 12 were transversions. Thirteen of the mutations were changes from A, G, or C nucleotides to T, with a G-to-T change being the most frequent (nine mutations). This suggests that a misincorporation of T was the most common mistake during in vitro transcription. Of the 17 nucleotide changes introduced by YFpol, including the two mistakes that occurred during plaque purifications plus beneficial mutations that accumulated during genetic stability passages of the clones, 11 changes were transitions and 6 were transversions. The most frequent changes were A to G (four mutations) and C to T (five mutations), which theoretically should be the two most common mistakes during late plus-strand RNA synthesis due to T · G and G · T mispairing, respectively. This observation may point toward the existence of distinct mechanisms for synthesis of double-stranded RNA (with high fidelity), representing most rounds of amplification in infected cells, and late transcription (with low fidelity). The majority of rounds of viral RNA synthesis should occur during the latent phase of infection accompanied by exponential duplication of the double-stranded replicative-form (RF) RNA. During the late linear phase, RNA synthesis is mostly single-round transcription on existing minus-strand template that is stable for at least 20 h and is renewed slowly (reviewed in references 1a and 30).
|
|
|---|
11,000 nt long), or 1.34 x 10-4 per nucleotide copied, which is comparable to the reported error rate of T7 RNA polymerase (0.5 x 10-4 [13]). However, it could be somewhat lower because all of the mutations may not be due to mistakes of SP6 polymerase and the accumulation of mutations in ChimeriVax-DEN1 may have occurred prior to the first plaque purification. On the other hand, if plaque purifications of ChimeriVax-DEN22001, -DEN3, and -DEN4 had been done at P1, rather than P2, more mistakes of SP6 polymerase may have been detected.
Two nucleotide changes were present in a total of 10.5 full genome equivalents that occurred during plaque purifications and therefore were clearly introduced by YFpol (Table 1). To estimate the YFpol error rate, we need to know how many times a genomic RNA molecule is copied during one act of plaque formation. Chimeras grow to titers of
107 PFU/ml in Vero cells; therefore, a plaque pick of roughly 100 µl in volume will contain
106 infectious particles. We assume that 100 to 1,000 times more plus strands are synthesized (108 to 109 molecules) and are not accounted for in infectious virus because of the loss of infectivity of mature secreted particles, as well as the accumulation of immature viral particles and possible RNA turnover inside the cells (reviewed in references 1a and 30). If most rounds of RNA amplification are semiconservative duplications of the RF RNA, as many as 25 to 30 rounds of synthesis can take place (e.g., 109 molecules
230, thus requiring 30 duplications) during one plaque formation, in
103 to 104 participating cells. In support of our estimate of the total number of plus strands synthesized, 20 to 40 µg of protein E per 107 cells is produced following infection with tick-borne encephalitis virus (1, 14a). Since each viral particle contains 180 protein E monomers (15), up to 2 x 105 viral particles (equal to secreted plus strands) are produced by each cell. This number does not include all plus strands synthesized in infected cell.
By dividing the number of YFpol mistakes (two mistakes) by the number of sequenced full genome equivalents (10.5 genomes), the estimated number of rounds of RNA synthesis during one plaque formation (25 to 30 rounds), and the number of plaque purification steps for each clone (three steps), the resulting error rate of YFpol is 2/[10.5 x (25 to 30) x 3] = 0.0021 - 0.0025 misincorporation per genomic RNA molecule (11,000 nt long) per synthesis, or 1.9 x 10-7 to 2.3 x 10-7 per copied nucleotide. With our system, we can detect only nonlethal mutations. We do not believe, however, that lethal mutations could constitute a percentage significant enough to affect the estimate, considering the overall low occurrence of mutations during plaque purifications. It is known that the distribution of mutations in RNA viruses is shifted toward neutral or slightly deleterious mutations rather than strongly deleterious or lethal mutations that are quickly eliminated by competition during virus replication; besides, a small fraction of mutations is beneficial (8, 9). To conclude that there could be a large proportion of undetected lethal mutations, we must observe many synonymous mutations, such as those that do not change the amino acid encoded, which is not the case.
Our data indicate that the mutation rate of YFpol is on the low end of the spectrum among RNA polymerases. This observation is in accord with previous studies demonstrating the high genetic stability of YF (12, 27, 33) and makes YF 17D an excellent vector for the construction of recombinant vaccines. It is believed that the evolution of YF and other flaviviruses in nature is constrained by the need of alternation of vertebrate and arthropod hosts. Our data extend this hypothesis further, suggesting that in the course of flavivirus evolution, variants with high fidelity of their RNA polymerases were selected to facilitate the necessary continuous host alternation.
|
|
|---|
![]() View larger version (12K): [in a new window] |
FIG. 2. Clustering of beneficial amino acid changes in the four ChimeriVax-DEN chimeras that accumulated in all cloned and uncloned viruses during genetic stability passages. The locations of mutations in the polyprotein and genomic RNA in each of the chimeras were plotted according to their positions in the genome.
|
If the latter hypothesis is correct, it is interesting that no mutations accumulated in the YF-specific capsid protein C. The mutations in the N-terminal part of prM could have functional significance in light of recent data indicating that prM functions as a chaperone necessary for proper folding of E (20). These mutations could facilitate a more appropriate folding of E such that the envelope more precisely conforms to the architecture of the nucleocapsid. All mutations in the cluster located in the middle of E protein are on the outward surface of the protein in the dimerization domain II (29). The effects of these changes on the structure and function of the E proteins of chimeric viruses need further investigation.
Most intriguing is the cluster of mutations in NS4B (Fig. 2), because it suggests that this protein plays a role in virus assembly. Flavivirus replication induces proliferation of smooth ER membranes and the accumulation of characteristic membranous structures described as convoluted membranes (CM), paracrystalline arrays (PC), and vesicle packets (VP) (21, 31, 32). Viral double-stranded RNA together with NS1, NS3, NS5, NS2A, and NS4A proteins colocalized with VP structures, indicating that they are the main sites of RNA synthesis, while CM and PC structures were proposed to be main sites of proteolytic cleavage (14, 21, 32). NS4B appears to have a distinct distribution, associated with ER membranes throughout the cytoplasm (31). Recently, NS2A was implicated in virus assembly, because specific mutations in this protein abolished production of infectious particles (16). In addition, it was proposed that NS3 also plays a role in assembly (14, 16, 19). Since NS2A and NS3 are both located in the sites of RNA synthesis and taking into account that both these proteins bind viral genomic RNA (21), it is tempting to speculate that NS2A and NS3 are responsible for genome encapsidation, while NS4B comes into play during the next step, budding, possibly by facilitating specific interaction of the nucleocapsid and envelope. NS4B could also facilitate particle formation indirectly, via its involvement in the formation of virus-induced membrane structures essential for the secretion pathways (14b, 20a).
Finally, the suggestive nature of these sequencing data should be recognized. The evolutionary, enzymological, and molecular aspects discussed above require further investigation. The process of accumulation of viral RNAs in infected cells needs to be precisely quantitated using modern approaches such as real-time PCR. This could establish the exact number of rounds of RNA synthesis, confirm or correct our estimate of YFpol fidelity, and add to a better understanding of the molecular mechanisms of flavivirus RNA synthesis. The high fidelity of YFpol requires confirmation by other available methods, including using wild-type YF. The hypothetical role of NS4B (as well as the elements in prM-E) in particle assembly requires direct demonstration.
Present address: DynPort Vaccine Co., LLC, Frederick, MD 21702. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»