Previous Article | Next Article ![]()
Journal of Virology, May 2009, p. 4051-4059, Vol. 83, No. 9
0022-538X/09/$08.00+0 doi:10.1128/JVI.02645-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Department of Immunology and Infectious Diseases,1 Harvard School of Public Health AIDS Initiative, Harvard School of Public Health, 651 Huntington Avenue, Boston, Massachusetts 021152
Received 23 December 2008/ Accepted 5 February 2009
|
|
|---|
|
|
|---|
Sub-Saharan Africa continues to bear the bulk of the burden of the HIV-1 epidemic, with 22.5 million infections out of the 33.2 million global infections (22). However, the developed world has begun to take notice by contributing unprecedented amounts of money and expertise to combat the epidemic. Much of the support for developing countries has gone toward drug treatment programs for those in need, making research into drug resistance of paramount importance. As more and more people get the treatment they require, drug resistance will increasingly play a major role in the epidemic in these hard-hit areas, where expertise and laboratory tests for clinical monitoring are often not available.
Some in vivo data point to different patterns of resistance between subtypes. Brenner et al. (3) have shown for subtype C that the 65R mutation occurs faster in cell culture in response to tenofovir pressure than it does in subtype B. A study from Brazil looked at 160 sequences from treatment-experienced HIV-1-infected individuals (21). It was found that subtype C viruses accumulated drug resistance mutations at a lower rate than subtype B isolates after at least 4 years of treatment.
Zidovudine (AZT) is a cheap antiretroviral and has relatively few side effects, making it a popular choice for a first-line regimen of HAART worldwide. While AZT drug resistance mutations have been extensively studied in subtype B, few studies have looked at non-B subtypes.
Resistance to AZT is conferred by thymidine analog mutations (TAMs). There are two pathways of TAMs that confer high levels of resistance with mutations in the viral reverse transcriptase (RT). TAM-1 is characterized by the 215Y mutation, along with the 41L and 210W mutations. TAM-2 is defined by the 215F mutation, along with 67N and 70R accessory mutations (14). In subtype B, TAM-1 occurs twice as often as TAM-2 (14, 24), and it is thought that viral fitness of the mutations plays a large part in this difference in frequency (10, 18). After years of treatment and different drug regimens, the mutations in the two pathways often become interchangeable due to the amount of pressure placed on the RT by continued antiretroviral therapy.
Viral fitness is a determinant of drug resistance mutation evolution. This is clearest when looking at the prevalence of nonnucleoside RT inhibitor drug resistance mutations. The most common nonnucleoside RT inhibitor mutations are 103N, 181C, 106 M, 188C, and 190A (6). Collins et al. (5) found that the fittest mutations are found more often on a population level, so that 103N is the fittest and most prevalent mutation. 181C, 190A, 188C, and 106M followed in order of fitness and prevalence. Several studies have looked at the viral fitness in different TAM pathways, as well. A different group found that the 41L215Y double mutant is much more fit than the 70R215Y double mutant, perhaps explaining the distinct separation of the two pathways in subtype B (10).
In subtype C, there is some evidence that the two pathways are not as definitive. Studies in our laboratory have shown that in a didanosine-containing regimen (along with AZT and efavirenz or nevirapine), 67N, 70R, and 215Y were the most common mutations seen together (17). The same study showed that the 41L mutation was rare compared to what was expected from the subtype B data; in fact, it did not occur in this study at all. In Brazil, where subtypes B and C cocirculate in the same population, a study found that the 41L and 210W TAM mutations occurred more frequently in subtype B than in subtype C (21).
These data led us to hypothesize that differences seen in the frequencies of TAMs in subtypes B and C are caused by different replicative costs associated with the pathways. We tested this by using Nef-tagged competing viruses in the same culture to determine the relative replication abilities of wild-type (WT) and mutant viruses. It was found that in subtype C, the 67N and 70R mutations enhanced HIV-1 replication compared to the WT, but no such enhancement was seen in subtype B. The TAM-1 double and triple mutants had no replicative cost for subtype C, but in subtype B, the TAM-1 mutations all had replication capacities similar to that of the single 215Y mutation alone. The TAM-2 pathway was the weakest set of mutations, with the triple mutant being the slowest virus in both subtypes. In subtype C, the mixed TAM pathway mutations, 67N70R215Y, were less fit than the WT but similar to the 215Y mutation alone; however, in subtype B, the triple mutant was less fit than the single mutation.
|
|
|---|
Cell culture. Peripheral blood mononuclear cells (PBMCs) were obtained from HIV-negative blood donors and isolated with a Ficoll gradient (MP Biomedicals, Solon, OH). The PBMCs were stimulated with phytohemagglutinin A (Sigma-Aldrich, St. Louis, MO) at a concentration of 10 µg/ml for 3 days in RPMI medium (Invitrogen, Carlsbad, CA) supplemented with 20% fetal bovine serum (FBS) (Invitrogen), 5 half-maximal units/ml human delectinized interleukin-2 (ABi, Columbia, MD), 100 U/ml penicillin G (Invitrogen), 100 U/ml streptomycin (Invitrogen), and 0.25 µg/ml amphotericin B (Fungizone) (Invitrogen). The PBMCs were maintained in RPMI medium supplemented with 20% FBS, 5 half-maximal units/ml interleukin-2, 100 U/ml penicillin G, 100 U/ml streptomycin, and 0.25 µg/ml amphotericin B.
293 cells were used to produce virus from DNA clones. The cells were maintained in Dulbecco's modified Eagle's medium (Invitrogen) supplemented with 10% FBS, 100 U/ml penicillin G, 100 U/ml streptomycin, and 0.25 µg/ml amphotericin B. Superfect transfection reagent (Qiagen, Valencia, CA) was used for DNA transfection with a modified protocol. Fifteen micrograms of plasmid DNA was diluted in a total volume of 150 µl of Dulbecco's modified Eagle's medium. Superfect (110 µl) was added to the DNA mixture and incubated for approximately 10 min; 1.25 ml of 293 maintenance medium was combined with the DNA, and the entire mixture was added to 5.28 x 106 cells in a T75 flask. The cells were incubated at 37°C for 3 hours and washed with Dulbecco's phosphate-buffered saline (Invitrogen). Fresh maintenance medium was added, and after 3 or 4 days, the viral titer was determined by p24 enzyme-linked immunosorbent assay (Perkin-Elmer, Waltham, MA). Cell-free viral stocks exceeding 50 ng/ml of p24 antigen were saved at –80°C.
Growth competition assay. Stimulated PBMCs (5 x 105) were infected at a total multiplicity of infection of less than 0.001. Mutant virus was competed with the subtype-specific WT virus at various initial concentrations because the replicative differences should not rely on the initial viral concentration in five different competitions. After an overnight infection, the supernatant was removed and 500 µl of fresh PBMC maintenance medium was added. On day 4 after infection, the culture was transferred to a new plate with 5 x 105 freshly stimulated PBMCs. The cells were fed on day 8 with one-half volume of fresh medium, and the culture was discarded on day 12. Viral RNA was isolated from the cell supernatant on days 0, 4, 8, and 12 with the Qiaamp Viral RNA Mini Kit (Qiagen) with the addition of an RNase-free DNase I (Qiagen) digestion step. All competitions were carried out in PBMCs from the same blood donor in five different concentrations. The WT, 67N, and 70R viruses were competed in five additional wells in PBMCs from a different blood donor for a total of 10 replicates.
AZT IC50. The AZT 50% inhibitory concentration (IC50) was determined by the protocol outlined previously (12) with some modification. PBMCs were placed in a 96-well plate at a concentration of 2 x 105 per well and stimulated with phytohemagglutinin A for 3 days. The cells were washed and resuspended in PBMC maintenance medium, and 50 µl of cells was added to each well. Fifty microliters of viral stock was also added to each well. After an overnight infection, 100 µl of AZT (Sigma-Aldrich) at concentrations from 0.002 to 200 µM was added to three wells for each dilution. The cells were fed on day 4 with fresh AZT solution, and on day 8, 50 µl of cell supernatant was saved at –80°C for p24 enzyme-linked immunosorbent assay determination. The IC50 was determined by nonlinear regression using SigmaPlot 2001 (Systat, San Jose, CA) by the Hill-Slope method. The expressed values are the percentages of growth at 100 µM AZT compared to growth with no AZT. This value was then normalized for each run with the WT, which was run with each set. Each mutant was tested in three different experiments in triplicate, except for the subtype C 67N70R215Y mutant, which was tested two times.
Quantitative PCR. Competing viral concentrations were determined with a quantitative PCR using sequence tags in nef, as described by Boutwell et al. (2). Briefly, the concentration of each viral species was determined with a unique forward primer and a common reverse primer using the Quantifast SYBR green RT-PCR master mix (Qiagen). The WT forward primer sequence is 5'-CAACACAGCCGCCAATA-3'; the Mut forward primer is CAACACTCCGGCGAACA-3'. The common reverse primer is 5'-CCCACAAATCAAGGATCT-3'. Relative amounts of virus were determined by comparison to a standard curve generated by a linearized plasmid of approximately full viral length, 9,078 bp.
Relative replication capacity estimation and statistical analysis.
The relative replication capacity was determined with the formula described by Maree et al. (15). Briefly,
![]() |
|
|
|---|
Quantification of competing viral species in the same culture. Placing two distinct viruses in the same culture is the best way to determine small differences in replication capacities. A dual-infection culture allows the control of well-to-well variation that confounds single-culture methods by establishing a control or WT virus. The two viral isolates are distinguished from each other with a sequence tag in nef. Four nucleotides are changed in the wobble positions of consecutive codons, which allows differential primer detection in the real-time PCR assay.
The replication cost of the sequence tag is determined by competing two viruses that differ only in the tag (Fig. 1A and B). Each growth competition is carried out in five different initial concentrations in PBMCs from the same blood donor. The populations show some variation but essentially remain the same over 12 days of infection, indicating that there is no advantage or disadvantage to the tag in nef in this system.
![]() View larger version (17K): [in a new window] |
FIG. 1. Representative growth curves of virus with sequence tags in nef. (A) The infecting ratio of virus is 80:20, with the native species (WT) at 80% and the tagged sequence (tagWT) at 20%. (B) Reciprocal infecting ratio of 20:80, with the WT at 20% and tagWT at 80% in the initial culture. (C) The ratio of WT virus to 215Y virus is 90% to 10%. (D) Reciprocal ratio of 10% WT to 90% 215Y virus. (E) Closer initial inoculum ratio of 60% WT to 40% 215Y. (F) Reciprocal of panel E, with 40% WT and 60% 215Y.
|
Single mutants in the TAM pathway. To investigate the replicative costs of TAMs, each mutation was added singly to both subtype B and C RTs. Each mutant infected PBMCs, along with its corresponding WT virus. Most single mutations behaved in the same manner regardless of subtype with two notable exceptions. In subtype B, the only mutations that were less fit than the WT were the 215Y and 215F mutations (P = 0.016 and P = 0.008, respectively); all other single mutations tested were similar to the WT. These values are similar to previously reported fitness indicators (18). In subtype C, the 215Y and 215F mutations were less fit than the WT (P = 0.008 for each), as well, and the 215Y results are shown in Fig. 2A and B. If the results for 215Y in subtype C are compared with the same graphs for subtype B (Fig. 1C and D), the results look similar, although the growth is somewhat delayed in subtype C.
![]() View larger version (25K): [in a new window] |
FIG. 2. (A) Competition with the subtype C WT (mjWT) and the 215Y mutant (mj215Y) with an initial inoculum of 90% WT to 10% mutant. (B) Reciprocal of panel A, with 10% WT to 90% 215Y. (C) Competition of the subtype C WT with the 67N mutant at a ratio of 90:10. (D) WT and 67N mutant at a ratio of 10:90. (E) 70R competition at a closer initial inoculum of 80% WT to 20% 70R. (F) Opposite of panel E, with 20% WT to 80% 70R. (G) Summary of all single-mutant competitions. Relative fitness is defined in Materials and Methods, and each mutant was competed with its subtype-specific WT. The bars represent the means, and the error bars are 1 standard deviation. The bars marked with asterisks are significantly different from the WT value according to the Mann-Whitney rank sum test, with a P value of less than 0.01.
|
For the purposes of comparing many different mutants, it is best to use a summary statistic (23). The formula defined by Maree et al. (15) is a comparison of the difference in the growth rates of two viral species, which they call the coefficient of selectivity (s). Relative viral fitness is s + 1. For the comparison between different mutants, the s + 1 value is used and graphed (Fig. 2G). The relative fitness of the WT is defined as 1 because there is no advantage or disadvantage to the tag. More fit viruses will have values greater than 1, and less fit viruses will have values less than 1. All values are averages of five replicates, except for the WT, 67N, and 70R which have 10 replicates, and significance is determined by a P value of less than 0.01 as determined by the Mann-Whitney rank sum test.
Double and triple mutations. Double mutants in the TAM pathway were tested as shown in Fig. 3. There are several interesting differences seen between the two subtypes in the relative costs of mutations along the TAM-1 pathway. The double mutant 41L215Y in subtype C does not significantly differ from WT levels of replication (Fig. 3A and B); however, in subtype B, it is less fit than the WT (P = 0.008). When the double mutant was compared to 215Y in subtype B, there was no significant change in fitness relative to 215Y alone (Fig. 3E).
![]() View larger version (25K): [in a new window] |
FIG. 3. (A) The subtype B WT competed with the 41L215Y double mutant at a ratio of 20:80, respectively. (B) The subtype C WT competed with the 41L215Y double mutant at the same initial ratio as in panel A. (C) The subtype B competition with the 70R215F mutant at a ratio of 40% WT to 60% mutant. (D) The 70R215F competition in subtype C at a ratio of 40% WT to 60% mutant. (E) Summary of all double-mutant competitions, with the 215Y and 215F single mutants included as a reference. The 41L215Y mutant is in the TAM-1 pathway. The 70R215F mutant is in the TAM-2 pathway, and the 70R215Y mutant is in the mixed TAM pathway. The bars are the means of five replicates, and the error bars are 1 standard deviation. All bars marked with asterisks have a P value of less than 0.01 according the Mann-Whitney rank sum test compared to the subtype-specific WT. The bars marked with a plus have a P value of less than 0.01 according to the Mann-Whitney rank sum test compared to each other.
|
The triple mutants show more differences between the two subtypes abilities to tolerate drug resistance mutations. The subtype C TAM-1 triple mutant, 41L210W215Y, has no replication capacity cost compared to the WT (Fig. 4B). In subtype B, the TAM-1 triple mutant is less fit than the WT (P = 0.032) (Fig. 4A) but has a relative fitness similar to that of the 215Y mutation alone. The TAM-2 triple mutant, 67N70R215F, confers a similar effect on each subtype. For both subtypes, the triple mutants are less fit than the WT and relatively less fit than the 215F mutation alone (P = 0.008 for each comparison). In fact, these mutations had the highest replicative cost of all mutants tested (Fig. 4E).
![]() View larger version (24K): [in a new window] |
FIG. 4. Relative fitness of triple mutants. (A) Subtype B competition of the TAM-1 mutant, 41L210W215Y, with the WT at a ratio of 20% WT to 80% mutant. (B) Subtype C TAM-1 competition with the WT at the same ratio of 20% WT to 80% mutant. (C) Mixed TAM, 67N70R215Y, competition with WT at an inoculum of 40% WT to 60% mutant in subtype B RT. (D) Subtype C competition of WT to mixed TAM at a ratio of 40:60, respectively. (E) Summary of triple-mutant competitions, with the 215Y and 215F mutations included for reference. The bars are the means of five different competitions, and the error bars are 1 standard deviation. All bars marked with asterisks differ significantly from WT, with a P value of less than 0.05 as determined by the Mann-Whitney rank sum test. All bars marked with +, #, or ^ are comparisons between the two marked bars with a P value of less than 0.01.
|
IC50. The IC50s for the WT, 215Y, 215F, and all triple mutants were determined using the standard protocol of comparing viral growth after 8 days with AZT to growth without the drug (12). Due to the nature of the assay, the values tended to be variable because AZT resistance mutations did not always shift the inhibition curve to the right (data not shown). However, the inhibition curve was shifted up, implying that AZT could not restrict the growth of the mutants completely, even at concentrations as high as 100 µM. Figure 5 illustrates the amount of viral growth expressed as a percentage of the growth if no AZT was present and was normalized by the WT value to give the resistance at 100 µM. Most mutants were different from the WT, according to the Mann-Whitney rank sum test, but there was no difference between individual mutants. The 41L215Y double mutants were not different from the WT, but this may be explained by the variability of the assay.
![]() View larger version (15K): [in a new window] |
FIG. 5. AZT IC50 results expressed as the resistance over WT at a 100 µM concentration. The bars represent the means, and the error bars are 1 standard deviation. TAM-1 is the 41L210W215Y mutant, TAM-2 is the 67N70R215F mutant, and mTAM is the 67N70R215Y mutant. The asterisks signify a P value of less than or equal to 0.001 as determined by the Mann-Whitney rank sum test compared to the subtype-specific WT.
|
![]() View larger version (5K): [in a new window] |
FIG. 6. Competition of the subtype B and C WTs against each other. (A) Competition with an inoculum of 80% subtype B, hxWT, and 20% subtype C, mjWT. (B) Opposite of panel A, with an inoculum of 20% subtype B and 80% subtype C.
|
|
|
|---|
TGG), which if done in the wrong order can result in a stop codon. The genotype seen by our group in Botswana (17), which we call the mixed TAM pathway of 67N70R215Y, showed a less striking difference between the subtypes. The triple mutants in subtypes B and C were each just as fit as the 215Y mutant alone; however, only in subtype C was there a benefit for this mixed pathway compared to TAM-2. When the single mutations were considered alone, subtype C showed enhanced replication of the 67N and 70R mutations compared to the WT. This led to the hypothesis that these mutations may occur first, and adding a 215Y mutation to either of the other residues results in a more competent virus than if a 215F mutation occurs instead. Mutations occur randomly, but the more fit amino acids are more likely to persist and be seen clinically due to their ability to outcompete the rest of the viral reservoir. These data do not explain why this pathway was seen more than TAM-1 or TAM-2 in Botswana, but the published results (17) are based on a restricted number of patients.
The TAM-2 pathway was the least fit of the triple mutants in both subtypes. It is somewhat puzzling that this was the second most common set of mutations seen in the Tshepo study (17), but the observed genotypes may be driven by the replicative fitness of the 67N and 70R mutations. These mutations may be selected early, thereby negating the possibility of the TAM-1 pathway.
The Novitsky et al. study is based on virological-failure cases from a clinical trial in which all patients described in the paper were on the AZT/didanosine/efavirenz or nevirapine arm of the trial. Their unusual genotype of the mixed TAM pathway, 67N70R215Y, could be a result of the drug combination. However, when the results were compared with subtype B sequences from patients on a similar regimen, there was a significant difference in the frequency of the 41L and 210W mutations between the two subtypes (17). Whether there will be a similar pattern of resistance on the most common drug regimen of AZT/lamivudine remains to be seen.
A study done in India found that 50% of HIV-1 subtype C strains containing nonnucleoside RT inhibitor drug resistance mutations included the 41L mutation (19). In addition, there was only one instance of the mixed pathway as seen in the African study, but it came from an isolate that included viral clones with all of the accessory TAMs, including the 219Q mutation. The cohort in India was made up of individuals on diverse kinds and numbers of drug regimens for different amounts of time, whereas the Novitisky et al. report was from a drug treatment study in which each participant received only one drug regimen for a prolonged period. The results discussed here may reconcile these two studies. While the TAM-1 pathway is clearly the fittest pathway, early in treatment, the 67N and 70R mutations may evolve first due to their enhanced replication capacities in HIV-1 subtype C, and it may take continued drug pressure for the TAM-1 mutations to evolve. It will be interesting to see if there is an increase in frequency relative to subtype B of the TAM-1 and mixed- pathway genotypes in response to the most common drug regimen of AZT and lamivudine.
A growing concern is the transmissibility of drug resistance mutations. With an increased replicative capacity of 67N and 70R compared to the WT and the absence of a fitness cost to the TAM-1 mutations, one would expect to see these mutations be easily transmissible through all routes of infection. Several studies have been carried out to determine the amount of transmitted resistance in areas where subtype C is endemic (1, 4, 11, 13, 20). They found that transmitted resistance is rare, even though our results predict that some mutations would be transmitted more easily; however, the use of HAART in southern Africa has lagged for at least a decade with respect to developed countries. With the increasing roll-out of antiretroviral drugs worldwide, it will be interesting to determine if the 67N and 70R single mutations and the fully resistant TAM-1 mutants are transmitted more and persist longer in a subtype C background.
As treatment continues in the developing world, it is important to understand how subtypes may differ in their responses to antiviral therapy. We have shown here that HIV-1 subtypes B and C differ in their abilities to tolerate drug resistance mutations. In the future, it will be important to learn if the differences between subtypes can be seen in clinical samples and with clinical consequences.
We thank C. Boutwell and M. F. McLane for scientific advice, R. Bosch for statistical advice, and L. Melton for editorial assistance.
Published ahead of print on 18 February 2009. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»