Drug resistance testing for HIV-1 infected individuals failing antiretroviral therapy (ART) can guide future therapies and improve treatment outcomes. Optimizing individual and population health outcomes in high HIV prevalence but resource-limited settings will ultimately require affordable and accessible drug resistance genotyping and interpretation methods.
HIV-1 drug resistance has the potential to seriously compromise the effectiveness and impact of antiretroviral therapy (ART). As ART programs in sub-Saharan Africa continue to expand, individuals on ART should be closely monitored for the emergence of drug resistance. Surveillance of transmitted drug resistance to track transmission of viral strains already resistant to ART is also critical. Unfortunately, drug resistance testing is still not readily accessible in resource limited settings, because genotyping is expensive and requires sophisticated laboratory and data management infrastructure. An open access genotypic drug resistance monitoring method to manage individuals and assess transmitted drug resistance is described. The method uses free open source software for the interpretation of drug resistance patterns and the generation of individual patient reports. The genotyping protocol has an amplification rate of greater than 95% for plasma samples with a viral load >1,000 HIV-1 RNA copies/ml. The sensitivity decreases significantly for viral loads <1,000 HIV-1 RNA copies/ml. The method described here was validated against a method of HIV-1 drug resistance testing approved by the United States Food and Drug Administration (FDA), the Viroseq genotyping method. Limitations of the method described here include the fact that it is not automated and that it also failed to amplify the circulating recombinant form CRF02_AG from a validation panel of samples, although it amplified subtypes A and B from the same panel.
The HIV epidemic in southern Africa has been evolving rapidly1 with a concomitant exponential increase in individuals on antiretroviral therapy (ART), especially in South Africa2,3. As evidence on the epidemiologic impact of large scale treatment programs in reducing incidence4 and increasing life expectancy in resource-limited settings (RLS)5 continues to accumulate, efforts to increase ART coverage will be intensified. The evolution of guidelines towards the use of treatment as a prevention tool6,7 under the test and treat programs means that the absolute number of individuals on treatment will further increase. Large numbers of individuals will be on ART for longer periods of time as the average life expectancy of individuals on ART nears that of the HIV uninfected population8. The development and transmission of HIV drug resistance has always been considered a threat to the achievements of ART9-12. Thus, there is a need for more rigorous surveillance and monitoring of drug resistance as more individuals are initiated onto ART.
Genotypic drug resistance testing (GRT) has been used successfully in developed countries, both for surveillance as well as monitoring of HIV-1 drug resistance in individuals receiving ART. In these settings, GRT has been integrated into the continuum of care for HIV-1 infected individuals. Most international guidelines recommend GRT for adult or pediatric patients failing ART (first-line and second-line)13-15, pediatric patients exposed to prevention of mother-to-child transmission (pMTCT) regimens but subsequently infected16, and in settings with high-levels of transmitted drug resistance among acutely infected individuals13-15. However, the cost, technology, and infrastructure requirements have limited the implementation of similar approaches to drug resistance monitoring in RLS.
The South African HIV treatment and monitoring guidelines do not currently recommend the use of GRT in guiding choice of ART for individuals failing first-line regimens17. Individuals are switched based primarily on virological (HIV-1 RNA viral load) parameters. However in 2012, the Southern African HIV Clinicians Society published the first Southern African ARV drug resistance testing guidelines18. These guidelines recommend GRT testing for all adults failing first-line and second-line ART and for infected infants and children exposed to pMTCT18. However, GRT is not recommended18 for acutely infected individuals because there is no current evidence for high levels of transmitted drug resistance in southern Africa19-29. It is expected that some of these recommendations will be integrated over time into the national treatment and monitoring guidelines of the various countries in the region. Already, in the 2013 South African treatment guidelines there is now recommendation of GRT at time of second-line failure for adults and at time of first- or second-line PI-based regimen failure for children30.
It has been shown that incorporating GRT into treatment guidelines in South Africa would be potentially cost-neutral. Considering the cost of the second line regimen drugs which are relatively more expensive than the first line drugs, using GRT to identify patients who truly need to be switched to second line therapy will not result any additional cost to the program. In addition, GRT can also identify other reasons for failure, conserve treatment options and generate information about emerging resistance patterns31. Therefore, it is necessary to reduce the cost of drug resistance monitoring methods even further in order to improve access, quality of care and outcomes.
Here, we present a GRT method designed to use generic (open source) primers for reverse transcription, polymerase chain reaction (PCR) and sequencing (Table 1), as well as mostly open source software for drug resistance interpretation. For clinical management, the protocol is complimented by a comprehensive review and reporting method with specialist interpretation of the laboratory drug resistance report with close adherence to the national treatment guidelines. The protocol is divided into four different components; 1) HIV Ribonucleic Acid (RNA) Extraction, 2) Reverse Transcription and Polymerase Chain Reaction (PCR) amplification of viral targets, 3) Sequencing and 4) Bioinformatics methods for analysis of chromatograms, alignment, curation and interpretation of sequence data.
1. Ethylenediaminetetraacetic Acid (EDTA) Whole Blood Processing
Note: Blood can be processed immediately after collection of can be stored at 4 °C for no more than 24 hr.
2. RNA Extraction
3. Reagent Preparation for Reverse Transcription
4. Reverse Transcription
5. Reagent Preparation for PCR
6. Nested PCR
Figure 1. Nested PCR cycling conditions. Click here to view larger image.
7. Gel Electrophoresis
8. PCR Product Cleanup
9. Sequencing Reactions
Figure 3. Scheme representation of a 96-well plate with 12 patient samples being sequenced with 4 primers each (RTC1F, RTC2R, RTC3F, and RTC4R). Click here to view larger image.
Figure 4. PCR cycling conditions for sequencing. Click here to view larger image.
10. Sequencing Cleanup
11. Bioinformatics
12. REGA DB Informatics
The method validated was a modification of a previously reported method20. The Viroseq genotyping method, which has been approved by FDA, was used as the reference method in the validation. A panel of proficiency testing samples obtained from the French National Agencies for Research on AIDS and Viral Hepatitis (ANRS) was used in the primary comparison between the two methods. The two genotyping methods were 100% concordant in identifying all clinically important drug resistance-associated mutations as interpreted by the HIVDB program for the samples that were successfully amplified by both methods. As shown in Table 6, the nucleotide sequences of the three pairs were 99.5% identical. The predicted amino acid sequences were 100% identical. One sample out of five could not be successfully amplified by Viroseq. In addition to the sample not amplified by Viroseq, the in-house method failed to amplify a second sample which was shown to be a circulating recombinant virus (CRF02_AG) by Viroseq. The three samples that amplified with both methodologies were subtype B (two samples) and subtype A (one sample).
Figure 5. Use of a HKY Neighbor Joining tree done as part of sequence quality assurance. There are four pairs/clusters of sequence with very short genetic distances. The genetic distance between RES655 and RES655_1 (same samples sequenced on different days) is 0.003. The is a potential error with the RES637_1/RES638 pair as their genetic distance is too short (0.075 ) for samples from different epidemiologically unlinked individuals. There is another RES637 on the tree with a distance of 0.075 when compared to RES638_1. The CQ01/CQ02 cluster suggests that the two samples from the panel are duplicates of the same sample. They cluster together with the subtype B reference sequence confirming the subtype assigned by the REGA Subtyping tool. CQ05 and CQ04 clustered with subtypes A and G respectively, whereas the REGA subtyping tool classified them as A and CRF02_AG respectively. Another useful tool for HIV subtyping and recombination is SCUEL, which is available at http://www.datamonkey.org. Click here to view larger image.
A panel of five samples was used to assess the precision of the in-house method. Ten replicate genotypes were generated for each of the five samples. Using the 16 Capillary 3130xl genetic analyzer, 48 of the 50 genotypes were generated from 24 runs, prepared on the same day. For all five samples, the predicted amino acid sequences were 100% concordant amongst replicates. For the nucleic acid sequences, there was >99% pairwise similarity.
During the first two years of the use of this method, sixty samples were repeated randomly from RNA extraction to sequencing. There were no statistically significant differences between the sequence quality score and the number of mixed bases between the replicates. Both the nucleotide and amino acid pairwise comparisons for the sixty pairs were greater than 99% identical. Thus the drug resistance mutations for all the pairs were 100% concordant.
Cost reduction
The reaction volumes for RT, PCR and sequencing were reduced by at least half, relative to the original method20,32, without compromising on the quality of the sequences generated. This enabled a reduction in cost of 50% for the RT and PCR stages.
The new method was originally designed to work with six sequencing primer to sequence all 99 codons of the protease gene and the first 300 codons of the reverse transcriptase gene20,32. Similar methods also use six to eight primers33,34. Some recently published methods have used less than six primers, although sometimes sequencing the protease and RT genes seprately35,36. We sought to reduce the number of sequencing primers from six to four, (Figure 6)
Figure 6. Comparison of contiguous sequences from six vs four sequencing primers for the generation of the 1197 bp pol sequence covering all 99 HIV-1 protease codons and the first 300 codons of the reverse transcriptase gene. Click here to view larger image.
Sequences from a set of 17 samples generated from six primers were compared to sequences generated after exclusion of two primers (MAW46 and RTY). The subtypes were 14 subtype C, two subtype B, and one subtype A. There were no significant differences in sequence quality scores. Again, the average pairwise identity between the 17 pairs of nucleic acid was 99% and 100% on the amino acid level. Thus, reducing the sequencing primers from six to four resulted in a reduction in the sequencing cost by almost a third.
The only proprietary software tool used in this protocol was Geneious for sequence assembly. The drug resistance interpretation tools, as well as the report generating tools are all free, open access tools. This reduces the cost further by eliminating the costs associated with the use of proprietary software. Further, collective negotiation allowed the reagents for this protocol to be packaged into a kit for easy access from Life Technologies and is available as the SATuRN/Life Technologies genotyping method 37. Furthermore, SATuRN members can access the reagents at a discounted price.
Clinical Setting
The described protocol has been implemented in the monitoring and surveillance of drug resistance in a rural community in KwaZulu-Natal. A total of 604 genotypes were generated from clinical samples between December 2010 and May 2013 at an amplification rate of 95% for samples with viral loads >1,000 RNA copies/ml. This clinical HIV drug resistance study was approved by the Biomedical Research Ethics Committee of the University of KwaZulu-Natal (ref. BF052/10) and the Health Research Committee of the KwaZulu-Natal Department of Health (ref. HRKM 176/10). Individual patient reports were generated and sent back to the clinics for patient management.
Seventy two (72) genotypes were also generated as part of a surveillance of transmitted drug resistance study, nested within a large prospective population-based HIV surveillance study. The primary samples were needle prick whole blood collected in EDTA microtubes. At genotyping there was an amplification rate of 79%19. Ethics approval for the genotyping of samples from the surveillance study was obtained from the University of KwaZulu-Natal Biomedical Research Ethics Committee (ref. BE066107).
Primer Name | Sequence | Length | Direction | HXB2 Position | |
MAW-26 | TTGGAAATGTGGAAA GGAAGGAC | 23 | Forward | 2028-2050 | 1st round PCR |
RT-21 | CTGTATTTCAGCTATC AAGTCCTTTGATGGG | 31 | Reverse | 3539-3509 | 1st round PCR |
Pro-1 | TAGAGCCAACAGCCC CACCA | 20 | Forward | 2147-2166 | 2nd round PCR |
RT-20 | CTGCCAATTCTAATTC TGCTTC | 22 | Reverse | 3462-3441 | 2nd round PCR |
RTC1F | ACCTACACCTGTCAA CATAATTG | 23 | Forward | 2486-2508 | Sequencing |
RTC2R | TGTCAATGGCCATTG TTTAACCTTTGG | 27 | Reverse | 2630-2604 | Sequencing |
RTC3F | CACCAGGGATTAGAT ATCAATATAATGTGC | 30 | Forward | 2956-2994 | Sequencing |
RTC4R | CTAAATCAGATCCTAC ATACAAGTCATCC | 29 | Reverse | 3129-3101 | Sequencing |
RT-y | GTGTCTCATTGTTTAT ACTAGG | 22 | Reverse | 2967-2946 | Sequencing |
MAW-46 | TCCCTCAGATCACTC TTTGGCAACGAC | 27 | Forward | 2251-2277 | Sequencing |
Table 1. Reverse transcription, PCR, and sequencing custom primers used in the generation of a 1197 bp pol fragment covering all the 99 HIV-1 Protease codons and the first 300 codons of the reverse trascriptase gene.
RT21 (5pmol/ml) | 0.5 | 0.2 |
dNTP (10 mM) | 0.5 | 0.4 |
Total | 1 |
Table 2. dNTP/Primer mix for the reverse transcription reaction.
Reagent | Volume (ml)/reaction | Concentration/reaction |
First Strand Buffer (10x) | 1 | 1 |
MgCl2 (25 mM) | 2 | 4 |
DTT (0.1 M) | 1 | 0.008 |
RNaseOUT (40 U/ml ) | 0.5 | 16 |
Superscript III Reverse Transcriptase (200U/ml ) | 0.5 | 8 |
Total | 5 |
Table 3. Enzyme mix for the reverse transcription reaction.
Reagent | Volume (ml)/reaction | Final Concentration/Reaction |
DEPC treated water | 18.4 | – |
PCR Buffer (10x) | 2.5 | 1 |
MgCI2 (50 mM) | 1 | 2 |
dNTP mix (10 mM) | 0.5 | 0.2 |
Froward primer (5 pmol/ml) | 0.25 | 0.05 |
Reverse primer (5 pmol/ml) | 0.25 | 0.05 |
Platinum Taq Polymerase (5 U/ml) | 0.1 | 0.02 |
Subtotal | 23 | – |
Table 4. Master mix for the nested PCR.
Reagent | Volume (ml)/reaction | Concentration/reaction |
DEPC treated water | 6.1 | |
Sequencing Buffer (5x) | 2 | 1 |
Primer (3.2 pmol/ml) | 0.5 | 0.16 |
Big Dye terminator Sequencing mix | 0.4 | – |
Total | 9 |
Table 5. Master mix for the sequencing reactions.
Viroseq | Inhouse | % NA Similarity | |||||||
Sample I.D | Subtype | Quality score | PR Mutations | RT mutations | Subtype | Quality score | PR mutations | RT Mutations | |
CQ01 | B | 99.9 | M46L, I54L, V82A, L90M | D67N, T69D, K70R, M184V, T215V, K219Q | B | 99.2 | M46L, I54L, V82A, L90M | D67N, T69D, K70R, M184V, T215V, K219Q | 100 |
CQ02 | B | 99.5 | M46L, I54L, V82A, L90M | D67N, T69D, K70R, M184V, T215V, K219Q | B | 99.5 | M46L, I54L, V82A, L90M | D67N, T69D, K70R, M184V, T215V, K219Q | 100 |
CQ03 | NA | NA | NA | NA | NA | NA | |||
CQ04 | CRF02_AG | 98.4 | I54V, V82F, I84V | M41L, L74I, L210W, T215Y, V108I, Y181C | NA | NA | NA | NA | NA |
CQ05 | A | 99.7 | K103N | A | 93 | K103N | 100 |
Table 6. Comparative results from a parallel analysis between the Viroseq genotyping method and the in-house method using a panel of samples provided by the ANRS.
Several low cost in-house methods have been described in efforts to try to make HIV drug resistance genotyping more affordable33,34,36. There is no doubt of the need to integrate drug resistance testing into the continuum of care for individuals on antiretroviral therapy in resource-limited settings. However, most of the reported methods focus on the application of drug resistance genotyping in the surveillance of drug resistance at a population level. The SATuRN/Life Technologies genotyping method is a fully integrated protocol for surveillance and monitoring of drug resistance. This method was designed to be an affordable protocol implementing mostly open source and open access bioinformatics resources for the interpretation of drug resistance and generation of reports for clinical management.
It was shown through comparison with the FDA approved Viroseq genotyping method to be accurate in identifying drug resistance mutations from a panel of ANRS proficiency testing samples, in 100% of laboratory panel samples that were successfully amplified. The accuracy was also assessed on clinical samples of subtype C viruses, the most dominant subtype in southern Africa. The method was as accurate on subtype C samples as it was on subtype A and B. However, if the method would be used in other parts of the world where CRF02_AG is prevalent, there is a need for the modification of the primers since the method failed to amplify one of the panel samples that was shown to have CRF02_AG. Alternatively, a degenerate set of primers sensitive to all group M viruses33,36 could be used in regions where the subtype distribution is more heterogeneous38.
The sensitivity of the reverse transcription and PCR can be increased by extracting RNA from higher volumes of plasma, such as 500 ml. The plasma can be centrifuged at 21,000 x g for 90 min to concentrate the viral particles before proceeding with the protocol as described by the QIAamp viral RNA extraction mini kit.
As shown, the new method has an additional advantage that it produces comprehensive reports for individual patient management. These reports are a consolidation of the genotype, the immunological and virologic monitoring data as well as clinical and treatment history from RegaDB. This is accompanied by a detailed laboratory interpretation of the resistance profile followed by an equally detailed review of the patient’s clinical history as well as treatment recommendations. The use of a specialist physicians to review the reports and provide treatment recommendations for the patients provides the much-needed mentorship for nurse practitioners as well as inexperienced clinicians, who are increasingly providing ART in South Africa as part of the WHO recommendations for task shifting. These clinical reports have been shown to be effective teaching aids for clinicians with little or no experience in drug resistance management. From a patient perspective, our method reduces the need to travel to centralized sites to access specialist HIV services.
Thus, the described protocol taken as a whole provides a good platform through which HIV drug resistance management can be integrated, at an affordable cost, into the continuum of care for HIV infected individuals failing ART. The data generated can be used for epidemiological purposes to assess the evolution and transmission of drug resistance in the community. The size of the pol fragment generated is good enough for more complex phylogenetic analysis which will produce better understanding of the epidemic at population level.
The authors have nothing to disclose.
The authors would like to acknowledge all colleagues who made this work possible, especially Maya Balamane, Elizabeth Johnston White, Sharon Sjoblom, Greg Ording Zakhona Gumede, Xolile Kineri, Phindile Mabaso, Lungisa Ndwandwe, James Garvey, Gavin Cobb, Senzo Maphanga, Terusha Chetty, Kevi Naidoo, Andrew Skingsley, Katharine Stott, and Lungani Ndwandwe. The authors would also like to thank all the personnel of the Department of Health and Africa Centre personnel who work the Hlabisa HIV Treatment and Care Programme.
Superscript III 1st strand Synthesis kit | Life Technologies | 18080051 | Reverse Transcription |
SATURN/LiFE Technologies Custom Primers | Life Technologies | 4473517 | PCR |
Platinum Taq | Life Technologies | 10966026 | PCR |
PureLink QUICK PCR Purification Kit | Life Technologies | K310002 | PCR |
Viroseq | ABBOTT | 4J94-20 | Reverse Transcription and PCR |
Agarose Tablets (Dnase/Rnase free) | BIOLINE | BIO-41027 | PCR |
TBE Buffer | MERCK | 1.06177.2500 | PCR |
O'Range Ruler 200bp DNA Ladder | Fermentas | FE SM0633 | PCR |
Novel Juice | GeneDireX | LD001-1000 | PCR |
MiniBis Bioimaging System | DNR Bioimaging Systems Ltd | Gel Documentation | |
Power Pac 300 | BIORAD | Gel Electrophoresis | |
Big Dye Terminator Kit ver 3.1 | Life Technologies | 4337456 | Sequencing |
Arrays | Life Technolgies | 4319899 | Sequencing |
PoP | Life technologies | 4363785 | Sequencing |
10 x EDTA Buffer | Life Technologies | 402824 | Sequencing |
Formamide | Life Technologies | 4311320 | Sequencing |
5 x Sequencing Buffer | Life Tecgnologies | 4336697 | Sequencing |
3130 xl Genetic Analyzer | Life Technologies | Sequencing | |
GeneAmp PCR System 9700 | Life Technologies | RT/PCR/Sequencing | |
Centrifuge 5804 | EPPENDORF | Sample Processing | |
Centrifuge 5415R | EPPENDORF | RNA Extraction | |
Centrifuge 5415R | EPPENDORF | RT and PCR | |
Centrifuge 5415D | EPPENDORF | PCR Product Clean up | |
Centrifuge 5810 | EPPENDORF | Sequencing Clean up | |
Picofuge | BIORAD | C1301-230V | RT and PCR |
Vortex Genius 3 | IKA | RNA Extraction | |
Vortex Genius 3 | IKA | reagent preparation | |
Vortex mixer | IKA | Sequencing Clean up | |
NanoDrop 2000 UV/VIS spectrophotometer | ThermoScientific | DNA quantification | |
3M Sodium Acetate | MERCK | 567422 | Sequencing Clean up |
Absolute Ethanol | MERCK | SAAR2233540LP | Sequencing Clean up |
1.5ml SARSTEDT Tubes | BIODEX | 72.692.005 | RNA Extraction |
2ml SARSTEDT Tubes | BIODEX | 72.693.005 | RNA Extraction |
2ml Collection tubes | SCIENTIFIC GROUP | MCT-200-NC/S | RNA Extraction |
Optical MicroAmp 96 Well reaction plates | Life Technologies | N8010560 | Sequencing |
200ul 8 Strip StarPCR Tubes with attached flat caps | STAR Lab – supplied by CELTIC | A1402-3700 | RT and PCR |
200ul PCR individual tubes | Scientific Group | CR/3745 | RT and PCR |
Geneious | Biomatters | Sequence analysis | |
Internet Access | Preferrably high speed | ||
Web resources | |||
hivdb.stanford.edu | Stanford University | Drug reistance analysis | |
http://bioafrica.mrc.ac.za:8080/regadb-ui/RegaDB | SATuRN | database | |
http://bioafrica.mrc.ac.za/tools/pppweb.html | SATuRN | Sequence quality tool |