Previous Article | Next Article ![]()
Journal of Virology, August 2005, p. 10451-10459, Vol. 79, No. 16
0022-538X/05/$08.00+0 doi:10.1128/JVI.79.16.10451-10459.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Gema Gómez-Mariano,1
Pedro R. Lowenstein,2 and
Esteban Domingo1*
Centro de Biología Molecular "Severo Ochoa" (CSIC-UAM), Cantoblanco, E-28049 Madrid, Spain,1 Department of Medicine, Department of Molecular and Medical Pharmacology, and Gene Therapeutics Research Institute, Cedars-Sinai Medical Center, and Department of Medicine, University of California, Los Angeles, California2
Received 28 February 2005/ Accepted 9 May 2005
|
|
|---|
|
|
|---|
The prototypic arenavirus lymphocytic choriomeningitis virus (LCMV) (12, 63, 64) showed extreme sensitivity to 5-fluorouracil (FU)-induced mutagenesis (41, 72) compared with other RNA viruses subjected to comparable doses of the mutagen (79). The extreme sensitivity of LCMV to FU offered an opportunity to analyze the capacity of LCMV to regain infectivity following FU mutagenesis as well as the modification of genomic nucleotide sequence variations as the virus moves toward or away from the error threshold and to explore the possible dominance of hypermutated genomes in the transition to extinction. The results reveal a remarkable capacity of LCMV populations to modify their fitness level while maintaining an invariant consensus sequence. Multiple molecular clones were analyzed to define a sequence at nucleotides 470 to 550 within the intergenic region of genomic segment L. A number of standard and mutated oligonucleotide primers have failed to produce evidence of hypermutated viral sequences in the L polymerase gene. The results suggest that limited mutagenesis may be sufficient to drive LCMV close to the error threshold and that LCMV can rapidly regain fitness in the absence of a mutagen. The invariance of the consensus genomic sequence, associated with the movements toward and away from the error threshold, suggests a decisive participation of the mutant spectrum in determining infectivity levels and supports a lethal defection model for virus extinction through enhanced mutagenesis.
|
|
|---|
Virus infections.
Semiconfluent (2.8 x 106 cells in 100-mm-diameter dishes [Falcon]) monolayers of BHK-21 were infected with 0.01 PFU of LCMV ARM 53b per cell in 5 ml of DMEM supplemented with 10% FCS, 10% tryptose phosphate broth, 2% L-glutamine, 0.52% glucose, and 50 µg/ml gentamicin. Viral supernatants were harvested 48 h postinfection, clarified by centrifugation at 2,500 rpm for 30 min at 4°C, and stored at 80°C. Supernatants were titrated on Vero cell monolayers under a semisolid agar medium, as described elsewhere (3). Viral yields usually were
108 PFU/ml.
Serial passages of LCMV ARM 53b were carried out by infecting control (D) or FU-treated (F) cell monolayers with the corresponding supernatants of the previous infection. For infections carried out with a virus obtained following a D passage, 2.8 x 106 BHK-21 cells were infected with about 3 x 104 PFU of LCMV, whereas for infections with a virus derived from an F passage virus, the input was
103 PFU of LCMV. Virus supernatants and cells were harvested 48 h postinfection. To obtain preextinction virus populations, semiconfluent BHK-21 monolayers (15 182-cm2 flasks with 1.6 x 107 cells per flask) were infected at an MOI of 0.01 PFU/cell and incubated for 48 h in 10 ml of infection medium (with 100 µg/ml of FU) per flask. No specific LCMV sequences could be amplified and no infectivity (<10 PFU/ml) obtained from cell culture supernatants after a second passage in 100 µg/ml FU, even after 25-fold concentration with a Centricon Plus-20 Biomax-100 (Millipore) filter. Virus titers are means of three titrations; standard deviations never exceeded 10% of the mean.
Drug treatment. Preparation of FU stock solutions, BHK-21 cell viability assays, and infections in the presence of the mutagenic analogue were performed as previously described (41, 79). LCMV was considered extinguished when no reverse transcriptase PCR (RT-PCR)-amplifiable material and no infectivity (<10 PFU/ml) could be detected after three blind passages of the undiluted viral population in BHK-21 cells in standard culture medium in the absence of mutagens.
RNA extraction, RT-PCR, nucleotide sequencing, molecular cloning, and calculation of mutant-spectrum complexity.
RNA was extracted from supernatants of infected cell cultures by using Trizol (Invitrogen) according to the manufacturer's protocol. RNAs were amplified by two-step RT-PCR using ThermoScript reverse transcriptase (Invitrogen) and a reverse (antisense) primer (http://www.cbm.uam.es/mkfactory.esdomain/webs/CBMSO/plt_LineasInvestigacion.aspx?IdObjeto=17) at 60°C for 45 min, followed by PCR with Pfu DNA polymerase (Promega). As a control to confirm that an excess of template was being amplified, 1/10 of the initial amount of template RNA was processed in parallel and yielded a positive amplification band (4, 41). cDNAs either were purified with a Wizard PCR purification kit (Promega) or were subjected to agarose gel electrophoresis and the cDNA band extracted from the gel using a QIAEX II gel extraction kit (QIAGEN). Purified cDNA was sequenced using Big Dye chemistry (ABI, Perkin-Elmer), and the products obtained were analyzed using an ABI 377 automated sequencer (Perkin-Elmer). Molecular clones corresponding to the intergenic regions of L RNA and of part of the GP-C-coding region in the S segment (residues 247 to 759) were obtained by ligation of cDNA into the pGEM-T Easy vector (Promega) and transformation into Escherichia coli DH5
cells. DNA from positive colonies was amplified with a TempliPhi amplification kit (Amersham; based on rolling-circle amplification by the high-fidelity
29 DNA polymerase) by following the manufacturer's protocol. Mutant-spectrum complexity was calculated by scoring mutations relative to the consensus sequence defined by all the clones analyzed. Mutation frequencies for mutant spectra of the GP-C-coding region were determined as the number of different mutations (repeated mutations were counted only once) divided by the total number of nucleotides sequenced (the length of the GP-C-coding region multiplied by the number of clones analyzed).
The entire LCMV genome was sequenced using the primers and the strategy given in the supplementary material (http://www.cbm.uam.es/mkfactory.esdomain/webs/CBMSO/plt_LineasInvestigacion.aspx?IdObjeto=17). For sampling of partial genomic sequences, three LCMV genomic regions were amplified by RT-PCR and sequenced, as follows: in the S genomic RNA, the glycoprotein-coding region between residues 268 and 738 was amplified using primers GP247F (forward) and GP759R (reverse); in the nucleoprotein-coding region, residues 2244 to 2723 were amplified with primers NP2223F (forward) and NP2743R (reverse); in the L RNA, the polymerase-coding region between nucleotides 2950 and 4240 was amplified with primers L2929F (forward) and L4260R (reverse) (this set amplified the region corresponding to amino acids 978 to 1422, which encompasses premotif A to motif E, including the catalytic domain of the viral polymerase [24, 52, 87]). Consensus nucleotide sequences were determined using entire viral RNA populations as templates, except for the L intergenic region, for which consensus sequences were derived from the alignment of sequences from multiple individual molecular clones. Most nucleotide sequences were determined two or more times using different sense and antisense primers (described at http://www.cbm.uam.es/mkfactory.esdomain/webs/CBMSO/plt_LineasInvestigacion.aspx?IdObjeto=17). Nucleotide positions are given in the viral (genomic) sense and refer to the consensus genomic sequence determined in the present study. LCMV genomic nucleotide sequences were retrieved from GenBank (74-76) (GenBank accession numbers BANKIT 691408 and AY 894816). Sequences were compared using SEQMAN from the DNASTAR software package (Lasergene), and alignments were done with the GeneDoc program (K. B. Nicholas and H. B. Nicholas, Jr., Pittsburgh Supercomputing Center, Pittsburgh, Pa.). Computer-based predictions of RNA secondary structure were carried out using the MFOLD program (92) with default folding.
Detection of hypermutated molecules.
Conditions for reverse transcription and PCR were set to increase the probability of amplification of putative hypermutated molecules. The RT reaction was done with avian myeloblastosis virus (AMV) reverse transcriptase (Promega) at 37°C for 45 min and then at 94°C for 5 min. This incubation temperature is that recommended for the synthesis of first-strand cDNA using random oligonucleotide primers. One microliter of RNA was thawed at 70°C for 3 min, and after cooling on ice, 8 pmol of reverse primer was added and the mixture incubated at 65°C for 5 min. A mixture of 2.5 U of AMV and 7.5 U of RNAsin RNA inhibitor (Promega) was added in a final reaction volume of 20 µl. Five microliters of the first-strand reaction was used in a touchdown PCR in a volume of 50 µl with 1.25 U of Pfu DNA polymerase and 16 pmol of reverse and forward primers. The reaction conditions were as follows: denaturation hold at 95°C for 2 min, followed by 20 cycles at 95°C for 1 min, 60°C for 30 s (starting at an annealing temperature of 60°C, with a decrease of 0.5°C per cycle), and 72°C for 2 min, and finally 10 additional cycles at an annealing temperature of 47°C. Assuming that the more-frequent mutations induced by FU are transitions of the types A
G and U
C (41, 65, 67, 72, 79), oligonucleotides with different degrees of degeneracy (which should favor hybridization with mutated templates) were synthesized. Oligonucleotides L3662F (forward) and L4260R (reverse) were modified to produce a battery of degenerate primers to amplify residues 3662 to 4260 of the L polymerase-coding region (Table 1). The cDNAs synthesized were purified with a Wizard PCR purification kit (Promega), ligated into the pGEM-T Easy vector (Promega), and used to transform E. coli DH5
; positive bacterial clones were amplified with the TempliPhi amplification kit (Amersham). Purified cDNA and molecular clones were subjected to cycle sequencing. Standard RNA consisted of an SP6 RNA polymerase runoff transcript of a molecular clone that contained the polymerase region amplified with primers L3662F and L4260R in the genomic sense. Quantification of RNAs of control and FU-treated virus populations was done by two-step RT-PCR with ThermoScript RT (Invitrogen) at 60°C for 45 min, followed by a hold at 85°C for 5 min; 2 µl of the first-strand reaction product was used in a 20-µl reaction mixture with the Fast Start DNA Master SYBR Green I kit (Roche) in a Light Cycler instrument (Roche). A standard curve employed serial 10-fold dilutions of the standard RNA from 1010 to 105 molecules. Mutation frequencies were calculated as described above.
|
View this table: [in a new window] |
TABLE 1. Oligonucleotide primers with position degeneracies used in the present study
|
|
|
|---|
![]() View larger version (25K): [in a new window] |
FIG. 1. Effect of repeated mutagenesis-recovery cycles on LCMV infectivity during serial passages of LCMV ARM 53b in BHK-21 cells. LCMV (parental LCMV ARM 53b stock p0 with 1 x 108 PFU/ml) was serially passaged by infecting BHK-21 monolayers (2.8 x 106 cells) at an MOI of 0.01 PFU per cell in the absence (solid lines) or presence (dashed lines) of 100 µg/ml FU. Virus supernatants were harvested at 48 postinfection, clarified by centrifugation, and titrated on Vero cell monolayers. (A) Virus infectivity reduction following a single passage in the presence of FU. (B) Variation of infectivity of LCMV subjected to FU mutagenesis and recovery in the absence of FU. Passages in the absence or presence of FU are indicated by D and F, respectively (i.e., F1D3F1D2 means that the parental LCMV stock p0 was subjected successively to one passage in the presence of FU, three passages in the absence of FU, one passage in the presence of FU, and finally two passages in the absence of FU). Virus extinction at passages 2, 3, and 6 was confirmed by the absence of infectivity and the absence of virus-specific RT-PCR amplification after two further passages in the absence of drug. Procedures are detailed in Materials and Methods.
|
Fitness variations with an invariant consensus sequence. To study whether large fluctuations in replication capacity were associated with modification of the consensus LCMV genomic sequence, the entire genomic nucleotide sequences of the virus after one passage in the absence of drug (D1), the virus subjected to one passage in the presence of FU (F1), and the virus subjected to three subsequent recovery passages in the absence of FU (F1D3) were determined (the nomenclature of LCMV populations that defines the passage history in the presence or absence of FU is described in the legend for Fig. 1B). No differences in the consensus nucleotide sequence of the three viruses were observed. Although some differences were seen relative to the previously published LCMV ARM 53b sequence (73-76), our sequence is in good agreement with the LCMV genomic sequence deposited in GenBank with accession number AY894816. Under the standard amplification conditions, it was not possible to determine a defined nucleotide sequence (i.e., the peak pattern suggested the presence of insertions and/or deletions) for the intergenic region of segment L (nucleotides 470 to 550). RT-PCR at a higher temperature (25, 87) in 5% dimethyl sulfoxide was used to obtain cDNA that was cloned into the pGEM-T Easy vector (Promega). The nucleotide sequences of 13 clones from viral populations D1, F1, and F1D3 indicated a remarkable genetic heterogeneity with G-rich and C-rich sequences. Despite this heterogeneity, no evidence of dominance of different subpopulations was obtained for viruses D1, F1, and F1D3 (Fig. 2A). Two of the clones harbored large deletions which likely correspond to stem-loop structures deleted during the amplification or during the sequencing reaction by polymerase slippage (71, 89); however, heterogeneity in the viral RNA cannot be excluded. We have not subjected in vitro-transcribed RNA containing this region to RT-PCR to assess the contribution to the observed heterogeneity of errors during the amplification procedure. Previous studies have reported deletions in the intergenic regions of LCMV (75) and Lassa fever virus (87).
![]() ![]() View larger version (33K): [in a new window] |
FIG. 2. Nucleotide sequence and predicted secondary structure of an LCMV genomic region belonging to the L intergenic region. (A) Nucleotide sequence of residues 470 to 550, corresponding to the intergenic L region of clones from populations D1, F1, and F1D3. Sequences of clones are aligned under our consensus LCMV ARM 53b; at the bottom, the sequence previously published (termed here ARM 53b*) (76) is given. Our sequence is in agreement with a recent sequence (deposited in GenBank with accession number AY894816), determined independently of ours. Deletions introduced for best alignment of heteropolymeric stretches are indicated by dashes. (B) Predicted secondary structure of intergenic L residues 365 to 466 of Lassa fever virus and predicted secondary structures of residues 470 to 550 of LCMV ARM 53b, clones D1-2, F1-4, and F1D3-4 analyzed in the present study. Procedures for nucleotide sequence determination and secondary-structure prediction are described in Materials and Methods.
|
To rule out the possibility that the invariance of the consensus sequence reflected a genetically homogeneous RNA population, the amplified DNAs of the GP-C-coding region between residues 247 and 759 of the S segment of populations D1, F1, and FD3 were cloned into pGEM-T, and 8 to 13 molecular clones from each population were sequenced. The mutation frequencies for D1, F1, and FD3 were 1.4 x 104, 2.3 x 104, and 3.9 x 103, respectively (basal mutation frequency levels, determined with a runoff transcript, are included in Fig. 3), indicating that despite the invariance of the consensus sequence, the quasispecies contain highly heterogeneous mutant spectra (see also Discussion).
![]() View larger version (21K): [in a new window] |
FIG. 3. Mutation frequencies in control and FU-treated LCMV populations. RNAs from control (D) and FU-treated (F) virus populations were copied to cDNA (residues 3662 to 4260 of the L polymerase-coding region) under standard amplification conditions (RT reaction with ThermoScript RT at 60°C followed by Pfu DNA polymerase amplification at an annealing temperature of 59°C) or under less stringent amplification conditions (RT reaction with AMV RT at 37°C followed by Pfu amplification in a touchdown PCR) as indicated. Combinations of degenerate and wild-type primers are indicated (primers used are listed in Table 1, although the first letter, L, and the nucleotide position are omitted here for simplification). As an internal control, two dilutions of an LCMV L RNA runoff transcript (St. RNA, standard RNA) with similar numbers of molecules as the D and F virus populations were processed in parallel under both amplification conditions. Criteria for the design of degenerate primers and RT-PCR amplification conditions are detailed in Materials and Methods.
|
For primers without degeneracy, lowering the temperature during first-strand synthesis, followed by touchdown amplification, allowed detection of genomes with a number of mutations that remained undetected under standard, more-stringent conditions of amplification. Under these new amplification conditions, the mutation frequency for the F population was 4.7-fold higher than the value obtained under standard amplification conditions whereas for the D population the mutation frequency was 2.0-fold lower (Fig. 3). Also, for the F population, the mutation frequencies scored were lower when any combination of degenerate oligonucleotides was used than with standard oligonucleotides representing the wild-type LCMV consensus sequence. Control amplifications employing similar numbers of molecules of L RNA runoff transcripts as those in the D and F populations indicate that the differences in the mutation frequency values were not due to the amplification conditions (Fig. 3). Thus, the results do not support our initial hypothesis that hypermutated molecules may be dominant during the transition into error catastrophe (41). Rather, the results suggest that mutant spectra, composed of genomes harboring modest numbers of fitness-decreasing mutations, dominate the molecular landscape in the LCMV transition into error catastrophe, supporting a lethal defection model for LCMV extinction (40).
|
|
|---|
LCMV shows higher sensitivity to FU mutagenesis than foot-and-mouth disease virus (FMDV) in cell culture when the two viruses are subjected to similar FU doses (41, 65, 67, 72, 79). The basis for this difference is not known, but it may relate either to the size of the essential polymerase gene, which occupies 62.6% and 17.4% of the total genetic information of LCMV and FMDV, respectively; to a larger number of essential genomic sites susceptible to FU-induced mutations in LCMV than in FMDV; or to other replicative differences between the two viruses still to be elucidated (different affinities of the two polymerases for FU-triphosphate, etc.). The comparison of complete genomic LCMV sequences reported here has indicated that neither an abrupt decrease in infectivity associated with FU treatment nor the ensuing recovery of the capacity to produce infectious progeny upon replication in the absence of mutagen modified the consensus nucleotide sequence of the virus. Likewise, no variations were detected among LCMV populations subjected to one or several rounds of replication in the absence or presence of FU, or among populations that underwent alternations between the two passage regimens. In these cases the consensus sequences analyzed involved three genomic regions, one of which was a part of the polymerase (L) that included premotif A and motifs A to E (11, 24, 52, 70, 86, 87). The possibility that FU, in addition to its mutagenic activity documented by increases in mutant-spectrum complexity, could exert some mutagenesis-independent, inhibitory activity (i.e., directly on the viral polymerase or indirectly through alteration of nucleotide pools, etc.) on LCMV replication cannot be excluded. This would mean that FU could contribute to decreases in viral load, thereby favoring mutagenesis-induced extinction, as documented with FMDV (65-67).
It could be argued that fitness variations despite an invariant consensus sequence may occur without modification of the mutant spectrum, simply by an increase of dominance of a preexisting subset of genomes with the same consensus sequence. A static mutant spectrum is highly unlikely given the number of LCMV progeny per cell (Fig. 1), since it would imply mutation rates much lower than those that have been calculated for RNA viruses (29). Low mutation rates are not supported by the estimates of the frequencies of cytotoxic T-lymphocyte (CTL)- and antibody-escape mutants in LCMV populations (1, 15, 43, 49; reviewed in reference 78). Fitness variations without modification of consensus genomic sequences have been reported for vesicular stomatitis virus clones in the process of fitness recovery after bottleneck passages (59).
The L RNA intergenic region may function as a transcription termination signal, as previously suggested for LCMV (69, 74), other arenaviruses (25, 34, 44), and the stem-loop present in the S RNA intergenic region of LCMV (55), Pichinde virus (7), and Lassa fever virus (8, 16). RNA secondary-structure predictions suggest that the L RNA intergenic region of LCMV ARM 53b forms a stem-loop stabilized by 15 GC pairs (Fig. 3B), similar to a loop predicted for the L RNA intergenic regions of Tacaribe and Lassa fever virus (25, 34, 44). We have not distinguished whether the nucleotide sequence heterogeneity seen in the L intergenic region (Fig. 3A) corresponds to a heterogeneity present in the LCMV used (perhaps with functional implications) or is due to the production of altered copies during RT-PCR, or both.
Although it cannot be excluded that highly mutated or hypermutated RNA molecules derived from LCMV RNA may be present during the transition to error catastrophe, attempts using standard and multiply substituted oligonucleotide primers have failed to reveal the dominance of such genomes (Table 1 and Fig. 3). In the amplification with primers L3662F and L4260 RY4, one sequence was found that had nine mutations, of which six were U
C transitions and the others were A
G, C
A, and G
A. The mutation frequency for this molecule relative to the consensus sequence (1.6 x 102) would lie within the lower limit in the range of 3.3 x 102 to 1.0 x 101 mutations per nucleotide for the U
C (or A
G) hypermutation, previously characterized in genomes of measles virus, human parainfluenza virus 3, vesicular stomatitis virus, and Rous-associated virus, which has been related to posttranscriptional unwinding and mutagenic activities (13, 36, 54, 56, 57, 61). Hypermutated genomes have not been identified, either, among multiple molecular clones copied from preextinction populations of FMDV (4, 67, 79). In the FU-mediated transition to extinction of LCMV in the course of a persistent infection of BHK-21 cells, decreases in viral infectivity preceded decreases in viral RNA levels (40). The specific infectivity of LCMV reached values 102- to 103-fold lower than those for LCMV from parallel cultures not subjected to FU mutagenesis (40). Similar decreases in specific infectivity were noted for mutagenized FMDV during persistent or cytolytic infections (4, 37, 65, 66). These quantifications of specific infectivity, together with the complexity of mutant spectra of mutagenized LCMV and FMDV populations (4, 37, 40, 41, 65-67, 72, 79), suggest that rather than requiring extensive mutagenesis, virus extinction may be mediated by the presence of defective genomes (which have been termed "defectors" in several studies of RNA virus evolution [85, 90]). Therefore, the present evidence suggests that low-activity mutagens capable of being incorporated specifically by viral polymerases may be sufficient to enrich the defector population to a level at which it may interfere with standard genome replication. A likely molecular mechanism mediating this type of interference is the effect of multiple, low-frequency dominant-negative mutants generated in mutagenized virus populations, an extension to an intrapopulation level of the mechanisms of interference between different viruses, well documented in classical viral genetics (reviewed in references 2, 17, and 91). Because such interfering genomesdue to their partial or complete defectivenesscannot dominate the population, no alteration of consensus sequences is observed despite profound fitness losses. Thus, the results reported here support a lethal defection model for the mutagenesis-induced genetic meltdown that results in viral extinction (40). Treatment with the mutagenic nucleoside analogue ribavirin of cell cultures infected with poliovirus or Hantaan virus resulted in decreases in infectivity that preceded decreases in viral RNA levels (19-21, 77), again suggesting a role of mutagenesis-induced defective genomes that maintain their replication competence prior to detectable losses of infectivity.
The results presented here underscore the relevance of mutant spectra in determining virus behavior, illustrated here by the capacity of LCMV to either decrease or increase its fitness without any alteration of the consensus nucleotide sequence. This observation must be added to results with several other virus-host systems showing suppressive or modulating effects of mutant spectra on specific types of variants with superior fitness (9, 14, 22, 37, 84).
This work was supported by grants BMC 2001-1823-C02-01 and CAM 08.2/0015/2001.1 and by Fundación R. Areces. A.G.-P. was supported by a postdoctoral contract from CAM and a RyC contract from MCyT.
Present address: Departamento de Genética, Universidad de Málaga, 29071 Málaga, Spain. ![]()
|
|
|---|
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»