Skip to main content

Genetic transmission networks of HIV-1 CRF07_BC strain among HIV-1 infections with virologic failure of ART in a minority area of China: a population-based study



The drug resistance and the virologic failure of antiretroviral therapy (ART) are quite severe in Liangshan. A better understanding of the virologic failure of ART and the HIV-1 transmission network dynamics is essential for the surveillance and prevention of HIV. Here, we analyzed the HIV-1 CRF07_BC strain genetic transmission networks and their associated factors among people living with HIV-1 (PLWH) who had virologic failure of ART by using close genetic links.


The drug-resistant mutations were determined using the Stanford University HIV Drug Resistance Database. HIV-1 pol genes sequences were used for phylogenetic and genotypic drug resistance analysis. The genetic transmission networks were performed by comparing sequences, constructing the phylogenetic tree, calculating the pairwise distance, and visualizing the network.


A total of 1050 PLWH with CRF07_BC pol sequences were finally identified and included in the genetic transmission network analysis from 2016 to 2017. Of the 1050 CRF07_BC pol sequences, 346 (32.95%) fell into clusters at a genetic distance of 0.006, resulting in 137 clusters ranging in size from 2 to 40 individuals. Subjects who were widowed or divorced were less likely to form a genetic transmission network (adjusted OR: 0.50), while subjects who had shared a needle ≥ five times were more likely to form a network (adjusted OR: 1.88).


The genetic transmission networks revealed the complex transmission pattern, highlighting the urgent need for transmission monitoring of virologic failure of ART and selection of more effective therapeutic regimens to promote viral suppression.

Peer Review reports


The Human Immunodeficiency Virus /Acquired Immune Deficiency Syndrome (HIV/AIDS) epidemic has become one of the most critical issues that seriously affect public health. Since the start of the epidemic, around 75 million people have been infected with HIV, and in 2018, 37.9 million people were living with HIV (PLWH) worldwide [1]. By the end of October 2019, approximately 958 thousand people lived with HIV in China [2]. Liangshan Prefecture, as the autonomous prefecture with the largest population of Yi people in southwest China, contains a huge of a proportion of PLWH in Sichuan Province, where the prevalence of PLWH was over 1% in Butuo County, Zhaojue County, Meigu county, Yuexi County, and Jinyang County, ranking among the top 5 counties in the country. According to previous sampling survey reports, the effective rate of antiretroviral therapy (ART) among PLWHin Liangshan Prefecture was only 55.9%, and the drug resistance rate of newly reported untreated patients was as high as over 5% [3]. The high virologic failure rate and drug resistance rate of ART in Liangshan Prefecture were mainly due to the poor adherence to ART, the high mobility, and the high lost-to-follow-up rate of PLWH [4, 5]. Besides, Linda et al. [6, 7] found a more unsatisfactory virologic response in patients with transmitted drug resistance than those without. Therefore, we hypothesized that the main reason for high drug resistance rate in Liangshan Prefecture might be because of transmitted drug resistance before treatment. Thus, effective intervention should be formatted to reduce drug resistance in this area.

Three circulating recombinant forms (CRFs) of HIV-1 predominate in China, namely CRF01_AE, CRF07_BC, and CRF08_BC [8]. Subtype CRF07_BC was first discovered in the population with intravenous drug abuse in Yunnan Province and then rapidly spread to Xinjiang, Guangxi, and Sichuan Provinces [9]. According to a nationwide molecular epidemiological survey in 2006, HIV-1 CRF07_BC had become one of the most widespread subtypes circulating in China [10]. In Liangshan Prefecture, CRF07_BC was the main subtype of drug users in the past two decades and was slowly introduced to heterosexual people [11]. Now, CRF07_BC has become the most common form in this area [12, 13]. This study intends to construct an HIV-1 CRF07_BC strain transmission network of PLWH with virologic failure to analyze the clustering of drug-resistant strains and observe the possibility of drug-resistant transmission.

The HIV-1 molecular transmission network refers to a group of sequences that are not randomly gathered and have a certain epidemiological correlation. It constructs a transmission network through the genetic information of people infected with HIV through similar viruses with similarities and connections, and restores the macro social network of infected people as much as possible, which aims to focus on analyzing the characteristics of infected people, and preventing and controlling the active and critical groups in the network [14, 15]. HIV is transmitted through networks formed by closely connected individuals who engage in injecting or sexual behaviours [16, 17]. Understanding the structure and features of the transmission networks among HIV individuals with virologic failure can be useful for identifying potential transmission partners and recognizing the links between different populations [15], which is essential for developing HIV intervention strategies.

To our knowledge, few studies were performed to construct the HIV genetic transmission network with pol sequences identified from individuals with virologic failure of ART in China. Through the analysis of molecular transmission networks, we could obtain drug-resistant transmission clusters to provide policymaking recommendations for precious prevention for the virologic failure of ART. In this study, we investigated the HIV-1 CRF07_BC strain transmission patterns among PLWH who had virologic failure by using close genetic links and explored the factors associated with genetic transmission networks in Liangshan Prefecture.


Ethics statement

All PLWH voluntarily participated in our study and signed informed consent forms before enrollment. The study protocol was approved by the Ethics Committee of the Sichuan Center for Disease Control and Prevention, and the study was carried out following the Helsinki Declaration of 1964.

Study participants

The inclusion criteria were: 1) being permanent residents or having stayed in the study sites for more than five years, 2) being confirmed with HIV-1 infection, 3) having received ART for more than six months, and 4) being confirmed to have a virologic failure of ART by the HIV RNA level.

Participants were recruited in Butuo county that had the most severe HIV epidemic of the Liangshan area. According to the basic information system for AIDS prevention and treatment, all PLWH of Butuo county in this system were included in our study.

A total of 5157 PLWH in Liangshan Prefecture were enrolled in this study. After excluding those with an HIV RNA level ≤ 1000 copies/ml, a total of 2156 PLWH showed virologic failure of ART. 1576 sequences were successfully obtained, and 1508 sequences of these were subtype CRF07_BC (95.69%). After excluding 69 other HIV-1 subtype sequences and 457 repeat sequences, a total of 1050 PLWH with CRF07_BC pol sequences were finally identified and included in the genetic transmission network analysis (Fig. 1).

Fig. 1

Sampling flowchart

Data collection

The collected data included participants’ socio-demographic information, HIV related behaviours information, and ART related information. Socio-demographic information included age, occupation, marital status, attained level of education, treatment time, and self-reported route of HIV transmission (e.g., injection drug use, casual sexual behaviour, and risk behaviours related to mother-to-child transmission, etc.). Information on HIV related behaviours included needle sharing, the number of casual sexual partners, and the history of sexually transmitted diseases apart from HIV. ART related information, including the time of ART, the initial regimen of ART, change of ART regimen, and drug resistance, were extracted from their medical records.

Laboratory tests

Each participant provided 5 mL of venous blood to detect the CD4+T cells and the viral load of HIV-1. The CD4+T cells were quantified within 12 h using flow cytometry at local county’s Center for Disease Control and Prevention (CDC). Plasma was isolated and sent on the cold chain to the Sichuan CDC to measure the HIV viral load. Virologic failure was defined as an HIV RNA level > 1000 copies/ml. For samples with an HIV RNA level > 1000 copies/ml, HIV drug resistance genotyping was conducted using an in-house PCR [18, 19].

Nucleic acid extraction, amplification, and sequencing

RNA was extracted by an automatic extraction machine (MagNA Pure LC 2.0 system, Roche, Branchburg, NJ, USA) and the MagNA Pure LC Total Nucleic Acid Isolation Kit-High Performance kit (Roche, Branchburg, NJ, USA). The HIV-1 Pol region was amplified by two round nested polymerase chain reactions (Nested-PCR) using RNA extracted as a template. The amplification length was 1200 bp, including protease and reverse transcriptase regions.

The first-round PCR was performed by the AccessQuick RT-PCR System (A1702) kit (Promega Biotechnology Co., Ltd., Beijing, China), and the second-round PCR was performed by the Perfect shot TM Ex Taq Mix kit (Tiangen Biochemical Technology Co., Ltd., Beijing, China). The PCR products were treated with electrophoresis with 1% agarose gel. The amplified products were sent to the Beijing Genomics Research Center Ltd. for purification and gene sequencing. The detailed amplification and sequencing were described previously [19].

For samples with an HIV RNA level > 1000 copies/ml, the evaluation of drug resistance mutations was performed by using the algorithm of the Stanford HIV drug resistance database ( All experimental protocols followed the manufacturer’s instructions.

HIV-1 genetic transmission network

The analysis of this step was described in our previous study [20]. By comparing sequences, constructing the phylogenetic tree, calculating the pairwise distance, and visualizing the network, the flow chart of the genetic transmission network was realized. To avoid the potential confounding effect of the convergent evolution of drug resistance, we constructed the genetic transmission network after excluding 32 codon positions in protease and reverse transcriptase associated with drug resistance.

After excluding repeat sequences, we performed the phylogenetic tree and molecular transmission networks based on pol sequences identified from individuals with virologic failure of ART. Calculations of the pairwise genetic distances of all sequences within the available clusters were performed by Run HyPhy 2.2.4 ( Within the genetic distances of 0.001–0.015, the gene distance with the largest number of clusters was selected as the linkage within a cluster. Finally, the network visualization analysis was based on the Cytoscape 3.5 network software ( [20].

Statistical analysis

Frequencies and percentages were used to represent the categorical variables, and the chi-square (χ2) test was used to compare the differences among groups. Taking the presence of individuals within the transmission networks as the dependent variable, unadjusted OR was estimated for the socio-demographic variables, HIV-related behaviours, ART-related information, the first CD4+ count after diagnosis, the stage of the diseases, and drug resistance. A binary regression model was fitted, and variables with significance in the univariate analysis were considered as candidate variables. SPSS version 21.0 for Windows (SPSS, Inc., Chicago, IL, USA) was used for data analysis, with p values< 0.05 taken as statistically significant.


Distribution of HIV-1 genotypes among PLWH with drug resistance

Most of the subtypes of HIV-1 were CRF07_BC (1508, 95.62%), followed by the subtypes CRF08_BC (42, 2.66%), and C (27, 1.71%). Among PLWH with virologic failure of ART, the drug-resistant rate was as high as 32.10% (507/1576). We found that 477 (30.25%) subjects showed drug resistance to non-nucleoside reverse transcriptase inhibitors (NNRTIs), 193 (12.24%) to nucleoside reverse transcriptase inhibitors (NRTIs), and 29 (1.84%) to protease inhibitors (PIs) in ART.

Characteristics of PLWH within the transmission network

A total of 1050 PLWH with CRF07_BC pol sequences were identified and included in the genetic transmission network analysis. A total of 346 PLWH were within the genetic transmission network with a genetic distance of 0.006. In the univariate analysis, using the presence of PLWH within the transmission networks as the dependent variable, PLWH with needle sharing behaviour were more likely to be included in the transmission networks than PLWH without drug use or sharing needles, and PLWH who were widowed or divorced were less likely to be included in the transmission networks than PLWH with marriage (Table 1).

Table 1 The distribution of PLWH with virologic failure in ART in the genetic transmission networks

In the summary of the binary logistic regression model, we found that PLWH who were widowed or divorced were less likely to form a genetic transmission network (adjusted OR: 0.50, 95%CI: 0.25–0.99; reference group: married), while PLWH who had shared a needle ≥5 times were more likely to be in a transmission network (adjusted OR: 1.88, 95%CI: 1.02–3.46; reference group: had never used illicit drugs or shared a needle) (Table 2). The ratio of men to women being widowed or divorced was 4:7, while the ratio was 6:1 among singles (Data not shown).

Table 2 Factors associated with the presence of PLWH within the genetic transmission networks

HIV-1 genetic transmission network

A total of 137 clusters were observed with the number of sequences per cluster ranging from 2 to 40. Of these clusters, 107 clusters had only two linked PLWH (one link), and 30 clusters had ≥2 potential transmission partners (Fig. 2).

Fig. 2

Genetic transmission networks by drug resistance. 1–40 are transmission clusters containing drug-resistant sequences

A total of 118 drug resistance sequences were observed in the transmission network. Among these sequences, 101 sequences (85.59%) were linked to at least one other drug resistance sequence, and 40 clusters had > 2 drug resistance sequences, indicating that PLWH with drug resistance were more likely to link with those with drug resistance (Fig. 2). Among these 40 clusters, 39 clusters (97.5%) had common drug resistance sites, 22 clusters (55%) shared the common site of K103N, 16 clusters (40%) shared the common site of M184V/I, and 9 clusters (22.5%) shared the common site of Y181C, which suggested that the linked sequences with drug resistance tended to have common drug resistance sites (Fig. 2, Additional file 1). On the HIV transmission routes, we found that drug injection and heterosexual contact were often cross-linked (Fig. 3). Moreover, 51.54% (50/97) of the clusters with two PLWH, 57.89% with 3–4 PLWH, and all clusters with more than 5 PLWH were from different villages or towns (Fig. 3).

Fig. 3

Genetic transmission networks by places of residence


Liangshan Prefecture is the largest traditional settlement of the Yi minority people in China and is one of the areas in China severely affected by HIV [21]. In Liangshan, the rapid scale-up of ART led to a significant decline in morbidity and mortality among PLWH [22]. However, the virologic failure of ART and drug resistance are still serious in this area. Therefore, monitoring the virologic failure of ART and understanding the HIV-1 transmission network dynamics among patients with virologic failure of ART are urgently needed.

We found that a large number of sequences of drug resistance to ART were identified in the genetic transmission networks through the genetic transmission network analysis. In terms of drug resistance-associated mutations in transmission clusters, 35.98% of the individuals carrying drug resistance were involved in the network. 101 PLWH (85.59%) with drug resistance were linked with each other, which suggested the possibility of drug resistance transmission. Moreover, most of the linked sequences had common drug resistance sites. Among 40 drug resistance clusters, K103N was observed in 22 clusters as the common drug resistance site and was the most frequent drug resistance-associated mutation. The possible reason might be that K103N mutation can persist without drug pressure and promote infection to another host [23, 24], which has been attributed to the widespread use of NNRTI in first-line ART [25]. Similar to K103N, M184V/I and Y181C can also be stably passaged without drug pressure [26] which may be the reason for the higher frequency of these drug resistance mutations in the transmission network. This finding was also observed in other undeveloped [25] and developed settings [27, 28]. The widespread use of ART may have led to a considerable increase in HIV drug resistance and drug resistance transmission in newly infected individuals. Drug resistance transmission may have potentially compromised the efficacy of the combination of ART regimens and may lead to failure of ART [29, 30]. Thus, there is an urgent need for routine resistance testing before starting treatment since the transmission of drug resistance was associated with virologic failure in patients [6, 7]. Also, close linkages between drug-resistance strains may be due to the high mutation of similar strains in the cluster, indicating that the ART regimen of PLWH within the cluster should be replaced if the drug-resistant mutations occurred in the same cluster.

Furthermore, our study indicated that most PLWH in the clusters were from different villages or towns. Unlike other minority groups in China, the Yi minority people in Liangshan Prefecture are more likely to have casual sexual behaviours [21, 31, 32]. Casual sex often involves condomless sexual behaviours and multiple sex partners, facilitating the transmission of HIV [31, 32]. With the progress of society, the network has become more developed, and the population mobility has increased, which may increase the risk of the spread of HIV-1 and drug resistance transmission. Therefore, we should enhance health education on HIV-related behaviours for patients, especially the floating population, to reduce HIV-1 transmission.

We also found that widowed or divorced PLWH were less likely to link with other partners with genetically linked viruses. Most of the widowed or divorced PLWH were women (the men to women ratio were 4:7). These women may have been infected with HIV by their husbands through heterosexual transmission [33]. Because of the local social norms, widowed or divorced Yi women rarely have casual sex and extramarital sexual behaviours, making it challenging to spread HIV to others [21]. Therefore, effective and tailor-made programs are needed to reduce sexual transmission among unmarried youth and HIV-1 discordant couples.


Four limitations should be highlighted. First, this study was a cross-sectional survey, time-based sequences and cause-effect relationships among these variables cannot be established. Second, many factors, such as the complication of ART, the number of sexual partners, insurance status, and socioeconomic status should be considered as potential factors for influencing drug resistance and transmission networks. However, we could not obtain detailed epidemiological information about the sequences that were in the clusters. Further, follow-up medical records are warranted to include more epidemiological information. Third, the phylogenetic analysis is limited by the impossibility of inferring the direction of transmission, and genetic linkage could not reflect direct transmission due to lack of behavioural information of participants [34]. Fourth, since we did not include the sequences of all PLWH for analysis, we could not draw a whole transmission network, which limited the representation of the PLWH in this area.


To our knowledge, we were the first to investigate the risk factors of transmission networks among HIV patients with virologic failure of ART in Liangshan. Our findings have enhanced the understanding of the HIV-1 genetic characteristics and transmission networks. Understanding the determinants of the genetic transmission network may contribute to designing preventive interventions for preventing HIV-1 transmission. However, it also highlighted the urgent need to monitor the transmission of virologic failure of ART and conduct effective therapeutic regimens to promote viral suppression.

Availability of data and materials

The quantitative data used to support the findings of this study were supplied by the Center for AIDS/STD Control and Prevention, Sichuan Center for Disease Control and Prevention under license, and so cannot be made freely available. Requests for access to these data should be made to Shujuan Yang,



Human Immunodeficiency Virus


Acquired Immune Deficiency Syndrome


People were living with HIV


Antiretroviral therapy


Circulating recombinant forms


Drug resistance


  1. 1.

    Union: AIDS [] Accessed 22 March 2020.

  2. 2.

    Bureau of disease prevention and control: In 2019, China's AIDS prevention and control work made new progress [] Accessed 1 August 2020.

  3. 3.

    Xu JL, Jike KCN, Ma Y, Yu G, Wang J, Wang K, Xu P, Liao Q, Gong YH. Zhong SY et al: [effects and associated factors of HIV/AIDS anti-retroviral therapy in Liangshan Yi autonomous prefecture, Sichuan Province]. Zhonghua Yu Fang Yi Xue Za Zhi. 2018;52(6):668–72.

  4. 4.

    Bian SC, Chen XR, Shen MG, Xiao L. Accepted HAART patients viral result of Liangshan Yi autonomous prefecture from 2011 to 2015. J Prev Med Inf. 2018;34(7):944–7.

    Google Scholar 

  5. 5.

    Miao FY, Wang QX, Jin CR. F.J. Z, R.D. M, Jimu ZS, Zhang TD, Yu HL, Ye SD, Gong YH et al: the analysis of reasons and impact factors about withdrawl from antiretroviral therapy among HIV/AIDS patients in Liangshan autonomous prefecture. Chin J AIDS STD. 2014;20(4):247–50.

    Google Scholar 

  6. 6.

    Wittkop L, Günthard HF, de Wolf F, Dunn D, Cozzi-Lepri A, de Luca A, Kücherer C, Obel N, von Wyl V, Masquelier B, et al. Effect of transmitted drug resistance on virological and immunological response to initial combination antiretroviral therapy for HIV (EuroCoord-CHAIN joint project): a European multicohort study. Lancet Infect Dis. 2011;11(5):363–71.

    Article  Google Scholar 

  7. 7.

    Bansi L, Geretti AM, Dunn D, Hill T, Green H, Fearnhill E, Gazzard B, Nelson M, Porter K, Phillips A, et al. Impact of transmitted drug-resistance on treatment selection and outcome of first-line highly active antiretroviral therapy (HAART). J Acquir Immune Defic Syndr. 2010;53(5):633–9.

    CAS  Article  Google Scholar 

  8. 8.

    Feng Y, Takebe Y, Wei H, He X, Hsi JH, Li Z, Xing H, Ruan Y, Yang Y, Li F, 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.

    CAS  Article  Google Scholar 

  9. 9.

    Li W, Li X, He Y, Ge Y, Ong JJ, Li X, Dong X, Chu J, Musa TH, Cao S, et al. The evolutionary and transmission characteristic of HIV-1 CRF07_BC in Nanjing. Journal of Medical Virology: Jiangsu; 2020.

    Google Scholar 

  10. 10.

    Jiang S, He X, Xing H, Ruan Y, Hong K, Cheng C, Hu Y, Xin R, Wei J, Feng Y 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(10).

  11. 11.

    Li L, Wei D, Hsu WL, Li T, Gui T, Wood C, Liu Y, Li H, Bao Z, Liu S, et al. CRF07_BC strain dominates the HIV-1 epidemic in injection drug users in Liangshan prefecture of Sichuan, China. AIDS Res Hum Retroviruses. 2015;31(5):479–87.

    CAS  Article  Google Scholar 

  12. 12.

    Vega Y, Delgado E, Fernández-García A, Cuevas MT, Thomson MM, Montero V, Sánchez M, Sánchez AM, Pérez-Álvarez L, Country. SGftSoNH-DiGaB: Epidemiological Surveillance of HIV-1 Transmitted Drug Resistance in Spain in 2004–2012: Relevance of Transmission Clusters in the Propagation of Resistance Mutations. PLoS One 2015, 10(5):e0125699.

  13. 13.

    Dong AB, Xiao L, Liang S, Liu L, Kang RH, Zhao S, Ruan YH, Qu SL. Xing H: [HIV-1 drug resistance and influencing factors among people living with HIV/AIDS before antiretroviral therapy in Liangshan Yi autonomous prefecture]. Zhonghua Liu Xing Bing Xue Za Zhi. 2019;40(6):648–53.

    CAS  PubMed  Google Scholar 

  14. 14.

    Gan MZ, Feng Y. Xing H: [Progress in research on the transmission characteristics of HIV-infected persons based on molecular network method]. Zhonghua Liu Xing Bing Xue Za Zhi. 2019;40(11):1487–91.

    CAS  PubMed  Google Scholar 

  15. 15.

    Oster AM, Wertheim JO, Hernandez AL, Ocfemia MC, Saduvala N, Hall HI. Using molecular HIV surveillance data to understand transmission between subpopulations in the United States. J Acquir Immune Defic Syndr. 2015;70(4):444–51.

    Article  Google Scholar 

  16. 16.

    Vasylyeva TI, Friedman SR, Paraskevis D, Magiorkinis G. Integrating molecular epidemiology and social network analysis to study infectious diseases: towards a socio-molecular era for public health. Infect Genet Evol. 2016;46:248–55.

    Article  Google Scholar 

  17. 17.

    Delva W, Leventhal GE, Helleringer S. Connecting the dots: network data and models in HIV epidemiology. AIDS. 2016;30(13):2009–20.

    Article  Google Scholar 

  18. 18.

    Wang J, He C, Hsi JH, Xu X, Liu Y, He J, Ling H, Ding P, Tong Y, Zou X, et al. Virological outcomes and drug resistance in Chinese patients after 12 months of 3TC-based first-line antiretroviral treatment, 2011-2012. PLoS One. 2014;9(2):e88305.

    Article  Google Scholar 

  19. 19.

    Xing H, Ruan Y, Li J, Shang H, Zhong P, Wang X, Liao L. L iH, Zhang M, Xue Y et al: HIV drug resistance and its impact on antiretroviral therapy in Chinese HIV-infected patients. PLoS One. 2013;8(2):e54917.

  20. 20.

    Yuan D, Du Z, Zhou J, Ye L, Su L, Yang H, Yuan F, Li Y, Liu H, Zhai W et al: HIV-1 subtype diversity, drug resistance, and genetic transmission networks in men who have sex with men with virologic failure in antiretroviral therapy in Sichuan, China, 2011 to 2017. Medicine 2019, 98(43).

  21. 21.

    Yang S, Luo M, Zhang S, Yao Y, Wang Q, Liao Q. Overview on culture and customs related to AIDS epidemic and prevention among Yi people in Liangshan. Chinese Journal of AIDS & STD. 2017;23(3):101–5.

    Google Scholar 

  22. 22.

    Günthard HF, Aberg JA, Eron JJ, Hoy JF, Telenti A, Benson CA, Burger DM, Cahn P, Gallant JE, Glesby MJ, et al. Antiretroviral treatment of adult HIV infection: 2014 recommendations of the international antiviral society-USA panel. JAMA. 2014;312(4):410–25.

    Article  Google Scholar 

  23. 23.

    Ruelle J, Ingels M-G, Jnaoui K, Ausselet N, Vincent A, Sasse A, Verhofstede C, Goubau P. Transmission network of an HIV type 1 strain with K103N in young Belgian patients from different risk groups. AIDS Res Hum Retrovir. 2013;29(10):1306–9.

    CAS  Article  Google Scholar 

  24. 24.

    Barth RE, Wensing AM, Tempelman HA, Moraba R, Schuurman R, Hoepelman AI. Rapid accumulation of nonnucleoside reverse transcriptase inhibitor-associated resistance: evidence of transmitted resistance in rural South Africa. AIDS. 2008;22(16):2210–2.

    CAS  Article  Google Scholar 

  25. 25.

    Hassan AS, Esbjornsson J, Wahome E, Thiong'o A, Makau GN, Price MA, Sanders EJ. HIV-1 subtype diversity, transmission networks and transmitted drug resistance amongst acute and early infected MSM populations from coastal Kenya. PLoS One. 2018;13(12):e0206177.

    Article  Google Scholar 

  26. 26.

    Jiao LY, Bao ZY, Li HP, Guo DX, Li JYJCJoM, Immunology: Evolvement of HIV-1 drug resistant mutations in vitro without drug pressure. 2010, 30(5):431–437.

  27. 27.

    Mbisa JL, Fearnhill E, Dunn DT, Pillay D, Asboe D, Cane PA. Evidence of self-sustaining drug resistant HIV-1 lineages among untreated patients in the United Kingdom. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America. 2015;61(5):829–36.

    Article  Google Scholar 

  28. 28.

    Brooks JI, Niznick H, Ofner M, Merks H, Angel JB. Local phylogenetic analysis identifies distinct trends in transmitted HIV drug resistance: implications for public health interventions. BMC Infect Dis. 2013;13:509.

    Article  Google Scholar 

  29. 29.

    HIVCC T. The effect of combined antiretroviral therapy on the overall mortality of HIV-infected individuals. Aids. 2010;24(1):123–37.

    Article  Google Scholar 

  30. 30.

    Yang J, Xing H, Niu J, Liao L, Ruan Y, He X, Feng Y, Li Z, Shao Y. The emergence of HIV-1 primary drug resistance genotypes among treatment-naïve men who have sex with men in high-prevalence areas in China. Arch Virol. 2013;158(4):839–44.

    CAS  Article  Google Scholar 

  31. 31.

    Yang Y, Luan RS, Liu P, Wu CL, Zhou Y, Chen W. Casual sex and concurrent sexual partnerships among young people from an Yi community with a high prevalence of HIV in China. Asian journal of andrology. 2012;14(5):758–65.

    Article  Google Scholar 

  32. 32.

    Qin S, Li AS, Ma MJ, Luo Q, Wang XL, Nan L, Zeng C, Xiao L, Li QQ, Zuo JB, et al. Unprotected sex with casual partners: a neglected source of HIV transmission among members of the Yi minority in southwestern China. Biomed Environ Sci. 2014;27(10):824–31.

    PubMed  Google Scholar 

  33. 33.

    Yu G, Yao YN, Yang SJ, Wang QX, Liao Q, Nan L, Xiao L, Gong YH, Liu ZX, Wang K, et al. HIV recent infection among women in Liangshan Yi autonomous prefecture during 2011-2014. Chin J AIDS STD. 2018;24(1):70–83.

    Google Scholar 

  34. 34.

    Levintow SN, Okeke NL, Hue S, Mkumba L, Virkud A, Napravnik S, Sebastian J, Miller WC, Eron JJ, Dennis AM: Prevalence and Transmission Dynamics of HIV-1 Transmitted Drug Resistance in a Southeastern Cohort. Open Forum Infect Dis 2018, 5(8):ofy178.

Download references


We thank the Center for AIDS/STD Control and Prevention, the Sichuan Center for Disease Control and Prevention, and the Butuo Center for Disease Control and Prevention to support data. We also thank the International Institute of Spatial Lifecourse Epidemiology (ISLE) for the research support.


This research was funded by the National Natural Science Foundation of China (81703279); the Sichuan Science and Technology Program (2019YJ0148); National Major Science and Technology Projects of China: Study on the assessment technology and prediction model of AIDS epidemic in China (2017ZX10201101); and Science and Technology Project of Sichuan Provincial Health Committee (20PJ121). The funding body had no role in the data collection, analysis, and interpretation of the verbal data, and writing of the manuscript.

Author information




SY and DY conducted the analysis. DY, ML, and SY interpreted the data and wrote the paper. YL, PJ, LS, LY, YZ, PL, CZ, SL, HY, and HL help collect and analyze the data. All authors have read and approved the manuscript.

Corresponding authors

Correspondence to Honglu Liu or Shujuan Yang.

Ethics declarations

Ethics approval and consent to participate

All subjects voluntarily participated in our study and signed informed consent forms before enrollment. The study protocol was approved by the Ethics Committee of the Sichuan Center for Disease Control and Prevention (201303), and the study was carried out following the Helsinki Declaration of 1964.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1.

Analysis of drug resistance sites in clusters containing two or more drug resistance sequences

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yuan, D., Liu, M., Li, Y. et al. Genetic transmission networks of HIV-1 CRF07_BC strain among HIV-1 infections with virologic failure of ART in a minority area of China: a population-based study. BMC Infect Dis 20, 631 (2020).

Download citation


  • HIV-1
  • Genetic transmission network
  • CRF07_BC
  • Virologic failure
  • Antiretroviral therapy