Mutagenicity and pausing of HIV reverse transcriptase during HIV plus-strand DNA synthesis

Removal of ribonucleotides by p53 protein incorporated during DNA synthesis by HIV-1 reverse transcriptase

OBJECTIVE(S)

HIV-1 reverse transcriptase frequently incorporates ribonucleotides into the proviral DNA in macrophages, but not in lymphocytes. The enzyme exerts an efficient ribonucleotide-terminated primer extension capacity. Furthermore, ribonucleotide-editing repair is attenuated in macrophages. Tumor suppressor p53 protein, displaying an intrinsic 3'→5' exonuclease activity, was found to be involved in efficient proofreading of base-base mismatches produced during DNA synthesis. As the presence of proofreading activity is cardinal for the DNA synthesis accuracy, it was of interest to assess whether p53 can serve as a trans-acting proofreader for HIV-1 reverse transcriptase during ribonucleotide incorporation.

DESIGN

We investigated the potential involvement of cytoplasmic p53 in error correction during insertion of ribonucleotides into DNA by recombinant HIV-1 reverse transcriptase in a p53-proficient and deficient background.

METHODS

Primer extension reactions were carried out to elucidate the incorporation and removal of ribonucleotides.

RESULTS

The biochemical studies suggest that p53 is involved in a ribonucleotide damage-associated repair mechanism through its capacity to remove preformed 3'-terminal ribonucleotides, to decrease ribonucleotide incorporation and to prevent the 3'-ribo-terminated primer extension during ongoing DNA synthesis by HIV-1 reverse transcriptase. A positive correlation exists between the presence of endogenous p53 and decrease in stable incorporation of ribonucleotides into DNA with p53-harboring lysates of HCT116 cells. p53, by preferential removal of purine over pyrimidine ribonucleotides, may affect the ribonucleotide mutation spectra produced by HIV-1 reverse transcriptase.

CONCLUSION

The data implies that p53 can excise incorrect sugar in addition to base mispairs, thereby expanding the role of p53 in the repair of nucleic acids replication errors.

Pre-steady state kinetic analysis of HIV-1 reverse transcriptase for non-canonical ribonucleoside triphosphate incorporation and DNA synthesis from ribonucleoside-containing DNA template

Non-dividing macrophages maintain extremely low cellular deoxyribonucleotide triphosphate (dNTP) levels, but high ribonucleotide triphosphate (rNTP) concentrations. The disparate nucleotide pools kinetically forces Human Immunodeficiency Virus 1 (HIV-1) reverse transcriptase (RT) to incorporate non-canonical rNTPs during reverse transcription. HIV-1 RT pauses near ribonucleoside monophosphates (rNMPs) embedded in the template DNA, which has previously been shown to enhance mismatch extension. Here, pre-steady state kinetic analysis shows rNTP binding affinity (Kd) of HIV-1 RT for non-canonical rNTPs was 1.4- to 43-fold lower, and the rNTP rate of incorporation (kpol) was 15- to 1551-fold slower than for dNTPs. This suggests that RT is more selective for incorporation of dNTPs rather than rNTPs. HIV-1 RT selectivity for dNTP versus rNTP is the lowest for ATP, implying that HIV-1 RT preferentially incorporates ATP when dATP concentration is limited. We observed that incorporation of a dNTP occurring one nucleotide before an embedded rNMP in the template had a 29-fold greater Kd and a 20-fold slower kpol as compared to the same template containing dNMP. This reduced the overall dNTP incorporation efficiency of HIV-1 RT by 581-fold. Finally, the RT mutant Y115F displayed lower discrimination against rNTPs due to its increase in binding affinity for non-canonical rNTPs. Overall, these kinetic results demonstrate that HIV-1 RT utilizes both substrate binding and a conformational change during: (1) enzymatic discrimination of non-canonical rNTPs from dNTPs and (2) during dNTP primer extension with DNA templates containing embedded rNMP.

Human immunodeficiency virus reverse transcriptase displays dramatically higher fidelity under physiological magnesium conditions in vitro

The fidelity of human immunodeficiency virus (HIV) reverse transcriptase (RT) has been a subject of intensive investigation. The mutation frequencies for the purified enzyme in vitro vary widely but are typically in the 10(-4) range (per nucleotide addition), making the enzyme severalfold less accurate than most polymerases, including other RTs. This has often been cited as a factor in HIV's accelerated generation of genetic diversity. However, cellular experiments suggest that HIV does not have significantly lower fidelity than other retroviruses and shows a mutation frequency in the 10(-5) range. In this report, we reconcile, at least in part, these discrepancies by showing that HIV RT fidelity in vitro is in the same range as cellular results from experiments conducted with physiological (for lymphocytes) concentrations of free Mg(2+) (~0.25 mM) and is comparable to Moloney murine leukemia virus (MuLV) RT fidelity. The physiological conditions produced mutation rates that were 5 to 10 times lower than those obtained under typically employed in vitro conditions optimized for RT activity (5 to 10 mM Mg(2+)). These results were consistent in both commonly used lacZα complementation and steady-state fidelity assays. Interestingly, although HIV RT showed severalfold-lower fidelity under high-Mg(2+) (6 mM) conditions, MuLV RT fidelity was insensitive to Mg(2+). Overall, the results indicate that the fidelity of HIV replication in cells is compatible with findings of experiments carried out in vitro with purified HIV RT, providing more physiological conditions are used.

IMPORTANCE

Human immunodeficiency virus rapidly evolves through the generation and subsequent selection of mutants that can circumvent the immune response and escape drug therapy. This process is fueled, in part, by the presumably highly error-prone HIV polymerase reverse transcriptase (RT). Paradoxically, results of studies examining HIV replication in cells indicate an error frequency that is ~10 times lower than the rate for RT in the test tube, which invokes the possibility of factors that make RT more accurate in cells. This study brings the cellular and test tube results in closer agreement by showing that HIV RT is not more error prone than other RTs and, when assayed under physiological magnesium conditions, has a much lower error rate than in typical assays conducted using conditions optimized for enzyme activity.

Effect of ribonucleotides embedded in a DNA template on HIV-1 reverse transcription kinetics and fidelity

HIV-1 reverse transcriptase (RT) frequently incorporates ribonucleoside triphosphates (rNTPs) during proviral DNA synthesis, particularly under the limited dNTP conditions found in macrophages. We investigated the mechanistic impacts of an rNMP embedded in DNA templates on HIV-1 RT-mediated DNA synthesis. We observed that the template-embedded rNMP induced pausing of RT and delayed DNA synthesis kinetics at low macrophage dNTP concentrations but not at high T cell dNTP concentrations. Although the binding affinity of RT to the rNMP-containing template-primer was not altered, the dNTP incorporation kinetics of RT were significantly reduced at one nucleotide upstream and downstream of the rNMP site, leading to pause sites. Finally, HIV-1 RT becomes more error-prone at rNMP sites with an elevated mismatch extension capability but not enhanced misinsertion capability. Together these data suggest that rNMPs embedded in DNA templates may influence reverse transcription kinetics and impact viral mutagenesis in macrophages.

Intracellular nucleotide levels and the control of retroviral infections

Retroviruses consume cellular deoxynucleoside triphosphates (dNTPs) to convert their RNA genomes into proviral DNA through reverse transcription. While all retroviruses replicate in dividing cells, lentiviruses uniquely replicate in nondividing cells such as macrophages. Importantly, dNTP levels in nondividing cells are extremely low, compared to dividing cells. Indeed, a recently discovered anti-HIV/SIV restriction factor, SAMHD1, which is a dNTP triphosphohydrolase, is responsible for the limited dNTP pool of nondividing cells. Lentiviral reverse transcriptases (RT) uniquely stay functional even at the low dNTP concentrations in nondividing cells. Interestingly, Vpx of HIV-2/SIVsm proteosomally degrades SAMHD1, which elevates cellular dNTP pools and accelerates lentiviral replication in nondividing cells. These Vpx-encoding lentiviruses rapidly replicate in nondividing cells by encoding both highly functional RTs and Vpx. Here, we discuss a series of mechanistic and virological studies that have contributed to conceptually linking cellular dNTP levels and the adaptation of lentiviral replication in nondividing cells.

The A-nucleotide preference of HIV-1 in the context of its structured RNA genome

A bipartition of HIV-1 RNA genome sequences into single- and double-stranded nucleotides is possible based on the secondary structure model of a complete 9 kb genome. Subsequent analysis revealed that the well-known lentiviral property of A-accumulation is profoundly present in single-stranded domains, yet absent in double-stranded domains. Mutational rate analysis by means of an unrestricted model of nucleotide substitution suggests the presence of an evolutionary equilibrium to preserve this biased nucleotide distribution.

The Impact of Macrophage Nucleotide Pools on HIV-1 Reverse Transcription, Viral Replication, and the Development of Novel Antiviral Agents

Macrophages are ubiquitous and represent a significant viral reservoir for HIV-1. Macrophages are nondividing, terminally differentiated cells, which have a unique cellular microenvironment relative to actively dividing T lymphocytes, all of which can impact HIV-1 infection/replication, design of inhibitors targeting viral replication in these cells, emergence of mutations within the HIV-1 genome, and disease progression. Scarce dNTPs drive rNTP incorporation into the proviral DNA in macrophages but not lymphocytes. Furthermore, the efficacy of a ribose-based inhibitor that potently inhibits HIV-1 replication in macrophages, has prompted a reconsideration of the previously accepted dogma that 2′-deoxy-based inhibitors demonstrate effective inhibition of HIV-1 replication. Additionally, higher levels of dUTP and rNTP incorporation in macrophages, and lack of repair mechanisms relative to lymphocytes, provide a further mechanistic understanding required to develop targeted inhibition of viral replication in macrophages. Together, the concentrations of dNTPs and rNTPs within macrophages comprise a distinctive cellular environment that directly impacts HIV-1 replication in macrophages and provides unique insight into novel therapeutic mechanisms that could be exploited to eliminate virus from these cells.

The origin of genetic diversity in HIV-1

One of the hallmarks of HIV infection is the rapid development of a genetically complex population (quasispecies) from an initially limited number of infectious particles. Genetic diversity remains one of the major obstacles to eradication of HIV. The viral quasispecies can respond rapidly to selective pressures, such as that imposed by the immune system and antiretroviral therapy, and frustrates vaccine design efforts. Two unique features of retroviral replication are responsible for the unprecedented variation generated during infection. First, mutations are frequently introduced into the viral genome by the error prone viral reverse transcriptase and through the actions of host cellular factors, such as the APOBEC family of nucleic acid editing enzymes. Second, the HIV reverse transcriptase can utilize both copies of the co-packaged viral genome in a process termed retroviral recombination. When the co-packaged viral genomes are genetically different, retroviral recombination can lead to the shuffling of mutations between viral genomes in the quasispecies. This review outlines the stages of the retroviral life cycle where genetic variation is introduced, focusing on the principal mechanisms of mutation and recombination. Understanding the mechanistic origin of genetic diversity is essential to combating HIV.

Frequent incorporation of ribonucleotides during HIV-1 reverse transcription and their attenuated repair in macrophages

Macrophages are well known long-lived reservoirs of HIV-1. Unlike activated CD4(+) T cells, this nondividing HIV-1 target cell type contains a very low level of the deoxynucleoside triphosphates (dNTPs) required for proviral DNA synthesis whereas the ribonucleoside triphosphate (rNTP) levels remain in the millimolar range, resulting in an extremely low dNTP/rNTP ratio. Biochemical simulations demonstrate that HIV-1 reverse transcriptase (RT) efficiently incorporates ribonucleoside monophosphates (rNMPs) during DNA synthesis at this ratio, predicting frequent rNMP incorporation by the virus specifically in macrophages. Indeed, HIV-1 RT incorporates rNMPs at a remarkable rate of 1/146 nucleotides during macrophage infection. This greatly exceeds known rates for cellular replicative polymerases. In contrast, little or no rNMP incorporation is detected in CD4(+) T cells. Repair of these rNMP lesions is also substantially delayed in macrophages compared with CD4(+) T cells. Single rNMPs embedded in a DNA template are known to induce cellular DNA polymerase pausing, which mechanistically contributes to mutation synthesis. Indeed, we also observed that embedded rNMPs in a dsDNA template also induce HIV-1 RT DNA synthesis pausing. Moreover, unrepaired rNMPs incorporated into the provirus during HIV-1 reverse transcription would be generally mutagenic as was shown in Saccharomyces cerevisiae. Most importantly, the frequent incorporation of rNMPs makes them an ideal candidate for development of a new class of HIV RT inhibitors.

HIV-1 nucleocapsid protein increases strand transfer recombination by promoting dimeric G-quartet formation

A preferred site for HIV-1 recombination was identified in vivo and in vitro surrounding the beginning of the HIV-1 gag gene. This G-rich gag hotspot for recombination contains three evenly spaced G-runs that stalled reverse transcriptase. Disruption of the G-runs suppressed both the associated pausing and strand transfer in vitro. Significantly, this same gag sequence was able to fold into a G-quartet monomer, dimer, and tetramer, depending on the cations employed. The pause band at the G-run (nucleotide (nt) 405-409), which was predicted to be involved in forming a G-quartet monomer, diminished with increased HIV-1 nucleocapsid (NC) protein. More NC induced stronger pauses at other G-runs (nt 363-367 and nt 382-384), a region that forms a G-quartet dimer, adhering the two RNA templates. We hypothesized that NC induces the unfolding of the monomeric G-quartet but stabilizes the dimeric interaction. We tested this by inserting a known G-quartet formation sequence, 5'-(UGGGGU)(4)-3', into a relatively structure-free template from the HIV-1 pol gene. Strand transfer assays were performed with cations that either encourage (K(+)) or discourage (Li(+)) G-quartet formation with or without NC. Strikingly, a G-quartet monomer was observed without NC, whereas a G-quartet dimer was observed with NC, both only in the presence of K(+). Moreover, the transfer efficiency of the dimerized template (with K(+) and NC) reached about 90%, approximately 2.5-fold of that of the non-dimerized template. Evidently, template dimerization induced by NC creates a proximity effect, leading to the unique high peak of transfer at the gag recombination hotspot.

Interplay between RNA structure and protein evolution in HIV-1

The genomes of many RNA viruses contain abundant secondary structures that have been shown to be important for understanding the evolution of noncoding regions and synonymous sites. However, the consequences for protein evolution are less well understood. Recently, the secondary structure of the HIV-1 RNA genome has been experimentally determined. Using this information, here we show that RNA structure and proteins do not evolve independently. A negative correlation exists between the extent of base pairing in the genomic RNA and amino acid variability. Relaxed RNA structures may favor the accumulation of genetic variation in proteins and, conversely, sequence changes driven by positive selection at the protein level may disrupt existing RNA structures. We also find that breakage of RNA base pairs might impose a fitness cost to drug resistance mutations in the protease and reverse transcriptase genes, thereby limiting their spread among untreated patients. Characterizing the evolutionary trade-offs between the selective pressures acting at the RNA and protein levels will help us to better understand the variability and evolution of HIV-1.

Single-molecule study of DNA polymerization activity of HIV-1 reverse transcriptase on DNA templates

HIV-1 RT (human immunodeficiency virus-1 reverse transcriptase) is a multifunctional polymerase responsible for reverse transcription of the HIV genome, including DNA replication on both RNA and DNA templates. During reverse transcription in vivo, HIV-1 RT replicates through various secondary structures on RNA and single-stranded DNA (ssDNA) templates without the need for a nucleic acid unwinding protein, such as a helicase. In order to understand the mechanism of polymerization through secondary structures, we investigated the DNA polymerization activity of HIV-1 RT on long ssDNA templates using a multiplexed single-molecule DNA flow-stretching assay. We observed that HIV-1 RT performs fast primer extension DNA synthesis on single-stranded regions of DNA (18.7 nt/s) and switches its activity to slow strand displacement synthesis at DNA hairpin locations (2.3 nt/s). Furthermore, we found that the rate of strand displacement synthesis is dependent on the GC content in hairpin stems and template stretching force. This indicates that the strand displacement synthesis occurs through a mechanism that is neither completely active nor passive: that is, the opening of the DNA hairpin is driven by a combination of free energy released during dNTP (deoxyribonucleotide triphosphate) hydrolysis and thermal fraying of base pairs. Our experimental observations provide new insight into the interchanging modes of DNA replication by HIV-1 RT on long ssDNA templates.

A recombination hot spot in HIV-1 contains guanosine runs that can form a G-quartet structure and promote strand transfer in vitro

The co-packaged RNA genomes of human immunodeficiency virus-1 recombine at a high rate. Recombination can mix mutations to generate viruses that escape immune response. A cell-culture-based system was designed previously to map recombination events in a 459-bp region spanning the primer binding site through a portion of the gag protein coding region. Strikingly, a strong preferential site for recombination in vivo was identified within a 112-nucleotide-long region near the beginning of gag. Strand transfer assays in vitro revealed that three pause bands in the gag hot spot each corresponded to a run of guanosine (G) residues. Pausing of reverse transcriptase is known to promote recombination by strand transfer both in vivo and in vitro. To assess the significance of the G runs, we altered them by base substitutions. Disruption of the G runs eliminated both the associated pausing and strand transfer. Some G-rich sequences can develop G-quartet structures, which were first proposed to form in telomeric DNA. G-quartet structure formation is highly dependent on the presence of specific cations. Incubation in cations discouraging G-quartets altered gel mobility of the gag template consistent with breakdown of G-quartet structure. The same cations faded G-run pauses but did not affect pauses caused by hairpins, indicating that quartet structure causes pausing. Moreover, gel analysis with cations favoring G-quartet structure indicated no structure in mutated templates. Overall, results point to reverse transcriptase pausing at G runs that can form quartets as a unique feature of the gag recombination hot spot.

Unequal distribution of RT-PCR artifacts along the E1-E2 region of Hepatitis C virus

Although viral variability studies have focused traditionally on consensus sequences, the relevance of molecular clone sequences for studying viral evolution at the intra-host level is being increasingly recognized. However, for this approach to be reliable, RT-PCR artifacts do not have to contribute excessively to the observed variability. Molecular clone sequences were obtained from an in vitro transcript to estimate the maximum error rate associated to RT-PCR for the Hepatitis C virus (HCV) E1-E2 region. On average, the frequency of RT-PCR errors was one order of magnitude lower than the level of intra-host genetic variability observed in samples from an HCV outbreak. However, RT-PCR errors were not distributed evenly along the E1-E2 region and were concentrated heavily in the hypervariable region 2 (HVR 2). Although it is concluded that RT-PCR molecular clone sequences are reliable, these results warn against extrapolation of RT-PCR error rates to different genome regions. The data suggest that the RNA sequence context or secondary structure can determine the fidelity of in vitro transcription or reverse transcription. Potentially, these factors might also modify the fidelity of the viral polymerase.

Identification of a preferred region for recombination and mutation in HIV-1 gag

We designed a cell culture-based system to test the hypothesis that recombination events during HIV-1 replication would be more frequent near the dimerization initiation sequence (DIS). A 459-bp region spanning the DIS through the 5'-end of gag was sequenced and analyzed to determine the frequency and distribution of crossover sites. We found a strong preference for recombination events occurring within a 112-nt-long region encompassing the gag AUG (64% of crossovers occurred in this region, compared to 10-14% in surrounding regions with similar lengths). Surprisingly, the region immediately surrounding the DIS was not a preferred site of recombination. Analysis of recombination events using RNA templates transcribed in vitro revealed a preference for crossover sites at the start of the gag coding region, similar to that observed in cell culture. This recombinogenic region was unusually G-rich and promoted extensive pausing by RT in vitro. Template features that induce RT pausing very likely contribute to the observed template switching events in gag during minus-strand synthesis. The region in gag that was a preferred site for recombination also had an approximately 2-fold higher mutation frequency compared to the rest of the region sequenced, but mutations were no more common in recombinant compared to non-recombinant clones, suggesting that recombination events were not mutagenic.

Resistance of virus to extinction on bottleneck passages: study of a decaying and fluctuating pattern of fitness loss

RNA viruses display high mutation rates and their populations replicate as dynamic and complex mutant distributions, termed viral quasispecies. Repeated genetic bottlenecks, which experimentally are carried out through serial plaque-to-plaque transfers of the virus, lead to fitness decrease (measured here as diminished capacity to produce infectious progeny). Here we report an analysis of fitness evolution of several low fitness foot-and-mouth disease virus clones subjected to 50 plaque-to-plaque transfers. Unexpectedly, fitness decrease, rather than being continuous and monotonic, displayed a fluctuating pattern, which was influenced by both the virus and the state of the host cell as shown by effects of recent cell passage history. The amplitude of the fluctuations increased as fitness decreased, resulting in a remarkable resistance of virus to extinction. Whereas the frequency distribution of fitness in control (independent) experiments follows a log-normal distribution, the probability of fitness values in the evolving bottlenecked populations fitted a Weibull distribution. We suggest that multiple functions of viral genomic RNA and its encoded proteins, subjected to high mutational pressure, interact with cellular components to produce this nontrivial, fluctuating pattern.

Molecular basis of fidelity of DNA synthesis and nucleotide specificity of retroviral reverse transcriptases

Reverse transcription involves the conversion of viral genomic RNAinto proviral double-stranded DNA that integrates into the host cell genome. Cellular DNA polymerases replicate the integrated viral DNA and RNA polymerase II transcribes the proviral DNA into RNA genomes that are packaged into virions. Although mutations can be introduced at any of these replication steps, reverse transcriptase (RT) errors play a major role in retroviral mutation. This review summarizes our current knowledge on fidelity of reverse transcriptases. Estimates of retroviral mutation rates or fidelity of retroviral RTs are discussed in the context of the different techniques used for this purpose (i.e., retroviral vectors replicated in culture, misinsertion and mispair extension fidelity assay, etc.). In vitro fidelity assays provide information on the RT's accuracy during the elongation reaction of DNA synthesis. In addition, other steps such as initiation of reverse transcription, or strand transfer, and factors including viral proteins such as Vpr [in the case of the human immunodeficiency virus type 1 (HIV-1)] have been shown to influence fidelity. A comprehensive description of the effect of amino acid substitutions on the fidelity of HIV-1 RT is presented. Published data point to certain dNTP-binding residues, as well as to various amino acids involved in interactions with the template or the primer strand, and to residues in the minor groove-binding track as major components of the fidelity center of retroviral RTs. Implications of these studies include the design of novel therapeutic strategies leading to virus extinction, by increasing the viral mutation rate beyond a tolerable threshold.

Viral RNA mutations are region specific and increased by ribavirin in a full-length hepatitis C virus replication system

High rates of genetic variation ensure the survival of RNA viruses. Although this variation is thought to result from error-prone replication, RNA viruses must also maintain highly conserved genomic segments. A balance between conserved and variable viral elements is especially important in order for viruses to avoid "error catastrophe." Ribavirin has been shown to induce error catastrophe in other RNA viruses. We therefore used a novel hepatitis C virus (HCV) replication system to determine relative mutation frequencies in variable and conserved regions of the HCV genome, and we further evaluated these frequencies in response to ribavirin. We sequenced the 5' untranslated region (5' UTR) and the core, E2 HVR-1, NS5A, and NS5B regions of replicating HCV RNA isolated from cells transfected with a T7 polymerase-driven full-length HCV cDNA plasmid containing a cis-acting hepatitis delta virus ribozyme to control 3' cleavage. We found quasispecies in the E2 HVR-1 and NS5B regions of untreated replicating viral RNAs but not in conserved 5' UTR, core, or NS5A regions, demonstrating that important cis elements regulate mutation rates within specific viral segments. Neither T7-driven replication nor sequencing artifacts produced these nucleotide substitutions in control experiments. Ribavirin broadly increased error generation, especially in otherwise invariant regions, indicating that it acts as an HCV RNA mutagen in vivo. Similar results were obtained in hepatocyte-derived cell lines. These results demonstrate the potential utility of our system for the study of intrinsic factors regulating genetic variation in HCV. Our results further suggest that ribavirin acts clinically by promoting nonviable HCV RNA mutation rates. Finally, the latter result suggests that our replication model may be useful for identifying agents capable of driving replicating virus into error catastrophe.

Resistance to extinction of low fitness virus subjected to plaque-to-plaque transfers: diversification by mutation clustering

Plaque-to-plaque transfers of RNA viruses lead to accumulation of mutations and fitness decrease. To test whether continuing plaque-to-plaque transfers would lead to viral extinction, we have subjected several low fitness foot-and-mouth disease virus (FMDV) clones to up to 130 successive plaque transfers, and have analyzed the evolution of plaque titers and genomic nucleotide sequences. No case of viral extinction could be documented. Some low fitness clones that posses an internal poly(A) tract evaded extinction by modifying the length or base composition of the poly(A) tract. The comparison of entire genomic sequences of FMDV clones at increasing plaque transfer number revealed that mutations accumulated at a uniform rate, and that they were distributed unevenly along the genome. Clusters of mutations were identified at different genomic sites in two plaque transfer lineages. Mutation clustering appears to occur stochastically and could not be related to fixation of compensatory mutations. The results document resistance of viral clones to extinction, and suggest that mutation clustering may be a mechanism of genetic diversification of low fitness virus.

Efficient virus extinction by combinations of a mutagen and antiviral inhibitors

The effect of combinations of the mutagenic base analog 5-fluorouracil (FU) and the antiviral inhibitors guanidine hydrochloride (G) and heparin (H) on the infectivity of foot-and-mouth disease virus (FMDV) in cell culture has been investigated. Related FMDV clones differing up to 10(6)-fold in relative fitness in BHK-21 cells have been compared. Systematic extinction of intermediate fitness virus was attained with a combination of FU and G but not with the mutagen or the inhibitor alone. Systematic extinction of high-fitness FMDV required the combination of FU, G, and H. FMDV showing high relative fitness in BHK-21 cells but decreased replicative ability in CHO cells behaved as a low-fitness virus with regard to extinction mutagenesis in CHO cells. This confirms that relative fitness, rather than a specific genomic sequence, determines the FMDV response to enhanced mutagenesis. Mutant spectrum analysis of several genomic regions from a preextinction population showed a statistically significant increase in the number of mutations compared with virus passaged in parallel in the absence of FU and inhibitors. Also, in a preextinction population the types of mutations that can be attributed to the mutagenic action of FU were significantly more frequent than other mutation types. The results suggest that combinations of mutagenic agents and antiviral inhibitors can effectively drive high-fitness virus into extinction.

Determination of the mutation rate of poliovirus RNA-dependent RNA polymerase

The fidelity of poliovirus RNA-dependent RNA polymerase (3D(pol)) was determined using a system based on the fidelity of synthesis of the alpha-lac gene which codes for a subunit of beta-galactosidase. Synthesis products are screened for mutations by an alpha-complementation assay, in which the protein product from alpha-lac is used in trans to complement beta-galactosidase activity in bacteria that do not express alpha-Lac. Several polymerases have been analyzed by this approach allowing comparisons to be drawn. The assay included RNA synthesis by 3D(pol) on an RNA template that coded for the N-terminal region of alpha-Lac. The product of this reaction was used as a template for a second round of 3D(pol) synthesis and the resulting RNA was reverse transcribed to DNA by MMLV-RT. The DNA was amplified by PCR and inserted into a vector used to transform Escherichia coli. The bacteria were screened for beta-galactosidase activity by blue-white phenotype analysis with white or faint blue colonies scored as errors made during synthesis on alpha-lac. Results showed a mutation rate for 3D(pol) corresponding to approximately 4.5x10(-4) errors per base (one error in approximately 2200 bases). Analysis of mutations showed that base substitutions occurred with greater frequency than deletions and insertions.

Generation of adenovirus vectors devoid of all viral genes by recombination between inverted repeats

Direct or inverse repeated sequences are important functional features of prokaryotic and eukaryotic genomes. Considering the unique mechanism, involving single-stranded genomic intermediates, by which adenovirus (Ad) replicates its genome, we investigated whether repetitive homologous sequences inserted into E1-deleted adenoviral vectors would affect replication of viral DNA. In these studies we found that inverted repeats (IRs) inserted into the E1 region could mediate predictable genomic rearrangements, resulting in vector genomes devoid of all viral genes. These genomes (termed DeltaAd.IR) contained only the transgene cassette flanked on both sides by precisely duplicated IRs, Ad packaging signals, and Ad inverted terminal repeat sequences. Generation of DeltaAd.IR genomes could also be achieved by coinfecting two viruses, each providing one inverse homology element. The formation of DeltaAd.IR genomes required Ad DNA replication and appeared to involve recombination between the homologous inverted sequences. The formation of DeltaAd. IR genomes did not depend on the sequence within or adjacent to the inverted repeat elements. The small DeltaAd.IR vector genomes were efficiently packaged into functional Ad particles. All functions for DeltaAd.IR replication and packaging were provided by the full-length genome amplified in the same cell. DeltaAd.IR vectors were produced at a yield of approximately 10(4) particles per cell, which could be separated from virions with full-length genomes based on their lighter buoyant density. DeltaAd.IR vectors infected cultured cells with the same efficiency as first-generation vectors; however, transgene expression was only transient due to the instability of deleted genomes within transduced cells. The finding that IRs present within Ad vector genomes can mediate precise genetic rearrangements has important implications for the development of new vectors for gene therapy approaches.

The effect of template RNA structure on elongation by HIV-1 reverse transcriptase

Reverse transcription of the RNA genome of retroviruses has to proceed through some highly structured regions of the template. The RNA genome of the human immunodeficiency virus type 1 (HIV-1) contains two hairpin structures within the repeat (R) region at the 5' end of the viral RNA (Fig. 1Fig. 1Template RNA structure of the HIV-1 R region and the position of reverse transcription pause sites. The HIV-1 R region (nucleotides +1/97) encodes two stable RNA structures, the TAR and polyA hairpins [5]. The latter hairpin contains the AAUAAA hexamer motif (marked by a box) that is involved in polyadenylation. The lower panel shows the predicted structures of the wild-type and two mutant forms of the polyA hairpin that were used in this study. Nucleotide substitutions are boxed, deletions are indicated by black triangle. The thermodynamic stability (free energy or DeltaG, in kcal/mol) was calculated according to the Zucker algorithm [71]. The TAR hairpin has a DeltaG of -24.8 kcal/mol. Minus-strand DNA synthesis on these templates was initiated by a DNA primer annealed to the downstream PBS. The position of reverse transcription pause sites observed in this study are summarized. All numbers refer to nucleotide positions on the wild-type HIV-1 transcript. Filled arrows represent stops observed on the wild-type template, and open arrows mark the pause sites that are specific for the structured A-mutant template. The sizes of the arrows correspond to the relative frequency of pausing. Little pausing was observed on the B-mutant template with the destabilized polyA hairpin.). These structures, the TAR and polyA hairpins, fulfil important functions in the viral life cycle. We analyzed the in vitro elongation properties of the HIV-1 reverse transcriptase (RT) enzyme on the wild-type RNA template and mutants thereof with either a stabilized or a destabilized polyA hairpin. Stable RNA structure was found to interfere with efficient elongation of the RT enzyme, as judged by the appearance of pause cDNA products. A direct relation was measured between the stability of template RNA structure and the extent of RT pausing. However, the position of structure-induced pause sites is rather diverse, with significant stops at a position approximately 6 nt ahead of the basepaired stem of the TAR and polyA hairpins. This suggests that the RT enzyme is stalled when its most forward domain contacts the RNA duplex. Addition of the viral nucleocapsid protein (NC) to the in vitro assay was found to overcome such structure-induced RT stops. These results indicate that the RT polymerase has problems penetrating regions of the template with stable RNA structure. This effect was more pronounced at high Mg2+ concentrations, which is known to stabilize RNA secondary structure. Such a structure-induced defect was not apparent in reverse transcription assays performed in virus-infected cells, which is either caused by the NC protein or other components of the virion particle. Thus, retroviruses can use relatively stable RNA structures to control different steps in the viral life cycle without interfering with the process of reverse transcription.

DNA secondary structure effects on DNA synthesis catalyzed by HIV-1 reverse transcriptase

The effect of DNA secondary structure on polymerization catalyzed by human immunodeficiency virus (HIV-1) reverse transcriptase (RT) was studied using a synthetic 66-nucleotide DNA template containing a stable hairpin structure. Four RT pause sites were identified within the first half of the hairpin stem. Additionally, five weak pause sites within the second half of the stem and the loop of the hairpin were identified at low temperatures. These weak pause sites were relocated to the site of the first few stem base pairs of two new hairpins formed due to a change in DNA secondary structure. Each pause site was correlated with a high free energy barrier of melting the stem base pair. Pre-steady state kinetic analysis of single nucleotide incorporation showed that polymerization at each pause site occurred by both a fast phase (10-20 s-1) and a slow phase (0. 02-0.07 s-1) during a single binding event. The reaction amplitudes of the fast phase were small (4-10% of enzyme sites), whereas the amplitudes of the slow phase were large (14-40%) at the pause sites. In contrast, only a single phase with a large reaction amplitude (32-50%) and a fast nucleotide incorporation rate (33-87 s-1) was observed at the non-pause sites. DNA substrates at all sites had similar dissociation rates (0.14-0.29 s-1) and overall binding affinity (16-86 nM). These results suggest that the DNA substrates at pause sites were bound in both productive and non-productive states at the polymerase site of RT. The non-productively bound DNA was slowly converted into a productive state upon melting of the next stem base pair without dissociation of the DNA from RT.

Evolution of envelope sequences from the genital tract and peripheral blood of women infected with clade A human immunodeficiency virus type 1

The development of viral diversity during the course of human immunodeficiency virus type 1 (HIV-1) infection may significantly influence viral pathogenesis. The paradigm for HIV-1 evolution is based primarily on studies of male cohorts in which individuals were presumably infected with a single virus variant of subtype B HIV-1. In this study, we evaluated virus evolution based on sequence information of the V1, V2, and V3 portions of HIV-1 clade A envelope genes obtained from peripheral blood and cervical secretions of three women with genetically heterogeneous viral populations near seroconversion. At the first sample following seroconversion, the number of nonsynonymous substitutions per potential nonsynonymous site (dn) significantly exceeded substitutions at potential synonymous sites (ds) in plasma viral sequences from all individuals. Generally, values of dn remained higher than values of ds as sequences from blood or mucosa evolved. Mutations affected each of the three variable regions of the envelope gene differently; insertions and deletions dominated changes in V1, substitutions involving charged amino acids occurred in V2, and sequential replacement of amino acids over time at a small subset of positions distinguished V3. The relationship among envelope nucleotide sequences obtained from peripheral blood mononuclear cells, plasma, and cervical secretions was evaluated for each individual by both phylogenetic and phenetic analyses. In all subjects, sequences from within each tissue compartment were more closely related to each other than to sequences from other tissues (phylogenetic tissue compartmentalization). At time points after seroconversion in two individuals, there was also greater genetic identity among sequences from the same tissue compartment than among sequences from different tissue compartments (phenetic tissue compartmentalization). Over time, temporal phylogenetic and phenetic structure was detectable in mucosal and plasma viral samples from all three women, suggesting a continual process of migration of one or a few infected cells into each compartment followed by localized expansion and evolution of that population.

Pausing of reverse transcriptase on retroviral RNA templates is influenced by secondary structures both 5' and 3' of the catalytic site

In the most extensive examination to date of the relationship between the pausing of reverse transcrip-tase (RT) and RNA secondary structures, pause events were found to be correlated to inverted repeats both ahead of, and behind the catalytic site in vitro. In addition pausing events were strongly associated with polyadenosine sequences and to a lesser degree diadenosines and monoadenosine residues. Pausing was also inversely proportional to the potential bond strength between the nascent strand and the template at the point of termination, for both mono and dinucleotides. A run of five adenosine and four uridine residues caused most pausing on the HIV-1 template, a region which is the site of much sequence heterogeneity in HIV-1. We propose that homopolyadenosine tracts can act as termination signals for RT in the context of inverted repeats as they do for certain RNA polymerases.

Insights into the interaction between tRNA and primer binding site from characterization of a unique HIV-1 virus which stably maintains dual PBS complementary to tRNA(Gly) and tRNA(His)

Previously our laboratory constructed an HIV-1 which stably maintained a primer binding site (PBS) complementary to tRNA(His) by mutating the region of the provirus within U5 postulated to interact with the anticodon of tRNA(His) (J. Wakefield, S-M Kang, and C. D. Morrow, 1996, J. Virol, 70, 966-975). From the analysis of the virus obtained after long-term culture, we identified an unusual proviral DNA in which the U5-PBS region contained a dual PBS complementary to tRNA(Gly) and tRNA(His), respectively, separated by a 21-nucleotide intervening sequence. To determine if this U5-PBS region containing the dual PBS would give rise to an infectious virus, the mutant U5-PBS containing the dual PBS was subcloned into an infectious HIV-1 proviral clone, pHXB2; the resultant proviral DNA was designated as pHXB2(Gly-His). Transfection of pHXB2(Gly-His) into cells gave rise to infectious virus. Analysis of the U5-PBS region revealed that the virus stably maintained the dual PBS rather than revert back to the wild-type PBS. In addition to genomes with the PBS complementary to tRNA(Gly) and tRNA(His), proviral genomes were identified after extended in vitro culture which contained dual PBS complementary to tRNA(Gly) and tRNA(Phe). To determine which PBS could be used for reverse transcription, we utilized an endogenous reverse transcription/PCR method which could discriminate (based on molecular size of the products) between the minus strand DNA initiated from the two PBSs. The results of this assay demonstrated that either the PBS complementary to tRNA(Gly) or tRNA(His) could be used for the initiation of reverse transcription. The results of our study highlight the complex interrelationship between U5-PBS and primer tRNA required for positioning the tRNA at the PBS and provides new insights into how the tRNA primer used to initiate reverse transcription is selected.

RNA virus mutations and fitness for survival

RNA viruses exploit all known mechanisms of genetic variation to ensure their survival. Distinctive features of RNA virus replication include high mutation rates, high yields, and short replication times. As a consequence, RNA viruses replicate as complex and dynamic mutant swarms, called viral quasispecies. Mutation rates at defined genomic sites are affected by the nucleotide sequence context on the template molecule as well as by environmental factors. In vitro hypermutation reactions offer a means to explore the functional sequence space of nucleic acids and proteins. The evolution of a viral quasispecies is extremely dependent on the population size of the virus that is involved in the infections. Repeated bottleneck events lead to average fitness losses, with viruses that harbor unusual, deleterious mutations. In contrast, large population passages result in rapid fitness gains, much larger than those so far scored for cellular organisms. Fitness gains in one environment often lead to fitness losses in an alternative environment. An important challenge in RNA virus evolution research is the assignment of phenotypic traits to specific mutations. Different constellations of mutations may be associated with a similar biological behavior. In addition, recent evidence suggests the existence of critical thresholds for the expression of phenotypic traits. Epidemiological as well as functional and structural studies suggest that RNA viruses can tolerate restricted types and numbers of mutations during any specific time point during their evolution. Viruses occupy only a tiny portion of their potential sequence space. Such limited tolerance to mutations may open new avenues for combating viral infections.

Mutations of the HIV type 1 V3 loop under selection pressure with neutralizing monoclonal antibody NM-01

Variants of human immunodeficiency virus type 1 (HIV-1) were selected for resistance to the neutralizing monoclonal antibody (MAb) NM-01. MAb NM-01 recognizes the center of the third hypervariable domain (V3 loop) of the envelope gp120, and neutralizes diverse HIV-1 strains. In the continuous presence of MAb NM-01, transmission and propagation of molecularly cloned HIV-1 were performed in vitro to isolate escape variants. The polymerase chain reaction-single-strand conformation polymorphism and sequence analyses of these variants indicated that the antigenic change against MAb NM-01 is due to a single base substitution resulting in one amino acid interchange within the recognition site of MAb NM-01 in the V3 loop. Mutational analyses also demonstrated a nonrandom event of variability and the existence of mutational hot spots in the V3 loop. The bias of variability could be interpreted by the specificity of error-prone replication by HIV-1 reverse transcriptase. Furthermore, the results suggest that distribution of mutability might correlate closely with the stability of the secondary structure of RNA encoding the V3 loop region.

Negative effects of chemical mutagenesis on the adaptive behavior of vesicular stomatitis virus

Changes in adaptability of vesicular stomatitis virus (VSV) upon treatment with chemical mutagens have been investigated. Results showed no improvement in virus viability or adaptability at any given level of mutagenesis. In fact, increasing inhibition of virus production and adaptability was observed with increasing levels of mutagenesis. This was true for all tested VSV variants replicating either in changing or constant host cell environments. Results also showed that mutagen-treated RNA virus populations which had undergone severe fitness declines were able to recover lost fitness completely after several large-population passages in BHK21, cells. The present findings illustrate the highly optimized states of RNA viruses and their potential to adapt readily. These results are significant for the possible development of specific antiviral agents designed to be mutagenic.

Wild-type and mutant HIV type 1 nucleocapsid proteins increase the proportion of long cDNA transcripts by viral reverse transcriptase

HIV-1 nucleocapsid, p7, contains two retroviral zinc fingers, which are both necessary for efficient packaging of genomic RNA and infectivity. The nucleocapsid protein is bound tightly to genomic RNA in the mature virion. In this study, the effect of p7 on polymerization of nascent cDNA by viral reverse transcriptase (RT) was examined. An 874-base RNA of HIV-1 was synthesized and used as a template in RT assays with varying concentrations of intact p7, mutants of p7 that have transposed or repeated zinc fingers, and several different peptides that represent various structural regions of p7. Results indicate that at greater than or equal to 50% saturation of p7-binding sites, with p7, there is up to a 90% reduction in total cDNA synthesis, as measured by nucleotide incorporation. However, the cDNA products that are made are almost exclusively full length. Three zinc finger mutants exhibited effects similar to those of wild-type p7. N-terminal and C-terminal halves of p7 inhibited total nucleotide incorporation, but also inhibited synthesis of long cDNA products by RT. In the absence of p7 an array of short transcripts (< 200 bases) was produced by RT. These studies show that full-length p7 is necessary to increase the proportion of long cDNA transcripts produced by RT. The relative position of the two zinc fingers is not critical for this effect.

Mechanisms of plant virus evolution

Plant viruses utilize several mechanisms to generate the large amount of genetic diversity found both within and between species. Plant RNA viruses and pararetroviruses probably have highly error prone replication mechanisms, that result in numerous mutations and a quasispecies nature. The plant DNA viruses also exhibit diversity, but the source of this is less clear. Plant viruses frequently use recombination and reassortment as driving forces in evolution, and, occasionally, other mechanisms such as gene duplication and overprinting. The amount of variation found in different species of plant viruses is remarkably different, even though there is no evidence that the mutation rate varies. The origin of plant viruses is uncertain, but several possible theories are proposed. The relationships between some plant and animal viruses suggests a common origin, possibly an insect virus. The propensity for rapid adaptation makes tracing the evolutionary history of viruses difficult, and long term control of virus disease nearly impossible, but it provides an excellent model system for studying general mechanisms of molecular evolution.

Misincorporation by HIV-1 reverse transcriptase promotes recombination via strand transfer synthesis

Genome heterogeneity in retroviruses derives from poor fidelity of the reverse transcriptase (RT) and recombination via RT-catalyzed strand transfer synthesis. RTs lack proofreading ability, and they proficiently extend primers with mismatched termini. Recombination reactions carried out in vitro are accompanied by a high frequency of base substitution errors, suggesting a relationship. Here we provide evidence that misincorporation during RNA-directed DNA synthesis promotes strand transfer recombination. Experiments involved measurement of DNA synthesis, RNase H-directed cleavage, and strand transfer synthesis from preformed mismatched primers on RNA templates by human immunodeficiency virus (HIV) RT in vitro. A significant pause in synthesis occurred from a G(primer). rA(template) mismatch compared to the synthesis from a correctly paired (T.A) primer. The misincorporation-induced pause allowed an unusually large area of RT-RNase H-directed cleavage of the template RNA beneath the primer. Strand transfer to an acceptor molecule with sequence identical to the template RNA was about 50% more efficient than if the primer had had a correctly paired terminus. Overall transfer was measured over a large region of homology. Assuming that enhanced transfer occurs primarily at the site of the mismatch, the actual increase in transfer at that site must have been 1-2 orders of magnitude. Inclusion of a different acceptor molecule with complete complementarity to the originally mismatched 3' primer terminus resulted in an additional 2-fold increase in strand transfer efficiency. Overall, these results suggest the mechanism by which misincorporation during minus strand DNA synthesis in retroviral replication would promote high frequency recombination.

Nucleotide substitution patterns can predict the requirements for drug-resistance of HIV-1 proteins

The enzyme reverse transcriptase (RT) plays a fundamental role in the replication of the human immunodeficiency virus type 1 (HIV-1) and several antiviral agents that target this key enzyme have been developed. Unfortunately, treatment of patients with RT inhibitors results in the appearance of drug-resistant variants with specific mutations in the RT protein. We hypothesized that if "difficult' resistance mutations (e.g. transversions/double-hits) are consistently observed at certain positions, it is likely that "easier' nucleotide substitutions (transitions/single-hits) at that codon do not result in a drug-resistant and/or active RT enzyme. In this study, we examined codon changes involved in RT drug resistance against nucleoside and non-nucleoside inhibitors and listed all easier substitutions, which apparently were not selected, either due to reduced enzyme RT activity or lack of drug resistance. These predictions on the requirements for resistance were confirmed by published mutational data on RT variants. We also propose that differences in mutation type can explain the order of appearance of substitutions in case multiple amino acid changes are required for optimal fitness. Differences in mutation pattern have been reported for drug-resistant HIV-1 variants selected in tissue culture compared with variants found in treated patients. In contrast to the in vivo situation, a relatively small population size is handled in in vitro tissue culture systems and this may limit the chances of creating a resistance mutation. Indeed, inspection of the codon changes indicates that the in vitro culture system is more strongly biased towards the relatively easy nucleotide substitutions. These results suggest that the nucleotide substitution pattern can provide important information on RT drug resistance.

Mechanisms of retroviral mutation

Retroviruses, like other RNA viruses, mutate at very high rates (0.05-1 mutations per genome per replication cycle) and exist as complex genetically heterogeneous populations ('quasispecies') that are ever changing. De novo mutations are generated by inherently error-prone steps in the retroviral life cycle that introduce base substitutions, frame shifts, genetic rearrangements and hypermutations.

Efficient extension of a misaligned tRNA-primer during replication of the HIV-1 retrovirus

The human immunodeficiency virus (HIV) and other retroviruses show extensive genomic variation, which is primarily due to error-prone replication by the viral reverse transcriptase (RT) enzymes. RT errors include misincorporation with subsequent extension of the mismatched terminal base, and extension of realigned primer-template duplexes. Whereas both RT-mediated mechanisms have been extensively studied in vitro, almost no in vivo experiments have been performed. In this work, we analyzed the ability of HIV-1 RT to extend a misaligned tRNA(Lys3) primer in vivo. This tRNA binds with its 3'-terminal 18 nt to a complementary sequence in the viral genome, referred to as the primer-binding site (PBS). We constructed a series of mutant viral genomes with small insertions or deletions in the PBS sequence, resulting in misalignment of the tRNA primer. Extension of the misaligned primer did occur with reasonable efficiency for some of the mutants, resulting in reversion to the wild-type viral sequence. The infectivity and reversion frequency of the PBS mutants is therefore a measure of the efficiency of extending a misaligned primer in vivo. Using virion-derived primer-template complexes, we also measured the tRNA-priming efficiency in vitro. The combined results show that HIV-1 RT can elongate a misaligned primer and that the efficiency of primer extension is determined by the extent of the mismatch.