Previous Article | Next Article ![]()
Journal of Clinical Microbiology, June 2006, p. 1930-1943, Vol. 44, No. 6
0095-1137/06/$08.00+0 doi:10.1128/JCM.02415-05
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Hulin Wu,4 and
Lisa M. Demeter1*
Department of Medicine, University of Rochester School of Medicine and Dentistry, Rochester, New York,1 Department of Pathology, University of Utah School of Medicine, Salt Lake City, Utah,2 Department of Hematology and Oncology, UCLA AIDS Institute, Los Angeles, California,3 Department of Biostatistics, University of Rochester School of Medicine and Dentistry, Rochester, New York4
Received 21 November 2005/ Returned for modification 20 January 2006/ Accepted 27 March 2006
|
|
|---|
|
|
|---|
All assays that measure HIV-1 replication efficiency in cell culture compare the clinical (or test) virus to a reference HIV-1 isolate. These assays vary widely, although they differ primarily in at least one of four different characteristics. First, such assays may either measure the replication efficiency of virus cultured directly from the patient sample or a recombinant virus constructed by combining a PCR-amplified segment of the patient virus genome with the sequence backbone from a laboratory strain of HIV-1. Second, these assays may measure the replication efficiency during a single virus replication cycle or over multiple cycles. Third, these assays may measure virus growth directly, using a viral gene or gene product, or measure virus growth indirectly by assaying a reporter gene substituted for a nonessential viral gene. Fourth, assays of HIV-1 replication efficiency may evaluate the growth of test and reference viruses in the same or separate cultures (referred to as growth competition and parallel infections, respectively) (28, 31).
Growth competition assays are considered to be the preferred method for measuring HIV-1 replication fitness in cell culture because they are more sensitive to subtle differences in replication efficiency than parallel infections and are not subject to artifact due to differences in culture conditions (10, 30). Growth competition assays require a method to distinguish the two virus variants, usually by measuring the relative proportions of each mutant (or a linked reporter gene) using bulk DNA sequencing, clonal analysis, heteroduplex tracking assay, or real-time PCR (16, 22, 26, 27, 32, 33, 36, 40). These methods make growth competition assays more labor-intensive than parallel infections, thus limiting their use in larger clinical trials.
We have designed a multiple-cycle, recombinant-virus, growth competition assay that obviates the need to purify and analyze HIV-1 RNA or DNA. Two recombinant HIV-1 virus constructs were designed to express either the mouse Thy1.1 or the mouse Thy1.2 genes in place of HIV-1 nef. The proportions of test and reference HIV-1 isolates were measured by using flow cytometry to detect the Thy1.1 or Thy1.2 gene products expressed by each virus on the surface of infected cells. Here we show that this assay can be used to reliably quantify the relative replication efficiencies of HIV-1 drug-resistant mutants, including those derived from clinical samples.
|
|
|---|
Cell culture. PM1 cells were grown in the presence of RPMI (Cellgro, Herndon, VA), supplemented with 10% FBS, L-glutamine (2 mM), penicillin (100 U/ml), and streptomycin (100 U/ml). 293 cells were grown in Dulbecco modified Eagle medium with 10% FBS, penicillin (100 U/ml), and streptomycin (100 U/ml).
Construction of pAT2 and pAT1 HIV-1 vectors.
Construction of the pAT1 and pAT2 vectors shown in Fig. 1 was accomplished in three separate steps. (i) The BamHI-NcoI fragment of the Thy1.2 gene of pNLThy was subcloned into pQE-31 and mutagenized at codon 108 (CAA
CGC) using PCR-mediated mutagenesis (QuikChange; Stratagene, La Jolla, CA). After we verified the absence of spurious mutations, this fragment was cloned back into the pNLThy vector to create the Thy1.1-expressing pNL4-3 vector pNLThy1.1. (ii) In order to facilitate the cloning of reverse transcriptase (RT) sequences amplified from patient plasma, we then introduced silent XmaI and XbaI sites near the beginning and end of the RT coding region by subcloning the SpeI-EcoRI fragment from the pNL4-3XX construct described previously (14). (iii) We next removed a second XbaI site which had been introduced between the 3' long terminal repeat and pUC19 bacterial plasmid sequence when Thy1.2 was originally introduced into the pNL4-3 vector (29). Removal of the additional XbaI site was accomplished by digesting the pNLThy vector with NcoI and XbaI, filling in the ends with T4 polymerase, and religating. The AatII-MluI fragment containing the deleted XbaI site was then subcloned into the pNLThy1.1 and pNLThy to generate pAT1 and pAT2, respectively (Fig. 1).
![]() View larger version (13K): [in a new window] |
FIG. 1. Maps of pAT1 (A) and pAT2 (B). The open reading frame locations for gag, pol, env, the accessory proteins, Thy, and the long terminal repeats (LTRs) are shown as black bars. The region containing the bacterial plasmid pUC19 sequence is marked as pUC. The 5' cellular flanking sequence is also marked by a black bar. The locations of restriction enzyme sites used to construct each vector are each marked with a vertical line and labeled with the nucleotide position. See Materials and Methods for details on the construction of each vector.
|
The antisense primers used for mutagenesis were the reverse and complement of each sense primer. The resultant clones were sequenced to confirm the presence of the mutation and the absence of spurious mutations. PR inhibitor resistance mutations were introduced into pAT1 and pAT2 by subcloning the ApaI-XmaI fragment from the mutant pRHAXX constructs, and the G190S mutation was introduced into pAT1 and pAT2 by subcloning the XmaI-XbaI fragment from the mutant pRHAXX. The V106A and P236L mutations were introduced separately into pAT1 and pAT2 by subcloning the SbfI-AgeI fragment from pNL4-3(V106A) (3) and the XmaI-XbaI fragment from pNL4-3XXP236L (14). All mutant pAT1 and pAT2 clones were sequenced to verify the presence of the appropriate mutation, the integrity of the cloning sites, and the presence of the correct Thy allele.
Amplification and cloning of HIV-1 PR and RT sequences from patient plasma. Plasma was obtained from HIV-infected patients seen at the University of Rochester Infectious Diseases Clinic, after obtaining informed consent, according to the guidelines of the University of Rochester Research Subjects Review Board. Viral RNA was purified from plasma, amplified by RT-PCR, and sequenced by using the Viroseq HIV-1 Genotyping System V2.0 (Celera, Rockville, MD). Amplification of viral RNA for cloning was performed with RNA purified from the Viroseq kit and the QIAGEN One-step RT-PCR kit (QIAGEN, Valencia, CA). The following primers were used, which contain ApaI and AgeI sites (underlined) to facilitate cloning: GagApaI-sense, 5'-GCAGGGCCCCTAGGAAAAAGGGCTGTTGGAAATGTGGAAAGGAAGG; and Age I-antisense, 5'-ACCGGTTCTTTTAGAATCTCCCTGTTTTCTGCC. The cycling conditions were as follows: reverse transcription (50°C for 50 min), PCR activation (95°C for 15 min), PCR cycling (40 cycles of 94°C for 30 s, 60°C for 30 s, and 72°C for 2 min), and a final extension (72°C for 10 min). The ApaI-AgeI product includes the 3' end of the gag coding region (codons 29 to 55 of p7 and all of p1 and p6), as well as the complete coding region of protease and the first 314 amino acids of the RT. The PCR product was gel isolated and cloned into pCR4-TOPO (Invitrogen, Carlsbad, CA). Clones were sequenced to ensure the absence of frameshift mutations and to characterize the PR and RT sequence. In order to generate recombinant infectious vectors, the ApaI-AgeI region was subcloned from the pCR4-TOPO constructs into pAT1. pAT1 clones were sequenced from the ApaI site to the AgeI site to be sure they were identical to the TOPO version of the clone. In addition, a panel of restriction digests were performed to ensure that no significant deletions or rearrangements occurred elsewhere in the genome due to a cloning artifact.
Generation of virus stocks. pAT1 and pAT2 containing wild-type NL4-3, mutant NL4-3, or patient viral sequences were purified from Escherichia coli and used to transiently transfect 293 cells (Superfect; QIAGEN). Because 293 cells have no CD4 receptor, progeny virus cannot reinfect these cells. This minimizes the introduction of spurious mutations and the loss of mutations due to reversion. Supernatants were harvested 72 h later and clarified by centrifugation at 500 x g. HIV-1 virus capsid protein (p24) quantitation was performed on virus stocks by using an HIV-1 p24 enzyme-linked immunosorbent assay (Perkin-Elmer, Norwalk, CT).
Growth competition assay. Seven million PM1 cells, pretreated with 10 µg of Polybrene (Sigma; St. Louis, MO)/ml in phosphate-buffered saline (PBS), were coinfected with a total of 300 ng of two virus stocks, one expressing Thy1.1 and the other expressing Thy1.2. The relative proportions were varied in different experiments and estimated based on the p24 antigen content of the virus stocks. Virus and cells were incubated at 37°C in a total volume of 3.5 ml for 1 h. To remove unbound virus, cells were washed with PBS and centrifuged at 500 x g. Cells were then seeded at 200,000 cells/ml in medium and cultured at 37°C. On days 3, 4, 5, and 6, half the culture was removed and replaced with fresh medium. A subset of cells removed during this process was stained with trypan blue and counted to determine the number of viable cells in the total culture. Cells were then stained with the Thy1.1 and Thy1.2 antibodies, either alone or in combination. The Thy1.1 and Thy1.2 antibodies were diluted 1:200 and 1:100, respectively, in wash buffer before use (Dulbecco PBS [Invitrogen], with 0.02% FBS [Valley], 0.02% sodium azide, 0.5 mM EDTA [Gibco]). Cells were fixed (4% paraformaldehyde [Sigma, St. Louis, MO] in Dulbecco PBS [Invitrogen]; pH 7.4 to 7.6) before analysis with a FACSCalibur (Becton Dickinson, San Jose, CA). A total of 100,000 events were collected from each stained sample. Cells were also collected and frozen as nonviable cell pellets for sequence analysis. Culture supernatant was collected and used for p24 quantitation.
Sequence analysis. Genomic DNA was extracted from 2 million cells by using the QIAGEN DNA Blood Minikit (QIAGEN). DNA was eluted with 200 µl of nuclease-free water. A 50-µl PCR was performed with 10 µl of genomic DNA and the following reagents: 1 U of Platinum Taq HiFi (Invitrogen, Carlsbad, CA), 0.2 mM deoxynucleoside triphosphates (PE Applied Biosystems, Foster City, CA), 2 mM MgSO4, and 0.2 µM primer RT18 and primer RT21 (the sequences are shown below). The following cycling conditions were used: 94°C for 2 min, 40 cycles of 94°C for 15 s, 55°C for 30 s, and 68°C for 2 min. PCR products were sequenced by using BigDye v3.1 (PE Applied Biosystems) and the primers RT18 and RT127 for protease codon 90, RT60 and 215D for RT codons 103 and 106, the primers RT162 and 215D for RT codon 190, and the primers RT162 and RT21 for RT codon 236.
The PCR and sequencing primers were as follows: RT18, 5'-GGAAACCAAAAATGATAGGGGGAATTGGAGG; RT21, 5'-CTGTATTTCTGCTATTAAGTCTTTTGATGGG; RT60, 5'-CTGAAAATCCATACAATACTCC; RT127, 5'-CTAGGTATGGTAAATGCAGT; RT162, 5'-AAAGGATCACCAGCAATATTCC; and 215D, 5'-TCTGTATGTCATTGACAGTCCAGC.
The relative prevalence of mutant variants at each time point was quantitated by averaging the relative peak heights from the sense and antisense sequencing reactions from a minimum of two replicate infections. We have found that this approach correlates with relative proportions as measured by clonal analysis (40).
Analysis of flow data. Flow cytometry data were analyzed by using CellQuest software (Becton Dickinson, San Jose, CA). Viable cells were gated based on their forward-scatter and side-scatter profiles. Gated cells were then analyzed by using a density plot of FL1 (FITC-Thy1.1) versus FL2 (PE-Thy1.2). Each culture was stained with either the anti-Thy1.1 or anti-Thy1.2 antibodies alone or with both antibodies. The gate defining positivity for the Thy1.1-stained sample was set using the Thy1.2 antibody-stained sample as a negative control, with a threshold of <0.05%. The same was done for the gate defining positivity for Thy1.2 using cells stained with the Thy1.1 antibody alone as a negative control. These gates were then used to determine the proportions of Thy1.1-infected cells, Thy1.2-infected cells, and dually infected cells in the dual-antibody-stained sample. The total percentage of Thy1.1-infected cells for each sample was determined by adding the percentage of singly stained (Thy1.1) positive cells and dually stained (Thy1.1 and Thy1.2) positive cells. The same was done for the total percentage of Thy1.2-positive cells. In order to minimize variability in setting gates, all gates were set by one observer, who was blinded as to the PR and RT genotype of the competing variants.
Fitness calculation and statistical analysis. The replication efficiency of each variant was estimated by using two time points according to the method of Maree et al. (25). We have shown that relative fitness, as estimated using this method is different than the selection coefficient (s), used in population genetics, which reflects the number of progeny contributed to the next generation by a mutant (41). We have therefore referred to the replication efficiency as estimated by Maree and coworkers as the production rate ratio or "p" (41).
We used the equation: p = ln[H(T)/H(0)]/(ln[W(T)/W(0)] +
T). W(T)/W(0) is the fold expansion of wild-type virus, which was calculated by dividing the number of wild-type infected cells in the culture at time T by the number at time zero. The number of viable wild-type or mutant infected cells was calculated by multiplying the number of viable cells in the total culture by the percent as determined by flow cytometry. H(T)/H(0) is the fold change in the mutant to wild-type ratio over T days of the experiment. This was calculated by determining the proportion of mutant to wild-type at time T (the number of cells infected by mutant virus divided by the sum of mutant and wild-type infected cells), divided by the same proportion at time zero. Delta (
) is defined as the life span of an infected cell and was assumed to be 0.5/day (25).
Three criteria were used to determine which two time points from each experiment to include in the calculation of p: (i) the percentage of dually infected cells, defined as cells which express both Thy1.1 and Thy1.2, should not be significant (defined as
20% of the less-well-replicating virus); (ii) measurable growth of the wild-type virus should be detected (defined as a positive slope in a plot of log percent infected cells versus time); and (iii) the percentage of cells infected by each variant should be above background (
0.05%). If more than one combination of time points satisfied these three criteria, then the longest and latest window of time was chosen.
Each mutant was competed against a reference strain in at least three separate infections, and the values of p were averaged. The averages of the production rate ratio, p, for each mutant were compared to K103N by using the t test. A significant P value was corrected for nonorthogonal multiple comparisons by using the Bonferroni adjustment.
|
|
|---|
Specificity and linearity of staining. Figure 2 shows the staining of cells coinfected with K103N-AT1 and wild-type-AT2 viruses at day 6, stained with either anti-Thy1.1 antibody alone (Fig. 2A), anti-Thy1.2 antibody alone (Fig. 2B), or both antibodies (Fig. 2C). Staining with the fluorescent-labeled Thy1.1 and Thy1.2 antibodies is specific: we found no staining of infected cells using isotype control antibodies (data not shown) and detected no staining when each Thy antibody was used to stain cells infected by the heterologous virus (Fig. 2A and B). Thresholds to distinguish positive- from negative-staining cells were set using the heterologous anti-Thy antibody as a control. Thresholds using singly stained samples were set so that <0.05% of Thy1.1-infected cells were positive when stained with the Thy1.2 antibody and vice versa. The percent infected cells, as determined by the dually stained sample, was used to calculate the proportions of the wild type and mutant. The percentage of cells infected by both viruses (upper right quadrant, Fig. 2C) was added to both the singly infected Thy1.1 positive and the Thy1.2 cells. All data presented here use gates set by one observer, who was blinded as to the genotype of the Thy1.1 and Thy1.2 viruses. We noted some minor differences in quantitation of the two viruses depending on whether single- or dual-stained samples were used. Therefore, we always compared the two variants using the dual-stained sample.
![]() View larger version (30K): [in a new window] |
FIG. 2. Density plots of infected cells. PM1 cells were coinfected with K103N-AT1 and wild-type AT2 viruses. Cells were stained at day 6 with either FITC-Thy1.1 antibody alone (A), PE-Thy1.2 antibody alone (B), or both antibodies (C). The log fluorescence intensities of FL1-H (FITC-Thy1.1) and FL2-H (PE-Thy1.2) are shown on the x axis and y axis, respectively. The horizontal and vertical lines represent gates that define positivity. The number of cells infected with K103N virus, wild-type virus, and both viruses is shown as a percentage in the upper left, lower right, and upper right quadrants, respectively.
|
![]() View larger version (15K): [in a new window] |
FIG. 3. Linearity of flow staining. PM1 cells were singly infected with Thy1.1- or Thy1.2-expressing wild-type virus. Cells from the Thy1.1-infected culture were counted and diluted with cells from the Thy1.2-infected culture, and the reciprocal was done for the Thy1.2-infected culture. The following dilutions of cells were then stained with either the Thy1.1 or the Thy1.2 antibodies: undiluted, 1:2, 1:3, 1:5, 1:10, 1:30, 1:50, 1:100, 1:500, and 1:1,000. The x axis represents the percentage of infected cells based on dilution. For the Thy1.1 data, this number was obtained by taking the percentage of Thy1.1-positive cells from the undiluted sample as measured by flow and dividing it by each dilution. The same was done for the Thy1.2 data. The y axis represents the percentage of infected cells as measured by flow cytometry. , Cells infected with wild-type Thy1.1; , cells infected with wild-type Thy1.2. The solid and dotted lines represent the best-fit linear curves for Thy1.1 and Thy1.2, respectively. The slopes and R2 values are shown next to each curve.
|
![]() View larger version (17K): [in a new window] |
FIG. 4. Correlation of p24 quantitation with flow analysis. Seven million PM1 cells were coinfected with the following combinations of virus: wild-type AT1/G190S-AT2, wild-type AT2/G190S-AT1, wild-type AT1/K103N-AT2, and wild-type AT2/K103N-AT1. A total of 300 ng of p24 was used at either equal ratios or ratios favoring the mutant. On days 3, 4, 5, and 6 cells were removed from the culture for staining, and supernatant was removed for p24 antigen quantitation. The x axis represents the p24 concentration in pg/ml. The y axis represents the percentage of infected cells obtained by flow analysis. The data from both Thy alleles were combined into one plot for each mutant competition. , G190S versus wild-type competitions; , K103N versus wild-type competitions. The solid and dotted lines represent the best-fit linear curve for the G190S and K103N competitions, respectively.
|
![]() View larger version (16K): [in a new window] |
FIG. 5. Competition assay of wild-type Thy1.1 versus wild-type Thy1.2. Seven million PM1 cells were coinfected with wild-type AT2 and AT1. A total of 300 ng of p24 was used at a ratio of 150 ng (50%) of Thy1.2 to 150 ng (50%) of Thy1.1. (A) The y axis represents the percentage of cells infected on a logarithmic scale (left side) or the corresponding total number of infected viable cells (right side) versus time after infection in days. (B) The percentage of the culture contributed by each genotype was calculated by dividing the percentage of infected cells (or the number of infected cells) determined for each virus and dividing it by the sum of the percent infected cells (or the number of infected cells) for both viruses and then multiplying that value by 100%.
|
Using these mutant viruses, this assay can distinguish the K103N mutant, which has a relatively preserved replication efficiency, from G190S, which has a replication efficiency substantially reduced compared to the wild type (Fig. 6). For example, the growth rate of G190S is substantially slower than the wild type in the same culture, even though the mutant was initially at a higher prevalence than wild-type virus (Fig. 6A and B). No such difference in replication rates can be seen with the K103N mutant over the 6-day course of the experiment (Fig. 6C and D).
![]() View larger version (14K): [in a new window] |
FIG. 6. Competition of G190S and K103N versus wild type. Seven million PM1 cells were coinfected with wild-type AT2 and G190S-AT1 or wild-type AT2 and K103N-AT1. A total of 300 ng of p24 was used at a ratio of 225 ng (75%) of mutant to 75 ng (25%) of wild type (WT). The x axis represents time in days. The y axis represents either the percentage of infected cells on a logarithmic scale for G190S versus the wild type (A) and K103N versus the wild type (C) or the proportions of mutant and wild type for G190S versus the wild type (B) and K103N versus the wild type (D). The dotted line represents the mutant, and the solid line represents the wild type. The proportion of each genotype was calculated by dividing the percent infected cells or the number of infected cells determined for each virus and dividing it by the sum of the percent infected cells or the number of infected cells for both viruses.
|
![]() View larger version (22K): [in a new window] |
FIG. 7. Proportion of mutant over time measured by flow or sequencing. Seven million PM1 cells were coinfected with wild-type and one of the following mutant viruses: L90M, K103N, G190S, or P236L. A total of 300 ng of p24 was used at a ratio equal to or favoring the mutant, either 50:50 or 75:25. The x axis represents time in days. The y axis represents the percent mutant. (A) L90M; (B) K103N; (C) G190S; (D) P236L. The dotted line represents the percent mutant as measured by sequencing, and the solid line represents the percent mutant as measured by flow.
|
![]() View larger version (10K): [in a new window] |
FIG. 8. Correlation of the proportion of infected cells, determined by flow cytometry, with sequence analysis. The proportions of mutant determined at each time point were combined into one plot of proportion measured by flow (x axis) versus proportion measured by sequencing (y axis).
|
We observed that certain experimental conditions led to higher variability in the production rate ratios for the same mutant. These conditions included the presence of significant proportions of cells that were dually infected with both viruses, significant lags in growth of wild-type virus, and levels of infected cells that were at or below our defined background threshold of 0.05%. As a result, we defined specific criteria for which datum points are included in the calculation of the production rate ratio. We observed that after day 6 the number of susceptible cells is limiting enough that significant numbers of dually infected cells occur. Since dually infected cells can produce recombinant viruses and since we have not evaluated the effect of dual infection on the competition between two viruses, we did not include time points in which there was significant dual infection (defined as
20% of the less-well-replicating virus). We also required that there should be detectable growth of the wild-type virus (defined as a positive slope in a plot of log percent infected cells versus time) and that the percentage of cells infected by each variant should be
0.05%. These guidelines were established during early experiments with the wild-type and K103N variants and were formulated before we initiated this analysis of the different mutants.
Using these criteria, we calculated the production rate ratio, p, for each of the mutants in both the Thy1.1 and the Thy1.2 backgrounds (Table 1). In order to limit the number of comparisons performed, we compared the value for each mutant to that of K103N, using it as an example of a mutant with relatively preserved fitness. The D30N, G190S, and P236L mutants all had significantly negative production rate ratios compared to K103N in both backgrounds, whereas L90M was not significantly different. The reduction in replication efficiency for V106A was significantly reduced compared to K103N when V106A was placed in the Thy1.2 background, but this reduction only trended toward statistical significance in the 1.1 background, when we corrected for multiple comparisons (P = 0.027). With the exception of the D30N mutant, all values for replication efficiency tended to be more negative when assayed in a Thy1.2 background compared to the Thy1.1 background. Thus, there may be a slight improvement in the relative fitness of the AT1 viruses relative to AT2, which we did not detect in our wild-type-wild-type competitions, that could have led to an obscuring of the fitness reduction conferred by V106A in the Thy1.1 background.
|
View this table: [in a new window] |
TABLE 1. Values for the production rate ratio (p) for site-directed mutants of HIV-1 PR and RT
|
Table 2 provides a summary for the patient clones tested. The sequences of p6, PR, and RT show that each clone had a unique genotype. Two of seventeen clones produced virus stocks with p24 concentrations too low to test in the competition assay (clones 3-1 and 5-1). These defects in p24 production may be due to one or more mutations in gag or protease that are unique to these clinical viral sequences or may reflect incompatibility between these clinical viral sequences and that of NL4-3. We think this is probably not due to a cloning artifact because we performed extensive restriction digests and sequenced the clinical isolate insert through the cloning sites, but we cannot completely rule out a small deletion or duplication in the viral vector sequence that may have been missed by our restriction enzyme digests and sequence analysis (see Materials and Methods). Of the 15 clones tested in culture, 7 did not establish a productive infection in the growth competition assay, with productive infection defined as a percentage of infected cells above baseline (
0.05%) at a minimum of 2 time points (clones 1-1, 2-1, 3-4, 4-1, 4-2, 5-2, and 5-3). Again, this may be due to deleterious mutations or, less likely, a cloning artifact. As can be seen in Table 2, when three to four clones per sample were tested, all patient samples yielded at least one clone that could be assayed and which resulted in a productive infection.
|
View this table: [in a new window] |
TABLE 2. Mutationsa p24 titer, and production rate ratios for clones derived from patient samples
|
![]() View larger version (17K): [in a new window] |
FIG. 9. Replication of patient clones. Seven million PM1 cells were coinfected with a total of 300 ng of p24 at a ratio of 225 ng (75%) of patient clone to 75 ng (25%) of wild type (WT). The x axis represents time in days. The y axis represents either the number of cells infected on a logarithmic scale for patient 2-clone B12 ( ) versus wild type or patient 2-clone C7 ( ) (A) or patient 1-clone 2 ( ) versus wild type or patient 1-clone 7 ( ) versus wild type (B). The dotted lines represent the patient clones, and the solid lines represent wild-type virus.
|
|
|
|---|
This assay gives results similar to a traditional growth competition assay when the relative replication efficiencies of single drug-resistant mutants are compared. Using this assay we confirmed studies by Martinez-Picado et al. showing that the PR inhibitor-resistant mutant D30N has reduced replication efficiency compared to L90M (26). Our laboratory has characterized the replication efficiencies of the non-nucleoside RT inhibitor-resistant mutants K103N, V106A, G190S, and P236L. We have consistently observed that the K103N mutant replicates more efficiently than V106A, G190S, and P236L (3, 15, 40). Although our early studies demonstrated a reduction in the replication kinetics of K103N relative to wild-type virus in parallel infections (15), our subsequent studies using growth competition assays have shown minimal reductions in the replication efficiency of K103N, which require 2 to 3 weeks of coculture to detect, using direct sequence analysis of PCR products (20). Another group, using a different cell line and approach to calculating relative fitness, has also found that K103N is more fit than other non-nucleoside RT inhibitor-resistant mutants tested, including G190A and V106A (10). The reduced fitness of G190S has also been noted by other investigators (19). The results produced with our flow cytometry-based assay are thus consistent with these previously published works.
In general, we found similar results for each mutant virus whether it was derived from the pAT1 or the pAT2 vector. We did note that the V106A production rate ratio, p, was significantly reduced compared to K103N when the V106A mutation was present in the Thy1.2 background but only showed a trend to being reduced in the Thy1.1 background. We have confirmed this result with independently generated virus stocks and independently cloned constructs (data not shown), so we do not believe this finding reflects a technical artifact. Although we have not definitively determined the cause for this finding, we postulate that there is a subtle increase in the relative replication efficiency of the AT1 viruses relative to AT2 that was not detected in our 6-day competition assay between the wild-type AT1 and AT2 viruses. If so, this could explain the observed trend toward slightly higher (i.e., less negative) relative fitness values for most of the AT1 viruses (Table 2) and may explain the lack of a significant reduction in the replication efficiency of the V106A-AT1 virus relative to the K103N mutant.
We did not observe significant reductions in the replication efficiency of the L90M and K103N mutants relative to wild-type virus using this assay. This finding is compatible with the relative preservation of replication efficiency of these mutants that has been observed in previous studies. We would not expect a growth competition assay that occurs over a week's time to reliably detect the subtle reductions in replication efficiency that occur with the K103N and L90M mutants; presumably, assaying replication over larger numbers of replication cycles would amplify the differences between the mutant and reference viruses and allow the detection of these subtle replication defects (9). The assay we present here is designed to detect larger reductions in replication efficiency, which we believe are more likely to be clinically significant. We believe that modifying the assay, so that culture supernatant from the day 6 culture can be used to reinfect fresh cells, could improve the sensitivity of this assay, although it would lengthen its turnaround time and increase its complexity.
It should be noted that the flow cytometry-based competition assay described here has some limitations. For example, this assay utilizes cell lines, rather than primary human peripheral blood mononuclear cells (PBMC), for both the production and the propagation of test and reference viruses. An advantage of using cell lines is reduced variability in infection efficiency; such variability can be observed in PBMC obtained from different donors. However, it is possible that the replication in these cell lines may not fully reflect the replication fitness of HIV-1 variants compared to the replication in PBMC. The lamivudine-resistant RT mutant M184V has been shown to have different replication properties in PBMC compared to MT2 cells (5). In contrast, similar relative replication efficiencies were observed in cell lines and PBMC for some non-nucleoside inhibitor- and enfuvirtide-resistant HIV variants (15, 23). Another potential problem that has been observed by others is the rapid reversion of a mutant over the course of a replication assay (37). To date, we have not observed reversion using direct sequence analysis (data not shown), but this could occur, particularly with poorly fit mutants, such as those studied by Sharma et al. (37). In addition, our assay is a recombinant virus assay, reflecting only determinants of replication efficiency that are present in the PR and RT regions of the viral genome. Thus, this assay may give different results than an assay that utilizes whole virus isolates. A whole virus growth competition assay has been shown to correlate with clinical progression in small numbers of patients (32) and may explain the relative prevalence of different HIV variants worldwide (4). Recombinant virus assays may prove to be more predictive of outcome in treatment-experienced, rather than untreated patients, although this remains to be demonstrated.
It is also not clear at present whether multiple cycle or single cycle assays best reflect replication fitness of virus isolates in patients. Much of the preliminary data correlating replication capacity with viral load and CD4 responses has been obtained using a commercially available recombinant-virus, single-cycle assay, in which test and reference strains are compared in parallel infections (7, 8, 13). Although extensive comparisons of single- and multiple-cycle assays have not been performed, it is quite possible that these assays could result in different assessments of relative fitness. Single-cycle assays primarily measure the infectivity of the virus, whereas multiple-cycle assays evaluate infectivity, as well as the rate and efficiency of virus production.
Some data suggest that replication efficiency, measured in cell culture by using different assays, correlates with CD4 count and viral load (6, 11, 13, 17, 38) and clinical prognosis (32). However, there is not yet a consensus on the clinical significance of replication fitness or the best way to measure this parameter. Further studies are thus needed to confirm these observations. Barbour et al. have also shown that the replication capacity conferred by clinical HIV-1 PR and RT sequences, as measured by a recombinant virus single-cycle parallel infection assay, varies dramatically (5 to 81% of the wild-type reference strain NL4-3) (7). It is not known whether differences such as these can influence the natural history of HIV-1 infection or response to therapy. These studies have primarily been performed by using parallel infections, in which the test and reference virus are cultured separately. Growth competition assays to measure replication fitness have not been used in large numbers of clinical samples, primarily because of the labor-intensive nature of these assays. This is unfortunate, because when directly compared, growth competition assays are more sensitive to subtle differences in replication fitness than are parallel cultures (10, 30). An important question is whether stronger correlations between fitness and outcome would be seen if test and reference viruses could be compared in the same culture. We have demonstrated with this assay that patient PR and RT sequences can be cloned into our vectors and tested against a wild-type reference strain to determine their fitness. Since our assay is less labor-intensive, we should be able to test a larger number of patient samples and address some of these questions.
This study was supported by National Institutes of Health grants R01-AI-041387, R01-AI-065217, and N01-AI-38858 (subcontract 200VC007) to L.M.D. and grant R01-AI-049075 to V.P.
Present address: Department of Epidemiology and Biostatistics, College of Public Health, University of South Florida, Tampa, FL 33612. ![]()
|
|
|---|
K reversion in human immunodeficiency virus type 1 reverse transcriptase double mutant K65R+L74V. Virology 321:222-234.[CrossRef][Medline]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»