This article has Open Peer Review reports available.
Diverse forms of HIV-1 among Burmese long-distance truck drivers imply their contribution to HIV-1 cross-border transmission
© Zhou et al.; licensee BioMed Central Ltd. 2014
Received: 19 February 2014
Accepted: 19 August 2014
Published: 26 August 2014
The China-Myanmar border is a particularly interesting region that has very high prevalence of and considerable diversity of HIV-1 recombinants. Due to the transient nature of their work, long-distance truck drivers (LDTDs) have a comparatively high potential to become infected with HIV-1 and further spread virus to other individuals in the area they travel within. In this study, we hypothesized that Burmese LDTDs crossing the China-Myanmar border frequently may potentially be involved in the cross-border transmission of HIV, and contribute to the extremely high prevalence of HIV-1 inter-subtype recombinants in this border region.
A molecular epidemiology study was conducted among 105 Burmese LDTDs between 2008 and 2010. HIV-1 genetic fragments including p17, pol, vif-vpr, vpr-env, and C2V3 were amplified and sequenced. The subtype characterization and HIV-1 transmission were determined by both phylogenetic and phylogeographic analyses.
Diverse forms of HIV-1, including subtypes CRF01_AE (41.9%), C (8.6%), B (4.8%), CRF02_AG (1.0%), and inter-subtype recombinants (33.3%), as well as dual infection (10.5%), were detected among the tested LDTDs. Phylogeographic analyses based on pure subtype revealed that 77.8% Burmese LDTDs acquired HIV-1 infection in Yunnan, and the others in Myanmar. Both the C-related and CRF01_AE-related recombinants from these LDTDs appeared to have close genetic relationship with those from IDUs in Myanmar and Dehong.
Burmese LDTDs may contribute to HIV-1 transmission along the China-Myanmar border. The results may provide some new perspective for understanding the on-going generation and prevalence of HIV-1 recombinants in the border region.
Myanmar has one of the highest rates of HIV-1 in southeast Asia, with CRF01_AE and B being the predominant subtypes among both injection drug users (IDUs) and heterosexual individuals in south and central Myanmar [1–3]. In northern Myanmar, CRF01_AE, B and C-involved unique recombinant forms (URFs) were the most prevalent strains among IDUs . China’s Yunnan province, which shares some 1977 kilometers of its border with Myanmar, has the highest number of HIV-1 infected individuals in China, accounting for 34.8% of national total cases as reported . In Yunnan, multiple HIV-1 genotypes, including CRF01_AE, B, C, CRF07_BC, CRF08_BC, as well as various URFs, have been identified [6, 7]. Within Yunnan, Dehong prefecture along the northern borders of Myanmar has considerably high HIV-1 prevalence among IDUs, and the proportion of HIV-1 recombinants reached 86.1%, very similar to that of northern Myanmar [7–9]. Why such an extremely high prevalence of HIV-1 recombinants occurred in both sides of the China-Myanmar border remains unclear.
Ruili city in Dehong prefecture is an important port linking China and Myanmar, serving as the major gateway for trade between both sides. Drugs (in particular heroin) from the infamous “Golden Triangle” are often illegally transferred to China via this port. Along the drug trafficking routes, HIV-1 has rapidly spread among IDUs. Unsurprisingly, Ruili became the first recorded site of HIV-1 outbreak in China . About 45.2% Burmese IDUs and 55.6% Chinese IDUs cross the China-Myanmar border each month to inject drugs, often with sharing needles and syringes [Dr. Lin Duo, personal communication]. This breadth of cross-border drug usage may partly explain the reason why both sides of the China-Myanmar border have similar HIV-1 epidemic and extremely high prevalence of HIV-1 recombinants.
Earlier studies supposed that HIV-1 expanded from Myanmar to Yunnan through these drug trafficking routes, but our recent phylogeographic analyses revealed that parental subtypes (e.g. B, C and CRF01_AE) of URFs_01BC circulating among Burmese IDUs were most likely transmitted from Dehong to Myanmar during 1990s . Interestingly, the prevalence of CRF01_AE-related URFs was remarkably higher in northern Myanmar than in Dehong, Yunnan, China. IDUs may be a main factor fueling HIV-1 bidirectional transmission and generation of many URFs, but could not simply explain this substantial difference in HIV-1 recombination patterns between both sides of the border. Therefore, we were curious as to whether or not other groups aside from IDUs may also have contributed to the rapid mixing of HIV-1 subtypes B, C and CRF01_AE, and the generation of the extremely high prevalence of HIV-1 recombinants in the China-Myanmar border region.
Long-distance truck drivers (LDTDs) are at an increased risk of infection with HIV, since they often stay away from their hometowns for long periods of time, and therefore less frequently visit a local, monogamous sexual partner and more frequently use commercial sex workers (CSWs) and have multiple sexual partners [12, 13]. A recent study showed that 5% of Burmese LDTDs crossing the China-Myanmar border had sex with occasional sexual partners during long-distance truck driving, and some LDTDs self-reported to inject drugs. Among this 5% population, 7.1% reported that they did not use condoms when engaging in sexual activity. HIV-1 prevalence was reported to be as high as 3.5-8.8% among Burmese LDTDs in general [14, 15]. Whether these LDTDs were involved in the high prevalence of HIV-1 recombinants in the China-Myanmar border region needs to be determined.
Here, we conducted a molecular epidemiological study of HIV-1 among Burmese LDTDs, and investigated their phylogenetic relationship with strains circulating among IDUs in the China-Myanmar border region.
Study population and specimens
This study was conducted in accordance with the principles of the Declaration of Helsinki and relevant international regulations, and approved by the Ethics Committee of Kunming Institute of Zoology, Chinese Academy of Sciences (Approval ID: SWYX2008010). Under the support of Australia HIV/AIDS Asia Regional Program (HAARP) and Yunnan HAARP project, 105 blood specimens were collected from HIV-positive Burmese LDTDs entering and exiting Ruili port between 2008 and 2010, after they provided written informed consent. Demographic data of all participants were obtained via an interview. The plasma were separated and stored at -70°C until use. All samples were serologically determined as HIV infection by an enzyme-linked immunosorbent assay (ELISA) kit (Wantai Biological Pharmacy Enterprise Co. Ltd, Beijing, China) and confirmed by another ELISA kit (Kehua bioengineering Co. Ltd, Shanghai, China). No HIV-2 infection was found among any of the studied participants.
HIV-1 amplification, cloning and sequencing
Five gene fragments including p17 (520 bp), pol (1273 bp), vif-vpr (702 bp), vpr-env (747 bp) and C2V3 (513 bp) were amplified from viral RNA by RT-nested PCR as described previously . The amplicons were purified and sent to sequence directly without prior cloning (Invitrogen Trading, Shanghai, China). If direct sequencing failed due to the presence of numerous double peaks in the chromatogram, the amplicons were cloned into pMD19-T vectors (TaKaRa Biotechnology) and several clones were selected randomly for re-sequencing. Furthermore, if one amplified fragment of certain samples had a subtyping result different from all other four fragments, this fragment was also subjected to TA-clone and 20–25 colonies were selected randomly for sequencing.
All obtained sequences were screened by the HIV BLAST tool to confirm lack of laboratory contamination and then were aligned with the HIV-1 reference strains obtained from the Los Alamos HIV database (http://www.hiv.lanl.gov/content/sequence/NEWALIGN/align.html), using BioEdit 7.0.9 (http://www.mbio.ncsu.edu/bioedit/bioedit.html) and MUSCLE (http://www.ebi.ac.uk/Tools/msa/muscle/). Phylogenetic trees were constructed by the neighbor-joining (NJ) method with 1000 bootstrap replicates, using Kimura’s two parameter model in MEGA 5.0 (http://www.megasoftware.net/). To detect HIV-1 recombinants, bootscan analysis was performed using SimPlot 3.5.1 (http://sray.med.som.jhmi.edu/SCRoftware/simplot/).
Phylogeographic analyses for CRF01_AE, B, C and URFs
To investigate the infection place of HIV-1 strains circulating among LDTDs, phylogeographic analyses were performed based on p17 sequences using BEAST 1.6.2 [16, 17]. Apart from p17 sequences obtained in this study, all available sequences with known sample year from Myanmar, China, Thailand, Vietnam, and India were also included in the analyses. The time scaled maximum clade credibility (MCC) trees were constructed for CRF01_AE, subtype B, C and URFs via the codon-based substitution model (SRD06), the uncorrelated log-normal relaxed clock model, and the constant population size. To ensure effective sample size (ESS) values >200, Markov Chain Monte Carlo (MCMC) chains for each subtype sequence set were run for 200 million generations. All trees were visualized using FigTree 1.3.1.
Nucleotide sequence accession numbers
All the sequences reported in this paper are available in GenBank under accession numbers from KC913597 to KC914092.
Demographics of the study population
All 105 Burmese LDTDs were males with an average age of 34.3 (range: 20–56 years). More than half of participants (58.1%) came from Mandalay, a city in central Myanmar. Sixty five participants (61.9%) self-reported the potential routes for HIV-1 infection including 59 (56.2%) via sexual contact, two (1.9%) via blood transfusion, two (1.9%) via sexual contact and/or transfusions, and two (1.9%) via sexual and/or injection drug use (IDU). None self-reported that IDU was the only risk behavior (see Additional file 1).
Subtype characterization of HIV-1 strains circulating among Burmese LDTDs
From 105 samples, we obtained 100 (95.2%) p17, 86 (81.9%) pol, 97 (92.4%) vif-vpr, 85 (81.0%) vpr-env and 90 (85.7%) C2V3 sequences. Because of an overlap of 156 nt between vif-vpr and vpr-env sequences, they were spliced into a 1328 nt sequence named vif-env when both were available for the same sample. As a result, 80 vif-env sequences were obtained (see Additional file 2).
High prevalence of HIV-1 dual infection among LDTDs
Origin of HIV-1 among LDTDs
Population movement is commonly identified as a key factor associated with transmission of human viruses, especially HIV-1 . LDTDs are a particular group with high potential to be infected with HIV-1 and further spread the virus to other individuals due to the long distances they travel, the long duration of staying away from hometowns, frequent use of CSWs, and some of them also being IDUs. Several previous studies have reported the association of the high risk behaviors of LDTDs with HIV-1 infection [15, 19]. However, there is less paper to characterize the molecular epidemiology of HIV-1 among LDTDs.
The China-Myanmar border region is a particularly interesting area that serves as a key channel for illegal drug trafficking from the “Golden triangle” to other countries/regions and has very high prevalence of HIV-1 and other viruses (e.g. HCV and HBV) [4, 9, 20]. In particular, over 86% HIV-1 strains circulating among IDUs in this region belong to URFs, which raises a question why such high proportion of HIV-1 recombinants occurred among them in this region. Furthermore, CRF01_AE-related URFs were detected among 45.6% of recombinants in northern Myanmar, significantly higher than 6.6% in its counterpart side (i.e. Dehong, Yunnan province of China) [4, 9]. IDUs were postulated as a main population responsible for the generation of so many URFs since they often injected drugs together across the border and shared needles/syringes. However, the behaviors of IDUs do not explain why such high proportion (45.6%) of CRF01_AE-related URFs occurs in northern Myanmar, implying the presence of other high-risk population responsible for this difference. In this study, we collected the samples from 105 Burmese LDTDs who drive between Mandalay, Myanmar and Ruili, China, and performed a molecular epidemiological investigation of HIV-1 among this cohort. The results showed that diverse forms of HIV-1, including subtypes CRF01_AE (41.9%), C (8.6%), B (4.8%), CRF02_AG (1.0%), and some inter-subtype recombinants (33.3%), were circulating among these LDTDs.
Phylogeographic analyses of the pure subtypes showed that majority (77.8%, 42/54) of Burmese LDTDs were infected with HIV-1 (five subtype B, six subtype C and 31 CRF01_AE) in Yunnan, and the others were infected in Myanmar (one subtype C and 9 CRF01_AE) (Figure 3). Likewise, MCC trees based on p17 regions derived from the URFs of LDTDs also supported that HIV-1 C and CRF01_AE in LDTDs have multiple origins including Dehong and Myanmar (Figure 4). On a broader scale, no matter where Burmese LDTDs acquired their infection, they have high potential to further spread the viruses to their wives and/or sexual partners in Myanmar and Yunnan. Therefore, Burmese LDTDs can mediate bi-directional transmission of HIV-1 between heterosexuals and IDUs and between China and Myanmar. This bi-directional transmission pattern of HIV-1 will facilitate the mixing of different HIV-1 forms among different high risk groups, and accelerate the generation of new recombinants, as observed among IDUs in the China-Myanmar border region. This concept was further supported by the detection of all CRF01_AE, B and C involved recombination forms (e.g. B/C, CRF01/B, CRF01/C, and CRF01/B/C) among Burmese LDTDs (Figure 1B). One major limitation is that because the samples from female sexual workers (FSWs) in the China-Myanmar border region were unavailable in this study, we were unable to obtain HIV-1 subtype characterization in this cohort, and could not find the direct evidence for the transmission association between LDTDs and FSWs. Hence we just further determined their role in HIV-1 cross-border transmission and their association with IDUs via phylogeographic analyses. To understand HIV genetic diversity and transmission complexity in the China-Myanmar border region, more concern should be paid to female sexual workers in future.
Dual infection is an essential condition for the generation of HIV-1 recombinants [21, 22]. Currently, inter-subtype recombinants are responsible for about 11.0% HIV-1 infection throughout the Asia (http://www.hiv.lanl.gov/), and dual/triple infections have been detected in Thailand and Taiwan, where multiple subtypes are co-existing [23–25]. In the China-Myanmar border area, 86.1-87.9% HIV-1 infected IDUs carried recombinants [4, 9]. Such high prevalence of recombinants implies a very frequent occurrence of dual infection. To date, however, no dual infection case was reported in this region. We found that 10.5% Burmese LDTDs were dually infected with various subtypes, implying frequent exposure of LDTDs to HIV-1. Among them, dual infection with CRF01_AE and C accounted for 72.7% (Figure 1C), which not only provided an explanation for why CRF01_AE/C recombinants were the most predominant recombination form (37.1%) among Burmese LDTDs, but also predicted an on-going generation of more new CRF01_AE/C recombinants among this population. In addition, CRF01_AE was involved in 90.9% cases of dual infection. These, together with the fact of high prevalence of CRF01_AE (41.9%), support the close association of Burmese LDTDs with sexual transmission group.Interestingly, two cases were found to be dual infection with certain pure subtype and its self-derived recombinants: one with CRF01_AE and CRF01_AE/B recombinants and the other with subtype C and B/C recombinant (Figure 2). These might represent a real-time observation on new recombinant generation. Because only the samples with sequences containing double peaks in Sanger sequencing were subjected to TA-cloning and in general, only 3–25 clones for each sample were selected for sequencing, this may lead to an underestimation of prevalence of dual infection among this cohort. The real dual infection rate may be higher when more samples are subjected to TA-cloning and sequencing, or ultra-deep sequencing, which is a robust tool to found low-frequency variants.
This study firstly represented the molecular epidemiological characterization of HIV-1 among Burmese LDTDs. High prevalence of HIV-1 CRF01_AE and inter-subtype recombinants, as well as dual infection, were detected among them. Phylogeographic analyses revealed that Burmese LDTDs may have made contribution in HIV-1 transmission between the heterosexuals and IDUs and between China and Myanmar in general. These findings may provide new insight into understanding of the extensive and complex of HIV-1 recombination in the China-Myanmar border area (Figure 5). In a broader sense, our results also provide new information relevant to policy-making decisions regarding the control of HIV/AIDS in southeast Asia.
The authors gratefully thank the participants in this study and the staff of the Australia HIV/AIDS Asia Regional Program (HAARP) and the Yunnan HAARP project helping collecting samples. We also thank Dr. Andrew Willden for editing the manuscript. This work was supported in part by grants from the NSFC (81271892, U1302224, U1202228, 81172876), Key Scientific and Technological Program of China (2012ZX10001-006, 2012ZX10001-007), National Basic Research Program of China (2012CBA01305), the Chinese Academy of Sciences (KSCX2-EW-R-13). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
- Kusagawa S, Sato H, Watanabe S, Nohtomi K, Kato K, Shino T, Thwe M, OO KY, Lwin S, Mra R, Kywe B, Yamazaki S, Takebe Y: Genetic and serologic characterization of HIV type 1 prevailing in Myanmar (Burma). AIDS Res Hum Retroviruses. 1998, 14: 1379-1385.View ArticlePubMedGoogle Scholar
- Motomura K, Kusagawa S, Kato K, Nohtomi K, Lwin HH, Tun KM, Thwe M, Oo KY, Lwin S, Kyaw O, Zaw M, Nagai Y, Takebe Y: Emergence of new forms of human immunodeficiency virus type 1 intersubtype recombinants in central Myanmar. AIDS Res Hum Retroviruses. 2000, 16: 1831-1843.View ArticlePubMedGoogle Scholar
- Takebe Y, Motomura K, Tatsumi M, Lwin HH, Zaw M, Kusagawa S: High prevalence of diverse forms of HIV-1 intersubtype recombinants in Central Myanmar: geographical hot spot of extensive recombination. AIDS. 2003, 17: 2077-2087.View ArticlePubMedGoogle Scholar
- Pang W, Zhang C, Duo L, Zhou YH, Yao ZH, Liu FL, Li H, Tu YQ, Zheng YT: Extensive and complex HIV-1 recombination between B, C and CRF01_AE among IDUs in southeast Asia. AIDS. 2012, 26: 1121-1129.View ArticlePubMedGoogle Scholar
- Jia M, Luo H, Ma Y, Wang N, Smith K, Mei J, Lu R, Lu J, Fu L, Zhang Q, Wu Z, Lu L: The HIV epidemic in Yunnan Province, China, 1989–2007. J Acquir Immune Defic Syndr. 2010, 53 (Suppl 1): S34-S40.View ArticlePubMedGoogle Scholar
- Zhang Y, Lu L, Ba L, Liu L, Yang L, Jia M, Wang H, Fang Q, Shi Y, Yan W, Chang G, Zhang L, Ho DD, Chen Z: Dominance of HIV-1 subtype CRF01_AE in sexually acquired cases leads to a new epidemic in Yunnan province of China. PLoS Med. 2006, 3: e443-View ArticlePubMedPubMed CentralGoogle Scholar
- Li L, Chen L, Yang S, Li T, Li J, Liu Y, Jia L, Yang B, Bao Z, Li H, Wang X, Zhuang D, Liu S, Li J: Recombination form and epidemiology of HIV-1 unique recombinant strains identified in Yunnan, China. PLoS One. 2012, 7: e46777-View ArticlePubMedPubMed CentralGoogle Scholar
- Yang R, Xia X, Kusagawa S, Zhang C, Ben K, Takebe Y: On-going generation of multiple forms of HIV-1 intersubtype recombinants in the Yunnan Province of China. AIDS. 2002, 16: 1401-1407.View ArticlePubMedGoogle Scholar
- Han X, An M, Zhao B, Duan S, Yang S, Xu J, Zhang M, McGoogan JM, Takebe Y, Shang H: High prevalence of HIV-1 intersubtype B’/C recombinants among injecting drug users in Dehong, China. PLoS One. 2013, 8: e65337-View ArticlePubMedPubMed CentralGoogle Scholar
- du Guerny J, Hsu L, Hong C: Population Movement and HIV/AIDS: the Case of Ruili, Yunnan, China. UNDP Southeast Asia HIV and Development Programme. 2003, [Available at: http://www.hivdevelopment.org/]Google Scholar
- Liu J, Jia Y, Xu Q, Zheng YT, Zhang C: Phylodynamics of HIV-1 unique recombinant forms in China-Myanmar border: implication for HIV-1 transmission to Myanmar from Dehong, China. Infect Genet Evol. 2012, 12: 1944-1948.View ArticlePubMedGoogle Scholar
- Manjunath JV, Thappa DM, Jaisankar TJ: Sexually transmitted diseases and sexual lifestyles of long-distance truck drivers: a clinico-epidemiologic study in south India. Int J STD AIDS. 2002, 13: 612-617.View ArticlePubMedGoogle Scholar
- Rao KS, Pilli RD, Rao AS, Chalam PS: Sexual lifestyle of long distance lorry drivers in India: questionnaire survey. BMJ. 1999, 318: 162-163.View ArticlePubMedPubMed CentralGoogle Scholar
- Liang YB, Li YP, Wang CM, Li XY, Yang RQ, Su H, Wang GL, Wang ZJ: An analysis of HIV detection results among the Burmese truck drivers entering and exiting Ruili port. J Chin AIDS/STD. 2010, 16: 372-374.Google Scholar
- Deng L, Zhang MY, Wang GL, Li SH, Xue HM, Duo L: Analysis on the effective elements relating Myanmar cross-border long-distance truck driver cross-border and risky behavior. Soft Sci Health. 2011, 25: 835-837.Google Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-View ArticlePubMedPubMed CentralGoogle Scholar
- Lemey P, Rambaut A, Drummond AJ, Suchard MA: Bayesian phylogeography finds its roots. PLoS Comput Bio. 2009, 5: e1000520-View ArticleGoogle Scholar
- Deane KD, Parkhurst JO, Johnston D: Linking migration, mobility and HIV. Trop Med Int Health. 2010, 15: 1458-1463.View ArticlePubMedGoogle Scholar
- Sanjeev K, S K G, S K B: A study of knowledge, sexual behaviour and practices regarding HIV/AIDS among long distance truck drivers. Ind J Public Health. 2009, 53: 243-245.Google Scholar
- Zhou YH, Liu FL, Yao ZH, Duo L, Li H, Sun Y, Zheng YT: Comparison of HIV-, HBV-HCV- and co-infection prevalence between Chinese and Burmese intravenous drug users of the China-Myanmar border region. PLoS One. 2011, 6: e16349-View ArticlePubMedPubMed CentralGoogle Scholar
- Nájera R, Delgado E, Pérez-Alvarez L, Thomson MM: Genetic recombination and its role in the development of the HIV-1 pandemic. AIDS. 2002, 16 (Suppl 4): S3-S16.View ArticlePubMedGoogle Scholar
- Templeton AR, Kramer MG, Jarvis J, Kowalski J, Gange S, Schneider MF, Shao Q, Zhang GW, Yeh MF, Tsai HL, Zhang H, Markham RB: Multiple-infection and recombination in HIV-1 within a longitudinal cohort of women. Retrovirology. 2009, 6: 54-View ArticlePubMedPubMed CentralGoogle Scholar
- Hu DJ, Subbarao S, Vanichseni S, Mock PA, Ramos A, Nguyen L, Chaowanachan T, Griensven F, Choopanya K, Mastro TD, Tappero JW: Frequency of HIV-1 dual subtype infections, including intersubtype superinfections, among injection drug users in Bangkok, Thailand. AIDS. 2005, 19: 303-308.PubMedGoogle Scholar
- Chen YJ, Lee CM, Chen M, Chuang SY, Liu HF, Wong WW, Lin YH, Tsai HC, Wang JH, Chen YM: Molecular epidemiology of HIV-1 infection in Taiwan from 2005 to 2008: further spread of CRF07_BC and emergence of CRF07_BC/subtype B dual infection. J Acquir Immune Defic Syndr. 2012, 59: 438-446.View ArticlePubMedGoogle Scholar
- van der Kuyl AC, Zorgdrager F, Jurriaans S, Back NK, Prins JM, Brinkman K, van Eeden A, Bakker M, Cornelissen M: Incidence of human immunodeficiency virus type 1 dual infections in Amsterdam, The Netherlands, during 2003–2007. Clin Infect Dis. 2009, 48: 973-978.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/14/463/prepub
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 (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.