Syphilis Testing as a Proxy Marker for a Subgroup of Men Who Have Sex With Men With a Central Role in HIV-1 Transmission in Guangzhou, China

Objectives: The objectives of this study were to distinguish the role of men who have sex with men (MSM) with or without syphilis testing in HIV-1 transmission and to provide molecular evidence of syphilis testing as a proxy marker for identifying the subgroup of MSM.

Methods: HIV-1 transmission clusters were constructed by HIV-TRACE and Cluster Picker using HIV-1 pol sequences from 729 newly diagnosed HIV-infected MSM from 2008 to 2012 in Guangzhou, China. The role of MSM in HIV-1 transmission networks was determined by a node influence measurement and centrality analysis. The association between syphilis testing and factors related to HIV-1 transmission and antiretroviral treatment (ART) were analyzed by the Cox regression model.

Results: Among HIV-infected MSM, 56.7% did not test for syphilis at the time of HIV-1 diagnosis. MSM without syphilis testing was a specific subgroup of MSM with a larger closeness centrality and clustering coefficient than the recipients of syphilis testing (P < 0.001), indicating their central position in the HIV-1 transmission networks. The median degree and radiality within HIV-1 transmission networks as well as the median K-shell scores were also greater for MSM without syphilis testing (P < 0.001), suggesting their relatively greater contribution in transmitting HIV-1 than the receipts of syphilis testing. MSM with syphilis testing usually did not disclose their occupation or were more likely to be unemployed or to take non-skilled jobs, to have a history of sexually transmitted infections (STIs), and to be AIDS patients when diagnosed with HIV-1 infection (P < 0.05). Multivariable Cox regression analysis indicated that syphilis testing per se did not promote the engagement of ART (P = 0.233) or affect the speed of CD4⁺ T cell count recovery after treatment (P = 0.256).

Conclusions: Our study identifies syphilis testing as a proxy marker of a specific subgroup of HIV-infected MSM who refuse syphilis testing during HIV-1 diagnosis with an important role in HIV-1 transmission. Specific prevention and intervention targeting MSM without syphilis testing during HIV-1 care are urgently needed.

Introduction

Over the past decades, syphilis has continued to spread globally, and disproportionately affects men who have sex with men (MSM) who are often co-infected with human immunodeficiency virus type one (HIV-1) (1–3). In light of the fact that syphilis and HIV-1 infection could facilitate the acquisition and transmission of each other, and even alter mutual disease progress (4–6), integrating syphilis testing into HIV-1 care has been recommended and promoted to curb the dual epidemic of syphilis and HIV-1 (7–9).

In China, despite the launch of a nationwide 10-year syphilis control plan in 2010 to promote the integration of HIV-1 and syphilis testing in the clinics of sexually transmitted infections (STIs) (10), there has been a continuous increase of syphilis cases and HIV-1/syphilis co-infections among MSM over the past decade (3, 11–13), due to low uptake of syphilis testing among HIV-infected MSM and other high-risk individuals in China (10, 13, 14). Furthermore, the current syphilis screening service might not fully cover all the risk groups. Hall et al. (13) have reported that MSM who did not engage in sex work were less likely to be tested for either HIV or syphilis. It is thus important to precisely determine the subgroup of MSM and to understand their difference in HIV-1 transmission (15). However, few studies have been published in China to investigate the role of MSM with or without syphilis testing in the control of HIV-1 epidemics (15), due to the lack of accurate and reliable information for MSM (16). Now, transmission network analysis using HIV-1 viral sequences has been revealed to be a valuable method to restore the social network of HIV-infected individuals in particular when epidemiological information, such as sexual behaviors, are not available or are inaccurate. The information obtained from HIV-1 transmission network topological structures may serve as a proxy for individuals with enhanced HIV-1 transmission risk (17, 18). For example, Ragonnet-Cronin et al. (19) adapted time-resolved phylogenetic analysis and HIV-1 transmission network analysis to identify 18.6% of all clustered heterosexual men as non-disclosed MSM in the UK. Based on HIV-1 transmission network analysis, we also found that non-disclosed MSM was a specific group and played a central role in HIV-1 transmission in China (20).

Therefore, in the current study, we adapted transmission network analysis to uncover the HIV-1 transmission risk among MSM with or without syphilis testing. We aimed to characterize the subgroup of MSM who refuse to test for syphilis and to provide molecular evidence of their important role in enhancing HIV-1 transmission in China.

Materials and Methods

Study Setting, Participants, and HIV-1 Sequences

Guangzhou is an international metropolis in the southern part of China, where the HIV-1 epidemic in MSM has been clustered for decades (21). The Center for Diseases Prevention and Control (CDC) is in charge of all HIV-1 cases in Guangzhou, including the diagnosis and confirmation of HIV-1 infection, CD4 cell measurement, and syphilis testing, etc. Eligibility criteria, recruitment, and measurement procedures for HIV-infected MSM have previously been described (22). Free consulting and testing for syphilis were available for all MSM at the time of HIV-1 testing or diagnosis, but they decided whether to take syphilis screening or not. Follow-up examinations were done every 3–6 months from June 2008 to February 2015. Demographic and laboratory testing data were collected from the medical charts of Guangzhou CDC. Blood samples were collected from 982 newly diagnosed and antiretroviral therapy (ART)-naïve HIV-1-infected MSM in Guangzhou CDC between January 2008 and December 2012. HIV-1 pol sequences (equal to the fragment of nucleotide 2,253–3,821 for the HIV-1 HXB2 strain) were successfully amplified by reverse transcription polymerase chain reaction (RT-PCR) and sequenced. The sequence data have been deposited in GenBank with the accession numbers listed in the Supplementary Material. Finally, a total of 729 (74.2%, 729/982) HIV-infected MSM were included in our study. The flow chart of the study is presented in Supplementary Figure 1. This study was approved by the Institutional Review Board of Guangzhou CDC (No. 2017030). Written informed consent was obtained from all the participants.

Phylogenetic Analysis

HIV-1 pol sequences were aligned using the reference sequences of HIV-1 group M and the circulating recombinant forms (CRFs) obtained from the Los Alamos HIV-1 database (http://www.hiv.lanl.gov). HIV-1 genotypes were determined by the phylogenetic trees constructed by IQ-TREE 1.6.9 with 1000 ultrafast bootstrap replicates and 1000 replicates of the Shimodaira–Hasegawa approximate likelihood-ratio test (SH-aLRT) (23).

Transmission Network Analysis

Transmission network analysis has been described previously (20). Briefly, all the HIV-1 sequences including a dataset of HIV-1 reference sequences from the Los Alamos National Laboratory HIV-1 database were aligned. The transmission clusters was identified by HIV Transmission Cluster Engine (HIV-TRACE) using a genetic distance threshold of 1.5% (24) or Cluster Picker with an intra-cluster genetic distance threshold of 4.5% (25, 26). The position of nodes in the transmission clusters were determined by nodal centrality indicators including degree, betweenness, closeness, clustering coefficient, and radiality, which were measured using NetworkAnalyzer 2.7 implemented in Cytoscape 3.7.0.13 (20, 27). K-shell (Ks) score, a measure of the cohesiveness of a subset of individuals among whom there are stronger, more direct, or more frequent ties than between other subgroups within the same network, was used to identify the hub node (28). Nodes with high Ks are more prone to infection and to be the most influential spreaders during epidemics (29–31). A novel graph theoretic clustering algorithm, Molecular Complex Detection (MCODE) implemented in Cytoscape 3.7.0.13, was adapted to identify the densely connected sub-networks in the large HIV-1 transmission clusters (32).

Statistical Analysis

To compare the difference of the centrality metrics and Ks score of nodes between HIV-infected MSM with or without syphilis testing, the Wilcoxon signed-rank test was adapted to estimate the statistical significance in all the HIV-infected MSM or in 1:1 matched pairs of MSM generated by the propensity score matching (PSM) method (33, 34). We compared their epidemiological data at HIV-1 diagnosis and further adapted the logistic regression analysis to explore the association between syphilis testing and factors related to HIV-1 transmission. All the candidate variables with a p-value of <0.1 in the univariate model were sequentially included in a binary multivariate model using the stepwise forward method to estimate the association with syphilis testing. Multivariate Cox proportional hazards regression analysis was used to identify the independent factors associated with the engagement of ART and the speed of ART-related CD4⁺ T cell recovery. SPSS version 25.0 (SPSS Inc., Chicago, IL, USA) and R software version 4.0.2 were used to perform the statistical analysis.

Results

MSM Without Syphilis Testing Were More Often Located in the Big HIV-1 Transmission Clusters

A total of 729 HIV-1 sequences obtained from newly diagnosed HIV-infected MSM were used to construct HIV-1 molecular transmission networks using HIV-TRACE. Among them, 638 HIV-1 sequences were located in 76 clusters (Supplementary Figure 2A) and 27 (35.5%) clusters contained only 2 individuals. The proportion of cluster sizes of 3–5, 6–10, and >10 was 30.3, 15.8, and 18.4%, respectively (Supplementary Figure 2B). For the small clusters with fewer than 5 members, 70–82% of them contained at least one recipient of syphilis testing while only 44–48% of these clusters contained MSM without syphilis testing. In contrast, for big clusters with more than 6 members, the proportion of the clusters that contained MSM without syphilis testing dramatically increased and further reached similar levels to the clusters containing the recipients of syphilis testing (Supplementary Figure 2B). Furthermore, in the small clusters, the subjects were dominated by MSM with syphilis testing (62–63%) while the number of MSM with or without syphilis testing were almost equal in the big clusters (Supplementary Figure 2C). These results indicated that MSM without syphilis testing were more often located in the big clusters rather than in the small clusters of HIV-1 transmission.

Furthermore, Cluster Picker analysis identified 45 HIV-1 transmission clusters (Table 1, Supplementary Figure 3). MSM infected with HIV-1 CRF55_01B, CRF 01_AE, and CRF 07_BC were more likely to be clustered than subtype B or other HIV-1 genotypes (p < 0.05, Table 1). However, the composition of HIV-1 genotypes and the demographic characteristics of the clusters identified by HIV-TRACE and Cluster Picker were not different (Table 1).

TABLE 1

Table 1. Demographic characteristics of HIV-infected MSM clustered by HIV-TRACE or Cluster Picker^&.

MSM Without Syphilis Testing Were More Efficient in HIV-1 Transmission

To further identify the key population driving HIV-1 transmission and to provide molecular evidence about the role of HIV-infected MSM without syphilis testing in HIV-1 transmission, we conducted centrality analysis and found that the following centrality indicators were significantly larger for MSM without syphilis testing than for the recipients of syphilis testing: the median degree (2.000, IQR 1.000–4.375 vs. 1.000, IQR 0.000–3.813, P = 0.004), the median radiality (1.000, 0.748–1.000 vs. 0.752, 0.000–1.000, P = 0.008), the clustering coefficient (0.645, 0.000–1.000 vs. 0.000, 0.000–0.792, P = 0.009), and the closeness centrality (1.000, 0.548–1.000 vs. 0.557, 0.000–1.000, P = 0.007), but no significant difference was observed for the betweenness centrality (P = 0.399, Table 2, Model 1). To reduce the potential bias due to the heterogeneity between the MSM with or without syphilis testing, we further analyzed the centrality indicators using 1:1 matched pairs of the MSM with or without syphilis testing based on PSM methodology to adjust the potential differences of the baseline characteristics (Supplementary Table 1). Similar results were obtained (Table 2, Model 2). These results indicated the central position of MSM without syphilis testing in the HIV-1 transmission networks and revealed their active role in HIV-1 transmission.

TABLE 2

Table 2. Centrality analysis of HIV-1 transmission networks for MSM with or without syphilis testing during 2008-2012 in Guangzhou, China.

Furthermore, we adapted a node influence measurement by calculating the Ks score of a node to identify the most influential spreaders in the HIV-1 transmission networks. Interestingly, we found that median Ks scores were statistically greater for the MSM without syphilis testing than those with syphilis testing (7.000, IQR 4.940–21.700 vs. 5.786, IQR 0.000–21.764, P = 0.022, Table 2, Figure 1). Then, we used MCODE analysis to analyze the densely connected sub-networks in the large transmission clusters and identified 33 key sub-networks and 8 seed nodes (Figure 1). Of note, the seed nodes were more frequently identified to be MSM without syphilis testing than those with syphilis screening (7/8 vs. 1/8) although the difference did not reach statistical significance (P = 0.069, Table 2, Model 1). These results reiterated that MSM without syphilis testing may be a specific subgroup of MSM with more efficient HIV-1 transmission.

FIGURE 1

Figure 1. Ks score-associated and group-associated genetic transmission networks for MSM with or without syphilis testing during 2008-2012 in Guangzhou, China. Ks score was measured by a K-shell decomposition algorithm. The colors and shapes represent different groups of the HIV-1 sequences from the MSM with syphilis screening (green circle) or without syphilis screening (red square) and the HIV-1 reference sequences (light yellow triangle). The density of the colors indicates the Ks score values. Seed is the highest scoring node in the densely connected sub-networks.

MSM With and Without Syphilis Testing Were Different Subgroups

Although free syphilis testing was provided upon the diagnosis of HIV-1 infection, in our study, 56.7% of HIV-infected MSM refused syphilis testing at the time of HIV-1 diagnosis. For those with or without syphilis testing, the significant differences were observed with regard to the occupation, the number of sex partners in the past 6 months, the history of STIs, and the stage of HIV-1 infection (Table 3). For example, 20.1% of the MSM without syphilis testing were unemployed or took non-skilled jobs, and 23.7% of them did not report their occupation while 77.5% of the recipients of syphilis testing had skilled/professional jobs or were students (P < 0.001). Furthermore, 42.4% of MSM without syphilis testing reported only one sex partner during the past 6 months whereas multiple sex partners were recorded in 88.9% of MSM with syphilis testing (p < 0.001). In addition, acquired immunodeficiency syndrome (AIDS) was recorded slightly more frequently in the MSM without syphilis testing than the recipients of syphilis testing (35.6 vs. 26.6%, p = 0.012). The percentage of past STIs was also different between MSM with or without syphilis testing (p < 0.001) in part due to the difference in regard to the awareness of their STIs (57.6 vs. 42.4%, p = 0.001, Supplementary Table 2). Multivariable logistic regression analysis also confirmed that MSM who did not report their occupations were most unlikely to receive syphilis testing [adjusted odds ratio (aOR) = 4.468, p < 0.001, Table 4]. In addition, syphilis testing was less likely to be accepted by those who disclosed their history of STIs than those who did not know their STIs (aOR = 0.466, p = 0.001, Table 4), Furthermore, MSM without syphilis testing were less likely to have multiple sex partners (p < 0.001, Table 4). Of note, the distribution of ethnicity (p = 0.129) and HIV-1 genotypes (p = 0.316) were not statistically different between MSM with and without syphilis testing (Table 3).

TABLE 3

Table 3. Comparison of demographic characteristics for HIV-infected MSM with or without syphilis screening from 2008 to 2012 in Guangzhou, China.

TABLE 4

Table 4. Factors associated with no syphilis testing among HIV-infected MSM from 2008 to 2012 in Guangzhou, China.

Syphilis Testing Was Not Associated With ART Initiation and ART-Related CD4 Cell Recovery

To investigate the difference between MSM with or without syphilis testing in the compliance of HIV-1 care, in particular the engagement of ART and the response to ART, we summarized the baseline characteristics between the subjects with or without ART (Supplementary Table 3), and found that the percentage of ART engagement was slightly higher in those without syphilis testing than the recipients of syphilis testing (69.4 vs. 55.5%, P = 0.001, Supplementary Table 3). Kaplan–Meier analysis also indicted that MSM without syphilis testing started ART earlier than the recipients of syphilis testing (P = 0.003, Figure 2A). To further explore the reasons for the difference of ART initiation, we adapted multivariate Cox proportional hazards regression analysis and found that syphilis testing was not an independent factor associated with ART initiation (adjusted p = 0.233, Table 5). The factors associated with the time to initiate ART included age, HIV-1 genotypes, the time of HIV-1 diagnosis, baseline CD4⁺ cell count, disease stage during HIV-1 diagnosis, and the disclosure of STIs (Table 5). Because STI history was associated with both syphilis testing with an unadjusted OR of 0.541 (Table 4) and ART initiation (unadjusted HR = 0.713, Table 5), we stratified the MSM population according to the disclosure of STIs in a time-to-event model to avoid confounding the STIs. The result indicated that ART engagement was not affected by syphilis testing in those aware (P = 0.546) or unaware of their STI history (P = 0.141, Supplementary Figure 4). These results indicated that syphilis testing per se did not affect the initiation of ART among HIV-infected MSM.

FIGURE 2

Figure 2. (A) Kaplan–Meier analysis of the time to initiate antiretroviral treatment between MSM with or without syphilis screening in Guangzhou, China, 2008–2012. The date of HIV-1 diagnosis was defined as time zero. The tick marks indicate the participants with censored data. P-values were calculated by the log-rank test. The median time to initiate ART was 2.7 years (range 2.3-3.0) for MSM without syphilis testing, and 3.6 years (range 3.1-4.1) for MSM with syphilis screening, respectively. MSM, men who have sex with men; HIV, human immunodeficiency virus. (B) Kaplan–Meier analysis of the time from CD4⁺ count <350–350 cells/μL following ART between MSM with or without syphilis screening. The date of ART initiation was defined as time zero. The tick marks indicate the participants with censored data. P-values were calculated by the log-rank test. The median time from CD4⁺ T-cell count <350–≥350 cells/μL was 9.1 months (range 4.3-13.9) for the patients without syphilis testing, and 5.1 months (range 1.8-8.5) for the patients with syphilis screening.

TABLE 5

Table 5. Factors associated with the time from HIV-1 diagnosis to antiretroviral treatment initiation between HIV-infected MSM with or without syphilis screening during 2008 and 2012 in Guangzhou, China.

Furthermore, the speed of CD4⁺ T cell recovery following ART was analyzed among the subjects with baseline CD4⁺ T cell count <350 cell/mm³ and at least two follow-up visit records. Among the 232 HIV-infected MSM analyzed, CD4⁺ cell count increased to ≥350 cells/mm³ in 84.8% (78/181) and 73.6% (103/181) of the subjects with or without syphilis testing, respectively. However, the difference was marginally significant (P = 0.044, Supplementary Table 4). Kaplan-Meier analysis indicated that HIV-infected MSM with syphilis testing achieved CD4⁺ cell count >350 cells/mm³ significantly earlier than the patients without syphilis testing. The median time was 5.1 months for the recipients of syphilis testing and 9.1 months for those without syphilis testing, respectively (P = 0.047, Figure 2B). Multivariate Cox proportional hazards regression analysis indicated that syphilis testing did not directly affect the outcome of CD4⁺ T cell recovery (P = 0.256, Table 6). In contrast, the speed of CD4⁺ T cell recovery was mainly associated with baseline CD4⁺ T cell counts at the time of HIV-1 diagnosis (P < 0.001) or ART initiation (P < 0.001, Table 6).

TABLE 6

Table 6. Factors associated with the time to CD4⁺ T-Cell recovery following initiation of ART between HIV-infected MSM with or without syphilis screening during 2008 and 2012 in Guangzhou, China.

Discussion

Syphilis screening is important for the control of both syphilis and HIV-1 infection. However, a large number of HIV-positive MSM do not receive syphilis testing as recommended (10, 14, 15). In our study, 56.7% of newly diagnosed HIV-infected MSM refused to take syphilis testing although it is provided for free in China. Our results are in line with the data obtained from HIV-infected MSM in Shanghai, China where the rate of syphilis testing was approximately 50% in 2010 (13). A recent online national survey in China also revealed that the acceptance rate of syphilis self-testing was 51.7% among MSM (35). Of note, the syphilis testing rate was suboptimal in the United States where only 68.3% of HIV-positive MSM took this testing during 2017–2018 (36). Though several studies have presented the socio-demographic characteristics of MSM with and without syphilis testing uptake (13, 35, 37), few studies investigated the impact of MSM who refuse to take syphilis testing on HIV-1 transmission and ART treatment partly because of the difficulty of identifying the subgroups of HIV-infected persons and precisely tracing their HIV-1 transmission.

Fortunately, HIV-1 genetic network analysis has been well-documented to be able to identify the subgroup of persons with specific significance in HIV-1 transmission, and to provide precise prevention measures targeting the specific sub-population (38, 39). In our study, we found that the MSM without syphilis testing were more likely to be found in the center of HIV-1 transmission networks and had a more intimate and direct relationship with other MSM according to their greater degree, radiality, and clustering coefficient within HIV-1 transmission networks. These results indicated that they were prone to contribute to the small world properties rather than in an intermediary form (40–42). Well-documented small world effect is characteristic of fast and efficient information transfer (40, 43). The presence of small world effect might be a mirror of efficient HIV-1 transmission among the HIV-infected MSM population who refuse to receive syphilis testing. Furthermore, coreness centrality obtained in the K-shell decomposition process is well-known to be a better measure than the degree to identify the potential influential spreaders in the transmission network (29–31). In our study, K-shell decomposition analysis also confirmed MSM who refuse syphilis testing as a specific subgroup of HIV-infected subjects with an active role in HIV-1 transmission. To our best knowledge, this is the first study to adapt HIV-1 transmission network analysis to characterize MSM without syphilis testing and to provide molecular evidence of their important role in enhancing HIV-1 transmission in China.

Our results indicated that such a specific group of MSM may lack knowledge about HIV-1 and syphilis or may not want to test for syphilis due to personal reasons (13). For example, in our study, MSM without syphilis testing were more likely to be unemployed or take non-skilled jobs, previous studies also indicated that non-testing of syphilis was more common among MSM with low HIV knowledge and lower income (13). Of note, we found that 49% of HIV-infected MSM did not know if they had STIs at the time of HIV-1 diagnosis. Several studies have shown that most persons with syphilis, including up to 44% of MSM, are unaware of their serostatus of syphilis as they may be asymptomatic for years (44, 45) or unwilling to disclose their STI diagnosis (46) and their sexual behaviors (15, 47). Our results also showed that MSM without syphilis screening were more likely to be highly connected in the transmission network, indicating that they may have more sexual partners. But unexpectedly, both a national online survey in China and our study found that MSM unwilling to screen for syphilis usually reported fewer sexual partners (35). These results suggested that the self-reported information from HIV-infected MSM may not be true. Therefore, we should interpret with caution the findings obtained through traditional epidemiological analysis according to the self-reported data. In contrast, the molecular epidemiological study was based on the HIV-1 sequences to construct the transmission networks. The results can be more objective and less likely to be affected by personal demographic information.

Our study is of great value for the precise prevention of a HIV-1 epidemic. (1) Our results further indicated that syphilis testing could be used as a proxy marker for the subgroup of MSM without syphilis testing during HIV-1 diagnosis. However, we would like to emphasize that syphilis testing per se was not associated with enhanced HIV-1 transmission and did not promote engagement of ART. As expected, we found that delayed diagnosis of HIV-1 and progressive disease stages were the key factors to affect the engagement of ART treatment. Furthermore, HIV-infected MSM without syphilis screening did not yield more from initiating ART earlier in our study. Consistent with the findings from previous studies, our study illustrated that the baseline CD4⁺ T cell count and early HIV-1 diagnosis were the major factors in predicting post-ART recovery of CD4⁺ cells (48), but not syphilis screening. (2) Molecular evolution and network analysis are important tools for the precise identification of HIV-1 transmission, which is usually difficult to define via traditional epidemiology surveys, due to the absence or inaccuracy of epidemiological data for HIV-infected individuals. (3) For the precise prevention and intervention of a HIV-1 epidemic, it is necessary to identify and distinguish the subgroups of HIV-infected subjects. We and others used a transmission network and cluster analysis to define potential non-disclosed MSM and their role in enhancing HIV-1 transmission.

Our study was subject to several limitations: (1) the study included 1 center, and was an observational study. Further investigations in multiple centers are needed to verify our findings. (2) The demographic information and the risk factors data such as the history of STIs, and the number of sexual partners in the past 6 months were self-reported although the information may not affect the phylogenetic analysis of HIV-1 sequences. (3) Syphilis screening data were captured from the HIV-1 primary care medical records only. If participants had gone elsewhere (i.e., hospital dermatology clinics) for syphilis screening, it would not be captured. (4) The construction of a HIV-1 transmission network based on only 1 HIV-1 pol gene may not be accurate enough and the data obtained were not authenticated using a wet lab procedure (20, 25).

Conclusions

HIV-1 transmission network analysis revealed that HIV-infected MSM who refuse to undertake syphilis testing at the time of HIV-1 diagnosis was a specific subgroup of MSM and played an important role in HIV-1 transmission. Therefore, syphilis testing may be a proxy marker for identifying the specific sub-population in HIV-1 transmission. Specific prevention and intervention targeting MSM without syphilis testing are urgently needed. Our results could contribute to evidence-based policy making for the precise intervention of the HIV-1 epidemic.

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 in the article/Supplementary Materials.

Ethics Statement

This study was approved by the Institutional Review Board of Guangzhou CDC (No. 2017030). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

ST and ZH conceived the study. HW, YL, QL and LH collected the data. LH and HY analyzed and interpreted the data. LH and ST wrote the manuscript. All authors reviewed, revised, and approved the final manuscript.

Funding

This work was supported by the National Key R&D Program of China (Project No. 2018ZX10732401-003) and the Collaborative Innovation Project for Health Care funded by the Bureau of Science and Information Technology of the Guangzhou Municipality (Project Nos. 201604020011, 201704020219, 201707010184, and 201803040002).

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.

Supplementary Material

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

Supplementary Figure 1. The flow chat of the study.

Supplementary Figure 2. The distribution of HIV-infected MSM with or without syphilis testing across and within the clusters of various sizes from 2008 to 2012 in Guangzhou, China. The HIV-1 transmission clusters identified by HIV-TRACE for MSM with or without syphilis testing during 2008–2012. The colors and shapes represent different groups of the HIV-1 sequences from the MSM with (green circle) or without (red square) syphilis testing and the HIV-1 reference sequences (light yellow triangle). The left and right y-axis are the total number of HIV-1 transmission clusters identified in this study and the proportion of these clusters, respectively. The x-axis shows the size of HIV-1 transmission clusters. The dusty blue bar represents the number of clusters while the lines indicate the clusters containing at least one HIV-infected MSM with (green) or without (orange) syphilis testing, respectively. The chart shows the coverage of syphilis testing among HIV-1 infected MSM clustered in 4 different cluster sizes. The blue and gray bar represent the proportion of HIV-1 infected MSM with and without syphilis testing in the corresponding clusters, respectively. MSM, men who have sex with men; HIV, human immunodeficiency virus.

Supplementary Figure 3. Phylogenetic transmission clusters for HIV-1 infected MSM in Guangzhou, China during 2008 and 2012. Phylogenetic HIV-1 transmission clusters identified by Cluster Picker. The outside circle displayed HIV-1 genotypes with different colors, i.e., HIV-1 CRF55_01B (yellow), CRF 01_AE (blue), CRF07_BC (green), subtype B (purple), and other HIV-1 genotypes (gray). The inside circle showed 46 transmission clusters with colored strips while the non-clustered sequences are displayed as white.

Supplementary Figure 4. Forest plot of the hazard ratio for the factors associated with antiretroviral initiation according to the disclosure of STIs history from 2008 to 2012 in Guangzhou, China. The forest plots summarize the hazard ratios (black squares) and their 95% confidential intervals (Cis, horizontal lines), as well as P-values for the interaction between ART initiation and the subgroup variables. A P-value < 0.05 is statistically significant. *P < 0.05; **P < 0.01; ***P < 0.001. MSM, men who have sex with men; STI, sexually transmitted infection; HIV, human immunodeficiency virus; AIDS, acquired immune deficiency syndrome.

Supplementary Table 1. Baseline characteristics of propensity-matched HIV-infected MSM diagnosed from 2008 to 2012 in Guangzhou, China.

Supplementary Table 2. Different percentage of STIs history between MSM with or without syphilis testing from 2008 to 2012 in Guangzhou, China.

Supplementary Table 3. Characteristics of HIV-infected MSM with or without antiviral treatment during 2008 and 2012 in Guangzhou, China.

Supplementary Table 4. Differences between ART-treated MSM with or without CD4⁺ T Cell recovery during 2008 and 2012 in Guangzhou, China.

References

1. Read P, Fairley CK, Chow EP. Increasing trends of syphilis among men who have sex with men in high income countries. Sex Health. (2015) 12:155–63. doi: 10.1071/SH14153

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Centers for Disease Control and Prevention. Sexually Transmitted Disease Surveillance, 2019. (2021). Atlanta, GA: Department of Health and Human Services. Available online at: https://www.cdc.gov/std/syphilis/stdfact-msm-syphilis.htm

Google Scholar

3. Wu Z, Xu J, Liu E, Mao Y, Xiao Y, Sun X, et al. HIV and syphilis prevalence among men who have sex with men: a cross-sectional survey of 61 cities in China. Clin Infect Dis. (2013) 57:298–309. doi: 10.1093/cid/cit210

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Buchacz K, Patel P, Taylor M, Kerndt P, Byers R, Holmberg S, et al. Syphilis increases HIV viral load and decreases CD4 cell counts in HIV-infected patients with new syphilis infections. AIDS. (2004) 18:2075–9. doi: 10.1097/00002030-200410210-00012

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Wu MY, Gong HZ, Hu KR, Zheng H-y, Wan X, Li J. Effect of syphilis infection on HIV acquisition: a systematic review and meta-analysis. Sex Transm Infect. (2020) 1–9. doi: 10.1136/sextrans-2020-054706

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Karp G, Schlaeffer F, Jotkowitz A, Riesenberg K. Syphilis and HIV co-infection. EurJ Internal Med. (2009) 20:9–13. doi: 10.1016/j.ejim.2008.04.002

CrossRef Full Text | Google Scholar

7. Tuite AR, Fisman DN, Mishra S. Screen more or screen more often? Using mathematical models to inform syphilis control strategies. BMC Public Health. (2013) 13:1–13. doi: 10.1186/1471-2458-13-606

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Fairley CK, Law M, Chen MY. Eradicating syphilis, hepatitis C and HIV in MSM through frequent testing strategies. Curr Opin Infect Dis. (2014) 27:56–61. doi: 10.1097/QCO.0000000000000020

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Tuite AR, Burchell AN, Fisman DN. Cost-effectiveness of enhanced syphilis screening among HIV-positive men who have sex with men: a microsimulation model. PLoS ONE. (2014) 9:e101240. doi: 10.1371/journal.pone.0101240

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Tucker JD, Cohen MS. China's syphilis epidemic: epidemiology, proximate determinants of spread, and control responses. Curr Opin Infect Dis. (2011) 24:50–5. doi: 10.1097/QCO.0b013e32834204bf

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Chow EPF, Wilson DP, Zhang L. HIV and syphilis co-infection increasing among men who have sex with men in China: a systematic review and meta-analysis. PLoS ONE. (2011) 6:e22768. doi: 10.1371/journal.pone.0022768

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Chen XS. Challenges in responses to syphilis epidemic. Lancet Infect Dis. (2017) 17:793–4. doi: 10.1016/S1473-3099(17)30327-4

CrossRef Full Text | Google Scholar

13. Hall CDX, Luu M, Nehl EJ, He N, Zheng T, Haardörfer R, et al. Do the same socio-demographic variables predict testing uptake and sero-status? HIV and syphilis among an observational sample of Chinese men who have sex with men in Shanghai, China. Int J STD AIDS. (2020) 31:939–49 doi: 10.1177/0956462420909733

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Ong JJ, Liao M, Lee A, Fu H, Pan SW, Tang W, et al. Bridging the HIV-syphilis testing gap: dual testing among men who have sex with men living in China. Sex Transm Infect. (2019) 95:251–3. doi: 10.1136/sextrans-2018-053527

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Ong JJ, Fu H, Pan S, Smith MK, Wu D, Wei C, et al. Missed opportunities for human immunodeficiency virus and syphilis testing among men who have sex with men in China: a cross-sectional study. Sex Transm Dis. (2018) 45:382–6. doi: 10.1097/OLQ.0000000000000773

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Shang H, Zhang L. MSM and HIV-1 infection in China. Nat Sci Rev. (2015) 2:388–91. doi: 10.1093/nsr/nwv060

CrossRef Full Text | Google Scholar

17. Campbell EM, Jia H, Shankar A, Hanson D, Luo W, Masciotra S, et al. Detailed transmission network analysis of a large opiate-driven outbreak of HIV infection in the United States. J Infect Dis. (2017) 216:1053–62. doi: 10.1093/infdis/jix307

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ferraz de Arruda G, Barbieri A, Rodríguez P, Moreno Y, da F. Costa L, Rodrigues F. Role of centrality for the identification of influential spreaders in complex networks. Phys Rev E. (2014) 90:032812. doi: 10.1103/PhysRevE.90.032812

CrossRef Full Text | Google Scholar

19. Ragonnet-Cronin M, Hué S, Hodcroft EB, Tostevin A, Dunn D, Fawcett T, et al. Non-disclosed men who have sex with men in UK HIV transmission networks: phylogenetic analysis of surveillance data. Lancet HIV. (2018) 5:e309–16. doi: 10.1016/S2352-3018(18)30062-6

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Yan H, He W, Huang L, Wu H, Liang Y, Li Q, et al. The central role of nondisclosed men who have sex with men in human immunodeficiency virus-1 transmission networks in Guangzhou, China. Open Forum Infect Dis. (2020) 7:ofaa154. doi: 10.1093/ofid/ofaa154

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Zhong F, Lin P, Xu H, Wang Y, Wang M, He Q, et al. Possible increase in HIV and syphilis prevalence among men who have sex with men in Guangzhou, China: results from a respondent-driven sampling survey. AIDS Behav. (2011) 15:1058–66. doi: 10.1007/s10461-009-9619-x

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Das A, Li J, Zhong F, Ouyang L, Mahapatra T, Tang W, et al. Factors associated with HIV and syphilis co-infection among men who have sex with men in seven Chinese cities. Int J STD AIDS. (2015) 26:145–55. doi: 10.1177/0956462414531560

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Trifinopoulos J, Nguyen L-T, von Haeseler A, Minh B. W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res. (2016) 44:gkw256. doi: 10.1093/nar/gkw256

PubMed Abstract | CrossRef Full Text

24. Wertheim J, Kosakovsky Pond S, Forgione L, Mehta S, Murrell B, Shah S, et al. Social and genetic networks of HIV-1 transmission in New York City. PLoS Pathog. (2017) 13:e1006000. doi: 10.1371/journal.ppat.1006000

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Ragonnet-Cronin M, Hodcroft E, Hué S, Fearnhill E, Delpech V, Brown AJ, et al. Automated analysis of phylogenetic clusters. BMC Bioinform. (2013) 14:317. doi: 10.1186/1471-2105-14-317

CrossRef Full Text

26. 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 Human Retroviruses. (2017) 33:211–8. doi: 10.1089/aid.2016.0205

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Wang J, Mo H, Wang F, Jin F. Exploring the network structure and nodal centrality of China's air transport network: a complex network approach. J Transp Geogr. (2011) 19:712–21. doi: 10.1016/j.jtrangeo.2010.08.012

CrossRef Full Text | Google Scholar

28. McGarry K, McDonald S. Complex network theory for the identification and assessment of candidate protein targets. Comput Biol Med. (2018) 97:113–23. doi: 10.1016/j.compbiomed.2018.04.015

PubMed Abstract | CrossRef Full Text | Google Scholar

29. De P, Singh AE, Wong T, Yacoub W, Jolly AM. Sexual network analysis of a gonorrhoea outbreak. Sex Transm Infect. (2004) 80:280–5. doi: 10.1136/sti.2003.007187

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Liu Y, Tang M, Zhou T, Younghae D. Core-like groups result in invalidation of identifying super-spreader by k-shell decomposition. Sci Rep. (2015) 5:9602. doi: 10.1038/srep09602

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Kitsak M, Gallos L, Havlin S, Liljeros F, Muchnik L, Stanley H, et al. Identification of influential spreaders in complex networks. Nat Phys. (2010) 6:888–93. doi: 10.1038/nphys1746

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Bader GD, Hogue CWV. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. (2003) 4:2. doi: 10.1186/1471-2105-4-2

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Rosenbaum PR, Rubin DB. The central role of the propensity score in observational studies for causal effects. Biometrika. (1983) 70:41–55. doi: 10.1093/biomet/70.1.41

CrossRef Full Text | Google Scholar

34. Walsh MC, Trentham-Dietz A, Newcomb PA, Gangnon R, Palta M. Using propensity scores to reduce case-control selection bias. Epidemiology. (2012) 23:772–3. doi: 10.1097/EDE.0b013e3182628365

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Wang C, Cheng W, Li C, Tang W, Ong JJ, Smith MK, et al. Syphilis Self-testing: a nationwide pragmatic study among men who have sex with men in China. Clin Infect Dis Offic Public Infect Dis Soc Am. (2020) 70:2178–86. doi: 10.1093/cid/ciz603

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Weiser J, Tie Y, Beer L, Pearson WS, Shouse RL. Receipt of prevention services and testing for sexually transmitted diseases among HIV-positive men who have sex with men, United States. Ann Internal Med. (2020) 173:162–4. doi: 10.7326/M19-4051

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Zhang T, Liu C, Han L, Tang W, Mao J, Wong T, et al. Community engagement in sexual health and uptake of HIV testing and syphilis testing among MSM in China: a cross-sectional online survey. J Int AIDS Soc. (2017) 20:21372. doi: 10.7448/IAS.20.01/21372

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Wang X, Wu Y, Mao L, Xia W, Zhang W, Dai L, et al. Targeting HIV prevention based on molecular epidemiology among deeply sampled subnetworks of men who have sex with men. Clin Infect Dis(2015) 61:1462–8. doi: 10.1093/cid/civ526

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Cao B, Saffer AJ, Yang C, Chen H, Peng K, Pan SW, et al. MSM behavior disclosure networks and HIV testing: an egocentric network analysis among MSM in China. AIDS Behav. (2019) 23:1368–74. doi: 10.1007/s10461-019-02404-z

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Watts DJ, Strogatz SH. Collective dynamics of ‘small-world’ networks. Nature. (1998) 393:440–2. doi: 10.1038/30918

CrossRef Full Text | Google Scholar

41. Achard S, Salvador R, Whitcher B, Suckling J, Bullmore E. A resilient, low-frequency, small-world human brain functional network with highly connected association cortical hubs. J Neurosci. (2006) 26:63–72. doi: 10.1523/JNEUROSCI.3874-05.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Xu T, Cullen KR, Mueller B, Schreiner MW, Lim KO, Schulz SC, et al. Network analysis of functional brain connectivity in borderline personality disorder using resting-state fMRI. Neuroimage Clin. (2016) 11:302–15. doi: 10.1016/j.nicl.2016.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Muldoon SF, Bridgeford EW, Bassett DS. Small-world propensity and weighted brain networks. Sci Rep. (2016) 6:22057. doi: 10.1038/srep22057

PubMed Abstract | CrossRef Full Text

44. CDC. The Center for Diseases Prevention and Control: prevalence and awareness of HIV infection among men who have sex with men — 21 cities, United States, 2008. MMWR Morb Mortal Wkly Rep. (2010) 59:1201–7.

Google Scholar

45. Denniston MM, Klevens RM, McQuillan GM, Jiles RB. Awareness of infection, knowledge of hepatitis C, and medical follow-up among individuals testing positive for hepatitis C: National Health and Nutrition Examination Survey 2001–2008. Hepatology. (2012) 55:1652–61. doi: 10.1002/hep.25556

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Dowshen N, Binns HJ, Garofalo R. Experiences of HIV-related stigma among young men who have sex with men. AIDS Patient Care STDS. (2009) 23:371–6. doi: 10.1089/apc.2008.0256

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Kumar N, Forastiere L, Zhang T, Yang F, Li KT, Tang W, et al. Lack of sexual behavior disclosure may distort STI testing outcomes. BMC Public Health. (2020) 20:616. doi: 10.1186/s12889-020-08768-5

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Stirrup OT, Copas AJ, Phillips AN, Gill MJ, Geskus RB, Touloumi G, et al. Predictors of CD4 cell recovery following initiation of antiretroviral therapy among HIV-1 positive patients with well-estimated dates of seroconversion. HIV Med. (2018) 19:184–94. doi: 10.1111/hiv.12567

PubMed Abstract | CrossRef Full Text | Google Scholar