Skip to main content

HIV-1 molecular transmission network and drug resistance in Chongqing, China, among men who have sex with men (2018–2021)



Over the past few years, HIV transmission among men who have sex with men (MSM) in China has increased significantly. Chongqing, located in the southwest of China, has the highest prevalence of HIV among MSM in the country.


Blood samples were taken from 894 MSM in Chongqing who had recently been diagnosed with HIV-1 infection and had not yet started getting treatment. In order to determine the distribution of HIV-1 subtypes, transmitted drug resistance, and assessments of molecularly transmitted clusters, we sequenced the Pol genes and employed them in phylogenetic analysis. The genetic distance between molecular clusters was 1.5%. To find potential contributing factors, logistic regression analyses were performed.


Of the 894 HIV-1 pol sequences acquired from study participants, we discovered that CRF07_BC (73.6%) and CRF01_AE (19.6%) were the two most prevalent HIV-1 genotypes in Chongqing among MSM, accounting for 93.2% of all infections. In addition, CRF08_BC (1.1%), B subtype (1.0%), CRF55_01B (3.4%), and URF/Other subtypes (1.3%) were less frequently observed. Among MSM in Chongqing, transmitted drug resistance (TDR) was reported to be present at a rate of 5.6%. 48 clusters with 600 (67.1%, 600/894) sequences were found by analysis of the molecular transmission network. The distributions of people by age, sexual orientation, syphilis, and genotype were significantly differentially related to being in clusters, according to the multivariable logistic regression model.


Despite the low overall prevalence of TDR, the significance of genotypic drug resistance monitoring needs to be emphasized. CRF07_BC and CRF01_AE were the two main genotypes that created intricate molecular transmission networks. In order to prevent the expansion of molecular networks and stop the virus’s spread among MSM in Chongqing, more effective HIV intervention plans should be introduced.


Acquired immune deficiency syndrome (AIDS), as a serious sexually transmitted disease, is induced by the human immunodeficiency virus (HIV) and poses a threat to global population health [1]. MSM have been identified as an important population and continue to represent 23.0% of new HIV infections globally [2]. There were more than one million people living with HIV in China in 2020 [3]. The transmission of HIV among MSM has shown rapid growth in China, increasing from 2.5% of newly reported cases in 2006 to 26% in 2014 [4].

Chongqing is located in southwestern China, with a total population of over 30 million. With approximately 8000 cases reported in each of the previous three years, Chongqing is one of the most severely affected areas by HIV-1. According to Chongqing’s system of AIDS reporting, MSM made up more than 13.0% of newly diagnosed cases with HIV in this city in 2021, while more than 24.5% of cases in the downtown region (containing 9 districts) did. The high prevalence of HIV-1 among MSM was found in Southwest China. And the city of Chongqing had the greatest prevalence of HIV (13.8%, 95% CI: 12.8-14.9%) due to acceptance of homosexuality and open attitudes toward sex [5, 6]. The great mobility of the MSM population in China has been the main factor contributing to the spread of HIV among various places and groups [7, 8].

The emergence of a large range of novel circulating recombinant forms (CRFs) and unique recombinant forms (URFs) was rapid across several risk groups and regions [9]. The diversity of HIV-1 across the nation is a significant obstacle for the development of an HIV vaccine and disease prevention [10]. According to a cross-sectional national survey that included 4704 partial pol sequences, CRF01_AE (39.0%) was the most prevalent subtype of HIV-1 in 2017, followed by CRF07_BC (35.6%), CRF08_BC (8.9%), and subtype B (5.5%) [11]. CRF01_AE has had a sharp increase across almost all risk groups and geographical areas, and it has taken over among MSM in China [12]. In contrast, the earlier study found that CRF07_BC was the dominant subtype in Chongqing, China, accounting for nearly 70% of the newly diagnosed HIV infected MSM [13].

Few studies have focused on sexual behavior among MSM but have not systematically analyzed the molecular epidemiology of incident infections in Chongqing. As the HIV-1 epidemic grows among MSM, CRFs will increase and heighten the possibility of drug resistance, resulting in treatment failure [14]. Studying the epidemic of HIV-1 genetics and looking into TDR are crucial [15]. Here we use phylogenetic analyses to infer the distribution and characteristics of the subtypes of HIV-1, the TDR, and the cluster of molecularly transmitted viruses from 2018 to 2021. The results of this study could help the MSM population in Chongqing, China, with precise preventative and control strategies.


Study population and samples

From January 2018 to December 2021, 894 patients with HIV-1 infection were enrolled after informed consent was obtained according to the following criteria: (1) population of MSM; (2) adult residents over the age of 18 living in Chongqing; (3) The newly diagnosed HIV-infected individual and has not received treatment. Prior to the blood draw, participants’ demographic information was obtained via face-to-face interviews at two volunteer counseling and testing clinics in Jiangbei district and Yuzhong district.

Extraction, amplification, and sequencing

According to the manufacturer’s instructions, HIV-1 RNA was isolated from a blood sample using the Shuoshi Viral RNA Kit (SDK60105, China). This HIV-1 RNA pol region encodes regions of the reverse transcriptase and protease genes (HBX2: 2147–3462), which were amplified by Reverse Transcription-Polymerase Chain Reaction (RT-PCR) technique with commercial primers according to previously published methods [16]. On a 1% agarose gel, PCR product identification was carried out using electrophoresis. The sequencing of positive products was then entrusted to Beijing Genomics Research Center Ltd. in Beijing, China.

Sequence analysis and determination of genotypes

The software Sequencher 5.4.6 was used to modify and combine the generated sequence fragments. After that, BioEdit 7.0.9 was used to align the sequences and compare them to the reference sequence from LANL’s HIV Sequence Database. With the aid of COMET online analysis tool (, genotype identification was carried out. Additionally, using the neighbor-joining approach with the Tamura-Nei 93 model in MEGA X [17], a phylogenetic tree was created. URFs were defined as sequences that could not be supported by the phylogenetic tree and COMET. We utilized the online tool ITOL ( to enhance the visualization of the tree and assign distinct subtype colors for labeling.

Analysis of drug-resistance mutations

Three kind of drug classes (NRTIs, NNRTIs, and PIs) were evaluated for resistance using the online tool ( [18]. At least one or more TDR mutations were required for HIV-1 strains to be considered resistant.

Construction of the transmission network

The sequences were entered into Hyphy 2.2.4 for computing the pairwise genetic distance with the Tamura-Nei 93 mode. The genetic distance between molecular clusters was 1.5% [19]. Cytoscape v3.9.0 software was then used to process and visualize the gratings.

Statistical analysis

IBM SPSS version 22.0 program was used for all statistical analyses. We describe the qualitative statistics using frequentists. Univariate and multivariable logistic regression analyses were conducted to find potential risk factors. Statistical significance was defined as a P value < 0.05.


Social-demographic characteristics

A total of 894 successfully amplified HIV pol sequences were obtained from 1010 research participants, with a success rate of 88.5%. These people were 31.3 years old on average (ranging from 16 to 82 years). 98.0% (876/894) were Han ethnicity, and 96.2% (860/894) were permanent residents (> 2 years). Of these subjects, 78.3% (700/894) were unmarried, 14.1% (126/894) were married or living with a spouse, and 7.6% (68/894) were divorced or widowed. The majority of the subjects (59.6%, 533/894) had a junior college or higher level of education. (Table 1)

Sexuality-related characteristics

Among the 894 participants, 80.1% were homosexual, 15.9% were bisexual, and 4.0% were uncertain of their sexual orientation. The sex role survey showed that 34.3% (307/894) were the “Top”, 31.0% (277/894) were the “Bottom”, and 34.7% (310/894) were the “Verse”. 61.5% (550/894) had their first same-sex sex at the age of 18–24 years old. Most (81.4%, 728/894) seek sexual partners via social software. Combined rates of HIV with syphilis positivity and hepatitis C positivity were found to be 17.7% (158/894) and 0.3% (3/894) respectively.

HIV‑1 genotypes distribution

We found that CRF07_BC (73.6%, 658/894), CRF01_AE (19.6%, 175/894) were the predominant HIV-1 genotypes circulating in Chongqing, which accounted for 93.2% of the total number of.

infections. And CRF55_01B (3.4%, 30/894), CRF08_BC (1.1%, 10/894), and subtype B (1.0%, 9/894) were less commonly observed in this study. The number of genotypes ≤ 3 (68_01B and 67_01B) and Other/URFs, were accounted for 1.3% (12/894). The color of the circle on the outside of the branch with the genotypes of CRF07_BC, CRF01_AE, CRF08_BC, CRF55_01B, B, and URF/Other shown in yellow, red, green, blue, purple, and gray, respectively, is indicated in the phylogenetic tree. (Fig. 1)

Table 1 Characteristics of participants
Fig. 1
figure 1

HIV‑1 genotypes distribution in the phylogenetic tree

Characteristics of TDR

The prevalence of TDR in MSM in Chongqing was found to be 5.6% (50/894). Of the 894 strains, Q58E (1.5%, 13/894) and M46L (0.5%, 4/894) were the most common mutations associated with PI. K103N (1.1%, 10/894), Y188L (0.6%, 5/894) and E138A (0.3%, 3/894) mutations were the most common mutations associated with NNRTI. And three strains (0.3%, 3/894) had drug resistance mutations in both PIs and NNRTIs. (Table 2)

Table 2 The prevalence of transmitted drug resistance mutations sites

Patients with the CRF07_BC strain (2.4%, 16/658) had the highest odds of acquiring a mutation associated with PI, followed by those with the CRF01_AE strain (2.3%, 4/175). Patients with the CRF01_AE strain (2.3%, 4/175) had the highest odds of acquiring a mutation associated with NRTI, followed by those with the CRF07_BC strain (0.2%, 1/658). The patients with the CRF55_01B and CRF08_BC strains were the most likely to acquire a mutation associated with NNRTI, followed by those who were infected with the CRF01_AE strain (5.1%, 9/175) and the CRF07_BC strain (2.7%, 18/658). Overall, the TDR prevalence rate of the CRF07_BC genotype was 7.4%, that of genotype CRF01_AE was 5.0%, that of genotype CRF08_BC was 10.0%, and that of genotype CRF55_01B was 10.0%. (Table 3)

Table 3 Transmitted drug resistance mutations sites in different subtypes

Characteristics of molecular transmission network

Analysis of the molecular transmission network identified 48 clusters containing 600 sequences (67.1%, 600/894), with 18 clusters of CRF07_BC, 25 clusters of CRF01_AE, 2 clusters of CRF55_01B, and 2 clusters of other genotypes. The largest molecular group was the CRF07_BC genotype consisting of 481 individuals and 5928 edges. (Fig. 2)

In the multivariable logistic regression model, age 23–34 years (AOR, 0.674; 95% CI 0.455–0.997) or age 35–44 years (AOR, 0.562; 95% CI 0.318–0.995), uncertain sexual orientation (AOR, 0.350; 95% CI 0.154–0.797), with no syphilis (AOR, 0.531; 95% CI 0.336–0.839), and different genotypes were significantly associated with being in clusters compared to those who were not in clusters. (Table 4)

Fig. 2
figure 2

The molecular transmission network. The color of blue represents CRF07_BC, red represents CRF01_AE, orange represents CRF55_01B, and black represents others/URF; The round node indicates that there is a drug resistant mutation site, while the triangular node with no drug resistant mutation site

Table 4 Factors associated with transmission within clusters


We performed a cross-sectional molecular epidemiologic study of HIV-1 to follow the distribution and characteristics of HIV-1 subtypes, TDR, and molecularly transmitted clusters in newly diagnosed infections in MSM in Chongqing, China. CRF07_BC was found to be the most prevalent subtype, which was in agreement with previous studies in 2012–2014 [13]. However, according to the national survey, the HIV-1 CRF01_AE strain is the most prevalent strain in China, particularly among MSM [20]. It is very odd that Chongqing has become the city with one of the highest prevalence of the CRF07_BC genotype subtype in China.

The CRF07_BC is a variant of an HIV-1 B’ and C recombination in persons who inject drugs (PWID) in southwestern China during the 1990s, it also spread to other cities and groups of heterosexuals at first via the drug trafficking route [21]. Then, because of the aggressive fight against drug trafficking and the introduction of methadone maintenance treatment programs [22], CRF07_BC has increasingly spread through sexual contact among MSM [5, 23]. Thus, we hypothesize that HIV flowed into Xinjiang and Yunnan Province through the transmission of drug use in the early years and then spread to become the dominant strain in the Chongqing region. There were also a few studies that suggested that CRF07_BC strains are commonly detected in MSM in certain regions of China, such as Jilin [24] and Shijiazhuang[25].

In our study, in addition to CRF07_BC, CRF01_AE also represented 19.6% of all participants which was less than that of the survey conducted in Shanghai [26] and Tianjin [27]. CRF01_AE transmitted primarily through the sexual route spread from the southeastern coast and southwestern border to the entire country [28]. It is also noteworthy that the CRF55_01B strain represents a higher percentage than CRF08_BC and the B subtype, despite being a relatively “young” HIV strain in MSM. CRF55_01B was first reported in Changsha, Hunan province, and Dongguan, Guangdong Province, China in 2013 [29] and has since spread rapidly across the country [30]. In a recent study, CRF55_01B, despite originating from MSM, was found to be transmitted among heterosexuals, indicating that heterosexual males played a critical role in the transmission and dissemination of this strain [31].

This study suggests that DRMs influencing the efficacy of NNRTIs, followed by NRTIs and PIs that were in line with other cities in China [11]. TDF, 3TC, and EFV are the most frequently utilized drugs in China’s current first-line therapy guidelines [32]. First-line treatment medications can still be used constantly in Chongqing because the aforementioned three drug kinds demonstrated that resistance was largely at a low potential level. K103N observed as the most frequent mutation in response to NNRTIs was reported in Anhui [33] and Guangdong [34] provinces and is also consistent with the results of our study.

HIV molecular transmission network can help us formulate precise interventions and improve the efficiency of public health [35]. For the purposes of this study, we identify all clusters with the genetic distance of ≤ 1.5%, which has been shown to be more appropriate for molecular networks of more than three years [36]. Within these transmission clusters, we found an enormous CRF07_BC cluster comprising 422 sequences suggestive of a strong association in MSM in Chongqing. Due to this complex transmission relationship, once the dominant strain has been identified, it will still continue, which to some degree, accounts for the consistently high proportion of CRF07_BC within Chongqing.

We discovered that being in a cluster was highly correlated with age, sexual orientation, syphilis, and genotypes. Age is one such factor, which is best recognized because MSM in their younger years engage in greater sexual activity and have more complicated molecular transmission networks [37]. Syphilis is a sexually transmitted illness that increases the chance of contracting HIV and, to some extent, reflects the state of the sexual activity network [38].

The MSM who self-identified as homosexual engaged in sexual conduct more frequently than those who were less convinced of their sexual orientation [39]. As a result, different sexual orientations displayed varied rates of network entry. The two main genotype transmission clusters in the network, CRF07_BC and CRF01_AE, need to be dynamically monitored in order to allow for precise preventive and control actions to stop their transmission.


Despite the low overall prevalence of TDR, the significance of genotypic drug resistance monitoring needs to be emphasized. CRF07_BC and CRF01_AE were the two main genotypes that created intricate molecular transmission networks. In order to prevent the expansion of molecular networks and stop the virus’s spread among MSM in Chongqing, more effective HIV intervention plans should be introduced.

Data availability

The corresponding author can provide all of the datasets used in this study upon reasonable request.



Acquired immunodeficiency syndrome


Human immunodeficiency viruses


Men who have sex with men


Sexually transmitted diseases


Surveillance/transmitted drug resistance mutations


Transmitted drug resistance


Tamura-Nei 93


Unique recombinant forms


Antiretroviral therapy


Circulating recombinant forms


Nonnucleoside reverse transcriptase inhibitors


Nucleoside reverse transcriptase inhibitors


Protease inhibitors


Unique recombinant forms


  1. GBD 2017 HIV collaborators. Global, regional, and national incidence, prevalence, and mortality, of HIV. 1980–2017, and forecasts to 2030, for 195 countries and territories: a systematic analysis for the Global Burden of Diseases, Injuries, and Risk Factors Study 2017. Lancet HIV. 2019;6:e831–59.

  2. Yu M, Song D, Zhang T, Yao T, Chen Y, Liu Y, et al. High risks of HIV transmission for men sex worker - a comparison of profile and risk factors of HIV infection between MSM and MSW in China. BMC Public Health. 2022;22:858.

    Article  PubMed  PubMed Central  Google Scholar 

  3. He N. Research Progress in the epidemiology of HIV/AIDS in China. China CDC Wkly. 2021;3:1022–30.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Zheng S. The growing threat of China’s HIV epidemic. Lancet Public Health. 2018;3:e311.

    Article  PubMed  Google Scholar 

  5. Dong M-J, Peng B, Liu Z-F, Ye Q-N, Liu H, Lu X-L, et al. The prevalence of HIV among MSM in China: a large-scale systematic analysis. BMC Infect Dis. 2019;19:1000.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Ouyang L, Guo-Hui WU, Zhou Y, Zheng JQ, Zhou QH. Sex behaviors and HIV infection status among men who have sex with men in Chongqing,2010–2013. Chin J Dis Surveill. 2014;29:124–9.

    Google Scholar 

  7. Feng Y, Wu Z, Detels R, Qin G, Liu L, Wang X et al. HIV/STD prevalence among men who have sex with men in Chengdu, China and associated risk factors for HIV infection. J Acquir Immune Defic Syndr 1999. 2010;53 Suppl 1:S74-80.

  8. Wu J, Shen Y, Zhong P, Feng Y, Xing H, Jin L, et al. The predominant cluster of CRF01_AE circulating among newly diagnosed HIV-1-positive people in Anhui Province, China. AIDS Res Hum Retroviruses. 2015;31:926–31.

    Article  CAS  PubMed  Google Scholar 

  9. Li X, Li W, Zhong P, Fang K, Zhu K, Musa TH, et al. Nationwide Trends in Molecular Epidemiology of HIV-1 in China. AIDS Res Hum Retroviruses. 2016;32:851–9.

    Article  PubMed  Google Scholar 

  10. Zhang R, Dong TL, Liang WL, Cao ZB, Xie Z, Liu KM, et al. [Analysis of HIV-1 genetic subtype and pretreatment drug resistance among men who have sex with men infected with HIV-1 from 19 cities of 6 provinces in China]. Chin J Epidemiol. 2022;43:523–7.

    CAS  Google Scholar 

  11. Zhao S, Feng Y, Hu J, Li Y, Zuo Z, Yan J, et al. Prevalence of transmitted HIV drug resistance in antiretroviral treatment naïve newly diagnosed individuals in China. Sci Rep. 2018;8:12273.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Shi P, Chen Z, Meng J, Su M, Yang X, Fan W, et al. Molecular transmission networks and pre-treatment drug resistance among individuals with acute HIV-1 infection in Baoding, China. PLoS ONE. 2021;16:e0260670.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Zeng P, Liu Y, He M, Wang J, Keating S, Mao W, et al. The infection staging and profile of genotypic distribution and drug resistance mutation among the human immunodeficiency virus-1 infected blood donors from five chinese blood centers, 2012–2014. PLoS ONE. 2017;12:e0179328.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Yu F, Li Q, Wang L, Zhao H, Wu H, Yang S, et al. Drug Resistance to HIV-1 integrase inhibitors among treatment-naive patients in Beijing, China. Pharmacogenomics Pers Med. 2022;15:195–203.

    CAS  Google Scholar 

  15. Zhou C, Liang S, Li Y, Zhang Y, Li L, Ye L, et al. Characterization of HIV-1 molecular epidemiology and transmitted drug-resistance in newly diagnosed HIV-infected patients in Sichuan, China. BMC Infect Dis. 2022;22:602.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Chen Y, Shen Z, Feng Y, Ruan Y, Li J, Tang S, et al. HIV-1 subtype diversity and transmission strain source among men who have sex with men in Guangxi, China. Sci Rep. 2021 Apr;15(1):8319.

  17. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol 2018 Jun 1;35(6):1547–9.

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

    Article  CAS  Google Scholar 

  19. Little SJ, Kosakovsky Pond SL, Anderson CM, Young JA, Wertheim JO, Mehta SR, et al. Using HIV networks to inform real time prevention interventions. PLoS ONE. 2014 Jun;5(6):e98443.

  20. Han X, An M, Zhang M, Zhao B, Wu H, Liang S et al. Identification of 3 distinct HIV-1 founding strains responsible for expanding epidemic among men who have sex with men in 9 Chinese cities. J Acquir Immune Defic Syndr 1999. 2013;64:16–24.

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Yin W, Hao Y, Sun X, Gong X, Li F, Li J, et al. Scaling up the national methadone maintenance treatment program in China: achievements and challenges. Int J Epidemiol. 2010;39(Suppl 2):29–37.

    Google Scholar 

  23. Ge Z, Feng Y, Zhang H, Rashid A, Zaongo SD, Li K, et al. HIV-1 CRF07_BC transmission dynamics in China: two decades of national molecular surveillance. Emerg Microbes Infect. 2021;10:1919–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Li X, Zang X, Ning C, Feng Y, Xie C, He X, et al. Molecular epidemiology of HIV-1 in Jilin province, northeastern China: emergence of a new CRF07_BC transmission cluster and intersubtype recombinants. PLoS ONE. 2014;9:e110738.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Li L, Lu X, Li H, Chen L, Wang Z, Liu Y, et al. High genetic diversity of HIV-1 was found in men who have sex with men in Shijiazhuang, China. Infect Genet Evol J Mol Epidemiol Evol Genet Infect Dis. 2011;11:1487–92.

    Google Scholar 

  26. Li X, Xue Y, Lin Y, Gai J, Zhang L, Cheng H, et al. Evolutionary Dynamics and complicated genetic transmission network patterns of HIV-1 CRF01_AE among MSM in Shanghai, China. Sci Rep. 2016;6:34729.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Zheng M, Yu M, Cheng S, Zhou N, Ning T, Li L, et al. Characteristics of HIV-1 molecular transmission networks and drug resistance among men who have sex with men in Tianjin, China (2014–2018). Virol J. 2020;17:169.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Xing H, Liang H, Wan Z-Y, Chen X, Wei M, Ma P-F, et al. Distribution of recombinant human immunodeficiency virus type-1 CRF01_AE strains in China and its sequence variations in the env V3-C3 region. Chin J Prev Med. 2004;38:300–4.

    CAS  Google Scholar 

  29. Han X, An M, Zhang W, Cai W, Chen X, Takebe Y, et al. Genome sequences of a Novel HIV-1 circulating recombinant form, CRF55_01B, identified in China. Genome Announc. 2013;1:e00050–12.

    PubMed  PubMed Central  Google Scholar 

  30. Han X, Takebe Y, Zhang W, An M, Zhao B, Hu Q, et al. A large-scale survey of CRF55_01B from men-who-have-sex-with-men in China: implying the Evolutionary History and Public Health Impact. Sci Rep. 2015;5:18147.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Gan M, Zheng S, Hao J, Ruan Y, Liao L, Shao Y, 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Zhang F, Liu L, Sun M, Sun J, Lu H. An analysis of drug resistance among people living with HIV/AIDS in Shanghai, China. PLoS ONE. 2017;12:e0165110.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Zhang D, Wu J, Zhang Y, Shen Y, Dai S, Wang X, et al. Genetic characterization of HIV-1 epidemic in Anhui Province, China. Virol J. 2020;17:17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Lan Y, Li L, He X, Hu F, Deng X, Cai W, et al. Transmitted drug resistance and transmission clusters among HIV-1 treatment-naïve patients in Guangdong, China: a cross-sectional study. Virol J. 2021;18:181.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Zhang Z, Dai L, Jiang Y, Feng K, Liu L, Xia W, et al. Transmission network characteristics based on env and gag sequences from MSM during acute HIV-1 infection in Beijing, China. Arch Virol. 2017;162:3329–38.

    Article  CAS  PubMed  Google Scholar 

  36. Rose R, Lamers SL, Dollar JJ, Grabowski MK, Hodcroft EB, Ragonnet-Cronin M, et al. Identifying transmission clusters with cluster picker and HIV-TRACE. AIDS Res Hum Retroviruses. 2017;33:211–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Jia D, Zhao J, Liu Y, Wang X, Jia L, Gui T, et al. Two-year cross-sectional studies reveal that single, young MSMs in Shenzhen, China are at high risk for HIV infection. Virol J. 2019;16:83.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Zhang Y, Wu G, Lu R, Xia W, Hu L, Xiong Y, et al. What has changed HIV and syphilis infection among men who have sex with men (MSM) in Southwest China: a comparison of prevalence and behavioural characteristics (2013–2017). BMC Public Health. 2019;19:1314.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Frye V, Nandi V, Egan J, Cerda M, Greene E, Van Tieu H, et al. Sexual orientation- and race-based discrimination and sexual HIV risk behavior among urban MSM. AIDS Behav. 2015;19:257–69.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


The authors are grateful to Jiangbei District Center for Disease Control and Prevention staff and all individuals who volunteered to take part in the study.


This study was supported by Chongqing Natural Science Foundation Project (cstc2021jcyj-msxmX1171), Chongqing Talents Program for Innovative and Entrepreneurial Pioneers (cstc2021ycjh-bgzxm0097), the first batch of key disciplines on public health in Chongqing (YWBF2022072), Chongqing medical scientific research project of Chongqing Health Commission and Science and Technology Bureau (2022GDRC017). Chongqing medical scientific research project (2023MSXM113).

Author information

Authors and Affiliations



The funds were provided by Guohui Wu, Chao Zhou. The study was conceived by Houlin Tang, Fangfang Chen, Rongrong Lu, and Xu Xiang. The samples and data were collected by Chongyang Bai, Tianyu Tan, Long Li, Wei Zhang, Ling Ouyang, Jing Xu; The experiments were conducted by Chongyang Bai, Tianyu Tan; Tianyu Tan analyzed the data and wrote the article.

Corresponding authors

Correspondence to Houlin Tang or Guohui Wu.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethics approval and consent to participate

This study was approved by the Ethics Review Committee of Chongqing Center for Disease Control and Prevention [CQCDCLS (2021)010]. Written informed consent was signed by all participants before donating samples.

Consent for publication

Not applicable.

Additional information

Publisher’s Note

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

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tan, T., Bai, C., Lu, R. et al. HIV-1 molecular transmission network and drug resistance in Chongqing, China, among men who have sex with men (2018–2021). Virol J 20, 147 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: