Front. Cell. Infect. Microbiol. Frontiers in Cellular and Infection Microbiology Front. Cell. Infect. Microbiol. 2235-2988 Frontiers Media S.A. 10.3389/fcimb.2024.1391215 Cellular and Infection Microbiology Original Research Men with a history of commercial heterosexual contact play essential roles in the transmission of HIV-1 CRF55_01B from men who have sex with men to the general population in Guangxi, China Luo Tong 1 2 † Zhang Fei 1 2 † Liang Huayue 1 2 3 † Yu Dee 1 2 4 Cen Ping 1 2 Zhong Shanmei 1 2 Qin Cai 1 2 Yang Yuan 1 2 Jiang Jiaxiao 1 2 Liao Yanyan 1 2 Li Mu 1 2 Zhang Rongjing 1 2 Li Zeyu 1 2 Lin Zhifeng 1 2 Ye Li 1 2 * Liang Hao 1 2 * Liang Bingyu 1 2 * 1 Guangxi Key Laboratory of AIDS Prevention and Treatment, School of Public Health, Guangxi Medical University, Nanning, China 2 Collaborative Innovation Centre of Regenerative Medicine and Medical Bioresource Development and Application Co-constructed by the Province and Ministry, Life Science Institute, Guangxi Medical University, Nanning, China 3 Science and Technology Department, The First People’s Hospital of Qinzhou, Qinzhou, China 4 International School of Public Health and One Health, Hainan Medical University, Haikou, China

Edited by: Curtis Brandt, University of Wisconsin-Madison, United States

Reviewed by: Antoinette van der Kuyl, University of Amsterdam, Netherlands

Bin Zhao, The First Affiliated Hospital of China Medical University, China

Mengze Gan, Huazhong University of Science and Technology, China, in collaboration with reviewer BZ

*Correspondence: Bingyu Liang, liangbingyu@gxmu.edu.cn; Hao Liang, lianghao@gxmu.edu.cn; Li Ye, yeli@gxmu.edu.cn

†These authors have contributed equally to this work

 23 08 2024 2024 14 1391215 25 02 2024 02 08 2024 Copyright © 2024 Luo, Zhang, Liang, Yu, Cen, Zhong, Qin, Yang, Jiang, Liao, Li, Zhang, Li, Lin, Ye, Liang and Liang 2024 Luo, Zhang, Liang, Yu, Cen, Zhong, Qin, Yang, Jiang, Liao, Li, Zhang, Li, Lin, Ye, Liang and Liang

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

 Background

There is increasing focus on HIV-1 CRF55_01B in China. However, there is limited information regarding the dissemination of CRF55_01B across different regions and populations in Guangxi. This study was performed to elucidate the evolutionary history of the introduction and dissemination of CRF55_01B in Guangxi.

 Methods

Molecular network and phylogenetic analyses were used to investigate the transmission characteristics of CRF55_01B in China. The analyses particularly focused on the cross-provincial spatial and temporal transmission patterns between Guangdong Province and Guangxi, as well as the transmission dynamics among different regions and populations within Guangxi.

 Results

In total, 2226 partial pol sequences of CRF55_01B strains sampled from 2007 to 2022 were collected, including 1895 (85.09%) sequences from Guangdong, 199 (8.94%) sequences from Guangxi, and 172 (7.59%) sequences from other provinces of China. Most people living with HIV in Guangxi were infected with HIV-1 through heterosexuals (52.76%). Among these, 19.10% had a history of commercial heterosexual contact (CHC) and 15.58% had a history of non-marital non-commercial heterosexual contact (NMNCHC). Overall, 1418 sequences were identified in the molecular network. Notably, the sequences from Guangdong Province were most closely linked to those from Guangxi. Phylogenetic analysis showed that CRF55_01B was first introduced from Shenzhen City to Nanning City around 2007. Subsequently, CRF55_01B established local transmission within Guangxi, with Nanning City serving as the transmission center from 2008 to 2017. After 2017, the CRF55_01B strain spread to other regions of Guangxi. Men who have sex with men (MSM) and men with a history of CHC have played a significant role in the transmission of CRF55_01B among different populations in Guangxi.

 Conclusions

This study provides evidence on the transmission trajectory of CRF55_01B among different regions and populations in Guangxi. Given the bridging role of men with a history of CHC in the dissemination of CRF55_01B from MSM to the general population, it is imperative to enhance surveillance among key populations to mitigate the secondary transmission of HIV-1.

 HIV-1 CRF55_01B MSM phylogenetic Bayesian section-in-acceptance Virus and Host Introduction

By the end of 2020, the number of people living with HIV (PLWH) in China was 1.22 million, with the vast majority infected with circulating recombinant forms (CRFs) (He, 2021). China is among the countries with a high prevalence of CRFs worldwide (Bbosa et al., 2019). According to the Los Alamos National Laboratory (LANL) HIV Database (https://www.hiv.lanl.gov/), more than 30 CRFs had been identified in China by the end of 2022. The high prevalence of CRFs will undoubtedly increase the complexity and challenges of HIV/AIDS prevention and control.

CRF55_01B is a novel HIV strain resulting from the recombination of CRF01_AE and B subtype variants. It was first identified among men who have sex with men (MSM) in Shenzhen City, China (Han et al., 2013). Although CRF55_01B has only been known for a decade and accounts for only 2.06% of the total number of HIV-1 infections in China (Wang et al., 2021), it has rapidly evolved and is considered the most important prevalent CRF by the Chinese Center for Disease Control and Prevention (Zhang et al., 2013), aside from CRF07_BC and CRF08_BC (Wei et al., 2021). Previous studies have indicated that the CRF55_01B strain originated in 2004 and was subsequently disseminated among MSM in 2013 (Zhao et al., 2014; Zai et al., 2020; Gan et al., 2021). Because of the rapidly developing railway system in China (Gan et al., 2021), the CRF55_01B strain has spread throughout the country through MSM, with a prevalence ranging from 1.5% to 12.5% in different regions (Wei et al., 2021). It is postulated that individuals infected with CRF55_01B may exhibit higher viral loads than those infected with CRF07_BC and CRF01_AE, and they may remain asymptomatic for an extended period because of a slower decline of their CD4+ T-cell count (Wei et al., 2021). Therefore, it is probable that CRF55_01B will exhibit more extensive transmission than CRF07_BC and CRF01_AE. CRF55_01B has reportedly spread from MSM to the heterosexual (HET) population and is progressively spreading within HET communities (Gan et al., 2021).

With the evolution of the HIV/AIDS epidemic, China has introduced two new categories of a non-marital heterosexual sexual contact (NMHC) history since 2014: commercial and non-commercial. In other words, a history of NMHC can now be further differentiated into a history of commercial heterosexual contact (CHC) or non-marital non-commercial heterosexual contact (NMNCHC). Previous research has suggested that HIV-1 transmission can occur among individuals with different infection routes. For instance, CRF55_01B can be transmitted from MSM to HETs, and CRF07_BC-N can be transmitted from MSM to injecting drug users and HETs (Gan et al., 2021; Gan et al., 2022). Few studies to date have been conducted to explore the relationships between HIV-1 transmission and diverse sexual contact histories.

Guangxi, located in the southwest region of China and sharing a border with Guangdong Province, was among the initial destination provinces to experience the transmission of CRF55_01B from Guangdong Province to other parts of China (Zai et al., 2020; Gan et al., 2021). Currently, there is no report on the transmission characteristics of CRF55_01B in different regions and populations with different sexual contact histories in Guangxi. In this study, data from the LANL HIV Database and the Guangxi Key Laboratory of AIDS Prevention and Control were used for molecular network and phylogenetic analyses to elucidate the transmission pathway of CRF55_01B in different regions and populations in Guangxi. The aim of this research is to provide valuable insights into local HIV/AIDS prevention and control efforts.

 Materials and methods Study population and sequences down load

From 1 January 2014 to 1 January 2023, we extracted viral nucleic acid from blood samples of HIV-infected individuals after their annual CD4+ T-cell count test. We then amplified the HIV pol partial region and sequenced the HIV pol sequence. Nearly 90% of the infected individuals in the sampled city undergo a CD4 + T-cell count test every year, allowing us to obtain these CRF55_01B pol sequences from 90% of the local HIV-infected individuals. Additional CRF55_01B pol sequences were obtained from the HIV database of the LANL. Each sequence represented an individual infected with HIV-1.

Information on the sample year, sample city, sex, and transmission route were collected for sequences downloaded from the LANL HIV Database. For the collected sequences, data on the sample year, sample city, sex, marital status, and sexual contact history were obtained from the National Free Antiretroviral Treatment Program Database. HET sexual contact history includes five categories: HC, MHC, NMHC, CHC, and NMNCHC. These are depicted in Figure 1 . Female sex workers were identified by staff from the local Center for Disease Control and Prevention.

Five categories of heterosexual contact history.

All sequences were aligned to the HXB2 (2253-3312nt, 1060 base pairs) reference sequence using the online tool HIVAlign from the LANL HIV Database via the multiple alignment fast Fourier transform method. Recombination analysis was performed using the online tools jpHMM and RIP, which are also available from the LANL HIV Database. The maximum likelihood (ML) phylogenetic tree was constructed using FastML v3.1.1. Duplicate sequences, sequences with >5% mixed bases, and sequences with less than 1060 base pairs were excluded.

 HIV-1 molecular network reconstruction

The genetic distance (GD) between two sequences was calculated using the Tamura-Nei 93 evolutionary model. The HIV-TRACE Transmission Cluster Engine (HIV-TRACE) was then employed to construct the molecular networks of CRF55_01B (Kosakovsky Pond et al., 2018), in accordance with the recommended GD threshold of 0.5% by the United States Centers for Disease Control and Prevention (National Center for HIV/AIDS VH, STD, and TB Prevention and Division of HIV/AIDS Prevention, 2018). The linkage degree is defined as the number of links between different regions in a molecular network.

 Spatiotemporal transmission dynamics of CRF55_01B

Three subsets were selected from the entire sequence dataset (N = 2226) for Bayesian analysis. The first dataset was designed to determine the origin and introduction time of CRF55_01B into Guangxi. It consisted of 271 sequences from various provinces in China. The second dataset included sequences from both Guangxi and Guangdong Province and was designed to investigate the origin and spatiotemporal dynamics of CRF55_01B transmission, with geographical information specified down to the city level. The sequences in the first and second datasets were filtered based on representative samples from different years and regions within the branches of the ML tree constructed using FastTree. The third dataset (n = 138) exclusively contained sequences from Guangxi and was utilized to determine the transmission trajectories of CRF55_01B among individuals with different sexual contact histories. The sequences selected for this dataset were chosen based on their representation across different years and sexual contact histories identified in the ML tree. The drug resistance mutation sites were edited and excluded before the Bayesian analysis.

The ML trees were constructed using FastML v3.1.1 (Ashkenazy et al., 2012). Outlier sequences were then examined and excluded by TempEst v1.5.3 (Rambaut et al., 2016). To identify the origin and transmission time of CRF55_01B in Guangxi, a Bayesian asymmetric discrete phylogeographic model in BEAST v1.10.4 was employed (Suchard et al., 2018). A general time-reversible substitution model was used to enhance the migration rates and nucleotide substitution rates. The uncorrelated relaxed clock and the Bayesian skyline were employed as the tree prior. Markov chain Monte Carlo analysis was performed on the empirical tree for 200 million iterations, with sampling every 20,000 steps. This process was run twice independently in BEAST. Markov chain Monte Carlo analysis was also employed to estimate the most recent common ancestor (tMRCA). The outputs were analyzed using Tracer (Rambaut et al., 2018), and effective sampling size (ESS) values of >200 indicated sufficient convergence. LogCombiner was used to merge the log.txt and trees.txt files from the two runs to obtain the optimal maximum clade credibility (MCC) tree. After excluding the first 10% of runs as burn-in, the maximum clade credibility trees were summarized by TreeAnnotator v1.10.4 (Suchard et al., 2018) and then visualized in FigTree v1.7.2. The spatiotemporal dissemination of CRF55_01B was presented using a Bayesian stochastic search variable selection (BSSVS) procedure (Lemey et al., 2009). Bayes factors (BFs) (Lemey et al., 2009) were used to quantify statistical support, and the results were summarized with the widely-used tool SpreaD3 (Bielejec et al., 2016). In this study, BF < 3 was considered anecdotal evidence, 3 < BF ≤ 10 as substantial evidence, 10 < BF ≤ 30 as strong evidence, 30 < BF ≤ 100 as very strong evidence, and BF > 100 as decisive evidence.

 Results Socio-demographic characteristics

In total, 2226 CRF55_01B sequences were collected from 2007 to 2022. Of these, 2077 sequences were obtained from the LANL HIV Database, while the remaining 189 sequences were amplified in the present study. These sequences covered 19 provinces in China. Most of the sequences were obtained from Guangdong Province (1894, 85.09%), followed by Guangxi (199, 8.94%). Most sequences lacked information on the transmission route (56.38%, 1255/2226); 34.46% (767/2226) were from MSM. Among the sequences from Guangxi, most individuals had been infected with HIV-1 through HET (52.76%, 105/199), of whom 19.10% had a CHC history and 15.58% had an NMNCHC history ( Table 1 ). More than half of the CRF55_01B-infected individuals in Guangxi came from Nanning City and Guigang City. A total of 90.45% of individuals infected with CRF55_01B were men, 69.85% had received more than 9 years of education, and 70.35% were divorced/widowed/unmarried. Further details are provided in Supplementary Table S1 and Supplementary Table S2 .

Demographic information of individuals infected with CRF55_01B in Guangxi.

Variables Total, n (%)
Sex
 Male 180 (90.45)
 Female 19 (9.55)
Ethnicity
 Han 101 (50.75)
 Other 98 (49.25)
Education, years
 ≤9 33 (16.58)
 >9 139 (69.85)
 Unknown 27 (13.57)
Marital status
 Married/cohabitated 59 (29.65)
 Divorced/widowed/unmarried 140 (70.35)
Sample year
 2009–2013 7 (3.52)
 2014–2016 3 (1.51)
 2017–2019 89 (44.72)
 2020–2022 100 (50.25)
Sexual contact history
 MSM 94 (47.24)
 HC 28 (14.07)
 MHC 5 (2.51)
 NMHC 3 (1.51)
 CHC 38 (19.10)
 NMNCHC 31 (15.58)
Sample location within Guangxi Province, region (cities)
 North (LZ) 12 (6.03)
 East (HZ, WZ, YL) 8 (4.02)
 Central (NN, GG) 114 (57.29)
 South (QZ, BH) 28 (14.07)
 West (BS, CZ, HC) 37 (18.59)
Total 199 (100)

MSM, men who have sex with men; HC, heterosexual contact; MHC, marital heterosexual contact; NMHC, non-marital heterosexual sexual contact; CHC, commercial heterosexual contact; NMNCHC, non-marital non-commercial heterosexual contact; LZ, Liuzhou; HZ, Hezhou; WZ, Wuzhou; YL, Yulin; NN, Nanning; GG, Guigang; QZ, Qinzhou; BH, Beihai; BS, Baise; CZ, Chongzuo; HC, Hechi.

 Molecular network of CRF55_01B in China

In total, 1418 CRF55_01B sequences were incorporated into the molecular network at a GD of 0.5%, forming 180 molecular clusters with a range of 2 to 715 nodes ( Supplementary Figure S1 ). The largest molecular cluster consisted of 715 nodes and contained sequences from 11 provinces. As shown in Figure 2 , the sequences from Guangdong Province played a pivotal role within the network by directly connecting with sequences from the other 11 provinces. Notably, the sequences from Guangxi, Beijing, Hebei Province, and Shanghai were highly correlated with those from Guangdong Province. The molecular network indicated that sequences from Guangdong Province, Hebei Province, and Shanghai were linked to Guangxi. Notably, sequences from Guangdong Province demonstrated the highest linkage degree with sequences from Guangxi.

Linkage degree between different provinces in the molecular network. The color of the line segment corresponds to the linkage degree, defined as the number of degrees in a molecular network.

 Spatiotemporal transmission characteristics of CRF55_01B

Phylogeographic analysis of the CRF55_01B sequences in China showed that this strain originated from Guangdong Province and that the estimated tMRCA was approximately 2003 ( Supplementary Figure S2 ). The CRF55_01B strain was first introduced into Guangxi in 2007, and it then spread to other provinces after 2010. Guangdong Province served as the gateway for the CRF55_01B strain to spread to other regions of China ( Figure 3 ). The regions in a transitional relationship with Guangxi were Guangdong Province, Shanghai, Anhui Province, and Zhejiang Province. Sequences from Guangdong Province and Guangxi exhibited a high degree of connectivity, with a BF of 74564 and a posterior probability of 1.0. Further details are shown in Supplementary Table S3 .

Transmission direction of CRF55_01B strain in China. The arrows represent the direction of transmission. The color of the line represents the size of the Bayes factor.

 Transmission dynamics across regions of Guangxi and/or Guangdong

To elucidate the origin and transmission dynamics of CRF55_01B in Guangxi, we conducted distinct phylogenetic analyses utilizing sequences from Guangxi and Guangdong Province. The findings from the phylogeographical inference revealed a complex history of viral migration, highlighting numerous interconnections between Guangdong Province and Guangxi as well as intricate interactions across diverse regions within Guangxi ( Figure 4 ).

Migration events of CRF55_01B between regions in Guangdong and Guangxi through time: (A) Migration events before 2008. (B) Migration events between 2008 and 2012. (C) Migration events between 2012 and 2017. (D) and after 2017. The arrows represent the direction of transmission. The color of the line represents the size of the Bayes factor.

These results indicate that the transmission of CRF55_01B was from Shenzhen City to Nanning City in 2007, substantiated by a BF of 221243 and a posterior probability of 1.0 ( Figure 4A ). Subsequently, the spread of CRF55_01B from Nanning City to other cities in Guangxi commenced in 2008 ( Figure 4B ), reaching its peak frequency of transmission between 2012 and 2017 ( Figure 4C ). After 2017, evidence of mutual viral transmission between distinct regions within Guangxi emerged. Notable migration events were supported by strong evidence, including transmissions from Nanning City to Chongzuo City, from Nanning City to Liuzhou City, and from Nanning City to Hechi City ( Figure 4D ). Nanning City emerged as a central hub facilitating the dissemination of CRF55_01B within Guangxi. Detailed data are provided in Supplementary Table S4 .

 Propagation characteristics of CRF55_01B among individuals with different sexual contact histories

A Bayesian model was employed to elucidate the correlation between the propagation of CRF55_01B and different sexual contact histories. As depicted in Figure 5 , MSM emerged as the primary population driving the dissemination of CRF55_01B in Guangxi. The findings of this study indicate a high probability that CRF55_01B was transmitted from MSM to men with an NMNCHC history (BF=278.93, posterior probability=0.97), to men with a CHC history (BF=16707.25, posterior probability=1.0), and to men with an HC history (BF = 1276.60, posterior probability = 0.99).

Transmission characteristics among groups with different sexual contact histories. The arrows represent the direction of transmission. The color of the line represents the size of the Bayes factor. NMHC.F, female with a history of non-marital heterosexual contact; CHC.M, men with a history of commercial heterosexual contact; HC.F, female with a history of heterosexual contact; MHC.F, female with a history of marital heterosexual contact; NMHC.M, men history of with non-marital heterosexual contact; HC.M, men with a history of heterosexual contact; NMNCHC.F, female with a history of non-marital non-commercial heterosexual contact.

Furthermore, our analysis revealed that women played a passive role in the transmission of CRF55_01B in Guangxi. Most infections in women could be traced back to men. Specifically, transmission to women with an NMNCHC history occurred through men with an HC history (BF = 1276.60, posterior probability = 0.99). Similarly, transmission to women with an NMHC history occurred through men with a CHC history (BF = 35.53, posterior probability = 0.79). Detailed data are provided in Supplementary Table S5 .

 Discussion

In this study, we elucidated for the first time the spatiotemporal transmission dynamics of CRF55_01B in different regions of Guangxi and the transmission patterns among populations with different sexual contact histories using a molecular network and phylogeographic and phylodynamic methods. Our findings demonstrated that CRF55_01B is extensively disseminated throughout Guangxi. The most probable source of transmission is Guangdong Province, with an estimated tMRCA of 2009. We determined that MSM and men with a history of CHC have played a key role in the dissemination of CRF55_01B in Guangxi.

We found that the CRF55_01B strain disseminated from MSM to other populations through diverse transmission routes, reflecting the intricate sexual contact relationship among individuals infected with CRF55_01B. Given that MSM and men with a CHC history have been identified as high-risk groups for second-generation transmission of HIV-1, it is imperative to strengthen the monitoring for MSM and men with a CHC history.

In addition, this study observed that CRF55_01B originated in 2003 from Guangdong Province and disseminated into Guangxi in 2007, which is consistent with previous findings (Zai et al., 2020; Gan et al., 2021). CRF55_01B was first identified in Guangdong, where it resulted in an outbreak and disseminated to other regions in China (Han et al., 2013; Zhao et al., 2014; Zai et al., 2020; Gan et al., 2021). Guangdong has the largest migrant population in China. Given that Guangxi shares a border with Guangdong Province and is the largest migrant population that has moved to Guangdong (Network, 2023), Guangdong is likely the primary source of the CRF55_01B strain prevalent in Guangxi. This study did not find evidence indicating that the CRF55_01B strain has spread from Guangxi to other provinces in China.

The Bayesian analysis suggested that CRF55_01B was first introduced into Nanning City from Shenzhen City in 2007 and subsequently spread from Nanning City to other cities within Guangxi. It is not uncommon for Chinese residents to travel to near metropolises for trade or work-related activities. Nanning City is a preferred destination for Guangxi residents engaged in trade or work because it serves as the economic and political center of Guangxi. Additionally, Nanning City is the region in Guangxi with the most severe HIV/AIDS epidemic (Lan et al., 2021). These factors have accelerated the dissemination of the CRF55_01B strain from Nanning City to other urban centers in Guangxi.

Previous studies have shown that MSM represent a key group affected by CRF55_01B infection (Zhao et al., 2014; Zai et al., 2020; Gan et al., 2021). However, our findings indicate that more than half of individuals infected with the CRF55_01B strain had a history of HC in Guangxi. Additionally, CHC and NCNMHC exhibited the highest prevalence within the HC group. There are two potential reasons for these findings. First, the primary mode of HIV-1 transmission among infected individuals in Guangxi is HC, particularly CHC and NCNMHC (Dong et al., 2021). MSM transmission accounts for only 7% of HIV-1 cases, which is significantly lower than the national average of 23% in China and 60% in Guangdong Province (Daily, 2016; He, 2021; Lan et al., 2021). Second, some men might conceal their same-sex behaviors because of a perceived higher level of discrimination, thus reporting a history of HC. Therefore, the proportion of individuals infected with CRF55_01B through HC in Guangxi might be overestimated.

In this study, we constructed a transmission pathway of the CRF55_01B strain among individuals with different sexual contact histories. The findings indicated that MSM played a predominant role in transmitting the CRF55_01B strain to other populations in Guangxi. Considering that the CRF55_01B strain originated from MSM and spread outward through this population (Gan et al., 2021), this finding is not unexpected. It is well acknowledged that MSM represent a high-risk group for the transmission of HIV-1. Previous studies have indicated that MSM in China exhibit lower rates of condom use, engage in multiple sexual partners, and might be identified as bisexual (Zhao et al., 2021). Because same-sex marriage is illegal in China, a significant number of MSM may engage in HET activity and marry women (i.e., MSMW). MSMW serve as a ‘bridge’ for the transmission of HIV-1 from MSM, a high-risk group, to women, a low-risk group (Tao et al., 2013; Oster et al., 2015). Furthermore, research indicates that compared with MSM only, MSMW are more likely to engage in unprotected and commercial sex (Tao et al., 2013; Wang et al., 2015). Therefore, it is essential to actively promote respect and social equality for MSM, with a particular focus on MSMW.

Our study also revealed that women represent the terminal end of the CRF55_01B transmission chain, whereas men engaged in commercial HET activities represent the main source of transmission of the CRF55_01B strain to women. This suggests a significant contribution of CHC to the dissemination of CRF55_01B in Guangxi. The primary mode of HIV-1 transmission among infected individuals in Guangxi is through CHC (Dong et al., 2021). Men with a history of CHC have a higher risk of HIV-1 infection and transmission (Kaufman, 2011). Our earlier research showed that men with CHC were responsible for the majority of local HIV-1 transmissions and were prone to transmitting the virus to other populations (Chen et al., 2021). Many patrons are reluctant to use male condoms because of difficulties with erection or discomfort, and this limits the promotion of condom use in sexual transactions. Prospective studies suggest that the provision of female condoms can reduce the incidence of unprotected sexual intercourse (Beksinska et al., 2019; Zhang et al., 2020). Therefore, promoting the use of female condoms may be an effective strategy for the prevention of HIV-1 transmission.

This study had several limitations. First, a relatively short sequence was used for the phylogenetic analysis. Shorter sequences cannot represent the dynamic evolution of the whole sequence, and this may have led to some bias of the results. Second, the non-Guangxi sequences were derived from the LANL HIV Database, which has been uploaded by other researchers. Therefore, these sequences may have been subject to publication bias, potentially leading to underestimation of the dissemination of CRF55_01B sequence in China. Third, the exclusive reliance on self-reported sexual contact histories may have led to recall bias and social desirability bias. Some individuals who reported having a history of HC were unable to further specify which type of HC they had (MHC, NMNCHC, or NMCHC). This might have diminished the reliability of the results regarding the transmission relationships among individuals with different exposure histories. Fourth, the accuracy of our results may have been affected by the quantity and quality of the sequence database. We downloaded all available pol sequences of the CRF55_01B strain, incorporated all CRF55_01B sequences sampled by our laboratory, and performed strict quality control prior to analysis.

 Conclusion

This study provides the first estimation of the dissemination of the CRF55_01B strain from high-risk groups to the general population in Guangxi. MSM and men with CHC serve as a ‘bridge’ for the transmission of CRF55_01B, and studies of these populations provide valuable insights into prevention and control efforts. Furthermore, Nanning City was found to be the center of the dissemination and proliferation of the CRF55_01B strain to other cities in Guangxi. Considering the complex relationship between HIV-infected individuals with different sexual contact histories in Guangxi, it is imperative to adopt a multifaceted approach to mitigate the second-generation transmission of CRF55_01B, especially among MSM and men with a history of CHC.

 Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/genbank/, 2hzDLVrHHT.

 Ethics statement

The studies involving humans were approved by Human Research Ethics Committee (IRB) of Guangxi Medical University (Ethical Review No.2019-SB-102). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

 Author contributions

TL: Writing – original draft. FZ: Writing – original draft. HYL: Writing – original draft. DY: Methodology, Software, Writing – original draft. PC: Visualization, Writing – original draft. SZ: Resources, Writing – original draft. CQ: Resources, Writing – original draft. YY: Resources, Writing – original draft. JJ: Resources, Writing – original draft. YL: Validation, Writing – original draft. ML: Investigation, Writing – original draft. RZ: Investigation, Writing – original draft. ZYL: Investigation, Writing – original draft. ZFL: Data curation, Writing – original draft. LY: Writing – review & editing. HL: Writing – review & editing. BL: Writing – review & editing.

 Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Natural Science Foundation of China (Grant Nos. 82103899, 82060610), the Guangxi Natural Science Foundation (Grant No. 2022JJA141110), the National Key R&D Program of China (Grant No. 2022YFC2305001), the Thousands of Young and Middle Age Key Teachers Training Program in Guangxi Colleges and Universities (to BL), and the Guangxi Bagui Young Top Scholar (to BL).

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2024.1391215/full#supplementary-material

References Ashkenazy H. Penn O. Doron-Faigenboim A. Cohen O. Cannarozzi G. Zomer O. . (2012). FastML: a web server for probabilistic reconstruction of ancestral sequences. Nucleic Acids Res. 40, W580–W584. doi: 10.1093/nar/gks498 Bbosa N. Kaleebu P. Ssemwanga D. (2019). HIV subtype diversity worldwide. Curr. Opin. HIV AIDS. 14, 153–160. doi: 10.1097/COH.0000000000000534 Beksinska M. Zulu B. Nkosi P. Mabude Z. Smit J. (2019). Offering choice in female condoms: The South African experience. Int. J. STD AIDS. 30, 985–990. doi: 10.1177/0956462419853376 Bielejec F. Baele G. Vrancken B. Suchard M. A. Rambaut A. Lemey P. (2016). SpreaD3: interactive visualization of spatiotemporal history and trait evolutionary processes. Mol. Biol. Evol. 33, 2167–2169. doi: 10.1093/molbev/msw082 Chen X. Qin C. Chen R. Huang Y. Xu Y. Tang Q. . (2021). Epidemiological profile and molecular genetic characterization of HIV-1 among female sex workers and elderly male clients in Guangxi, China. Emerg. Microbes Infect. 10, 384–395. doi: 10.1080/22221751.2021.1888659 Daily N. (2016). Guangdong: Rapid increase in HIV transmission among men who have sex with men (Central Government of the People's Republic of China). Available at: https://www.gov.cn/xinwen/2016-11/30/content_5140382.htm. Dong Z. Pan X. Cai C. Qin Q. Gao G. F. Lyu F. (2021). Spatial analysis of people living with HIV/AIDS transmitted through commercial heterosexual contact or non-marital non-commercial heterosexual contact - China, 2018. China CDC Wkly. 3, 316–319. doi: 10.46234/ccdcw2021.086 Gan M. Zheng S. Hao J. Ruan Y. Liao L. Shao Y. . (2021). The prevalence of CRF55_01B among HIV-1 strain and its connection with traffic development in China. Emerg. Microbes Infect. 10, 256–265. doi: 10.1080/22221751.2021.1884004 Gan M. Zheng S. Hao J. Ruan Y. Liao L. Shao Y. . (2022). Spatiotemporal patterns of CRF07_BC in China: A population-based study of the HIV strain with the highest infection rates. Front. Immunol. 13, 824178. doi: 10.3389/fimmu.2022.824178 Han X. An M. Zhang W. Cai W. Chen X. Takebe Y. . (2013). Genome Sequences of a Novel HIV-1 Circulating Recombinant Form, CRF55_01B, Identified in China. Genome Announc 1. doi: 10.1128/genomeA.00050-12 He N. (2021). Research progress in the epidemiology of HIV/AIDS in China. China CDC Wkly. 3, 1022–1030. doi: 10.46234/ccdcw2021.249 Kaufman J. H. I. V. (2011). sex work, and civil society in China. J. Infect. Dis. 204 Suppl 5, S1218–S1222. doi: 10.1093/infdis/jir538 Kosakovsky Pond S. L. Weaver S. Leigh Brown A. J. Wertheim J. O. (2018). HIV-TRACE (TRAnsmission cluster engine): a tool for large scale molecular epidemiology of HIV-1 and other rapidly evolving pathogens. Mol. Biol. Evol. 35, 1812–1819. doi: 10.1093/molbev/msy016 Lan G. Shen Z. Ge X. Zhu Q. Zhou Y. Liang S. . (2021). Review and prospect of comprehensive AIDS prevention, control and management in Guangxi. Chin. J. New Clin. Med. 14, 951–955. doi: 10.3969/j.issn.1674-3806.2021.10.01 Lemey P. Rambaut A. Drummond A. J. Suchard M. A. (2009). Bayesian phylogeography finds its roots. PLoS Comput. Biol. 5, e1000520. doi: 10.1371/journal.pcbi.1000520 National Center for HIV/AIDS VH, STD, and TB Prevention Division of HIV/AIDS Prevention (2018). DETECTING AND RESPONDING TO HIV transmission clusters A GUIDE FOR HEALTH DEPARTMENTS (National Center for HIV/AIDS, Viral Hepatitis, STD, and TB Prevention, Division of HIV/AIDS Prevention). Network G. S. I. (2023). In 2022, Guangdong's permanent population will continue to rank first in the country and grow steadily 2023. Available online at: http://stats.gd.gov.cn/tjkx185/content/post_4147054.html. (Accessed: June 30, 2023). Oster A. M. Wertheim J. O. Hernandez A. L. Ocfemia M. C. Saduvala N. Hall H. I. (2015). Using molecular HIV surveillance data to understand transmission between subpopulations in the United States. J. Acquir. Immune Defic. Syndr. 70, 444–451. doi: 10.1097/QAI.0000000000000809 Rambaut A. Drummond A. J. Xie D. Baele G. Suchard M. A. (2018). Posterior summarization in bayesian phylogenetics using tracer 1.7. Syst. Biol. 67, 901–904. doi: 10.1093/sysbio/syy032 Rambaut A. Lam T. T. Max Carvalho L. Pybus O. G. (2016). Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol. 2, vew007. doi: 10.1093/ve/vew007 Suchard M. A. Lemey P. Baele G. Ayres D. L. Drummond A. J. Rambaut A. (2018). Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 4, vey016. doi: 10.1093/ve/vey016 Tao J. Ruan Y. Yin L. Vermund S. H. Shepherd B. E. Shao Y. . (2013). Sex with women among men who have sex with men in China: prevalence and sexual practices. AIDS Patient Care STDS. 27, 524–528. doi: 10.1089/apc.2013.0161 Wang H. Y. Xu J. J. Zou H. C. Reilly K. H. Zhang C. M. Yun K. . (2015). Sexual risk behaviors and HIV infection among men who have sex with men and women in China: evidence from a systematic review and meta-analysis. BioMed. Res. Int. 2015, 850132. doi: 10.1155/2015/850132 Wang X. Zhang Y. Liu Y. Li H. Jia L. Han J. . (2021). Phylogenetic analysis of sequences in the HIV database revealed multiple potential circulating recombinant forms in China. AIDS Res. Hum. Retroviruses 37, 694–705. doi: 10.1089/aid.2020.0190 Wei L. Li H. Lv X. Zheng C. Li G. Yang Z. . (2021). Impact of HIV-1 CRF55_01B infection on the evolution of CD4 count and plasma HIV RNA load in men who have sex with men prior to antiretroviral therapy. Retrovirology. 18, 22. doi: 10.1186/s12977-021-00567-z Zai J. Liu H. Lu Z. Chaillon A. Smith D. Li Y. . (2020). Tracing the transmission dynamics of HIV-1 CRF55_01B. Sci. Rep. 10, 5098. doi: 10.1038/s41598-020-61870-x Zhang H. Hsieh E. Wang L. Liao S. (2020). HIV/AIDS among female sex workers in China: epidemiology and recent prevention strategies. Curr. HIV/AIDS Rep. 17, 151–160. doi: 10.1007/s11904-019-00477-y Zhao J. Cai W. Zheng C. Yang Z. Xin R. Li G. . (2014). Origin and outbreak of HIV-1 CRF55_01B among MSM in Shenzhen, China. J. Acquir. Immune Defic. Syndr. 66, e65–e67. doi: 10.1097/QAI.0000000000000144 Zhang L. Chow E. P. Jing J. Zhuang X. Li X. He M. . (2013). HIV prevalence in China: integration of surveillance data and a systematic review. Lancet Infect. Dis. 13, 955–963. doi: 10.1016/S1473-3099(13)70245-7 Zhao P. Yang Z. Li B. Xiong M. Zhang Y. Zhou J. . (2021). Simple-to-use nomogram for predicting the risk of syphilis among MSM in Guangdong Province: results from a serial cross-sectional study. BMC Infect. Dis. 21, 1199. doi: 10.1186/s12879-021-06912-z