Previous Article | Next Article ![]()
Journal of Clinical Microbiology, October 2007, p. 3342-3351, Vol. 45, No. 10
0095-1137/07/$08.00+0 doi:10.1128/JCM.00814-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

International Centre for Reproductive Health, Ghent University, Ghent, Belgium,1 International Centre for Reproductive Health, Mombasa, Kenya,2 AIDS Reference Laboratory, Ghent University, Ghent, Belgium,3 Coast Province General Hospital, Mombasa, Kenya,4 Department of Medical Microbiology, University of Nairobi, Nairobi, Kenya5
Received 17 April 2007/ Returned for modification 10 June 2007/ Accepted 17 July 2007
|
|
|---|
|
|
|---|
In a previous study, we presented a sensitive method for combined detection of the presence of HIV-1 in plasma and subsequent sequencing of the protease (PR) and reverse transcriptase (RT) genes (36). Although this technique allows a significant reduction in the price of follow-up testing, the need for a proper infrastructure for collection and transportation of blood samples remains.
The aim of the current study was to evaluate the possibility of using cellular DNA instead of viral RNA for the detection of drug resistance mutations. As DNA is more stable than RNA, the use of DNA would eliminate the necessity for fast processing of the blood samples after collection and reduce the complexity of the manipulations. The fact that a reverse transcription step is no longer needed for the amplification of the genes of interest also reduces the overall cost of the procedure. Blood spots collected on filter paper have been used for DNA amplification purposes, but reports on the feasibility of performing genotypic analysis for HIV-1 drug resistance remain scarce (24, 33, 43). In resource-poor settings, dried blood spots (DBS) are the preferred method of sample collection in the field, as they do not require electricity or a cold chain. Moreover, once the samples are air dried on the paper, they can be handled as noninfectious and transported to the laboratory without the risk of infection transmission (19).
Whether the same resistance information can be obtained from DNA sequencing and RNA sequencing is still unclear. Results of published reports in this regard are very inconsistent (6, 10, 14, 21, 24, 34). The present study aims to investigate both the possibility of DNA sequencing from whole blood cells or from DBS and the usefulness of DNA sequencing for determination of drug resistance in patients from a resource-poor setting with a high level of HIV-1 subtype diversity.
|
|
|---|
Ten milliliters of EDTA blood was collected, and 50 µl of whole blood was immediately spotted four times on filter paper (Schleicher and Schuell 903). Filter papers were air dried overnight and stored at –20°C until processing. CD4 cell counts were performed (FACScount; Becton Dickinson Immunocytometry, Oxford, United Kingdom), and the remainder of the EDTA blood was centrifuged to collect the plasma and the buffy coat cells. Plasma and whole blood cells were stored at –80°C and –20°C, respectively. All samples were shipped to the AIDS Reference Laboratory at the University Hospital in Ghent, Belgium.
RNA viral load testing. HIV RNA quantification was performed using the ultrasensitive Cobas Amplicor HIV-1 Monitor test, version 1.5 (Roche Molecular Systems, Branchburg, NJ), with a detection limit of 50 copies/ml.
DNA extraction. DNA was extracted from the blood cells by use of the QIAamp DNA blood mini kit (QIAGEN, Venlo, The Netherlands), following the instructions of the manufacturer, and eluted in 50 µl elution buffer.
The procedure to extract DNA from the DBS was adopted from that of Fisher et al. (15). Briefly, one whole DBS was cut into small pieces and transferred to a 2-ml tube. One milliliter of phosphate-buffered saline solution (Cambrex Bio Science, Verviers, Belgium) with 0.1% Tween 20 (Sigma-Aldrich, Bornem, Belgium) was added. The mixture was subjected to a vortex briefly, and after 10 min of incubation at room temperature, the supernatant was removed. This wash step was repeated twice. After the final wash step, 200 µl of 5% Chelex-100 resin solution (Bio-Rad, Nazareth, Belgium) was added and the sample was incubated for 30 min at 56°C and another 30 min at 95°C. The supernatant was then transferred to a clean tube and stored at –20°C until processing.
Amplification of RNA and DNA extracts. Amplification of HIV RNA was performed by nested RT-PCR as previously described (36). DNA extracts from the blood cells and the DBS were amplified by nested PCR using the same primer sets as the ones used for RNA amplification. The first round of amplification was carried out using 5 U/µl of Taq polymerase (Applied Biosystems Inc., Foster City, CA) in a 25-µl reaction volume containing buffer II (Applied Biosystems) with 25 mM MgCl2, 1.25 mM deoxynucleoside triphosphates (Amersham Biosciences, Buckinghamshire, United Kingdom), 0.1% bovine serum albumin (Roche Diagnostics, GmbH Mannheim, Germany), and 0.05 µM each of two forward primers (GAG2 [5'-GAGGAAGCTGCAGAATGGG-3'] and PR1 [5'-ATGATGCAGAGAGGCAATTT-3']) and two reverse primers (RT137 [5'-TTCTGTATGTCATTGACAGTCCAGC-3'] and RT3303 [5'-TAAYTTYTGTATRTCATTGAC-3']). Primers were selected to have a broad specificity for different HIV-1 subtypes.
Cycling conditions were 94°C for 5 min followed by 35 cycles at 94°C for 30 s, 50°C for 30 s, and 72°C for 2 min and a final extension for 7 min at 72°C. Two-microliter volumes of the outer RT-PCR products were used in two separate inner PCRs, one with nested PR primers (PR3 [5'-AGAGCCAACAGCCCCACCA-3'] and PR4 [5'-GGGCCATCCATTCCTGGCTT-3']) and one with nested RT primers (RT1 [5'-CCAAAAGTTAAACAATGGCCATTGACAGA-3'] and RT4 [5'-AGTTCATAACCCATCCAAAG-3']). The second-round PCR was carried out using Taq polymerase (Applied Biosystems Inc., Foster City, CA) in a 50-µl reaction volume containing Taq buffer II with 5 U/µl of Taq polymerase, 25 mM MgCl2, 1.25 mM deoxynucleoside triphosphates (Amersham Biosciences, Buckinghamshire, United Kingdom), and 0.05 µM each of the specific forward and reverse primers. The second-round PCR included an initial denaturation step at 94°C for 5 min followed by 35 cycles at 94°C for 30 s, 50°C for 30 s, and 72°C for 30 s and a final extension for 7 min at 72°C.
In each PCR run, a sample containing 10 copies of HIV-1 DNA (DNA equivalent of 10 8E5 cells) was included as a positive control. Ultrapure water was used as a negative control. Results were accepted only if the results of both controls were correct.
The PR primers amplify a 458-bp fragment spanning the whole PR gene. The RT primers amplify a 646-bp fragment spanning nucleotides corresponding to amino acids 27 to 227 of the RT gene. The amplification products were visualized after electrophoresis in a 1% agarose gel with ethidium bromide.
Genotyping of PR and RT genes. Genotyping was performed using a homemade sequencing assay as described earlier (36). Direct sequencing of both sense and antisense strands was done with the dRhodamine Terminator Cycle Sequencing Ready Reaction kit (Applied Biosystems Inc., Foster City, CA). Sequencing reaction products were analyzed on an ABI310 or ABI3130XL genetic analyzer (Applied Biosystems Inc.). All validations and subsequent manipulations of the sequencing results, as well as the interpretations of the genotyping data and the subtyping, were performed using the Smartgene HIV software packages (Integrated Database Network System; Smartgene, Zug, Switzerland). The selection of drug resistance mutations was based on the Stanford algorithm (http://hivdb.stanford.edu).
Phylogenetic analyses. Phylogenetic analyses were performed using version 3.6 of the PHYLIP package (http://evolution.genetics.washington.edu/phylip.html) with a maximum-likelihood distance matrix and a transition-to-transversion ratio of 2.0. Tree diagrams were plotted with Treeview, version 1.4 (included in the PHYLIP package). Nucleotide differences were calculated using the Emboss infoalign software (http://embossgui.sourceforge.net/demo/infoalign.html).
Statistical analysis. All statistical analyses were performed using SPSS 15.0 (SPSS, IL). Pearson's chi-square test was used to detect possible statistically significant differences between two groups of samples.
|
|
|---|
Sequencing of the viral RNA was attempted on all 29 of these samples and was successful for PR in 28 (96.6%) and for RT in 26 (89.7%). The three patients for whom one or both failed had very low viral loads (73, 74, and 81 copies/ml, respectively). The viral load was undetectable in 31 patients, all treated.
Efficiency of DNA sequencing on whole blood cells and DBS.
Sequencing of the PR and RT genes from whole-blood-cell DNA was successful for 28 (96.6%) and 27 (93.1%), of the 29 blood-cell samples collected from patients with a detectable viral load, respectively, and for 21 (67.7%) and 19 (61.3%) of the 31 samples from patients with an undetectable viral load, respectively. For the DBS, the success rate was slightly lower, with 27 (93.1%) successful sequencing reactions for both PR and RT from the patients with a detectable viral load and 17 (54.8%) PR and 18 (58.1%) RT sequences from the patients with undetectable viral load. No statistically significant difference between the efficiency of genotyping on whole blood DNA and the efficiency of genotyping on DBS DNA was found (P values of 0.2743 for PR and 0.8311 for RT [Pearson's
2 test]).
Subtype distribution. For the 56 patients from whom PR and RT sequences were available, from DNA, RNA, or both, subtyping of the HIV-1 strains was performed by similarity search using the Smartgene subtyping tool. Results showed that the majority of patients was infected with a subtype A virus (n = 31); in order of prevalence, subtype A was followed by subtypes D (n = 7), CRF16_AD (n = 7), C (n = 6), and G (n = 1). Discordances between the subtype determined from the PR sequence and the subtype determined from the RT sequence were seen for four patients (PR/RT subtype combinations D/C, A/C, CRF_AE/C, and C/A). The results of the subtype determination were not influenced by the sample type.
Comparison of the nucleotide sequences from whole-blood DNA, DBS DNA, and plasma RNA. A phylogenetic tree was constructed using 100 PR and RT gene sequences from 38 patients (Fig. 1). The tree revealed a close clustering of the DNA and RNA sequences from the same patients in all cases. The mean nucleotide difference, expressed as the percentage of the total nucleotide length (857 bp), between paired DNA sequences of blood cells and DBS was 2.10% (range, 0.35% to 4.00%). The mean nucleotide difference between the sequences obtained from plasma RNA and blood-cell DNA and the mean nucleotide difference between the sequences obtained from plasma RNA and DBS DNA were 1.82% (range, 0.35% to 3.73%) and 1.90% (range, 0.58% to 3.50%), respectively. Mixtures of nucleotides (International Union of Pure and Applied Chemistry nomenclature) were considered differences. For comparison, the observed mean nucleotide differences between duplicate sequencing reactions on RNA and DNA were 1.45% (range, 0% to 3.37%) and 1.61% (range, 0% to 3.61%), respectively.
![]() View larger version (15K): [in a new window] |
FIG. 1. A phylogenetic tree was constructed based on 100 PR and RT gene sequences from 38 patients. Reference strains for subtypes A, C, D, and CRF16_AD were included in the phylogenetic analysis. The tree was rooted with a subtype B strain. BC, blood cells.
|
The mutation resulting in L33F was detected consistently in all three samples from the same treatment-naïve patient (BV021). The mutation resulting in D30N was detected as part of a mixed population in one of the DNA sequences of patient BV011 (Table 1).
|
View this table: [in a new window] |
TABLE 1. Overview of all patients with mutations at resistance-related positions in the PR gene
|
No LPV-related PI resistance-associated mutations were observed in the two patients on a PI-based regimen (BV114 and BV120) (Table 1).
No NRTI or NNRTI resistance-associated mutations were observed in samples from the 10 treatment-naïve patients. For the 43 of 50 treated patients for whom at least one RT sequence (RNA and/or DNA) was available, NRTI or NNRTI resistance-associated mutations were detected in 13 samples (Table 2). A total of 34 nucleotide substitutions were seen; 10 of these were detected consistently in the RNA and in both DNA sequences, 5 were detected in the RNA and in one of the DNA sequences, 3 were detected only in the DNA from blood cells, and 2 were detected only in the DBS DNA (Table 2).
|
View this table: [in a new window] |
TABLE 2. Overview of all patients with mutations at resistance-related positions in the RT gene
|
The number of resistance mutations detected in the DNA samples of five treated patients with a viral load below 5,000 copies/ml (n = 4) was clearly lower than the number of resistance mutations detected in the DNA samples of six treated patients with a viral load of more than 5,000 copies/ml (n = 17), but the difference did not reach statistical significance (Pearson's chi-square P = 0.11). The numbers of resistance mutations in the RNA in both groups, on the other hand, were comparable, with 12 and 16 mutations found in the RNA from patients with a viral load below 5,000 copies/ml and from patients with a viral load above 5,000 copies/ml (Pearson's chi-square P = 0.83), respectively.
Reproducibility of the DNA sequencing. Replicate amplification and sequencing reactions were performed on the DNA extracts of all individuals listed in Table 2, with the exception of one ART-naïve patient. The results showed the inconsistent detection of some mutations (Table 3). Five resistance mutations that originally were not detected in the DNA but were present in the RNA sample were picked up in at least one of the replicate DNA-sequencing reactions. However, despite replicate analysis, several important resistance mutations remained undetectable even after triplicate sequencing of the DNA: mutations resulting in the presence of 103N and 184V in patient BV141, mutations resulting in the presence of 103N and 184V in patient BV123, mutations resulting in the presence of 67N and 184V in patient BV113, and a mutation resulting in the presence of 65R in patient BV106.
|
View this table: [in a new window] |
TABLE 3. Results of replicate amplification and sequencing reactions performed on DNA extracts
|
Reproducibility of the RNA-sequencing reactions. To evaluate the influence of sampling on RNA sequencing, replicate reactions were performed on four RNA samples (from patients BV113, BV123, BV132, and BV141). Results were much more consistent than those for DNA. In samples from patients BV123 and BV141, the same resistance mutations were detected in all reactions. For patient BV113, the mutation resulting in a D67N substitution was detected in only two of the four replicates. For patient BV132, the mutation resulting in a K103N substitution was detected in only two of the five replicates (results not shown).
|
|
|---|
The results of the study presented confirmed the feasibility of extracting and subsequently sequencing HIV-1 DNA from DBS as shown by others (24, 43). The sensitivity of our in-house genotyping assay for sequencing DNA from DBS was only slightly lower than the sensitivity for sequencing DNA extracted from whole blood cells. Besides allowing the successful sequencing of the PR and RT gene in 93.1% of the patients with a detectable viral load, the DBS also allowed PR and RT gene sequencing of samples from more than 50% of the patients with an undetectable viral load. The reduced sensitivity compared to that of the patients with a detectable viral load suggests a correlation between the viral DNA contributions in DBS and the plasma viral load, as was reported by McNulty et al. (24). Although one might expect that in patients with a low CD4 count, the total number of HIV-infected cells per volume of blood is limited, no relation was seen between low CD4 counts and failure to amplify viral DNA from DBS or failure to detect resistant strains in the DNA (results not shown). The possibility of obtaining sequence information from patients with an undetectable viral load can be valuable for epidemiological studies and subtype determination and might allow the detection of mutations selected during previous suboptimal treatment regimens, as is discussed below.
If all nucleotides are considered, irrespective of the importance of their position for drug resistance, good correlations between the nucleotide sequences of DBS DNA, whole-blood-cell DNA, and RNA were obtained, with differences that only slightly exceeded the differences observed between replicate sequencing reactions starting from the same DNA or RNA extract. For the detection of mutations at positions associated with drug resistance, however, superiority of RNA sequencing over DNA sequencing was observed. Relying on only the DNA sequence would lead to a misinterpretation of the assumed phenotypic resistance in 6 of the 11 patients with significant resistance in the virus in the plasma, irrespective of whether blood cells or DBS sequences were considered.
The feasibility of using viral DNA for the establishment of drug resistance is still controversial. Some studies reported more mutations in the plasma RNA (6, 10, 21, 24), and others reported more mutations in the cellular DNA (14, 34). Most of these discordances may be attributed to differences in the selection of the samples used for the evaluation. While plasma RNA represents the population of short-lived actively replicating virus, viral DNA from infected cells is composed of a heterogeneous mix of DNA from acutely infected, actively virus-producing cells as well as quiescent cells that form the viral reservoir (13, 27, 40). Due to the limited capacities of current population-based sequencing genotyping methods, minor variants can remain undetectable (38). The virus replicating in plasma in the early stage of resistance development most probably derives from a small number of cells in the reservoir, and therefore, detection of resistant virus in plasma will precede detection in the cells. Since all patients in this study initiated HAART very recently and about half of the 11 patients in whom drug resistance mutations were detected still had a viral load of less than 5,000 copies/ml at the time of sampling, they can be considered early treatment failures. Although the number of patients is limited, a correlation between low viral load and a reduced capacity for the detection of resistance mutations in DNA compared to that of resistance mutations in RNA was observed, further strengthening the hypothesis of a relationship between duration of treatment failure and the ability to detect drug resistance in DNA, as was suggested by Bi et al. (6).
Despite a lack of sensitivity in the early phase of treatment failure, sequencing of viral DNA can be valuable, especially for patients with previous episodes of treatment failure that were not documented by resistance analysis. Interruption of a failing regimen or a failing drug will result in an overgrowth of the drug-resistant variant by the more fit wild type (13, 27, 40). In the cellular reservoir on the other hand, the resistant strains will remain detectable for longer time. Venturi et al. (39) demonstrated the kinetics of resistant and wild-type virus in the RNA and DNA compartment very clearly. They found more mutations in the RNA samples than in the DNA samples in the patients on a failing treatment, while those undergoing treatment interruption showed more mutations in the DNA samples. The same higher level of sensitivity of DNA for the detection of drug resistance mutations in a group of patients undergoing treatment interruption was reported by Devereux et al. (14). Previous exposure to EFV explains the mutation resulting in the presence of 190A that was found in the DNA but not in the RNA of patient BV114 one month after a switch from EFV to LPV. All the other patients that were selected for this study were assumed to be on their first line of HAART. Nevertheless, programs aimed at evaluating the influence of a short course of NVP or AZT on mother-to-child HIV-1 transmission have run in this region. Previous exposure to NVP therefore seems the most plausible explanation for the observed 103N mutation in the DNA of patient BV112. This, however, cannot explain the resistance pattern for BV117, a male patient. Surprisingly, both patients BV112 and BV117 showed undetectable viral loads while on NVP-based HAART. The possibility of obtaining an undetectable viral load with an NNRTI-based regimen despite the presence of the mutation resulting in K103N has been observed before (11).
The mutation resulting in the presence of 65R seen in samples from patients BV106 and BV132 is typically known as a tenofovir-associated mutation, but its selection by d4T, especially in patients infected with non-B subtypes, has been reported (37). For patient BV132, the 65R variant was seen in one of the DNA samples only. From the sequencing results for the plasma RNA isolated from this patient, it appears that a mixture of wild-type and mutant strains is present, and this may be the cause for the mutation resulting in the presence of 65R to be missed. The other NRTI or NNRTI resistance-associated mutations that were detected in the DNA sequences only (those resulting in the presence of 70R, 75I, 108I, 219E, and 225H) are considered to have only a minor impact on the drug resistance (http://hivdb.stanford.edu).
The finding of PI-associated mutations in the DNA of two patients who were not on a PI regimen was unexpected. Although the occurrence of 33F in the absence of other PI resistance mutations is rare, we assume that its presence in an isolate from patient BV021 is due to a natural polymorphism. The mutations resulting in the presence of 30N and 46I are seen as part of a mixed population in one of the DNA sequences of patient BV011. We cannot exclude the possibility that this was due to sequencing errors.
The discordances seen between the replicate DNA amplification and sequencing reactions indicate that DNA genotyping, especially genotyping performed on DNA preparations containing small amounts of viral DNA, is prone to a founder effect.
In conclusion, the results of this study demonstrate that whole-blood DNA or DBS DNA sequencing can be a very useful tool for epidemiological studies aimed at analyzing the HIV subtype distribution or the overall distribution of drug resistance mutations in a population. For the follow-up of individual patients, on the other hand, replacement of RNA sequencing by DNA sequencing cannot be recommended, due to the high number of missed mutations, especially in the early phase of treatment failure. The use of DNA and DBS sequencing for the detection of drug resistance mutations selected during previous drug exposure deserves further study.
Kim Steegen is supported by the Flemish Interuniversity Council.
Published ahead of print on 1 August 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»