Point mutant frequencies in the pol gene of human immunodeficiency virus type 1 are two- to threefold lower than those of env

Characterization of minority HIV-1 drug resistant variants in the United Kingdom following the verification of a deep sequencing-based HIV-1 genotyping and tropism assay

BACKGROUND

The widespread global access to antiretroviral drugs has led to considerable reductions in morbidity and mortality but, unfortunately, the risk of virologic failure increases with the emergence, and potential transmission, of drug resistant viruses. Detecting and quantifying HIV-1 drug resistance has therefore become the standard of care when designing new antiretroviral regimens. The sensitivity of Sanger sequencing-based HIV-1 genotypic assays is limited by its inability to identify minority members of the quasispecies, i.e., it only detects variants present above ~ 20% of the viral population, thus, failing to detect minority variants below this threshold. It is clear that deep sequencing-based HIV-1 genotyping assays are an important step change towards accurately monitoring HIV-infected individuals.

METHODS

We implemented and verified a clinically validated HIV-1 genotyping assay based on deep sequencing (DEEPGEN™) in two clinical laboratories in the United Kingdom: St. George's University Hospitals Healthcare NHS Foundation Trust (London) and at NHS Lothian (Edinburgh), to characterize minority HIV-1 variants in 109 plasma samples from ART-naïve or -experienced individuals.

RESULTS

Although subtype B HIV-1 strains were highly prevalent (44%, 48/109), most individuals were infected with non-B subtype viruses (i.e., A1, A2, C, D, F1, G, CRF02_AG, and CRF01_AE). DEEPGEN™ was able to accurately detect drug resistance-associated mutations not identified using standard Sanger sequencing-based tests, which correlated significantly with patient's antiretroviral treatment histories. A higher proportion of minority PI-, NRTI-, and NNRTI-resistance mutations was detected in NHS Lothian patients compared to individuals from St. George's, mainly M46I/L and I50 V (associated with PIs), D67 N, K65R, L74I, M184 V/I, and K219Q (NRTIs), and L100I (NNRTIs). Interestingly, we observed an inverse correlation between intra-patient HIV-1 diversity and CD4⁺ T cell counts in the NHS Lothian patients.

CONCLUSIONS

This is the first study evaluating the transition, training, and implementation of DEEPGEN™ between three clinical laboratories in two different countries. More importantly, we were able to characterize the HIV-1 drug resistance profile (including minority variants), coreceptor tropism, subtyping, and intra-patient viral diversity in patients from the United Kingdom, providing a rigorous foundation for basing clinical decisions on highly sensitive and cost-effective deep sequencing-based HIV-1 genotyping assays in the country.

Deep sequencing: becoming a critical tool in clinical virology

Population (Sanger) sequencing has been the standard method in basic and clinical DNA sequencing for almost 40 years; however, next-generation (deep) sequencing methodologies are now revolutionizing the field of genomics, and clinical virology is no exception. Deep sequencing is highly efficient, producing an enormous amount of information at low cost in a relatively short period of time. High-throughput sequencing techniques have enabled significant contributions to multiples areas in virology, including virus discovery and metagenomics (viromes), molecular epidemiology, pathogenesis, and studies of how viruses to escape the host immune system and antiviral pressures. In addition, new and more affordable deep sequencing-based assays are now being implemented in clinical laboratories. Here, we review the use of the current deep sequencing platforms in virology, focusing on three of the most studied viruses: human immunodeficiency virus (HIV), hepatitis C virus (HCV), and influenza virus.

Variation of human immunodeficiency virus type-1 reverse transcriptase within the simian immunodeficiency virus genome of RT-SHIV

RT-SHIV is a chimera of simian immunodeficiency virus (SIV) containing the reverse transcriptase (RT)-encoding region of human immunodeficiency virus type 1 (HIV-1) within the backbone of SIV_mac239. It has been used in a non-human primate model for studies of non-nucleoside RT inhibitors (NNRTI) and highly active antiretroviral therapy (HAART). We and others have identified several mutations that arise in the "foreign" HIV-1 RT of RT-SHIV during in vivo replication. In this study we catalogued amino acid substitutions in the HIV-1 RT and in regions of the SIV backbone with which RT interacts that emerged 30 weeks post-infection from seven RT-SHIV-infected rhesus macaques. The virus set points varied from relatively high virus load, moderate virus load, to undetectable virus load. The G196R substitution in RT was detected from 6 of 7 animals at week 4 post-infection and remained in virus from 4 of 6 animals at week 30. Virus from four high virus load animals showed several common mutations within RT, including L74V or V75L, G196R, L214F, and K275R. The foreign RT from high virus load isolates exhibited as much variation as that of the highly variable envelope surface glycoprotein, and 10-fold higher than that of the native RT of SIV_mac239. Isolates from moderate virus load animals showed much less variation in the foreign RT than the high virus load isolates. No variation was found in SIV_mac239 genes known to interact with RT. Our results demonstrate substantial adaptation of the foreign HIV-1 RT in RT-SHIV-infected macaques, which most likely reflects selective pressure upon the foreign RT to attain optimal activity within the context of the chimeric RT-SHIV and the rhesus macaque host.

Detection of SARS coronavirus in patients with severe acute respiratory syndrome by conventional and real-time quantitative reverse transcription-PCR assays

BACKGROUND

A novel coronavirus (CoV) was recently identified as the agent for severe acute respiratory syndrome (SARS). We compared the abilities of conventional and real-time reverse transcription-PCR (RT-PCR) assays to detect SARS CoV in clinical specimens.

METHODS

RNA samples isolated from nasopharyngeal aspirate (NPA; n = 170) and stool (n = 44) were reverse-transcribed and tested by our in-house conventional RT-PCR assay. We selected 98 NPA and 37 stool samples collected at different times after the onset of disease and tested them in a real-time quantitative RT-PCR specific for the open reading frame (ORF) 1b region of SARS CoV. Detection rates for the conventional and real-time quantitative RT-PCR assays were compared. To investigate the nature of viral RNA molecules in these clinical samples, we determined copy numbers of ORF 1b and nucleocapsid (N) gene sequences of SARS CoV.

RESULTS

The quantitative real-time RT-PCR assay was more sensitive than the conventional RT-PCR assay for detecting SARS CoV in samples collected early in the course of the disease. Real-time assays targeted at the ORF 1b region and the N gene revealed that copy numbers of ORF 1b and N gene sequences in clinical samples were similar.

CONCLUSIONS

NPA and stool samples can be used for early diagnosis of SARS. The real-time quantitative RT-PCR assay for SARS CoV is potentially useful for early detection of SARS CoV. Our results suggest that genomic RNA is the predominant viral RNA species in clinical samples.

Characterization of the mutant spectra of a fish RNA virus within individual hosts during natural infections

Infectious hematopoietic necrosis virus (IHNV) is an RNA virus that causes significant mortalities of salmonids in the Pacific Northwest of North America. RNA virus populations typically contain genetic variants that form a heterogeneous virus pool, referred to as a quasispecies or mutant spectrum. This study characterized the mutant spectra of IHNV populations within individual fish reared in different environmental settings by RT-PCR of genomic viral RNA and determination of partial glycoprotein gene sequences of molecular clones. The diversity of the mutant spectra from ten in vivo populations was low and the average mutation frequencies of duplicate populations did not significantly exceed the background mutation level expected from the methodology. In contrast, two in vitro populations contained variants with an identical mutational hot spot. These results indicated that the mutant spectra of natural IHNV populations is very homogeneous, and does not explain the different magnitudes of genetic diversity observed between the different IHNV genogroups. Overall the mutant frequency of IHNV within its host is one of the lowest reported for RNA viruses.

Common genetic arrangements among human immunodeficiency virus type 1 subtype A and D recombinant genomes vertically transmitted in Tanzania

Human immunodeficiency virus type 1 (HIV-1) subtypes A, C, and D are cocirculating in Tanzania, and large numbers of recombinant genomes have been reported from this region. Here we describe full-length sequences of six unlinked HIV-1 subtype A and D recombinants. The samples came from newborns, indicating that the recombination patterns were vertically transmitted and were functionally competent. All six genomes had different recombination patterns with one to eight cross-over points frequently located at the beginning or end of functionally defined regions. In five of the six viruses most of gag, pol, tat, and rev and the intracytoplasmic domain of gp41 were classified as subtype D. In all but one genome, the external domain of gp41 and the majority of gp120 belonged to subtype A. A recombination site common to four of the six genomes was located at the transmembrane domain of gp41, at the end of the rev response element. The identification of subtype patterns among intersubtype recombinant genomes from recently infected individuals may reveal genetic determinants of improved viral fitness or advantage for transmission.

Survey of reverse transcriptase from the heterosexual epidemic of human immunodeficiency virus type 1 CRF01_AE in Thailand from 1990 to 2000

Genetic diversity of the HIV-1 envelope gene has shown a steady increase over time in the Thai and other regional epidemics. A serial survey of subtype CRF01_AE polymerase gene (RT) diversity in Thailand was performed, using 48 novel and 15 reported sequences covering the period 1990--2000. These sequences were gathered from individuals whose sole risk factor for infection was heterosexual contact. By contrast to envelope, diversity was low and, despite a 40% increase early in the epidemic, has remained static since 1996. These results indicate that epidemic HIV-1 may be constrained within defined limits of genetic diversity at least in some genomic regions.

Unusual distribution of mutations associated with serial bottleneck passages of human immunodeficiency virus type 1

Repeated bottleneck passages result in fitness losses of RNA viruses. In the case of human immunodeficiency virus type 1 (HIV-1), decreases in fitness after a limited number of plaque-to-plaque transfers in MT-4 cells were very drastic. Here we report an analysis of entire genomic nucleotide sequences of four HIV-1 clones derived from the same HIV-1 isolate and their low-fitness progeny following 7 to 15 plaque-to-plaque passages. Clones accumulated 4 to 28 mutations per genome, with dominance of A --> G and G --> A transitions (57% of all mutations) and 49% nonsynonymous replacements. One clone-but not three sibling clones-showed an overabundance of G --> A transitions, evidencing the highly stochastic nature of some types of mutational bias. The distribution of mutations along the genome was very unusual in that mutation frequencies in gag were threefold higher than in env. Particularly striking was the complete absence of replacements in the V3 loop of gp120, confirmed with partial nucleotide sequences of additional HIV-1 clones subjected to repeated bottleneck passages. The analyses revealed several amino acid replacements that have not been previously recorded among natural HIV-1 isolates and illustrate how evolution of an RNA virus genome, with regard to constant and variable regions, can be profoundly modified by alterations in population dynamics.

Identification of biased amino acid substitution patterns in human immunodeficiency virus type 1 isolates from patients treated with protease inhibitors

Human immunodeficiency virus type 1 (HIV-1) amino acid substitutions observed during antiretroviral drug therapy may be caused by drug selection, non-drug-related evolution, or sampling error introduced by the sequencing process. We analyzed HIV-1 sequences from 371 untreated patients and from 178 patients receiving a single protease inhibitor. Amino acid substitution patterns during treatment were compared with inferred substitution patterns arising evolutionarily without treatment. Our results suggest that most treatment-associated amino acid substitutions are caused by selective drug pressure, including substitutions not previously associated with drug resistance.

Human immunodeficiency virus type 1 genetic evolution in patients with prolonged suppression of plasma viremia

Treatment of human immunodeficiency virus type 1 (HIV-1)-infected patients with combination drug regimens results in a reduction of plasma viral load to levels below the limit of detection. To investigate the genomic fluctuations in HIV-1 populations from long-term responders to antiviral therapies we analyzed the viral sequence evolution of env and pol genes from sequential peripheral blood mononuclear cell (PBMC) DNA samples of three infected patients. Analyses of sequences covering the V3 and flanking env regions obtained from blood samples at the beginning of the therapy and at 14 or 24 months from baseline revealed that HIV-1 quasispecies continue to evolve in the three patients following combination antiretroviral therapy. Minor drug-resistant mutant subpopulations were also searched for and found in one patient. Interestingly, no minor resistant subpopulations were found in the other two patients despite the fact that they showed evidence of ongoing viral replication. Finally, the genetic analysis of the env gene shows a reduction in PBMC env viral population diversity after long-term response to the therapy in all the patients analyzed.

Genetic variation and susceptibilities to protease inhibitors among subtype B and F isolates in Brazil

The genetic variation of the human immunodeficiency virus type 1 (HIV-1) protease gene (prt) permits the classification of HIV-1 strains into five distinct protease subtypes, which follow the gag subtyping patterns. The susceptibilities of non-B-subtype strains to protease inhibitors (PIs) and other antiretroviral drugs remain largely unknown. Subtype F is the main non-B strain contributing to the Brazilian epidemic, accounting for 15 to 20% of these infections. In this work, we report the findings on 81 isolates from PI-naive Brazilian patients collected between 1993 and 1997. In addition, the relevant PI resistance mutations and their phenotypes were determined in vitro for 15 of these patients (B = 9 and F = 6). Among these, the subtype F samples evidenced high sensitivities in vitro to ritonavir and indinavir, with MICs at which 50 and 90% of the isolates are inhibited similar to those of both the Brazilian and the U.S. subtype B isolates. Analysis of the 81 Brazilian prt sequences demonstrated that the subtype F consensus sequence differs from the U.S. and Brazilian subtype B consensus in eight positions (I15V, E35D, M36I, R41K, R57K, Q61N, L63P, and L89M). The frequency of critical PI resistance substitutions (amino acid changes D30N, V82A/F/T, I84V, N88D, and L90M) among Brazilian isolates is very low (mean, 2.5%), and the associated secondary substitutions (amino acid positions 10L, 20K, 36M, 46M, 48G, 54I, 63P, 71A, and 77A) are infrequent. These observations document the relative rarity of resistance to PIs in the treatment of patients infected with HIV-1 subtype F in South America.

Sequence analysis of the polymerase domain of HIV-1 reverse transcriptase in naive and zidovudine-treated individuals reveals a higher polymorphism in alpha-helices as compared with beta-strands

We report a statistical analysis of genetic heterogeneity of the reverse transcriptase (RT)-coding region of human immunodeficiency virus type 1. Both newly determined sequences and sequences contained in the data banks have been examined. For the calculations, the viral samples and the regions within the RT molecule were divided in two groups. The viral samples were split into those from patients not subjected to antiretroviral therapy and those from patients treated with zidovudine (AZT, 3'-azido-3'-deoxythymidine) alone or in combination with other RT inhibitors. The RT-coding region was divided into segments encoding beta-strands and segments encoding alpha-helices. A significantly lower heterogeneity was observed in beta-strands relative to the alpha-helix coding segments. Application of the D test of Tajima has provided evidence of operation of negative (or purifying) selection in sequences from viruses of patients not subjected to antiretroviral treatment as well as in treated patients. In the group of untreated individuals, regions encoding beta-strands are subjected to stronger negative selection than those encoding alpha-helices. It is likely that the observed differences reflect stronger functional constraints in beta-strands than in alpha-helices of RT.

Analysis of pol gene heterogeneity, viral quasispecies, and drug resistance in individuals infected with group O strains of human immunodeficiency virus type 1

Nucleotide sequences of the reverse transcriptase (RT) coding region have been compared in four new human immunodeficiency virus type 1 (HIV-1) group O isolates. Phylogenetic analysis of this pol region highlights a cluster of these four HIV-1 group O sequences with seven other group O isolates (5% intracluster nucleotide sequence diversity) similar to clusters classified as subtypes in HIV-1 group M (an average of 4.9% intrasubtype sequence diversity). Based on these analyses, this group O cluster has been designated subtype A-O. A longitudinal study of a heterosexual couple infected with group O (ESP1 and ESP2) allowed a detailed analysis of RT sequences (amino acids 28 to 219). Directed evolution and a slightly higher mutation frequency was observed in the RT sequences of patient ESP2, treated with antiretroviral drugs, than that from the untreated patient ESP1. Antiretroviral treatment also selected for specific substitutions, M184V and T215Y in the RT coding region, conferring resistance to 3'-dideoxy-3'-thiacytidine and zidovudine, respectively. A Gly98 to Glu RT substitution identified in the treated patient suggests a possible reversion of a nonnucleoside RT inhibitor-resistant phenotype. Using RT clones from this longitudinal study, both heteroduplex tracking assay and cloning-sequencing techniques were employed for an extensive genetic analysis of pol gene quasispecies. Amino acid substitutions (i.e., Phe-77 to Leu, Lys-101 to Glu, and Val-106 to Iso) associated with antiretroviral resistance were identified in RT clones from HIV-1 group O-infected patients not subjected to drug therapy or treated with unrelated drugs. Finally, phylogenetic relationships between RT clones of the treated ESP2 patient and those of the untreated ESP1 patient show how drug pressure can direct evolution of viral pol gene quasispecies independently of direct drug-resistant substitutions.

LTR and tat variability of HIV-1 isolates from patients with divergent rates of disease progression

The genetic heterogeneity and transcription activity of the human immunodeficiency virus type 1 (HIV-1) LTR region and tat gene have been examined. Comparison involved the relevant genomic regions of viruses isolated from twenty long-term survivors and from ten typical progressors. No significant differences were observed in mutation frequencies among the two groups, although there was a significant higher proportion of synonymous substitutions in the tat gene of viruses from typical progressors. Four LTR sequences showed an insertion of 20-31 residues at the junction between the LTR Nef-coding and the LTR noncoding region. Neither these insertions nor other genetic changes found in these sequences affected the LTR transcription function, as measured in transient expression assays using transfection of both established cell lines and peripheral blood lymphocytes with plasmid DNA. The results did not allow the association of structural or functional alterations in LTR or tat with a degree of disease progression. The results reinforce the concepts of complexity of HIV-1 evolution in infected individuals, and the multifactorial nature of progression to AIDS.

Human immunodeficiency virus type 1 subtype F reverse transcriptase sequence and drug susceptibility

We sequenced and phylogenetically analyzed the reverse transcriptase (RT) regions of the pol genes of 14 human immunodeficiency virus type 1 (HIV-1) isolates from Romanian patients, which were classified as subtype F on the basis of env gene structure. The RT sequences showed that the strains clustered phylogenetically and were equidistant from other HIV-1 subtypes as shown by the neighbor-joining and maximum-likelihood methods, allowing us to define HIV-1 subtype F according to the pol classification. The subtype F RT sequences differed from reported group M RT sequences by 10.94% (for nucleotides) and 7.6% (for amino acids). Phenotypic analysis of subtype F susceptibility to three classes of antiretroviral compounds showed an increase in the 50% inhibitory concentration of the tetrahydroimidazo[4,5,1-jk] [1,4]-benzodiazepin-2-(1H)-one and -thione (TIBO) derivate R82913 for one strain which was naturally resistant to this compound. This first report of subtype F pol sequences confirms the perfect correlation between the phylogenetic positions determined by env and pol analyses and suggests that virus variability might influence the efficacy of antiretroviral treatments. This finding warrants a global evaluation of the phenotypic and genotypic susceptibility of HIV-1 subtypes to antiretroviral drugs.

Switch to unusual amino acids at codon 215 of the human immunodeficiency virus type 1 reverse transcriptase gene in seroconvertors infected with zidovudine-resistant variants

Sequences of the human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) domain were determined by direct sequencing of HIV-1 RNA in successive plasma samples from eight seroconverting patients infected with virus bearing the T215Y/F amino acid substitution associated with zidovudine (ZDV) resistance. At baseline, additional mutations associated with ZDV resistance were detected. Three patients had the M41L amino acid change, which persisted. Two patients had both the D67N and the K70R amino acid substitutions; reversion to the wild type was seen at both positions in one of these patients and at codon 70 in the other one. Reversion to the wild type at codon 215 was observed in only one of eight patients. Unusual amino acids, such as aspartic acid (D) and cysteine (C), appeared at position 215 in four patients during follow-up. These variants isolated by coculturing were sensitive to ZDV. Overgrowth of these variants suggests that they have better fitness than the original T215Y variant. Intraindividual nucleoside substitutions over time were 10 times more frequent in codons associated with ZDV resistance (41, 67, 70, 215, and 219) than in other codons of the RT domain. The predominance of nonsynonymous substitutions observed over time suggests that most changes reflect adaptation of the RT function. The variance in sequence evolution observed among patients, in particular at codon 215, supports a role for chance in the evolution of the RT domain.

Nonsynonymous mutations within the human immunodeficiency virus type 1 p17 gene are clustered to sequences binding to the host human leukocyte antigen class I molecules

We have analyzed the relation between intrapatient variabilities of the p17 gene and the location of known host p17 cytotoxic T lymphocyte (CTL) epitopes in five patients infected with human immunodeficiency virus type 1 (HIV-1). All patients were typed with respect to the human leukocyte antigen (HLA) class I type. One to seven previously fine-mapped p17 CTL epitopes corresponded to the HLA class I restriction elements of each patient. An average of 28+/-16% of the p17 gene of each patient encoded CTL epitopes corresponding to the HLA restriction elements of the host. Twenty full-length p17 gene clones were sequenced from each patient. The intrapatient homology between the p17 sequences ranged from 96.4 to 98.9%. The interpatient homology between the consensus sequences of each patient ranged from 83.1 to 91.6%. A total of 246 nucleotide differences within the 100 p17 clones was noted. Fifteen (16%) of 96 synonymous substitutions were found within host CTL epitopes, whereas 72 (48%) of 150 nonsynonymous nucleotide changes were found within CTL epitopes corresponding to the HLA restriction elements of the host (p < 0.0001; Fisher's exact test). Subsequently, variable residues indicating the evolution of at least two major p17 species (i.e., >20% of the clones) were determined to be more common at positions contained within these CTL epitopes (p < 0.01). The present data suggest that the evolution of the p17 gene is influenced by contact areas with the host HLA class I molecules.

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.

Characterization of the reverse transcriptase of a human immunodeficiency virus type 1 group O isolate

The catalytic properties and sensitivity to different inhibitors have been determined for the reverse transcriptase (RT) of group O human immunodeficiency virus type 1 (HIV-1). The RT-coding region was cloned from a new HIV-1 group O isolate from Spain, expressed in Escherichia coli, and purified by affinity chromatography. This new RT showed 79% amino acid sequence identity with the corresponding enzyme of group M subtype B strain BH10. The two enzymes showed very similar kinetics of RNA-dependent DNA polymerization using homopolymeric template-primers and RNase H specific activity. Inhibitor sensitivity to ddTTP and 3'-azido-2',3'-dideoxythymidine triphosphate (AZTTP) was also similar for both enzymes. However, the two enzymes differed dramatically in their sensitivity to several inhibitors. While the RT of the BH10 isolate was sensitive to nevirapine and loviride (IC50 ranged from 0.16 to 8.2 microM, depending on the substrates used), the enzyme of the Spanish HIV-1 group O isolate showed high-level resistance to those compounds (IC50 > 200 microM). The amino acid sequence of the RT of group O HIV-1 contains three amino acids (Cys-181, Glu-179, and Gly-98), which are found in group M subtype B strains resistant to nonnucleoside RT inhibitors. The recombinant group O HIV-1 RT should be useful for studies aimed at discovering and designing drugs directed toward group O isolates of HIV-1.

Sequence and drug susceptibility of subtype C reverse transcriptase from human immunodeficiency virus type 1 seroconverters in Zimbabwe

Naturally occurring human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) variability has implications for the success of antiretroviral therapy. We determined the sequence of the polymerase-coding region of RT from virus isolates from 12 Zimbabwean individuals recently infected with HIV-1. The 12 RT sequences differed from the consensus B RT sequence at 10.5% of nucleotides and 5.8% of amino acids. Susceptibility testing of five isolates to zidovudine, didanosine, lamivudine, and nevirapine demonstrated susceptibilities similar to those of wild-type subtype B isolates. Phylogenetic analysis of 40 HIV-1 RT sequences, including the 12 Zimbabwean subtype C sequences, 11 subtype B sequences, and the 17 remaining published non-subtype B sequences showed sufficient intrasubtype RT sequence variation to differentiate subtype A, B, C, and D isolates. Five recently reported subtype C RT sequences from India grouped with the Zimbabwean subtype C sequences but had significantly less intraisolate sequence variation. Both intra- and intersubtype RT comparisons were notable for extraordinarily high ratios of synonymous to nonsynonymous differences. Although substitutions in the HIV-1 RT gene are limited by functional constraints, variation between RT sequences demonstrates phylogenetic relationships that parallel env and gag gene variation.