ANTIMICROBIAL AGENTS AND CHEMOTHERAPY, Mar. 2000, p. 568 573 Vol. 44, No. 3 0066-4804/00/$04.00 0 Copyright 2000, American Society for Microbiology. All Rights Reserved. A Novel Human Immunodeficiency Virus Type 1 Reverse Transcriptase Mutational Pattern Confers Phenotypic Lamivudine Resistance in the Absence of Mutation 184V KURT HERTOGS, 1 * STUART BLOOR, 2 VERONIQUE DE VROEY, 1 CHRISTEL VAN DEN EYNDE, 1 PASCALE DEHERTOGH, 1 ANJA VAN CAUWENBERGE, 1 MARTIN STÜRMER, 3 TIMOTHY ALCORN, 4 SCOTT WEGNER, 5 MARGRIET VAN HOUTTE, 1 VERONICA MILLER, 3 AND BRENDAN A. LARDER 2 Virco NV, Mechelen, Belgium 1 ; Virco Ltd., Cambridge, United Kingdom 2 ; Klinikum der J.W. Goethe-Universität, Zentrum der Inneren Medizin, Frankfurt, Germany 3 ; LabCorp, Center for Molecular Biology and Pathology, Research Triangle Park, North Carolina 4 ; and U.S. Military HIV Research Program, Rockville, Maryland 5 Received 1 July 1999/Returned for modification 22 September 1999/Accepted 9 December 1999 We describe a new human immunodeficiency virus type 1 (HIV-1) mutational pattern associated with phenotypic resistance to lamivudine (3TC) in the absence of the characteristic replacement of methionine by valine at position 184 (M184V) of reverse transcriptase. Combined genotypic and phenotypic analyses of clinical isolates revealed the presence of moderate levels of phenotypic resistance (between 4- and 50-fold) to 3TC in a subset of isolates that did not harbor the M184V mutation. Mutational cluster analysis and comparison with the phenotypic data revealed a significant correlation between moderate phenotypic 3TC resistance and an increased incidence of replacement of glutamic acid by aspartic acid or alanine and of valine by isoleucine at residues 44 and 118 of reverse transcriptase, respectively. This occurred predominantly in those isolates harboring zidovudine resistance-associated mutations (41L, 215Y). The requirement of the combination of mutations 41L and 215Y with mutations 44D and 44A and/or 118I for phenotypic 3TC resistance was confirmed by site-directed mutagenesis experiments. These data support the assumption that HIV-1 may have access to several different genetic pathways to escape drug pressure or that the increase in the frequency of particular mutations may affect susceptibility to drugs that have never been part of a particular regimen. The emergence of drug-resistant human immunodeficiency virus type 1 (HIV-1) variants is almost always observed during the course of treatment of patients with antiretroviral drugs (3, 10, 14 16, 18, 21, 27; L. T. Bacheler, E. Anton, S. Jeffrey, H. George, G. Hollis, K. Abremski, and the Sustiva Resistance Study Team, Abstr. 2nd Int. Workshop HIV Drug Resistance and Treatment Strategies, abstr. 19, p. 15, 1998). The mutational profile of the resistant viruses generally is characteristic for the particular drug(s) taken. For example, mutations at codons 41, 67, 70, 210, 215, and 219 of reverse transcriptase (RT) typically confer resistance to zidovudine (ZDV) (6, 12, 13, 27). Similarly, mutation M184V in RT has been shown to be specifically associated with high-level ( 50-fold) phenotypic resistance to lamivudine (3TC) (1, 22, 28). No specific mutation(s) associated with moderate levels of phenotypic resistance (4- to 50-fold) to 3TC has been described before. Those mutations that confer moderate (4- to 50-fold) levels of phenotypic resistance to 3TC reported previously always appeared in the context of a constellation of mutations that confer resistance to multiple nucleoside analogues or as a cross-resistance phenomenon that appears with the emergence of resistance to another nucleoside analogue. This has been the case for the nucleoside multidrug resistance complex of mutations Q151M, F77L, F116Y, A62V, and V75I, although the increase in the level of phenotypic resistance to 3TC in viruses that harbor those mutations is slight (9, 20, 24, 25). In the case * Corresponding author. Mailing address: Virco NV, Generaal De Wittelaan 11 B4, B-2800 Mechelen, Belgium. Phone: 32 15 28 63 20. Fax: 32 15 28 63 46. E-mail: kurt.hertogs@virco.be. of the insertion mutations near position 69 of RT, a notable increase in the frequency of 3TC resistance has been reported together with an increased frequency of phenotypic resistance to other nucleosides (2, 17, 29). The K65R mutation appears infrequently during the course of treatment with dideoxyinosine (ddi; didanosine) and dideoxycytosine and confers concomitant low-level 3TC resistance (4, 5). In another study, phenotypic 3TC resistance (from 2.2-fold to 8.6-fold) was observed in the absence of M184V in patient isolates with various levels of phenotypic resistance to ZDV (from 4-fold to 378- fold) (26). However, no specific 3TC resistance-associated mutation was reported in this patient group. Mutations present simultaneously may act synergistically, causing drug resistance to increase as the number of resistance-associated mutations increases (19). On the other hand, a mutation may reverse the effect of another mutation occurring concurrently in a viral variant. This is the case for the 3TC resistance mutation M184V, which can reverse the effect conferred by the ZDV mutation T215Y (1, 14, 28). Thus, although resistance to a particular antiretroviral drug is conferred by a specific mutation(s) associated with a particular drug, this can also be modulated by other mutations present in the background. The implementation of high-throughput phenotyping and genotyping assays has allowed us to establish a database containing phenotypic resistance data for and the genotypic sequences of over 6,000 clinical isolates. Correlative data analysis and mutational cluster analysis then enabled us to search for mutational patterns with the accompanying phenotypic resistance profile. 568
VOL. 44, 2000 PHENOTYPIC RESISTANCE TO 3TC IN HIV-1 569 We describe here a novel mutational pattern in HIV-1 RT associated with a moderate level of phenotypic resistance (4- to 50-fold) to 3TC that is distinguishable from the characteristically high levels of phenotypic resistance ( 50-fold) associated with the M184V mutation and whose expression is modulated by a particular genetic background. MATERIALS AND METHODS Plasma samples. Plasma samples obtained from HIV-1-infected individuals from routine clinical practices in Europe and the United States were shipped to the laboratory (Virco NV, Virco Ltd.) on dry ice and were stored at 70 C until genotypic and phenotypic analyses. Since the research at this stage is concerned with the discovery and description of the mutational patterns that underlie resistance and with virological and molecular biological aspects of resistance, selection of plasma samples was not necessarily based upon the availability of the therapeutic histories for each subject who provided a sample. Phenotypic resistance testing. Phenotypic analysis was performed by the recombinant virus assay (11) approach described by Hertogs et al. (7) with the modifications described by Pauwels et al. (Abstr. 2nd Int. Workshop HIV Drug Resistance and Treatment Strategies, Abstr. 51, p. 35, 1998) (Antivirogram; Virco, Mechelen, Belgium). Briefly, protease (PR)- and RT-coding sequences were amplified from patient-derived viral RNA with HIV-1-specific primers. After homologous recombination of amplicons into a proviral clone from which the PR- and RT-coding sequences were deleted, the resulting recombinant viruses were harvested, titrated, and used for testing of in vitro susceptibility to antiretroviral drugs. The results of this analysis are expressed as fold change values, which reflect the fold increase in the mean 50% inhibitory concentration (IC 50 ) (micromolar) of a particular drug when it was tested with patient-derived recombinant virus isolates relative to the mean IC 50 (in micromolar) of the same drug obtained when it was tested with a reference wild-type virus isolate (strain IIIB/LAI). The 3TC resistance values generated by the phenotypic assay were grouped into one of three susceptibility categories, susceptible, moderately resistant, and high-level resistant, corresponding to fold changes in the IC 50 of between 0 and 4-fold, 4-fold and 50-fold, and 50-fold, respectively. Sequence analysis and genotypic-phenotypic correlative data analysis. Genotypic analysis was performed by automated population-based full-sequence analysis (ABI). The results of the genotypic analysis are reported as amino acid changes at positions along the RT gene compared to the RT gene sequence of wild-type reference strain (strain HXB2). The database containing the genetic sequences was searched to determine the occurrence and frequencies of occurrence of mutational patterns present in the sequences of the clinical isolates. The corresponding phenotypic resistance profiles of those isolates were retrieved from the phenotypic database for comparison. Site-directed mutagenesis. Mutations were generated in the RT gene of HXB2, a wild-type laboratory HIV-1 strain, with the QuickChange Site-Directed Mutagenesis Kit (Stratagene; Stratagene Cloning Systems, La Jolla, Calif.). RESULTS Analysis of the clinical isolates. The phenotypic database was interrogated for isolates susceptible or resistant to ZDV and/or 3TC. The corresponding genotypes were retrieved from the database. These isolates were drawn from a pool of 1,083 clinical samples for which both genotypic and phenotypic resistance data were available. In this group 42.2, 17.5, and 40.3% were found to have phenotypic susceptibility, moderate resistance, and high-level resistance to 3TC, respectively. Table 1 reports the frequencies of occurrence of the mutations 41L, 215Y, 184V, 44D and 44A, and 118I in RT in isolates with various levels of phenotypic resistance to ZDV and 3TC. Throughout our analysis we found that the ZDV and 3TC phenotypic resistance profile, as well as the co-occurrence of ZDV and 3TC phenotypic resistance with any of the ZDV resistance-associated genetic changes (see below), did not differ between isolates with 44D and 44A mutations (data not shown). Therefore, the data for both mutations (referred to as 44D/A) are pooled in Tables 1 to 3 and Fig. 1. Isolates that are phenotypically susceptible (wild type) to both ZDV and 3TC (n 314). As expected, the frequency of occurrence of any of the mutations listed above was low for the subset of isolates that are phenotypically susceptible to both ZDV and 3TC. TABLE 1. Frequency of ZDV and 3TC resistance-associated mutations in clinical isolates susceptible or resistant to ZDV and/or 3TC compared to a sample of fully susceptible isolates Resistance class a Frequency (%) of mutations ZDV resistanceassociated mutations 3TC resistanceassociated mutations 41L 215Y 184V 44D/A 118I No. of samples ZDV ( 4), 3TC ( 4) 4.5 4.8 0 1.3 3.1 314 ZDV ( 4), 3TC ( 10) 18.3 18.8 90 1.3 6.3 240 ZDV ( 10), 3TC ( 4) 59.5 68.9 0 14.9 18.9 74 ZDV ( 10), 3TC ( 4, 50) 77 72.2 4 30.2 39.7 126 ZDV ( 10), 3TC ( 50) 77.5 66.9 84.1 28.5 37.8 151 a Resistance (in parentheses) is expressed as the fold increase in the mean IC 50 of the drug relative to the mean IC 50 of the same drug for a wild-type reference laboratory HIV-1 strain. Isolates that are phenotypically susceptible to ZDV but for which the 3TC IC 50 is greater than 10-fold (n 240). Table 1 shows that the frequency of occurrence of 3TC resistanceassociated mutation 184V is high among isolates that are phenotypically susceptible to ZDV but for which the 3TC IC 50 is greater than 10-fold. The frequencies of occurrence of mutations 44D/A and 118I are low. About 18% of the isolates in this group harbor the major ZDV resistance-associated mutations 41L and 215Y, yet the level of ZDV resistance is low among these isolates, presumably due to the reversal of phenotypic ZDV resistance by the 184V mutation. Isolates that are phenotypically resistant to ZDV (>10-fold; n 351). The group of isolates that are phenotypically resistant to ZDV ( 10-fold) was divided into the three 3TC resistance categories as detailed above. The mutation frequencies indicate that the rate of occurrence of ZDV resistance-associated mutations 41L and 215Y is high among all three 3TC resistance groups. As expected, the 184V mutation was almost exclusively present in the subset of isolates with high-level phenotypic 3TC resistance ( 50-fold), whereas the same mutation was low in frequency (4%) in the group of ZDV-resistant isolates with a moderate level of phenotypic 3TC resistance. This mutation was absent from the subset of isolates susceptible to 3TC. Mutations 44D/A and 118I were present in isolates in all 3TC resistance categories, and their incidence was higher among isolates with moderate and high-level phenotypic 3TC resistance than among isolates that were phenotypically susceptible to 3TC. These results showed that moderate phenotypic 3TC resistance (4- to 50-fold) was not correlated with the presence of 184V. Indeed, this mutation was present in low numbers in isolates in these categories. Table 1 further indicates that mutations 44D/A and 118I were present at high frequencies only in the presence of the ZDV resistance mutations 41L and 215Y. In the isolates that were phenotypically susceptible to ZDV, the frequency of occurrence of ZDV resistance-associated mutations was low, as were the frequencies of occurrence of the 44D/A and 118I mutations (even though 3TC resistance was greater than 10-fold). In this group, the high frequency of the 184V mutation accounted for the resistance to 3TC. Figure 1 graphically presents the frequencies of occurrence and the distributions of mutations 184V, 44D/A, and 118I, in addition to the 65R mutation, the two-amino-acid insertion at position 69, and the 151M mutation, among the groups susceptible to 3TC and with moderate and high-level resistance to 3TC. The frequencies of occurrence of the 65R mutation, the
570 HERTOGS ET AL. ANTIMICROB. AGENTS CHEMOTHER. FIG. 1. Frequency and distribution of mutation 65R, the two-amino-acid insertion at position 69, and mutations 151M, 184V, 44D/A, and 118I among the groups susceptible to 3TC and with intermediate and high-level phenotypic resistance to 3TC (n 1,083). two-amino-acid insertion at position 69, and the 151M mutation were very low among isolates in all three 3TC resistance categories, whereas the frequency of occurrence of the 184V mutation was very high among isolates in the high-level 3TC resistance category. The 44D/A and 118I changes were present in substantial numbers of isolates, especially those in the moderate and high-level 3TC resistance categories. Next, we queried the genotypic database for the co-occurrence of mutations 44D/A and 118I, any of the ZDV resistance-associated mutations 41L, 67N, 70R, 210W, 215Y/F, and 219Q/E, and the 3TC resistance-associated mutation 184V. Table 2 shows the frequency with which mutations 44D/A and 118I occurred together with any of the mutations 41L, 67N, 70R, 210W, 215Y/F, and 219Q/E, as well as with 184V. These data extracted from the database confirmed the results obtained earlier and presented above, which indicated that a ZDV resistance background is required for the mutations 44A/D and 118I to accumulate in substantial numbers. These data also confirm that this requirement is less stringent for 118I than for 44D/A. Statistical analysis showed that the co-occurrence of 44D/A and 118I with ZDV resistance-associated mutations was significant (P 0.0001). However, no such association was seen between the occurrence of 44D/A and 118I and the occurrence of 184V (P 0.15) in these clinical isolates. Mutations 118I, 44D, and 44A occurred with overall frequencies of 15.07, 9.23, and 1.46%, respectively, in our database. Forty-eight percent of the clinical isolates with mutation 118I simultaneously harbored the 44D/A mutations. In contrast, 65 and 81% of the clinical isolates carrying the 44D and 44A mutations, respectively, simultaneously harbored the 118I mutation. Analysis of mutant viruses generated by site-directed mutagenesis. In order to investigate the relationship between phenotypic ZDV and/or 3TC resistance and the presence of mutations at positions 41, 215, 184, 44, and 118, different combinations of genotypic profiles were generated by site-directed mutagenesis. Table 3 lists a series of mutants that carried codon changes introduced into the wild-type HXB2 background together with the corresponding fold change in IC 50 obtained by a drug susceptibility assay. Three mutants with a change at codon 44, three mutants with a change at codon 118, and three mutants with changes at both codons 44 and 118 were generated. Within each of these three groups two mutants also had changes at positions associated with resistance to ZDV, whereas one mutant had the wild-type sequence at those codons. The drug resistance data for those mutants showed that the presence of mutations at codons 44 and 118, singly or together, could cause moderate phenotypic resistance to 3TC (7- to 32-fold). However, the moderate resistance to 3TC was observed only when mutations at positions 44 and/or 118 were in a ZDV resistance background (41L, 67N, 210W, 215Y). Six mutants with a change at codon 184 were generated, with one of those retaining the wild-type sequences at the other positions, whereas three others carried a series of ZDV resistanceassociated mutations and two had changes at codons 44 and 118 as well as a series of ZDV resistance-associated mutations. All six of those mutants had high-level resistance to 3TC ( 50- fold) that was distinguishable from the moderate level of phenotypic resistance seen in the mutants with changes at codons 44 and 118. In all of the mutants with the 184V mutation, resistance to 3TC was not related to the presence of the ZDV resistance-associated mutations. All mutants showed the predicted ZDV resistance or susceptibility pattern. The presence of the 184V mutation had a ZDV resensitizing effect on the mutants that carried ZDV resistance-associated mutations. This was not the case for the mutants with changes at codons 44 and 118, as phenotypic ZDV resistance remained unchanged when the mutations were introduced into the ZDV resistance background (Table 3). Relationship between presence of changes at RT positions 44 or 118 in clinical samples and antiretroviral therapy. As can be deduced from Tables 1 and 2, changes at positions 44 and 118 can occur in viruses with or without the M184V substitution, but their incidence appeared to be the highest in viruses with ZDV resistance. It was therefore of interest to examine the antiretroviral treatments of patients infected with HIV isolates containing 44D or 118I. We identified a subset of 86 samples with isolates with the 44D mutation and 88 samples with isolates with the 118I mutation that originated from patients for whom antiretroviral treatment histories were available. It was not possible to draw conclusions regarding the incidence of changes at position 44 or 118 from this subset according to treatment history, as this was not a randomized study. However, this analysis sheds light on the requirements that lead to mutations at these positions. For the subset of samples with isolates with the 44D mutation, 50 of 86 of the samples originated from patients who were receiving 3TC at the sampling data; 5 samples in this subset were from patients who had never received 3TC prior to the sampling date. All five patients had received ZDV-ddI at some TABLE 2. Incidence of co-occurrence of mutations 44D/A and 118 (either singly or together in an isolate) with ZDV resistanceassociated mutations 41L, 67N, 70R, 210W, 215Y/F, and 219Q/E and with 3TC resistance-associated mutation 184V in RT in clinical HIV-1 isolates a Status of ZDV resistance-associated positions Status of 3TC resistance-associated position 184 44D/A No. (%) of isolates with the following mutation(s): 118I 44D/A 118I Wild type Wild type 3 (1.54) 33 (7.55) 1 (0.25) Mutant Wild type 104 (53.33) 197 (45.08) 188 (47.36) Mutant Mutant 87 (44.62) 191 (43.71) 175 (44.08) Wild type Mutant 1 (0.51) 16 (3.66) 33 (8.31) Total 195 (100) 437 (100) 397 (100) a A total of 5,536 genetic sequences were queried. The frequencies of mutations 44D, 44A, and 118I in the genetic database were 9.23, 1.46, and 15.07%, respectively.
VOL. 44, 2000 PHENOTYPIC RESISTANCE TO 3TC IN HIV-1 571 TABLE 3. 3TC and ZDV resistance-associated mutations and phenotypic resistance in mutants with site-directed mutations Mutant 3TC resistance profile AZT resistance profile Mutation(s) Fold resistance a n b Mutations Fold resistance n SDM05 4 (0.3) c 4 41L, 210W,211K,214F,215Y 64 (15.0) 5 SDM18 2 (0.6) 3 41L, 67N,210W,211K,214F,215Y 45 (13.4) 3 SDM19 4 (0.7) 3 41L, 67N,69D,210W,211K,214F,215Y 46 (18.2) 2 SDM28 44D 1 (0.1) 6 2 (0.3) 6 SDM31 44D 22 (2.5) 5 41L, 67N,210W,211K,214F,215Y 48 (11.8) 4 SDM32 44D 8 (2.2) 4 41L, 67N,69D,210W,211K,214F,215Y 49 (5.9) 6 SDM29 118I 2 (0.2) 6 2 (0.4) 6 SDM33 118I 7 (1.0) 5 41L, 67N,210W,211K,214F,215Y 49 (8.0) 6 SDM34 118I 32 (3.9) 5 41L, 67N,69D,210W,211K,214F,215Y 34 (14.4) 5 SDM30 44D,118I 3 (0.3) 6 1 (0.4) 6 SDM35 44D,118I 14 (1.4) 6 41L, 67N,210W,211K,214F,215Y 49 (9.6) 5 SDM36 44D,118I 15 (2.1) 6 41L, 67N,69D,210W,211K,214F,215Y 49 (10.6) 5 SDM22 184V 78 (16.3) 5 41L,210W,211K,214F,215Y 7 (0.9) 5 SDM23 184V 82 (13.6) 6 2 (0.5) 6 SDM24 184V 85 (14.2) 6 69S-S-S, 210W,211K,214F,215Y 27 (16.2) 5 SDM26 184V 72 (13.8) 5 41L, 67N,210W,211K,214F,215Y 25 (1.4) 5 SDM38 184V,44D,118I 82 (13.5) 6 41L, 67N,210W,211K,214F,215Y 20 (4.1) 6 SDM39 184V,44D,118I 84 (13.9) 6 41L, 67N,69D,210W,211K,214F,215Y 21 (5.3) 5 a Fold increase in the mean IC 50 of the drug relative to the mean IC 50 of the same drug for a wild-type reference laboratory HIV-1 strain. b n, number of replicate tests run for each phenotypic drug resistance determination. c Standard errors are indicated in parentheses. time, and all virus isolates from these five patients had the wild-type sequence at codon 184. The ZDV treatment experiences of these patients were extensive. All except one of the patients had received ZDV in combination with other nucleoside RT inhibitors, and 70 of 86 had also received ZDV monotherapy at some time previously. The one patient reported to be ZDV naive had received stavudine. The virus in the sample from this patient nevertheless contained mutations 41L and 215Y. Results for the subset of isolates with the 118I mutation were similar: 55 of 88 samples with isolates with this mutation originated from patients who were receiving 3TC at the sampling date. Two patients had never received 3TC (both had received ZDV plus ddi). Eighty-three of 88 of the patients had received ZDV in combination with other nucleoside RT inhibitors, and 70 had also received ZDV monotherapy. The five ZDV-naive patients had received stavudine. DISCUSSION The results presented here indicate that the E44D/A and V118I substitutions in HIV-1 RT confer a moderate level (4- to 50-fold) of phenotypic resistance to 3TC when they occur together with a ZDV resistance background. The presence of the M184V mutation is not required. The results obtained from the site-directed mutagenesis experiments confirmed the phenotypic data for the clinical isolates. This resistance-conferring mutational pattern has not previously been described and was discovered through mutational cluster analysis and correlative data analysis of a database containing phenotypic drug resistance data and the corresponding genetic sequences of these clinical isolates. The study by Skowron et al. (26) described low-level resistance to 3TC in the absence of the M184V mutation in patient isolates that were also resistant to ZDV. Those investigators noted that the 3TC resistance levels were highly correlated with the ZDV resistance levels. Mutations at positions 41, 67, 70, 215, and 219 were present in the majority of the patient isolates. The investigators reported that patients whose isolates had mutations at codons 70 and 215 experienced a lesser response to treatment with 3TC than patients who did not harbor virus with those mutations, implying that mutations at codons 70 and 215 might interfere with the action of 3TC. In our study, the frequency of the mutation at position 215 (as well as the mutation at position 41) was high in all isolates with intermediate and high-level ZDV resistance, regardless of the level of 3TC resistance. This leads us to conclude that other mutations are responsible for the low and intermediate levels of 3TC resistance that we observed among the isolates in these samples. The cluster analysis of the approximately 1,000 genotypically and phenotypically characterized clinical isolates and the results from the site-directed mutagenesis experiment confirm that mutations at codons 44 and 118 are indeed associated with moderate levels of phenotypic resistance to 3TC when they are present with ZDV resistance-associated mutations. The analysis of the clinical samples from patients for whom therapeutic histories were available and for which prior ZDV exposure was shown to be extensive confirmed the results obtained from our large clinical data set in that mutations 44D/A and 118I appeared in the context of ZDV mutations. The appearance of a mutation at position 118 in a ZDV resistance background has been documented previously (23). In a study of 39 patients who received ZDV and ddi combination therapy for 2 years, isolates from 5 patients were reported to have a mutation at position 118 at the end of the 2-year period. All five patients had received ZDV prior to the start of the combination therapy. The isolates from the five patients all had mutations at positions 41 and 215. Interestingly, isolates from one of the patients in this group of five patients harbored the mutation at position 118 at the start of combination therapy and also had a mutation at position 215 at the start of combination therapy. Mutations 44D/A and 118I in the context of ZDV mutations are both capable of independently generating resistance to 3TC. The site-directed mutagenesis experiments did not indicate the existence of synergistic effects between the two mutations with respect to their phenotypic effect on 3TC resistance. The 3TC resistance mutation 184V appears independently
572 HERTOGS ET AL. ANTIMICROB. AGENTS CHEMOTHER. of a ZDV resistance background and clearly appears to be due to 3TC selection pressure, whereas the accumulation of 44D/A and 118I appears to be driven by ZDV and, fortuitously, confers resistance to 3TC. The mutational patterns of the mutants with site-directed mutations and the clinical isolates showed that the 44D/A and 118I mutations can occur with or without the 184V mutation. This indicates that 44D/A and 118I are not necessarily alternative mutations to 184V but suggests that 184V is a more dominant mutation with respect to 3TC resistance. We are not able to differentiate the effects of 44D/A and 118I from the effect of 184V in the isolates in which they occur together. Thus, we do not know whether 44D/A and 118I contribute to phenotypic resistance when 184V is simultaneously present. The data presented here thus point to the existence of one possible genetic pathway for intermediate phenotypic 3TC resistance in clinical isolates of HIV-1. The data in Table 1 show that not all recombinant clinical isolates with an intermediate level of phenotypic resistance to 3TC and with high-level resistance to ZDV possess mutations at positions 44 and 118. This implies that other polymorphisms with an effect on phenotypic 3TC resistance may exist. These polymorphisms will need to be identified through further study. The clinical relevance of the accumulation of 44D/A and 118I in a ZDV resistance background is under investigation by means of a longitudinal follow-up study. Further studies are also required to establish whether the 3TC resistance-associated changes that occur in the absence of 3TC treatment may provide an explanation for the possible delayed appearance of 184V during subsequent 3TC treatment (data not shown). It is of interest to examine the positions 44 and 118 in the three-dimensional structures of the RT enzyme. The recently described structure of the catalytic complex of RT with a deoxynucleoside triphosphate template and primer sheds new insight on the possible way that nucleoside analog mutations confer resistance (8). Residue 215 is located near the incoming nucleotide-binding site. Of interest, residue 118 is located close to residue 215 in the three-dimensional structure of the catalytic complex. Furthermore, residue 44 is located in close proximity to residue 41. It is conceivable that changes at positions 215 and 41 of RT facilitate the accumulation of changes at positions 118 and 44 because of three-dimensional structural constraints. However, the mutations at positions 44 and 118 do not affect phenotypic ZDV resistance. Whether the changes at positions 118 and 44 influence the replicative fitness of ZDVresistant viruses remains to be determined. Our analysis shows that the pathway toward resistance to a particular drug may critically depend on the sequence in which drugs are administered in the course of therapy. As new drugs are approved for treatment, physicians will not only need to continue to evaluate the effects of the drugs in the context of the concurrently administered drugs but will also need to consider previously taken drugs. This not only is because of the possibility of cross-resistance but is also because of phenomena such as the one discussed here. This is also true in situations in which recycling of drugs is being considered. The results presented in this report indicate that the significant differences in phenotype which will probably be reflected in the response to treatment may not be readily predicted by the genotype. This points to the importance of combining the knowledge about both the phenotype and the genotype when therapeutic decisions are made. ACKNOWLEDGMENTS We thank the staff at Virco in Belgium, Virco Ltd. in the United Kingdom, the J.W. Goethe-Universität in Frankfurt, Germany, the U.S. Military HIV Research Program in Rockville, Md., and LabCorp at Research Triangle Park, N.C., for assistance with the phenotypic and genotypic analyses. REFERENCES 1. Boucher, C. A. B., N. Cammack, P. Schipper, R. Schuurman, P. Rouse, M. A. Wainberg, and J. M. Cameron. 1993. High-level resistance to ( ) enantiomeric 2 -deoxy-3 -thiacytidine in vitro is due to one amino acid substitution in the catalytic site of human immunodeficiency virus type 1 reverse transcriptase. Antimicrob. Agents Chemother. 37:2231 2234. 2. de Jong, J. J., J. Goudsmit, V. V. Lukashov, M. E. Hillebrand, E. Baan, R. Huismans, S. A. Danner, J. H. ten Veen, F. de Wolf, and S. Jurriaans. 1999. Insertion of two amino acids combined with changes in reverse transcriptase containing tyrosine-215 of HIV-1 resistant to multiple nucleoside analogs. AIDS 13:75 80. 3. Drusano, G. L., J. A. Bilello, D. S. Stein, M. Nessly, A. Meibohm, E. A. Emini, P. Deutsch, J. Condra, J. Chodakewitz, and D. J. Holder. 1998. Factors influencing the emergence of resistance to indinavir: role of virologic, immunologic, and pharmacologic variables. J. Infect. Dis. 178:360 367. 4. Gu, Z., Q. Gao, H. Fang, H. Salomon, M. A. Partiak, E. Goldberg, J. Cameron, and M. A. Wainberg. 1994. Identification of a mutation at codon 65 in the IKKK motif of reverse transcriptase that encodes human immunodeficiency virus resistance to 2,3 -dideoxycytidine and 2,3 -dideoxy-3 thiacytidine. Antimicrob. Agents Chemother. 38:275 281. 5. Gu, Z., R. S. Fletcher, E. J. Arts, M. A. Wainberg, and M. A. Parniak. 1994. The K65R mutant reverse transcriptase of HIV-1 cross-resistant to 2,3 -dideoxycytidine, 2,3 -dideoxy-3 -thiacytidine, and 2,3 -dideoxyinosine shows reduced sensitivity to specific dideoxynucleoside triphosphate inhibitors in vitro. J. Biol. Chem. 269:28118 28122. 6. Harrigan, P. R., I. Kinghorn, S. Bloor, S. D. Kemp, I. Nájera, A. Kohli, and B. A. Larder. 1996. Significance of amino acid variation at human immunodeficiency virus type 1 reverse transcriptase residue 210 for zidovudine susceptibility. J. Virol. 70:5930 5934. 7. Hertogs, K., M.-P. de Béthune, V. Miller, T. Ivens, P. Schel, A. Van Cauwenberge, C. Van Den Eynde, V. Van Gerwen, H. Azijn, M. Van Houtte, F. Peeters, S. Staszewski, M. Conant, S. Bloor, S. Kemp, B. Larder, and R. Pauwels. 1998. A rapid method for simultaneous detection of phenotypic resistance to inhibitors of protease and reverse transcriptase in recombinant human immunodeficiency virus type 1 isolates from patients treated with antiretroviral drugs. Antimicrob. Agents Chemother. 42:269 276. 8. Huang, H., R. Chopra, G. L. Verdine, and S. C. Harrison. 1998. Structure of a covalently trapped catalytic complex of HIV-1 reverse transcriptase: implications for drug resistance. Science 282:1669 1675. 9. Iversen, A. K., R. W. Shafer, K. Wehrly, M. A. Winters, J. Mullins, B. Cheseboro, and T. C. Merigan. 1996. Multidrug-resistant human immunodeficiency virus type 1 strains resulting from combination antiretroviral therapy. J. Virol. 70:1086 1090. 10. Jacobsen, H., M. Hänggi, M. Ott, I. B. Duncan, S. Owen, M. Andreoni, S. Vella, and J. Mous. 1996. In vivo resistance to a human immunodeficiency virus type 1 proteinase inhibitor: mutations, kinetics, and frequencies. J. Infect. Dis. 173:1379 1387. 11. Kellam, P., and B. A. Larder. 1994. Recombinant virus assay: a rapid, phenotypic assay for assessment of drug susceptibility of human immunodeficiency virus type 1 isolates. Antimicrob. Agents Chemother. 38:23 30. 12. Kellam, P., C. A. B. Boucher, and B. A. Larder. 1992. Fifth mutation in human immunodeficiency virus type 1 reverse transcriptase contributes to the development of high-level resistance to zidovudine. Proc. Natl. Acad. Sci. USA 89:1934 1938. 13. Larder, B. A., and S. D. Kemp. 1989. Multiple mutations in HIV-1 reverse transcriptase confer high-level resistance to zidovudine (AZT). Science 246: 1155 1158. 14. Larder, B. A., S. D. Kemp, and P. R. Harrigan. 1995. Potential mechanism for sustained antiretroviral efficacy of AZT-3TC combination therapy. Science 269:696 699. 15. Miller, V., M.-P. de Béthune, A. Kober, M. Stürmer, K. Hertogs, R. Pauwels, P. Stoffels, and S. Staszewski. 1998. Patterns of resistance and cross-resistance to human immunodeficiency virus type 1 reverse transcriptase inhibitors in patients treated with the nonnucleoside reverse transcriptase inhibitor loviride. Antimicrob. Agents Chemother. 42:3123 3129. 16. Montaner, J. S. G., P. Reiss, D. Cooper, S. Vella, M. Harris, B. Conway, M. A. Wainberg, D. Smith, P. Robinson, D. Hall, M. Myers, and J. M. A. Lange for the INCAS Study Group. 1998. A randomized, double-blind trial comparing combinations of nevirapine, didanosine and zidovudine for HIVinfected patients. The INCAS trial. JAMA 279:930 937. 17. Rakik, A., M. Ait-Khaled, P. Griffin, D. Thomas, J.-P. Kleim, and M. Tisdale, for the Abacavir CNA2007 International Study Group. 1998. A novel genotype encoding a single amino acid insertion and five other substitutions between residues 64 and 74 of the HIV-1 reverse transcriptase confers cross-resistance to NRTIs. AIDS 12(Suppl. 4):S25. 18. Richman, D. D., D. Havlir, J. Corbeil, D. Looney, C. Ignacio, S. A. Spector, J. Sullivan, S. Cheeseman, K. Barringer, D. Pauletti, C.-K. Shih, M. Myers, and J. Griffin. 1994. Nevirapine resistance mutations of human immunode-
VOL. 44, 2000 PHENOTYPIC RESISTANCE TO 3TC IN HIV-1 573 ficiency virus type 1 selected during therapy. J. Virol. 68:1660 1666. 19. Richman, D. D., J. M. Grimes, and S. W. Lagakos. 1990. Effect of stage of disease and drug dose on zidovudine susceptibilities of isolates of human immunodeficiency virus. J. Acquir. Immune Defic. Syndr. 3:743 746. 20. Schmit, J. C., K. Van Laethem, L. Ruiz, P. Hermans, S. Sprecher, A. Sonnerborg, M. Leal, T. Harrer, B. Clotet, V. Arendt, E. Lissen, M. Witvrouw, J. Desmyter, E. De Clercq, and A.-M. Vandamme. 1998. Multiple dideoxynucleoside analogue-resistant isolated from patients from different European countries. AIDS 12:2007 2015. 21. Schmit, J.-C., L. Ruiz, B. Clotet, A. Raventos, J. Tor, J. Leonard, J. Desmyter, E. De Clercq, and A.-M. Vandamme. 1996. Resistance-related mutations in the HIV-1 protease gene of patients treated for 1 year with the protease inhibitor ritonavir (ABT-538). AIDS 10:995 999. 22. Schuurman, R., M. Nijhuis, R. van Leeuwen, P. Schipper, D. de Jong, P. Collis, S. A. Danner, J. Mulder, C. Loveday, C. Christopherson, S. Kwok, J. Sninsky, and C. A. B. Boucher. 1995. Rapid changes in human immunodeficiency virus type 1 RNA load and appearance of drug-resistant populations in persons treated with lamivudine (3TC). J. Infect. Dis. 171:1411 1419. 23. Shafer, R. W., A. K. Iversen, M. A. Winters, E. Aguiniga, D. A. Katzenstein, and T. C. Merigan. 1995. Drug resistance and heterogeneous long-term virologic responses of human immunodeficiency virus type 1-infected subjects to zidovudine and didanosine combination therapy. The AIDS Clinical Trials Group 143 Virology Team. J. Infect. Dis. 172:70 78. 24. Shirasaka, T., M. F. Kavlick, T. Ueno, W.-Y. Gao, E. Kojima, M. L. Alcaide, S. Chokekijchai, B. M. Roy, E. Arnold, R. Yarchoan, and H. Mitsuya. 1995. Emergence of human immunodeficiency virus type 1 variants with resistance to multiple dideoxynucleosides in patients receiving therapy with dideoxynucleosides. Proc. Natl. Acad. Sci. USA 92:2398 2402. 25. Shirasaka, T., R. Yarchoan, M. O Brien, R. Husson, B. Anderson, E. Kojima, S. Broder, and H. Mitsuya. 1993. Changes in drug sensitivity of human immunodeficiency virus type 1 during therapy with azidothymidine, dideoxycytidine, and dideoxyinosine: an in vitro comparative study. Proc. Natl. Acad. Sci. USA 90:562 566. 26. Skowron, G., C. Petropoulos, M. Holodniy, M. Wesley, L. Ferrigno, and K. Frost for the American Foundation for AIDS Research Community Based Clinical Trials Network, New York, NY, USA. 1998. Reduced nucleoside analogue susceptibility patterns and correlation with proportion of mutant virus detected by differential hybridization in patients receiving AZT and ddi. AIDS 12(Suppl. 4):S19. 27. St. Clair, M. H., J. L. Martin, G. Tudor-Williams, M. C. Bach, C. L. Vavro, D. M. King, P. Kellam, S. D. Kemp, and B. A. Larder. 1991. Resistance to ddi and sensitivity to AZT induced by a mutation in HIV-1 reverse transcriptase. Science 253:1557 1559. 28. Tisdale, M., S. D. Kemp, N. R. Parry, and B. A. Larder. 1993. Rapid in vitro selection of human immunodeficiency type 1 resistant to 3 -thiacytidine inhibitors due to a mutation in the YMDD region of reverse transcriptase. Proc. Natl. Acad. Sci. USA 90:5653 5656. 29. Winters, M. A., K. L. Coolley, Y. A. Girard, D. J. Levee, H. Hamdan, R. W. Shafer, D. A. Katzenstein, and T. C. Merigan. 1998. A 6-basepair insert in the reverse transcriptase gene of human immunodeficiency virus type 1 confers resistance to multiple nucleoside inhibitors. J. Clin. Invest. 102:1769 1775. Downloaded from http://aac.asm.org/ on June 15, 2018 by guest