Previous Article | Next Article ![]()
Journal of Virology, February 2007, p. 1879-1887, Vol. 81, No. 4
0022-538X/07/$08.00+0 doi:10.1128/JVI.01899-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
Centro de Biología Molecular Severo Ochoa (CSIC-UAM), Universidad Autónoma de Madrid, Cantoblanco, 28049 Madrid, Spain
Received 31 August 2006/ Accepted 23 November 2006
|
|
|---|
|
|
|---|
The rapid replication, high mutation rates, and large population sizes of RNA viruses lead to genetically highly heterogeneous populations, termed viral quasispecies (for some reviews, see references 17, 18, 19, and 20). Viral quasispecies have an enormous potential for rapid evolution through adaptation to new environments. As quasispecies contain not only many different single mutants but also many different multiple mutants, evolution through the simultaneous or sequentially rapid fixation of combinations of mutations, including compensatory mutations, may occur very frequently in RNA viruses. Indeed, many documented cases of compensatory mutations (see references mentioned above) involve proteins from RNA viruses.
The establishment and maintenance of a persistent infection of mammalian cells by a picornavirus, foot-and-mouth disease virus (FMDV), involved the coevolution of cells and virus (13, 14). The determinant element in the initiation of the persistent infection was the selection of cell variants with a slightly increased resistance to the FMDV clone (C-S8c1) used to initiate the infection (34). The surviving cells became increasingly resistant to infection, while the virus shed evolved to become progressively more virulent to the original cells and acquired, in addition, other phenotypic traits. These traits included a smaller plaque size, temperature sensitivity, and a reduced stability of the virion during purification (13, 16, 49). Nucleotide sequencing revealed that several mutations had been fixed in the C-S8c1 genome during the persistent infection (16, 54). In particular, nine amino acid substitutions were identified in the capsid of the virus (R100) rescued after 100 passages of persistently infected cells (16). Most of these substitutions had not occurred in previously sequenced FMDV variants, which were generally isolated from animals acutely infected in the field or from cytolytic infections in the laboratory, or even in related picornaviruses. These results led those authors to conclude that if provided with a new suitable environment, like residence in carrier cells, even positions in the viral genome that are regarded as highly invariant can be subjected to rapid change (16).
A major interest of our group is the study of the molecular determinants of assembly and stability of viruses. In a previous analysis, we provided evidence that in FMDV, most residues at the capsid interpentamer interfaces, and their interactions, are important for the infectivity, and stability, of the virion (35). In addition, a few mutations outside these interfaces have been shown to have an effect on the conformational stability of picornaviruses, including FMDV (2, 24, 50, 55). Four (nearly half) of the mutations in the capsid of FMDV R100 involve identical clusters of residues that surround the pores at each capsid fivefold symmetry axis (16) (Fig. 1) and are also a part of the interfaces between the capsid protomers. The low stability of R100 during purification, and its temperature sensitivity in cytolytic infections, led us to investigate whether this cluster of mutations in the R100 capsid and, specifically, the loss of interprotomer disulfide bonds and salt bridges around those pores (because of C3007V and D3009A mutations) could be responsible for a low conformational stability. We were also interested in determining whether these clustered mutations could act through nonadditive ("epistatic") and perhaps compensatory effects. The mutations have been introduced into an infectious FMDV clone, and their effects on infectivity and particle stability and on the biological fitness of the virus in competition experiments have been analyzed. We have found unexpected, complex mutation dynamics that involved the repeated fixation in a cytolytic infection context of some of the mutations previously found in R100, a virus that emerged in a persistent infection context. In addition, combinations of these mutations led to complex compensatory effects and an increase in the biological fitness of cytolytic infections.
![]() View larger version (28K): [in a new window] |
FIG. 1. Location in the virus capsid of a cluster of residues found to be mutated in FMDV R100. (A) C tracing model of a pentameric subunit in the FMDV C-S8c1 capsid viewed from the exterior of the virion. Five copies of each of the capsid proteins VP1, VP2, VP3, and VP4 are depicted in different colors. The residues around the fivefold axis pore (center) that were found to be mutated in FMDV R100 (16) are represented as space-filling models and are color coded (green, Cys3007; red, Asp3009; yellow, Asn3013; magenta, Met3014). (B) Close-up view of the N termini of five symmetry-related copies of VP3 around the fivefold axis pore in a pentamer shown in the same orientation depicted in A. The tracing of the main-chain backbone is indicated by ribbons. The residues that were found to be mutated in FMDV R100 are represented as in A. The five cysteine residues (green) constrict the pore (center). (C). Side view of the fivefold axis region in a pentamer. The outer surface is up. The ß-annulus formed by ß-strands belonging to the five VP3 subunits (colored arrows) is shown below the residues found mutated in FMDV R100 (shown as colored space-filling models as in A and B).
|
|
|
|---|
2292 corresponds to pO1K/C-S8c1 with a segment of 2,292 bp (between the AflII and AvrII restriction sites) deleted. This segment includes part of the L protease and the entire capsid region coding for the capsid proteins, except for a short C-terminal segment of VP1 (identical in sequence between C-S8c1 and O1K).
Site-directed mutagenesis, subcloning, and DNA sequencing.
Site-directed mutagenesis of FMDV capsid residues was carried out on plasmid p3242/C-S8c1 by an inverse PCR method using the QuikChange site-directed mutagenesis kit (Stratagene). The mutations introduced were confirmed by automated sequencing using dideoxynucleotide terminators. The sequences were obtained by the Genomic Unit of the Madrid Science Park at the Universidad Autónoma de Madrid. The ABI PRISM BigDye Terminator kit and an ABI PRISM 3730 or ABI PRISM 3700 automated sequencer (Applied Biosystems) were used. A region of about 600 nucleotides around the introduced mutations was sequenced in each case. For some mutants, indicated in Results, the entire sequence of the 2,292-bp segment to be subcloned was determined. No second-site mutations were found in the mutated plasmid. The mutagenized 2,292-bp segments were subcloned into pO1K/
2292 to obtain pO1K/C-S8c1 plasmids, and the presence of the introduced mutation was again checked by sequencing.
Transcription of viral RNA and electroporation of eukaryotic cells. Transcription and transfection of infectious viral RNA have been previously described (35). In each experiment, the same amount of mutated and nonmutated FMDV RNA was used. A negative control containing no viral RNA was also included.
Titration, amplification of viruses, and extraction of viral RNA. Virus titers were always determined at least in duplicate in standard plaque assays. When necessary, mutant virions obtained were amplified by a single passage in BHK-21 cell monolayers at the highest possible multiplicity of infection (MOI). Extraction of viral RNA was performed by using TRIzol (Invitrogen) essentially as described previously (35).
Competition assays. Confluent BHK-21 cell monolayers cultured in B25 flasks (3 x 106 cells) were infected with a mixture containing equal amounts of the competing viruses at an MOI of 0.1. Adsorption of viruses was allowed to proceed for 60 min at 37°C. Fresh Dulbecco's modified Eagle medium (DMEM) supplemented with 2% fetal calf serum was then added. Viruses were harvested when a complete cytopathic effect was observed, clarified by centrifugation, and diluted 10 times with fresh DMEM. A fraction of the supernatant was stored at 70°C and used for RNA extraction. Another fraction of the supernatant was used to infect new cell monolayers. This cycle of infection and harvest was repeated up to seven times. For each pair of competing viruses, the experiment was independently repeated two or three times. When eight passages had been completed, the RNA of the viral populations obtained at different passages was extracted, reverse transcribed and amplified by reverse transcription-PCR, and sequenced as described above to estimate the degree of dominance of each of the two competing viruses at the passage considered and whether other mutations had been fixed in a detectable proportion of the viral RNA molecules during replication. In general, the entire capsid-coding region was sequenced. In some cases, as indicated in Results, the entire genome was sequenced. The approximate proportions of the competing genomes and of any other mutation that happened to arise during replication were estimated by densitometry of the chromatogram peaks corresponding to the nucleotide present in the mutant virus and the nucleotide present in the nonmutated virus. The ratio of the integrated areas under each peak was used to estimate the ratio of mutant genomes to nonmutant genomes at the position considered.
Kinetic assays of virus production. For each assay, two confluent BHK-21 cell monolayers in P60 petri dishes (2.5 x 106 cells) were infected with the chosen mutant FMDV or with FMDV C-S8c1 at an MOI of 10. The progression of the cytopathic effect was compared at different times by examination under a microscope. The kinetics of extracellular viral production were determined by the titration of aliquots taken at different times postinfection. Production assays were always carried out in duplicate. For each assay, titers were also always calculated in duplicate.
Thermal inactivation assays. Two types of assays were used. In the first type, 1-ml aliquots of cell culture supernatants obtained from cytolytic infections or after viral RNA transfection of cell monolayers containing viral populations with an approximate titer of 2 x 107 to 4 x 107 PFU/ml were clarified and subsequently heated for different lengths of time at 42°C. The percentage of infectious particles remaining after heating was determined by use of standard plaque assays, and the values obtained at different times were fitted to a first-order exponential decay model by using the program Kaleidagraph (Abelbeck Software) to obtain the inactivation rate constant (35). Alternatively, the logarithmic values of the viral titers were linearly fitted using the same program. In a second type of assay, the same cell culture supernatants containing virus populations were diluted in DMEM plus 2% fetal calf serum to a final concentration of about 1,000 PFU/ml. Aliquots (100 µl) of the diluted samples were incubated for different lengths of time at 42°C, and the percentage of infectious particles remaining after heating was determined by plaque assays as described above.
|
|
|---|
Analyses of 211 FMDV sequences currently available from the GenBank database confirmed a previously reported observation (16) that Asp3009 is highly conserved (98%) in FMDV from acute/lytic infections and that Cys3007 is also remarkably conserved (87%). To determine whether mutations of residues Cys3007 or Asp3009 and the respective removal of disulfide bonds or buried salt bridges had an effect on the infectivity and/or stability of FMDV, we introduced them individually into the capsid of C-S8c1, the virus used to establish a persistent infection. Viral RNAs carrying the C3007V or D3009A mutation and nonmutated viral RNA were obtained from an infectious DNA clone with the corresponding mutation and were used to electroporate susceptible cells. The progeny virus populations obtained were collected and titrated (Table 1). Relative to the nonmutated control, no significant differences in viral titer at 45 h or 55 h postelectroporation, or in plaque size, were observed for the C3007V mutant in repeated electroporation assays. In contrast, for the D3009A mutant, nine independent electroporations yielded widely different virus titers at 45 h postelectroporation, with values that ranged from undetectable to a titer close to that of the nonmutated control. On most occasions, plaque size heterogeneity was also observed. In one experiment, D3009A populations were recovered at 45 h, 55 h, and 95 h postelectroporation, and a gradual increase in titer, up to nearly wild-type levels, was observed (Table 1). These data suggested that the D3009A mutation had a drastic negative effect on infectivity and that reversion or new mutations had occasionally occurred during the replication of D3009A, leading to heterogeneous populations with variable amounts of virus genotypes that differed drastically in infectivity.
|
View this table: [in a new window] |
TABLE 1. Values of viral titer and plaque size for a set of viruses with mutations around the fivefold axis of the FMDV capsid
|
To summarize the above-described results, of the mutations that were fixed around the fivefold axes in the capsid of FMDV R100 during a persistent infection, C3007V and M3014L were individually tolerated without any substantial effect on the viral yield in cytolytic infections. In contrast, D3009A had a drastic negative effect on viability but strongly favored the fixation of the M3014L mutation at a neighboring position, leading to a compensatory effect that allowed the recovery of infectious virions in high yields. In turn, the D3009A/M3014L double mutation facilitated the partial fixation of the N3013H mutation, again at a neighboring position.
Effect of mutations around the fivefold axes of the R100 capsid on the thermostability of the FMDV virion. The thermostabilities of the viable C3007V mutants were determined in heat inactivation assays (Fig. 2). The inactivation rate constants at 42°C (average ± standard deviation of two or three determinations from independent experiments) were, respectively, 0.023 ± 0.010 min1, 0.013 ± 0.001 min1, and 0.023 ± 0.04 min1. By comparison, the inactivation rate constant of the nonmutated virus (obtained from 55 independent determinations as a measure of the reproducibility of the assay) was 0.019 ± 0.006 min1. In the case of the C3007V mutant and the nonmutated control, inactivation analyses of diluted or undiluted virus were carried out and yielded similar results (Fig. 2A and B, respectively). No substantial differences in thermostabilities were found between any of those three mutants and the nonmutated infectious virion (Fig. 2B). These results reveal that contrary to what we suspected, Cys3007 and Met3014 alone or in combination with Asp3009 are not involved in the stabilization of the FMDV particle against thermal inactivation.
![]() View larger version (15K): [in a new window] |
FIG. 2. Thermal inactivation kinetics of the parental and mutant FMDVs. (A) Inactivation of undiluted parental virus (circles) and mutant C3007V (squares) (using the first type of inactivation assay described in Materials and Methods) at 42°C. Linear fitting of the logarithmic values is indicated. (B) Inactivation of the diluted parental virus (circles), C3007V (squares), M3014L (triangles), and D3009A/M3014L (inverted triangles) (using the second type of inactivation assay described in Materials and Methods) at 42°C. Fitting of the data to a single exponential is indicated.
|
In the case of the M3014L mutant against the nonmutated virus, serial passaging of the virus mixture led to a gradual increase in the nonmutated genotype that was highly dominant from passage 6 onwards in two independent competition experiments (Fig. 3A). Sequencing of the entire capsid region of the resulting populations (passage 8) showed no mutation present in a detectable proportion at any position.
![]() View larger version (22K): [in a new window] |
FIG. 3. Approximate percentage of FMDV genomes carrying a specific residue at a defined position during serial passages in competition experiments. The global sequence of the viral populations obtained after each passage was determined by automated sequencing, and the relative abundance of each residue type at the positions of interest was estimated from the ratio between the integrated areas of the corresponding bands in the chromatograms, as described in Materials and Methods. (A) Competition between the M3014L mutant and the parental virus. Circles, Met; squares, Leu at position 3014. (B) Competition between the D3009A/M3014L mutant and the parental virus. (B.1) Circles, Ala; squares, Asp at position 3009; black inverted triangles, Leu; black triangles, Met at position 3014. (B.2) Squares, Thr; circles, Pro at position 1104. (C) competition between the C3007V mutant and the parental virus. (C.1) Squares, Val; circles, Cys at position 3007. (C.2) Squares, Ala; circles, Asp at position 3009. All three competition assays were carried out in duplicate or triplicate with very similar results.
|
The above-mentioned observations indicate that in cytolytic infections, the M3014L mutation is biologically disadvantageous and that D3009A causes strongly reduced viability when introduced in isolation; however, the introduction of D3009A leads to the fixation of M3014L. The two disadvantageous mutations together could yield a near neutral genotype relative to the nonmutated genotype, because the relative abundance of both viruses was not significantly altered during the first passages in competition assays. The final takeover by the mutant could have been caused exclusively by the fixation of an additional mutation, P1104T. However, the takeover and the fixation of this mutation were not concomitant (see above). We consider it more likely that the 3009/3014 double mutation conferred a slight advantage that may have been potentiated by the effect of the P1104T mutation.
In the case of the C3007V mutant against the nonmutated virus, contrary to what could be expected from the high conservation of Cys3007 in viruses from acute infections, the nonmutated virus did not reach dominance. Instead, in three independent competition experiments, the proportion of viruses carrying the C3007V mutation increased from passage 3 upwards, and in the last passages, they had reached full dominance (Fig. 3C.1). Remarkably, in all three competitions, the D3009A mutation became gradually fixed in the population together with C3007V (Fig. 3C.2). The entire genome of the progeny virus at passage 8 of the three series was sequenced, and no other mutations in addition to C3007V and D3009A were found. Thus, the presence of Val instead of Cys at position 3007 is biologically not advantageous in cytolytic infections but allows the fixation of Ala instead of Asp at position 3009, leading to the dominance of the double mutant.
Finally, we investigated the functional reason why those virus variants with compensatory mutations, specifically C3007V/D3009A, may show increased biological fitness. The viral progeny yields obtained upon infection with this mutant, the C3007V single mutant, or the nonmutated virus were similar, as were their plaque sizes. However, infection kinetic assays showed that the double mutant produced extracellular progeny viruses significantly faster, reaching the maximum amount in about 8 h, while the C3007V mutant and the nonmutated control reached similar titers but only after 24 h (Fig. 4). This suggests that the increased fitness of the C3007V/D3009A mutant is due, at least in part, to a faster multiplication rate.
![]() View larger version (27K): [in a new window] |
FIG. 4. Kinetics of extracellular virus production of nonmutated and mutant FMDVs. (A and B) Absolute virus titers at different times after infection. (A) Circles, nonmutated virus; squares, C3007V mutant. (B) Circles, nonmutated virus; squares, C3007V/D3009A mutant. (C and D) Relative infectivity of mutant viruses with respect to the infectivity of the nonmutated virus at different times after infection. The infectivity of the nonmutated virus at each time after infection has been given the reference value of 100. (C) White bars, nonmutated virus; black bars, C3007V mutant. (D) White bars, nonmutated virus; black bars, C3007V/D3009A mutant. Standard deviations are indicated in each case.
|
|
|
|---|
Unlike D3009A, other mutations fixed around the capsid pores in R100 (C3007V, N3013H, and M3014L) may have not occurred as a necessary consequence of a modified environment in persistently infected cells. These mutations were not fixed in other virus populations rescued from cells that had been persistently infected (22, 54). In addition, sequence analysis showed that Val or Leu was present at position 3007 or 3014 in as many as 11% or 54% of the more than 200 FMDV isolates from acute/lytic infections, respectively. His at position 3013 was not observed in viruses not related to persistent infections, but Asn at this latter position was not conserved either, as it was present in only 19% of the isolates. More importantly, in the present study, we showed that some of these mutations in R100 also became repeatedly fixed during cytolytic infections in response to the introduction of a previous mutation, as discussed below.
Linked mutational events may occur in response to the introduction of other mutations in the FMDV capsid. Our results also suggest that D3009A, irrespective of its role in a persistent infection, may be the key mutation for a series of linked mutational events in the FMDV capsid that involved a cluster of residues around the fivefold axis pores. The C3007V mutant outcompeted C-S8c1 but only through the fixation of the D3009A mutation. The D3009A mutant obtained by site-directed mutagenesis (mimicking the possible natural fixation of this mutation during a persistent infection) yielded high titers only through the fixation of the M3014L mutation. The introduction of the D3009A/M3014L double mutation facilitated the fixation of the N3013H mutation. Finally, the D3009A/M3014L mutation led, in competition experiments, to the fixation of a new mutation, P1104T, in a residue located close to the other mutated residues. In addition, every FMDV isolate analyzed that carried the D3009A mutation also carried an accompanying mutation of another residue that is spatially very close to residue 3009. These second mutations included C3007V or C3007A, S3008T or S3008A, N3013H, M3014L, and V3015Q in VP3 and H1108R, G1110R, and V1112L in VP1. Four of these mutations (S3008T, N3013H, H1108R, and G1110R) were very rarely observed in isolates that did not have the D3009A mutation, suggesting a linkage with D3009A. The introduction of a mutation like D3009A in a specific region of the virus capsid (either as a natural consequence of a modified environment, for example, during persistence, or by force using site-directed mutagenesis) may facilitate, or even require, the fixation of other mutations at neighboring locations.
Mutually compensatory effects of individually disadvantageous mutations may occur in the FMDV capsid. The results have shown that the D3009A and M3014L mutations are, per se, biologically disadvantageous and that C3007V may be slightly disadvantageous or, at least, does not provide a selective advantage. However, the simultaneous presence of D3009A and M3014L or of C3007V and D3009A led to an increase of the biological fitness relative to that of the nonmutated virus. The possibility of a putative substitution outside the capsid region in the viral RNA cannot be completely excluded in some instances. However, we consider this possibility to be unlikely. First, the same substitutions in the capsid repeatedly occurred in response to the same previous mutations; second, many compensating molecular events involve spatially close residues, as observed in our study; finally, in one case, the entire mutant genome was sequenced, and no other mutation was found. These observations provide clear evidence of second-site mutations that, even though intrinsically disadvantageous, may not only compensate for the detrimental effects of the primary mutations but also lead to the acquisition of new phenotypic characteristics that may provide a selective advantage. These effects open new possibilities for virus evolution.
An understanding of the structural basis of these compensatory effects must await further studies. However, note that the mutations fixed in response to other mutations (D3009A in response to C3007V, M3014L in response to D3009A, or N3013H in response to D3009A/M3014L) involved residues that are spatially in contact with or very close to each other, form a part of the interprotomer interfaces, and closely surround the pores at the capsid fivefold axes (Fig. 1). A likely possibility, which awaits testing, is that some of these compensation effects could be related to the kinetics of virion assembly. They could be also related to translocation events that might occur through permanent or transient capsid openings. Such translocation events have been detected in icosahedral viruses including picornaviruses. For example, in the minute virus of mice, a parvovirus, mutations of residues surrounding either the fivefold axis pores or buried cavities located close to the pores were shown to affect a conformational rearrangement associated with the translocation of the VP2 N terminus through the pores, an event that is needed for infectivity (9, 47). An impaired or untimely translocation event was also proposed to explain these effects (9). In coxsackievirus A9, a picornavirus, mutation of some residues at interprotomer interfaces located close to the capsid pores could result in the untimely or hindered externalization of the VP1 N terminus together with VP4 (2). In poliovirus, recent results indicate that the externalization of the myristoylated capsid protein VP4 and the N terminus of VP1, which are required for infectivity, takes place not through the fivefold axis pores but through transient openings at the interface between VP1 and VP3 that occur after a conformational change of the capsid (8). Further studies are required to determine whether the complex effect of mutations like C3007V, D3009A, and M3014L in the capsid of FMDV R100 could be related to a putative conformational rearrangement of the capsid and a translocation event, through either the pores or other transient openings, on the virion (under study).
To conclude, from the point of view of viral genetics, this study (i) indicates that the clustering of mutations at conserved residues around the fivefold axes in the capsid of an FMDV variant rescued from a persistent infection may also occur in cytolytic infections in response to changes in the local structural context, rather than being a direct consequence of adaptation to a different environment, and (ii) provides evidence for a chain of deterministic, mutually compensatory effects of individually disadvantageous or neutral mutations in the FMDV capsid that also led to the acquisition of new phenotypic traits.
This work was supported by grants from the Spanish Ministerio de Ciencia y Tecnología (BIO2003-04445 and BIO2006-00793) and Comunidad Autónoma de Madrid to M.G.M. and by an institutional grant from the Fundación Ramón Areces.
Published ahead of print on 6 December 2006. ![]()
|
|
|---|
X174. Genetics 170:989-999.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»