Skip to main content

Comparison of predicted susceptibility between genotype and virtual phenotype HIV drug resistance interpretation systems among treatment-naive HIV-infected patients in Asia: TASER-M cohort analysis



Accurate interpretation of HIV drug resistance (HIVDR) testing is challenging, yet important for patient care. We compared genotyping interpretation, based on the Stanford University HIV Drug Resistance Database (Stanford HIVdb), and virtual phenotyping, based on the Janssen Diagnostics BVBA’s vircoTYPE™ HIV-1, and investigated their level of agreement in antiretroviral (ARV) naive patients in Asia, where non-B subtypes predominate.


Sequences from 1301 ARV-naive patients enrolled in the TREAT Asia Studies to Evaluate Resistance – Monitoring Study (TASER-M) were analysed by both interpreting systems. Interpretations from both Stanford HIVdb and vircoTYPE™ HIV-1 were initially grouped into 2 levels: susceptible and non-susceptible. Discrepancy was defined as a discordant result between the susceptible and non-susceptible interpretations from the two systems for the same ARV. Further analysis was performed when interpretations from both systems were categorised into 3 levels: susceptible, intermediate and resistant; whereby discrepancies could be categorised as major discrepancies and minor discrepancies. Major discrepancy was defined as having a susceptible result from one system and resistant from the other. Minor discrepancy corresponded to having an intermediate interpretation in one system, with a susceptible or resistant result in the other. The level of agreement was analysed using the prevalence adjusted bias adjusted kappa (PABAK).


Overall, the agreement was high, with each ARV being in “almost perfect agreement”, using Landis and Koch’s categorisation. Highest discordance was observed for efavirenz (75/1301, 5.8%), all arising from susceptible Stanford HIVdb versus non-susceptible vircoTYPE™ HIV-1 predictions. Protease Inhibitors had highest level of concordance with PABAKs all above 0.99, followed by Nucleoside Reverse Transcriptase Inhibitors with PABAKs above 0.97 and non-NRTIs with the lowest PABAK of 0.88. The 68/75 patients with discordant efavirenz results harboured the V179D/E mutations compared to 7/1226 with no efavirenz discrepancy (p-value <0.001). In the 3-level comparison, all but one of the discrepancies was minor.


The two systems agreed well with lowest concordance observed for efavirenz. When interpreting HIVDR, especially in non-B subtypes, clinical correlation is crucial, in particular when efavirenz resistance is interpreted based on V179D/E.


In recent years, developing countries have experienced a rapid expansion of antiretroviral drugs (ARVs) in the treatment of HIV-1 infection[1]. Treatment-naive HIV-infected patients may harbour drug-resistance-associated mutations (RAMs) prior to their initial treatment through infection from pre-treated patients. The prevalence of RAMs in treatment-naive patients in resource limited settings was found to vary between different settings[25], mainly below 10%. This implies that drug resistance monitoring prior to treatment may aid the selection of appropriate ARVs.

The aim of the Therapeutics, Research, Education and AIDS Training in Asia (TREAT Asia) Studies to Evaluate Resistance – Monitoring Study (TASER–M) is to monitor HIV-1 drug resistance (HIVDR) and its effects in HIV-infected patients in Asia. The TASER-M cohort consists mainly of ARV-naive Asian individuals infected with HIV-1 AE circulating recombinant form (CRF01_AE)[6]. This is in contrast to developed countries where HIV-1 subtype B predominates. HIV-1 subtypes differ in their gene sequences and this may influence susceptibility to ARVs[7, 8]. Some non-B subtypes have natural polymorphisms at resistance associated positions[9], leading to possible misinterpretation of drug susceptibility.

There are several software-based drug resistance interpretation systems currently available, which are used in patient care to design subsequent regimens. The two systems adopted in TASER-M are (i): the Stanford University HIV Drug Resistance Database (Stanford HIVdb) which has a publicly available tool that provides genotypic analysis and interpretation via web-based sequence submission. The Stanford HIVdb uses a rules based system in which all mutation scores are added to derive a level of predicted viral resistance to each ARV. These rules are pre-determined based on research findings of published studies[10, 11]; (ii) Janssen Diagnostics BVBA’s vircoTYPE™ HIV-1 which is a proprietary software that provides a predicted or ‘virtual’ phenotype of a sequence, based on a large genotype-phenotype database. The tool provides a calculated fold change (FC) predicted phenotype of the patient’s genotype, compared to a reference sequence, derived from linear regression modelling[12, 13]. The Stanford HIVdb and vircoTYPE™ HIV-1 utilise consensus B and HXB2, respectively, as the reference strains[10, 11, 13, 14].

Because different interpretation algorithms use different rules to predict drug susceptibility, results may differ. Studies have shown that differences do exist with varying degree of discordances[1521], and the proportion of discordances for each ARV could be subtype dependent[22, 23].

The objective of this study was to determine the yet unstudied level of agreement between genotype and virtual phenotype from the Stanford HIVdb and vircoTYPE™ HIV-1 resistance interpretation in the TASER-M cohort, where non-B subtypes predominate.


TASER-M recruitment began in 2007, with a total of 12 participating sites from Thailand, Hong Kong, Malaysia, Philippines and Indonesia. The TASER-M cohort comprises of 95% treatment naive patients at time of enrolment with pre-treatment sequences contributing to more than 90% of total available sequences. Ethics approvals were obtained from the University of New South Wales Ethics Committee and institutional review boards at the participating clinical sites. Informed consent was obtained from participants prior to enrolment. ARV treatment-naive patients were required to have viral load and resistance testing performed prior to initiation of ARVs and at yearly intervals post initiation[24]. Genotyping was performed in local TAQAS[25] laboratories.

Sequences were included in this analysis if they were pre-treatment sequences derived from HIV-infected patients enrolled as treatment-naive in the TASER-M cohort. FASTA sequence files were analysed by both the Stanford HIVdb and vircoTYPE™ HIV-1. FASTA files were submitted to the Stanford HIVdb (Version 6.0.11) through Sierra – The Stanford HIV Web Service (Version 1.0) for genotypic interpretation. Subtyping was analysed using REGA HIV-1 Subtyping Tool - Version 2.0.

For virtual phenotyping, FASTA files were submitted to Janssen Diagnostics. The vircoTYPE™ HIV-1 (VPT 4.3.01) results were then transferred to Advanced Biological Laboratories (ABL) TherapyEdge 3.9.2 web based software, and retrieved directly from TherapyEdge 3.9.2. The interpretation of the FC is made in reference to the clinical cut-offs (CCOs) or biological cut-offs (BCOs)[26, 27]. If the CCOs are used, the resistance interpretations are categorised as maximal response if the FC is less than the lower CCO (CCO1), reduced response if the FC is between CCO1 and upper CCO (CCO2), and minimal response if the FC is greater than CCO2. For ARVs utilising the BCO, the interpretations are either susceptible if the FC is below the BCO, or resistant if the FC is above the BCO[14]. For lamivudine and etravirine, a combination of both BCO and CCO are used in vircoTYPE™ HIV-1 predictions. The lower CCO (CCO1) is replaced by the BCO, resulting in two cut-off points: BCO and CCO2. If the FC value is below the BCO, the interpretation will be susceptible. If the FC is between the BCO and CCO2, the result will be classified as reduced response, whilst minimal response is the resulting interpretation when FC is greater than CCO2 for these three drugs.

Two level interpretations

A total of 15 ARV interpretations (7 Nucleoside Reverse Transcriptase Inhibitors (NRTIs), 3 non-NRTIs (NNRTIs) and 5 Protease Inhibitors (PIs)) and were compared. The Stanford HIVdb provides 5 levels of susceptibility prediction for each ARV (susceptible, potential low-level resistance, low-level resistance, intermediate resistance and high-level resistance). The vircoTYPE™ HIV-1 utilises either 2 or 3-level resistance prediction for any given ARV, depending on the cut-offs used.

To be able to compare both systems effectively, the 5-level interpretations from the Stanford HIVdb needed to be combined in order to match the levels of interpretation provided by vircoTYPE™ HIV-1. For comparison of all ARVs, 2-level comparisons were performed to enable the inclusion of ARVs where virtual phenotyping only provided predictions at 2 levels. That is, the 3-level interpretations available from the Stanford’s Sierra database (susceptible, intermediate and resistant), initially grouped from the 5-level interpretation, were combined further into susceptible (susceptible) and non-susceptible (intermediate and resistant). The vircoTYPE™ HIV-1’s susceptible and maximal response were categorised as susceptible whilst reduced response, minimal response and resistant were grouped as non-susceptible. Discrepancy was defined as a discordant drug resistant interpretation between the Stanford HIVdb and the vircoTYPE™ HIV-1 report, using the 2-level interpretations, i.e. susceptible versus non-susceptible, and vice versa.

The ARV that exhibited the highest discordances was chosen for further investigation into the mutation pattern and predicted FC values associated with the discrepancies. The mutations were analysed using the Stanford’s HIV Drug Resistance Mutations by Drug Class (November 6, 2009)[28] (Stanford HIVDR list) only, as vircoTYPE™ HIV-1’s list of relevant mutations was not available. Differences in FC values were analysed by comparing the median FC between the discordant and non-discordant groups.

Three level interpretations

To evaluate whether the discrepancies were major or minor in nature, a sub comparison was also made using 3-level susceptibility predictions. In this comparison however, efavirenz, nevirapine and emtricitabine were excluded as they did not have 3 levels of virtual phenotype interpretations available. As a consequence, the available prediction output for all other ARVs were susceptible, maximal response, reduced response and minimal response. The susceptible category here belongs to lamivudine and etravirine as a result of the BCO being used as the lower cut off point, instead of the CCO1, thus replacing “maximal response” with “susceptible” interpretation. The vircoTYPE™ HIV-1’s susceptible and maximal response were then categorised as susceptible, while reduced response was intermediate, and minimal response was defined as resistant. The interpretations were then compared to the 3-level Stanford HIVdb results. Major discrepancy was defined as having a susceptible result from one system and resistant from the other. Minor discrepancy corresponded to having an intermediate interpretation in one system, with a susceptible or resistant result in the other.

Statistical analysis

To take into account agreement due to chance, the prevalence adjusted bias adjusted kappa coefficient (PABAK) was used to analyse the level of agreement between the two interpretation systems[29]. The interpretation of PABAK was made based on Landis and Koch’s[30] categorisation: 0 to 0.20 slight agreement; 0.21 to 0.40 fair agreement; 0.41 to 0.60 moderate agreement; 0.61 to 0.80 substantial agreement; and 0.81 to 1.00 almost perfect agreement.

Median values were analysed using Wilcoxon rank sum test. Categorical data was analysed using Chi-square, Fisher’s exact test or McNemar’s test.

Analyses were performed using SAS software version 9.2 (SAS Institute Inc., Cary, NC, USA) and STATA software version 10.1 (STATA Corp., College Station, TX, USA). Ordinal scale PABAK (PABAK-OS) for 3-level interpretation was analysed using PABAK-OS web-based tool[31].



A total of 1301 eligible naive patients consented between 2007 and September 2010 were included. Out of the 1301 baseline sequence files analysed, most were from male patients (66%) with a median age of 36 years (IQR 31–44). The main ethnicity was Thai (71%) followed by Chinese (20%). Heterosexual contact was the predominant mode of exposure (72%). The median pre-treatment CD4 and RNA counts were 100 cells/μL (IQR 34–199) and 100,000 copies/mL (IQR 41,345 – 230,000), respectively. The main subtype was CRF01_AE (78%). Subtype B accounted for 15% and other non-B subtypes were 7%.

Two level interpretations

The number of discrepant results, using the 2-level interpretation method, between Stanford HIVdb and the vircoTYPE™ HIV-1 for each ARV are illustrated in Figure1. The number of discrepancies ranged from 0 to 75. The ARV with the highest number of discordances was efavirenz (75/1301, 5.8%), followed by zidovudine (21/1301, 1.6%) and etravirine and abacavir (18/1301, 1.4% each). All of the 75 discrepancies for efavirenz arose from susceptible interpretation from the Stanford HIVdb versus non-susceptible interpretation from vircoTYPE™ HIV-1. The observed efavirenz agreement was 94% with 1199 patients being susceptible and 27 non-susceptible. For zidovudine, the observed agreement was 98% with 1274 susceptible and 6 non susceptible patients. Emtricitabine and darunavir/r both had 100% agreement with no discrepancy, however 1293/1301 agreements for emtricitabine were from susceptible prediction and 8/1301 were non-susceptible. For darunavir/r, 100% of the agreements were susceptible from both systems.

The PABAK values (Table1) ranged from 0.88 to 1.00. The lowest PABAK (0.88) was observed for efavirenz. This is consistent with efavirenz having the lowest number of concordances. The PABAK for etravirine, abacavir, didanosine and zidovudine was second lowest at 0.97. For the two drugs with no discrepancy, the corresponding PABAK for both drugs was 1.00 with 95% confidence interval (95%CI) of (1.00-1.00). Overall, PIs had highest level of concordance with individual PABAK values all above 0.99, followed by NRTIs with PABAKs above 0.97 and NNRTIs with lowest PABAK of 0.88. Using Landis and Koch’s categorisation, the agreement between the two interpreting systems for all ARVs considered could be classified as “almost perfect agreement” since the PABAK values all lie within the 0.81-1.00 interval.

Figure 1

Number and percentage of discordant drug resistance interpretations for each ARV – 2-level interpretation.

Table 1 PABAK values for each ARV – 2-level interpretation

As efavirenz had significantly higher number of discordances compared to other ARVs (p-value<0.001 for all comparisons) and the lowest PABAK, it was chosen for further investigation. The relevant mutations were extracted from the Stanford HIVdb and compared against the Stanford HIVDR list only. 68/75 (91%) patients with discrepant efavirenz results harboured the V179D/E mutations (V179D (83%), V179E (17%)) compared to only 7/1226 (0.6%) of patients with no efavirenz discrepancy (p-value <0.001). Apart from V179D/E, these 68 patients did not present any other efavirenz RAMs. The other 7/75 (9%) patients who had the discrepancy did not harbour any efavirenz-associated mutations. These 7 discrepant patients may harbour other RAMs that were found to significantly affect susceptibility to efavirenz by the vircoTYPE™ HIV-1’s algorithm, which were not present in the Stanford HIVDR list. Of note, T39K, K43S, K122E, D123S, D177E, Q207A and R211S were the most common RT-mutations (86%) found in these 7 patients with respect to the reference wildtype virus. However no common RT-mutation was identified across all 7 patients. Out of the 75 patients with efavirenz discrepancy, 79% were non-B subtype compared to 85% in patients showing concordance (p=0.135). Among the 1226 patients with no discrepancy, the efavirenz-associated mutations were A98G (0.3%), K103N (0.8%), K103S (0.2%), V106M (0.2%), V179D (0.6%), Y181C (0.6%), G190A (0.3%), G190E (0.1%), P225H (0.2%), F227C (0.1%) and K238T (0.1%). RAMs found in discrepant patients for other NNRTIs were: etravirine (18 discordances): V179D (38.9%), V179E (11.1%), G190A (5.6%) and G190E (5.6%); nevirapine (12 discordances) V179D (75%); V179E (8.3%) and F227C (8.3%). Although G190A/E and F227 were found in a proportion of these patients, a similar pattern to efavirenz could be seen whereby discrepant patients with V179D/E did not harbour any additional NNRTI RAMs.

In vircoTYPE™ HIV-1, the BCO was used to determine the level of susceptibility of efavirenz. Therefore only one cut-off point (BCO = 3.3) and two interpretations were available. There was no discrepancy found in patients with FC less than the assigned BCO of 3.3 for efavirenz. That is, all patients with the discrepancies had reported FC greater than 3.3, which resulted in a non-susceptible or resistant virtual phenotype interpretation.

In light of 100% efavirenz discrepant patients possessing non-susceptible virtual phenotype results, median FC values were compared between the discordant and concordant group for those with non-susceptible vircoTYPE™ HIV-1 interpretation only. This is because by restricting the comparison to those that resulted in non-susceptible interpretation (FC>3.3), it was possible to make direct comparison without diluting the results with the low FC values whereby no discrepancy was found. There was a total of 102 patients with FC above the BCO of 3.3 for efavirenz. Out of these 102 patients, 27 had no discrepant result, with median FC of 31.90 and inter quartile range (IQR) (7.10-87.10). The median FC for those with discrepancies was 5.10 (p-value <0.001; IQR 4.30-6.10). This suggests that for patients with non-susceptible virtual phenotype results for efavirenz, the group with the discordances had a median FC that is lower and nearer to the cut-off point than those without discordances.

Three level interpretations

Interpretations were regrouped into susceptible, intermediate and resistant in order to determine the number of minor and major discrepancies. Efavirenz, emtricitabine and nevirapine were excluded. Figure2 and Table2 indicate that all of the discordances, but one, were classified as minor discrepancies. The proportion of exact agreement ranged from 98% to 100%. As with the two level interpretations, darunavir/r had 100% agreement with no discrepancy.

Figure 2

Number and percentage of major and minor discrepancies for each ARV– 3-level interpretation, excluding efavirenz, nevirapine and emtricitabine.

Table 2 PABAK – OS – 3-level interpretation

The results of Table2 suggest that almost all of the discrepancies that exist between Stanford HIVdb and vircoTYPE™ HIV-1 were minor in nature. That is, one system produced an intermediate interpretation while the other produced either a susceptible or resistant interpretation. Etravirine was the only ARV that possessed a major discrepancy where the Stanford HIVdb’s interpretation was susceptible and vircoTYPE™ HIV-1 was resistant. Without the inclusion of efavirenz, the number of minor discrepancies in Table2 ranged from 0 to 21. By comparing the 2 and 3-level discrepancies, it was found that the total number of discrepancies from the 3-level comparison was slightly greater than that from the 2-level comparison for a number of ARVs. For example, abacavir had 20 and 18 discrepancies in the 3 and 2-level comparisons, respectively. This is due to the fact that the differences between intermediate and resistant interpretations would not have been captured in the 2-level categorisation, as both interpretations were grouped as non-susceptible. However, the total discrepancies for each ARV between the 2 and 3-level groupings differed by a maximum of 2, and the PABAK and PABAK-OS agreement interpretations remained unchanged, therefore the simplification of the interpretations into two categories did not greatly alter the discordance outcome.


The overall results suggest that both the Stanford HIVdb and the vircoTYPE™ HIV-1 provide comparable drug resistance interpretations. In the 2-level comparison, efavirenz showed the highest number of disagreements, all of which were Stanford HIVdb susceptible - vircoTYPE™ HIV-1 resistant pairs. The PABAK values for all ARVs considered remained in the “almost perfect agreement” category. The results from the 2-level comparison suggest that the lowest level of agreement belongs to the NNRTI drug class. This is in contrast to previous studies[1519] where NRTIs were found to have lowest concordances. Of note, the majority of the patients in this cohort harboured CRF01_AE subtype virus, whereas B subtype was predominant in preceding studies[16, 17, 19]. Studies comparing non-B sequences reported highest discordances in either PI or NRTI drug class[20, 21].

The RAMs associated with efavirenz discrepancy were V179D/E. These are likely to be natural polymorphisms in HIV-1 RT. V179D/E are often found in treatment-naive individuals infected with non-B subtypes[6, 32, 33]. V179D is associated with low-level resistance to NNRTIs and is considered to have no significant effect on the efficacy of NNRTI–containing regimen. The 5-level Stanford HIVdb predicted “potential low-level resistance” in all sequences with discordant efavirenz interpretation containing V179D/E. A combination of V106I and V179D has been shown to confer significant resistance to efavirenz[34]. Among TASER-M patients, only 1 had a combination of V106I and V179D. V179D also confers significant resistance to NNRTIs when presented together with K103R[35]. Only 2 patients had a combination of K103R and V179D.

Since the mutations were extracted and compared against the Stanford HIVDR list only, it should be noted that the RAMs found to be associated with the discrepancies are those based on the Stanford HIVDR list alone. This is one limitation of the study. The likely presence of additional mutations in vircoTYPE™ HIV-1 coupled with the differences in mutation weights or scores are the most likely causes of the discrepancies. Therefore, the significance of V179D/E found in this study should not be interpreted as the cause of the differences between the two interpreting systems, but merely reflects the mutation pattern based on the Stanford HIVDR list. However, it has been reported that natural polymorphisms found in non-B subtype have contributed to the majority of the discordances, although this was mainly associated with the protease region of the sequence[20].

The 3-level comparison showed that most of the discrepancies found were minor in nature. These may have less clinical importance than the major discordances. The availability of only two virtual phenotype predictions for efavirenz, nevirapine and emtricitabine, another limitation of this study, did now allow for further break down of discordant groups. As such, discrepancies between intermediate and resistant categorisation for these three drugs could not be captured. Nevertheless, the minor discrepancies illustrated in this study suggest that the results did not differ greatly from one system to the other.

The differences in median FC between those with and without efavirenz discrepancies suggest that one should take into consideration the actual predicted FC values in addition to the resistance interpretation when evaluating efavirenz resistance. Although discrepancies only occurred in sequences with resistant interpretation, those who were resistant with lower FC tended to result in discordant interpretation compared to those with high FC. Additionally, by simplifying resistance calls into categories, the true value of predicted phenotype is diminished. In the Stanford HIVdb, the 3 or 5-level classification should also be interpreted together with the resistance scores in order to take into account the variation of susceptibility within each classification. It would be of interest to re investigate when the CCOs become available for efavirenz, nevirapine and emtricitabine. It is likely that the FC values associated with efavirenz discrepancies would lie in the middle region between CCO1 and CCO2, resulting in minor discrepancies.


This study has shown high concordance between the Stanford HIVdb and vircoTYPE™ HIV-1 in naive patients. The differences in the proportions of discordances from different patient groups reported in previous studies[20, 22, 23] and the higher efavirenz discrepancies compared to all other ARVs analysed indicate that more data is required to enable the algorithms to adapt appropriately to different sequence variations. When interpreting HIVDR, especially in non-B subtypes, careful mutation-specific, interpretation of genotyping and virtual phenotyping results are recommended to identify potential discordances.

Members of the TASER study include

· PCK Li† and MP Lee, Queen Elizabeth Hospital and KH Wong, Integrated Treatment Centre, Hong Kong, China;

· N Kumarasamy and S Saghayam, YRG Centre for AIDS Research and Education, Chennai, India;

· S Pujari and K Joshi, Institute of Infectious Diseases, Pune, India;

· TP Merati‡ and F Yuliana, Faculty of Medicine, Udayana University & Sanglah Hospital, Bali, Indonesia;

· A Kamarulzaman and LY Ong, University Malaya Medical Center, Kuala Lumpur, Malaysia;

· C KC Lee and B HL Sim, Hospital Sungai Buloh, Sungai Buloh, Malaysia;

· M Mustafa and N Nordin, Hospital Raja Perempuan Zainab II, Kota Bharu, Malaysia;

· R Ditangco and RO Bantique, Research Institute for Tropical Medicine, Manila, Philippines;

· YMA Chen, YJ Chen and YT Lin, Taipei Veterans General Hospital and AIDS Prevention and Research Centre, National Yang-Ming University, Taipei, Taiwan;

· P Kantipong and P Kambua, Chiang Rai Regional Hospital, Chiang Rai, Thailand;

· P Phanuphak and S Sirivichayakul, HIV-NAT/Thai Red Cross AIDS Research Centre, Bangkok, Thailand;

· W Ratanasuwan and R Sriondee, Faculty of Medicine, Siriraj Hospital, Mahidol University, Bangkok, Thailand;

· T Sirisanthana and J Praparattanapan, Research Institute for Health Sciences, Chiang Mai University, Chiang Mai, Thailand;

· S Sungkanuparph, S Kiertiburanakul, and L Chumla, Faculty of Medicine, Ramathibodi Hospital, Mahidol University, Bangkok, Thailand;

· R Kantor, Brown University, Rhode Island, U.S.A.;

· AH Sohn, N Durier and T Singtoroj, TREAT Asia, amfAR -- The Foundation for AIDS Research, Bangkok, Thailand;

· DA Cooper, MG Law, and A Jiamsakul, The Kirby Institute, University of New South Wales, Sydney, Australia.

† Steering Committee Chair, ‡ Co-Chair

The content of this publication is solely the responsibility of the authors and does not necessarily represent the official views of any of the institutions mentioned above.


  1. 1.

    Global HIV/AIDS response: WHO Library Cataloguing-in-Publication Data. epidemic update and health sector progress towards universal access: progress report 2011. 2011

    Google Scholar 

  2. 2.

    Liao L, Xing H, Shang H, Li J, Zhong P, Kang L, Cheng H, Si X, Jiang S, Li X: The prevalence of transmitted antiretroviral drug resistance in treatment-naive HIV-infected individuals in China. J Acquir Immune Defic Syndr. 2010, 53 (Suppl 1): S10-14.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  3. 3.

    Dean J, Ta Thi TH, Dunford L, Carr MJ, Nguyen LT, Coughlan S, Connell J, Nguyen HT, Hall WW, Nguyen Thi LA: Prevalence of HIV type 1 antiretroviral drug resistance mutations in Vietnam: a multicenter study. AIDS Res Hum Retroviruses. 2011, 27 (7): 797-801. 10.1089/aid.2011.0013.

    PubMed  CAS  Article  Google Scholar 

  4. 4.

    Hamers RL, Siwale M, Wallis CL, Labib M, van Hasselt R, Stevens WS, Schuurman R, Wensing AM, Van Vugt M, Rinke de Wit TF: HIV-1 drug resistance mutations are present in six percent of persons initiating antiretroviral therapy in Lusaka, Zambia. J Acquir Immune Defic Syndr. 2010, 55 (1): 95-101. 10.1097/QAI.0b013e3181e544e0.

    PubMed  CAS  Article  Google Scholar 

  5. 5.

    Apisarnthanarak A, Jirayasethpong T, Sa-nguansilp C, Thongprapai H, Kittihanukul C, Kamudamas A, Tungsathapornpong A, Mundy LM: Antiretroviral drug resistance among antiretroviral-naive persons with recent HIV infection in Thailand. HIV Med. 2008, 9 (5): 322-325. 10.1111/j.1468-1293.2008.00562.x.

    PubMed  CAS  Article  Google Scholar 

  6. 6.

    Sungkanuparph S, Oyomopito R, Sirivichayakul S, Sirisanthana T, Li PC, Kantipong P, Lee CK, Kamarulzaman A, Messerschmidt L, Law MG: HIV-1 drug resistance mutations among antiretroviral-naive HIV-1-infected patients in Asia: results from the TREAT Asia Studies to Evaluate Resistance-Monitoring Study. Clin Infect Dis. 2011, 52 (8): 1053-1057. 10.1093/cid/cir107.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  7. 7.

    Holguin A, Ramirez de Arellano E, Rivas P, Soriano V: Efficacy of antiretroviral therapy in individuals infected with HIV-1 non-B subtypes. AIDS Rev. 2006, 8 (2): 98-107.

    PubMed  Google Scholar 

  8. 8.

    Kantor R: Impact of HIV-1 pol diversity on drug resistance and its clinical implications. Curr Opin Infect Dis. 2006, 19 (6): 594-606. 10.1097/QCO.0b013e3280109122.

    PubMed  CAS  Article  Google Scholar 

  9. 9.

    Kantor R, Katzenstein DA, Efron B, Carvalho AP, Wynhoven B, Cane P, Clarke J, Sirivichayakul S, Soares MA, Snoeck J: Impact of HIV-1 subtype and antiretroviral therapy on protease and reverse transcriptase genotype: results of a global collaboration. PLoS Med. 2005, 2 (4): e112-10.1371/journal.pmed.0020112.

    PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Liu TF, Shafer RW: Web resources for HIV type 1 genotypic-resistance test interpretation. Clin Infect Dis. 2006, 42 (11): 1608-1618. 10.1086/503914.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  11. 11.

    Shafer RW: Rationale and uses of a public HIV drug-resistance database. J Infect Dis. 2006, 194 (Suppl 1): S51-58.

    PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Vermeiren H, Van Craenenbroeck E, Alen P, Bacheler L, Picchio G, Lecocq P: Prediction of HIV-1 drug susceptibility phenotype from the viral genotype using linear regression modeling. J Virol Methods. 2007, 145 (1): 47-55. 10.1016/j.jviromet.2007.05.009.

    PubMed  CAS  Article  Google Scholar 

  13. 13.

    Pattery T, Verlinden Y, De Wolf H, Nauwelaers D, Van Baelen K, Van Houtte M, Mc Kenna P, Villacian J: Development and performance of conventional HIV-1 phenotyping (Antivirogram(R)) and genotype-based calculated phenotyping assay (virco(R)TYPE HIV-1) on protease and reverse transcriptase genes to evaluate drug resistance. Intervirology. 2012, 55 (2): 138-146. 10.1159/000332013.

    PubMed  CAS  Article  Google Scholar 

  14. 14.

    Steegen K, Bronze M, Van Craenenbroeck E, Winters B, Van der Borght K, Wallis CL, Stevens W, de Wit TF R, Stuyver LJ: A comparative analysis of HIV drug resistance interpretation based on short reverse transcriptase sequences versus full sequences. AIDS Res Ther. 2010, 7: 38-10.1186/1742-6405-7-38.

    PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Muñoz M, Carmona R, Pérez-Alvarez L, Cilla G, Suarez MD, Delgado E, Contreras G, Corral J, de Goicoetxea MJL, Medrano L: Analysis of discrepancies in the interpretation of antiretroviral drug resistance results in HIV-1 infected patients of Basque Country, Spain. Journal of Clinical Virology. 2005, 33 (3): 224-229. 10.1016/j.jcv.2004.06.014.

    PubMed  Article  Google Scholar 

  16. 16.

    Sturmer M, Doerr HW, Staszewski S, Preiser W: Comparison of nine resistance interpretation systems for HIV-1 genotyping. Antivir Ther. 2003, 8 (3): 239-244.

    PubMed  Google Scholar 

  17. 17.

    Puchhammer-Stockl E, Steininger C, Geringer E, Heinz FX: Comparison of virtual phenotype and HIV-SEQ program (Stanford) interpretation for predicting drug resistance of HIV strains. HIV Med. 2002, 3 (3): 200-206. 10.1046/j.1468-1293.2002.00116.x.

    PubMed  CAS  Article  Google Scholar 

  18. 18.

    Zazzi M, Romano L, Venturi G, Shafer RW, Reid C, Dal Bello F, Parolin C, Palu G, Valensin PE: Comparative evaluation of three computerized algorithms for prediction of antiretroviral susceptibility from HIV type 1 genotype. J Antimicrob Chemother. 2004, 53 (2): 356-360. 10.1093/jac/dkh021.

    PubMed  CAS  Article  Google Scholar 

  19. 19.

    Liu L, May S, Richman DD, Hecht FM, Markowitz M, Daar ES, Routy JP, Margolick JB, Collier AC, Woelk CH: Comparison of algorithms that interpret genotypic HIV-1 drug resistance to determine the prevalence of transmitted drug resistance. AIDS. 2008, 22 (7): 835-839. 10.1097/QAD.0b013e3282f5ff71.

    PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Vergne L, Snoeck J, Aghokeng A, Maes B, Valea D, Delaporte E, Vandamme AM, Peeters M, Van Laethem K: Genotypic drug resistance interpretation algorithms display high levels of discordance when applied to non-B strains from HIV-1 naive and treated patients. FEMS Immunol Med Microbiol. 2006, 46 (1): 53-62. 10.1111/j.1574-695X.2005.00011.x.

    PubMed  CAS  Article  Google Scholar 

  21. 21.

    Poonpiriya V, Sungkanuparph S, Leechanachai P, Pasomsub E, Watitpun C, Chunhakan S, Chantratita W: A study of seven rule-based algorithms for the interpretation of HIV-1 genotypic resistance data in Thailand. J Virol Methods. 2008, 151 (1): 79-86. 10.1016/j.jviromet.2008.03.017.

    PubMed  CAS  Article  Google Scholar 

  22. 22.

    Snoeck J, Kantor R, Shafer RW, Van Laethem K, Deforche K, Carvalho AP, Wynhoven B, Soares MA, Cane P, Clarke J: Discordances between interpretation algorithms for genotypic resistance to protease and reverse transcriptase inhibitors of human immunodeficiency virus are subtype dependent. Antimicrob Agents Chemother. 2006, 50 (2): 694-701. 10.1128/AAC.50.2.694-701.2006.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  23. 23.

    Yebra G, de Mulder M, del Romero J, Rodriguez C, Holguin A: HIV-1 non-B subtypes: High transmitted NNRTI-resistance in Spain and impaired genotypic resistance interpretation due to variability. Antiviral Res. 2010, 85 (2): 409-417. 10.1016/j.antiviral.2009.11.010.

    PubMed  CAS  Article  Google Scholar 

  24. 24.

    Hamers RL, Oyomopito R, Kityo C, Phanuphak P, Siwale M, Sungkanuparph S, Conradie F, Kumarasamy N, Botes ME, Sirisanthana T: Cohort Profile: The PharmAccess African (PASER-M) and the TREAT Asia (TASER-M) Monitoring Studies to Evaluate Resistance--HIV drug resistance in sub-Saharan Africa and the Asia-Pacific. Int J Epidemiol. 2011

    Google Scholar 

  25. 25.

    Land S, Cunningham P, Zhou J, Frost K, Katzenstein D, Kantor R, Chen YM, Oka S, DeLong A, Sayer D: TREAT Asia Quality Assessment Scheme (TAQAS) to standardize the outcome of HIV genotypic resistance testing in a group of Asian laboratories. J Virol Methods. 2009, 159 (2): 185-193. 10.1016/j.jviromet.2009.03.016.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  26. 26.

    Winters B, Van Craenenbroeck E, Van der Borght K, Lecocq P, Villacian J, Bacheler L: Clinical cut-offs for HIV-1 phenotypic resistance estimates: update based on recent pivotal clinical trial data and a revised approach to viral mixtures. J Virol Methods. 2009, 162 (1–2): 101-108.

    PubMed  CAS  Article  Google Scholar 

  27. 27.

    Winters B, Montaner J, Harrigan PR, Gazzard B, Pozniak A, Miller MD, Emery S, van Leth F, Robinson P, Baxter JD: Determination of clinically relevant cutoffs for HIV-1 phenotypic resistance estimates through a combined analysis of clinical trial and cohort data. J Acquir Immune Defic Syndr. 2008, 48 (1): 26-34. 10.1097/QAI.0b013e31816d9bf4.

    PubMed  Article  Google Scholar 

  28. 28.

    HIV Drug Resistance Mutations by Drug Class: 2009, [],

  29. 29.

    Byrt T, Bishop J, Carlin JB: Bias, prevalence and kappa. J Clin Epidemiol. 1993, 46 (5): 423-429. 10.1016/0895-4356(93)90018-V.

    PubMed  CAS  Article  Google Scholar 

  30. 30.

    Landis JR, Koch GG: The measurement of observer agreement for categorical data. Biometrics. 1977, 33 (1): 159-174. 10.2307/2529310.

    PubMed  CAS  Article  Google Scholar 

  31. 31.

    Single Case Research: web based calculators for SCR analysis. (Version 1.0). College Station, TX: Texas A&M University, [],

  32. 32.

    Sagoe KW, Dwidar M, Lartey M, Boamah I, Agyei AA, Hayford AA, Mingle JA, Arens MQ: Variability of the human immunodeficiency virus type 1 polymerase gene from treatment naive patients in Accra, Ghana. J Clin Virol. 2007, 40 (2): 163-167. 10.1016/j.jcv.2007.07.016.

    PubMed  CAS  Article  Google Scholar 

  33. 33.

    Oumar AA, Jnaoui K, Kabamba-Mukadi B, Yombi JC, Vandercam B, Goubau P, Ruelle J: Genotypic evaluation of etravirine sensitivity of clinical human immunodeficiency virus type 1 (HIV-1) isolates carrying resistance mutations to nevirapine and efavirenz. Acta Clin Belg. 2010, 65 (4): 242-244.

    PubMed  CAS  Article  Google Scholar 

  34. 34.

    Gatanaga H, Ode H, Hachiya A, Hayashida T, Sato H, Oka S: Combination of V106I and V179D polymorphic mutations in human immunodeficiency virus type 1 reverse transcriptase confers resistance to efavirenz and nevirapine but not etravirine. Antimicrob Agents Chemother. 2010, 54 (4): 1596-1602. 10.1128/AAC.01480-09.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  35. 35.

    Parkin NT, Gupta S, Chappey C, Petropoulos CJ: The K101P and K103R/V179D mutations in human immunodeficiency virus type 1 reverse transcriptase confer resistance to nonnucleoside reverse transcriptase inhibitors. Antimicrob Agents Chemother. 2006, 50 (1): 351-354. 10.1128/AAC.50.1.351-354.2006.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

Download references


The TREAT Asia Studies to Evaluate Resistance is an initiative of TREAT Asia, a program of amfAR, The Foundation for AIDS Research, with major support provided by the Dutch Ministry of Foreign Affairs through a partnership with Stichting Aids Fonds, and with additional support from amfAR and the National Institute of Allergy and Infectious Diseases (NIAID) of the U.S. National Institutes of Health (NIH) and the National Cancer Institute (NCI) as part of the International Epidemiologic Databases to Evaluate AIDS (IeDEA) (grant no. U01AI069907). Queen Elizabeth Hospital and the Integrated Treatment Centre are supported by the Hong Kong Council for AIDS Trust Fund. The Kirby Institute is funded by the Australian Government Department of Health and Ageing, and is affiliated with the Faculty of Medicine, The University of New South Wales. The use of vircoTYPE™ HIV-1 is supported by Janssen Diagnostics BVBA.

Author information



Corresponding author

Correspondence to Awachana Jiamsakul.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AJ performed data manipulation, analysis plan, statistical analysis, interpretation of results and drafted the manuscript. SS originated the study concept, participated in the analysis plan, data collection, helped interpret results and edited the manuscript. RK participated in the analysis plan, helped interpret results and edited the manuscript. PL, SS, TS, PK, CL, AK, WR, RD, TS performed data collection and revised the final manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Jiamsakul, A., Kantor, R., Li, P.C. et al. Comparison of predicted susceptibility between genotype and virtual phenotype HIV drug resistance interpretation systems among treatment-naive HIV-infected patients in Asia: TASER-M cohort analysis. BMC Res Notes 5, 582 (2012).

Download citation


  • Asia
  • HIV
  • Resistance
  • Interpretation
  • Algorithm