Consistency of spatial dynamics of HIV-1 and HCV among HIV-1/HCV coinfected drug users in China
BMC Infectious Diseases volume 21, Article number: 1001 (2021)
As the transmission routes of human immunodeficiency virus type 1 (HIV-1) and hepatitis C virus (HCV) are similar, previous studies based on separate research on HIV-1 and HCV assumed a similar transmission pattern. However, few studies have focused on the possible correlation of the spatial dynamics of HIV-1 and HCV among HIV-1/HCV coinfected patients.
A total of 310 HIV-1/HCV coinfected drug users were recruited in Yingjiang and Kaiyuan prefectures, Yunnan Province, China. HIV-1 env, p17, pol and HCV C/E2, NS5B fragments were amplified and sequenced from serum samples. The genetic characteristics and spatial dynamics of HIV-1 and HCV were explored by phylogenetic, bootscanning, and phylogeographic analyses.
Among HIV-1/HCV coinfected drug users, eight HCV subtypes (1a, 1b, 3a, 3b, 6a, 6n, 6v, and 6u) and two HIV-1 subtypes (subtype B and subtype C), three HIV-1 circulating recombinant forms (CRF01_AE, CRF07_BC and CRF08_BC), and four unique recombinant forms (URF_BC, URF_01B, URF_01C and URF_01BC) were identified. HCV subtype 3b was the most predominant subtype in both Yingjiang and Kaiyuan prefectures. The dominant circulating HIV-1 subtypes for drug users among the two areas were CRF08_BC and URF_BC. Maximum clade credibility trees revealed that both HIV-1 and HCV were transmitted from Yingjiang to Kaiyuan.
The spatial dynamics of HIV-1 and HCV among HIV-1/HCV coinfected drug users seem to have high consistency, providing theoretical evidence for the prevention of HIV-1 and HCV simultaneously.
At the end of 2019, the total numbers of people living with human immunodeficiency virus type 1 (HIV-1) and hepatitis C virus (HCV) worldwide were estimated to be 38 and 71 million, respectively (http://www.who.int). As both HIV-1 and HCV can transmit through blood, sexual contact and mother to child, HIV-1/HCV coinfection is common worldwide. A global systematic analysis showed that in HIV-infected individuals, HIV-1/HCV coinfection was 82.4% in people who injected drugs (PWIDs) . In China, the prevalence of HIV-1/HCV coinfection among PWIDs varied from 3.4% to 21.4% . Previous studies have found that the prevalence of HIV-1/HCV coinfection is 15% among PWIDs in Yunnan province, and a high proportion can be seen in Yingjiang and Kaiyuan prefectures .
In China, current circulating HIV-1 subtypes comprise subtype B, subtype C, circulating recombinant form (CRF) 01_AE, CRF07_BC, CRF08_BC and other CRFs and unique recombinant forms (URFs). HIV-1 CRF01_AE is mainly transmitted through sexual contact and originates from Thailand, while HIV-1 subtypes C and B are mainly transmitted through injection of drugs and originate from India and Thailand, respectively. The three subtypes all spread to Yunnan province and further to other provinces [4, 5]. HIV-1 CRF07_BC and CRF08_BC were first identified among drug users in Yunnan province , and the two HIV-1 recombinants spread rapidly among drug users in China. The most prevalent HCV isolates in China are genotypes 1, 2, 3 and 6 [7,8,9]. Different HCV subtypes had diverse transmission routes, such as subtypes 3b, 6n, and 6u originating from Yunnan and Guangxi and subtypes 2a and 6a originating from southern China (e.g., Guangdong), which further spread to Jiangsu . Zhang and colleagues found that HCV 3a spread from Xinjiang to Jiangsu and from Yunnan to Guangxi; these routes were similar to the routes of HIV-1 CRF07_BC and CRF08_BC . The study indicated that there may be a correlation between HIV-1 and HCV, but it did not focus on the same population. A previous study identified the origin and transmission patterns of HIV-1 and HCV in former blood donor (FBD) patients with coinfections. The results showed that HIV-1 infections in FBDs were introduced from Thailand and HCV infections were introduced from Japan ; however, this study did not analyze the correlation between HIV-1 and HCV. Here, our study documented that eight HCV subtypes and nine HIV-1 subtypes were found among HIV-1/HCV coinfected drug users and that there was a possible association in transmission patterns between HIV-1 and HCV.
A cross-sectional study was conducted among drug users recruited from communities and methadone maintenance treatment programs with the assistance of local Centers for Disease Control and Prevention (CDC) in the Yingjiang and Kaiyuan prefectures of Yunnan province between 2009 and 2011 (Fig. 1). The study recruited 310 HIV-1/HCV coinfected drug users (183 from Kaiyuan prefecture and 127 from Yingjiang prefecture). Information regarding demographic characteristics, such as age, gender, and ethnicity, was collected from the interview. Venous blood was collected, centrifuged, and stored in a − 80 °C freezer.
Amplification of HIV-1 and HCV gene fragments
Viral RNA was extracted from the plasma of HIV-1/HCV coinfected drug users using a High Pure Viral RNA Kit (Roche, Mannheim, Germany) and then subjected to amplification of HIV-1 env, p17, pol and HCV C/E2 and NS5B genomic fragments. Reverse transcription PCRs were performed using the PrimeScript™ II 1st Strand cDNA Synthesis Kit (TaKaRa Biotechnology, Dalian, China), and then the PCR products were subjected to nested PCR using TransTaq DNA Polymerase High Fidelity (Beijing TransGen Biotech Co., Ltd., Beijing, China). The primer pairs used in the nested PCR were modified from previous papers [11,12,13]. The nested PCR products were sequenced with an ABI PRISM 377XL DNA sequencer (Applied BioSystems, California, America).
Phylogenetic analyses of HIV-1 and HCV sequences
The HIV-1 and HCV sequences were aligned with a set of reference sequences available at Los Alamos National Laboratory HIV-1 Database and HCV Database using ClustalW. The alignment was edited using BioEdit 5.0.9. The gaps were removed manually, and the sequences were trimmed to obtain fragments of equivalent length. Phylogenetic trees were generated using the maximum-likelihood method implemented in MEGA 7.0 software, and the branch significance was analyzed by bootstrapping with 1000 replicates. To determine HIV-1 recombination, bootscan analyses were performed using SimPlot 3.5.1 software.
Phylogeographic analysis of HIV-1 and HCV sequences
To investigate the possible correlation of HIV-1 and HCV among HIV-1/HCV coinfected drug users, reference HIV-1 p17 and HCV NS5B sequences from prefectures outside Yingjiang and Kaiyuan in Yunnan were downloaded from the HIV-1 and HCV databases. Subsequently, HIV-1 subtype C (including the sequences of HIV-1 CRF07_BC and CRF08_BC, as their p17 fragments were subtype C) and HCV 3b sequences were selected according to the results of the maximum-likelihood tree analysis and the online tool “Recombinant Identification Program” available in HIV-1 database. Bayesian phylogeographic analysis was then performed using BEAST 1.6.2 as previously described [11, 14, 15]. The accession numbers of reference sequences downloaded from GenBank and used for Bayesian phylogeographic analysis are listed in Additional file 1: Table S1 and Additional file 2: Table S2.
The sequences reported in this article are available in GenBank under accession numbers MF990907-MF991024. MG334013-MG334120, MG385895-MG386024, MG432013-MG432106, MG450396-MG450541, MG461702-MG461853, MG494384-MG494479, MG549076-MG549188, MG549193-MG549305, MG763455-MG763560.
Most of the 183 HIV-1/HCV coinfected samples from Kaiyuan prefecture were from men who were of Han ethnicity, unemployed and single and who had an education level of secondary school (Table 1). Most of the 127 HIV-1/HCV coinfected samples from Yingjiang prefecture were from men who were of Dai ethnicity, farmers and single and who had an education level of primary school (Table 1).
Of the participants in Kaiyuan, 42 were self-reported as undergoing antiretroviral therapy (ART), 62 were ART-naïve, and 79 had not answered the questions regarding the status of ART; the numbers in Yingjiang were 26, 33, and 68, respectively. Information on patients receiving anti-HCV therapy was not collected in the present study.
Subtype characterization of HIV-1 and HCV strains
Among the 183 HIV-1/HCV coinfected samples from Kaiyuan prefecture, 84.2% and 70.0% of them were successfully amplified with at least one fragment of the HCV C/E2, NS5B and HIV-1 p17, pol, env genes, respectively (Additional file 3: Table S3). Maximum-likelihood trees showed that HCV 3b was predominantly based on HCV C/E2 and NS5B, followed by 3a, 6n, 1b, 6a and 6v (Fig. 2). HIV-1 subtype C was predominantly based on HIV-1 p17 and env fragments, followed by subtype B and CRF01_AE (Additional file 4: Figure S1). While for the HIV-1 pol fragments, CRF08_BC was predominant (Fig. 3). Taken together, as the p17 fragments of CRF07_BC and CRF08_BC is subtype C, the results showed that HCV 3b (45.4%) and HIV-1 CRF08_BC (66.1%) were the most predominant subtypes among HIV-1/HCV coinfected drug users from Kaiyuan prefecture.
Among 127 HIV-1/HCV coinfected samples from Yingjiang prefecture, 90.6% and 89.8% of them were successfully amplified with at least one fragment of the HCV C/E2, NS5B and HIV-1 p17, pol, env genes, respectively (Additional file 3: Table S3). Maximum-likelihood trees showed that HCV 3b was predominantly based on HCV C/E2 and NS5B, followed by 6u, 6n, 1a, 3a, 6v and 1b (Fig. 2). HIV-1 subtype C was predominantly based on HIV-1 p17 and env fragments, followed by subtype B and CRF01_AE (Additional file 1: Figure S1). While for the HIV-1 pol fragments, URFs comprising HIV-1 subtypes B and C (URF_BC) were predominant (Fig. 3). Taken together, as the p17 fragments of CRF07_BC and CRF08_BC is subtype C, the results showed that HCV 3b (43.4%) and URF_BC (53%) were the most predominant subtypes among HIV-1/HCV coinfected drug users from Yingjiang prefecture.
Spatial dynamics of HIV-1 and HCV
A maximum clade credibility tree based on HIV-1 subtype C p17 fragments showed that HIV-1 among drug users in Kaiyuan was introduced from Yingjiang in 2004 and then spread to other prefectures of Yunnan Province (Fig. 4). The results based on HCV 3b NS5B fragments showed that HCV among drug users in Kaiyuan was introduced from Yingjiang through multiple lineages in 2007 and then spread to other prefectures of Yunnan Province (Fig. 4). These results indicated that the transmission of HIV-1 and HCV among HIV-1/HCV coinfected patients had high consistency.
Due to high-risk behaviors such as drug injection and sharing of needles, drug users, especially those who inject drugs, are at high risk of HIV-1 transmission in China. The first cases of HIV-1 in China were found in drug users of Yunnan province , and then HIV-1 spread in China rapidly. By the end of 2017, the cumulative number of HIV-1-positive individuals reported in China was 718,270. Among these HIV-1-infected cases, a high proportion are transmitted through drug users. Data from the 1995 to 2011 China National AIDS Case Report and Sentinel Surveillance System showed that the rate of HIV-1 infection among drug users increased rapidly before 2000. After 2004, the rate of HIV-1 infection among drug addicts declined, with an infection rate of 6.4% in 2011 . In China, there are differences in HIV-1 infection rates in different regions. In areas close to the “Golden Triangle” and along drug trafficking routes, such as Xinjiang and Yunnan, HIV-1 infection among drug users will reach 40% [18, 19]. Similar to HIV-1, HCV can also transmit among drug users and is considered to be the most common viral infection among drug users . Systematic analysis showed that the prevalence of HCV in Chinese drug addicts reached 50% [21, 22], while the HCV infection rate in Yunnan was 90% . As HIV-1 and HCV share similar transmission routes, HIV-1/HCV coinfection is more common among drug users, while HIV-1/HCV coinfection makes the management of infection in drug users more complicated. The coinfection rate of HIV-1/HCV is not consistent in different regions. The data showed that the highest rate of coinfection was 46.3% in Yunnan, followed by 36.2% in Xinjiang .
Yunnan province is located in southwestern China and situated along drug trafficking routes that channel heroin into China. Laos, Vietnam and Thailand all border this province. Dehong prefecture is bordered by Myanmar, and Honghe prefecture borders Vietnam. In this study, the prevalence of HIV-1/HCV coinfection was 32% and 56% among drug users of Yingjiang and Kaiyuan prefectures, Yunnan province, respectively.
The main subtypes of HIV-1 circulating in drug users in 1990 were B, C and CRF01_AE [24, 25]. From 2000 to 2001, a study analyzed HIV-1 subtypes of PWIDs of Dehong, Honghe and Wenshan prefectures, Yunnan province. The results showed that the main HIV-1 subtype in Dehong was URFs (71%), and the other subtype was B, while in Honghe and Wenshan, CRF 08_BC was the most prevalent . At the same time, a new recombinant comprising CRF07_BC and CRF08_BC was also found in Honghe prefecture . A study also showed that there are more URFs among drug users . In this study, we found a high prevalence of HIV-1 CRF08_BC and URF_BC recombinants among HIV-1/HCV coinfected drug users in Kaiyuan province and Yingjiang prefecture, Yunnan province, respectively, which was similar to previous studies. HCV genotypes 3 and 6 were common among drug users in Yunnan province, and a higher rate of HCV genotype 6 was seen in recent years [23, 29, 30]. We also concluded that HCV 3b was the most predominant among HIV-1/HCV coinfected drug users in Yingjiang and Kaiyuan prefectures, Yunnan province. In addition, a high proportion of HCV genotype 6 was found in this population.
A previous study discovered that HCV prevalence may predict the HIV-1 epidemic , and another study illustrated that individuals with HIV-1 virus subtypes that clustered with those of HIV-1/HCV coinfected patients had a higher risk for acquiring HCV . These studies all indicate that HIV-1 and HCV may have some correlation. Our study further investigated the spatial dynamics of HIV-1 and HCV by performing phylogeographic analyses. The results showed that the transmission of the two viruses was the same. This finding can help us predict the transmission route of one virus based on the known route of another virus to some extent.
In summary, the present study explored the subtypes and spatial dynamics of HIV-1 and HCV among HIV-1/HCV coinfected drug users from Yingjiang and Kaiyuan prefectures, Yunnan province. Phylogenetic analysis based on partial HCV C/E2, NS5B and HIV-1 p17, pol and env genes revealed that HCV 3b and HIV-1 CRF08_BC and URF_BCs were predominant in this risk group. Furthermore, phylogeographic analysis confirmed that the transmission patterns of HIV-1 and HCV were highly consistent, providing new insights into the prevention of HIV-1 and HCV.
Availability of data and materials
The datasets generated during the current study are available in the GenBank, https://www.ncbi.nlm.nih.gov/, with accession numbers MF990907-MF991024. MG334013-MG334120, MG385895-MG386024, MG432013-MG432106, MG450396-MG450541, MG461702-MG461853, MG494384-MG494479, MG549076-MG549188, MG549193-MG549305, MG763455- MG763560. All the original data are available from the corresponding author as required.
Platt L, Easterbrook P, Gower E, McDonald B, Sabin K, McGowan C, Yanny I, Razavi H, Vickerman P. Prevalence and burden of HCV co-infection in people living with HIV: a global systematic review and meta-analysis. Lancet Infect Dis. 2016;16(7):797–808.
Xie J, Han Y, Qiu Z, Li Y, Li Y, Song X, Wang H, Thio CL, Li T. Prevalence of hepatitis B and C viruses in HIV-positive patients in China: a cross-sectional study. J Int AIDS Soc. 2016;19(1):20659.
Zhou YH, Yao ZH, Liu FL, Li H, Jiang L, Zhu JW, Zheng YT. High prevalence of HIV, HCV, HBV and co-infection and associated risk factors among injecting drug users in Yunnan province, China. PLoS ONE. 2012;7(8):e42937.
Beyrer C, Razak MH, Lisam K, Chen J, Lui W, Yu X-F. Overland heroin trafficking routes and HIV-1 spread in south and south-east Asia. AIDS. 2000;14(1):75–83.
Deng X, Liu H, Shao Y, Rayner S, Yang R. The epidemic origin and molecular properties of B′: a founder strain of the HIV-1 transmission in Asia. AIDS. 2008;22(14):1851–8.
Takebe Y, Liao H, Hase S, Uenishi R, Li Y, Li X-J, Han X, Shang H, Kamarulzaman A, Yamamoto N. Reconstructing the epidemic history of HIV-1 circulating recombinant forms CRF07_BC and CRF08_BC in East Asia: the relevance of genetic diversity and phylodynamics for vaccine strategies. Vaccine. 2010;28:B39–44.
Zhang Y, Chen L-M, He M. Hepatitis C Virus in mainland China with an emphasis on genotype and subtype distribution. Virol J. 2017;14(1):1–14.
Ye M, Chen X, Wang Y, Duo L, Zhang C, Zheng Y-T. Identification of a new HCV subtype 6xg among injection drug users in Kachin, Myanmar. Front Microbiol. 2019;10:814.
Chen M, Ma Y, Chen H, Dai J, Luo H, Jia M, Song Z. Complete genome sequencing and evolutionary analysis of HCV subtype 6xg from IDUs in Yunnan, China. PLoS ONE. 2019;14(5):e0217010.
Zhang C, Wu N, Liu J, Ge Q, Huang Y, Ren Q, Feng Q, He G. HCV subtype characterization among injection drug users: implication for a crucial role of Zhenjiang in HCV transmission in China. PLoS ONE. 2011;6(2):e16817.
Liu J, Zhang C. Phylogeographic analyses reveal a crucial role of Xinjiang in HIV-1 CRF07_BC and HCV 3a transmissions in Asia. PLoS ONE. 2011;6(8):e23347.
Du L, Wu J, Qian P, Xin R, Ni Y, Han R, Meng Z, Xia J. Phylogeographical analysis reveals distinct sources of HIV-1 and HCV transmitted to former blood donors in China. AIDS Res Hum Retroviruses. 2017;33(3):284–9.
Pang W, Zhang C, Duo L, Zhou Y-H, Yao Z-H, Liu F-L, Li H, Tu Y-Q, Zheng Y-T. Extensive and complex HIV-1 recombination between B’, C and CRF01_AE among IDUs in south-east Asia. AIDS. 2012;26(9):1121–9.
Liu J, Jia Y, Xu Q, Zheng Y-T, 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(8):1944–8.
Chen X, Zhou Y-H, Ye M, Wang Y, Duo L, Pang W, Zhang C, Zheng Y-T. Burmese injecting drug users in Yunnan play a pivotal role in the cross-border transmission of HIV-1 in the China-Myanmar border region. Virulence. 2018;9(1):1195–204.
Zheng X, Tian C, Choi K-H, Zhang J, Cheng H, Yang X, Li D, Lin J, Qu S, Sun X. Injecting drug use and HIV infection in southwest China. AIDS. 1994;8(8):1141–7.
Wang L, Guo W, Li D, Ding Z, McGoogan JM, Wang N, Wu Z, Wang L, Program CNHAS. HIV epidemic among drug users in China: 1995–2011. Addiction. 2015;110:20–8.
Bao Y, Liu Z. Systematic review of HIV and HCV infection among drug users in China. Int J STD AIDS. 2009;20(6):399–405.
Yao Y, Wang N, Chu J, Ding G, Jin X, Sun Y, Wang G, Xu J, Smith K. Sexual behavior and risks for HIV infection and transmission among male injecting drug users in Yunnan, China. Int J Infect Dis. 2009;13(2):154–61.
Agdamag DM, Kageyama S, Alesna ET, Solante RM, Leaño PS, Heredia AML, Abellanosa-Tac-An IP, Vibal ET, Jereza LD, Ichimura H. Rapid spread of hepatitis C virus among injecting-drug users in the Philippines: implications for HIV epidemics. J Med Virol. 2005;77(2):221–6.
Bennett H, Waser N, Johnston K, Kao J-H, Lim Y-S, Duan Z-P, Lee Y-J, Wei L, Chen C-J, Sievert W. A review of the burden of hepatitis C virus infection in China, Japan, South Korea and Taiwan. Hepatol Int. 2015;9(3):378–90.
Xia X, Luo J, Bai J, Yu R. Epidemiology of hepatitis C virus infection among injection drug users in China: systematic review and meta-analysis. Public Health. 2008;122(10):990–1003.
Xia X, Lu L, Tee KK, Zhao W, Wu J, Yu J, Li X, Lin Y, Mukhtar MM, Hagedorn CH. The unique HCV genotype distribution and the discovery of a novel subtype 6u among IDUs co-infected with HIV-1 in Yunnan, China. J Med Virol. 2008;80(7):1142–52.
Cheng H, Zhang J, Capizzi J, Young N, Mastro T. HIV-1 subtype E in Yunnan, China. Lancet. 1994;344(8927):953–4.
Luo C-C, Tian C, Hu D, Kai M, Dondero T, Zheng X. HIV-1 subtype C in China. Lancet. 1995;345(8956):1051–2.
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(10):1401–7.
Yang R, Kusagawa S, Zhang C, Xia X, Ben K, Takebe Y. 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, China. J Virol. 2003;77(1):685–95.
Li L, Chen L, Yang S, Li T, Li J, Liu Y, Jia L, Yang B, Bao Z, Li H, et al. Recombination form and epidemiology of HIV-1 unique recombinant strains identified in Yunnan, China. PLoS ONE. 2012;7(10):e46777.
Zhang Z, Yao Y, Wu W, Feng R, Wu Z, Cun W, Dong S. Hepatitis C virus genotype diversity among intravenous drug users in Yunnan Province, Southwestern China. PLoS ONE. 2013;8(12):e82598.
Wan Z, Chen Q, Chen X, Duo L, Li P, Zheng Y-T, Zhang C. HCV diversity among Chinese and Burmese IDUs in Dehong, Yunnan, China. PLoS ONE. 2016;11(9):e0163062.
Akbarzadeh V, Mumtaz GR, Awad SF, Weiss HA, Abu-Raddad LJ. HCV prevalence can predict HIV epidemic potential among people who inject drugs: mathematical modeling analysis. BMC Public Health. 2016;16(1):1–17.
Kouyos RD, Rauch A, Böni J, Yerly S, Shah C, Aubert V, Klimkait T, Kovari H, Calmy A, Cavassini M. Clustering of HCV coinfections on HIV phylogeny indicates domestic and sexual transmission of HCV. Int J Epidemiol. 2014;43(3):887–96.
We thank the staff of the local CDC for the assistance in collecting samples. We also thank all the participants who participated in this research.
This work was supported by grants from the National Natural Science Foundation of China (U1302224, 81601802, 81271892, and U1802284), Project of Innovative Research Team of Yunnan Province (202005AE160005). The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Ethical approval for the study and the informed consent process were approved by the Ethics Committee of Kunming Institute of Zoology, Chinese Academy of Sciences (approval ID: SWYX-2008010; approval date: 9 January 2008). The research was conducted in accordance with basic principles of the Helsinki declaration and the relevant international rules. Written informed consent was obtained from all participants.
Consent for publication
The authors declare that there are no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The reference sequences of HCV NS5B fragments that downloaded from GenBank and used for Bayesian phylogeographic analysis.
The reference sequences of HIV-1 p17 fragments that downloaded from GenBank and used for Bayesian phylogeographic analysis.
The amplification results of HIV-1 and HCV fragments among HIV-1/HCV coinfected drug users in Yingjiang and Kaiyuan prefectures, Yunnan, China.
The maximum-likelihood trees based on env and p17 fragments of HIV-1 among HIV-1/HCV coinfected drug users in Yunnan Province, China. The red triangles and the blue circles indicate the sequences that amplified from drug users in Yingjiang and Kaiyuan Prefectures, respectively. The different colored sectors indicate the sequences with different subtypes.
About this article
Cite this article
Wang, Y., Chen, X., Ye, M. et al. Consistency of spatial dynamics of HIV-1 and HCV among HIV-1/HCV coinfected drug users in China. BMC Infect Dis 21, 1001 (2021). https://doi.org/10.1186/s12879-021-06711-6