Human Immunodeficiency Virus Type 1 Phylodynamics to Detect and Characterize Active Transmission Clusters in North Carolina

Molecular Epidemiology of HIV-1 Subtype B Infection across Florida Reveals Few Large Superclusters with Metropolitan Origin

Florida is considered an epicenter of HIV in the United States. The U.S. federal plan for Ending the HIV Epidemic (EHE) within 10 years prioritizes seven of Florida's 67 counties for intervention. We applied molecular epidemiology methods to characterize the HIV infection networks in the state and infer whether the results support the EHE. HIV sequences (N = 34,446) and associated clinical/demographic metadata of diagnosed people with HIV (PWH), during 2007 to 2017, were retrieved from the Florida Department of Health. HIV genetic networks were investigated using MicrobeTrace. Associates of clustering were identified through boosted logistic regression. Assortative trait mixing was also assessed. Bayesian phylogeographic methods were applied to evaluate evidence of imported HIV-1 lineages and illustrate spatiotemporal flows within Florida. We identified nine large clusters spanning all seven EHE counties but little evidence of external introductions, suggesting-in the absence of undersampling-an epidemic that evolved independently from the rest of the country or other external influences. Clusters were highly assortative by geography. Most of the sampled infections (82%) did not cluster with others in the state using standard molecular surveillance methods despite satisfactory sequence sampling in the state. The odds of being unclustered were higher among PWH in rural regions, and depending on demographics. A significant number of unclustered sequences were observed in counties omitted from EHE. The large number of missing sequence links may impact timely detection of emerging transmission clusters and ultimately hinder the success of EHE in Florida. Molecular epidemiology may help better understand infection dynamics at the population level and underlying disparities in disease transmission among subpopulations; however, there is also a continuous need to conduct ethical discussions to avoid possible harm of advanced methodologies to vulnerable groups, especially in the context of HIV stigmatization. IMPORTANCE The large number of missing phylogenetic linkages in rural Florida counties and among women and Black persons with HIV may impact timely detection of ongoing and emerging transmission clusters and ultimately hinder the success of epidemic elimination goals in Florida.

Optimized phylogenetic clustering of HIV-1 sequence data for public health applications

Clusters of genetically similar infections suggest rapid transmission and may indicate priorities for public health action or reveal underlying epidemiological processes. However, clusters often require user-defined thresholds and are sensitive to non-epidemiological factors, such as non-random sampling. Consequently the ideal threshold for public health applications varies substantially across settings. Here, we show a method which selects optimal thresholds for phylogenetic (subset tree) clustering based on population. We evaluated this method on HIV-1 pol datasets (n = 14, 221 sequences) from four sites in USA (Tennessee, Washington), Canada (Northern Alberta) and China (Beijing). Clusters were defined by tips descending from an ancestral node (with a minimum bootstrap support of 95%) through a series of branches, each with a length below a given threshold. Next, we used pplacer to graft new cases to the fixed tree by maximum likelihood. We evaluated the effect of varying branch-length thresholds on cluster growth as a count outcome by fitting two Poisson regression models: a null model that predicts growth from cluster size, and an alternative model that includes mean collection date as an additional covariate. The alternative model was favoured by AIC across most thresholds, with optimal (greatest difference in AIC) thresholds ranging 0.007-0.013 across sites. The range of optimal thresholds was more variable when re-sampling 80% of the data by location (IQR 0.008 - 0.016, n = 100 replicates). Our results use prospective phylogenetic cluster growth and suggest that there is more variation in effective thresholds for public health than those typically used in clustering studies.

Undiagnosed HIV Infections May Drive HIV Transmission in the Era of "Treat All": A Deep-Sampling Molecular Network Study in Northeast China during 2016 to 2019

Universal antiretroviral therapy (ART, “treat all”) was recommended by the World Health Organization in 2015; however, HIV-1 transmission is still ongoing. This study characterizes the drivers of HIV transmission in the “treat all” era. Demographic and clinical information and HIV pol gene were collected from all newly diagnosed cases in Shenyang, the largest city in Northeast China, during 2016 to 2019. Molecular networks were constructed based on genetic distance and logistic regression analysis was used to assess potential transmission source characteristics. The cumulative ART coverage in Shenyang increased significantly from 77.0% (485/630) in 2016 to 93.0% (2598/2794) in 2019 (p < 0.001). Molecular networks showed that recent HIV infections linked to untreated individuals decreased from 61.6% in 2017 to 28.9% in 2019, while linking to individuals with viral suppression (VS) increased from 9.0% to 49.0% during the same time frame (p < 0.001). Undiagnosed people living with HIV (PLWH) hidden behind the links between index cases and individuals with VS were likely to be male, younger than 25 years of age, with Manchu nationality (p < 0.05). HIV transmission has declined significantly in the era of “treat all”. Undiagnosed PLWH may drive HIV transmission and should be the target for early detection and intervention.

Exploring the Evolutionary History and Phylodynamics of Human Immunodeficiency Virus Type 1 Outbreak From Unnao, India Using Phylogenetic Approach

Certain rural and semiurban settings in the Unnao district, Uttar Pradesh, India observed an unprecedented increase in the detection of HIV cases during July 2017. Subsequent investigations through health camps and a follow-up case-control study attributed the outbreak to the unsafe injection exposures during treatment. In this study, we have undertaken a secondary analysis to understand the phylogenetic aspects of the outbreak-associated HIV-1 sequences along with the origin and phylodynamics of these sequences. The initial phylogenetic analysis indicated separate monophyletic grouping and there was no mixing of outbreak-associated sequences with sequences from other parts of India. Transmission network analysis using distance-based and non-distance-based methods revealed the existence of transmission clusters within the monophyletic Unnao clade. The median time to the most recent common ancestor (tMRCA) for sequences from Unnao using the pol gene region was observed to be 2011.87 [95% highest posterior density (HPD): 2010.09–2013.53], while the estimates using envelope (env) gene region sequences traced the tMRCA to 2010.33 (95% HPD: 2007.76–2012.99). Phylodynamics estimates demonstrated that the pace of this local epidemic has slowed down in recent times before the time of sampling, but was certainly on an upward track since its inception till 2014.

Statewide Longitudinal Trends in Transmitted HIV-1 Drug Resistance in Rhode Island, USA

Background

HIV-1 transmitted drug resistance (TDR) remains a global challenge that can impact care, yet its comprehensive assessment is limited and heterogenous. We longitudinally characterized statewide TDR in Rhode Island.

Methods

Demographic and clinical data from treatment-naïve individuals were linked to protease, reverse transcriptase, and integrase sequences routinely obtained over 2004-2020. TDR extent, trends, impact on first-line regimens, and association with transmission networks were assessed using the Stanford Database, Mann-Kendall statistic, and phylogenetic tools.

Results

In 1123 individuals, TDR to any antiretroviral increased from 8% (2004) to 26% (2020), driven by non-nucleotide reverse transcriptase inhibitor (NNRTI; 5%-18%) and, to a lesser extent, nucleotide reverse transcriptase inhibitor (NRTI; 2%-8%) TDR. Dual- and triple-class TDR rates were low, and major integrase strand transfer inhibitor resistance was absent. Predicted intermediate to high resistance was in 77% of those with TDR, with differential suppression patterns. Among all individuals, 34% were in molecular clusters, some only with members with TDR who shared mutations. Among clustered individuals, people with TDR were more likely in small clusters.

Conclusions

In a unique (statewide) assessment over 2004-2020, TDR increased; this was primarily, but not solely, driven by NNRTIs, impacting antiretroviral regimens. Limited TDR to multiclass regimens and pre-exposure prophylaxis are encouraging; however, surveillance and its integration with molecular epidemiology should continue in order to potentially improve care and prevention interventions.

Interplay Between Geography and HIV Transmission Clusters in Los Angeles County

BACKGROUND

Clusters of HIV diagnoses in time and space and clusters of genetically linked cases can both serve as alerts for directing prevention and treatment activities. We assessed the interplay between geography and transmission across the Los Angeles County (LAC) HIV genetic transmission network.

METHODS

Deidentified surveillance data reported for 8186 people with HIV residing in LAC from 2010 through 2016 were used to construct a transmission network using HIV-TRACE. We explored geographic assortativity, the tendency for people to link within the same geographic region; concordant time-space pairs, the proportion of genetically linked pairs from the same geographic region and diagnosis year; and Jaccard coefficient, the overlap between geographical and genetic clusters.

RESULTS

Geography was assortative in the genetic transmission network but less so than either race/ethnicity or transmission risk. Only 18% of individuals were diagnosed in the same year and location as a genetically linked partner. Jaccard analysis revealed that cis-men and younger age at diagnosis had more overlap between genetic clusters and geography; the inverse association was observed for trans-women and Blacks/African Americans.

CONCLUSIONS

Within an urban setting with endemic HIV, genetic clustering may serve as a better indicator than time-space clustering to understand HIV transmission patterns and guide public health action.

The Relationship Between the Human Immunodeficiency Virus-1 Transmission Network and the HIV Care Continuum in Los Angeles County

BACKGROUND

Public health action combating human immunodeficiency virus (HIV) includes facilitating navigation through the HIV continuum of care: timely diagnosis followed by linkage to care and initiation of antiretroviral therapy to suppress viral replication. Molecular epidemiology can identify rapidly growing HIV genetic transmission clusters. How progression through the care continuum relates to transmission clusters has not been previously characterized.

METHODS

We performed a retrospective study on HIV surveillance data from 5226 adult cases in Los Angeles County diagnosed from 2010 through 2014. Genetic transmission clusters were constructed using HIV-TRACE. Cox proportional hazard models were used to estimate the impact of transmission cluster growth on the time intervals between care continuum events. Gamma frailty models incorporated the effect of heterogeneity associated with genetic transmission clusters.

RESULTS

In contrast to our expectations, there were no differences in time to the care continuum events among individuals in clusters with different growth dynamics. However, upon achieving viral suppression, individuals in high growth clusters were slower to experience viral rebound (hazard ratio 0.83, P = .011) compared with individuals in low growth clusters. Heterogeneity associated with cluster membership in the timing to each event in the care continuum was highly significant (P < .001), with and without adjustment for transmission risk and demographics.

CONCLUSIONS

Individuals within the same transmission cluster have more similar trajectories through the HIV care continuum than those across transmission clusters. These findings suggest molecular epidemiology can assist public health officials in identifying clusters of individuals who may benefit from assistance in navigating the HIV care continuum.

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.

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

Precise identification of HIV transmission among populations is a key step in public health responses. However, the HIV transmission network is usually difficult to determine. HIV molecular networks can be determined by phylogenetic approach, genetic distance-based approach, and a combination of both approaches. These approaches are increasingly used to identify transmission networks among populations, reconstruct the history of HIV spread, monitor the dynamics of HIV transmission, guide targeted intervention on key subpopulations, and assess the effects of interventions. Simulation and retrospective studies have demonstrated that these molecular network-based interventions are more cost-effective than random or traditional interventions. However, we still need to address several challenges to improve the practice of molecular network-guided targeting interventions to finally end the HIV epidemic. The data remain limited or difficult to obtain, and more automatic real-time tools are required. In addition, molecular and social networks must be combined, and technical parameters and ethnic issues warrant further studies.

HIV Outbreak Control With Effective Access to Care and Harm Reduction in North Carolina, 2017-2018

Objectives. To assess and control a potential outbreak of HIV among people who inject drugs in Western North Carolina.Methods. Disease intervention specialists offered testing for hepatitis B and hepatitis C, harm reduction materials, and linkage to care to 7 linked people recently diagnosed with HIV who also injected drugs. Contacts were offered the same services and HIV testing. HIV genotype analysis was used to characterize HIV spread. We assessed testing and care outcomes by using state surveillance information.Results. Disease intervention specialists contacted 6 of 7 linked group members and received information on 177 contacts; among 96 prioritized contacts, 42 of 96 (44%) were exposed to or diagnosed with hepatitis C, 4 of 96 (4%) had hepatitis B, and 14 of 96 (15%) had HIV (2 newly diagnosed during the investigation). HIV genotype analysis suggested recent transmission to linked group members and 1 contact. Eleven of 14 with HIV were virally suppressed following the outbreak response.Conclusions. North Carolina identified and rapidly responded to an HIV outbreak among people reporting injecting drugs. Effective HIV care, the availability of syringe exchange services, and the rapid response likely contributed to controlling this outbreak.