Molecular network-based intervention brings us closer to ending the HIV pandemic

Evolving molecular HIV clusters revealed genotype-specific dynamics in Guangzhou, China (2008-2020)

OBJECTIVES

This study investigated the genotype-specific dynamics of molecular HIV clusters (MHCs) in Guangzhou, China, aiming to enhance HIV control.

METHODS

HIV pol sequences from people with HIV (PWH) in Guangzhou (2008-2020) were obtained for genotyping and molecular network creation. MHCs were identified and categorized into three types: emerging, growing, or stable. Clustering rates, proportions of cluster types, and members within each type were calculated and their trends were assessed using joinpoint regression.

RESULTS

Among 8395 PWH, the most prevalent HIV-1 genotypes were CRF07_BC (39.7%) and CRF01_AE (32.6%). The genotype composition has been stable since 2012 (Ps > 0.05). The overall clustering rate was 43.3%, with significant variations across genotypes (P < 0.001), indicating genotype-specific transmission fitness. Significant declines in overall and genotype-specific clustering rates toward the end of 2020 (Ps < 0.05), potentially offer support for HIV control efforts in reducing local infections. The continuously increasing proportions of stable clusters and the gradually decreasing proportions of emerging and growing clusters (either Ps < 0.05 or Ps > 0.05) suggest a trend toward stable molecular network structure. However, growing clusters exhibited CRF55_01B, CRF07_BC, and CRF59_01B dominance that indicate their priority for interventions.

CONCLUSION

The evolving MHCs highlight the genotype-specific cluster dynamics, providing fresh insights for enhanced prevention and control strategies.

Interpreting the Epidemiological Characteristics of HIV-1 in Heterosexually Transmitted Population Based on Molecular Transmission Network in Kunming, Yunnan: A Retrospective Cohort Study

Heterosexuals have become the most prevalent group of HIV-1 in Kunming, Yunnan Province. Utilizing the principle of genetic similarity between their gene sequences, we built a molecular transmission network by gathering data from earlier molecular epidemiological studies. This allowed us to analyze the epidemiological features of this group and offer fresh concepts and approaches for the prevention and management of HIV-1 epidemics. Cytoscope was used to visualize and characterize the network following the processing of the sample gene sequences by BioEdit and HyPhy. The number of possible links and the size of the clusters were investigated as influencing factors using a zero-inflated Poisson model and a logistic regression model, respectively. A scikit-learn-based prediction model was developed to account for the dynamic changes in the HIV-1 molecular network. Six noteworthy modular clusters with network scores ranging from 4 to 9 were found from 150 clusters using Molecular Complex Detection analysis at a standard genetic distance threshold of 0.01. The size of the number of possible links and the network's clustering rate were significantly impacted by sampling time, marital status, and CD4+ T lymphocytes (all p < 0.05). The gradient boosting machine (GBM) model had the highest area under the curve value, 0.884 ± 0.051, according to scikit-learn. Though not all cluster subtypes grew equally, the network clusters were relatively specific and aggregated. The largest local transmission-risk group for HIV-1CRF08_BC is now the heterosexual transmission population. The most suitable model for constructing the HIV-1 molecular network dynamics prediction model was found to be the GBM model.

Combining molecular transmission network analysis and spatial epidemiology to reveal HIV-1 transmission pattern among the older people in Nanjing, China

BACKGROUND

In China, the problem of HIV infection among the older people has become increasingly prominent. This study aimed to analyze the pattern and influencing factors of HIV transmission based on a genomic and spatial epidemiological analysis among this population.

METHODS

A total of 432 older people who were aged ≥ 50 years, newly diagnosed with HIV-1 between January 2018 and December 2021 and without a history of ART were enrolled. HIV-1 pol gene sequence was obtained by viral RNA extraction and nested PCR. The molecular transmission network was constructed using HIV-TRACE and the spatial distribution analyses were performed in ArcGIS. The multivariate logistic regression analysis was performed to analyze the factors associated with clustering.

RESULTS

A total of 382 sequences were successfully sequenced, of which CRF07_BC (52.3%), CRF01_AE (32.5%), and CRF08_BC (6.8%) were the main HIV-1 strains. A total of 176 sequences entered the molecular network, with a clustering rate of 46.1%. Impressively, the clustering rate among older people infected through commercial heterosexual contact was as high as 61.7% and three female sex workers (FSWs) were observed in the network. The individuals who were aged ≥ 60 years and transmitted the virus by commercial heterosexual contact had a higher clustering rate, while those who were retirees or engaged other occupations and with higher education degree were less likely to cluster. There was a positive spatial correlation of clustering rate (Global Moran I = 0.206, P < 0.001) at the town level and the highly aggregated regions were mainly distributed in rural area. We determined three large clusters which mainly spread in the intra-region of certain towns in rural areas. Notably, 54.5% of cases in large clusters were transmitted through commercial heterosexual contact.

CONCLUSIONS

Our joint analysis of molecular and spatial epidemiology effectively revealed the spatial aggregation of HIV transmission and highlighted that towns of high aggregation were mainly located in rural area. Also, we found vital role of commercial heterosexual contact in HIV transmission among older people. Therefore, health resources should be directed towards highly aggregated rural areas and prevention strategy should take critical persons as entry points.

Characteristics of Subtype and Molecular Transmission Networks among Newly Diagnosed HIV-1 Infections in Patients Residing in Taiyuan City, Shanxi Province, China, from 2021 to 2023

The HIV-1 pandemic, spanning four decades, presents a significant challenge to global public health. This study aimed to understand the molecular transmission characteristics of newly reported HIV infections in Taiyuan, Shanxi Province, China, to analyze the characteristics of subtypes and the risk factors of the transmission network, providing a scientific basis for precise prevention and intervention measures. A total of 720 samples were collected from newly diagnosed HIV-1 patients residing in Taiyuan between 2021 and 2023. Sequencing of partial genes of the HIV-1 pol gene resulted in multiple sequence acquisitions and was conducted to analyze their subtypes and molecular transmission networks. Out of the samples, 584 pol sequences were obtained, revealing 17 HIV-1 subtypes, with CRF07_BC (48.29%), CRF01_AE (31.34%), and CRF79_0107 (7.19%) being the dominant subtypes. Using a genetic distance threshold of 1.5%, 49 molecular transmission clusters were generated from the 313 pol gene sequences. Univariate analysis showed significant differences in the HIV transmission molecular network in terms of HIV subtype and household registration (p < 0.05). Multivariate logistic regression analysis showed that CRF79_0107 subtype and its migrants were associated with higher proportions of sequences in the HIV transmission network. These findings provide a scientific foundation for the development of localized HIV-specific intervention strategies.

HIV-1 Molecular Networks and Pretreatment Drug Resistance at the Frontier of Yunnan Province, China

The border areas of Yunnan Province in China are severely affected by human immunodeficiency virus (HIV). To investigate the risk of HIV transmission and assess the prevalence of pretreatment drug resistance (PDR) in the border area, blood samples were collected from individuals with newly reported HIV in 2021 in three border counties (Cangyuan, Gengma, and Zhenkang) in Yunnan Province. Among the 174 samples successfully genotyped, eight circulating recombinant forms (CRFs), two subtypes, and several unique recombinant forms (URFs) were identified. CRF08_BC (56.9%, 99/174), URFs (14.4%, 25/174), CRF01_AE (10.9%, 19/174), and CRF07_BC (8.0%, 14/174) were the main genotypes. CRF08_BC and URFs were detected more frequently in Chinese and Burmese individuals, respectively. CRF07_BC was found more frequently in men who have sex with men. The proportion of individuals detected in HIV-1 networks was only associated with case-reporting counties. When stratified by county, individuals aged ≤40 years in Cangyuan and ≥41 years in Gengma were more likely to be found in these networks. Furthermore, 93.8% (15/16) of the links in Cangyuan and 79.4% (50/63) of those in Gengma were located within their own counties. The prevalence of PDR to any antiretroviral drug, nucleoside reverse transcriptase inhibitors (NRTIs), and non-nucleoside reverse transcriptase inhibitors (NNRTIs) were 10% (17/170), 0.6% (1/170), and 9.4% (16/170), respectively. The most frequent resistance-associated mutations (RAMs) were V179D/VD/E/T (22.9%, 39/170) and E138A/G/K/R (13.5%, 23/170). In the molecular networks, six clusters shared common RAMs. HIV-1 genetics has become more diverse in border areas. HIV-1 molecular network analysis revealed the different characteristics of the HIV-1 epidemic in the border counties. The prevalence of PDR showed an upward trend, and the PDR to NNRTIs was close to the public response threshold. These findings provide information for the development of AIDS prevention and treatment strategies.

Molecular transmission network analysis of newly diagnosed HIV-1 infections in Nanjing from 2019 to 2021

OBJECTIVE

The objective of this study was to conduct a comprehensive analysis of the molecular transmission networks and transmitted drug resistance (TDR) patterns among individuals newly diagnosed with HIV-1 in Nanjing.

METHODS

Plasma samples were collected from newly diagnosed HIV patients in Nanjing between 2019 and 2021. The HIV pol gene was amplified, and the resulting sequences were utilized for determining TDR, identifying viral subtypes, and constructing molecular transmission network. Logistic regression analyses were employed to investigate the epidemiological characteristics associated with molecular transmission clusters.

RESULTS

A total of 1161 HIV pol sequences were successfully extracted from newly diagnosed individuals, each accompanied by reliable epidemiologic information. The analysis revealed the presence of multiple HIV-1 subtypes, with CRF 07_BC (40.57%) and CRF01_AE (38.42%) being the most prevalent. Additionally, six other subtypes and unique recombinant forms (URFs) were identified. The prevalence of TDR among the newly diagnosed cases was 7.84% during the study period. Employing a genetic distance threshold of 1.50%, the construction of the molecular transmission network resulted in the identification of 137 clusters, encompassing 613 nodes, which accounted for approximately 52.80% of the cases. Multivariate analysis indicated that individuals within these clusters were more likely to be aged ≥ 60, unemployed, baseline CD4 cell count ≥ 200 cells/mm3, and infected with the CRF119_0107 (P < 0.05). Furthermore, the analysis of larger clusters revealed that individuals aged ≥ 60, peasants, those without TDR, and individuals infected with the CRF119_0107 were more likely to be part of these clusters.

CONCLUSIONS

This study revealed the high risk of local HIV transmission and high TDR prevalence in Nanjing, especially the rapid spread of CRF119_0107. It is crucial to implement targeted interventions for the molecular transmission clusters identified in this study to effectively control the HIV epidemic.

Exploring Dynamic Changes in HIV-1 Molecular Transmission Networks and Key Influencing Factors: Cross-Sectional Study

BACKGROUND

The HIV-1 molecular network is an innovative tool, using gene sequences to understand transmission attributes and complementing social and sexual network studies. While previous research focused on static network characteristics, recent studies' emphasis on dynamic features enhances our understanding of real-time changes, offering insights for targeted interventions and efficient allocation of public health resources.

OBJECTIVE

This study aims to identify the dynamic changes occurring in HIV-1 molecular transmission networks and analyze the primary influencing factors driving the dynamics of HIV-1 molecular networks.

METHODS

We analyzed and compared the dynamic changes in the molecular network over a specific time period between the baseline and observed end point. The primary factors influencing the dynamic changes in the HIV-1 molecular network were identified through univariate analysis and multivariate analysis.

RESULTS

A total of 955 HIV-1 polymerase fragments were successfully amplified from 1013 specimens; CRF01_AE and CRF07_BC were the predominant subtypes, accounting for 40.8% (n=390) and 33.6% (n=321) of the specimens, respectively. Through the analysis and comparison of the basic and terminal molecular networks, it was discovered that 144 sequences constituted static molecular networks, and 487 sequences contributed to the formation of dynamic molecular networks. The findings of the multivariate analysis indicated that the factors occupation as a student, floating population, Han ethnicity, engagement in occasional or multiple sexual partnerships, participation in anal sex, and being single were independent risk factors for the dynamic changes observed in the HIV-1 molecular network, and the odds ratio (OR; 95% CIs) values were 2.63 (1.54-4.47), 1.83 (1.17-2.84), 2.91 (1.09-7.79), 1.75 (1.06-2.90), 4.12 (2.48-6.87), 5.58 (2.43-12.80), and 2.10 (1.25-3.54), respectively. Heterosexuality and homosexuality seem to exhibit protective effects when compared to bisexuality, with OR values of 0.12 (95% CI 0.05-0.32) and 0.26 (95% CI 0.11-0.64), respectively. Additionally, the National Eight-Item score and sex education experience were also identified as protective factors against dynamic changes in the HIV-1 molecular network, with OR values of 0.12 (95% CI 0.05-0.32) and 0.26 (95% CI 0.11-0.64), respectively.

CONCLUSIONS

The HIV-1 molecular network analysis showed 144 sequences in static networks and 487 in dynamic networks. Multivariate analysis revealed that occupation as a student, floating population, Han ethnicity, and risky sexual behavior were independent risk factors for dynamic changes, while heterosexuality and homosexuality were protective compared to bisexuality. A higher National Eight-Item score and sex education experience were also protective factors. The identification of HIV dynamic molecular networks has provided valuable insights into the characteristics of individuals undergoing dynamic alterations. These findings contribute to a better understanding of HIV-1 transmission dynamics and could inform targeted prevention strategies.

Optimization of genetic distance threshold for inferring the CRF01_AE molecular network based on next-generation sequencing

Introduction

HIV molecular network based on genetic distance (GD) has been extensively utilized. However, the GD threshold for the non-B subtype differs from that of subtype B. This study aimed to optimize the GD threshold for inferring the CRF01_AE molecular network.

Methods

Next-generation sequencing data of partial CRF01_AE pol sequences were obtained for 59 samples from 12 transmission pairs enrolled from a high-risk cohort during 2009 and 2014. The paired GD was calculated using the Tamura-Nei 93 model to infer a GD threshold range for HIV molecular networks.

Results

2,019 CRF01_AE pol sequences and information on recent HIV infection (RHI) from newly diagnosed individuals in Shenyang from 2016 to 2019 were collected to construct molecular networks to assess the ability of the inferred GD thresholds to predict recent transmission events. When HIV transmission occurs within a span of 1-4 years, the mean paired GD between the sequences of the donor and recipient within the same transmission pair were as follow: 0.008, 0.011, 0.013, and 0.023 substitutions/site. Using these four GD thresholds, it was found that 98.9%, 96.0%, 88.2%, and 40.4% of all randomly paired GD values from 12 transmission pairs were correctly identified as originating from the same transmission pairs. In the real world, as the GD threshold increased from 0.001 to 0.02 substitutions/site, the proportion of RHI within the molecular network gradually increased from 16.6% to 92.3%. Meanwhile, the proportion of links with RHI gradually decreased from 87.0% to 48.2%. The two curves intersected at a GD of 0.008 substitutions/site.

Discussion

A suitable range of GD thresholds, 0.008-0.013 substitutions/site, was identified to infer the CRF01_AE molecular transmission network and identify HIV transmission events that occurred within the past three years. This finding provides valuable data for selecting an appropriate GD thresholds in constructing molecular networks for non-B subtypes.

Using the Molecular Transmission Networks to Analyze the Epidemic Characteristics of HIV-1 CRF08_BC in Kunming, Yunnan

HIV-1CRF08_BC is the most prevalent epidemic subtype among heterosexual (HET) and intravenous drug users (IDUs) in Kunming, Yunnan. Using the pol region of gene sequences derived from molecular epidemiological surveys, we developed a molecular transmission network for the purpose of analyzing its epidemiological characteristics, assessing its epidemiological trends, identifying its potential transmission relationships, and developing targeted interventions. HyPhy 2.2.4 was used to calculate pairwise genetic distances between sequences; GraphPad-Prism 8.0 was employed to determine the standard genetic distance; and Cytoscope 3.7.2 was applied to visualize the network. We used the network analysis tools to investigate network characteristics and the Molecular Complex Detection (MCODE) tool to observe the growth of the network. We utilized a logistic regression model to examine the factors influencing clustering and a zero-inflated Poisson model to investigate the factors influencing potential transmission links. At the standard genetic distance threshold of 0.008, 406 out of 858 study participants were clustered in 132 dissemination networks with a total network linkage of 868, and the number of links per sequence ranged from 1 to 19. The MCODE analysis identified three significant modular clusters in the networks, with network scores ranging from 4.9 to 7. In models of logistic regression, HET, middle-aged and elderly individuals, and residents of northern and southeastern Kunming were more likely to enter the transmission network. According to the zero-inflated Poisson model, age, transmission category, sampling year, marital status, and CD4+ T level had a significant effect on the size of links. The molecular clusters in Kunming's molecular transmission network are specific and aggregate to a certain extent. HIV-1 molecular network analysis provided information on local transmission characteristics, and these findings helped to determine the priority of transmission-reduction interventions.

HIV transmission dynamics and population-wide drug resistance in rural South Africa

Despite expanded antiretroviral therapy (ART) in South Africa, HIV-1 transmission persists. Integrase strand transfer inhibitors (INSTI) and long-acting injectables offer potential for superior viral suppression, but pre-existing drug resistance could threaten their effectiveness. In a community-based study in rural KwaZulu-Natal, prior to widespread INSTI usage, we enroled 18,025 individuals to characterise HIV-1 drug resistance and transmission networks to inform public health strategies. HIV testing and reflex viral load quantification were performed, with deep sequencing (20% variant threshold) used to detect resistance mutations. Phylogenetic and geospatial analyses characterised transmission clusters. One-third of participants were HIV-positive, with 21.7% having detectable viral loads; 62.1% of those with detectable viral loads were ART-naïve. Resistance to older reverse transcriptase (RT)-targeting drugs was found, but INSTI resistance remained low (<1%). Non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance, particularly to rilpivirine (RPV) even in ART-naïve individuals, was concerning. Twenty percent of sequenced individuals belonged to transmission clusters, with geographic analysis highlighting higher clustering in peripheral and rural areas. Our findings suggest promise for INSTI-based strategies in this setting but underscore the need for RPV resistance screening before implementing long-acting cabotegravir (CAB) + RPV. The significant clustering emphasises the importance of geographically targeted interventions to effectively curb HIV-1 transmission.

Molecular network analysis of 308 newly diagnosed HIV infection and 210 ART failure patients from rural counties in Sichuan

BACKGROUND

Few studies on molecular epidemiology have studied people with newly diagnosed HIV infection and ART Failure Patients at the same time in rural China. With more serious HIV epidemic than in other provinces in China, Sichuan is an area suitable for this study.

OBJECTIVE

To analyze the characteristics of HIV-1 molecular networks and factors related to network entry among newly diagnosed HIV infection and ART Failure Patients in three county-level cities (A, B, C) in Sichuan Province, to provide scientific basis for accurate prevention and control.

METHODS

Nested PCR amplification method was used to amplify HIV-1 pol gene region of 530 blood samples, Sequencer 4.9 was used to edit, clean and splice the gene sequence, Bioedit correction, Fastree 2.1.8 and Figtree 1.4.2 to construct evolutionary tree and determine genotype. HyPhy2.2.4 and Cytoscape 3.6.1 software were used to construct molecular network. Logistic regression analysis was applied.

RESULTS

523(98.68%) pol sequences were obtained, and a total of 518 valid sequences with basic information came into the final analyses. A total of 6 genotypes were detected, namely CRF01_AE (320,61.78%), CRF07_BC (149,28.76%), B (30,5.79%), CRF08_BC (11, 2.12%), CRF55_01B (6, 1.16%) and C (2, 0.39%). 186 of 518(35.91%) sequences entered the network at a genetic distance of 0.8%, forming 42 propagation clusters. "High-risk transmitters"(connected with two and more) accounted for 21.62%. Logistic regression showed that≥50 years old (OR = 2.474) were more risky than 18-49 years old, CRF07_BC sub-type (OR = 0.174) were less risky than CRF01_AE sub-type, B sub-type (OR = 6.698) is higher risky than CRF01_AE sub-type, and District B (OR = 0.077) less risky than that of A city.

CONCLUSION

The sources of HIV infection in rural Sichuan are diversified and complicated. The prevention and control of HIV infection in Sichuan Province should focus on strengthening the long-term dynamic detection of elderly population, B strain sub-type, and in City A.

Transmission network and phylogenetic analysis reveal older male-centered transmission of CRF01_AE and CRF07_BC in Guangxi, China

In China, the number of newly reported HIV infections in older people is increasing rapidly. However, clear information on the impact of older people on HIV transmission is limited. This study aims to reveal the local HIV transmission patterns, especially how older people affect virus transmission. Subtype analysis based on available pol sequences obtained from HIV patients revealed that CRF01_AE and CRF08_BC were predominant in patients aged <50 years, whereas CRF01_AE was predominant in older people aged ≥50 years (χ²= 29.299, P < 0.001). A total of 25 patients (5.2%, 25/484) were identified with recent HIV infection (RHI). Transmission network analysis found 267 genetically linked individuals forming 55 clusters (2-63 individuals), including 5 large transmission clusters and 12 transmission clusters containing RHI. Bayesian phylogenetic analysis suggested that transmission events in CRF01_AE and CRF07_BC were centred on older males, while transmission events in CRF08_BC were centred on younger males. Multivariable logistic regression analysis showed that older people were more likely to cluster within networks (AOR = 2.303, 95% CI: 1.012-5.241) and that RHI was a significant factor associated with high linkage (AOR = 3.468, 95% CI: 1.315-9.146). This study provides molecular evidence that older males play a central role in the local transmission of CRF01_AE and CRF07_BC in Guangxi. Given the current widespread of CRF01_AE and CRF07_BC in Guangxi, there is a need to recommend HIV screening as part of free national medical examinations for older people to improve early detection, timely treatment, and further reduce second-generation transmission.

Characteristics of the different HIV-1 risk populations based on the genetic transmission network of the newly diagnosed HIV cases in Jiangsu, Eastern China

Introduction

The HIV-1 prevalence has been steadily increasing in Jiangsu, China. HIV-1 genetic transmission network can be used to explore the transmission kinetics and precision intervention in high-risk populations. Thus, we generated an HIV-1 genetic transmission network, explored key risk populations based on different risk factors and found out the risk factors for HIV-1 prevention and control among the newly-diagnosed HIV-1 cases from 2017 to 2018.

Method

We amplified the HIV-1 pol sequences from the plasma samples of the newly-diagnosed HIV-1 cases from 2017 to 2018 and obtained the infection data from The National HIV/AIDS Surveillance System. HIV-Trace and Cytoscape Software were both used to construct the HIV-1 genetic network with a gene distance of <0.005. The R software was used to analyze the risk factors for inclusion into the network.

Results

We obtained 3362 sequences with the pol gene region, of which 3316 contained detailed individual information. CRF01_AE accounted for 42.3 % of the HIV-1 subtypes in the samples. The median CD4+T lymphocyte count was 329 cells/μL in 2017 and 313 cells/μL in 2018. At the gene distance threshold of 0.005, 481 sequences were incorporated into the HIV-1 gene network, constructing 202 clusters. Age over 60 years old, heterosexual transmission route, subtype (CRF105_0107, CRF55_01 B, and CRF67_01 B) and CD4+T lymphocyte count (>200) were the risk factors influencing inclusion into the HIV-1 gene network. Moreover, south Jiangsu cities had higher inclusion in the network. Thus, key risk populations in the clusters with different transmission routes, new emerging subtypes, drug resistance nodes, and individuals above 60 years of age in the network represented the critical risk populations that should be focused more on for intervention.

Conclusion

The HIV-1 genetic transmission network is adept at discovering undiagnosed HIV-infected cases and linking all diagnosed cases for determination of risk infections. Therefore we should pay more attention to these risk infections with further investigation and intervention, helping to achieve the goal of 95 % use combination prevention from the World Health Organization, and push to end AIDS epidemic.

HIV transmission dynamics and population-wide drug resistance in rural South Africa

Despite the scale-up of antiretroviral therapy (ART) in South Africa, HIV-1 incidence remains high. The anticipated use of potent integrase strand transfer inhibitors and long-acting injectables aims to enhance viral suppression at the population level and diminish transmission. Nevertheless, pre-existing drug resistance could impede the efficacy of long-acting injectable ART combinations, such as rilpivirine (an NNRTI) and cabotegravir (an INSTI). Consequently, a thorough understanding of transmission networks and geospatial distributions is vital for tailored interventions, including pre-exposure prophylaxis with long-acting injectables. However, empirical data on background resistance and transmission networks remain limited. In a community-based study in rural KwaZulu-Natal (2018-2019), prior to the widespread use of integrase inhibitor-based first-line ART, we performed HIV testing with reflex HIV-1 RNA viral load quantification on 18,025 participants. From this cohort, 6,096 (33.9%) tested positive for HIV via ELISA, with 1,323 (21.7%) exhibiting detectable viral loads (> 40 copies/mL). Of those with detectable viral loads, 62.1% were ART-naïve, and the majority of the treated were on an efavirenz + cytosine analogue + tenofovir regimen. Deep sequencing analysis, with a variant abundance threshold of 20%, revealed NRTI resistance mutations such as M184V in 2% of ART-naïve and 32% of treated individuals. Tenofovir resistance mutations K65R and K70E were found in 12% and 5% of ART-experienced individuals, respectively, and in less than 1% of ART-naïve individuals. Integrase inhibitor resistance mutations were notably infrequent (< 1%). Prevalence of pre-treatment drug resistance to NNRTIs was 10%, predominantly consisting of the K103N mutation. Among those with viraemic ART, NNRTI resistance was 50%, with rilpivirine-associated mutations observed in 9% of treated and 6% of untreated individuals. Cluster analysis revealed that 20% (205/1,050) of those sequenced were part of a cluster. We identified 171 groups with at least two linked participants; three quarters of clusters had only two individuals, and a quarter had 3-6 individuals. Integrating phylogenetic with geospatial analyses, we revealed a complex transmission network with significant clustering in specific regions, notably peripheral and rural areas. These findings derived from population scale genomic analyses are encouraging in terms of the limited resistance to DTG, but indicate that transitioning to long-acting cabotegravir + rilpivirine for transmission reduction should be accompanied by prior screening for rilpivirine resistance. Whole HIV-1 genome sequencing allowed identification of significant proportions of clusters with multiple individuals, and geospatial analyses suggesting decentralised networks can inform targeting public health interventions to effectively curb HIV-1 transmission.

Epidemiological characteristics of HIV transmission in southeastern China from 2015 to 2020 based on HIV molecular network

Objective

HIV/AIDS remains a global public health problem, and understanding the structure of social networks of people living with HIV/AIDS is of great importance to unravel HIV transmission, propose precision control and reduce new infections. This study aimed to investigate the epidemiological characteristics of HIV transmission in Fujian province, southeastern China from 2015 to 2020 based on HIV molecular network.

Methods

Newly diagnosed, treatment-naive HIV/AIDS patients were randomly sampled from Fujian province in 2015 and 2020. Plasma was sampled for in-house genotyping resistance test, and HIV molecular network was created using the HIV-TRACE tool. Factors affecting the inclusion of variables in the HIV molecular network were identified using univariate and multivariate logistic regression analyses.

Results

A total of 1,714 eligible cases were finally recruited, including 806 cases in 2015 and 908 cases in 2020. The dominant HIV subtypes were CRF01_AE (41.7%) and CRF07_BC (38.3%) in 2015 and CRF07_BC (53. 3%) and CRF01_AE (29.1%) in 2020, and the prevalence of HIV drug resistance was 4.2% in 2015 and 5.3% in 2020. Sequences of CRF07_BC formed the largest HIV-1 transmission cluster at a genetic distance threshold of both 1.5 and 0.5%. Univariate and multivariate logistic regression analyses showed that ages of under 20 years and over 60 years, CRF07_BC subtype, Han ethnicity, sampling in 2015, absence of HIV drug resistance, married with spouse, sampling from three cities of Jinjiang, Nanping and Quanzhou resulted in higher proportions of sequences included in the HIV transmission molecular network at a genetic distance threshold of 1.5% (p < 0.05).

Conclusion

Our findings unravel the HIV molecular transmission network of newly diagnosed HIV/AIDS patients in Fujian province, southeastern China, which facilitates the understanding of HIV transmission patterns in the province.

HIV-1 drug resistance and genetic transmission network among newly diagnosed people living with HIV/AIDS in Ningbo, China between 2018 and 2021

BACKGROUND

As the HIV epidemic continues to grow, transmitted drug resistance(TDR) and determining relationship of HIV transmission are major barriers to reduce the risk of HIV transmissions.This study aimed to examine the molecular epidemiology and TDR and evaluated the transmission pattern among newly diagnosed people living with HIV/AIDS(PLWHA) in Ningbo city, which could contribute to the development of targeted precision interventions.

METHODS

Consecutive cross-sectional surveys were conducted in Ningbo City between January 2018 and December 2021. The HIV-1 pol gene region was amplified and sequenced for drug resistance and genetic transmission network analysis. TDR was determined using the Stanford University HIV Drug Resistance Database. Genetic transmission network was visualized using Cytoscape with the genetic distance threshold of 0.013.

RESULTS

A total of 1006 sequences were sequenced successfully, of which 61 (6.1%) showed evidence of TDR. The most common mutations were K103N (2.3%), E138A/G/Q (1.7%) and V179D/E (1.2%). 12 HIV-1 genotypes were identified, with CRF07_BC being the major genotype (43.3%, 332/767), followed by CRF01_AE (33.7%, 339/1006). 444 (44.1%) pol sequences formed 856 links within 120 transmission clusters in the network. An increasing trend in clustering rate between 2018 and 2021(χ2 = 9.546, P = 0.023) was observed. The odds of older age (≥ 60 years:OR = 2.038, 95%CI = 1.072 ~ 3.872, compared to < 25 years), HIV-1 genotypes (CRF07_BC: OR = 2.147, 95%CI = 1.582 ~ 2.914; CRF55_01B:OR = 2.217, 95%CI = 1.201 ~ 4.091, compared to CRF01_AE) were significantly related to clustering. Compared with CRF01_AE, CRF07_BC were prone to form larger clusters. The largest cluster with CRF07_BC was increased from 15 cases in 2018 to 83 cases in 2021.

CONCLUSIONS

This study revealed distribution of HIV-1 genotypes, and genetic transmission network were diverse and complex in Ningbo city. The prevalence of TDR was moderate, and NVP and EFV were high-level NNRTI resistance. Individuals aged ≥ 60 years old were more easily detected in the networks and CRF07_BC were prone to form rapid growth and larger clusters. These date suggested that surveillance and comprehensive intervention should be designed for key rapid growth clusters to reduce the potential risk factors of HIV-1 transmission.

Using molecular network analysis to understand current HIV-1 transmission characteristics in an inland area of Yunnan, China

HIV-1 molecular surveillance provides a new approach to explore transmission risks and targeted interventions. From January to June 2021, 663 newly reported HIV-1 cases were recruited in Zhaotong City, Yunnan Province, China. The distribution characteristics of HIV-1 subtypes and HIV-1 molecular network were analysed. Of 542 successfully subtyped samples, 12 HIV-1 strains were identified. The main strains were CRF08_BC (47.0%, 255/542), CRF01_AE (17.0%, 92/542), CRF07_BC (17.0%, 92/542), URFs (8.7%, 47/542), and CRF85_BC (6.5%, 35/542). CRF08_BC was commonly detected among Zhaotong natives, illiterates, and non-farmers and was mostly detected in Zhaoyang County. CRF01_AE was frequently detected among married and homosexual individuals and mostly detected in Weixin and Zhenxiong counties. Among the 516 pol sequences, 187 (36.2%) were clustered. Zhaotong natives, individuals aged ≥60 years, and illiterate individuals were more likely to be found in the network. Assortativity analysis showed that individuals were more likely to be genetically associated when stratified by age, education level, occupation, and reporting area. The genetic diversity of HIV-1 reflects the complexity of local HIV epidemics. Molecular network analyses revealed the subpopulations to focus on and the characteristics of the risk networks. The results will help optimise local prevention and control strategies.

Combining molecular network analysis and field epidemiology to quantify local HIV transmission and highlight ongoing epidemics

OBJECTIVES

This study aimed to establish a collaborative approach to quantify local HIV transmission, which is an issue of great concern to public health.

METHODS

We linked HIV-1 pol gene sequences to demographic information and epidemiological investigations in Hangzhou (a central city in East China). We estimated local acquisition rates from a collaboration of molecular network analysis (with a distance-based approach) and epidemiological investigations.

RESULTS

Among 1064 newly diagnosed patients with HIV, 857 pol sequences were acquired and subsequently analyzed. Multiple subtypes were identified, with circulating recombinant form (CRF)07_BC (42.5%) and CRF01_AE (39.2%) predominating, followed by 13 other subtypes and 26 unique recombinant forms. By integrating the molecular network analysis and epidemiological investigations, we estimated that the proportion of local infection was 63.2%. The multivariable analyses revealed that individuals in clusters were more likely to be local residents, be aged 50 years or older, work as farmers, and have a higher first cluster of differentiation 4 count level (P <0.05). The proportions of local acquisitions over 70% were observed in local residents (79.9%, 242/303), individuals aged 50 years or older (73.6%, 181/246), and farmers (75.6%, 99/131).

CONCLUSION

The molecular network analysis can augment traditional HIV epidemic surveillance. This study establishes a paradigm for quantifying local HIV transmission for generalization in other areas.

Ending the HIV Epidemic: Identifying Barriers and Facilitators to Implement Molecular HIV Surveillance to Develop Real-Time Cluster Detection and Response Interventions for Local Communities

The rapid implementation of molecular HIV surveillance (MHS) has resulted in significant challenges for local health departments to develop real-time cluster detection and response (CDR) interventions for priority populations impacted by HIV. This study is among the first to explore professionals' strategies to implement MHS and develop CDR interventions in real-world public health settings. Methods: Semi-structured qualitative interviews were completed by 21 public health stakeholders in the United States' southern and midwestern regions throughout 2020-2022 to identify themes related to the implementation and development of MHS and CDR. Results for the thematic analysis revealed (1) strengths and limitations in utilizing HIV surveillance data for real-time CDR; (2) limitations of MHS data due to medical provider and staff concerns related to CDR; (3) divergent perspectives on the effectiveness of partner services; (4) optimism, but reluctance about the social network strategy; and (5) enhanced partnerships with community stakeholders to address MHS-related concerns. Conclusions: Enhancing MHS and CDR efforts requires a centralized system for staff to access public health data from multiple databases to develop CDR interventions; designating staff dedicated to CDR interventions; and establishing equitable meaningful partnerships with local community stakeholders to address MHS concerns and develop culturally informed CDR interventions.

Molecular Transmission Network of Newly Reported HIV Infections in Pengzhou, Sichuan Province: A Study Based on Genomics and Spatial Epidemiology

OBJECTIVE

The objective of this study was to understand the molecular transmission characteristics of newly reported HIV infections in the city of Pengzhou, Sichuan Province, to analyze the risk factors of transmission network and spatial clustering and the transmission characteristics, and to provide a scientific basis for precision prevention and intervention.

METHODS

Anticoagulated whole blood was collected from newly reported HIV infections in Pengzhou from March 2019 to August 2021. After the plasma was isolated, the HIV-1 pol gene was amplified and sequenced by reverse transcriptase polymerase chain reaction (PCR). The obtained gene sequences were used to construct a maximum likelihood phylogenetic tree for the analysis of virus subtypes, and a molecular transmission network was constructed using the genetic distance method to evaluate the transmission pattern of people living with HIV/AIDS in Pengzhou. A logistic regression model was used to find out the potential risk factors for entering the molecular transmission network with the number of nodes ≥ 2. Spatial analysis is used to show the geographical pattern of the proportion of newly reported HIV infections entering the molecular transmission network, and a flow map is used to show the intensity of transmission within and between townships.

RESULTS

A total of 463 newly reported HIV-infection sequences were obtained in this study, including 237 cases (51.19%) of CRF01_ AE, 159 cases (34.34%) of CRF07_BC, 45 cases (9.72%) of B, 15 cases (3.24%) of CRF08_BC and 7 cases (1.5%) of others. The number of clusters was the highest when the gene distance was 0.009, with a total of 246 sequences entering the network, forming 54 clusters, and the network entry rate was 55.36%. There were 170 sequences with more than two nodes in the network sequence. The logistic regression showed that compared with age < 50 years old, age ≥ 50 years old has a higher risk of transmission (OR = 3.43, 95% CI = 2.06-5.71); compared with farmers, the risk of transmission within industry is lower (OR = 0.046, 95% CI = 0.25-0.87); and compared with CRF07_BC, CRF01_AE (OR = 6.09, 95% CI = 3.60-10.30) and B (OR = 20.31, 95% CI = 8.94-46.13) had a higher risk of transmission. Men aged ≥ 50 years are mainly clustered with women between 50 and 70 years of age. In addition to being clustered with gay men, there are nine (50%) and three (16.7%) chains of transmission between gay men and heterosexual men and women, respectively. In the geographical space, there is no hot spot clustering of the molecular propagation network. The subtype B was mainly distributed in the town of Tianpeng and formed transmission networks in eastern Pengzhou;0020CRF01_AE is mainly distributed in the town of Lichun and formed transmission networks in the west and north of Pengzhou.

CONCLUSION

This study reveals the characteristics and influencing factors of molecular network transmission in the region, as well as the spatial transmission characteristics of newly reported HIV infections in recent years, and reveals the geographical differences in HIV-1 transmission. The results provide a scientific basis for the development of local AIDS-specific intervention measures.

Community Insights in Phylogenetic HIV Research: The CIPHR Project Protocol

Inferring HIV transmission networks from HIV sequences is gaining popularity in the field of HIV molecular epidemiology. However, HIV sequences are often analyzed at distance from those affected by HIV epidemics, namely without the involvement of communities most affected by HIV. These remote analyses often mean that knowledge is generated in absence of lived experiences and socio-economic realities that could inform the ethical application of network-derived information in 'real world' programmes. Procedures to engage communities are noticeably absent from the HIV molecular epidemiology literature. Here we present our team's protocol for engaging community activists living in Nairobi, Kenya in a knowledge exchange process - The CIPHR Project (Community Insights in Phylogenetic HIV Research). Drawing upon a community-based participatory approach, our team will (1) explore the possibilities and limitations of HIV molecular epidemiology for key population programmes, (2) pilot a community-based HIV molecular study, and (3) co-develop policy guidelines on conducting ethically safe HIV molecular epidemiology. Critical dialogue with activist communities will offer insight into the potential uses and abuses of using such information to sharpen HIV prevention programmes. The outcome of this process holds importance to the development of policy frameworks that will guide the next generation of the global response.

Non-disclosed men who have sex with men within local MSM HIV-1 genetic transmission networks in Guangyuan, China

Background

Most men who have sex with men (MSM), especially those with HIV infection, do not disclose their same-sex behaviors in China due to Chinese family values and fear of stigmatization, rejection, or prejudice. However, disclosure of same-sex behaviors to healthcare providers (HCPs) can be beneficial for reducing viral transmission and promoting their physical and mental health. In this study, by combining phylogenetic analysis with traditional epidemiological approaches, we tried to identify the MSM who do not disclose to HCPs in transmission networks and explored the factors related to the non-disclosed behaviors.

Method

Phylogenetic analysis was conducted using HIV pol sequences obtained from the drug-resistant surveillance program, which was collected as part of routine clinical care since 2012. Sequences were linked to the demographic data collected in the Chinese HIV/AIDS Comprehensive Response Information Management System (CRIMS). First, male patients in whom genetic sequences were within the molecular transmission clusters involving self-reported MSM were identified as potential MSM (pMSM). Then, a cross-sectional survey was conducted to supplement behavioral information and attitudes toward MSM.

Results

Our sample consisted of 190 pMSM patients. In total, 43.16% of the patients were likely to conceal same-sex behaviors during the first-self-report, and 14.73% of patients might continue to conceal a history of same-sex behaviors even after receiving medical care. The pMSM who concealed their same-sex behaviors were reluctant to accept medical services such as Voluntary Counseling and Testing (VCT) and had a lower likelihood of condom use. In addition, the related factors for non-disclosed behavior were associated with current address, income before diagnosis, and attitudes toward MSM.

Conclusion

Non-disclosure of same-sex behaviors to HCPs may be a major obstacle for certain medical services for MSM who exhibit risky sexual behaviors. The pMSM from developing areas, with high monthly income, and with neutral or un-supportive attitudes toward MSM may represent non-disclosure of their same-sex behaviors. Thus, policies facilitating MSM to disclose their same-sex behaviors are recommended, such as legislations protecting homosexual rights on employment, education, marriage, and so on.

Non-student young men put students at high risk of human immunodeficiency virus acquisition in Guangxi, China: a phylogenetic analysis of surveillance data

Background

We sought to identify students and their sexual partners in a molecular transmission network.

Methods

We obtained 5996 HIV protease and reverse transcriptase gene sequences in Guangxi (165 from students and 5831 from the general populations) and the relevant demographic data. We constructed a molecular transmission network and introduced a permutation test to assess the robust genetic linkages. We calculated the centrality measures to describe the transmission patterns in clusters.

Results

At the network level, 68 (41.2%) students fell within the network across 43 (8.1%) clusters. Of 141 genetic linkages between students and their partners, only 25 (17.7%) occurred within students. Students were more likely than random permutations to link to other students (odds ratio [OR], 7.2; P < .001), private company employees aged 16–24 years (OR, 3.3; P = .01), private company or government employees aged 25–49 years (OR, 1.7; P = .03), and freelancers or unemployed individuals aged 16–24 years (OR, 5.0; P < .001). At the cluster level, the median age of nonstudents directly linked to students (interquartile range) was 25 (22–30) years, and 80.3% of them had a high school or higher education background. Compared with students, they showed a significantly higher median degree (4.0 vs 2.0; P < .001) but an equivalent median Eigenvector Centrality (0.83 vs 0.81; P = .60).

Conclusions

The tendency of genetic linkage between students and nonstudent young men and their important position in the HIV transmission network emphasizes the urgent need for 2-pronged public health interventions based on both school and society.

Integrated molecular and affiliation network analysis: Core-periphery social clustering is associated with HIV transmission patterns

This study investigates the two-mode core-periphery structures of venue affiliation networks of younger Black men who have sex with men (YBMSM). We examined the association between these structures and HIV phylogenetic clusters, defined as members who share highly similar HIV strains that are regarded as a proxy for sexual affiliation networks. Using data from 114 YBMSM who are living with HIV in two large U.S. cities, we found that HIV phylogenetic clustering patterns were associated with social clustering patterns whose members share affiliation with core venues that overlap with those of YBMSM. Distinct HIV transmission patterns were found in each city, a finding that can help to inform tailored venue-based and network intervention strategies.

Epidemiology and Molecular Transmission Characteristics of HIV in the Capital City of Anhui Province in China

Hefei, Anhui province, is one of the cities in the Yangtze River Delta, where many people migrate to Jiangsu, Zhejiang and Shanghai. High migration also contributes to the HIV epidemic. This study explored the HIV prevalence in Hefei to provide a reference for other provinces and assist in the prevention and control of HIV in China. A total of 816 newly reported people with HIV in Hefei from 2017 to 2020 were recruited as subjects. HIV subtypes were identified by a phylogenetic tree. The most prevalent subtypes were CRF07_BC (41.4%), CRF01_AE (38.1%) and CRF55_01B (6.3%). Molecular networks were inferred using HIV-TRACE. The largest and most active transmission cluster was CRF55_01B in Hefei’s network. A Chinese national database (50,798 sequences) was also subjected to molecular network analysis to study the relationship between patients in Hefei and other provinces. CRF55_01B and CRF07_BC-N had higher clustered and interprovincial transmission rates in the national molecular network. People with HIV in Hefei mainly transmitted the disease within the province. Finally, we displayed the epidemic trend of HIV in Hefei in recent years with the dynamic change of effective reproductive number (Re). The weighted overall Re increased rapidly from 2012 to 2015, with a peak value of 3.20 (95% BCI, 2.18–3.85). After 2015, Re began to decline and remained stable at around 1.80. In addition, the Re of CRF55_01B was calculated to be between 2.0 and 4.0 in 2018 and 2019. More attention needs to be paid to the rapid spread of CRF55_01B and CRF07_BC-N strains among people with HIV and the high Re in Hefei. These data provide necessary support to guide the targeted prevention and control of HIV.

HIV-1 CRF07_BC transmission dynamics in China: two decades of national molecular surveillance

By analyzing an unprecedentedly large, longitudinal HIV-1 CRF07_BC sequence dataset collected from China in the past two decades, we sought to build CRF07_BC lengthwise transmission networks, and understand its transmission dynamics. We divided CRF07_BC into two clusters based on phylogenetic analysis and an estimation of the pairwise genetic distance at 0.7%. Of 6213 sequences, 3607 (58.1%) linked to ≥1 other sequence. CRF07_BC was divided into two clusters: 07BC_O and 07BC_N. The 07BC_O is the original CRF07_BC, circulating in people who inject drugs (PWID) and heterosexuals, predominantly in southwestern and northwestern provinces of China. The 07BC_N is a new cluster, identified mostly in men having sex with men (MSM) in the northern provinces of China. Bayesian analysis indicates that CRF07_BC has experienced two phases of exponential growth, which was first driven by 07BC_O then 07BC_N. Compared to 07BC_O, the proportion of the parameter of population transmission risk (TR) of 07BC_N has risen constantly. The power-law function analyses reveal that 07BC_N has increased over years with higher degree. In 07BC_N, only 13.16% of MSM were linked to other risk groups, but these links represent 41.45%, 54.25%, and 55.07% of links among heterosexual females, heterosexual males, and male PWID respectively. This study indicates that CRF07_BC has evolved into two clusters in China, and their distributions are distinct across risk groups and geographical regions. 07BC_N shows a greater risk of transmission, and has gradually replaced 07BC_O. Furthermore, the results show that strengthening the MSM interventions could lower the rapidity of 07BC_N transmission in all risk groups.

Analysis of the Driving Factors of Active and Rapid Growth Clusters Among CRF07_BC-Infected Patients in a Developed Area in Eastern China

Background

The purpose of this study was to research the molecular transmission and genetic evolutionary characteristics among CRF07_BC-infected patients in a developed area in Eastern China.

Methods

Plasma samples from newly diagnosed HIV-1-positive patients from 2015–2018 and basic demographic and epidemiological information were obtained. Pol sequences from CRF07_BC-infected patients were selected for phylogenetic, molecular transmission network, and Bayesian evolutionary analyses.

Results

Pol sequences were successfully obtained from 258 samples of CRF07_BC. Phylogenetic analysis revealed 2 distinct lineages: lineage 1 (66.3%, 171/258), primarily from men who have sex with men (MSM) and some heterosexual individuals, and lineage 2 (33.7%, 87/258), primarily from heterosexual individuals. Under an optimal genetic distance of 0.01 substitutions/site, 163 individuals (63.2%, 163/258) formed 23 groups comprising 6 clusters and 17 dyads in the networks. A distinctly large and rapidly growing cluster (C1) containing 105 individuals was identified, in which MSM with ≥4 links had quite a high transmission risk (low educational background, active sexual behavior, low sexual protection awareness, etc.). According to Bayesian analyses, most C1 clades formed from 2005 to 2009, most of which were closely geographically related to CRF07_BC epidemic strains from Anhui province.

Conclusions

Here, we elucidated the local transmission characteristics and epidemic pattern of HIV-1 CRF07_BC, revealing that MSM (especially with ≥4 links) may be a significant driver in the formation of active and rapid growth networks in regional CRF07_BC epidemics. Thus, unique region– and risk group–specific transmission network analysis based on a molecular approach can provide critical and insightful information for more effective intervention strategies to limit future HIV-1 transmission.

Advances in research of HIV transmission networks

Transmission network analysis is a crucial evaluation tool aiming to explore the characteristics of the human immunodeficiency virus epidemic, develop evidence-based prevention strategies, and contribute to various areas of human immunodeficiency virus/acquired immunodeficiency syndrome prevention and control. Over recent decades, transmission networks have made tremendous strides in terms of modes, methods, applications, and various other aspects. Transmission network methods, including social, sexual, and molecular transmission networks, have played a pivotal role. Each transmission network research method has its advantages, as well as its limitations. In this study, we established a systematic review of these aforementioned transmission networks with respect to their definitions, applications, limitations, recent progress, and synthetic applications.

Dynamics of HIV-1 Molecular Networks Reveal Effective Control of Large Transmission Clusters in an Area Affected by an Epidemic of Multiple HIV Subtypes

This study reconstructed molecular networks of human immunodeficiency virus (HIV) transmission history in an area affected by an epidemic of multiple HIV-1 subtypes and assessed the efficacy of strengthened early antiretroviral therapy (ART) and regular interventions in preventing HIV spread. We collected demographic and clinical data of 2221 treatment-naïve HIV-1–infected patients in a long-term cohort in Shenyang, Northeast China, between 2008 and 2016. HIV pol gene sequencing was performed and molecular networks of CRF01_AE, CRF07_BC, and subtype B were inferred using HIV-TRACE with separate optimized genetic distance threshold. We identified 168 clusters containing ≥ 2 cases among CRF01_AE-, CRF07_BC-, and subtype B-infected cases, including 13 large clusters (≥ 10 cases). Individuals in large clusters were characterized by younger age, homosexual behavior, more recent infection, higher CD4 counts, and delayed/no ART (P < 0.001). The dynamics of large clusters were estimated by proportional detection rate (PDR), cluster growth predictor, and effective reproductive number (R_e). Most large clusters showed decreased or stable during the study period, indicating that expansion was slowing. The proportion of newly diagnosed cases in large clusters declined from 30 to 8% between 2008 and 2016, coinciding with an increase in early ART within 6 months after diagnosis from 24 to 79%, supporting the effectiveness of strengthened early ART and continuous regular interventions. In conclusion, molecular network analyses can thus be useful for evaluating the efficacy of interventions in epidemics with a complex HIV profile.