Previous Article | Next Article ![]()
Journal of Virology, May 2005, p. 5421-5427, Vol. 79, No. 9
0022-538X/05/$08.00+0 doi:10.1128/JVI.79.9.5421-5427.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
,
Antonio V. Bordería,2,
Esteban Domingo,1 and
Cecilio López-Galíndez2*
Centro de Biología Molecular "Severo Ochoa," CSIC-UAM, Universidad Autónoma, Cantoblanco,1 Servicio de Virologia Molecular, Centro Nacional de Microbiología, Instituto de Salud Carlos III, Majadahonda, Madrid, Spain2
Received 8 October 2004/ Accepted 15 December 2004
|
|
|---|
|
|
|---|
The measurement of fitness values under different experimental conditions in several RNA viruses has led to the observation that fitness variations are extremely dependent on the population size of the evolving RNA virus population. In particular, repeated bottleneck passages (the multiplication of viral genomes from only one or a small number of genomes) result in average fitness losses in several different RNA viruses (7, 11, 12). The decrease in fitness mediated by repeated bottlenecks has been interpreted as the result of an accentuation of Muller's ratchet effect (30) due to the incorporation of deleterious mutations and the lack of compensatory mechanisms (26) such as recombination (10-12, 27, 35). The operation of Muller's ratchet effect in viruses has been demonstrated with bacteriophage
6 (7), vesicular stomatitis virus (11), foot-and-mouth disease virus (13), and HIV-1 (47, 48). In contrast, massive infections of cells with large numbers of infectious viruses lead to rapid fitness gains (10, 14, 32). As the population size of an RNA virus increases, the repertoire of multiple mutants that may be subjected to positive selection also increases. This process leads to the accumulation of more fit genomes and a loss or decrease in the frequency of less adapted genomes. For retroviruses, decreases in fitness as a result of serial bottleneck passages were first documented for HIV-1 after plaque-to-plaque passages on MT-4 cells (48). For this virus, fitness losses were unexpectedly drastic and rapid compared with those experienced by other RNA viruses, such as the bacteriophage
6 (7), vesicular stomatitis virus (11), and foot-and-mouth disease virus (13), which were subjected to similar passage regimens. Only 4 of 10 HIV-1 clones produced viable progeny after 15 plaque-to-plaque transfers, and three of the four survivors displayed important decreases in fitness which were associated with an unusual distribution of mutations (47). The mutation frequencies in gag were threefold higher than those in env, and the complete absence of replacements in the V3 loop of gp120 was remarkable. The amino acid replacements observed were also atypical, and several of the amino acid changes have not been previously recorded among natural HIV-1 isolates (47).
Here we report an increase in fitness and describe the associated mutations that occur when debilitated HIV-1 clones are subjected to large population passages. The results document a surprising preferential location of mutations in the 5' untranslated leader sequence. The analysis also revealed that few mutations (an average of five mutations per clone) were associated with fitness recovery and that most of the mutations in the coding regions were amino acid replacements. The large numbers of reversions (25%) and amino acid substitutions to regain the most common amino acids found for the corresponding positions among natural HIV-1 isolates were particularly striking. These observations reflect a high selective pressure during fitness recovery.
|
|
|---|
Population passages of virus. Adaptive passages of large virus populations were performed as follows (Fig. 1). MT-4 cells (5 x 106) were infected with 5 x 105 PFU of each clone (multiplicity of infection of 0.1 PFU/cell) in a final volume of 5 ml. The virus was recovered in the supernatant when cytopathology was complete (at 5 to 7 days postinfection). For each passage, fresh MT-4 cells were infected with a 1:10 dilution from the supernatant of the previous passage, taking into account that viral titers do not change significantly during serial passages of HIV-1 (37). The cells were cultured in closed bottles, and all of the procedures were aimed at minimizing the possibility of cross-contamination, which was monitored by the use of uninfected control cultures which were maintained in parallel. No evidence of contamination was observed at any time.
![]() View larger version (15K): [in a new window] |
FIG. 1. Scheme of passages of HIV-1 clones subjected to plaque-to-plaque transfers (filled squares) and large population passages (open circles) in MT-4 cells. The experimental procedures and the origins of the natural isolate s61 and all of the different clones were previously reported (41, 48) and are given in Materials and Methods. HIV-1 clones are indicated by a letter followed by a number that gives the total number of plaque-to-plaque transfers and p1 to p10, indicating the number of large population transfers.
|
RNA extraction, PCR amplification, and nucleotide sequencing. RNAs were extracted with a guanidinium isothiocyanate lysis buffer and glass milk as described previously (5). HIV-1 RNAs were amplified by reverse transcription-PCR and nested PCR. The first amplification was performed by use of the Access RT-PCR system (Promega) and six different pairs of primers that covered the entire genome of HIV-1. Internal amplifications were performed with EcoTaq DNA polymerase (Ecogen) and internal primers appropriate for each amplification. The primer sequences used for all amplifications are available upon request. Both external and internal amplifications involved 35 cycles, with temperatures chosen according to the compositions of the oligonucleotide primers. Negative amplification controls were run in parallel to monitor the absence of contamination.
Complete genomic nucleotide sequences were determined for the initial and final viral populations of the plaque-to-plaque passages and were reported in a previous work (47). The complete genomic nucleotide sequences of the final viral populations after the large population passages were determined for this study. Sequences were determined for the two cDNA strands by use of an ABI 373 automatic sequencer. Multiple sequence alignments were performed with the CLUSTAL W program (45). The numbering of nucleotide changes detected is reported relative to the numbering of isolate HXB2 included in the Los Alamos database alignment (21).
HTA.
HTAs were performed as previously reported (48) by use of a cDNA copy of a 650-bp fragment corresponding to the V1-V2 region of the env gene. For amplification of this DNA, a nested PCR was performed with the primers 91ECU (5' CTTAGGCATCTCCTATGGC 3'; positions 5956 to 5974 [numbering relative to HXB2]) and 92ECD (5' GGAGCAGCAGGAAGCAC 3'; complementary to positions 7810 to 7827) and the nested primers 99ECU (5' AGAGCAGAAGACAGTGGC 3'; positions 6205 to 6223) and 52EV1D (5' TAATGTATGGGAATTGGCTCAA 3'; complementary to positions 6874 to 6896). Amplifications were carried out for 35 cycles of 94°C for 1 min, 55°C for 55 s, and 72°C for 1 min, with a final extension cycle of 10 min at 72°C. In order to determine the proportion of each molecular species during competition passages, we labeled the amplified cDNA of clone J1 with [
-32P]dCTP (3,000 Ci/mmol) by PCR amplification. About 10,000 cpm of this radioactive PCR probe was mixed with 50 to 100 ng of the unlabeled second-round PCR product from the competing virus in annealing buffer (0.1 M NaCl, 10 mM Tris-HCl [pH 7.8], and 2 mM EDTA). The DNA mixtures were denatured at 94°C for 2 min and then quickly cooled. Heteroduplexes were resolved in denaturing 8% polyacrylamide-15% urea gels in TBE (0.1 mM Tris-borate, 89 mM boric acid, and 2 mM EDTA, pH 8). Autoradiograms were obtained by exposing gels on a Fuji 2000 instrument for 2 h. Quantification of the ratio of the J1 cDNA (homoduplex) to the cDNA of the competing virus (heteroduplex) was performed by densitometry with the help of the PCBAS program.
|
|
|---|
![]() View larger version (18K): [in a new window] |
FIG. 2. Fitness evolution of each clone during plaque-to-plaque and large population passages. Fitness values were calculated from the exponential slopes of the corresponding fitness vectors (48) as described previously (19). Plaque-to-plaque transfers are indicated by filled squares, and large population passages are indicated by open circles. The fitness of clone D15p10 could not be calculated, as explained in Table 1. The asterisk indicates a debilitated clone (G7) that could not produce plaques after plaque-to-plaque transfer 7. The fitness of this clone could not be determined by competition passages, as described in Materials and Methods.
|
|
View this table: [in a new window] |
TABLE 1. Fitness values of clones after 15 plaque-to-plaque transfers and 1, 5, and 10 large population passages
|
G, C
T, and G
A changes accounted for 35, 20, and 15%, respectively, of the total mutations (Table 2). Mutation frequencies ranged from 2.1 x 104 to 7.5 x 104 substitutions per nucleotide (mean, 5.4 x 104 substitutions per nucleotide). Thirteen substitutions were in untranslated regulatory regions, and eight were in coding regions. Among the latter mutations, six were nonsynonymous (75% of the total) and two were synonymous. Most of these amino acid changes replaced residues which are present with low frequencies in the database with amino acids that are present in >90% of HIV-1 isolates (21), and the degree of acceptability is high in most of the amino acid replacements (15) (Table 3).
![]() View larger version (11K): [in a new window] |
FIG. 4. Distribution in the complete genome of all mutations appearing during large population passages. Horizontal bars below the gene organization scheme represent the genomes of the four clones studied. Vertical lines illustrate nonsynonymous mutations, whereas vertical arrows indicate synonymous mutations. The 5' U5 region was expanded in order to more precisely map the mutations accumulated in the clones.
|
|
View this table: [in a new window] |
TABLE 2. Numbers and types of mutations in the genomes of HIV-1 clones subjected to large population transfers
|
|
View this table: [in a new window] |
TABLE 3. Amino acid replacements and synonymous mutations during large population passages of HIV-1
|
![]() View larger version (31K): [in a new window] |
FIG. 3. Localization of mutations accumulated in the U5 PBS loop during large population passages. (A) Sequence alignment of viruses before and after fitness recovery in the U5 PBS region in comparison with the HXB2 and consensus B sequences, which are marked in bold. Boxes indicated nucleotide changes associated with large population passages. Asterisks indicate reversions to the sequence present before the fitness loss (47, 48). The PAS and PBS sequences are shaded. (B) Secondary structure model of the U5 PBS leader region of HIV-1 (LAI sequence). The locations of the PBS and the PAS are indicated. The positions of the mutations accumulated during fitness recovery in all final p10 viral populations are marked. Watson-Crick base pairing is marked with dashes. Clones that are underlined and marked with asterisks indicate reversions to the sequence that was present before the fitness loss (47, 48). Modified from reference 3.
|
Accumulation of mutations in the primer binding site hairpin. An analysis of the mutations that appeared during fitness recovery revealed that they were not randomly distributed along the genome, as 12 substitutions (60% of the total) were localized in the 5' untranslated region, one was in the 3' U3 region, and the remaining ones were situated within the coding regions. The 12 mutations in the 5' untranslated leader region mapped within the primer binding site (PBS) loop (Fig. 3 and 4). The majority of these mutations correspond to conserved nucleotides in the sequence database (20). The 335-nucleotide 5' untranslated leader region is structurally conserved in HIV-1 (21) and contains multiple sequence and structural motifs that play important roles in viral replication (4). Several mutations were found in more than one population. The G570A mutation was found in populations D15p10 and K15p10, and the A584G and T584G mutations were found in populations D15p10 and G7p10, although the A587G mutation which occurred in all four clones was more intriguing. The observation of the same mutation in independent populations is an example that illustrates evolutionary convergence and may support the deterministic behavior of RNA viruses when selection predominates (38, 39). Two of the mutations found in this region (G570A in K15p10 and T584G in G7p10) were reversions to the sequence present in the virus before fitness loss. The T584G reversion was also found to be the only mutation incorporated into clone G7p1 during fitness recovery. The mutation at position 584 is located in an eight-nucleotide motif in the U5 region (positions 576 to 584 [numbering according to HXB2]) (21) (Fig. 3) which is critical for the tRNA3Lys-mediated initiation of reverse transcription in vitro (3). This motif is referred to as the primer activation signal (PAS) and is important for the activation of the PBS-bound tRNA primer (2).
|
|
|---|
The large majority of mutations associated with fitness losses through repeated bottlenecks occurred in the gag gene (47), with mutation frequencies in gag being threefold higher than those in env. Surprisingly, most of the mutations associated with large population passages and fitness recovery were not located in any of the open reading frames, but were in the 5' untranslated region. An analysis of the coding region mutations occurring during large population passages revealed that only one of the amino acid replacements had a low frequency in the HIV-1 sequence database (21), and four of the six amino acid replacements had a frequency above 90% (21) (Table 3).
The observation in our analysis of a preferential location of mutations in the 5' leader region, with 12 of the 20 mutations associated with fitness recovery, was particularly interesting. The 5' untranslated leader region of the viral genome is a highly structured RNA domain which is conserved among isolates of the same subtype (21). This region includes the PBS hairpin, displaying the PBS motif, which does not constitute a hot spot for genetic variation (21). Within this region, several interactions are required for efficient binding of the tRNA to the PBS in order to initiate reverse transcription, which is a critical step for viral replication (2, 3, 17). First, a cellular tRNA3Lys, which is used as a primer, binds the PBS motif in the 5' end of the viral RNA (4). Second, a new critical motif, the PAS, was recently identified for tRNA3Lys as a mediator of the initiation of reverse transcription (3). This motif interacts with the anti-PAS motif of the tRNA and triggers the initiation of retrotranscription. Mutations in the PAS motif have been associated with a reduced efficiency of tRNA-primed reverse transcription (1, 2). Interestingly, two of the low-fitness viruses in our study (D15 and G7) had mutations in the PAS motif. In clone G7, the mutation reverted after a fitness gain (T
G), and in D15 it changed to the common nucleotide (A
G), rendering in both cases the PAS wild-type sequence complementary to the anti-PAS sequence of the tRNA. RNA structure analysis performed with MFOLD (http://bioweb.pasteur.fr/seqanal/interfaces/mfold-simple.html) indicated a decrease in
G between the RNA structure for large populations and that for plaque-to-plaque passaged virus, showing that the T/A584G (in clones D15 and G7) and A587G (present in all recovered clones) mutations conferred more stability to the 5' U5 structure (data not shown). Moreover, the accumulation of mutations in different positions of the PBS loop and around the PAS motif (Fig. 3) supports the proposed hypothesis (2) that additional interactions within this region are also required to assemble an initiation-competent and processive reverse transcription complex. The preferential location of mutations associated with fitness recovery in the PBS loop suggests that a more efficient initiation of retrotranscription is a critical factor for the optimization of in vitro HIV-1 viral replication.
Both fitness loss upon plaque-to-plaque transfers and fitness recovery upon large population passages were more rapid and larger for HIV-1 than for other RNA viruses analyzed (13, 14, 47, 48). During plaque-to-plaque transfers (44), 6 of 10 HIV-1 clones became extinct, in contrast to the results of other Muller's ratchet studies with viruses for which extinction was absent or rare (6, 10). For large population passages, the HIV-1 viral populations showed very significant changes in fitness values, ranging from 1.28- to 1,780-fold, which were much higher than those with other systems (14, 33) and those obtained with fewer serial passages.
In the first five large population passages, clones G and I showed slight decreases in fitness (Table 1). Although this decrease was not significant compared to the global fitness gains observed, this observation may be related to the increase in the viral population size and also the viral diversity occurring during the transition from homogeneous clonal populations to large population passages.
It is remarkable that five mutations associated with the recovery of fitness were reversions, with three located in the coding region (Table 3) and two located in the PBS loop. Reversions in the PBS loop were recently associated with an increase in replication capacity when low-replication viruses that employed the nonself tRNAPro or tRNA1,2Lys recovered the ability to use the tRNA3Lys molecule (1). The occurrence of reversions in the PAS motif indicates a high degree of functional constraint, consistent with the high structural conservation observed in this region (2, 3). Also, the abundance of reversions suggests that in HIV-1, after drastic losses, reversion is the preferred mechanism for the in vitro recovery of viral fitness, at least in this genomic region. In contrast to what occurs in foot-and-mouth disease virus, in which reversions are virtually absent during fitness recovery (<10%) (14, 22), compensatory mutations were favored in HIV-1 and also in vesicular stomatitis virus (33).
The natural history of HIV-1 infection is characterized by fluctuations in the viral population size, with transitions from extremely large populations in many infected patients to periods with drastic reductions in population sizes such as occurs during transmission. The changes in viral populations when patients undergo antiretroviral therapy, from undetectable levels for successful potent antiviral therapy to exponential growth during treatment failures, or during structured treatment interruptions are also significant. This study gives a better understanding of in vitro HIV-1 population dynamics and of the dramatic consequences of bottleneck and large population passages on HIV-1 fitness, which may be relevant to in vivo infections.
Work at Centro de Biología Molecular "Severo Ochoa" was supported by grant BMC 2001-1823-C02-01 and by Fundación Ramón Areces. Work in CNM was supported by grants MPY 1359/02 and MPY 1028/04 and in part by the "Red Tematica Cooperativa de Investigación en SIDA (Red de grupos 173) del FISss." A.V.B. was supported by an Instituto Carlos III fellowship. E.Y. was supported by a postdoctoral fellowship from Comunidad Autonoma de Madrid.
E.Y. and A.V.B. contributed equally to this work. ![]()
Present address: New England Primate Research Center, Harvard Medical School, Southborough, MA 01772-9102. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»