Skip to main content

High prevalence of HIV-1 transmitted drug resistance among therapy-naïve Burmese entering travelers at Dehong ports in Yunnan, China



The overall success of Human immunodeficiency virus type 1 (HIV-1) antiretroviral therapy (ART) was heavily challenged upon the occurrence of drug resistance. Dehong Prefecture witnessed not only the first report of HIV-1 infection but also the experimental adoption of antiviral treatment in China. The transmission and epidemic of HIV-1 in Dehong is impacted by cross-border activities. The characteristics of HIV-1 drug resistance among therapy-naïve Burmese entering travelers in Yunnan and their speculated origin are still not clarified.


Two hundred ninety-eight HIV-1 infected Burmese entering travelers at Dehong ports were recruited between 2003 and 2012. The partial HIV-1 pol gene fragments were amplified and sequenced for the analysis of drug-resistance mutations (DRMs). Phylogenetic analysis on gag-pol gene was conducted to elucidate phylogenetic and evolutionary characteristics of these drug resistant strains.


It was figured out that the occurrence ratio of HIV-1 drug resistance among HIV-1 infected entering travelers from Myanmar was up to 12.8%. The resistant mutations covered several types, including one type of PI mutations (L33F), six types of NRTI mutations and seven types of NNRTI. Close genetic relationship was observed in the phylogenetic analysis on gag-pol gene among the drug resistant strains respectively from Dehong, other Yunnan areas, neighboring provinces (Guangxi) and neighboring countries (Thailand and Myanmar).


The findings in this study revealed that HIV drug resistant locus is spreading from the population who is receiving drug-resistance treatment to the new infectors, which indicates the urgency of surveillance work on drug resistance among the migrant population with high risks of HIV infection.

Peer Review reports


Antiretroviral therapy (ART) has greatly contributed to reducing HIV-1 viral load, prolonging patients’ lives, lowered the AIDS morbidity and mortality, and also prevented HIV-1 transmission [1,2,3]. However, HIV-1 drug resistance occured in a substantial proportion of treated patients and accumulated over time with the long-term use of ART [4]. Drug resistance has become the magnificent challenge for current antiretroviral therapy, especially in developing countries.

In the wake of implementation of ART, resistant strains carrying drug-resistant mutations have been increasingly prevalent [5]. The new infected individuals have greater chance to be infected directly by existing drug resistant strains instead of the natural generation of strains during the antiretroviral therapy procedure. The transmitted drug resistance (TDR), also known as acquired drug resistance, makes strains get more advantages in circuiting due to its insensitive response to various antiviral drugs. Hence, TDR has been a major difficulty and an intensive focus in AIDS treatment and HIV surveillance.

Yunnan province has been worst hit by HIV-1 in China [6]. The multi-channel introduction from neighboring southeast Asian countries, coexistence of multiple subtypes and continuous occurrence of new recombinants all lead to the most complex genetic characteristics of HIV-1 in Yunnan. As the pioneering ART implemented region, with a long history of ART use, Dehong Prefecture, undoubtedly, had a larger proportion of resistant strains among all epidemic ones in Yunnan. It has been well documented that the cross-border travelers were proved to be a key factor of the high infection rate and complex distribution of HIV in Yunnan [7, 8]. HIV-1 cross-border transmission and exchanges of various circulating subtypes were considered to be accelerated by this floating population. Moreover, migrations of travelers carrying resistant HIV-1 strains also bring great challenges to TDR surveillance work and AIDS treatment.

Located in the most southwestern Yunnan Province, Dehong prefecture has a shared border line of 503.8 km with Myanmar and is considered as the first station of illegal drug traffic from golden triangle region into China. Due to the special geographic location, the first endemic of HIV-1 in China was identified in Dehong, with the most HIV-1 prevalence and the complex distribution of HIV-1 subtypes or recombinants [9,10,11]. With the implementation of reform and opening-up policy in China, the cross-border travel and business became more frequent. More than 23 million travelers entered Yunnan province from Myanmar via Dehong port in 2017 (

Besides the continuous cross-border transmission of HIV-1 between China and Myanmar through travelers’ immigration, the prevalence of TDR HIV-1 strains among this huge cross-border population should be investigated too. Although it is obvious that HIV-1 infection rate of this population is much higher than other general groups based on our previous reports [8], the frequency of occurrence and types of transmitted drug resistance mutations (TDRMs) are still unidentified in this population. In current study, the plasma samples were collected randomly from entering travelers at Dehong port between 2003 and 2012. HIV-1 infection was determined and TDR mutations were further analyzed by sequencing on entire pol gene of HIV-1. The occurrence, types and molecular characteristics of these TDR strains will be further explored.



From 2003 to 2012, one of every five-hundred Myanmar travelers entering China at land ports of Dehong prefecture, Yunnan, was randomly sampled for HIV-1 test. The processes and methods of HIV-1 detection and laboratory confirmation have been described in our previous reports [8]. Totally, 22,699 travelers were required for HIV test, and 1163 were identified as HIV-1 infectors. Among them, 298 travelers who never received any antivirus therapy before (ART-naïve), were selected in this study. Those cases fall intoeach year irregularly, ranging from 27 (year 2011) to 56 (year 2008). The collected plasma samples were stored at − 80 °C until RNA extraction. Data of the social-demographic characteristics and risk behaviors for HIV infection were obtained from the enrolled participants via face-to-face interviews with trained medical personnel by administering a structured questionnaire.

HIV-1 pol gene amplification and sequencing

HIV-1 RNA was extracted from plasma using the High Pure Viral RNA kit (Qiagen, Valencia, CA, USA). The partial HIV-1 pol gene fragment (HXB2:2147–3462) was amplified using the One-Step reverse transcription PCR kits (TaKaRa, Dalian, China) with primers of MAW26 (5′-TTGGAAATGTGGAAAGGAAGGAC-3′) and RT21 (5′-CTGTATTTCTGCTATTAAGTCTTTTGATGGG-3′). RT-PCR procedure was as follows: HIV-1 RNA denaturation at 65 °C for 30 s, addition of the reaction mixtures at 4 °C, incubation at 50 °C for 30 min, 94 °C for 2 min, then 35 cycles of 94 °C for 30 s, 55 °C for 30 s, and 72 °C for 2 min 30 s. Nested PCR was performed using 2× Taq PCR MasterMix (TaKaRa) with primers PRO-1 (5′-CAGAGCCAACAGCCCCACCA-3′) and RT20 (5′-CTGCCAGTTCTAGCTCTGCTTC-3′). PCR conditions were: 94 °C for 5 min, then 35 cycles of 94 °C for 30 s, 63 °C for 30 s, and then 72 °C for 2 min 30 s. Positive PCR products were visualized in agarose gel electrophoresis, and commercially sequenced by Biomed company (Beijing, China).

Genotyping and drug-resistance mutations (DRMs) analysis

The obtained sequence of pol gene was assembled, edited and aligned as previously described [12, 13], and then submitted to NCBI GenBank to get accession numbers as MG787428-MG787459. The Stanford University HIV Drug Resistance Database HIVdb Program (version 8.4, was used to investigate the presence of drug resistance mutations conferring resistance to Protease Inhibitors (PIs), Nucleoside Reverse Transcriptase Inhibitors (NRTIs), and Non-Nucleoside Reverse Transcriptase Inhibitors (NNRTIs) [14, 15]. Genotypes were determined in reconstructed maximum likelihood (M-L) tree using MEGA 7.0 software in the Kimura 2 parameter model with gamma distribution and invariant sites.

Phylogenetic and recombination analyses of HIV drug resistant strains

To elucidate phylogenetic characteristics of these drug resistant strains, gag-pol gene sequences were obtained by using reverse transcription nested PCR as described before [12]. Phylogenetic analyses were performed using the M-L (Maximum likelihood) tree method based on the Kimura 2-parameter model with 1000 bootstrap replicates in MEGA version 7.0 [16]. The HIV-1 recombinant structure was screened using the Recombinant Identification Program available from the HIV database ( The suspected novel recombinant strains were subjected to boot scanning and informative-site analysis using SimPlot version 3.5.1 with boot-scan window sizes of 200 bases, a step size of 20 bases and 100 replicates [17].


Demographic characteristics of recruited patients

Between 2003 and 2012, a total of 298 Burmese ART-naïve travelers with HIV-1 infection, who entered China at Dehong port in Yunnan province, China, were recruited in this study. Most of their age varieded from 21 to 40 (accounted for 76%) at their sampling time and 26 (10.44%) were female. Their occupations were long-distance drivers, businessman, migrant rural workers, hair salon workers, and tourists. Long-distance drivers, who transported goods shuttling between China and Myanmar, accounted for the majority in this population (Table 1). Unsafe sexual contact (12.80%) and intravenous drug use (33.20%) were major informed risk factors for HIV-1 infection; however, 135 infected travelers did not know how they got infected. Pol gene sequences (1.3 kb) were successfully obtained from 250 travelers. Phylogenetic analysis showed that CRF01_AE was the predominant HIV-1 subtype/CRF (142 cases, 56.8%), followed by subtype C (49 cases, 19.60%), subtype B (31 cases, 12.40%), CRF08_BC (9 cases, 3.60%), CRF07_BC (5 cases, 2.00%). 14 cases identified other CRFs, including CRF01B, CRF64_BC, and other types of B/C recombination (Additional file 1: Figure S1).

Table 1 Demographic characteristics of all recruited subjects for drug-resistance mutations detection

Prevalence of HIV drug resistant mutations

According to HIV Drug Resistance Database, 45 mutations associated with drug resistance were defined in 32 therapy-naïve Burmese travelers. Statistically, the rate of TDRs in this population was 12.80% (32/250). The annual TDRs rate was distributed irregularly and ranged from 0 (2005) to 21.88% (2004). Of 32 HIV-1 infected travelers with TDRs mutation, CRF_01AE (20 cases) was identified as predominant subtype/CRF, followed by subtype C (4 cases), CRF08_BC (2 cases), subtype B (2 cases) and other undetermined subtypes (4 cases). The TDRs rate was significantly different from the subtypes/CRF. Although 9 cases were infected with CRF08_BC in this population, 2 of them (22.22%) were proved to carry drug resistant sites. There was no TDRs observed with CRF07_BC.

In detail, mutations associated with resistance to NNRTI, NRTI, or PI were found in 29 (64.44%), 11 (24.44%) and 5 (11.11%) cases respectively (Fig. 1). For NNRTI resistance-associated mutations, V179D/T was identified as the major mutation (20 cases), and the other 6 types of mutation were E138A, A98G, K101P, K103 N, Y181H and Y188H. NRTI resistance-associated mutations found in this population included M184 V (n = 3, 6.67%), M41 L (n = 2, 4.44%), T215F (n = 2, 4.44%), L74F/V (n = 2, 4.44%), D67G (n = 1, 2.22%) and K70Q (n = 1, 2.22%) (Figure1). In addition, there was only one PI resistance-associated mutation found (L33F) in 5 cases. Furthermore, four (12.50%) subjects displaying DRMs had dual-class drug resistance mutations, including 2 cases to NRTI and NNRTI and other 2 cases to NNRTI and PI. No obvious annual change of DRMs occurrence was found from 2003 to 2012, and also no different DRMs distribution was observed in various subtypes/CRF (Fig. 2).

Fig. 1
figure 1

Forty-five mutations associated with transmitted resistance to any drug, defined according to the World Health Organization

Fig. 2
figure 2

Trend of DRMs rate of Therapy-naïve Burmese entering travelers at Dehong port in the past decade (a) and the TDR frequency in different subtype (b)

Phylogenetic characteristics of the drug resistant strains

To elucidate phylogenetic characteristics of these drug resistant strains, sequences of gag-pol gene were used for phylogenetic analysis. As a result, 7 subtypes/CRFs of HIV-1 were found among these travelers infected with drug resistant strains (Fig. 3). All the used reference CRF08_BC strains were identified in China. The strains of subtype B and C in the current study were clustered together with the strains from China and neighboring countries, Thailand and Myanmar. The CRF01_AE strains, the most prevalent subtype among this population, showed close genetic relationship with reference strains from Southeast Asian countries (Thailand and Vietnam) and China. Two independent Chinese clusters were observed in the M-L tree, including 7 strains from this study. A Drug resistant strain, TDR-12018, was classified as CRF64_BC, which was firstly identified in IDUs (Injection drug users) from Dehong. The other two strains, TDR-06017 and TDR-08278, formed a tight cluster with the reference CRF01_B strains from Guangxi and CRF0107 from Jiangsu Province. However, TDR-08269 strain couldn’t be classified into any subtype/CRF. Further restructuring analysis showed it was inter-subtype of CRF01_AE and subtype C (Fig. 3) with small fragment gene of C inserted in CRF01_AE gag-pol gene.

Fig. 3
figure 3

Maximum likelihood phylogenetic tree of the drug resistant strains in this study based on the sequence of HIV-1 gag-pol fragments and the recombinant structures of TDR-08269. The strains labeled with black circles were collected in this study. Strains that could not be precisely classified as known subtypes, CRFs, or URFs are shown in dashed red squares in the neighbor-joining tree


Yunnan Province has been the key midpoint of traffic and business between south Asian countries and China. Dehong has been an important port in the business with Myanmar. The numerous cross-border population not only boosted rapid economic development but worked as carriers for HIV transmission. Thus, they were considered a ‘bridge population’ who may accelerate the exchange of virus genes between Yunnan and Myanmar. In our previous research, we have confirmed that HIV infection ratio (5.12%) of entering travelers at Dehong port is higher than any other ports in Yunnan [8]. This population, shaped by the specific characteristics of migration and HIV-1 infection susceptibility, has become a critical medium in the spreading of drug resistant strains. Therefore, it is of great necessity to investigate the prevalence of HIV infection and drug resistance strains among travelers from Myanmar.

In our previous study, a total of 1163 HIV infected Burmese travelers entering Yunnan at Dehong port have been detected. To address the prevalence of drug-resistance mutations among this bridge population, 298 travelers who never received any anti-virus therapy before, were randomly selected for this study. The young long-distance drivers and businessmen dominated in this population. Although the infection route of 173 cases is still unidentified, face-to-face interviewing has guaranteed the authenticity of obtained epidemiologic information and privacy of participants. Among the informed infection route, drug injection is still the main risk behavior due to the common syringe sharing. Phylogenetically, more than 6 subtypes/CRFs have been discovered, including the main CRF01_AE (142 cases, 56.8%), following by subtype C, subtype B. The HIV-1 subtype distribution is similar to that of Yunnan Province in previous reports [18, 19].

With the implementation of active ART, the occurrence of drug resistant mutation becomes more frequent, which results in more resistant strains transmitting among new HIV-1 infectors and therapy-naïve infectors [20]. A global study on TDR frequency shows the highest (12.9%) TDR rate is in North America and the lowest (4.2%) in Asia [21]. It has been documented that the prevalence of drug-resistance was 45.1% in the ART-failure individuals and 2.1% in the controlled therapy-naïve individuals in Yunnan province [22]. However, the prevalence of drug-resistance is still controversial in the controlled therapy-naïve individuals due to small recruited participants (2/94). In current study, the TDR rate is 12.8% among the therapy-naïve Burmese travelers entering Dehong port in Yunnan Province, which is far higher than the earlier reported level. The high prevalence of TDR mutation was also discovered in other regions of China, like Henan Province (14.3%) [23] and Beijing (11.96%) [24]. More than two decades of ART treatment in these areas could explain the higher TDR prevalence in Yunnan, Beijing and Henan, etc. However, the prevalence of TDR among the bridge population is obviously higher than any other groups in Yunnan, even compared with homosexual group in Kunming (4.6%) [25], and general population in Dehong (< 5%) [26].

One type of PI mutations (L33F), six types of NRTI mutations and seven types of NNRTI mutations have been identified in this population. The major TDRMs observed in this study have been reported in a previous study of HIV Drug Resistance among ART-failure individuals in Yunnan Province, which figured out mutations such as M184 V/I, K103 N, V106A, Y181C and G190A were common [22]. Notably, the similar prevalence of TDRMs has also been found among therapy-naïve HIV patients in neighboring countries Vietnam [20] and Thailand [27]. It suggested a close correlation between the TDRMs from this population and the drug resistant mutants generated during the course of ART. Furthermore, M184 V is the most frequently detected NRTI resistance mutation, which was also found in previous studies [28, 29]. L33F and V179D/T are the types of mutation that confer potential low-level resistance. Even considering the previous report of occurrence in Hebei and Beijing, the mutation rate of V179D/T up to 6.92% is still an astonishing level [30, 31]. Despite of the previous report stating the possible relationship between TDRMs and HIV-1 subtype [32], there is little evidence shown in this study supporting this statement, possibly due to the relatively smaller sampled group in the investigation. However, DRMs analysis was not carried out on the mutations of integrase inhibitors in this study, which were also significant class of resistant mutations.

Phylogenetic analysis on gag-pol gene of all resistant strains has also been performed to identify their possible origins in this study. The results demonstrate that the Dehong epidemic strains share a very close genetic relationship with the strains from Yunnan, neighboring provinces (Guangxi) and countries (Thailand, Myanmar, etc). The high prevalence of TDR strains in these regions is reasonable considering the long and earliest history of ART implementation. The floating population, whose inherent mobility determines the multi-regional sources of the epidemic resistant strains, is considered as the highly possible TDR strains carrier in further spreading to broader range. Thus, it is critical and fundamental to monitor their drug-resistance proceeding to prevent them from further transmission of drug resistant strains.


This study focuses on the occurrence of TDR and its phylogenetic characteristics in Burmese travelers. The prevalence of TDR among this group is as high as to 12.8%, far in excess of general population and reported level in China. The high similarities of drug resistant mutations between this population and ART failure individuals, imply a very close genetic distance between our resistant strains and epidemic strains from Yunnan and its surrounding areas. Based on these similarities, we speculate that HIV-1 drug resistant mutants are transmitted and spreading from the ART-undergoing populations to the new infectors.



Acquired immune deficiency syndrome


antiretroviral therapy


drug-resistance mutations


Human immunodeficiency virus type 1


Injection drug user


Maximum likelihood


Non-Nucleoside Reverse Transcriptase Inhibitors


Nucleoside Reverse Transcriptase Inhibitors


Protease Inhibitors


transmitted drug resistance


transmitted drug resistance mutations


  1. Cohen MS, Chen YQ, McCauley M, Gamble T, Hosseinipour MC, Kumarasamy N, et al. Prevention of HIV-1 infection with early antiretroviral therapy. N Engl J Med. 2011;365:493–505.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Attia S, Egger M, Müller M, Zwahlen M, Low N. Sexual transmission of HIV according to viral load and antiretroviral therapy: systematic review and meta-analysis. AIDS. 2009;23:1397–404.

    Article  PubMed  Google Scholar 

  3. Porco TC, Martin JN, Page-Shafer KA, Cheng A, Charlebois E, Grant RM, et al. Decline in HIV infectivity following the introduction of highly active antiretroviral therapy. AIDS. 2004;18:81–8.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Pham QD, Wilson DP, Law MG, Kelleher AD, Zhang L. Global burden of transmitted HIV drug resistance and HIV-exposure categories: a systematic review and meta-analysis. AIDS. 2014;28(18):2751–62.

    Article  PubMed  Google Scholar 

  5. WHO HIV/AIDS Programme. The HIV drug resistance report-2012.

  6. Sun X, Nan J, Guo Q. AIDS and HIV infection in China. AIDS. 1994;8(Suppl 2):S55–9.

    PubMed  Google Scholar 

  7. Wang B, Liang Y, Wang Y, Feng Y, Li Y, et al. HIV prevalence and phylogenetic characteristics among entry travelers in Xishuangbanna prefecture, Yunnan province, between 2003 and 2012. J Med Virol. 2017;89(6):1112–5.

    Article  CAS  PubMed  Google Scholar 

  8. Wang B, Liang Y, Feng Y, Li Y, Wang Y, et al. Prevalence of human immunodeficiency virus 1 infection in the last decade among entry travelers in Yunnan Province, China. BMC Public Health. 2015;15:362.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Shan D, Sun J, Khoshnood K, Fu J, Duan S, Jiang C, et al. The impact of comprehensive prevention of mother-to-child HIV transmission in Dehong prefecture, Yunnan province, 2005-2010: a hard-hit area by HIV in southern China. Int J STD AIDS. 2014;25(4):253–60.

    Article  PubMed  Google Scholar 

  10. Zheng XW, Zhang JP, Tian CQ, Cheng HH, Yang XZ, Duan S, et al. Cohort study of HIV infection among drug users in Ruili, Longchuan and Luxi of Yunnan Province, China`. Biomed Environ Sci. 1993;6(4):348–51.

    CAS  PubMed  Google Scholar 

  11. Zheng X, Tian C, Choi KH, Zhang J, Cheng H, Yang X, et al. Injecting drug use and HIV infection in Southwest China. AIDS. 1994;8(8):1141–7.

    Article  CAS  PubMed  Google Scholar 

  12. Wang Y, Liang Y, Feng Y, Wang B, Li Y, et al. HIV-1 prevalence and subtype/recombinant distribution among travelers entering China from Vietnam at the HeKou port in the Yunnan province, China, between 2003 and 2012. J Med Virol. 2015;87(9):1500–9.

    Article  CAS  PubMed  Google Scholar 

  13. Aiyar A. The use of CLUSTAL W and CLUSTAL X for multiple sequence alignment. Methods Mol Biol. 2000;132:221–41.

    CAS  PubMed  Google Scholar 

  14. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Shafer RW, Rhee SY, Pillay D, Miller V, Sandstrom P, Schapiro JM, et al. HIV-1 protease and reverse transcriptase mutations for drug resistance surveillance. AIDS. 2007;21:215–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–1874.

    Article  CAS  PubMed  Google Scholar 

  17. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, et al. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999;73:152–60.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. Su YZ, Ma YL, Jia MH, He X, Yang L, Cheng CL, et al. Update on diversity and distribution of HIV-1 subtypes in Yunnan province. Epidemiol Infect. 2013;141(11):2418–27.

    Article  PubMed  Google Scholar 

  19. Li L, Sun G, Liang S, Li J, Li T, Wang Z, et al. Different distribution of HIV-1 subtype and drug resistance were found among treatment naïve individuals in Henan, Guangxi, and Yunnan province of China. PLoS One. 2013;8(10):e75777.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Pham HV, Ishizaki A, Nguyen CH, Saina MC, Hoang HT, Tran VT, et al. Change in the prevalence of HIV-1 and the rate of transmitted drug-resistant HIV-1 in Haiphong, northern Vietnam. AIDS Res Hum Retrovir. 2015;31(7):757–9.

    Article  PubMed  Google Scholar 

  21. Frentz D, Boucher CA, van de Vijver DA. Temporal changes in the epidemiology of transmission of drug-resistant HIV-1 across the world. AIDS Rev. 2012;14(1):17–27.

    PubMed  Google Scholar 

  22. Li H, Zhong M, Guo W, Zhuang D, Li L, Liu Y, et al. Prevalence and mutation patterns of HIV drug resistance from 2010 to 2011 among ART-failure individuals in the Yunnan Province, China. PLoS One. 2013;8(8):e72630.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Li L, Sun B, Zeng H, Sun Z, Sun G, Yang R. Relatively high prevalence of drug resistance among antiretroviral-naive patients from Henan, Central China. AIDS Res Hum Retrovir. 2014;30(2):160–4.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Liu J, Jiao Y, Li L, Zeng Z, Li C, Li S, et al. Research on the mutations of relevant drug resistance gene sites in newly confirmed HIV-1 infectors. Chinese J Health Lab Technol. 2016;5:753–5. (in Chinese)

    Google Scholar 

  25. Chen M, Ma Y, Su Y, Yang L, Zhang R, et al. HIV-1 genetic characteristics and transmitted drug resistance among men who have sex with men in Kunming, China. PLoS One. 2014;9(1):e87033.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Chen M, Wang J, Xing H, Ma Y, Yao S, Chen H, et al. Study on the threshold of HIV-1 drug resistance in Dehong prefecture of Yunnan province in 2014. Zhonghua Liu Xing Bing Xue Za Zhi. 2014;35(8):909–12. (in Chinese)

    CAS  PubMed  Google Scholar 

  27. Sanguansittianant S, Nooroon N, Phaengchomduan P, Ammaranond P. Trends in prevalence of HIV-1 drug resistance in Thailand 2009-2010. J Clin Lab Anal. 2013;27(5):346–53.

    Article  CAS  PubMed  Google Scholar 

  28. Barth RE, van der Loeff MF, Schuurman R, Hoepelman AI, Wensing AM. Virological follow-up of adult patients in antiretroviral treatment programmes in sub-Saharan Africa: a systematic review. Lancet Infect Dis. 2010;10(3):155–66.

    Article  CAS  PubMed  Google Scholar 

  29. Hosseinipour MC, Gupta RK, Van Zyl G, Eron JJ, Nachega JB. Emergence of HIV drug resistance during frst- and second-line antiretroviral therapy in resource-limited settings. J Infect Dis. 2013;207(Suppl 2):S49–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Jiao Y, Li S, Li Z, Zhang Z, Zhao J, Li L, et al. HIV-1 transmitted drug resistance-associated mutations and mutation co-variation in HIV-1 treatment-naïve MSM from 2011 to 2013 in Beijing, China. BMC Infect Dis. 2014;14:689.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Lu X, Zhao H, Zhang Y, Wang W, Zhao C, Li Y, et al. HIV-1 drug-resistant mutations and related risk factors among HIV-1-positive individuals experiencing treatment failure in Hebei Province, China. AIDS Res Ther. 2017;14(1):4.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Sui H, Gui T, Jia L, Guo W, Han J, Liu Y, et al. Different frequencies of drug resistance mutations among HIV-1 subtypes circulating in China: a comprehensive study. PLoS One. 2014;9(3):e91803.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


The authors would like to express immense gratitude to all patients associated with this study.


This study was supported by grants from the National Natural Science Foundation of China (81460509); the Science and Technology Project (2016RA014), the Key Science and Technology Planning Project of Yunnan Provincial Science and Technology Department (2016FC005 and 2016 BC005).

Availability of data and materials

The sequences used in this study were submitted to NCBI GenBank (accession numbers: MG787428-MG787459).

Author information

Authors and Affiliations



XX and WB designed the study; XQ, WB, XZ and LS performed statistical analysis and drafted the first version of this manuscript; YS, ZT and ZA participated in the writing process; SH and QW collected sample and the demographic data. All other authors read and approved the final manuscript.

Corresponding authors

Correspondence to Binghui Wang or Xueshan Xia.

Ethics declarations

Ethics approval and consent to participate

All the participants were informed at enrolment and written consent was received before sample collection. The research was approved by the Institutional Ethical Committee of Kunming University of science and technology (NO.2014sk026). Consent was obtained from a parent or guardian on behalf of any participants under the age of 16.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Figure S1. Maximum likelihood method was chosen for phylogenetic analysis and constructed phylogenetic tree using all the pol gene sequences. (PDF 133 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xuan, Q., Liang, S., Qin, W. et al. High prevalence of HIV-1 transmitted drug resistance among therapy-naïve Burmese entering travelers at Dehong ports in Yunnan, China. BMC Infect Dis 18, 211 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: