Public health interventions guided by clustering of HIV-1 molecular sequences may be impacted by choices of analytical approaches. We identified commonly-used clustering analytical approaches, applied them to 1886 HIV-1 Rhode Island sequences from 2004–2018, and compared concordance in identifying molecular HIV-1 clusters within and between approaches. We used strict (topological support ≥ 0.95; distance 0.015 substitutions/site) and relaxed (topological support 0.80–0.95; distance 0.030–0.045 substitutions/site) thresholds to reflect different epidemiological scenarios. We found that clustering differed by method and threshold and depended more on distance than topological support thresholds. Clustering concordance analyses demonstrated some differences across analytical approaches, with RAxML having the highest (91%) mean summary percent concordance when strict thresholds were applied, and three (RAxML-, FastTree regular bootstrap- and IQ-Tree regular bootstrap-based) analytical approaches having the highest (86%) mean summary percent concordance when relaxed thresholds were applied. We conclude that different analytical approaches can yield diverse HIV-1 clustering outcomes and may need to be differentially used in diverse public health scenarios. Recognizing the variability and limitations of commonly-used methods in cluster identification is important for guiding clustering-triggered interventions to disrupt new transmissions and end the HIV 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, t ₀ and t ₁, representing adjacent years and categorized clusters as stable (cluster in t ₀ phylogeny = cluster in t ₁ phylogeny) or unstable (cluster in t ₀ ≠ cluster in t ₁). Unstable clusters were further categorized as emerging (t ₁ phylogeny only) or growing (larger in t ₁ 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.

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.

Molecular HIV cluster data can guide public health responses towards ending the HIV epidemic. Currently, real-time data integration, analysis, and interpretation are challenging, leading to a delayed public health response. We present a comprehensive methodology for addressing these challenges through data integration, analysis, and reporting. We integrated heterogeneous data sources across systems and developed an open-source, automatic bioinformatics pipeline that provides molecular HIV cluster data to inform public health responses to new statewide HIV-1 diagnoses, overcoming data management, computational, and analytical challenges. We demonstrate implementation of this pipeline in a statewide HIV epidemic and use it to compare the impact of specific phylogenetic and distance-only methods and datasets on molecular HIV cluster analyses. The pipeline was applied to 18 monthly datasets generated between January 2020 and June 2022 in Rhode Island, USA, that provide statewide molecular HIV data to support routine public health case management by a multi-disciplinary team. The resulting cluster analyses and near-real-time reporting guided public health actions in 37 phylogenetically clustered cases out of 57 new HIV-1 diagnoses. Of the 37, only 21 (57%) clustered by distance-only methods. Through a unique academic-public health partnership, an automated open-source pipeline was developed and applied to prospective, routine analysis of statewide molecular HIV data in near-real-time. This collaboration informed public health actions to optimize disruption of HIV transmission.

Drug resistance remains a global challenge in children and adolescents living with HIV (CALWH). Characterizing resistance evolution, specifically using next generation sequencing (NGS) can potentially inform care, but remains understudied, particularly in antiretroviral therapy (ART)-experienced CALWH in resource-limited settings. We conducted reverse-transcriptase NGS and investigated short-and long-term resistance evolution and its predicted impact in a well-characterized cohort of Kenyan CALWH failing 1st-line ART and followed for up to ~8 years. Drug resistance mutation (DRM) evolution types were determined by NGS frequency changes over time, defined as evolving (up-trending and crossing the 20% NGS threshold), reverting (down-trending and crossing the 20% threshold) or other. Exploratory analyses assessed potential impacts of minority resistance variants on evolution. Evolution was detected in 93% of 42 participants, including 91% of 22 with short-term follow-up, 100% of 7 with long-term follow-up without regimen change, and 95% of 19 with long-term follow-up with regimen change. Evolving DRMs were identified in 60% and minority resistance variants evolved in 17%, with exploratory analysis suggesting greater rate of evolution of minority resistance variants under drug selection pressure and higher predicted drug resistance scores in the presence of minority DRMs. Despite high-level pre-existing resistance, NGS-based longitudinal follow-up of this small but unique cohort of Kenyan CALWH demonstrated continued DRM evolution, at times including low-level DRMs detected only by NGS, with predicted impact on care. NGS can inform better understanding of DRM evolution and dynamics and possibly improve care. The clinical significance of these findings should be further evaluated.

Background: Characterizing HIV clustering rates and their trends over time can improve understanding a local epidemic and enhance its control. Methods: Leveraging an academic–public health partnership in Rhode Island, we explored longitudinal dynamics of statewide clustering rates among key populations from 1991 to 2023. Partial HIV-1 pol sequences were grouped by year of HIV-1 diagnosis. Molecular clusters were identified in cumulative annual phylogenies. Overall clustering rates, and clustering rates of newly diagnosed and prevalent infections, and of specific sociodemographic characteristics of key populations over time were determined. Mann–Kendall statistics were used to estimate clustering rate trends and relationships among groups. Results: By the end of 2023, 2630 individuals with sequences represented the statewide epidemic in Rhode Island. Overall clustering rates increased from 7% in 1991 to 46% in 2023, correlating with cumulative sequence increase. Clustering rates of newly diagnosed and prevalent infections significantly increased over time, higher in newly diagnosed individuals since the early 2000s. Increases were also observed among groups defined by gender, age, transmission risks, race, mental illness, HIV-1 subtypes, and country of birth, with some crossovers and divergence patterns over time. Conclusion: Exploring dynamics of HIV clustering rates over three decades in a statewide HIV-1 epidemic expanded its characterization and provided insight into its evolving changes. These dynamics may indicate a gradual shift towards a more concentrated and localized HIV-1 epidemic, highlighting important opportunities for targeted interventions to effectively prevent new HIV transmissions.