Transmitted Drug Resistance Mutations in Antiretroviral-Naïve Injection Drug Users with Chronic HIV-1 Infection in Iran

Abstract

The growing incidence and transmission of drug resistant HIV-1 strains due to widespread use of antiretroviral therapy (ART) can jeopardize the success of first-line ART. While there is a known moderate prevalence of transmitted drug resistance (TDR) among newly infected Iranians, no data exist about the rate of these primary resistance mutations among the ART-naïve, chronically infected individuals who are, in fact, the main candidates for ART initiation. To address this issue, we collected blood samples from 40 ART-naïve injection drug-users (IDUs) with chronic HIV-1 infection (seroconversion time ranging from 2 to 9 years) living in Sanandaj, Iran, followed by sequencing of the protease and reverse-transcriptase regions from their HIV-1 genome. Phylogenetic analyses of the sequenced regions revealed that all samples were CRF35_AD. Transmitted resistance mutations were interpreted as surveillance drug-resistant mutations (SDRMs) based on the world health organization (WHO) algorithm. The frequency of SDRMs to any class of antiretroviral drugs was 15%, which included mutations to nucleoside reverse transcriptase inhibitors (NRTIs, 10%), with M41L and M184V as the most common (5%), and non-nucleoside reverse transcriptase inhibitors (NNRTIs, 5%), with K103N as the only detected mutation (5%). Although not in the WHO SDRMs list, several minor protease inhibitor resistant mutations listed in the International Antiviral Society-USA panel were identified, of which M36I, H69K, L89M/V/I (each one 100%) and K20R/T (92.5%) can be considered as polymorphic signatures for CRF35_AD.The relatively high rate of TDR mutations in our study raises concerns about the risk of treatment failure in chronically infected IDUs of Sanandaj city. These results suggest that routine resistance testing should be considered before the therapy initiation in this area. Additional surveillance studies are required to generalize this deduction to other cities of Iran.

About the Authors:

Arash Memarnejadian

* E-mail: [email protected] (AM); [email protected] (SM)

Current address: Department of Microbiology and Immunology, Schulich School of Medicine and Dentistry, The University of Western Ontario, London, Ontario, Canada

Affiliation: Department of Hepatitis and AIDS, Pasteur Institute of Iran, Tehran, Iran

Shahoo Menbari

* E-mail: [email protected] (AM); [email protected] (SM)

Affiliation: Department of Pathology and Medical Laboratory Sciences, Faculty of Paramedicine, Kurdistan University of Medical Sciences, Sanandaj, Iran

Seyed Ali Mansouri

Contributed equally to this work with: Seyed Ali Mansouri, Leila Sadeghi

Affiliation: Department of Microbiology, Faculty of Sciences, Islamic Azad University, Qom, Iran

Leila Sadeghi

Contributed equally to this work with: Seyed Ali Mansouri, Leila Sadeghi

Affiliation: Department of Microbiology, Faculty of Sciences, Islamic Azad University, Karaj, Iran

Rouhollah Vahabpour

Affiliation: Department of Virology, Faculty of Public Health, Tehran University of Medical Sciences, Tehran, Iran

Mohammad Reza Aghasadeghi

Affiliation: Department of Hepatitis and AIDS, Pasteur Institute of Iran, Tehran, Iran

Ehsan Mostafavi

Affiliation: Department of Epidemiology, Pasteur Institute of Iran, Tehran, Iran

Mohammad Abdi

Affiliation: Department of Pathology and Medical Laboratory Sciences, Faculty of Paramedicine, Kurdistan University of Medical Sciences, Sanandaj, Iran

Introduction

Incorporation of highly active antiretroviral therapy (HAART) into clinical practice has resulted in a 60% to 80% decline in rates of Acquired Immunodeficiency Syndrome (AIDS) death and hospitalization [1]. However, an adverse consequence of antiretroviral therapy (ART) is the emergence and selection of antiretroviral resistant mutant variants, the major cause of ART failure in the treatment of AIDS[2–3]. These variants have become widespread, in drug-treated and untreated individuals infected with human immunodeficiency virus (HIV), and have compromised the therapeutic options in drug-naïve infected persons [4]. Transmission of resistant mutants from drug-experienced patients to newly infected drug-naïve individuals was initially noted in developed countries with good access to antiretroviral drugs [5–8]. More recently, ART scaling-up in resource-limited settings is resulting in the occurrence of primary mutations in developing countries, as well [2, 9]. The world health organization (WHO) recommends periodic surveillance of transmitted drug resistance (TDR) mutations in drug-naïve, recently infected individuals in distinct geographical areas[10]. Furthermore, current treatment guidelines recommend routine laboratory testing to assess drug resistance-associated mutations (DRAMs) in patients with acute and chronic infections prior to ART initiation to optimize the treatment regimen [11–13]. These are particularly recommended in countries scaling up ART and in areas where primary resistance has been consistently documented[12, 14]. Lack of testing for baseline resistance, in addition to, other factors including interruption in treatment due to disruption in drug supply, or as a result of financial restrictions and improper administration of drug regimens are the major causes for the occurrence and expansion of drug resistance in developing countries [2].

Since 1986, when the first HIV-1 positive case was reported in Iran, the number of people living with HIV/AIDS has been increased to 27041 cases that have been registered by the end of 2013. In 2015, the number of Iranians with HIV/AIDS is estimated to be more than 120,000[15]. Although an increasing incidence of sexual transmission of HIV has been recently observed, the main route of HIV transmission in Iran is still through injection, so that injection drug users (IDUs) comprise the majority of the HIV-infected population[15–16]. ART in Iran was initiated in2004; however, drug resistance testing is not a prerequisite for ART initiation yet, likely due to limited financial resources and lack of accumulated data on primary resistance. Two studies that used the WHO surveillance drug-resistant mutation(SDRM) list [17] found TDR frequencies of 4.3% [18] and 5.1% [19] among 47 and 39 newly infected ART-naïve Iranian subjects, respectively, with one noting a moderate predicted prevalence (5–15%) throughout the country [18]. Another study using the Stanford drug resistance database[20]noted a TDR mutation frequency of 6.7% among30 ART-naïve subjects living in Tehran city, irrespective of the length of infection [21]. Apart from these, we are aware of no other study addressing the existence of primary drug resistance in Iran and, more importantly, no data is available regarding the prevalence of baseline resistance in ART-naïve Iranians with long-term established (chronic) HIV-1 infection. This lack of information is highlighted in relation to the fact that people with chronic HIV-1 infection are typically candidates for initiation of ART. The objective of this study, therefore, is to determine the frequency of TDR mutations in a chronically infected ART-naïve population in Iran.

Methods

Study subjectsand sample collection

For this cross-sectional study, we enrolled forty HIV-1 positive IDUs who visited Voluntary Counseling and Testing (VCT) centers of Sanandaj, a western city of Iran, in 2011. All participants had been diagnosed as HIV-1 positive before 2009 (seroconversion time ranging from 2 to 9 years), were antiretroviral-naïve and based upon the attending physician’s decision became candidates to initiate HAART at the time of enrollment. Blood (5 mL) was collected from the participants and prepared plasma samples were stored at -70°C. Participants were provided with a detailed questionnaire including socio-demographic questions and blood samples for laboratory results obtained on the same day.

Ethics statement

Study procedure was approved by institutional Ethics Committee of the Pasteur Institute of Iran. All 40 participants signed a written informed consent before blood sample collection.

Determination of drug resistance mutations and HIV-1 subtype

The drug resistance genotypic testing protocol previously described [22]was used with some modifications. Briefly, viral RNA was extracted from 140 μl of plasma by QIA amp viral RNA mini kit (Qiagen, Düsseldorf, Germany), according to the manufacturer’s protocol. WholeHIV-1 protease fragment (PR, nucleotide positions of 2285 to 2550 in HIV-1 HXB2 sequence) and N-terminal fragment of reverse transcriptase (RT, nucleotide positions of 2673 to 3269 in HIV-1 HXB2 sequence)were amplified using specific primers[22]and SuperScript III one-step RT-PCR kit with platinum Taq (Life Technologies, Carlsbad, CA), and further re-amplified in nested PCR using i-Taq Maxime PCR Premix (Intron Biotechnologies, Kyungki-Do, Korea). Final PCR products were purified by QIA quick PCR purification kit (Qiagen, Düsseldorf, Germany) and sequenced using BigDye terminator cycle sequencing kit v3.1and ABI PRISM 3130 Genetic Analyzer (Applied Biosystems, Foster City, CA). Obtained sequences of each sample were assembled to HXB2 PR-RT (GenBank accession number: K03455) as the reference sequence by CLC Main Workbench 5.5 software (CLC bio, Boston, MA). TDR mutations were determined based on the WHO-recommended surveillance drug-resistant mutations (SDRMs) list [17]. To determine the subtypes, PR-RT sequences (1017 bp) were aligned against the subtype and CRF reference sequences retrieved from the Los Alamos HIV-1 database (http://www.hiv.lanl.gov/) by Clustal W. A phylogenetic tree was constructed using the neighbor-joining method with 1000 bootstraps, as implemented in MEGA 5 software [23]. To validate the results, we further submitted PR-RT sequences to the REGA HIV-1 Subtyping Tool V3 [24]. PR-RT sequences obtained in this study (1017 bp) were deposited in the GenBank database and are available under the accession numbers [GenBank: KF544024-KF544063].

Results

The baseline characteristics of the study cases, mean age 35.9 years of age with most cases being male (85%) and single (62.5%) are shown in Table 1. CD4+ T cell counts ranged from 50 to 682 cells/mm3 (median 428), and duration of HIV infection since the first detected seroconversion ranged from 2 to 9 years(median 4.5 years), and 37.5% of the cases were co-infected with HCV. Phylogenetic analysis of PR-RT sequences against the reference sequences of 9 subtypes and 72 CRFs revealed that all 40 samples (100%) were CRF35_AD (the summarized tree is shown in Fig 1). This was also supported by results obtained from analysis by the REGA Subtyping Tool (data not shown).

[Figure omitted. See PDF.]

Table 1. Demographics and laboratory data of study samples.

https://doi.org/10.1371/journal.pone.0126955.t001

[Figure omitted. See PDF.]

Fig 1. Summarized phylogenetic tree of 40 Sanandaj (SDJ) samples based on protease and reverse-transcriptase (PR-RT, 1017 bp) sequences.

Tree was constructed using neighbor-joining method with 1000 replicates. As a summarized representation, only reference sequences from HIV-1 subtypes/CRFs previously reported in Iran including B (●; closed circles), D (□; open squares), A1 (■; closed squares), A2 (○; open circles), CRF01_AE (◆; closed diamonds) and CRF35_AD (▲; closed triangle) are included. GenBank accession numbers for the reference sequences, as well as, bootstrap values over 70% are shown. The scale bar represents nucleotidesubstitutions per site.

https://doi.org/10.1371/journal.pone.0126955.g001

Among the detected DRAMs (Table 2), M41L and M184V that confer the clinical resistance to zidovudine (AZT), lamivudine (3TC), stavudine (d4T) and abacavir (ABC) were the most common to nucleoside reverse transcriptase inhibitors (NRTIs; 5%). K103N associated with resistance to efavirenz (EFV) and nevirapine (NVP) was the only mutation to non-nucleoside reverse transcriptase inhibitors (NNRTIs; 5%). Detection of these mutations is important in view of the fact that AZT, 3TC, EFV and NVP are among the most common first-line regimens in Iran. We found no SDRMs to protease inhibitors (PIs); however, based on the mutations listed in the International Antiviral Society-USA (IAS-USA) panel, 13minor mutations were detected in the PR region, of which M36I, H69K, and L89M/V/I were found in all 40 (100%) andK20R/T was detected in 37 (92.5%) samples. V11I, L33V and D60E that were already reported as minor PI DRAMs in CRF35_AD [18] were not detected in Sanandaj samples. Instead, T74S (2.5%) that is a common accessory mutation and V77I (2.5%) that confers moderate resistance to indinavir (IDV), nelfinavir (NFV) and saquinavir (SQV) were identified in our samples. Based on these data, the overall frequency of WHO-recommended SDRMs in the present population was 15%, which was classified to NRTIs (10%), NNRTIs (5%) and PIs (0%).

[Figure omitted. See PDF.]

Table 2. Frequency of HIV-1 drug resistance-associated mutations among 40 sequenced samples.

https://doi.org/10.1371/journal.pone.0126955.t002

Discussion

In this study we report the HIV-1 circulating subtype and the prevalence of TDR among 40 chronically infected cases in a Western city of Iran. Considering that IDUs comprise the majority of HIV-infected population in Iran [15–16]our samples were selected from a cohort of ART-naïve IDUs. All study cases were determined to be infected with the CRF35_AD variant. This finding is in support of previous studies reporting predominance of CRF35_AD among Iranian HIV-infected IDUs [18–19, 25]. CRF35_AD was originally identified in 2007 in Afghanistan[26] and since then has been reported as a circulating variant in Iran, Afghanistan and Pakistan [27–28].

Genotypic HIV-1 drug resistance can be interpreted using various available databases. Comparison of these algorithms for the interpretation of TDR has shown minor difference for resistance to PIs and NNRTIs, and more pronounced variations for resistance to NRTIs [29]. To address this concern, WHO recommends the use of the SDRMs standard list to make TDR results obtained from different regions, different times and different groups comparable [14]. Accordingly, the overall frequency of SDRMs we found in chronically infected subjects was 15% (NRTIs 10%, NNRTIs 5%, and PIs 0%). Compared to the previously reported TDR frequencies of 5.1% [19]and 4.2% [18] among Iranian newly infected subjects, the frequency of transmitted mutations in our samples appears relatively high although it still remains within the moderate prevalence (5–15%) previously calculated for the SDRMs in Iran [18].

In contrast to these earlier reports that identified a limited number of NRTI mutations (T215D, K219Q [18] and D67G, V75A [19]) with a low overall frequency (4.2% and 5.1%, respectively) in newly infected Iranian cases, we detected a variety of NRTI SDRMs (M41L, D67N, K70R, V75M, F116Y, M184V, L210W, T215Y, and K219E) with a higher overall frequency (10%) in chronically infected IDUs in the city of Sanandaj (Table 2). Furthermore, results from the two mentioned studies found no NNRTI SDRMs, while K103N was identified in 5% of our samples.

We did not find any WHO-recommended SDRMs within the protease region of our samples. However, based on the IAS-USA algorithm, several minor mutations were detected. Moreover, in agreement with previous reports[18–19] we identified high frequencies of M36I, H69K, L89M/V/I and K20R/T mutations that may be considered as polymorphic signatures within the protease region of CRF35_AD. The absence of major mutations to PIs was anticipated and previously reported in newly infected drug-naïve Iranian subjects, as well. The general consensus is that limited access to PIs in Iran has caused an uncommon transmission of PR mutations [18].

It is believed that to facilitate viral replication, initial mutations present in recent infections tend to revert to the wild-type overtime, and hence, in the absence of drug selection, a lower frequency of resistance-associated mutations are expected to be found in subjects with chronic infections compared to those recently infected[30–31]. On the contrary, there are studies reporting no substantial differences between subjects with acute and long-lasting infection [4, 32]. Herein, compared with two previous Iranian studies on recent infections we identified a higher frequency of SDRMs among individuals with chronic infection. Differences in cohort characteristics, as well as, time and region of sampling may have contributed to this discrepancy. Of note, previous studies worked on heterogeneous HIV infected populations from multiple geographical areas, while we focused on IDU samples from a single city. Furthermore, when disregarding the duration of infection, the fact that we selected our samples merely from a cohort of IDUs living in one city may explain the relatively high frequency of mutations we detected. It has been shown that high infection risk behaviors such as sharing syringes and needles, as well as, participating in unprotected sexual practice are frequent among Iranian IDUs [16]. These common risk behaviors in addition to living in one geographical area could result in circulation of mutant viruses in this population. Given this possibility, epidemiological linkage of the study participants may have increased the rate of TDR mutations detected in our study.

In conclusion, our survey shows that the transmission of RT inhibitor resistance occurs in a relatively high frequency in chronically infected HIV-1 positive IDUs in the city of Sanandaj raising the concern of treatment failure in this population. These results are further highlighted when compared to first-line ART in Iran, which since 2008 has been mainly based on the co-administration of NRTIs and NNRTIs. Our findings also emphasize the requirement for routine drug resistance testing to optimize first-line treatment in HIV-1 chronic infections. Our study was geographically limited to one city in the West of Iran. Therefore, in order to generalize these results, prospective drug resistance surveys should be considered in drug-naïve individuals with long-term established HIV-1 infection in other cities throughout Iran.

Acknowledgments

Authors would like to thank Prof. W. Sugiura, and Drs. S. Ibe, J. Hattori and K. Azadmanesh for their valuable supports. They also thank D. Mazzuca for editing the manuscript. AM received “HIV Research Trust Scholarship 2010” to be trained on HIV molecular epidemiology and drug resistance detection in Center for AIDS Research and Education (CARE) in Nagoya Medical Center, Japan.

Author Contributions

Conceived and designed the experiments: AM SM. Performed the experiments: AM SAM LS RV. Analyzed the data: AM SAM LS EM. Contributed reagents/materials/analysis tools: AM SM MRA MA EM. Wrote the paper: AM.

Citation: Memarnejadian A, Menbari S, Mansouri SA, Sadeghi L, Vahabpour R, Aghasadeghi MR, et al. (2015) Transmitted Drug Resistance Mutations in Antiretroviral-Naïve Injection Drug Users with Chronic HIV-1 Infection in Iran. PLoS ONE 10(5): e0126955. https://doi.org/10.1371/journal.pone.0126955

References

1. Palella FJ Jr., Delaney KM, Moorman AC, Loveless MO, Fuhrer J, Satten GA, et al. Declining morbidity and mortality among patients with advanced human immunodeficiency virus infection. HIV Outpatient Study Investigators. N Engl J Med. 1998;338(13):853–60. Epub 1998/03/27. pmid:9516219.

2. Bakhouch K, Oulad-Lahcen A, Bensghir R, Blaghen M, Elfilali KM, Ezzikouri S, et al. The prevalence of resistance-associated mutations to protease and reverse transcriptase inhibitors in treatment-naive (HIV1)-infected individuals in Casablanca, Morocco. J Infect Dev Ctries. 2009;3(5):380–91. Epub 2009/09/18. pmid:19759509.

3. Novak RM, Chen L, MacArthur RD, Baxter JD, Huppler Hullsiek K, Peng G, et al. Prevalence of antiretroviral drug resistance mutations in chronically HIV-infected, treatment-naive patients: implications for routine resistance screening before initiation of antiretroviral therapy. Clin Infect Dis. 2005;40(3):468–74. Epub 2005/01/26. CID33817 [pii] pmid:15668873.

4. Re MC, Monari P, Bon I, Borderi M, Chiodo F. Conflicting interpretations of the prevalence of mutations associated with drug resistance in antiviral naive HIV-1 patients with acute and chronic infection. Int J Antimicrob Agents. 2004;23(2):164–8. Epub 2004/03/12. S0924857903003704 [pii]. pmid:15013042.

5. Re MC, Monari P, Bon I, Gibellini D, Vitone F, Borderi M, et al. Analysis of HIV-1 drug resistant mutations by line probe assay and direct sequencing in a cohort of therapy naive HIV-1 infected Italian patients. BMC Microbiol. 2001;1:30. Epub 2001/12/12. pmid:11737863; PubMed Central PMCID: PMC60646.

6. Sullivan PS, Buskin SE, Turner JH, Cheingsong R, Saekhou A, Kalish ML, et al. Low prevalence of antiretroviral resistance among persons recently infected with human immunodeficiency virus in two US cities. Int J STD AIDS. 2002;13(8):554–8. Epub 2002/08/27. pmid:12194739.

7. Wensing AM, van de Vijver DA, Angarano G, Asjo B, Balotta C, Boeri E, et al. Prevalence of drug-resistant HIV-1 variants in untreated individuals in Europe: implications for clinical management. J Infect Dis. 2005;192(6):958–66. Epub 2005/08/19. JID34030 [pii] pmid:16107947.

8. Manosuthi W, Thongyen S, Nilkamhang S, Manosuthi S, Sungkanuparph S. HIV-1 drug resistance-associated mutations among antiretroviral-naive Thai patients with chronic HIV-1 infection. J Med Virol. 2013;85(2):194–9. Epub 2012/11/20. pmid:23161095.

9. Toni T, Masquelier B, Minga A, Anglaret X, Danel C, Coulibaly A, et al. HIV-1 antiretroviral drug resistance in recently infected patients in Abidjan, Cote d'Ivoire: A 4-year survey, 2002–2006. AIDS Res Hum Retroviruses. 2007;23(9):1155–60. Epub 2007/10/09. pmid:17919113.

10. Bennett DE, Bertagnolio S, Sutherland D, Gilks CF. The World Health Organization's global strategy for prevention and assessment of HIV drug resistance. Antivir Ther. 2008;13 Suppl 2:1–13. Epub 2008/06/26. pmid:18578063.

11. Adolescents PoAGfAa. Guidelines for the use of antiretroviral agents in HIV-1-infected adults and adolescents.: Department of Health and Human Services; 2014.

12. Thompson MA, Aberg JA, Hoy JF, Telenti A, Benson C, Cahn P, et al. Antiretroviral treatment of adult HIV infection: 2012 recommendations of the International Antiviral Society-USA panel. JAMA. 2012;308(4):387–402. Epub 2012/07/24. 1221704 [pii] pmid:22820792.

13. Vandamme AM, Camacho RJ, Ceccherini-Silberstein F, de Luca A, Palmisano L, Paraskevis D, et al. European recommendations for the clinical use of HIV drug resistance testing: 2011 update. AIDS Rev. 2011;13(2):77–108. Epub 2011/05/19. pmid:21587341.

14. Bennett DE, Myatt M, Bertagnolio S, Sutherland D, Gilks CF. Recommendations for surveillance of transmitted HIV drug resistance in countries scaling up antiretroviral treatment. Antivir Ther. 2008;13 Suppl 2:25–36. Epub 2008/06/26. pmid:18575189.

15. National AIDS Committee Secretariat MoHaME. Islamic Re-public of Iran AIDS Progress Report On Monitoring of the United Nations General Assembly Special Session (UNGASS) on HIV and AIDS. March 2014. Report No.

16. Khajehkazemi R, Osooli M, Sajadi L, Karamouzian M, Sedaghat A, Fahimfar N, et al. HIV prevalence and risk behaviours among people who inject drugs in Iran: the 2010 National Surveillance Survey. Sex Transm Infect. 2013.

17. Bennett DE, Camacho RJ, Otelea D, Kuritzkes DR, Fleury H, Kiuchi M, et al. Drug resistance mutations for surveillance of transmitted HIV-1 drug-resistance: 2009 update. PLoS One. 2009;4(3):e4724. Epub 2009/03/07. pmid:19266092; PubMed Central PMCID: PMC2648874.

18. Jahanbakhsh F, Hattori J, Matsuda M, Ibe S, Monavari SH, Memarnejadian A, et al. Prevalence of transmitted HIV drug resistance in Iran between 2010 and 2011. PLoS One. 2013;8(4):e61864. Epub 2013/04/30. PONE-D-12-40726 [pii]. pmid:23626742; PubMed Central PMCID: PMC3633940.

19. Mousavi SM, Hamkar R, Gouya MM, Safaie A, Zahraei SM, Yazdani Z, et al. Surveillance of HIV drug resistance transmission in Iran: experience gained from a pilot study. Arch Virol. 2010;155(3):329–34. Epub 2010/01/07. pmid:20052598.

20. Rhee SY, Gonzales MJ, Kantor R, Betts BJ, Ravela J, Shafer RW. Human immunodeficiency virus reverse transcriptase and protease sequence database. Nucleic Acids Res. 2003;31(1):298–303. Epub 2003/01/10. pmid:12520007; PubMed Central PMCID: PMC165547.

21. Baesi K, Ravanshad M, Ghanbarisafari M, Saberfar E, Seyedalinaghi S, Volk JE. Antiretroviral drug resistance among antiretroviral-naive and treatment experienced patients infected with HIV in Iran. J Med Virol. 2014;86(7):1093–8. Epub 2014/04/18. pmid:24740443.

22. Sugiura W, Matsuda M, Abumi H, Yamada K, Taki M, Ishikawa M, et al. Prevalence of drug resistance-related mutations among HIV-1s in Japan. Jpn J Infect Dis. 1999;52(1):21–2. Epub 2000/05/16. pmid:10808258.

23. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9. Epub 2011/05/07. msr121 [pii] pmid:21546353; PubMed Central PMCID: PMC3203626.

24. Pineda-Pena AC, Faria NR, Imbrechts S, Libin P, Abecasis AB, Deforche K, et al. Automated subtyping of HIV-1 genetic sequences for clinical and surveillance purposes: performance evaluation of the new REGA version 3 and seven other tools. Infect Genet Evol. 2013;19:337–48. Epub 2013/05/11. S1567-1348(13)00181-0 [pii] pmid:23660484.

25. Hamkar R, Mohraz M, Lorestani S, Aghakhani A, Truong HH, McFarland W, et al. Assessing subtype and drug-resistance-associated mutations among antiretroviral-treated HIV-infected patients. AIDS. 2010;24 Suppl 2:S85–91. Epub 2010/07/17. 00002030-201007002-00011 [pii]. pmid:20610954.

26. Sanders-Buell E, Saad MD, Abed AM, Bose M, Todd CS, Strathdee SA, et al. A nascent HIV type 1 epidemic among injecting drug users in Kabul, Afghanistan is dominated by complex AD recombinant strain, CRF35_AD. AIDS Res Hum Retroviruses. 2007;23(6):834–9. Epub 2007/07/03. pmid:17604548.

27. Mumtaz G, Hilmi N, Akala FA, Semini I, Riedner G, Wilson D, et al. HIV-1 molecular epidemiology evidence and transmission patterns in the Middle East and North Africa. Sex Transm Infect. 2011;87(2):101–6. Epub 2010/11/03. sti.2010.043711 [pii] pmid:21036790.

28. Shah S, Xing H, Altaf A, Chen B, Liao L, Jia Y, et al. Antiretroviral drug resistance mutations among treated and treatment-naive patients in Pakistan: diversity of the HIV type 1 pol gene in Pakistan. AIDS Res Hum Retroviruses. 2011;27(12):1277–82. Epub 2011/05/20. pmid:21591988; PubMed Central PMCID: PMC3227241.

29. Liu L, May S, Richman DD, Hecht FM, Markowitz M, Daar ES, et al. Comparison of algorithms that interpret genotypic HIV-1 drug resistance to determine the prevalence of transmitted drug resistance. AIDS. 2008;22(7):835–9. Epub 2008/04/23. 00002030-200804230-00005 [pii]. pmid:18427201; PubMed Central PMCID: PMC2716722.

30. Perno CF, Cozzi-Lepri A, Balotta C, Bertoli A, Violin M, Monno L, et al. Low prevalence of primary mutations associated with drug resistance in antiviral-naive patients at therapy initiation. AIDS. 2002;16(4):619–24. Epub 2002/03/02. pmid:11873006.

31. Salomon H, Wainberg MA, Brenner B, Quan Y, Rouleau D, Cote P, et al. Prevalence of HIV-1 resistant to antiretroviral drugs in 81 individuals newly infected by sexual contact or injecting drug use. Investigators of the Quebec Primary Infection Study. AIDS. 2000;14(2):F17–23. Epub 2000/03/09. pmid:10708278.

32. Kiertiburanakul S, Chaiwarith R, Sirivichayakul S, Ditangco R, Jiamsakul A, Li PC, et al. Comparisons of Primary HIV-1 Drug Resistance between Recent and Chronic HIV-1 Infection within a Sub-Regional Cohort of Asian Patients. PLoS One. 2013;8(6):e62057. Epub 2013/07/05. PONE-D-12-35130 [pii]. pmid:23826076; PubMed Central PMCID: PMC3694952.