Identifying Clusters of Recent and Rapid HIV Transmission Through Analysis of Molecular Surveillance Data

Enhancements to the National HIV Surveillance System, United States, 2013-2023

HIV infection is monitored through the National HIV Surveillance System (NHSS) to help improve the health of people with HIV and reduce transmission. NHSS data are routinely used at federal, state, and local levels to monitor the distribution and transmission of HIV, plan and evaluate prevention and care programs, allocate resources, inform policy development, and identify and respond to rapid transmission in the United States. We describe the expanded use of HIV surveillance data since the 2013 NHSS status update, during which time the Centers for Disease Control and Prevention (CDC) coordinated to revise the HIV surveillance case definition to support the detection of early infection and reporting of laboratory data, expanded data collection to include information on sexual orientation and gender identity, enhanced data deduplication processes to improve quality, and expanded reporting to include social determinants of health and health equity measures. CDC maximized the effects of federal funding by integrating funding for HIV prevention and surveillance into a single program; the integration of program funding has expanded the use of HIV surveillance data and strengthened surveillance, resulting in enhanced cluster response capacity and intensified data-to-care activities to ensure sustained viral suppression. NHSS data serve as the primary source for monitoring HIV trends and progress toward achieving national initiatives, including the US Department of Health and Human Services' Ending the HIV Epidemic in the United States initiative, the White House's National HIV/AIDS Strategy (2022-2025), and Healthy People 2030. The NHSS will continue to modernize, adapt, and broaden its scope as the need for high-quality HIV surveillance data remains.

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.

DeepDynaForecast: Phylogenetic-informed graph deep learning for epidemic transmission dynamic prediction

In the midst of an outbreak or sustained epidemic, reliable prediction of transmission risks and patterns of spread is critical to inform public health programs. Projections of transmission growth or decline among specific risk groups can aid in optimizing interventions, particularly when resources are limited. Phylogenetic trees have been widely used in the detection of transmission chains and high-risk populations. Moreover, tree topology and the incorporation of population parameters (phylodynamics) can be useful in reconstructing the evolutionary dynamics of an epidemic across space and time among individuals. We now demonstrate the utility of phylodynamic trees for transmission modeling and forecasting, developing a phylogeny-based deep learning system, referred to as DeepDynaForecast. Our approach leverages a primal-dual graph learning structure with shortcut multi-layer aggregation, which is suited for the early identification and prediction of transmission dynamics in emerging high-risk groups. We demonstrate the accuracy of DeepDynaForecast using simulated outbreak data and the utility of the learned model using empirical, large-scale data from the human immunodeficiency virus epidemic in Florida between 2012 and 2020. Our framework is available as open-source software (MIT license) at github.com/lab-smile/DeepDynaForcast.

Simulation of Full HIV Cluster Networks in a Nationally Representative Model Indicates Intervention Opportunities

BACKGROUND

Clusters of rapid HIV transmission in the United States are increasingly recognized through analysis of HIV molecular sequence data reported to the National HIV Surveillance System. Understanding the full extent of cluster networks is important to assess intervention opportunities. However, full cluster networks include undiagnosed and other infections that cannot be systematically observed in real life.

METHODS

We replicated HIV molecular cluster networks during 2015-2017 in the United States using a stochastic dynamic network simulation model of sexual transmission of HIV. Clusters were defined at the 0.5% genetic distance threshold. Ongoing priority clusters had growth of ≥3 diagnoses/year in multiple years; new priority clusters first had ≥3 diagnoses/year in 2017. We assessed the full extent, composition, and transmission rates of new and ongoing priority clusters.

RESULTS

Full clusters were 3-9 times larger than detected clusters, with median detected cluster sizes in new and ongoing priority clusters of 4 (range 3-9) and 11 (range 3-33), respectively, corresponding to full cluster sizes with a median of 14 (3-74) and 94 (7-318), respectively. A median of 36.3% (range 11.1%-72.6%) of infections in the full new priority clusters were undiagnosed. HIV transmission rates in these clusters were >4 times the overall rate observed in the entire simulation.

CONCLUSIONS

Priority clusters reflect networks with rapid HIV transmission. The substantially larger full extent of these clusters, high proportion of undiagnosed infections, and high transmission rates indicate opportunities for public health intervention and impact.

AUTO-TUNE: SELECTING THE DISTANCE THRESHOLD FOR INFERRING HIV TRANSMISSION CLUSTERS

Molecular surveillance of viral pathogens and inference of transmission networks from genomic data play an increasingly important role in public health efforts, especially for HIV-1. For many methods, the genetic distance threshold used to connect sequences in the transmission network is a key parameter informing the properties of inferred networks. Using a distance threshold that is too high can result in a network with many spurious links, making it difficult to interpret. Conversely, a distance threshold that is too low can result in a network with too few links, which may not capture key insights into clusters of public health concern. Published research using the HIV-TRACE software package frequently uses the default threshold of 0.015 substitutions/site for HIV pol gene sequences, but in many cases, investigators heuristically select other threshold parameters to better capture the underlying dynamics of the epidemic they are studying. Here, we present a general heuristic scoring approach for tuning a distance threshold adaptively, which seeks to prevent the formation of giant clusters. We prioritize the ratio of the sizes of the largest and the second largest cluster, maximizing the number of clusters present in the network. We apply our scoring heuristic to outbreaks with different characteristics, such as regional or temporal variability, and demonstrate the utility of using the scoring mechanism's suggested distance threshold to identify clusters exhibiting risk factors that would have otherwise been more difficult to identify. For example, while we found that a 0.015 substitutions/site distance threshold is typical for US-like epidemics, recent outbreaks like the CRF07_BC subtype among men who have sex with men (MSM) in China have been found to have a lower optimal threshold of 0.005 to better capture the transition from injected drug use (IDU) to MSM as the primary risk factor. Alternatively, in communities surrounding Lake Victoria in Uganda, where there has been sustained hetero-sexual transmission for many years, we found that a larger distance threshold is necessary to capture a more risk factor-diverse population with sparse sampling over a longer period of time. Such identification may allow for more informed intervention action by respective public health officials.

Molecular epidemiology to aid virtual elimination of HIV transmission in Australia

The Global UNAIDS 95/95/95 targets aim to increase the percentage of persons who know their HIV status, receive antiretroviral therapy, and have achieved viral suppression. Achieving these targets requires efforts to improve the public health response to increase access to care for those living with HIV, identify those yet undiagnosed with HIV early, and increase access to prevention for those most at risk of HIV acquisition. HIV infections in Australia are among the lowest globally having recorded significant declines in new diagnoses in the last decade. However, the HIV epidemic has changed with an increasing proportion of newly diagnosed infections among those born outside Australia observed in the last five years. Thus, the current prevention efforts are not enough to achieve the UNAIDS targets and virtual elimination across all population groups. We believe both are possible by including molecular epidemiology in the public health response. Molecular epidemiology methods have been crucial in the field of HIV prevention, particularly in demonstrating the efficacy of treatment as prevention. Cluster detection using molecular epidemiology can provide opportunities for the real-time detection of new outbreaks before they grow, and cluster detection programs are now part of the public health response in the USA and Canada. Here, we review what molecular epidemiology has taught us about HIV evolution and spread. We summarize how we can use this knowledge to improve public health measures by presenting case studies from the USA and Canada. We discuss the successes and challenges of current public health programs in Australia, and how we could use cluster detection as an add-on to identify gaps in current prevention measures easier and respond quicker to growing clusters. Lastly, we raise important ethical and legal challenges that need to be addressed when HIV genotypic data is used in combination with personal data.

Joint modeling HIV and HPV using a new hybrid agent-based network and compartmental simulation technique

Persons living with human immunodeficiency virus (HIV) have a disproportionately higher burden of human papillomavirus infection (HPV)-related cancers. Causal factors include both behavioral and biological. While pharmaceutical and care support interventions help address biological risk of coinfection, as social conditions are common drivers of behaviors, structural interventions are key part of behavioral interventions. Our objective is to develop a joint HIV-HPV model to evaluate the contribution of each factor, to subsequently inform intervention analyses. While compartmental modeling is sufficient for faster spreading HPV, network modeling is suitable for slower spreading HIV. However, using network modeling for jointly modeling HIV and HPV can generate computational complexities given their vastly varying disease epidemiology and disease burden across sub-population groups. We applied a recently developed mixed agent-based compartmental (MAC) simulation technique, which simulates persons with at least one slower spreading disease and their immediate contacts as agents in a network, and all other persons including those with faster spreading diseases in a compartmental model, with an evolving contact network algorithm maintaining the dynamics between the two models. We simulated HIV and HPV in the U.S. among heterosexual female, heterosexual male, and men who have sex with men (men only and men and women) (MSM), sub-populations that mix but have varying HIV burden, and cervical cancer among women. We conducted numerical analyses to evaluate the contribution of behavioral and biological factors to risk of cervical cancer among women with HIV. The model outputs for HIV, HPV, and cervical cancer compared well with surveillance estimates. Model estimates for relative prevalence of HPV (1.67 times) and relative incidence of cervical cancer (3.6 times), among women with HIV compared to women without, were also similar to that reported in observational studies in the literature. The fraction attributed to biological factors ranged from 22-38% for increased HPV prevalence and 80% for increased cervical cancer incidence, the remaining attributed to behavioral. The attribution of both behavioral and biological factors to increased HPV prevalence and cervical cancer incidence suggest the need for behavioral, structural, and pharmaceutical interventions. Validity of model results related to both individual and joint disease metrics serves as proof-of-concept of the MAC simulation technique. Understanding the contribution of behavioral and biological factors of risk helps inform interventions. Future work can expand the model to simulate sexual and care behaviors as functions of social conditions to jointly evaluate behavioral, structural, and pharmaceutical interventions for HIV and cervical cancer prevention.

Molecular HIV Clustering Among Individuals with Mpox and HIV Co-Morbidity in New York State, Excluding New York City

The 2022 global mpox outbreak created an opportunity to test the utility of molecular HIV surveillance (MHS) to identify high-risk transmission networks. Individuals diagnosed with mpox in New York State (NYS) outside New York City-[Rest of State (ROS)] were matched to the NYS HIV and sexually transmitted infection registries. The demographic characteristics of individuals diagnosed with mpox in ROS mirror national trends. HIV-mpox comorbid individuals were more likely to be included in HIV molecular clusters compared to persons living with diagnosed HIV in ROS overall, men who have sex with men (MSM) in ROS, and age-adjusted MSM (to match individuals with mpox diagnosis) in ROS. For the 3-year 0.5% clusters, which are used to define national priority clusters, the HIV-mpox comorbid individuals clustered 2.4 times more frequently than the age/risk-adjusted control group. This study supports the use of HIV MHS to identify populations for priority public health interventions.

Modeling homophily in dynamic networks with application to HIV molecular surveillance

BACKGROUND

Efforts to control the HIV epidemic can benefit from knowledge of the relationships between the characteristics of people who have transmitted HIV and those who became infected by them. Investigation of this relationship is facilitated by the use of HIV genetic linkage analyses, which allows inference about possible transmission events among people with HIV infection. Two persons with HIV (PWH) are considered linked if the genetic distance between their HIV sequences is less than a given threshold, which implies proximity in a transmission network. The tendency of pairs of nodes (in our case PWH) that share (or differ in) certain attributes to be linked is denoted homophily. Below, we describe a novel approach to modeling homophily with application to analyses of HIV viral genetic sequences from clinical series of participants followed in San Diego. Over the 22-year period of follow-up, increases in cluster size results from HIV transmissions to new people from those already in the cluster-either directly or through intermediaries.

METHODS

Our analytical approach makes use of a logistic model to describe homophily with regard to demographic, clinical, and behavioral characteristics-that is we investigate whether similarities (or differences) between PWH in these characteristics are associated with their sequences being linked. To investigate the performance of our methods, we conducted on a simulation study for which data sets were generated in a way that reproduced the structure of the observed database.

RESULTS

Our results demonstrated strong positive homophily associated with hispanic ethnicity, and strong negative homophily, with birth year difference. The second result implies that the larger the difference between the age of a newly-infected PWH and the average age for an available cluster, the lower the odds of a newly infected person joining that cluster. We did not observe homophily associated with prior diagnosis of sexually transmitted diseases. Our simulation studies demonstrated the validity of our approach for modeling homophily, by showing that the estimates it produced matched the specified values of the statistical network generating model.

CONCLUSIONS

Our novel methods provide a simple and flexible statistical network-based approach for modeling the growth of viral (or other microbial) genetic clusters from linkage to new infections based on genetic distance.

Network typologies predict future molecular linkages in the network of HIV transmission

OBJECTIVE

HIV molecular transmission network typologies have previously demonstrated associations to transmission risk; however, few studies have evaluated their predictive potential in anticipating future transmission events. To assess this, we tested multiple models on statewide surveillance data from the Florida Department of Health.

DESIGN

This was a retrospective, observational cohort study examining the incidence of new HIV molecular linkages within the existing molecular network of persons with HIV (PWH) in Florida.

METHODS

HIV-1 molecular transmission clusters were reconstructed for PWH diagnosed in Florida from 2006 to 2017 using the HIV-TRAnsmission Cluster Engine (HIV-TRACE). A suite of machine-learning models designed to predict linkage to a new diagnosis were internally and temporally externally validated using a variety of demographic, clinical, and network-derived parameters.

RESULTS

Of the 9897 individuals who received a genotype within 12 months of diagnosis during 2012-2017, 2611 (26.4%) were molecularly linked to another case within 1 year at 1.5% genetic distance. The best performing model, trained on two years of data, was high performing (area under the receiving operating curve = 0.96, sensitivity = 0.91, and specificity = 0.90) and included the following variables: age group, exposure group, node degree, betweenness, transitivity, and neighborhood.

CONCLUSIONS

In the molecular network of HIV transmission in Florida, individuals' network position and connectivity predicted future molecular linkages. Machine-learned models using network typologies performed superior to models using individual data alone. These models can be used to more precisely identify subpopulations for intervention.

Assessing the reliability of the CD4 depletion model in the presence of Ending the HIV Epidemic initiatives

BACKGROUND

Accurate estimates of HIV incidence are necessary to monitor progress towards Ending the HIV Epidemic (EHE) initiative targets (90% decline by 2030). U.S. incidence estimates are derived from a CD4 depletion model (CD4 model). We performed simulation-based analyses to investigate the ability of this model to estimate HIV incidence when implementing EHE interventions that have the potential to shorten the duration between HIV infection and diagnosis (diagnosis delay).

METHODS

Our simulation study evaluates the impact of three parameters on the accuracy of incidence estimates derived from the CD4 model: rate of HIV incidence decline, length of diagnosis delay, and sensitivity of using CD4 + cell counts to identify new infections (recency error). We model HIV incidence and diagnoses after the implementation of a theoretical prevention intervention and compare HIV incidence estimates derived from the CD4 model to simulated incidence.

RESULTS

Theoretical interventions that shortened the diagnosis delay (10-50%) result in overestimation of HIV incidence by the CD4 model (10-92%) in the first year and by more than 10% for the first 6 years after implementation of the intervention. Changes in the rate of HIV incidence decline and the presence of recency error had minimal impact on the accuracy of incidence estimates derived from the CD4 model.

CONCLUSION

In the setting of EHE interventions to identify persons with HIV earlier during infection, the CD4 model overestimates HIV incidence. Alternative methods to estimate incidence based on objective measures of incidence are needed to assess and monitor EHE interventions.

Public Health Response to Clusters of Rapid HIV Transmission Among Hispanic or Latino Gay, Bisexual, and Other Men Who Have Sex with Men - Metropolitan Atlanta, Georgia, 2021-2022

During February 2021-June 2022, the Georgia Department of Public Health (GDPH) detected five clusters of rapid HIV transmission concentrated among Hispanic or Latino (Hispanic) gay, bisexual, and other men who have sex with men (MSM) in metropolitan Atlanta. The clusters were detected through routine analysis of HIV-1 nucleotide sequence data obtained through public health surveillance (1,2). Beginning in spring 2021, GDPH partnered with health districts with jurisdiction in four metropolitan Atlanta counties (Cobb, DeKalb, Fulton, and Gwinnett) and CDC to investigate factors contributing to HIV spread, epidemiologic characteristics, and transmission patterns. Activities included review of surveillance and partner services interview data,† medical chart reviews, and qualitative interviews with service providers and Hispanic MSM community members. By June 2022, these clusters included 75 persons, including 56% who identified as Hispanic, 96% who reported male sex at birth, 81% who reported male-to-male sexual contact, and 84% of whom resided in the four metropolitan Atlanta counties. Qualitative interviews identified barriers to accessing HIV prevention and care services, including language barriers, immigration- and deportation-related concerns, and cultural norms regarding sexuality-related stigma. GDPH and the health districts expanded coordination, initiated culturally concordant HIV prevention marketing and educational activities, developed partnerships with organizations serving Hispanic communities to enhance outreach and services, and obtained funding for a bilingual patient navigation program with academic partners to provide staff members to help persons overcome barriers and understand the health care system. HIV molecular cluster detection can identify rapid HIV transmission among sexual networks involving ethnic and sexual minority groups, draw attention to the needs of affected populations, and advance health equity through tailored responses that address those needs.

Not all clusters are equal: dynamics of molecular HIV-1 clusters in a statewide Rhode Island epidemic

OBJECTIVES

Molecular epidemiology is a powerful tool to characterize HIV epidemics and prioritize public health interventions. Typically, HIV clusters are assumed to have uniform patterns over time. We hypothesized that assessment of cluster evolution would reveal distinct cluster behavior, possibly improving molecular epidemic characterization, towards disrupting HIV transmission.

DESIGN

Retrospective cohort.

METHODS

Annual phylogenies were inferred by cumulative aggregation of all available HIV-1 pol sequences of individuals with HIV-1 in Rhode Island (RI) between 1990 and 2020, representing a statewide epidemic. Molecular clusters were detected in annual phylogenies by strict and relaxed cluster definition criteria, and the impact of annual newly-diagnosed HIV-1 cases to the structure of individual clusters was examined over time.

RESULTS

Of 2153 individuals, 31% (strict criteria) - 47% (relaxed criteria) clustered. Longitudinal tracking of individual clusters identified three cluster types: normal, semi-normal and abnormal. Normal clusters (83-87% of all identified clusters) showed predicted growing/plateauing dynamics, with approximately three-fold higher growth rates in large (15-18%) vs. small (∼5%) clusters. Semi-normal clusters (1-2% of all clusters) temporarily fluctuated in size and composition. Abnormal clusters (11-16% of all clusters) demonstrated collapses and re-arrangements over time. Borderline values of cluster-defining parameters explained dynamics of non-normal clusters.

CONCLUSIONS

Comprehensive tracing of molecular HIV clusters over time in a statewide epidemic identified distinct cluster types, likely missed in cross-sectional analyses, demonstrating that not all clusters are equal. This knowledge challenges current perceptions of consistent cluster behavior over time and could improve molecular surveillance of local HIV epidemics to better inform public health strategies.

HIV-1 Transmission Cluster in Injection Drug Users in Nairobi City, Kenya

Background

While there is a striking increase in the prevalence of HIV in injection drug users, information on envelope-gene subtypes and transmission clusters in injection drug users is scarce.

Method

In a cross-sectional study, 247 injection drug users were recruited via out-rich method. Deoxyribonucleic acid was extracted from dry blood spot samples, amplified by Polymerase Chain Reaction and sequenced. Subtyping was performed using COntext-based Modeling for Expeditious Typing (COMET) and Recombinant Identification Program (RIP) tools. Phylogenetic diversity and Transmission clusters were identified using MEGA version 6.0 and TreeLink, respectively.

Results

Overall, 42 (17.0%) injection drug users were sero-positive for HIV-1. Of the 37 samples successfully sequenced, 29 (78.4%) sequences were identified as A1, 6 (16.2%) as AG while 1 (2.7%) as A1/G/AE and A1/C recombinants. The HIV subtypes formed clusters with little genetic diversity.

Conclusion

The high HIV prevalence was associated with transmission clusters and diversity in subtypes indicating ongoing local transmission. Therefore, there is need for comprehensive HIV care tailored to this population.

Quantifying and Predicting Ongoing Human Immunodeficiency Virus Type 1 Transmission Dynamics in Switzerland Using a Distance-Based Clustering Approach

BACKGROUND

Despite effective prevention approaches, ongoing human immunodeficiency virus 1 (HIV-1) transmission remains a public health concern indicating a need for identifying its drivers.

METHODS

We combined a network-based clustering method using evolutionary distances between viral sequences with statistical learning approaches to investigate the dynamics of HIV transmission in the Swiss HIV Cohort Study and to predict the drivers of ongoing transmission.

RESULTS

We found that only a minority of clusters and patients acquired links to new infections between 2007 and 2020. While the growth of clusters and the probability of individual patients acquiring new links in the transmission network was associated with epidemiological, behavioral, and virological predictors, the strength of these associations decreased substantially when adjusting for network characteristics. Thus, these network characteristics can capture major heterogeneities beyond classical epidemiological parameters. When modeling the probability of a newly diagnosed patient being linked with future infections, we found that the best predictive performance (median area under the curve receiver operating characteristic AUCROC = 0.77) was achieved by models including characteristics of the network as predictors and that models excluding them performed substantially worse (median AUCROC = 0.54).

CONCLUSIONS

These results highlight the utility of molecular epidemiology-based network approaches for analyzing and predicting ongoing HIV transmission dynamics. This approach may serve for real-time prospective assessment of HIV transmission.

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.

Insights into HIV-1 Transmission Dynamics Using Routinely Collected Data in the Mid-Atlantic United States

Background: Molecular epidemiological approaches provide opportunities to characterize HIV transmission dynamics. We analyzed HIV sequences and virus load (VL) results obtained during routine clinical care, and individual’s zip-code location to determine utility of this approach. Methods: HIV-1 pol sequences aligned using ClustalW were subtyped using REGA. A maximum likelihood (ML) tree was generated using IQTree. Transmission clusters with ≤3% genetic distance (GD) and ≥90% bootstrap support were identified using ClusterPicker. We conducted Bayesian analysis using BEAST to confirm transmission clusters. The proportion of nucleotides with ambiguity ≤0.5% was considered indicative of early infection. Descriptive statistics were applied to characterize clusters and group comparisons were performed using chi-square or t-test. Results: Among 2775 adults with data from 2014−2015, 2589 (93%) had subtype B HIV-1, mean age was 44 years (SD 12.7), 66.4% were male, and 25% had nucleotide ambiguity ≤0.5. There were 456 individuals in 193 clusters: 149 dyads, 32 triads, and 12 groups with ≥ four individuals per cluster. More commonly in clusters were males than females, 349 (76.5%) vs. 107 (23.5%), p < 0.0001; younger individuals, 35.3 years (SD 12.1) vs. 44.7 (SD 12.3), p < 0.0001; and those with early HIV-1 infection by nucleotide ambiguity, 202/456 (44.3%) vs. 442/2133 (20.7%), p < 0.0001. Members of 43/193 (22.3%) of clusters included individuals in different jurisdictions. Clusters ≥ four individuals were similarly found using BEAST. HIV-1 viral load (VL) ≥3.0 log10 c/mL was most common among individuals in clusters ≥ four, 18/21, (85.7%) compared to 137/208 (65.8%) in clusters sized 2−3, and 927/1169 (79.3%) who were not in a cluster (p < 0.0001). Discussion: HIV sequence data obtained for HIV clinical management provide insights into regional transmission dynamics. Our findings demonstrate the additional utility of HIV-1 VL data in combination with phylogenetic inferences as an enhanced contact tracing tool to direct HIV treatment and prevention services. Trans-jurisdictional approaches are needed to optimize efforts to end the HIV epidemic.

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.

Impact of molecular sequence data completeness on HIV cluster detection and a network science approach to enhance detection

Detection of viral transmission clusters using molecular epidemiology is critical to the response pillar of the Ending the HIV Epidemic initiative. Here, we studied whether inference with an incomplete dataset would influence the accuracy of the reconstructed molecular transmission network. We analyzed viral sequence data available from ~ 13,000 individuals with diagnosed HIV (2012-2019) from Houston Health Department surveillance data with 53% completeness (n = 6852 individuals with sequences). We extracted random subsamples and compared the resulting reconstructed networks versus the full-size network. Increasing simulated completeness was associated with an increase in the number of detected clusters. We also subsampled based on the network node influence in the transmission of the virus where we measured Expected Force (ExF) for each node in the network. We simulated the removal of nodes with the highest and then lowest ExF from the full dataset and discovered that 4.7% and 60% of priority clusters were detected respectively. These results highlight the non-uniform impact of capturing high influence nodes in identifying transmission clusters. Although increasing sequence reporting completeness is the way to fully detect HIV transmission patterns, reaching high completeness has remained challenging in the real world. Hence, we suggest taking a network science approach to enhance performance of molecular cluster detection, augmented by node influence information.

A deep learning approach to real-time HIV outbreak detection using genetic data

Pathogen genomic sequence data are increasingly made available for epidemiological monitoring. A main interest is to identify and assess the potential of infectious disease outbreaks. While popular methods to analyze sequence data often involve phylogenetic tree inference, they are vulnerable to errors from recombination and impose a high computational cost, making it difficult to obtain real-time results when the number of sequences is in or above the thousands. Here, we propose an alternative strategy to outbreak detection using genomic data based on deep learning methods developed for image classification. The key idea is to use a pairwise genetic distance matrix calculated from viral sequences as an image, and develop convolutional neutral network (CNN) models to classify areas of the images that show signatures of active outbreak, leading to identification of subsets of sequences taken from an active outbreak. We showed that our method is efficient in finding HIV-1 outbreaks with R0 ≥ 2.5, and overall a specificity exceeding 98% and sensitivity better than 92%. We validated our approach using data from HIV-1 CRF01 in Europe, containing both endemic sequences and a well-known dual outbreak in intravenous drug users. Our model accurately identified known outbreak sequences in the background of slower spreading HIV. Importantly, we detected both outbreaks early on, before they were over, implying that had this method been applied in real-time as data became available, one would have been able to intervene and possibly prevent the extent of these outbreaks. This approach is scalable to processing hundreds of thousands of sequences, making it useful for current and future real-time epidemiological investigations, including public health monitoring using large databases and especially for rapid outbreak identification.

Characterising persons diagnosed with HIV as either recent or long-term using a cross-sectional analysis of recent infection surveillance data collected in Malawi from September 2019 to March 2020

OBJECTIVES

In Malawi, a recent infection testing algorithm (RITA) is used to characterise infections of persons newly diagnosed with HIV as recent or long term. This paper shares results from recent HIV infection surveillance and describes distribution and predictors.

SETTING

Data from 155 health facilities in 11 districts in Malawi were pooled from September 2019 to March 2020.

PARTICIPANTS

Eligible participants were ≥13 years, and newly diagnosed with HIV. Clients had RITA recent infections if the rapid test for recent infection (RTRI) test result was recent and viral load (VL) ≥1000 copies/mL; if VL was <1000 copies/mL the RTRI result was reclassified as long-term. Results were stratified by age, sex, pregnancy/breastfeeding status and district.

RESULTS

13 838 persons consented to RTRI testing and 12 703 had valid RTRI test results and VL results after excluding clients not newly HIV-positive, RTRI negative or missing data (n=1135). A total of 12 365 of the 12 703 were included in the analysis after excluding those whose RTRI results were reclassified as long term (n=338/784 or 43.1%). The remainder, 446/12 703 or 3.5%, met the definition of RITA recent infection. The highest percentage of recent infections was among breastfeeding women (crude OR (COR) 3.2; 95% CI 2.0 to 5.0), young people aged 15-24 years (COR 1.6; 95% CI 1.3 to 1.9) and persons who reported a negative HIV test within the past 12 months (COR 3.3; 95% CI 2.6 to 4.2). Factors associated with recent infection in multivariable analysis included being a non-pregnant female (adjusted OR (AOR) 1.4; 95% CI 1.2 to 1.8), a breastfeeding female (AOR 2.2; 95% CI 1.4 to 3.5), aged 15-24 years (AOR 1.6; 95% CI 1.3 to 1.9) and residents of Machinga (AOR 2.0; 95% CI 1.2 to 3.5) and Mzimba (AOR 2.4; 95% CI 1.3 to 4.5) districts.

CONCLUSIONS

Malawi's recent HIV infection surveillance system demonstrated high uptake and identified sub-populations of new HIV diagnoses with a higher percentage of recent infections.

Genetic network analysis of human immunodeficiency virus sexual transmission in rural Southwest China after the expansion of antiretroviral therapy: A population-based study

Background

This study is used to analyze the genetic network of HIV sexual transmission in rural areas of Southwest China after expanding antiretroviral therapy (ART) and to investigate the factors associated with HIV sexual transmission through the genetic network.

Materials and methods

This was a longitudinal genetic network study in Guangxi, China. The baseline survey and follow-up study were conducted among patients with HIV in 2015, and among those newly diagnosed from 2016 to 2018, respectively. A generalized estimating equation model was employed to explore the factors associated with HIV transmission through the genetic linkage between newly diagnosed patients with HIV (2016-2018) and those at baseline (2015-2017), respectively.

Results

Of 3,259 identified HIV patient sequences, 2,714 patients were at baseline, and 545 were newly diagnosed patients with HIV at follow-up. A total of 8,691 baseline objectives were observed by repeated measurement analysis. The prevention efficacy in HIV transmission for treated HIV patients was 33% [adjusted odds ratio (AOR): 0.67, 95% confidence interval (CI): 0.48-0.93]. Stratified analyses indicated the prevention efficacy in HIV transmission for treated HIV patients with a viral load (VL) of <50 copies/ml and those treated for 4 years with a VL of <50 copies/ml to be 41 [AOR: 0.59, 95% CI: 0.43-0.82] and 65% [AOR: 0.35, 95% CI: 0.24-0.50], respectively. No significant reduction in HIV transmission occurred among treated HIV patients with VL missing or treated HIV patients on dropout. Some factors were associated with HIV transmission, including over 50 years old, men, Zhuang and other nationalities, with less than secondary schooling, working as a farmer, and heterosexual transmission.

Conclusion

This study reveals the role of ART in reducing HIV transmission, and those older male farmers with less than secondary schooling are at high risk of HIV infection at a population level. Improvements to ART efficacy for patients with HIV and precision intervention on high-risk individuals during the expansion of ART are urgently required.

Use of Injection Drugs and Any Form of Methamphetamine in the Portland, OR Metro Area as a Driver of an HIV Time-Space Cluster: Clackamas, Multnomah, and Washington Counties, 2018-2020

We describe the response to detection of a time-space cluster of new HIV infection in the Portland, OR metro area among people who inject drugs (PWID) and/or people who use any form of methamphetamine. This time-space cluster took place in a region with a syndemic of homelessness and drug use. The investigation included new HIV diagnoses in 2018, 2019, and 2020 in Clackamas, Multnomah, and Washington Counties. Public health response included activating incident command, development and implementation of an enhanced interview tool, outreach testing, and stakeholder engagement. We identified 396 new cases of HIV infection, 116 (29%) of which met the cluster definition. Most cluster cases had no molecular relationships to previous cases. Persons responding to the enhanced interview tool reported behaviors associated with HIV acquisition. Field outreach testing did not identify any new HIV cases but did identify hepatitis C and syphilis infections. We show the importance of a robust public health response to a time-space cluster of new HIV infections in an urban area.

Using molecular network analysis to explore the characteristics of HIV-1 transmission in a China-Myanmar border area

Background

The China-Myanmar border area is considered a hot spot of active HIV-1 recombination in Southeast Asia. To better understand the characteristics of HIV-1 transmission in this area, a cross-sectional HIV-1 molecular epidemiological survey was conducted in Baoshan Prefecture of Yunnan Province.

Methods

In total, 708 newly reported HIV-1 cases in Baoshan Prefecture from 2019 to 2020 were included in this study. HIV-1 gag, pol and env genes were sequenced, and the spatial and demographic distributions of HIV-1 genotypes were analyzed. The characteristics of HIV-1 transmission were investigated using the HIV-1 molecular network method.

Results

In the 497 samples with genotyping results, 19 HIV-1 genotypes were found, with URFs being the predominant strains (30.2%, 150/497). The main circulating HIV-1 strains were mostly distributed in the northern area of Baoshan. URFs were more likely identified in Burmese individuals, intravenous drug users and those younger than 50 years old. CRF08_BC was more likely detected in farmers and those of Han ethnicity, CRF01_AE in the young and those of Han ethnicity, and CRF07_BC in the subpopulation with junior middle school education and higher. Moreover, CRF118_BC and CRF64_BC were more likely found in the subpopulation aged ≥40 years and ≥50 years, respectively. Among 480 individuals with pol sequence detection, 179 (37.3%) were grouped into 78 clusters, with Baoshan natives being more likely to be in the network. The proportion of the linked individuals showed significant differences when stratified by the regional origin, marital status, age and county of case reporting. In the molecular network, recent infections were more likely to occur among nonfarmers and individuals aged below 30 years.

Conclusions

HIV-1 genetics has become complex in Baoshan. HIV-1 molecular network analysis provided transmission characteristics in the local area, and these findings provided information to prioritize transmission-reduction interventions.

Do partner services linked to molecular clusters yield people with viremia or new HIV?

OBJECTIVES

We examined whether molecular cluster membership was associated with public health identification of HIV transmission potential among named partners in Chicago.

DESIGN

Historical cohort study.

METHODS

We matched and analyzed HIV surveillance and partner services data from HIV diagnoses (2012-2016) prior to implementation of cluster detection and response interventions. We constructed molecular clusters using HIV-TRACE at a pairwise genetic distance threshold of 0.5% and identified clusters exhibiting recent and rapid growth according to the Centers for Disease Control and Prevention definition (three new cases diagnosed in past year). Factors associated with identification of partners with HIV transmission potential were examined using multivariable Poisson regression.

RESULTS

There were 5208 newly diagnosed index clients over this time period. Average age of index clients in clusters was 28; 47% were Black, 29% Latinx/Hispanic, 6% female and 89% MSM. Of the 537 named partners, 191 (35.6%) were linked to index cases in a cluster and of those 16% were either new diagnoses or viremic. There was no statistically significant difference in the probability of identifying partners with HIV transmission potential among index clients in a rapidly growing cluster versus those not in a cluster [adjusted relative risk 1.82, (0.81-4.06)].

CONCLUSION

Partner services that were initiated from index clients in a molecular cluster yielded similar new HIV case finding or identification of those with viremia as did interviews with index clients not in clusters. It remains unclear whether these findings are due to temporal disconnects between diagnoses and cluster identification, unobserved cluster members, or challenges with partner services implementation.

Beyond HIV outbreaks: protocol, rationale and implementation of a prospective study quantifying the benefit of incorporating viral sequence clustering analysis into routine public health interventions

INTRODUCTION

HIV continues to have great impact on millions of lives. Novel methods are needed to disrupt HIV transmission networks. In the USA, public health departments routinely conduct contact tracing and partner services and interview newly HIV-diagnosed index cases to obtain information on social networks and guide prevention interventions. Sequence clustering methods able to infer HIV networks have been used to investigate and halt outbreaks. Incorporation of such methods into routine, not only outbreak-driven, contact tracing and partner services holds promise for further disruption of HIV transmissions.

METHODS AND ANALYSIS

Building on a strong academic-public health collaboration in Rhode Island, we designed and have implemented a state-wide prospective study to evaluate an intervention that incorporates real-time HIV molecular clustering information with routine contact tracing and partner services. We present the rationale and study design of our approach to integrate sequence clustering methods into routine public health interventions as well as related important ethical considerations. This prospective study addresses key questions about the benefit of incorporating a clustering analysis triggered intervention into the routine workflow of public health departments, going beyond outbreak-only circumstances. By developing an intervention triggered by, and incorporating information from, viral sequence clustering analysis, and evaluating it with a novel design that avoids randomisation while allowing for methods comparison, we are confident that this study will inform how viral sequence clustering analysis can be routinely integrated into public health to support the ending of the HIV pandemic in the USA and beyond.

ETHICS AND DISSEMINATION

The study was approved by both the Lifespan and Rhode Island Department of Health Human Subjects Research Institutional Review Boards and study results will be published in peer-reviewed journals.

Early Emergence Phase of SARS-CoV-2 Delta Variant in Florida, US

SARS-CoV-2, the causative agent of COVID-19, emerged in late 2019. The highly contagious B.1.617.2 (Delta) variant of concern (VOC) was first identified in October 2020 in India and subsequently disseminated worldwide, later becoming the dominant lineage in the US. Understanding the local transmission dynamics of early SARS-CoV-2 introductions may inform actionable mitigation efforts during subsequent pandemic waves. Yet, despite considerable genomic analysis of SARS-CoV-2 in the US, several gaps remain. Here, we explore the early emergence of the Delta variant in Florida, US using phylogenetic analysis of representative Florida and globally sampled genomes. We find multiple independent introductions into Florida primarily from North America and Europe, with a minority originating from Asia. These introductions led to three distinct clades that demonstrated varying relative rates of transmission and possessed five distinct substitutions that were 3-21 times more prevalent in the Florida sample as compared to the global sample. Our results underscore the benefits of routine viral genomic surveillance to monitor epidemic spread and support the need for more comprehensive genomic epidemiology studies of emerging variants. In addition, we provide a model of epidemic spread of newly emerging VOCs that can inform future public health responses.

Transmission Clusters, Predominantly Associated With Men Who Have Sex With Men, Play a Main Role in the Propagation of HIV-1 in Northern Spain (2013–2018)

Viruses of HIV-1-infected individuals whose transmission is related group phylogenetically in transmission clusters (TCs). The study of the phylogenetic relations of these viruses and the factors associated with these individuals is essential to analyze the HIV-1 epidemic. In this study, we examine the role of TCs in the epidemiology of HIV-1 infection in Galicia and the Basque County, two regions of northern Spain. A total of 1,158 HIV-1-infected patients from both regions with new diagnoses (NDs) in 2013–2018 were included in the study. Partial HIV-1 pol sequences were analyzed phylogenetically by approximately maximum-likelihood with FastTree 2. In this analysis, 10,687 additional sequences from samples from HIV-1-infected individuals collected in Spain in 1999–2019 were also included to assign TC membership and to determine TCs’ sizes. TCs were defined as those which included viruses from ≥4 individuals, at least 50% of them Spaniards, and with ≥0.95 Shimodaira-Hasegawa-like node support in the phylogenetic tree. Factors associated to TCs were evaluated using odds ratios (OR) and their 95% CI. Fifty-one percent of NDs grouped in 162 TCs. Male patients (OR: 2.6; 95% CI: 1.5–4.7) and men having sex with men (MSM; OR: 2.1; 95% CI: 1.4–3.2) had higher odds of belonging to a TC compared to female and heterosexual patients, respectively. Individuals from Latin America (OR: 0.3; 95% CI: 0.2–0.4), North Africa (OR: 0.4; 95% CI: 0.2–1.0), and especially Sub-Saharan Africa (OR: 0.02; 95% CI: 0.003–0.2) were inversely associated to belonging to TCs compared to native Spaniards. Our results show that TCs are important components of the HIV-1 epidemics in the two Spanish regions studied, where transmission between MSM is predominant. The majority of migrants were infected with viruses not belonging to TCs that expand in Spain. Molecular epidemiology is essential to identify local peculiarities of HIV-1 propagation. The early detection of TCs and prevention of their expansion, implementing effective control measures, could reduce HIV-1 infections.

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.

Dating the Origin and Estimating the Transmission Rates of the Major HIV-1 Clusters in Greece: Evidence about the Earliest Subtype A1 Epidemic in Europe

Our aim was to estimate the date of the origin and the transmission rates of the major local clusters of subtypes A1 and B in Greece. Phylodynamic analyses were conducted in 14 subtype A1 and 31 subtype B clusters. The earliest dates of origin for subtypes A1 and B were in 1982.6 and in 1985.5, respectively. The transmission rate for the subtype A1 clusters ranged between 7.54 and 39.61 infections/100 person years (IQR: 9.39, 15.88), and for subtype B clusters between 4.42 and 36.44 infections/100 person years (IQR: 7.38, 15.04). Statistical analysis revealed that the average difference in the transmission rate between the PWID and the MSM clusters was 6.73 (95% CI: 0.86 to 12.60; p = 0.026). Our study provides evidence that the date of introduction of subtype A1 in Greece was the earliest in Europe. Transmission rates were significantly higher for PWID than MSM clusters due to the conditions that gave rise to an extensive PWID HIV-1 outbreak ten years ago in Athens, Greece. Transmission rate can be considered as a valuable measure for public health since it provides a proxy of the rate of epidemic growth within a cluster and, therefore, it can be useful for targeted HIV prevention programs.

Phylogenetic prioritization of HIV-1 transmission clusters with viral lineage-level diversification rates

Background and objectives

Public health officials faced with a large number of transmission clusters require a rapid, scalable and unbiased way to prioritize distribution of limited resources to maximize benefits. We hypothesize that transmission cluster prioritization based on phylogenetically derived lineage-level diversification rates will perform as well as or better than commonly used growth-based prioritization measures, without need for historical data or subjective interpretation.

Methodology

9822 HIV pol sequences collected during routine drug resistance genotyping were used alongside simulated sequence data to infer sets of phylogenetic transmission clusters via patristic distance threshold. Prioritized clusters inferred from empirical data were compared to those prioritized by the current public health protocols. Prioritization of simulated clusters was evaluated based on correlation of a given prioritization measure with future cluster growth, as well as the number of direct downstream transmissions from cluster members.

Results

Empirical data suggest diversification rate-based measures perform comparably to growth-based measures in recreating public heath prioritization choices. However, unbiased simulated data reveals phylogenetic diversification rate-based measures perform better in predicting future cluster growth relative to growth-based measures, particularly long-term growth. Diversification rate-based measures also display advantages over growth-based measures in highlighting groups with greater future transmission events compared to random groups of the same size. Furthermore, diversification rate measures were notably more robust to effects of decreased sampling proportion.

Conclusions and implications

Our findings indicate diversification rate-based measures frequently outperform growth-based measures in predicting future cluster growth and offer several additional advantages beneficial to optimizing the public health prioritization process.

Forecasting HIV-1 Genetic Cluster Growth in Illinois,United States

BACKGROUND

HIV intervention activities directed toward both those most likely to transmit and their HIV-negative partners have the potential to substantially disrupt HIV transmission. Using HIV sequence data to construct molecular transmission clusters can reveal individuals whose viruses are connected. The utility of various cluster prioritization schemes measuring cluster growth have been demonstrated using surveillance data in New York City and across the United States, by the Centers for Disease Control and Prevention (CDC).

METHODS

We examined clustering and cluster growth prioritization schemes using Illinois HIV sequence data that include cases from Chicago, a large urban center with high HIV prevalence, to compare their ability to predict future cluster growth.

RESULTS

We found that past cluster growth was a far better predictor of future cluster growth than cluster membership alone but found no substantive difference between the schemes used by CDC and the relative cluster growth scheme previously used in New York City (NYC). Focusing on individuals selected simultaneously by both the CDC and the NYC schemes did not provide additional improvements.

CONCLUSION

Growth-based prioritization schemes can easily be automated in HIV surveillance tools and can be used by health departments to identify and respond to clusters where HIV transmission may be actively occurring.

Performance evaluation of the Asante Rapid Recency Assay for verification of HIV diagnosis and detection of recent HIV-1 infections: Implications for epidemic control

We previously described development of a rapid test for recent infection (RTRI) that can diagnose HIV infection and detect HIV-1 recent infections in a single device. This technology was transferred to a commercial partner as Asante Rapid Recency Assay (ARRA). We evaluated performance of the ARRA kits in the laboratory using a well-characterized panel of specimens. The plasma specimen panel (N = 1500) included HIV-1 (N = 570), HIV-2 (N = 10), and HIV-negatives (N = 920) representing multiple subtypes and geographic locations. Reference diagnostic data were generated using the Bio-Rad HIV-1-2-O EIA/Western blot algorithm with further serotyping performed using the Multispot HIV-1/2 assay. The LAg-Avidity EIA was used to generate reference data on recent and long-term infection for HIV-1 positive specimens at a normalized optical density (ODn) cutoff of 2.0 corresponding to a mean duration of about 6 months. All specimens were tested with ARRA according to the manufacturer's recommendations. Test strips were also read for line intensities using a reader and results were correlated with visual interpretation. ARRA's positive verification line (PVL) correctly classified 575 of 580 HIV-positive and 910 of 920 negative specimens resulting in a sensitivity of 99.1% (95% CI: 98.0-99.6) and specificity of 98.9% (95% CI: 98.1-99.4), respectively. The reader-based classification was similar for PVL with sensitivity of 99.3% (576/580) and specificity of 98.8% (909/920). ARRA's long-term line (LTL) classified 109 of 565 HIV-1 specimens as recent and 456 as long-term compared to 98 as recent and 467 as long-term (LT) by LAg-Avidity EIA (cutoff ODn = 2.0), suggesting a mean duration of recent infection (MDRI) close to 6 months. Agreement of ARRA with LAg recent cases was 81.6% (80/98) and LT cases was 93.8% (438/467), with an overall agreement of 91.7% (kappa = 0.72). The reader (cutoff 2.9) classified 109/566 specimens as recent infections compared to 99 by the LAg-Avidity EIA for recency agreement of 81.8% (81/99), LT agreement of 9% (439/467) with overall agreement of 91.9% (kappa = 0.72). The agreement between visual interpretation and strip reader was 99.9% (95% CI: 99.6-99.9) for the PVL and 98.1% (95% CI: 96.6-98.9) for the LTL. ARRA performed well with HIV diagnostic sensitivity >99% and specificity >98%. Its ability to identify recent infections is comparable to the LA-Avidity EIA corresponding to an MDRI of about 6 months. This point-of-care assay has implications for real-time surveillance of new infections among newly diagnosed individuals for targeted prevention and interrupting ongoing transmission thus accelerating epidemic control.

HIV Drug Resistance and Transmission Networks Among a Justice-Involved Population at the Time of Community Reentry in Washington, D.C

Justice-involved (JI) populations bear a disproportionate burden of HIV infection and are at risk of poor treatment outcomes. Drug resistance prevalence and emergence, and phylogenetic inference of transmission networks, understudied in vulnerable JI populations, can inform care and prevention interventions, particularly around the critical community reentry period. We analyzed banked blood specimens from CARE+ Corrections study participants in Washington, D.C. (DC) across three time points and conducted HIV drug resistance testing using next-generation sequencing (NGS) at 20% and 5% thresholds to identify prevalent and evolving resistance during community reentry. Phylogenetic analysis was used to identify molecular clusters within participants, and in an extended analysis between participants and publicly available DC sequences. HIV sequence data from 54 participants (99 specimens) were analyzed. The prevalence of transmitted drug resistance was 14% at both thresholds, and acquired drug resistance was 47% at 20%, and 57% at 5% NGS thresholds, respectively. The overall prevalence of drug resistance was 43% at 20%, and 52% at 5% NGS thresholds, respectively. Among 34 participants sampled longitudinally, 21%–35% accumulated 10–17 new resistance mutations during a mean 4.3 months. In phylogenetic analysis within the JI population, 11% were found in three molecular clusters. The extended phylogenetic analysis identified 46% of participants in 22 clusters, of which 21 also included publicly-available DC sequences, and one JI-only unique dyad. This is the first study to identify a high prevalence of HIV drug resistance and its accumulation in a JI population during community reentry and suggests phylogenetic integration of this population into the non-JI DC HIV community. These data support the need for new, effective, and timely interventions to improve HIV treatment during this vulnerable period, and for JI populations to be included in broader surveillance and prevention efforts.

Inferring population HIV incidence trends from surveillance data of recent HIV infection among HIV testing clients

BACKGROUND

Measuring recent HIV infections from routine surveillance systems could allow timely and granular monitoring of HIV incidence patterns. We evaluated the relationship of two recent infection indicators with alternative denominators to true incidence patterns.

METHODS

We used a mathematical model of HIV testing behaviours, calibrated to population-based surveys and HIV testing services programme data, to estimate the number of recent infections diagnosed annually from 2010 to 2019 in Côte d'Ivoire, Malawi, and Mozambique. We compared two different denominators to interpret recency data: those at risk of HIV acquisition (HIV-negative tests and recent infections) and all people testing HIV positive. Sex and age-specific longitudinal trends in both interpretations were then compared with modelled trends in HIV incidence, testing efforts and HIV positivity among HIV testing services clients.

RESULTS

Over 2010-2019, the annual proportion of the eligible population tested increased in all countries, while positivity decreased. The proportion of recent infections among those at risk of HIV acquisition decreased, similar to declines in HIV incidence among adults (≥15 years old). Conversely, the proportion of recent infections among HIV-positive tests increased. The female-to-male ratio of the proportion testing recent among those at risk was closer to 1 than the true incidence sex ratio.

CONCLUSION

The proportion of recent infections among those at risk of HIV acquisition is more indicative of HIV incidence than the proportion among HIV-positive tests. However, interpreting the observed patterns as surrogate measures for incidence patterns may still be confounded by different HIV testing rates between population groups or over time.

Characteristics and growth of the genetic HIV transmission network of Mexico City during 2020

INTRODUCTION

Molecular surveillance systems could provide public health benefits to focus strategies to improve the HIV care continuum. Here, we infer the HIV genetic network of Mexico City in 2020, and identify actively growing clusters that could represent relevant targets for intervention.

METHODS

All new diagnoses, referrals from other institutions, as well as persons returning to care, enrolling at the largest HIV clinic in Mexico City were invited to participate in the study. The network was inferred from HIV pol sequences, using pairwise genetic distance methods, with a locally hosted, secure version of the HIV-TRACE tool: Seguro HIV-TRACE. Socio-demographic, clinical and behavioural metadata were overlaid across the network to design focused prevention interventions.

RESULTS

A total of 3168 HIV sequences from unique individuals were included. One thousand and one-hundred and fifty (36%) sequences formed 1361 links within 386 transmission clusters in the network. Cluster size varied from 2 to 14 (63% were dyads). After adjustment for covariates, lower age (adjusted odds ratio [aOR]: 0.37, p<0.001; >34 vs. <24 years), being a man who has sex with men (MSM) (aOR: 2.47, p = 0.004; MSM vs. cisgender women), having higher viral load (aOR: 1.28, p<0.001) and higher CD4+ T cell count (aOR: 1.80, p<0.001; ≥500 vs. <200 cells/mm3 ) remained associated with higher odds of clustering. Compared to MSM, cisgender women and heterosexual men had significantly lower education (none or any elementary: 59.1% and 54.2% vs. 16.6%, p<0.001) and socio-economic status (low income: 36.4% and 29.0% vs. 18.6%, p = 0.03) than MSM. We identified 10 (2.6%) clusters with constant growth, for prioritized intervention, that included intersecting sexual risk groups, highly connected nodes and bridge nodes between possible sub-clusters with high growth potential.

CONCLUSIONS

HIV transmission in Mexico City is strongly driven by young MSM with higher education level and recent infection. Nevertheless, leveraging network inference, we identified actively growing clusters that could be prioritized for focused intervention with demographic and risk characteristics that do not necessarily reflect the ones observed in the overall clustering population. Further studies evaluating different models to predict growing clusters are warranted. Focused interventions will have to consider structural and risk disparities between the MSM and the heterosexual populations.

Brief Report: HIV Diagnoses Among Persons Who Inject Drugs by the Urban-Rural Classification-United States, 2010-2018

BACKGROUND

After many years of decline, HIV diagnoses attributed to injection drug use in the United States increased in 2015, the year of a large outbreak among persons who inject drugs (PWIDs) in Indiana. We assessed trends in HIV diagnoses among PWID across the urban-rural continuum.

METHODS

We conducted national and county-level analyses of diagnoses among persons aged ≥13 years with HIV attributed to injection drug use only and reported to the National HIV Surveillance System through December 2019; county of residence at diagnosis was classified according to the Centers for Disease Control and Prevention's National Center for Health Statistics Urban-Rural Classification Scheme. National trends for diagnoses occurring during 2010-2014 and 2014-2018 were assessed by the estimated annual percentage change (EAPC). Counties were considered to have an "alert" (ie, an increase above baseline) if the number of 2019 diagnoses among PWID was >2 SDs and >2 diagnoses greater than the mean of annual diagnoses during 2016-2018.

RESULTS

Nationally, HIV diagnoses among PWID declined 33% during 2010-2014 from 3314 to 2220 (EAPC: -9.7%; 95% confidence interval: -10.8 to -8.6); EAPCs declined significantly in 5 of 6 urban-rural strata. During 2014-2018, diagnoses increased 11% to 2465 (EAPC: 2.4%; 95% confidence interval: 1.1 to 3.8); EAPCs were >0 for all urban-rural strata, although most were nonsignificant. Alerts were detected in 23 counties, representing 5 urban-rural strata.

CONCLUSIONS

Vigilance is needed for increases in HIV among PWID in counties across the urban-rural continuum, particularly those with indicators of increased drug use. Prompt detection, investigation, and response are critical for stemming transmission.

Interdisciplinary HIV Sentinel Case Review: Identifying Practices to Prevent Outbreaks in Philadelphia

INTRODUCTION

The Ending the HIV Epidemic in the U.S. initiative considers cluster and outbreak response essential. This article describes the design, implementation, and early findings of a Philadelphia-based project to systematically assess sentinel cases among priority populations for improving public health infrastructure and preventing future outbreaks.

METHODS

Sentinel HIV cases (i.e., early-stage or acute infection or molecular cluster cases) were identified among priority populations (Black and Hispanic/Latino men who have sex with men, youth aged 18-24 years, and transgender people who have sex with men). Chart abstraction and structured interview data were reviewed to determine themes and service gaps and to identify, prioritize, and implement recommendations. Interdisciplinary review teams included individuals with lived experience, frontline staff, and local agency leadership.

RESULTS

Data were collected during July 2019-December 2020 and analyzed for 53 of 126 sentinel cases of HIV diagnosed since July 1, 2018. The majority were men who have sex with men (79.3%), those aged 18-24 years (67.9%), and non-Hispanic Black (67.9%). More than half received sexually transmitted infection and HIV testing ≤3 years preceding HIV diagnosis (56.6% and 54.7%, respectively), had a healthcare visit within 12 months before diagnosis (64.2%), and had no evidence of pre-exposure prophylaxis awareness (58.5%). Project recommendations effectuated actions to improve pre-exposure prophylaxis provision, integrate sexually transmitted infection and HIV testing, and educate primary care providers.

CONCLUSIONS

HIV sentinel case review is a model for health departments to rapidly respond to recent transmission, identify missed HIV prevention opportunities, strengthen community partnerships, and implement programmatic and policy changes. Such efforts may prevent outbreaks and inform longer-term strategies.

Response to a Large HIV Outbreak, Cabell County, West Virginia, 2018-2019

INTRODUCTION

In January 2019, the West Virginia Bureau for Public Health detected increased HIV diagnoses among people who inject drugs in Cabell County. Responding to HIV clusters and outbreaks is 1 of the 4 pillars of the Ending the HIV Epidemic in the U.S. initiative and requires activities from the Diagnose, Treat, and Prevent pillars. This article describes the design and implementation of a comprehensive response, featuring interventions from all pillars.

METHODS

This study used West Virginia Bureau for Public Health data to identify HIV diagnoses during January 1, 2018-October 9, 2019 among (1) people who inject drugs linked to Cabell County, (2) their sex or injecting partners, or (3) others with an HIV sequence linked to Cabell County people who inject drugs. Surveillance data, including HIV-1 polymerase sequences, were analyzed to estimate the transmission rate and timing of infections using molecular clock phylogenetic analysis. Federal, state, and local partners designed and implemented a comprehensive response during January 2019-October 2019.

RESULTS

Of 82 people identified in the outbreak, most were male (60%), were White (91%), and reported unstable housing (80%). In a large molecular cluster containing 56 of 60 (93%) available sequences, 93% of inferred transmissions occurred after January 1, 2018. HIV testing, HIV pre-exposure prophylaxis, and syringe services were rapidly expanded, leading to improved linkage to HIV care and viral suppression.

CONCLUSIONS

Evidence of rapid transmission in this outbreak galvanized robust collaboration among federal, state, and local partners, leading to critical improvements in HIV prevention and care services. HIV outbreak response requires increased coordination and creativity to improve service delivery to people affected by rapid HIV transmission.

HIV Cluster and Outbreak Detection and Response: The Science and Experience

The Respond pillar of the Ending the HIV Epidemic in the U.S. initiative, which consists of activities also known as cluster and outbreak detection and response, offers a framework to guide tailored implementation of proven HIV prevention strategies where transmission is occurring most rapidly. Cluster and outbreak response involves understanding the networks in which rapid transmission is occurring; linking people in the network to essential services; and identifying and addressing gaps in programs and services such as testing, HIV and other medical care, pre-exposure prophylaxis, and syringe services programs. This article reviews the experience gained through 30 HIV cluster and outbreak responses in North America during 2000-2020 to describe approaches for implementing these core response strategies. Numerous jurisdictions that have implemented these response strategies have demonstrated success in improving outcomes related to HIV care and viral suppression, testing, use of prevention services, and reductions in transmission or new diagnoses. Efforts to address important gaps in service delivery revealed by cluster and outbreak detection and response can strengthen prevention efforts broadly through multidisciplinary, multisector collaboration. In this way, the Respond pillar embodies the collaborative, data-guided approach that is critical to the overall success of the Ending the HIV Epidemic in the U.S. initiative.

Incorporating metadata in HIV transmission network reconstruction: A machine learning feasibility assessment

HIV molecular epidemiology estimates the transmission patterns from clustering genetically similar viruses. The process involves connecting genetically similar genotyped viral sequences in the network implying epidemiological transmissions. This technique relies on genotype data which is collected only from HIV diagnosed and in-care populations and leaves many persons with HIV (PWH) who have no access to consistent care out of the tracking process. We use machine learning algorithms to learn the non-linear correlation patterns between patient metadata and transmissions between HIV-positive cases. This enables us to expand the transmission network reconstruction beyond the molecular network. We employed multiple commonly used supervised classification algorithms to analyze the San Diego Primary Infection Resource Consortium (PIRC) cohort dataset, consisting of genotypes and nearly 80 additional non-genetic features. First, we trained classification models to determine genetically unrelated individuals from related ones. Our results show that random forest and decision tree achieved over 80% in accuracy, precision, recall, and F1-score by only using a subset of meta-features including age, birth sex, sexual orientation, race, transmission category, estimated date of infection, and first viral load date besides genetic data. Additionally, both algorithms achieved approximately 80% sensitivity and specificity. The Area Under Curve (AUC) is reported 97% and 94% for random forest and decision tree classifiers respectively. Next, we extended the models to identify clusters of similar viral sequences. Support vector machine demonstrated one order of magnitude improvement in accuracy of assigning the sequences to the correct cluster compared to dummy uniform random classifier. These results confirm that metadata carries important information about the dynamics of HIV transmission as embedded in transmission clusters. Hence, novel computational approaches are needed to apply the non-trivial knowledge collected from inter-individual genetic information to metadata from PWH in order to expand the estimated transmissions. We note that feature extraction alone will not be effective in identifying patterns of transmission and will result in random clustering of the data, but its utilization in conjunction with genetic data and the right algorithm can contribute to the expansion of the reconstructed network beyond individuals with genetic data.

Characterization of HIV Risk Behaviors and Clusters Using HIV-Transmission Cluster Engine Among a Cohort of Persons Living with HIV in Washington, DC

Molecular epidemiology (ME) is one tool used to end the HIV epidemic in the United States. We combined clinical and behavioral data with HIV sequence data to identify any overlap in clusters generated from different sequence datasets; to characterize HIV transmission clusters; and to identify correlates of clustering among people living with HIV (PLWH) in Washington, District of Columbia (DC). First, Sanger sequences from DC Cohort participants, a longitudinal HIV study, were combined with next-generation sequences (NGS) from participants in a ME substudy to identify clusters. Next, demographic and self-reported behavioral data from ME substudy participants were used to identify risks of secondary transmission. Finally, we combined NGS from ME substudy participants with Sanger sequences in the DC Molecular HIV Surveillance database to identify clusters. Cluster analyses used HIV-Transmission Cluster Engine to identify linked pairs of sequences (defined as distance ≤1.5%). Twenty-eight clusters of ≥3 sequences (size range: 3-12) representing 108 (3%) participants were identified. None of the five largest clusters (size range: 5-12) included newly diagnosed PLWH. Thirty-four percent of ME substudy participants (n = 213) reported condomless sex during their last sexual encounter and 14% reported a Syphilis diagnosis in the past year. Seven transmission clusters (size range: 2-19) were identified in the final analysis, each containing at least one ME substudy participant. Substudy participants in clusters from the third analysis were present in clusters from the first analysis. Combining HIV sequence, clinical and behavioral data provided insights into HIV transmission that may not be identified using traditional epidemiological methods alone. Specifically, the sexual risk behaviors and STI diagnoses reported in the substudy survey may not have been disclosed during Partner Services activities and the survey data complemented clinical data to fully characterize transmission clusters. These findings can be used to enhance local efforts to interrupt transmission and avert new infections.

Longitudinal typing of molecular HIV clusters in a statewide epidemic

BACKGROUND

HIV molecular epidemiology is increasingly integrated into public health prevention. We conducted cluster typing to enhance characterization of a densely sampled statewide epidemic towards informing public health.

METHODS

We identified HIV clusters, categorized them into types, and evaluated their dynamics between 2004 and 2019 in Rhode Island. We grouped sequences by diagnosis year, assessed cluster changes between paired phylogenies, t0 and t1, representing adjacent years and categorized clusters as stable (cluster in t0 phylogeny = cluster in t1 phylogeny) or unstable (cluster in t0 ≠ cluster in t1). Unstable clusters were further categorized as emerging (t1 phylogeny only) or growing (larger in t1 phylogeny). We determined proportions of each cluster type, of individuals in each cluster type, and of newly diagnosed individuals in each cluster type, and assessed trends over time.

RESULTS

A total of 1727 individuals with available HIV-1 subtype B pol sequences were diagnosed in Rhode Island by 2019. Over time, stable clusters and individuals in them dominated the epidemic, increasing over time, with reciprocally decreasing unstable clusters and individuals in them. Conversely, proportions of newly diagnosed individuals in unstable clusters significantly increased. Within unstable clusters, proportions of emerging clusters and of individuals in them declined; whereas proportions of newly diagnosed individuals in growing clusters significantly increased over time.

CONCLUSION

Distinct molecular cluster types were identified in the Rhode Island epidemic. Cluster dynamics demonstrated increasing stable and decreasing unstable clusters driven by growing, rather than emerging clusters, suggesting consistent in-state transmission networks. Cluster typing could inform public health beyond conventional approaches and direct interventions.

The Role of Phylogenetics in Unravelling Patterns of HIV Transmission towards Epidemic Control: The Quebec Experience (2002-2020)

Phylogenetics has been advanced as a structural framework to infer evolving trends in the regional spread of HIV-1 and guide public health interventions. In Quebec, molecular network analyses tracked HIV transmission dynamics from 2002-2020 using MEGA10-Neighbour-joining, HIV-TRACE, and MicrobeTrace methodologies. Phylogenetics revealed three patterns of viral spread among Men having Sex with Men (MSM, n = 5024) and heterosexuals (HET, n = 1345) harbouring subtype B epidemics as well as B and non-B subtype epidemics (n = 1848) introduced through migration. Notably, half of new subtype B infections amongst MSM and HET segregating as solitary transmissions or small cluster networks (2-5 members) declined by 70% from 2006-2020, concomitant to advances in treatment-as-prevention. Nonetheless, subtype B epidemic control amongst MSM was thwarted by the ongoing genesis and expansion of super-spreader large cluster variants leading to micro-epidemics, averaging 49 members/cluster at the end of 2020. The growth of large clusters was related to forward transmission cascades of untreated early-stage infections, younger at-risk populations, more transmissible/replicative-competent strains, and changing demographics. Subtype B and non-B subtype infections introduced through recent migration now surpass the domestic epidemic amongst MSM. Phylodynamics can assist in predicting and responding to active, recurrent, and newly emergent large cluster networks, as well as the cryptic spread of HIV introduced through migration.