Volume 36 Issue 5
May  2023
Turn off MathJax
Article Contents

CEN Ping, LIANG Hua Yue, YANG Yuan, ZHANG Fei, YANG Shi Xiong, MO Ju Cong, FENG Yi, HUANG Jie Gang, NING Chuan Yi, HUANG Chun Yuan, YANG Yao, LIANG Na, LIANG Bing Yu, YE Li, LIANG Hao. HIV-1 Transmission among Injecting Drug Users is Principally Derived from Local Circulating Strains in Guangxi, China[J]. Biomedical and Environmental Sciences, 2023, 36(5): 418-430. doi: 10.3967/bes2023.051
Citation: CEN Ping, LIANG Hua Yue, YANG Yuan, ZHANG Fei, YANG Shi Xiong, MO Ju Cong, FENG Yi, HUANG Jie Gang, NING Chuan Yi, HUANG Chun Yuan, YANG Yao, LIANG Na, LIANG Bing Yu, YE Li, LIANG Hao. HIV-1 Transmission among Injecting Drug Users is Principally Derived from Local Circulating Strains in Guangxi, China[J]. Biomedical and Environmental Sciences, 2023, 36(5): 418-430. doi: 10.3967/bes2023.051

HIV-1 Transmission among Injecting Drug Users is Principally Derived from Local Circulating Strains in Guangxi, China

doi: 10.3967/bes2023.051
Funds:  National Natural Science Foundation of China [grant no. 82060610]; Guangxi Scientific and Technological Key Project [Gui Ke AB19245038]; National Natural Science Foundation of China [grant no. 82103899]; Guangxi Scientific and Technological Key Project [Guike 2022AC23005, 2022AC20031,2022JJA141110]; and Science and Technology Project of Nanning [20223051]
More Information
  • Author Bio:

    CEN Ping, female, born in 1982, PhD Candidate, majoring in epidemiology of AIDS

    LIANG Hua Yue, male, born in 1995, Master, majoring in epidemiology of AIDS

  • Corresponding author: LIANG Hao, Professor, PhD, Tel: 86-771-5323713, E-mail: lianghao@gxmu.edu.cn; YE Li, Professor, PhD, Tel: 86-771-5344361, E-mail: yeli@gxmu.edu.cn; LIANG Bing Yu; Professor, PhD, Tel: 86-771-5334215, E-mail: liangbingyu@gxmu.edu.cn
  • LIANG Hao, YE Li, and LIANG Bing Yu conceived and designed the study. YANG Shi Xiong and MO Ju Cong collected the specimen information. YANG Yuan, ZHANG Fei, HUANG Chun Yuan, YANG Yao, and LIANG Na performed experimental work. FENG Yi, HUANG Jie Gang, and NING Chuan Yi provided insight into the experimental design and data analysis. CEN Ping and LIANG Hua Yue performed data analysis and wrote the first draft. All authors approved the manuscript.
  • The authors declare that there are no conflicts of interest.
  • &These authors contributed equally to this work.
  • Received Date: 2022-08-01
  • Accepted Date: 2022-10-09
  •   Objective   The mode of human immunodeficiency virus (HIV) transmission via injection drug use (IDU) still exists, and the recent shift in IDU-related transmission of HIV infection is largely unknown. The purpose of this study was to analyze the spatiotemporal sources and dynamics of HIV-1 transmission through IDU in Guangxi.   Methods   We performed a molecular epidemiological investigation of infections across Guangxi from 2009 to 2019. Phylogenetic and Bayesian time-geographic analyses of HIV-1 sequences were performed to confirm the characteristics of transmission between IDUs in combination with epidemiological data.   Results   Among the 535 subjects, CRF08_BC (57.4%), CRF01_AE (28.4%), and CRF07_BC (10.7%) were the top 3 HIV strains; 72.6% of infections were linked to other provinces in the transmission network; 93.6% of sequence-transmitted strains were locally endemic, with the rest coming from other provinces, predominantly Guangdong and Yunnan; 92.1% of the HIV transmission among people who inject drugs tended to be transmitted between HIV-positive IDUs.   Conclusion  HIV recombinants were high diversity, and circulating local strains were the transmission sources among IDUs in Guangxi. However, there were still cases of IDUs linked to other provinces. Coverage of traditional prevention strategies should be expanded, and inter-provincial collaboration between Guangxi, Yunnan, and Guangdong provinces should be strengthened.
  • LIANG Hao, YE Li, and LIANG Bing Yu conceived and designed the study. YANG Shi Xiong and MO Ju Cong collected the specimen information. YANG Yuan, ZHANG Fei, HUANG Chun Yuan, YANG Yao, and LIANG Na performed experimental work. FENG Yi, HUANG Jie Gang, and NING Chuan Yi provided insight into the experimental design and data analysis. CEN Ping and LIANG Hua Yue performed data analysis and wrote the first draft. All authors approved the manuscript.
    The authors declare that there are no conflicts of interest.
    &These authors contributed equally to this work.
  • 加载中
  • [1] Sheng L, Cao WK. HIV/AIDS epidemiology and prevention in China. Chin Med J (Engl), 2008; 121, 1230−6. doi:  10.1097/00029330-200807010-00016
    [2] Des Jarlais DC, Johnston P, Friedmann P, et al. Patterns of HIV prevalence among injecting drug users in the cross-border area of Lang Son Province, Vietnam, and Ning Ming County, Guangxi Province, China. BMC Public Health, 2005; 5, 89. doi:  10.1186/1471-2458-5-89
    [3] Peng ZH, Cheng YJ, Reilly KH, et al. Spatial distribution of HIV/AIDS in Yunnan province, People's Republic of China. Geospat Health, 2011; 5, 177−82. doi:  10.4081/gh.2011.169
    [4] Ma SZ, Chen Y, Lai XL, et al. Predicting the HIV/AIDS epidemic and measuring the effect of AIDS Conquering Project in Guangxi Zhuang Autonomous Region. PLoS One, 2022; 17, e0270525. doi:  10.1371/journal.pone.0270525
    [5] Chen HH, Luo LH, Pan SW, et al. HIV epidemiology and prevention in southwestern China: trends from 1996-2017. Curr HIV Res, 2019; 17, 85−93. doi:  10.2174/1570162X17666190703163838
    [6] He X, Xing H, Ruan YH, et al. A comprehensive mapping of HIV-1 genotypes in various risk groups and regions across China based on a nationwide molecular epidemiologic survey. PLoS One, 2012; 7, e47289. doi:  10.1371/journal.pone.0047289
    [7] Tang HL, Cai C, Cui Y, et al. Epidemiological characteristics of newly diagnosed cases of HIV through injection drug use - China, 2012-2019. China CDC Wkly, 2020; 2, 920−3. doi:  10.46234/ccdcw2020.251
    [8] Zhang B, Yan XY, Li YJ, et al. Epidemic of HIV infection among persons who inject drugs in mainland China: a series, cross-sectional study. Harm Reduct J, 2021; 18, 63. doi:  10.1186/s12954-021-00511-6
    [9] Hemelaar J, Gouws E, Ghys PD, et al. Global trends in molecular epidemiology of HIV-1 during 2000-2007. AIDS, 2011; 25, 679−89. doi:  10.1097/QAD.0b013e328342ff93
    [10] Chen M, Ma YL, Chen HC, et al. Spatial clusters of HIV-1 genotypes in a recently infected population in Yunnan, China. BMC Infect Dis, 2019; 19, 669. doi:  10.1186/s12879-019-4276-9
    [11] Li K, Liu ML, Chen HH, et al. Using molecular transmission networks to understand the epidemic characteristics of HIV-1 CRF08_BC across China. Emerg Microbes Infect, 2021; 10, 497−506. doi:  10.1080/22221751.2021.1899056
    [12] Feng Y, Takebe Y, Wei HM, et al. Geographic origin and evolutionary history of China's two predominant HIV-1 circulating recombinant forms, CRF07_BC and CRF08_BC. Sci Rep, 2016; 6, 19279. doi:  10.1038/srep19279
    [13] Ge ZW, Feng Y, Zhang H, et al. HIV-1 CRF07_BC transmission dynamics in China: two decades of national molecular surveillance. Emerg Microbes Infect, 2021; 10, 1919−30. doi:  10.1080/22221751.2021.1978822
    [14] Nikolopoulos GK, Kostaki EG, Paraskevis D. Overview of HIV molecular epidemiology among people who inject drugs in Europe and Asia. Infect Genet Evol, 2016; 46, 256−68. doi:  10.1016/j.meegid.2016.06.017
    [15] Fan Q, Zhang JF, Luo MY, et al. Molecular genetics and epidemiological characteristics of HIV-1 epidemic strains in various sexual risk behaviour groups in developed eastern China, 2017-2020. Emerg Microbes Infect, 2022; 11, 2326−39. doi:  10.1080/22221751.2022.2119167
    [16] Liao HN, Tee KK, Hase S, et al. Phylodynamic analysis of the dissemination of HIV-1 CRF01_AE in Vietnam. Virology, 2009; 391, 51−6. doi:  10.1016/j.virol.2009.05.023
    [17] Wang XL, He X, Zhong P, et al. Phylodynamics of major CRF01_AE epidemic clusters circulating in mainland of China. Sci Rep, 2017; 7, 6330. doi:  10.1038/s41598-017-06573-6
    [18] Piyasirisilp S, McCutchan FE, Carr JK, et al. A recent outbreak of human immunodeficiency virus type 1 infection in southern China was initiated by two highly homogeneous, geographically separated strains, circulating recombinant form AE and a novel BC recombinant. J Virol, 2000; 74, 11286−95. doi:  10.1128/JVI.74.23.11286-11295.2000
    [19] Su L, Graf M, Zhang YZ, et al. Characterization of a virtually full-length human immunodeficiency virus type 1 genome of a prevalent intersubtype (C/B') recombinant strain in China. J Virol, 2000; 74, 11367−76. doi:  10.1128/JVI.74.23.11367-11376.2000
    [20] Kato K, Kusagawa S, Motomura K, et al. Closely related HIV-1 CRF01_AE variant among injecting drug users in northern Vietnam: evidence of HIV spread across the Vietnam-China border. AIDS Res Hum Retroviruses, 2001; 17, 113−23. doi:  10.1089/08892220150217201
    [21] Jiang Y, Zhang LC, Hou ZJ, et al. Prevalence of different genotypes of HIV-1 in injection drug users in China: a systematic review and meta-analysis. Curr HIV Res, 2019; 17, 240−57. doi:  10.2174/1570162X17666190919115036
    [22] Su L, Feng Y, Liang S, et al. The origin and spread of CRF85_BC, driven by heterosexual transmission among older people in Sichuan, China. BMC Infect Dis, 2020; 20, 772. doi:  10.1186/s12879-020-05488-4
    [23] Li JJ, Feng Y, Shen ZY, et al. HIV-1 transmissions among recently infected individuals in southwest china are predominantly derived from circulating local strains. Sci Rep, 2018; 8, 12831. doi:  10.1038/s41598-018-29201-3
    [24] Wang XF, Liu XS, Li F, et al. Epidemiological surveillance of HIV-1 transmitted drug resistance among newly diagnosed individuals in Shijiazhuang, northern China, 2014-2015. PLoS One, 2018; 13, e0198005. doi:  10.1371/journal.pone.0198005
    [25] Pang XW, Wei H, Huang JH, et al. Patterns and risk of HIV-1 transmission network among men who have sex with men in Guangxi, China. Sci Rep, 2021; 11, 513. doi:  10.1038/s41598-020-79951-2
    [26] Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol, 2007; 7, 214. doi:  10.1186/1471-2148-7-214
    [27] Gill MS, Lemey P, Faria NR, et al. Improving Bayesian population dynamics inference: a coalescent-based model for multiple loci. Mol Biol Evol, 2013; 30, 713−24. doi:  10.1093/molbev/mss265
    [28] Drummond AJ, Ho SY, Phillips MJ, et al. Relaxed phylogenetics and dating with confidence. PLoS Biol, 2006; 4, e88. doi:  10.1371/journal.pbio.0040088
    [29] Beyrer C, Razak MH, Lisam K, et al. Overland heroin trafficking routes and HIV-1 spread in south and south-east Asia. AIDS, 2000; 14, 75−83. doi:  10.1097/00002030-200001070-00009
    [30] Yang RG, Kusagawa S, Zhang CY, et al. Identification and characterization of a new class of human immunodeficiency virus type 1 recombinants comprised of two circulating recombinant forms, CRF07_BC and CRF08_BC, in China. J Virol, 2003; 77, 685−95. doi:  10.1128/JVI.77.1.685-695.2003
    [31] Bao YP, Liu ZM, Lu L. Review of HIV and HCV infection among drug users in China. Curr Opin Psychiatry, 2010; 23, 187−94. doi:  10.1097/YCO.0b013e328338658b
    [32] Feng Y, He X, Hsi JH, et al. The rapidly expanding CRF01_AE epidemic in China is driven by multiple lineages of HIV-1 viruses introduced in the 1990s. AIDS, 2013; 27, 1793−802. doi:  10.1097/QAD.0b013e328360db2d
    [33] Zhang J, Shen ZY, Li Z, et al. Genetic characteristics of CRF01_AE among newly diagnosed HIV-1-infected 16- to 25-year olds in 3 geographic regions of Guangxi, China. Medicine (Baltimore), 2015; 94, e894. doi:  10.1097/MD.0000000000000894
    [34] Chen Y, Shen ZY, Feng Y, et al. HIV-1 subtype diversity and transmission strain source among men who have sex with men in Guangxi, China. Sci Rep, 2021; 11, 8319. doi:  10.1038/s41598-021-87745-3
    [35] Baeten JM, Chohan B, Lavreys L, et al. HIV-1 subtype D infection is associated with faster disease progression than subtype A in spite of similar plasma HIV-1 loads. J Infect Dis, 2007; 195, 1177−80. doi:  10.1086/512682
    [36] Stephenson KE, Barouch DH. A global approach to HIV-1 vaccine development. Immunol Rev, 2013; 254, 295−304. doi:  10.1111/imr.12073
    [37] Hué S, Pillay D, Clewley JP, et al. Genetic analysis reveals the complex structure of HIV-1 transmission within defined risk groups. Proc Natl Acad Sci USA, 2005; 102, 4425−9. doi:  10.1073/pnas.0407534102
    [38] Gan MZ, Zheng S, Hao JJ, et al. The prevalence of CRF55_01B among HIV-1 strain and its connection with traffic development in China. Emerg Microbes Infect, 2021; 10, 256−65. doi:  10.1080/22221751.2021.1884004
    [39] Alirol E, Getaz L, Stoll B, et al. Urbanisation and infectious diseases in a globalised world. Lancet Infect Dis, 2011; 11, 131−41. doi:  10.1016/S1473-3099(10)70223-1
    [40] Hammett TM, Des Jarlais DC, Kling R, et al. Controlling HIV epidemics among injection drug users: eight years of Cross-Border HIV prevention interventions in Vietnam and China. PLoS One, 2012; 7, e43141. doi:  10.1371/journal.pone.0043141
    [41] Shao YM, Jia ZW. Challenges and opportunities for HIV/AIDS control in China. Lancet, 2012; 379, 804.
    [42] Region TPsGoGZA. Implementation plan of Guangxi AIDS conquering project (2010-2014). 2010.
  • 22260+Supplementary Materials.pdf
  • 加载中
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(6)  / Tables(6)

Article Metrics

Article views(650) PDF downloads(51) Cited by()

Proportional views
Related

HIV-1 Transmission among Injecting Drug Users is Principally Derived from Local Circulating Strains in Guangxi, China

doi: 10.3967/bes2023.051
Funds:  National Natural Science Foundation of China [grant no. 82060610]; Guangxi Scientific and Technological Key Project [Gui Ke AB19245038]; National Natural Science Foundation of China [grant no. 82103899]; Guangxi Scientific and Technological Key Project [Guike 2022AC23005, 2022AC20031,2022JJA141110]; and Science and Technology Project of Nanning [20223051]
  • Author Bio:

  • Corresponding author: LIANG Hao, Professor, PhD, Tel: 86-771-5323713, E-mail: lianghao@gxmu.edu.cn YE Li, Professor, PhD, Tel: 86-771-5344361, E-mail: yeli@gxmu.edu.cn LIANG Bing Yu; Professor, PhD, Tel: 86-771-5334215, E-mail: liangbingyu@gxmu.edu.cn
  • LIANG Hao, YE Li, and LIANG Bing Yu conceived and designed the study. YANG Shi Xiong and MO Ju Cong collected the specimen information. YANG Yuan, ZHANG Fei, HUANG Chun Yuan, YANG Yao, and LIANG Na performed experimental work. FENG Yi, HUANG Jie Gang, and NING Chuan Yi provided insight into the experimental design and data analysis. CEN Ping and LIANG Hua Yue performed data analysis and wrote the first draft. All authors approved the manuscript.
  • The authors declare that there are no conflicts of interest.
  • &These authors contributed equally to this work.

Abstract:   Objective   The mode of human immunodeficiency virus (HIV) transmission via injection drug use (IDU) still exists, and the recent shift in IDU-related transmission of HIV infection is largely unknown. The purpose of this study was to analyze the spatiotemporal sources and dynamics of HIV-1 transmission through IDU in Guangxi.   Methods   We performed a molecular epidemiological investigation of infections across Guangxi from 2009 to 2019. Phylogenetic and Bayesian time-geographic analyses of HIV-1 sequences were performed to confirm the characteristics of transmission between IDUs in combination with epidemiological data.   Results   Among the 535 subjects, CRF08_BC (57.4%), CRF01_AE (28.4%), and CRF07_BC (10.7%) were the top 3 HIV strains; 72.6% of infections were linked to other provinces in the transmission network; 93.6% of sequence-transmitted strains were locally endemic, with the rest coming from other provinces, predominantly Guangdong and Yunnan; 92.1% of the HIV transmission among people who inject drugs tended to be transmitted between HIV-positive IDUs.   Conclusion  HIV recombinants were high diversity, and circulating local strains were the transmission sources among IDUs in Guangxi. However, there were still cases of IDUs linked to other provinces. Coverage of traditional prevention strategies should be expanded, and inter-provincial collaboration between Guangxi, Yunnan, and Guangdong provinces should be strengthened.

LIANG Hao, YE Li, and LIANG Bing Yu conceived and designed the study. YANG Shi Xiong and MO Ju Cong collected the specimen information. YANG Yuan, ZHANG Fei, HUANG Chun Yuan, YANG Yao, and LIANG Na performed experimental work. FENG Yi, HUANG Jie Gang, and NING Chuan Yi provided insight into the experimental design and data analysis. CEN Ping and LIANG Hua Yue performed data analysis and wrote the first draft. All authors approved the manuscript.
The authors declare that there are no conflicts of interest.
&These authors contributed equally to this work.
CEN Ping, LIANG Hua Yue, YANG Yuan, ZHANG Fei, YANG Shi Xiong, MO Ju Cong, FENG Yi, HUANG Jie Gang, NING Chuan Yi, HUANG Chun Yuan, YANG Yao, LIANG Na, LIANG Bing Yu, YE Li, LIANG Hao. HIV-1 Transmission among Injecting Drug Users is Principally Derived from Local Circulating Strains in Guangxi, China[J]. Biomedical and Environmental Sciences, 2023, 36(5): 418-430. doi: 10.3967/bes2023.051
Citation: CEN Ping, LIANG Hua Yue, YANG Yuan, ZHANG Fei, YANG Shi Xiong, MO Ju Cong, FENG Yi, HUANG Jie Gang, NING Chuan Yi, HUANG Chun Yuan, YANG Yao, LIANG Na, LIANG Bing Yu, YE Li, LIANG Hao. HIV-1 Transmission among Injecting Drug Users is Principally Derived from Local Circulating Strains in Guangxi, China[J]. Biomedical and Environmental Sciences, 2023, 36(5): 418-430. doi: 10.3967/bes2023.051
    • Following China’s open-door policy, the use of illicit substances emerged, which was associated with global drug-trafficking activities in the 1980s[1]. The proximity of the “Golden Triangle” (the border areas of Thailand, Myanmar, and Laos) facilitates drug trafficking from the origin of the drug to Yunnan and Guangxi, two provinces in southwest China, or from Vietnam to Guangxi, and then to other neighboring provinces[2,3]. In 1996, the first HIV case in Guangxi was identified among people who inject drugs (PWID)[4]. Since then, Guangxi has become an HIV-1 hotspot, with many HIV-infected PWID[5]. Until 2006, IDUs were a major contributor to the spread of HIV in China, although sexual transmission was subsequently identified as the main route of HIV transmission[6]. From 2012 to 2019, 28.0% of newly diagnosed cases of HIV infected via IDU had baseline CD4 cell counts greater than or equal to 500 cells/μL, indicating that these individuals may have been infected within the past 2 years[7]. Since 2010, the HIV prevalence among new entrants tested for HIV at detoxification centers has been increasing[8], which means that HIV transmission via IDU still exists[7].

      Previous evidence has said that HIV is characterized by mutation and recombination, which leads to high genetic diversity and complex heterogeneity[9]. HIV-1 subtype B’ (Tai B) strains and C were co-circulated among IDU in Yunnan Province in the early 1990s[10-12]. Subsequently, CRF07_BC, CRF08_BC, and a large number of unique recombinant forms began to circulate widely among IDUs in China[13,14]. CRF01_AE was found in populations at risk of sexual transmission and IDUs in Yunnan and Guangxi[15-17], bordering Thailand, Myanmar, and Vietnam, and quickly became the most widespread HIV-1 strain in the national molecular epidemiological survey[6]. Guangxi’s HIV epidemic had become one of the worst in China as the initial epidemic was driven by IDUs in the 1990s. In Guangxi, HIV subtypes B′ (Tai B), C, D, and CRF01_AE were identified in 1996[17]. HIV subtype B/C recombinants were found among IDUs in Yunnan Province and Baise City, Guangxi[5]. During this period, the CRF01_AE strain was found among IDUs in Pingxiang and Nanning cities in Guangxi, which are located along the main drug-trafficking route from Vietnam to China[18]. CRF07_BC and CRF08_BC, originally formed in Yunnan Province, were first detected among IDUs in Guangxi in 1997[19]. In contrast, the first CRF01_AE strain identified among IDUs in Guangxi had a significant clustering with strains found in Northern Vietnam[20]. Recent studies have found that HIV dominant subtypes among IDUs in China are CRF07-BC, followed by CRF01-AE and CRF08-BC[21].

      In Guangxi, HIV-1 transmission via IDU has not yet been eliminated, and the epidemic is characterized by substantial viral diversity and distinct molecular properties. However, there have been few studies on the regional HIV epidemic among PWIDs in recent years. A detailed and comprehensive analysis of the geographic distribution of viral genotypes is needed. In this study, to determine the source of infections, we provide evidence of the genotypes and the temporal-geographical origins of HIV transmission among IDUs in Guangxi. Furthermore, we spatialized the HIV-1 genetic transmission network, with an explanation of spatial distribution and spatial clustering patterns of phylogenetic clustering rates and spatial links of HIV-1 among cities, which helps to better understand the dynamic characteristics toward the IDU-related transmission of HIV-1 subtypes at the regional level and may contribute to developing strategies to reduce the risk of HIV transmission via IDU.

    • The blood samples from 535 confirmed HIV-1-positive IDUs were continuously collected between January 2009 and December 2019 through voluntary counseling and testing sites in Guangxi, China. We performed a BLAST search using the sequences and selected the top five sequences with the highest homology from publicly available HIV databases at Los Alamos National Laboratory[22]. Seven hundred fifty reference sequences from the HIV database were selected for analysis after excluding repeated sequences (accessed May 2021). Overall, 1,285 sequences were included for subsequent analyses.

    • Sequencher 5.1 software was used to edit all sequences, and BioEdit 7.1 software was used to align them. The reference sequences were obtained from the Los Alamos National Laboratory database and covered the major HIV-1 subtypes and circulating recombinant forms (CRFs). Phylogenetic tree analyses were performed using the neighbor-joining method based on the Kimura2-parameter model with 1,000 bootstrap replicates by using MEGA7.0 software[23].

    • One thousand two hundred eighty-five aligned sequences, including reference and subjects’ sequences, were entered into the HyPhy software to calculate genetic distances, and Tamura-Nei93 pairwise genetic distances were calculated for all pairs of sequences. Two sequences with a genetic distance of 0.6% were identified as potential transmission partners. Cytoscape V3.5.1 was used to process the network data and construct the transmission network. As previously mentioned, the characteristics of the network, including nodes (individuals in the network), edges (the link between two nodes, representing the potential transmission relationship between the two individuals), degrees (the number of edges linking one node to other nodes), network sizes (the number of individuals in a cluster), and clusters (groups of linked sequences), were determined[24]. The Cytoscape V3.5.1 software was developed by Institute for Systems Biology. The network data were visualized using ArcGIS Desktop 10.7, Version: 10.7.0.10450 (Environmental Systems Research Institute)[25].

    • For accurate analysis and better display of phylogenetic trees, 516 pol sequences were divided into three data sets according to genotypes CRF01_AE, CRF07_BC, and CRF08_BC, with 152, 57, and 307 sequences, respectively. Since only one sequence belonged to subtype C and CRF57_BC, no phylogenetic tree was required to determine the transmission source. Furthermore, additional 666 unique pol sequences identified by BLAST as CRF01_AE, CRF07_BC, and CRF08_BC, were divided into data sets with 368, 134, and 164 sequences, respectively. Thus, 1,182 (516 + 666) pol sequences were used to determine the sources of HIV-1 transmission strains.

      To reconstruct the spatial dynamics and estimate the sources of the strains, a Bayesian discrete phylogeographic approach was performed using Markov chain Monte Carlo (MCMC) runs of 300 million generations with BEAST v.1.8.4 under a Bayesian skygrid demographic model. The first 10%–30% of the states from each run were discarded as the burn-in[26,27]. All three data sets were analyzed using a general time-reversible model specifying a gamma distribution based on a neighbor-joining tree as a prior on each relative substitution rate and a relaxed uncorrelated lognormal molecular clock model to infer the timescale of HIV evolution with a gamma distribution prior on the mean clock rate (shape = 0.001, scale = 1,000)[28]. The Bayesian MCMC output was analyzed using Tracer v1.6 (http://beast.bio.ed.ac.uk/Tracer). Maximum clade credibility (MCC) trees were generated in TreeAnnotator v1.8.4. Finally, the most probable temporal, geographical origin, and risk category of HIV transmission strains for the study participants were inferred according to the output of the posterior of Bayesian estimation and visualized on MCC trees using the FigTree software v1.4.3 (http://beast.bio.ed.ac.uk). The tree figure editing was conducted using Adobe Illustrator 2020 (V 24.0.1.341).

      Sequences with a posterior probability over 0.7 were defined as local circulating strains and risk categories in Guangxi, and the posterior probability refers to the support of locations at the nodes. The provinces of origin of HIV transmission strains in our study samples were then tabulated according to different HIV subtypes.

    • The edges between the two nodes were calculated to explore factors associated with potential transmission. Comparisons between individuals were based on the following factors: year of diagnosis, residence region, age, marital status, educational level, ethnicity, HIV subtype, and time of infection. Significant differences in categorical variables were analyzed using chi-square or Fisher’s exact tests. P < 0.05 was considered statistically significant.

    • Table 1 provides a summary of the 535 HIV-positive IDUs samples. The average age of these individuals was 42.53 ± 8.13 years, of which 66.7% (357/535) were over 35 years old. There were 489 (91.4%) women and 46 (8.6%) men. Of these IDUs, 38.9% (208/535) were permanent residents in Nanning, 33.6% (327/535) were from Qinzhou, and 27.5% were from other cities in Guangxi. The ethnicities of the participants were mainly Han (411, 76.8%) and Zhuang (105, 19.6%), and the education level was mainly junior school or below (484, 90.5%). In terms of marital status, 61.3% (328/535) were single, and 35.7% (191/535) were married. For CD4+ T-cell count, 110 (20.6%) individuals had a CD4+ T-cell count ≥ 500 cells/μL, 163 (30.5%) had a CD4+ T-cell count < 250 cells/μL, and 194 (36.3%) between 250 and 500 cells/μL. Of all participants, 320 subjects (59.8%) were ART naïve (Table 1).

      VariablesOverall01AE07BC08BCOther subtypescP-value*
      Overall5351525730719
      Age group (years), n (%)
       ≤ 253 (0.6)1 (0.7)0 (0.0)2 (0.7)0< 0.001
       26–35114 (21.3)20 (13.2)6 (10.5)86 (28.0)2 (10.5)
       > 36357 (66.7)112 (73.6)50 (87.7)182 (59.2)13 (68.4)
       Unknown61 (11.4)19 (12.5)1 (1.8)37 (12.1)4 (21.1)
      Gender, n (%)
       Female46 (8.6)17 (11.2)4 (7.0)24 (7.8)1 (5.3)0.576
       Male489 (91.4)135 (88.8)53 (93.0)283 (92.2)18 (94.7)
      Marital status, n (%)
       Married191 (35.7)58 (38.2)17 (29.8)110 (35.8)6 (31.6)0.639
       Single328 (61.3)87 (57.2)39 (68.4)189 (61.6)13 (68.4)
       Unknown16 (3.0)7 (4.6)1 (1.8)8 (2.6)0
      Region, n (%)
       Nanning208 (38.9)65 (42.7)4 (7.0)130 (42.4)9 (47.4)< 0.001
       Liuzhou98 (18.3)57 (37.5)40 (70.1)1 (0.3)0 (0.0)
       Qinzhou180 (33.6)22 (14.5)12 (21.1)136 (44.3)10 (52.6)
       Othera49 (9.2)8 (5.3)1 (1.8)40 (13.0)0
      Ethnicity, n (%)
       Han411 (76.8)106 (69.8)47 (82.4)240 (78.2)18 (94.7)0.28
       Zhuang105 (19.6)37 (24.3)9 (15.8)58 (18.9)1 (5.3)
       Otherb3 (0.6)2 (1.3)0 (0.0)1 (0.3)0
       Unknown16 (3.0)7 (4.6)1 (1.8)8 (2.6)0
      Education, n (%)
       Junior school or less484 (90.5)129 (84.9)46 (80.7)291 (94.8)18 (94.7)0.001
       High school34 (6.3)14 (9.2)10 (17.5)9 (2.9)1 (5.3)
       College or higher2 (0.4)2 (1.3)0 0 0
       Unknown15 (2.8)7 (4.6)1 (1.8)7 (2.3)0
      HIV diagnosed year, n (%)
       2009–2016274 (51.2)67 (44.1)45 (78.9)155 (50.5)7 (36.8)< 0.001
       2017–2019261 (48.8)85 (55.9)12 (21.1)152 (49.5)12 (63.2)
      CD4+ (cells/mm3), n (%)
       < 250163 (30.5)59 (38.8)18 (31.6)80 (26.1)6 (31.6)0.179
       250–500194 (36.3)43 (28.3)24 (42.1)119 (38.8)8 (42.1)
       > 500110 (20.6)29 (19.1)11 (19.3)68 (22.1)2 (10.5)
       Unknown68 (12.6)21 (13.8)4 (7.0)40 (13.0)3 (15.8)
      Treatment, n (%)
       Treat215 (40.2)94 (61.8)27 (47.4)91 (29.6)3 (15.8)< 0.001
       Untreated320 (59.8)58 (38.2)30 (52.6)216 (70.4)16 (84.2) 
        Note. Unknown data are not available. aOther regions consist of Baise, Chongzuo, and Laibin. bOther minority ethnicity consists of Yao, Hui, and Shui. cOther subtypes consist of C (n = 1), CRF_57BC (n = 1), and unique recombinant forms (n = 17). *Chi-square trend test.

      Table 1.  Demographic and clinical characteristics of the study population

    • As shown in Figure 1, 535 samples were sequenced and genotyped. HIV-1 subtypes were mainly CRF08_BC, CRF01_AE, and CRF07_BC. A few other subtypes/CRFs were C and CRF57_BC. The 152 CRF01_AE sequences were divided into four clusters: CRF01_AE Cluster 1 (82/152, 53.9%), CRF01_AE Cluster 2 (33/152, 21.7%), CRF01_AE Cluster 3 (3/152, 2.0%), and CRF01_AE Cluster 4 (2/152, 1.3%) (Figure 1, Table 2). Seventeen (3.2%) URFs were discovered, mainly between CRF_01AE and CRF_07/08_BC (Supplementary Figure S1, Supplementary Table S1, available in www.besjournal.com). There were significant differences in the distribution of genotypes by the time of HIV diagnosis, participants’ region, age, education, and treatment (Table 1).

      Figure 1.  Maximum-likelihood trees of 535 HIV-1 pol sequences sampled from IDUs in Guangxi, China. The stability of each node was assessed by bootstrap analyses with 1,000 replicates. The red spots and green spots indicate the sequences of this study and references from the HIV database, respectively.

      SubtypeThis studyHIV databaseTotal
      CRF01_AE152368519
      CRF07_BC57134191
      CRF08_BC307164471
      CRF57_BC101
      CRF52_01B022
      CRF85_BC011
      A1088
      C11213
      URF176178
      Total5357501,285

      Table 2.  List of HIV-1 sequences used in the study

    • At a genetic distance threshold of 0.6%, transmission network analysis found 34.3% (441/1,285) genetically linked individuals, forming 128 clusters and ranging in size from 2 to 80. The largest cluster consisting of 595 links was comprised predominantly of subtype CRF08_BC, which clustered more easily than other subtypes, indicating that CRF08_BC tended to group into a large cluster. In the network, 121 subjects (27.4%) were from Guangxi, and 320 (72.6%) were from other provinces (Figure 2). To explore factors associated with the potential transmission, 535 Guangxi HIV-positive IDUs samples were included in the analysis. The chi-square test revealed significant differences in the transmission network between subtypes (P = 0.005), age (P < 0.001), region (P < 0.001), year of diagnosis (P < 0.001), and treatment status (P < 0.001) (Table 3).

      Figure 2.  Region and transmission route associated with HIV-1 transmission clusters. The network was constructed using R script with the network package in R version 4.0.2 software. Shape codes the reported transmission route, and colors denote regions for each individual. HET, heterosexual; IDU, intravenous drug user; MSM, men who have sex with men.

      VariablesOverallClusteredNon-clusteredP-value*
      Overall535142393
      Subtype, n (%)
       01AE152 (28.4)31 (21.8)121 (30.8)0.005
       07BC57 (10.6)10 (7.0)47 (12.0)
       08BC307 (57.4)99 (69.6)208 (52.9)
       Otherc19 (3.6)2 (1.4)17 (4.3)
      Age group (years), n (%)
       ≤ 253 (0.6)2 (1.4)1 (0.3)< 0.001
       26–35114 (21.3)51 (35.9)63 (16.0)
       > 36357 (66.7)78 (54.9)279 (71.0)
       Unknown61 (11.4)11 (7.8)50 (12.7)
      Gender, n (%)
       Female46 (8.6)13 (9.2)33 (8.4)0.919
       Male489 (91.4)129 (90.8)360 (91.6)
      Marital status, n (%)
       Married191 (35.7)41 (28.9)150 (38.2)0.041
       Single328 (61.3)99 (69.7)229 (58.2)
       Unknown16 (3.0)2 (1.4)14 (3.6)
      Region, n (%)
       Nanning208 (38.9)43 (30.4)165 (42.0)< 0.001
       Liuzhou98 (18.3)11 (7.7)87 (22.1)
       Qinzhou180 (33.6)81 (57.0)99 (25.2)
       Othera49 (9.2)7 (4.9)42 (10.7)
      Ethnicity, n (%)
       Han411 (76.8)119 (83.8)292 (74.3)0.131
       Zhuang105 (19.6)19 (13.4)86 (21.9)
       Otherb3 (0.6)1 (0.7)2 (0.5)
       Unknown16 (3.0)3 (2.1)13 (3.3)
      Education, n (%)
       Junior school or less484 (90.4)133 (93.7)351 (89.3)0.407
       High school34 (6.4)7 (4.9)27 (6.9)
       College or higher2 (0.4)0 (0.0)2 (0.5)
       Unknown15 (2.8)2 (1.4)13 (3.3)
      HIV diagnosed year, n (%)
       2009–2016274 (51.2)91 (64.1)183 (46.6)< 0.001
       2017–2019261 (48.8)51 (35.9)210 (53.4)
      CD4+ (cells/μL), n (%)
       < 250163 (30.5)42 (29.5)121 (30.8)0.059
       250–500194 (36.2)49 (34.5)145 (36.9)
       > 500110 (20.6)39 (27.5)71 (18.1)
       Unknown68 (12.7)12 (8.5)56 (14.2)
      Treatment, n (%)
       Treat215 (40.2)38 (26.8)177 (45.0)< 0.001
       Untreated320 (59.8)104 (73.2)216 (55.0) 
        Note. Unknown data are not available. aOther regions consist of Baise, Chongzuo, and Laibin. bOther minority ethnicity consists of Yao, Hui, and Shui. cOther subtypes consist of C (n = 1) and unique recombinant forms (n = 18). *Chi-square trend test.

      Table 3.  Factors associated with the probability of belonging to a transmission cluster of the study population

    • To explore the geographic dimension of HIV-1 transmission among IDUs, we constructed transmission routes using the R language. Guangxi was linked with Vietnam (2.7%) and 11 provinces in China, including Guangdong (43.5%), Yunnan (16.7%), and others (37.1%). There were 99 (22.3%) intercity links from Guangxi to other provinces, with Guangdong (44, 9.9%) having the strongest link. In addition, 248 (55.8%) intercity links were from other provinces to Guangxi, of which Guangdong (128, 28.8%) had the strongest link, followed by Yunnan (69, 15.5%) (Figure 3).

      Figure 3.  The geographic dimension of HIV-1 transmission. HIV-1 transmission between Guangxi and other provinces. The line width indicates the number of individuals.

    • Considering that CRF08_BC, CRF01_AE, and CRF07_BC were the dominant strains, other subtypes were excluded from the analysis due to small sample sizes (n < 5).

      The sources of HIV-1 transmission strains in recently infected individuals in Guangxi were estimated by the reconstruction of Bayesian discrete phylogeographic approaches under a Bayesian skygrid demographic model. Based on the criteria of posterior probability over 0.7, the transmission source of 516 samples could be determined. As shown in Table 4, 483 (93.6%) samples were derived from strains already circulating in Guangxi, with the remaining 33 (6.4%) derived from strains circulating outside Guangxi. HIV transmission sources varied by subtypes and education levels. Subtype CRF01_AE had the most diversified source (Table 5).

      Posterior probabilityProvinceTotalCRF01_AECRF07_BCCRF08_BC
      > 0.7aGuangxi48313353297
      Guangdong191027
      Zhejiang4310
      Yunnan3111
      Anhui2200
      Hongkong, China2002
      Beijing1100
      Shenzhen1100
      Philippines1100
      Total51615257307
        Note. aSequences with posterior probability over 0.7 were defined as Guangxi local circulating strains.

      Table 4.  Distribution of HIV transmission strains circulating provinces

      VariablesOverallGuangxiOutside GuangxiP-value*
      Overall51648333
      Subtype, n (%)
       01AE152 (29.5)133 (27.5)19 (57.5)< 0.001
       07BC57 (11.0)48 (9.9)9 (27.3)
       08BC307 (59.5)302 (62.6)5 (15.2)
      Age group (years), n (%)
       ≤ 253 (0.6)3 (0.6)0 0.815
       26–35112 (21.7)106 (21.9)6 (18.2)
       > 36344 (66.7)322 (66.7)22 (66.6)
       Unknown57 (11.0)52 (10.8)5 (15.2)
      Gender, n (%)
       Female45 (8.7)41 (8.5)4 (12.1)0.692
       Male471 (91.3)442 (91.5)29 (87.9)
      Marital status, n (%)
       Married185 (35.9)171 (35.4)14 (42.4)0.716
       Single315 (61.0)297 (61.5)18 (54.6)
       Unknown16 (3.1)15 (3.1)1 (3.0)
      Region, n (%)
       Nanning199 (38.6)186 (38.5)13 (39.4)0.62
       Liuzhou98 (19.0)91 (18.8)7 (21.2)
       Qinzhou170 (32.9)158 (32.8)12 (36.4)
       Othera49 (9.5)48 (9.9)1 (3.0)
      Ethnicity, n (%)
       Han393 (76.1)366 (75.8)27 (81.9)0.178
       Zhuang104 (20.2)100 (20.7)4 (12.1)
       Otherb3 (0.6)2 (0.4)1 (3.0)
       Unknown16 (3.1)15 (3.1)1 (3.0)
      Education, n (%)
       Junior school or less466 (90.3)441 (91.3)25 (75.7)< 0.001
       High school33 (6.4)28 (5.8)5 (15.2)
       College or higher2 (0.4)0 2 (6.1)
       Unknown15 (2.9)14 (2.9)1 (3.0)
      HIV diagnosed year, n (%)
       2009–2016267 (51.7)248 (51.3)19 (57.6)0.608
       2017–2019249 (48.3)235 (48.7)14 (42.4)
      CD4+ (cells/μL), n (%)
       < 250157 (30.5)144 (29.9)13 (39.4)0.329
       250–500186 (36.0)176 (36.4)10 (30.3)
       > 500108 (20.9)104 (21.5)4 (12.1)
       Unknown65 (12.6)59 (12.2)6 (18.2)
      Treatment, n (%) 
       Treat212 (41.1)201 (41.6)11 (33.3)0.452
       Untreated304 (58.9)282 (58.4)22 (66.7)
        Note. Unknown data are not available. aOther regions consist of Baise, Chongzuo, and Laibin. bOther minority ethnicity consists of Yao, Hui, and Shui. *Chi-square trend test.

      Table 5.  Factors associated with the probability of Guangxi local transmission source of the study population

      The introduction time for CRF01_AE strains in Guangxi IDUs with effective population size circulating may be 1994.25 (1992.97–1995.58), from Vietnam but originated in Thailand. Two distinct phylogenetic clusters of CRF01_AE were identified. Cluster strains were placed within a large phylogenetic group consisting of CRF01_AE strains that originated in Thailand and then spread to Guangxi through Vietnam. This cluster contained CRF01_AE strains circulating among IDUs in Guangxi, Guangdong, and Yunnan. Another cluster originated in Thailand and then spread to Guangxi through Guangdong and formed three clusters of Guangxi local circulating strains. CRF07_BC strains probably came from Guangdong in 1995.23 (1993.64–1996.64). CRF07_BC strains in Guangxi may originate in Yunnan and then spread to Guangxi through Guangdong. It is worth noting that some individuals may be associated with the epidemic among MSMs in Guangdong. CRF08_BC strains, mainly directly originating from Yunnan, probably spread to Guangxi in 1992.29 (1989.05–1995.36), forming large local clusters among IDUs (Supplementary Figure S2, available in www.besjournal.com).

      Figure S2.  This figure illustrated time-scaled phylogeographic history of three HIV-1 subtypes (CRF01_AE, CRF08_BC, CRF07_BC) for IDU in Guangxi. Background colours around branches represented the most probable origin province of the parental node of each branch. The names of the origin provinces were indicated by the colour of the legend in the upper left corner. The values of probability inferred by Bayesian skygrid demographic model estimation for the most probable ancestral location were specified at the corresponding nodes. The probable time of introduction for these three key subtypes with an effective population size circulating in Guangxi IDU was indicated at key nodes. Red and blue dots at the ends of branches denoted the sequences of study participants and references, respectively. The Bayesian MCC trees were visualised through the FigTree software v1.4.3 (http://beast.bio.ed.ac.uk). The figure editing was using Adobe Illustrator 2020 (V 24.0.1.341).

    • HIV transmission risk categories for 516 sequences with CRF08_BC, CRF01_AE, and CRF07_BC were visualized on the MCC tree based on a posterior probability criterion exceeding 0.7 (Supplementary Figure S2). Of these, 92.1% of sequences were more frequently linked to others reporting IDU risk. The percentage of CRF08_BC sequences linked to IDU risk was as high as 97.7%, followed by CRF07_BC (94.7%) and CRF01_AE (80.1%). In addition, some were also linked to reporting heterosexual (7.2%) and MSM (0.6%) risks. Compared with the CRF08_BC and CRF07_BC groups, the CRF01_AE group had a higher proportion of IDUs (18.5%) linked to those reporting heterosexual risk (Supplementary Figure S3, available in www.besjournal.com).

      Figure S3.  This figure illustrated HIV transmission risk category link of three HIV-1 subtypes (CRF01_AE, CRF07_BC, CRF08_BC) for IDU in Guangxi. The transmission risk category of the strains is color-coded as shown in the inset. Red and blue dots at the ends of branches denoted sequences of study participants and references, respectively. Background colors around branches represent the most probable origin of the parental node of each branch. The Bayesian MCC trees were visualized through the FigTree software v1.4.3 (http://beast.bio.ed.ac.uk). The figure editing was using Adobe Illustrator 2020 (V 24.0.1.341).

    • This study investigated the source and genetic transmission patterns of HIV-1 strains circulating among IDUs in Guangxi using molecular epidemiological approaches. Although a previous study discovered that HIV-1 transmission in Guangxi predominantly occurs from strains already circulating locally, the small sample size in a previous study failed to indicate the current HIV/AIDS status of IDUs in Guangxi[23]. In recent years, HIV infection due to drug use still disproportionately affects the AIDS epidemic, and the transmission patterns of this high-risk group still need to be researched. To the best of our knowledge, this is the first large-scale study of HIV genetic transmission among drug users in Guangxi using molecular epidemiological methods. We found that HIV-1 transmission among recently diagnosed IDUs was predominantly derived from locally circulating strains and formed epidemic clusters in Guangxi. The major HIV-1 subtype strains among IDUs in Guangxi were CRF08_BC, CRF01_AE, and CRF07_BC, which came from different regions, respectively, rather than a single source. The most probable source of CRF01_AE strains was either directly from Vietnam to Guangxi in 1994.25 (1992.97–1995.58) or indirectly from Thailand to Guangxi through Guangdong and then spread as a single lineage throughout the region. This result is in agreement with previous studies, suggesting that cross-border transmission of CRF01_AE between Vietnam and Guangxi occurs[16]. In contrast, CRF08_BC strains, mainly originated in Yunnan, spread to Guangxi in 1992.29 (1989.05–1995.36) and formed the main epidemic strain among IDUs. CRF07_BC strains may originate in Yunnan and spread to Guangxi via Guangdong in 1995.23 (1993.64–1996.64). In contrast to previous studies[18,29,30], our findings showed the source provinces of the majority of HIV-1 strains, except CRF08_BC, in recently infected IDUs were Guangdong, although we also observed some strains from Vietnam and Yunnan. Changes in the source provinces of these strains illustrate a shift in HIV-1 transmission drivers. They may have contributed to a new emergence of CRFs and URFs of HIV-1 and increased the complexity of the HIV epidemic among IDUs in Guangxi.

      We identified three major CRFs (CRF08_BC, CRF01_AE, and CRF07_BC) in this study, which are consistent with previous results among IDUs in southwest provinces in China[21] but were different from those in the Southeast coast and the Central Plains in China[6,31]. Our study further confirmed that CRF08_BC is the dominant subtype present in IDUs in Yunnan and Guangxi and is rarely distributed in other places in the country. A previous study in China categorized CRF01_AE into seven independent lineages, namely, CRF01_AE Clusters 1–7[32]. Our study identified four lineages (Clusters 1–4)[23] among CRF01_AE strains circulating in Guangxi, but no strains from CRF01_AE Clusters 5–7 were found in this study. In addition, CRF01_AE Cluster 1 accounted for more than half of the total CRF01_AE strains, which is consistent with the findings of some recent studies[33]. In our study, we observed a low prevalence of CRF07_BC strains and a high prevalence of CRF08_BC strains among IDUs in Guangxi, which is consistent with a systematic review reporting that CRF08_BC is the most prevalent genotype in the southwest region of China[21]. Notably, CRF07_BC may be experiencing an increasing epidemic trend, especially among MSMs in Guangxi[34], and we need to be alert to the possibility of a new wave of HIV infection in IDUs and MSMs driven by subtype CRF07_BC strains. Furthermore, we found 17 recombinant subtypes formed of inter-subtypes and inter-CRF, which was higher than those in previous studies[6] and increased the complexity of HIV molecular patterns. This is also a major issue of a public health concern as novel strains may increase transmission, promote disease progression, or change resistance patterns to antiretroviral drugs[35,36].

      In the networks, we found that two large networks containing more than 100 IDUs infected with CRF01_AE or CRF_08BC were linked to heterosexuals from outside Guangxi and abroad, supporting the importance of identifying and controlling this “bridge population” to curb HIV transmission[37]. The geographic dimension analysis (Figure 3) indicated that Guangdong and Yunnan were the major provinces for HIV-1 cross-regional transmission in Guangxi, accounting for 60.2% of all cases. High-speed railway and aviation may have facilitated the migration between these regions and promoted HIV-1 cross-regional transmission[38]. Of course, there may be more fundamental social factors associated with the spatial mobility of HIV-1. As revealed by a previous study, the heterogeneity in the health status of urban dwellers could increase exposure and mobility and result in a high risk of disease transmission in large urban populations[39]. Cities become incubators and provide conditions for outbreaks, and future migration and urbanization may strengthen disease transmission between cities[39]. We also note that there was still an isolated cross-border transmission link to Vietnam. Although cross-border interventions have been effective in controlling the HIV epidemic among IDUs in recent years[40], the risk of cross-border transmission still needs to be considered.

      Although three common CRF strains in recently infected IDUs in Guangxi could be traced back to Vietnam, Guangdong, or Yunnan, the epidemic of HIV in Guangxi is mainly driven by local transmission. Specifically, in our study, more than 93% of HIV infections in the past 10 years were derived from HIV-1 strains already circulating in local Guangxi. Three common HIV CRF strains were most likely introduced to Guangxi IDUs in the early 1990s (1992–1995). Nonetheless, 6.4% IDUs were observed to be infected with HIV-1 from non-local circulating strains, which is comparable with the data for heterosexuals in Guangxi (5.9%)[23]. Non-local circulating strains were mainly transmitted from Guangdong, consistent with our findings in geographic dimension analysis. Two putative reasons could explain the shift in HIV transmission sources among IDUs found in this study. First, in the mid-1990s, Guangxi experienced a high HIV prevalence among IDUs, which was derived from two emerging epidemics: one from Vietnam and another from Yunnan Province[5]. In 2002, a Cross-Border Project for HIV prevention was well implemented among IDUs on both sides of the international border between Northern Vietnam and southern China[40], playing an important role in reducing HIV epidemics for nearly a decade. Second, the high population mobility between Guangxi and Guangdong and high urbanized migration in Guangxi may contribute to the shift in HIV-1 transmission source. Accordingly, this study highlights the importance of understanding the role of external-province population mobility and intracity urbanized migration in HIV transmission.

      An interesting finding of our study is that the primary HIV-1 transmission risk population among IDUs tended to be limited to PWID in recent years in Guangxi. Fortunately, these changes have not yet caused any rebound in HIV prevalence or incidence among IDUs. In fact, Guangxi suffered a massive HIV-1 epidemic among PWID between 1996 and 2003[4]. Following the PWID epidemic in the 1990s, the predominant mode of HIV transmission shifted to heterosexuality through the bridge role of FSWs in HIV transmission from IDUs to the general population[41]. In response to the HIV epidemic, the Guangxi government launched the Guangxi AIDS Conquering Project (GACP) in 2010[42]. From 2009 to 2015, HIV/AIDS prevention efforts in Guangxi has been significantly strengthened, the scope of prevention services, MMT or needle exchange, peer education, and HIV testing has been expanded, and the awareness of HIV/AIDS knowledge among at-risk groups has been improved over time. These structural interventions have substantially altered HIV-related risk behavior and reduced HIV morbidity. Data from this study indicate that the local GACP initiative combined with the National AIDS Control Policy helps to stabilize HIV incidence in Guangxi.

      Our study has several limitations. First, our study had limited samples and a convenience sampling method, both of which may lead to bias. Nevertheless, to the best of our knowledge, this study contains the largest number of sequences among IDUs in Guangxi to date to analyze the source of HIV strain transmission. Second, the missing variables of some classification indicators may have a potential impact on the analysis results. Third, most of the sequences from the other provinces were downloaded from the HIV database, as not all individuals living with HIV in other provinces have had their HIV-1 sequenced and submitted to the database. Last, we cannot identify the direct source of each transmission because there may be an intermediary between a case and its source. However, our analysis could identify closely linked infections.

      In conclusion, we found a high diversity of HIV-1 strains among IDUs in Guangxi over the past decade and determined that HIV-1 transmission among IDUs in Guangxi mainly came from local strains that were already circulating. In addition to implementing traditional HIV prevention models, it is necessary to continue to strengthen the coverage of NSP and MMT. Furthermore, we have identified complicated sources of HIV transmission in Guangxi, suggesting appropriate measures to limit the spread of HIV along illicit drug-trafficking routes in the region.

    • The individual participant-level data that underlie the results reported in this article will be shared after de-identification (text, tables, figures, and appendices).

    • Sequence nameEthnicitySexFragment start positionBreakpoint intervalFragment end positionFragment subtype
      2015QZ0137HanMale22532253-2713271308_BC
      27132713-3309330907_BC
      2019NN0153ZhuangMale22532253-2663266308_BC
      26632663-2973297301_AE
      29732973-3309330908_BC
      2016NND0016HanFemale22532253-2513251301_AE
      25132513-3053305307_BC
      30533053-3309330901_AE
      2018NN0035HanMale22532253-2793279308_BC
      27932793-3023302307_BC
      30233023-3309330908_BC
      2016NND0037HanMale22532253-2542254201_AE
      25422542-3034303407_BC
      30343034-3309330901_AE
      2018NN0038HanMale22532253-2543254301_AE
      25432543-3033303307_BC
      30333033-3309330901_AE
      2018NN0054HanMale22532253-2633263301_AE
      26332633-3309330907_BC
      2016NND0057HanMale22532253-2993299301_AE
      29932993-3309330908_BC
      2017NND0087HanMale22532253-2703270308_BC
      27032703-3309330907_BC
      2018QZ0437HanMale22532253-2693269308_BC
      26932693-3309330907_BC
      2015QZ0187HanMale22532253-2933293357_BC
      29332933-3309330907_BC
      2015QZ0197HanMale22532253-2723272308_BC
      27232723-3309330907_BC
      2015QZ0253HanMale22532253-2953295307_BC
      29532953-3309330901_AE
      2018QZ0564HanMale22532253-2713271308_BC
      27132713-3309330907_BC
      2018QZ0567HanMale22532253-2723272308_BC
      27232723-3309330907_BC
      2018QZ0571HanMale22532253-2728272808_BC
      27282728-3309330907_BC
      2018QZ0627HanMale22532253-2653265301_AE
      26532653-3309330907_BC

      Table S1.  Detail information of HIV-1 unique recombinant forms of 17 sample from 535 IDUs in Guangxi, China

      Figure S1.  Recombination patterns and parental origin in 1.2 kb pol region of 17 sample from 535 IDUs in Guangxi, China. The genome map of the pol gene DNA sequence of 17 recombinant sample was generated using the Online Recombinant HIV-1 Drawing Tool (https://www.hiv.lanl.gov/content/sequence/DRAW_CRF recom_mapper.html).

Reference (42)
Supplements:
22260+Supplementary Materials.pdf

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return