Molecular epidemiology of a primarily MSM acute HIV-1 cohort in Bangkok, Thailand and connections within networks of transmission in Asia

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.

Characteristics of genotype, drug resistance, and molecular transmission network among newly diagnosed HIV-1 infections in Shenzhen, China

The HIV-1 pandemic has persisted for four decades, and poses a major challenge to global public health. Shenzhen, a city with large number of migrant populations in China, is suffering HIV-1 epidemic. It is necessary to continuously conduct the molecular surveillance among newly diagnosed HIV-1 patients in these migrant population. In this study, plasma samples of newly diagnosed and ART-naive HIV-1 infections were collected from Shenzhen city in China. The partial genes of HIV-1 gag and pol were amplified and sequenced for the analysis of genotype, drug resistance, and molecular transmission network. Ninety-one sequences of pol gene were obtained from newly diagnosed HIV-1 infections in Shenzhen, and seven HIV-1 subtypes were revealed in this investigation. Among them, the circulating recombinant form (CRF) 07_BC was the mostly frequent subtype (53.8%, 49/91), followed by CRF01_AE (20.9%, 19/91), CRF55_01B (9.9%, 9/91), unique recombinant forms (URFs) (8.8%, 8/91), B (3.3%, 3/91), CRF59_01B (2.2%, 2/91), and CRF08_BC (1.1%, 1/91). The overall prevalence of pretreatment drug resistance (PDR) was 23.1% (21/91), and 52.38% (11/21) of the PDR was specific for the nonnucleotide reverse transcriptase inhibitors (NNRTIs). Furthermore, a total of 3091 pol gene sequences were used to generate 19 molecular transmission clusters, and then one growing cluster, a new cluster, and a cluster with growth reactivation were identified. The result revealed that more sexual partner, CRF_07BC subtype, and seven amino acid deletions in gag p6 region might be the influencing factors associated with the high risk of transmission behavior. Compared with CRF01_AE subtype, CRF07_BC subtype strains were more likely to form clusters in molecular transmission network. This suggests that long-term surveillance of the HIV-1 molecular transmission should be a critical measure to achieve a precise intervention for controlling the spread of HIV-1 in China.

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.

Evolution of Antibody Responses in HIV-1 CRF01_AE Acute Infection: Founder Envelope V1V2 Impacts the Timing and Magnitude of Autologous Neutralizing Antibodies

Understanding the dynamics of early immune responses to HIV-1 infection, including the evolution of initial neutralizing and antibody-dependent cellular cytotoxicity (ADCC)-mediating antibodies, will inform HIV vaccine design. In this study, we assess the development of autologous neutralizing antibodies (ANAbs) against founder envelopes (Envs) from 18 participants with HIV-1 CRF01_AE acute infection. The timing of ANAb development directly associated with the magnitude of the longitudinal ANAb response. Participants that developed ANAbs within 6 months of infection had significantly higher ANAb responses at 1 year (50% inhibitory concentration [IC50] geometric mean titer [GMT] = 2,010 versus 184; P = 0.001) and 2 years (GMT = 3,479 versus 340; P = 0.015), compared to participants that developed ANAb responses after 6 months. Participants with later development of ANAb tended to develop an earlier, potent heterologous tier 1 (92TH023) neutralizing antibody (NAb) response (P = 0.049). CRF01_AE founder Env V1V2 loop lengths correlated indirectly with the timing (P = 0.002, r = -0.675) and directly with magnitude (P = 0.005, r = 0.635) of ANAb responses; Envs with longer V1V2 loop lengths elicited earlier and more potent ANAb responses. While ANAb responses did not associate with viral load, the viral load set point correlated directly with neutralization of the heterologous 92TH023 strain (P = 0.007, r = 0.638). In contrast, a striking inverse correlation was observed between viral load set point and peak ADCC against heterologous 92TH023 Env strain (P = 0.0005, r = -0.738). These data indicate that specific antibody functions can be differentially related to viral load set point and may affect HIV-1 pathogenesis. Exploiting Env properties, such as V1V2 length, could facilitate development of subtype-specific vaccines that elicit more effective immune responses and improved protection. IMPORTANCE Development of an effective HIV-1 vaccine will be facilitated by better understanding the dynamics between the founder virus and the early humoral responses. Variations between subtypes may influence the evolution of immune responses and should be considered as we strive to understand these dynamics. In this study, autologous founder envelope neutralization and heterologous functional humoral responses were evaluated after acute infection by HIV-1 CRF01_AE, a subtype that has not been thoroughly characterized. The evolution of these humoral responses was assessed in relation to envelope characteristics, magnitude of elicited immune responses, and viral load. Understanding immune parameters in natural infection will improve our understanding of protective responses and aid in the development of immunogens that elicit protective functional antibodies. Advancing our knowledge of correlates of positive clinical outcomes should lead to the design of more efficacious vaccines.

Tracking coreceptor switch of the transmitted/founder HIV-1 identifies co-evolution of HIV-1 antigenicity, coreceptor usage and CD4 subset targeting

The CCR5 (R5) to CXCR4 (X4) coreceptor switch in natural HIV-1 infection is associated with faster progression to AIDS, but the underlying mechanisms remain unclear. The difficulty in capturing the earliest moment of coreceptor switch in vivo limits our understanding of this phenomenon. Here, by tracking the evolution of the transmitted/founder (T/F) HIV-1 in a prospective cohort of individuals at risk for HIV-1 infection identified very early in acute infection, we investigated this process with high resolution. The earliest X4 variants evolved from the R5 tropic T/F strains. Strong X4 usage can be conferred by a single mutation. The mutations responsible for coreceptor switch can confer escape to neutralization and drive X4 variants to replicate mainly in the central memory and naïve CD4+ T cells. We propose a novel concept to explain the co-evolution of virus antigenicity and entry tropism termed "escape by shifting". This concept posits that for viruses with receptor or coreceptor flexibility, entry tropism alteration represents a mechanism of immune evasion in vivo .

Using Molecular Transmission Networks to Reveal the Epidemic of Pretreatment HIV-1 Drug Resistance in Guangxi, China

Introduction

Pretreatment drug resistance (PDR) is becoming an obstacle to the success of ART. This study investigated the prevalence of PDR and the transmission clusters (TCs) of drug resistance mutations (DRMs) in two cities where drug abuse used to be high to describe the local HIV-1 transmission dynamics.

Methods

Plasma samples were obtained from 1,027 ART-naïve patients in Guangxi. Viral subtypes and DRMs were identified. Transmission network and related factors were also determined.

Results

A total of 1,025 eligible sequences were obtained from Qinzhou (65.8%) and Baise (34.2%) cities. The predominant HIV-1 genotype was CRF08_BC (45.0%), followed by CRF01_AE (40.9%). The overall prevalence of PDR was 8.3%, and resistance to NNRTI was the most common. Putative links with at least one other sequence were found in 543/1,025 (53.0%) sequences, forming 111 clusters (2–143 individuals). The most prevalent shared DRMs included V106I (45.35%), V179D (15.1%), and V179E (15.1%). Clusters related to shared DRMs were more frequent and larger in CRF08_BC. The prevalence of shared DRMs increased with time, while the proportion of PDR gradually decreased. Age > 50 years was associated with clustering. Subtype CRF08_BC was more likely to have DRMs, PDR propagation, and DRM sharing.

Conclusion

PDR prevalence is moderate in this region. The association between PDR and subtype CRF08_BC suggested that DRMs spreading from injection drug users (IDUs) to heterosexuals (HETs) might be the major source of PDR in this region. Our findings highlight the significance of continuous surveillance of PDR.

The Emerging of CRF01_AE: A Clinical Story and Future HIV/AIDS Situation in Thailand

The HIV epidemic in Thailand in the 1980's compromised the country's socio-economic development. The epidemic first became evident in the community of men with male sexual partners (MSM), and subsequently spread to intravenous drug users (IVDU), female commercial sex workers (CSW) and their male clients, and, ultimately, to their partners and children. The HIV epidemic has devastated the country's working-age population. The extensive negative impact and social stigma associated with the disease do not only have an impact on the victims of HIV but also on their descendants and relatives. An epicenter of the HIV epidemic has been in the northern provinces of Thailand. An HIV-1 subtype CRF01_AE, a complex chimeric virus composed of both A and E subtypes, is prevalent in Northern Thailand. The virus has quickly become a predominant viral strain circulating in Thailand, other neighboring Southeast Asian countries, and China as well as some other countries throughout the world. The epidemiology, evolution, and biology of CRF01_AE offer a unique model for further scientific investigations which would advance the knowledge of and curative strategies against HIV. In addition, Thailand has developed suitable national guidelines on HIV/AIDS treatment and prevention in order to control the epidemic. Effective antiretroviral drugs are, therefore, able to be made available to those who live with HIV. The national surveillance system has also been effective. The great efforts and resources which Thailand has dedicated to the fight against the epidemic have eventually paid off. In 2010, a plan was proposed to eliminate mother-to-child HIV transmission and Thailand has become the first country to be effective in this objective. Thailand therefore has become recognized as being the global leader in HIV prevention and treatment. The experience which Thailand has gained from the past and the current research and management strategies of the HIV epidemic has prepared the country for emerging strains of HIV-1 in the future.

HIV-1 subtype diversity and transmission strain source among men who have sex with men in Guangxi, China

With the rapid increase in HIV prevalence of men who have sex with men (MSM) in recent years and common human migration and travelling across different provinces in China, MSM are now finding it easier to meet each other, which might contribute to local HIV epidemics as well as fueling cross-province transmission. We performed a cross-sectional survey in 2018-2019 to investigate the current HIV subtype diversity and inferred HIV strain transmission origin among MSM in Guangxi province, China based on a phylogenetic analysis. Based on 238 samples, we found that the HIV-1 subtype diversity was more complicated than before, except for three major HIV subtypes/circulating recombinant forms (CRFs): CRF07_BC, CRF01_AE, CRF55_01B, five other subtypes/CRFs (CRF59_01B, B, CRF08_BC, CRF67_01B, CRF68_01B) and five unique recombinant forms (URFs) were detected. In total, 76.8% (169/220) of samples were infected with HIV from local circulating strains, while others originated from other provinces, predominantly Guangdong and Shanghai. The high diversity of HIV recombinants and complicated HIV transmission sources in Guangxi MSM indicates that there has been an active sexual network between HIV positive MSM both within and outside Guangxi without any effective prevention. Inter-province collaboration must be enforced to provide tailored HIV prevention and control services to MSM in China.

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.

Dynamics of Human Immunodeficiency Virus-1 Genetic Diversification During Acute Infection

We analyzed human immunodeficiency virus envelope diversity in 98 acute infections. The within-host genetic diversity, divergence from transmitted/founder (T/F) strain, and the observed frequency of multiple T/F infections increased with Fiebig stage. These data identify rapid viral dynamics during acute infection with implications for clinical trials conducted in this setting.

A Geographic Hotspot and Emerging Transmission Cluster of the HIV-1 Epidemic Among Older Adults in a Rural Area of Eastern China

Despite the implementation of health education and free condom distribution for decades, the HIV/AIDS epidemic among older adults in China shows no sign of declining. This study aims to identify HIV transmission patterns and pathways in a rural county area and provide insight for developing effective HIV prevention strategies among older adults. Epidemiological field surveys combined with phylogenetic analysis were used to identify potential HIV transmission linkage in one rural county with a rapidly increasing HIV epidemic among older adults. A total of 160 HIV-positive individuals and their HIV-positive sexual partners diagnosed between 2015 and 2018 were recruited. Among them, 69.4% (n = 111) were diagnosed at age 50 or older, 68.1% (n = 109) were men, and 55.0% (n = 88) were married. Ninety-six participants self-reported as having been infected with HIV through commercial heterosexual behavior and 24 were infected from their HIV-positive spouse. CRF08_BC (67.0%, n = 73) and CRF07_AE (19.3%, n = 21) were the prevalent HIV stains and formed 15 clusters at distance cutoff of 0.0025. Participants diagnosed with HIV at age ≥50, those who were divorced or bereft their spouse, or infected with HIV locally were more prevalent within the 15 molecular clusters. We identified specific venues as the primary hotspot and found that commercial heterosexual contact between older men and commercial sex workers contributed to the local HIV epidemic. In rural areas with an increasing HIV epidemic, detecting HIV transmission clusters through epidemiological and phylogenetic analysis as well as designing localized HIV prevention strategies should be prioritized.

Incorporation of information diffusion model for enhancing analyses in HIV molecular surveillance

Molecular surveillance of infections is essential in monitoring their transmission in the population. In this study, newly diagnosed HIV patients' phylogenetic, clinical and behavioural data were integrated, and an information diffusion model was incorporated in analysing transmission dynamics. A genetic network was constructed from HIV sequences, from which transmission cascades were extracted. From the transmission cascades, CRF01_AE had higher values of information diffusion metrics, including scale, speed and range, than that of B, signifying the distinct transmission patterns of two circulating subtypes in Hong Kong. Patients connected in the network, were more likely male, younger, of main circulating subtypes, to have acquired HIV infection locally, and a higher CD4 level at diagnosis. Genetic connections varied among men who have sex with men (MSM) who used different channels of sex networking and varied in their engagement in risk behaviours. MSM using recreational drugs for sex held positions of greater importance within the network. Significant differences in network metrics were observed among MSM as differentiated by their mobile apps usage patterns, evidencing the impact of social network on transmission networks. The applied model in the presence of consistently collected longitudinal data could enhance HIV molecular epidemiologic surveillance for informing future intervention planning.