Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, November 2003, p. 3478-3484, Vol. 47, No. 11
0066-4804/03/$08.00+0 DOI: 10.1128/AAC.47.11.3478-3484.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Institute of Clinical and Molecular Virology, German National Reference Centre for Retroviruses, University of Erlangen-Nürnberg, Erlangen,1 Max Planck Institute for Informatics, Saarbrücken,2 Institute for Virology, University of Cologne, Cologne,4 Center of Advanced European Studies and Research, Bonn,5 Max Planck Institute of Molecular Plant Physiology, Golm, Germany,6 Virology Networks and University Medical Center Utrecht, Utrecht, The Netherlands,3 Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium7
Received 26 February 2003/ Returned for modification 12 May 2003/ Accepted 22 July 2003
|
|
|---|
|
|
|---|
In vitro experiments (20, 25) and the administration of tenofovir to SIV-infected rhesus macaques (24) have shown the emergence of the K65R mutation, which confers three- to fourfold reduced susceptibility to tenofovir. Furthermore, tenofovir resistance has been observed for multiple dideoxynucleoside-resistant viruses with insertions between RT amino acids 67 and 70, while tenofovir was still active in the presence of Q151M (10). K65R as well as insertions have only rarely been observed in clinical isolates (5). Likewise, resistance to tenofovir was detected at only a low frequency in isolates from treatment-experienced patients, with the greatest degree of cross-resistance observed in the context of reduced susceptibility to thymidine analogues (5). Furthermore, the M184V mutation was described to confer increased susceptibility to tenofovir (25) and to reduce the level of K65R-induced resistance (28). In the presence of thymidine analogue mutations (TAMs), M184V was shown to increase the susceptibility of viruses to tenofovir, at least in part (5); however, the clinical relevance of this effect is still unclear (8). Recently, 4.0-fold reduced susceptibility to tenofovir has been defined as the cutoff for the failure of therapy (8). More recent data suggest that the efficacy of tenofovir can be reduced at even lower levels of resistance (7).
This study was performed to focus on three aspects of tenofovir resistance: (i) to characterize the resistance profile of multiple dideoxynucleoside-resistant viral strains, (ii) to analyze TAM-induced cross-resistance and the degree of resensitization by M184V, and (iii) to develop computational models to predict tenofovir resistance from the genotypes by using the recently suggested clinically relevant cutoffs for a bioinformatics approach which has previously been described for other antiretroviral drugs (2, 3).
(This study was presented in part at the 11th International HIV Drug Resistance Workshop Basic Principles and Clinical Implications, Seville, Spain, July 2002 [K. Wolf, H. Walter, T. Schnell, W. Keulen, N. Beerenwinkel, J. Selbig, A.-M. Vandamme, K. Korn, and B. Schmidt, Abstr. 11th Int. HIV Drug Resist. Workshop Basic Principles Clin. Implications, abstr. 20, 2002].)
|
|
|---|
Resistance testing. Phenotypic resistance testing was performed by a recombinant virus assay as described previously (26). In brief, the complete protease gene and the first 900 bp of the RT gene were amplified from patient plasma and cloned into a matched deletion mutant of pNL4-3. Transient transfection of 293T cells resulted in recombinant viral stocks, which were analyzed on a CEMx174-derived cell line containing the gene for the secreted alkaline phosphatase (SEAP) under the control of the SIV long terminal repeat (9). Tenofovir resistance was determined in triplicate by measuring secreted alkaline phosphatase activity after 3 days by using six tenofovir concentrations ranging from 0.45 to 150 µM. The fold change in susceptibility was calculated by dividing the 50% inhibitory concentration (IC50) for the respective recombinant virus by the IC50 for nonresistant reference strain NL4-3, which was included in each independent assay. The fold reductions in the susceptibilities of viruses resistant to multiple dideoxynucleosides as well as the viral clones were determined in at least three independent runs. Genotypic analysis was performed by population sequencing of the amplification product described above by using dye terminators (Amersham, Cleveland, Ohio). The sequences were aligned by using the Wisconsin package (version 10.0; Genetics Computer Group, Madison, Wis.). The detection limit for minority species was about 30%. Resistance-associated mutations were analyzed as described by Schinazi et al. (17).
Mutagenesis. Viral clones were constructed by site-directed mutagenesis (6). Viral clones were constructed by insertion of the following mutations into the HXB2 backbone: M41L; K70R; M184I; M184V; T215F; T215Y; M41L and K70R; M41L and M184V; M41L and T215F; M41L and T215Y; K70R and M184V; K70R and T215Y; M184V and T215Y; M41L, K70R, and T215Y; M41L, M184V, and T215Y; K70R, M184V, and T215Y; and M41L, K70R, M184V, and T215Y. Construction of viral clones with mutations within the NL4-3 backbone (mutations M184V; M184V and T215Y; L210W and T215Y; M184V, L210W, and T215Y; M41L, L210W, R211K, and T215Y; and M41L, M184V, L210W, R211K, and T215Y) was described previously (27). For the mutagenesis of positions 62 and 65 of the RT, the following sense (s) and antisense (as) primers were used: A65A-K65K-s (5'-CCAGTATTTGCCATAAAGAAAAAAAATA-3'), A62A-K65K-as (5'-ATTTTTTTTCTTTATGGCAAATACTGGA-3'), A62V-K65K-s (5'-CCAGTATTTGTAATAAAGAAAAAAAA-3'), A62V-K65K-as (5'-TTTTTTTCTTTATTACAAATACTGGA-3'), A62A-K65R-s (5'-CCAGTATTTGCCATAAAGAGAAAAAA-3'), and A62A-K65R-as (5'-TTTTTCTCTTTATGGCAAATACTGGA-3').
Bioinformatics analysis. The DNA sequences were aligned with the pol gene of reference strain HXB2CG (GenBank accession number K03455) and translated into amino acids. One attribute was defined for each of the first 250 amino acids of the RT, with the attributes classified as an amino acid, a deletion, or unknown for ambiguous sequence information or a lack of sequence information. Additionally, one binary attribute was defined for the occurrence of an insertion. The viruses in the sample set were divided into susceptible and resistant subgroups according to the two phenotypic cutoffs which have recently been proposed as clinically relevant for a reduced response to therapy and treatment failure (1.5- and 4.0-fold, respectively) (7, 8). Mutual information profiles were created. These profiles quantify the information content (in bits) of each sequence position for the discrimination between tenofovir-susceptible and -resistant strains. Decision trees were generated by recursively splitting the data set according to the amino acid position with the largest amount of normalized mutual information plus the ratio between the amount of mutual information and its average amount over all RT sequence positions and by repeating this splitting for the respective subsets. Pruning was used to avoid overfitting.
Statistical analysis. For statistical evaluation of the data, Fisher's exact test and the Wilcoxon rank-sum test were applied as appropriate. Linear regression analysis was used to investigate the relationship between the number of TAMs and the median fold change in susceptibility. Significant values in the mutual information profiles were calculated by using 5,000 permutations for each position, which allowed the detection of significant differences from background at a 95% confidence level with Bonferroni's correction for multiple comparisons (see the legend to Fig. 2). The expected prediction error of the decision trees for unseen sequences was estimated by the leave-one-out method, a common cross-validation technique (3).
![]() View larger version (24K): [in a new window] |
FIG. 2. Mutual information profiles and decision tree models for tenofovir resistance obtained by using the two cutoffs for which clinical significance has been suggested (7, 8). (A) Mutual information profiles for 1.5- and 4.0-fold reduced susceptibility (upper and lower panels, respectively). The leftmost bar (position 0) on the x axis corresponds to the attribute indicating an insertion; positions 1 through 250 represent the first 250 amino acid positions of the HIV type 1 RT. The y axis gives the mutual information (in bits), which represents the information content of each sequence position that can be used to discriminate between tenofovir-susceptible and -resistant viruses. All amino acid positions with an information content significantly different from the background level are labeled. (B) Decision tree for 1.5-fold reduced susceptibility, which requires two branches with at least four viruses, resulting in a leave-one-out error of 19.1%. (C) Decision tree for 4.0-fold reduced susceptibility, which requires two branches with at least seven viruses, resulting in a leave-one-out error of 22.2%. In panels B and C, each number below an "S" (susceptible) or "r" (resistant) indicates the number of viruses classified along this path and the number in parentheses indicates the number of misclassifications. Fractional numbers are due to sequence heterogeneities resulting in more than one amino acid at the respective position.
|
|
|
|---|
Tenofovir resistance in multiple-dideoxynucleoside-resistant viruses. Among the 321 clinical samples analyzed, the viruses in 27 displayed a pattern of multiple dideoxynucleoside resistance. Of these, six samples contained viruses with an insertion of two amino acids (SV and SS in three viruses each) between RT positions 67 and 70, one sample contained a virus with a deletion of one amino acid in this region (M41L, deletion of amino acid 67, T69G, K70R, L74I, A98G, V108I, M184V, L214F, T215F, K219E), and 19 isolates displayed the Q151M mutation pattern (18). The virus in one sample contained a deletion, in addition to the Q151M complex (A62V, S68G, deletion of amino acid 69, mixture of 75V and 75I, F77L, Q151M, M184V, R211K, L214F). All viruses with insertions were resistant to tenofovir, with a median reduction in susceptibility of 16.8-fold (range, 10.0- to 19.3-fold). The isolate exhibiting a deletion in the context of the Q151M complex was susceptible to tenofovir (1.7-fold), whereas the isolate displaying the deletion without Q151M showed a 6.3-fold reduced susceptibility to tenofovir. The group of viruses with the Q151M complex could be divided into two subsets: the viruses from 12 samples showed a median reduction in susceptibility to tenofovir of 2.0-fold (range, 0.9- to 3.3-fold), whereas the viruses from 8 samples were highly resistant, with a median reduction in susceptibility of 13.5-fold (range, 9.3- to 17.0-fold). These subsets differed in the presence of K65R: none of the 12 susceptible viruses but 7 of 8 resistant viruses exhibited K65R (P < 0.0001, Fisher's exact test). Since the second subset contained five and two consecutive samples from two patients, respectively, statistics were recalculated by including the data for only one sample from each patient. The difference between the subsets was still significant (P < 0.05, Fisher's exact test).
Clonal analysis of K65R in the presence of Q151M. Since all tenofovir-resistant viruses carrying the Q151M mutation also exhibited mutations at position 62 (V in viruses from six samples, F in the virus from one sample, and P in the virus from one sample), in addition to K65R, we performed a clonal analysis to confirm the role of K65R and to verify the contribution of a mutation at position 62. An individual clone was isolated from the plasmid preparation of one of the five consecutive samples described above which had been generated during phenotypic resistance testing (the RT profile was A62V, K65R, D67N, S68K, T69N, V75I, F77L, K103N, V108I, Y115F, F116Y, Q151M, I178L, M184V, R211K, L214F, K219E). Mutations K65R and A62V were consecutively removed from this clone by site-directed mutagenesis (Fig. 1). The clones with K65R were highly resistant to tenofovir, irrespective of the presence of A62V, whereas the clones with the wild-type sequence at position 65 were susceptible. Resistance could thus be attributed to the presence of K65R.
![]() View larger version (40K): [in a new window] |
FIG. 1. Tenofovir resistance of a viral clone derived from a patient sample (with RT mutations A62V, K65R, D67N, S68K, T69N, V75I, F77L, K103N, V108I, Y115F, F116Y, Q151 M, I178L, M184V, R211K, L214F, and K219E) and three clones whose sequences are identical at all positions except positions 62 and 65. Data are given as the medians of at least three independent runs for each phenotypic resistance test. IQR, interquartile range (the interquartile range was not available for mutant A62V-K65R).
|
|
View this table: [in a new window] |
TABLE 1. Tenofovir resistance of individual clones in the context of TAMs in the absence or presence of M184Va
|
|
View this table: [in a new window] |
TABLE 2. Tenofovir resistance of clinical isolates in the context of TAMs in the absence or presence of M184Va
|
|
|
|---|
However, these highly resistant variants are rare even in populations of patients who have been heavily treated in the past. Therefore, frequent mutational patterns conferring lower levels of tenofovir resistance may have a much greater influence on the success of tenofovir therapy in patients previously treated with an NRTI. An analysis of more than 5,000 clinical viruses (5) revealed that only 12, 4, and 1% exhibited more than 3.0-, 5.0-, and 10.0-fold resistance to tenofovir, respectively; and thus, it was argued that the degree of cross-resistance to tenofovir is low. The corresponding values for our data set are higher (34, 21, and 7%, respectively). The high prevalence of viruses with >10.0-fold resistance is caused by the preselection of viruses resistant to multiple dideoxynucleosides. The higher frequency of viruses with 4.0- to 10.0-fold resistance may be explained by the extensive treatment with NRTIs that our clinical cohort had previously received; i.e., more than half of the patients had received four or more NRTIs. Furthermore, resistance factors may not be directly comparable between different phenotypic platforms, although the IC50 for our nonresistant reference virus was very similar to those in three previous studies (12, 25, 28), and clones with the M184V mutation showed the same extent of hypersusceptibility determined in our assay (10, 12, 25).
In addition, previous data and the present data should be interpreted in light of the clinically relevant cutoffs for tenofovir resistance. Thus, a durable treatment response with a decline in the plasma viral load of more than 0.5 log10 copies/ml was observed in a cohort of HIV-infected patients for more than 48 weeks if <4.0-fold resistance to tenofovir was present at the baseline, whereas tenofovir therapy failed if >4.0-fold resistance was present (8). However, another recent study showed that the response to tenofovir was impaired even with >1.5-fold resistance (7), a cutoff which interferes with the interassay variabilities of phenotypic assays. In consideration of these cutoffs, the level of clinically relevant cross-resistance to tenofovir may be broader than was originally thought, in particular in patients previously heavily treated with NRTIs. This may be a reason for the lower than expected level of reduction of the viral load with tenofovir treatment in clinical studies (19). However, other mechanisms that could contribute to the durability of the antiretroviral response should be considered, such as less efficient removal of tenofovir by pyrophosphorolysis and nucleotide-dependent chain terminator removal (11).
Our data allow some interesting conclusions about TAM-induced tenofovir resistance and the clinical role of M184V-induced resensitization to be made. First, the level of tenofovir resistance increased gradually with the number of TAMs. Without the M184V mutation, the median resistance of the clinical samples rose to >4.0-fold for viruses with five or six TAMs, whereas one clone with only three TAMs (M41L, L210W, and T215Y) also reached this threshold. Thus, not all TAMs are of equal importance for tenofovir resistance, which should be evaluated in further studies. Concerning the effect of M184V in NRTI-treated patients, recent data showed a more extensive viral load reduction in patients carrying viruses with M184V (8). After the study had been completed, it became obvious that the same effect was observed with tenofovir and placebo, so the authors concluded that the occurrence of M184V indicated good compliance rather than a resensitizing effect (8). Our study adds some information to this point of interest. First, data on the tenofovir hypersusceptibilities of viral clones carrying M184V could be confirmed. Second, data from clonal studies also support the fact that M184V resensitizes clones with TAM-induced resistance (5), although the resistance of viruses with the M41L, L210W, and T215Y mutations was only partially reversed to a level which may still be clinically relevant (Table 1). Third, stratified analysis of clinical isolates showed significantly lower levels of tenofovir resistance in viruses carrying TAMs with 184V than in viruses carrying TAMs with 184M (Table 2) but not in viruses with two TAMs. This subgroup contained four virus isolates from one patient with two TAMs (D67N and K219E) plus M184V, but also K65R and Q151M, which explains the higher median level of resistance. Overall, the resensitizing effect of M184V seemed to be less pronounced in the clinical isolates than in the clones. This was further supported by the fact that M184V did not play an important role in the decision trees. A possible explanation may be the accumulation of additional mutations such as those indicated in the mutual information profiles, e.g., V118I. Additionally, population sequencing of clinical isolates showed mutations from a pool of viruses, and this approach possibly failed to detect minorities of resistant variants.
The decision trees confirmed the central role of mutations at position 215 and 65 for 1.5- and 4.0-fold reduced susceptibilities to tenofovir, respectively. Importantly, these data were derived from viruses from tenofovir-naïve patients, thus describing cross-resistance rather than resistance to tenofovir in tenofovir-treated patients. Since the K65R mutation is rare in clinical samples, the bulk of cross-resistance to tenofovir appears to be due to TAMs. When the in vitro data and clinical data for tenofovir are combined, the resistance-response relationship appears to be more continuous, ranging from full response to an intermediate response to no response. The phenotypic cutoffs of 1.5- and 4.0-fold approximately define these response zones. A similar type of continuous response pattern has been observed for abacavir and will likely hold true for other antiretrovirals, which should be evaluated in further clinical studies.
Financial support for this study was obtained through grants from the Bayerische Staatsministerium für Kultus, Erziehung und Wissenschaft (to K. Korn), the Robert Koch-Institute, Berlin (National Reference Centre for Retroviruses), the Deutsche Forschungsgemeinschaft (to D. Hoffmann, R. Kaiser, and J. Selbig), and the Federal Ministry of Education and Research (HIV Competence Network grant AZ 01 KI 0211). This work was supported in part by the AIDS Reference Laboratory of Leuven, Belgium, which received funding from the Belgian Ministry of Social Affairs through the Health Insurance System.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»