Phylogenetic studies of transmission dynamics in generalized HIV epidemics: an essential tool where the burden is greatest?

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

BACKGROUND

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Assessing transmission attribution risk from simulated sequencing data in HIV molecular epidemiology

BACKGROUND

HIV molecular epidemiology (ME) is the analysis of sequence data together with individual-level clinical, demographic, and behavioral data to understand HIV epidemiology. The use of ME has raised concerns regarding identification of the putative source in direct transmission events. This could result in harm ranging from stigma to criminal prosecution in some jurisdictions. Here we assessed the risks of ME using simulated HIV genetic sequencing data.

METHODS

We simulated social networks of men-who-have-sex-with-men, calibrating the simulations to data from San Diego. We used these networks to simulate consensus and next-generation sequence (NGS) data to evaluate the risks of identifying direct transmissions using different HIV sequence lengths, and population sampling depths. To identify the source of transmissions, we calculated infector probability and used phyloscanner software for the analysis of consensus and NGS data, respectively.

RESULTS

Consensus sequence analyses showed that the risk of correctly inferring the source (direct transmission) within identified transmission pairs was very small and independent of sampling depth. Alternatively, NGS analyses showed that identification of the source of a transmission was very accurate, but only for 6.5% of inferred pairs. False positive transmissions were also observed, where one or more unobserved intermediaries were present when compared to the true network.

CONCLUSION

Source attribution using consensus sequences rarely infers direct transmission pairs with high confidence but is still useful for population studies. In contrast, source attribution using NGS data was much more accurate in identifying direct transmission pairs, but for only a small percentage of transmission pairs analyzed.

Reconstructing relative transmission rates in Bayesian phylodynamics: Two-fold transmission advantage of Omicron in Berlin, Germany during December 2021

Phylodynamic methods have lately played a key role in understanding the spread of infectious diseases. During the coronavirus disease (COVID-19) pandemic, large scale genomic surveillance has further increased the potential of dynamic inference from viral genomes. With the continual emergence of novel severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) variants, explicitly allowing transmission rate differences between simultaneously circulating variants in phylodynamic inference is crucial. In this study, we present and empirically validate an extension to the BEAST2 package birth-death skyline model (BDSKY), BDSKY[Formula: see text], which introduces a scaling factor for the transmission rate between independent, jointly inferred trees. In an extensive simulation study, we show that BDSKY[Formula: see text] robustly infers the relative transmission rates under different epidemic scenarios. Using publicly available genome data of SARS-CoV-2, we apply BDSKY[Formula: see text] to quantify the transmission advantage of the Omicron over the Delta variant in Berlin, Germany. We find the overall transmission rate of Omicron to be scaled by a factor of two with pronounced variation between the individual clusters of each variant. These results quantify the transmission advantage of Omicron over the previously circulating Delta variant, in a crucial period of pre-established non-pharmaceutical interventions. By inferring variant- as well as cluster-specific transmission rate scaling factors, we show the differences in transmission dynamics for each variant. This highlights the importance of incorporating lineage-specific transmission differences in phylodynamic inference.

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

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.

Characterization of HIV-1 Transmission Clusters Inferred from the Brazilian Nationwide Genotyping Service Database

The study of HIV-1 transmission networks inferred from viral genetic data can be used to clarify important factors about the dynamics of HIV-1 transmission, such as network growth rate and demographic composition. In Brazil, HIV transmission has been stable since the early 2000s and the study of transmission clusters can provide valuable data to understand the drivers of virus spread. In this work, we analyzed a nation-wide database of approximately 53,000 HIV-1 nucleotide pol sequences sampled from genotyped patients from 2008-2017. Phylogenetic trees were reconstructed for the HIV-1 subtypes B, C and F1 in Brazil and transmission clusters were inferred by applying genetic distances thresholds of 1.5%, 3.0% and 4.5%, as well as high (>0.9) cluster statistical support. An odds ratio test revealed that young men (15-24 years) and individuals with more years of education presented higher odds to cluster. The assortativity coefficient revealed that individuals with similar demographic features tended to cluster together, with emphasis on features, such as place of residence and age. We also observed that assortativity weakens as the genetic distance threshold increases. Our results indicate that the phylogenetic clusters identified here are likely representative of the contact networks that shape HIV transmission, and this is a valuable tool even in sites with low sampling density, such as Brazil.

Using phylogenetics to infer HIV-1 transmission direction between known transmission pairs

Identifying the transmission direction between individuals provides unparalleled power to understand infectious disease epidemiology. With epidemiological and clinical information typically unavailable to infer transmission direction, phylogenetic analysis of pathogen sequence data offers an alternative approach. While the success of this phylogenetic analysis varies, the reasons remain unknown. We analyze sequence data from over 100 transmission pairs for which both the transmission direction of HIV is known and detailed additional information is available. We find that easily quantifiable phylogenetic and sampling characteristics discriminate whether a phylogenetically inferred transmission direction is correct. Our analysis highlights that while phylogenetic approaches to infer transmission direction are unsuitable for individual-level analysis, such as forensic investigations, confidence in source attribution can be incorporated in population-level analyses.

Inferring the transmission direction between linked individuals living with HIV provides unparalleled power to understand the epidemiology that determines transmission. Phylogenetic ancestral-state reconstruction approaches infer the transmission direction by identifying the individual in whom the most recent common ancestor of the virus populations originated. While these methods vary in accuracy, it is unclear why. To evaluate the performance of phylogenetic ancestral-state reconstruction to determine the transmission direction of HIV-1 infection, we inferred the transmission direction for 112 transmission pairs where transmission direction and detailed additional information were available. We then fit a statistical model to evaluate the extent to which epidemiological, sampling, genetic, and phylogenetic factors influenced the outcome of the inference. Finally, we repeated the analysis under real-life conditions with only routinely available data. We found that whether ancestral-state reconstruction correctly infers the transmission direction depends principally on the phylogeny's topology. For example, under real-life conditions, the probability of identifying the correct transmission direction increases from 32%—when a monophyletic–monophyletic or paraphyletic–polyphyletic tree topology is observed and when the tip closest to the root does not agree with the state at the root—to 93% when a paraphyletic–monophyletic topology is observed and when the tip closest to the root agrees with the root state. Our results suggest that documenting larger differences in relative intrahost diversity increases our confidence in the transmission direction inference of linked pairs for population-level studies of HIV. These findings provide a practical starting point to determine our confidence in transmission direction inference from ancestral-state reconstruction.

Deep-sequence phylogenetics to quantify patterns of HIV transmission in the context of a universal testing and treatment trial - BCPP/ Ya Tsie trial

Background: Mathematical models predict that community-wide access to HIV testing-and-treatment can rapidly and substantially reduce new HIV infections. Yet several large universal test-and-treat HIV prevention trials in high-prevalence epidemics demonstrated variable reduction in population-level incidence.

Methods: To elucidate patterns of HIV spread in universal test-and-treat trials we quantified the contribution of geographic-location, gender, age and randomized-HIV-intervention to HIV transmissions in the 30-community Ya Tsie trial in Botswana. We sequenced HIV viral whole genomes from 5,114 trial participants among the 30 trial communities.

Results: Deep-sequence phylogenetic analysis revealed that most inferred HIV transmissions within the trial occurred within the same or between neighboring communities, and between similarly-aged partners. Transmissions into intervention communities from control communities were more common than the reverse post-baseline (30% [12.2 - 56.7] versus 3% [0.1 - 27.3]) than at baseline (7% [1.5 - 25.3] versus 5% [0.9 - 22.9]) compatible with a benefit from treatment-as-prevention.

Conclusion: Our findings suggest that population mobility patterns are fundamental to HIV transmission dynamics and to the impact of HIV control strategies.

Funding: This study was supported by the National Institute of General Medical Sciences (U54GM088558); the Fogarty International Center (FIC) of the U.S. National Institutes of Health (D43 TW009610); and the President's Emergency Plan for AIDS Relief through the Centers for Disease Control and Prevention (CDC) (Cooperative agreements U01 GH000447 and U2G GH001911).

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 Discerning HIV-1 Mixing among Vulnerable Populations and Geographic Regions in Sub-Saharan Africa: A Systematic Review

To reduce global HIV-1 incidence, there is a need to understand and disentangle HIV-1 transmission dynamics and to determine the geographic areas and populations that act as hubs or drivers of HIV-1 spread. In Sub-Saharan Africa (sSA), the region with the highest HIV-1 burden, information about such transmission dynamics is sparse. Phylogenetic inference is a powerful method for the study of HIV-1 transmission networks and source attribution. In this review, we assessed available phylogenetic data on mixing between HIV-1 hotspots (geographic areas and populations with high HIV-1 incidence and prevalence) and areas or populations with lower HIV-1 burden in sSA. We searched PubMed and identified and reviewed 64 studies on HIV-1 transmission dynamics within and between risk groups and geographic locations in sSA (published 1995-2021). We describe HIV-1 transmission from both a geographic and a risk group perspective in sSA. Finally, we discuss the challenges facing phylogenetic inference in mixed epidemics in sSA and offer our perspectives and potential solutions to the identified challenges.

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

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

Transmission dynamics and molecular characterization of HIV-1 epidemic among therapeutic failure patients in Santa Catarina state, southern Brazil

The HIV-1 epidemic in southern Brazil is mostly caused by subtype C, which contrasts the dominance of subtype B in the other regions of the country. Santa Catarina (SC), although the smallest state in the southern region, presents one of the highest incidences and mortality rates in Brazil due to AIDS. This work investigated the HIV-1 molecular diversity and phylogenetic transmission networks in SC state by analyzing a database of 3070 sequences of the national genotyping service. HIV-1C proved to be the most frequent subtype, with a significant increase in prevalence over time. HIV-1B was observed to be associated with highly educated men, suggesting a compartmentalization from other subtypes. Such observation was confirmed by the high frequency of HIV-1B circulating in MSM transmission networks. Identified transmission clusters were majority composed by individuals living up to 25 km away and interstate linkages were mainly between southern neighbor states. In general, individuals between 25 and 40 years old and sequences sampled after 2014 were more likely to be in transmission chains, in agreement with the universal treatment protocol launched in 2014. The present study brings new insights about HIV-1 transmission dynamics in southern Brazil.

Evaluation of Phylogenetic Methods for Inferring the Direction of Human Immunodeficiency Virus (HIV) Transmission: HIV Prevention Trials Network (HPTN) 052

BACKGROUND

Phylogenetic analysis can be used to assess human immunodeficiency virus (HIV) transmission in populations. We inferred the direction of HIV transmission using whole-genome HIV sequences from couples with known linked infection and known transmission direction.

METHODS

Complete next-generation sequencing (NGS) data were obtained for 105 unique index-partner sample pairs from 32 couples enrolled in the HIV Prevention Trials Network (HPTN) 052 study (up to 2 samples/person). Index samples were obtained up to 5.5 years before partner infection; partner samples were obtained near the time of seroconversion. The bioinformatics method, phyloscanner, was used to infer transmission direction. Analyses were performed using samples from individual sample pairs, samples from all couples (1 sample/person; group analysis), and all available samples (multisample group analysis). Analysis was also performed using NGS data from defined regions of the HIV genome (gag, pol, env).

RESULTS

Using whole-genome NGS data, transmission direction was inferred correctly (index to partner) for 98 of 105 (93.3%) of the individual sample pairs, 99 of 105 (94.3%) sample pairs using group analysis, and 31 of the 32 couples (96.9%) using multisample group analysis. There were no cases where the incorrect transmission direction (partner to index) was inferred. The accuracy of the method was higher with greater time between index and partner sample collection. Pol region sequences performed better than env or gag sequences for inferring transmission direction.

CONCLUSIONS

We demonstrate the potential of a phylogenetic method to infer the direction of HIV transmission between 2 individuals using whole-genome and pol NGS data.

Identification of Hidden Population Structure in Time-Scaled Phylogenies

Population structure influences genealogical patterns, however, data pertaining to how populations are structured are often unavailable or not directly observable. Inference of population structure is highly important in molecular epidemiology where pathogen phylogenetics is increasingly used to infer transmission patterns and detect outbreaks. Discrepancies between observed and idealized genealogies, such as those generated by the coalescent process, can be quantified, and where significant differences occur, may reveal the action of natural selection, host population structure, or other demographic and epidemiological heterogeneities. We have developed a fast non-parametric statistical test for detection of cryptic population structure in time-scaled phylogenetic trees. The test is based on contrasting estimated phylogenies with the theoretically expected phylodynamic ordering of common ancestors in two clades within a coalescent framework. These statistical tests have also motivated the development of algorithms which can be used to quickly screen a phylogenetic tree for clades which are likely to share a distinct demographic or epidemiological history. Epidemiological applications include identification of outbreaks in vulnerable host populations or rapid expansion of genotypes with a fitness advantage. To demonstrate the utility of these methods for outbreak detection, we applied the new methods to large phylogenies reconstructed from thousands of HIV-1 partial pol sequences. This revealed the presence of clades which had grown rapidly in the recent past and was significantly concentrated in young men, suggesting recent and rapid transmission in that group. Furthermore, to demonstrate the utility of these methods for the study of antimicrobial resistance, we applied the new methods to a large phylogeny reconstructed from whole genome Neisseria gonorrhoeae sequences. We find that population structure detected using these methods closely overlaps with the appearance and expansion of mutations conferring antimicrobial resistance. [Antimicrobial resistance; coalescent; HIV; population structure.].

Factors Associated With Human Immunodeficiency Virus Infections Linked in Genetic Clusters But Disconnected in Partner Tracing

BACKGROUND

Successful partner notification can improve community-level outcomes by increasing the proportion of persons living with human immunodeficiency virus (HIV) who are linked to HIV care and virally suppressed, but it is resource intensive. Understanding where HIV transmission pathways may be undetected by routine partner notification may help improve case finding strategies.

METHODS

We combined partner notification interview and HIV sequence data for persons diagnosed with HIV in Wake County, NC in 2012 to 2013 to evaluate partner contact networks among persons with HIV pol gene sequences 2% or less pairwise genetic distance. We applied a set of multivariable generalized estimating equations to identify correlates of disparate membership in genetic versus partner contact networks.

RESULTS

In the multivariable model, being in a male-male pair (adjusted odds ratio [AOR], 16.7; P = 0.01), chronic HIV infection status (AOR, 4.5; P < 0.01), and increasing percent genetic distance between each dyad member's HIV pol gene sequence (AOR, 8.3 per each 1% increase, P < 0.01) were all associated with persons with HIV clustering but not being identified in the partner notification network component. Having anonymous partners or other factors typically associated with risk behavior were not associated.

CONCLUSIONS

Based on genetic networks, partnerships which may be stigmatized, may have occurred farther back in time or may have an intervening partner were more likely to be unobserved in the partner contact network. The HIV genetic cluster information contributes to public health understanding of HIV transmission networks in these settings where partner identifying information is not available.

Persistence of HIV transmission clusters among people who inject drugs

OBJECTIVE

We investigated the duration of HIV transmission clusters.

DESIGN

Fifty-four individuals newly infected at enrollment in the ALIVE cohort were included, all of whom had sequences at an intake visit (T1) and from a second (T2) and/or a third (T3) follow-up visit, median 2.9 and 5.4 years later, respectively.

METHODS

Sequences were generated using the 454 DNA sequencing platform for portions of HIV pol and env (HXB2 positions 2717-3230; 7941-8264). Genetic distances were calculated using tn93 and sequences were clustered over a range of thresholds (1--5%) using HIV-TRACE. Analyses were performed separately for individuals with pol sequences for T1 + T2 (n = 40, 'Set 1') and T1 + T3 (n = 25; 'Set 2'), and env sequences for T1 + T2 (n = 47, 'Set 1'), and T1 + T3 (n = 30; 'Set 2').

RESULTS

For pol, with one exception, a single cluster contained more than 75% of samples at all thresholds, and cluster composition was at least 90% concordant between time points/thresholds. For env, two major clusters (A and B) were observed at T1 and T2/T3, although cluster composition concordance between time points/thresholds was low (<60%) at lower thresholds for both sets 1 and 2. In addition, several individuals were included in clusters at T2/T3, although not at T1.

CONCLUSION

Caution should be used in applying a single threshold in population studies where seroconversion dates are unknown. However, the retention of some clusters even after 5 + years is evidence for the robustness of the clustering approach in general.

Brief Report: Temporal Changes in HIV Transmission Patterns Among Young Men Who Have Sex With Men, United States, 2009-2016

BACKGROUND

In the United States, young men (aged 13-24 years) who have sex with men (MSM) bear a disproportionate burden of HIV. Transmission among MSM has been found to be disassortative by age.

METHODS

We analyzed HIV-1 pol sequences reported to the US National HIV Surveillance System from MSM with HIV diagnosed during 2009-2016. Using an HIV genetic transmission network, we identified persons with closely related viruses (ie, genetic distance ≤1.5%) and used multivariable logistic regression to examine changes from 2009-2012 to 2013-2016 in proportions of MSM linked to young MSM who were >5 years older or of the same race/ethnicity.

RESULTS

Among 9510 young MSM linked to another MSM with a closely related virus, 37% linked to an older MSM and 62% linked to an MSM of the same race/ethnicity. Comparing 2013-2016 with 2009-2012, we found increases in linkage of older MSM to young MSM, with the most substantial increases seen in Hispanic/Latinos aged 13-19 [adjusted prevalence ratio (APR) = 1.31, 95% confidence interval (CI) = 1.11 to 1.56] and blacks aged 13-19 (APR = 1.23, CI = 1.06 to 1.41) and 20-24 years (APR = 1.14, CI = 1.02 to 1.28). By contrast, change in linkage patterns among racial/ethnic groups was unremarkable.

CONCLUSIONS

We found evidence of increased age mixing among MSM with respect to HIV transmission over time, which coincides temporally with changes in partner-seeking behavior such as increased use of mobile applications. These findings indicate the importance of social factors on HIV sexual and transmission networks and suggest that prevention efforts need to effectively reach MSM of all ages.

A Comprehensive Genomics Solution for HIV Surveillance and Clinical Monitoring in Low-Income Settings

Viral genetic sequencing can be used to monitor the spread of HIV drug resistance, identify appropriate antiretroviral regimes, and characterize transmission dynamics. Despite decreasing costs, next-generation sequencing (NGS) is still prohibitively costly for routine use in generalized HIV epidemics in low- and middle-income countries. Here, we present veSEQ-HIV, a high-throughput, cost-effective NGS sequencing method and computational pipeline tailored specifically to HIV, which can be performed using leftover blood drawn for routine CD4 cell count testing. This method overcomes several major technical challenges that have prevented HIV sequencing from being used routinely in public health efforts; it is fast, robust, and cost-efficient, and generates full genomic sequences of diverse strains of HIV without bias. The complete veSEQ-HIV pipeline provides viral load estimates and quantitative summaries of drug resistance mutations; it also exploits information on within-host viral diversity to construct directed transmission networks. We evaluated the method's performance using 1,620 plasma samples collected from individuals attending 10 large urban clinics in Zambia as part of the HPTN 071-2 study (PopART Phylogenetics). Whole HIV genomes were recovered from 91% of samples with a viral load of >1,000 copies/ml. The cost of the assay (30 GBP per sample) compares favorably with existing VL and HIV genotyping tests, proving an affordable option for combining HIV clinical monitoring with molecular epidemiology and drug resistance surveillance in low-income settings.

Human Immunodeficiency Virus (HIV)-1 Transmission Among Persons With Acute HIV-1 Infection in Malawi: Demographic, Behavioral, and Phylogenetic Relationships

BACKGROUND

Understanding sexual networks involving acute human immunodeficiency virus (HIV)-1 infections (AHI) may lead to prevention opportunities to mitigate high rates of onward transmission. We evaluated HIV-1 phylogenetic and behavioral characteristics among persons with AHI and their referred partners.

METHODS

Between 2012 and 2014, 46 persons with AHI in Malawi participated in a combined behavioral and biomedical intervention. Participants referred sexual partners by passive referral. Demographics and sexual behaviors were collected through interviews and HIV-1 genetic relationships were assessed with phylogenetics.

RESULTS

Among 45 AHI participants with HIV-1 sequences, none was phylogenetically-linked with another AHI index. There were 19 (42%) AHI participants who referred a single partner that returned for testing. Most partners (n = 17) were HIV-infected, with 15 (88%) presenting with an established infection. There were 14 index-partner pairs that had sequences available; 13 (93%) pairs were phylogenetically-linked dyads. The AHI index was female in 7/13 (54%) dyads. Age-disparate relationships among dyads were common (≥5-year age difference in 67% of dyads), including 3/6 dyads involving a male index and a younger woman. Index participants with a referred partner were more likely to report no casual partners and to be living with their current partner than participants not in dyads.

CONCLUSIONS

Passive-partner referral successfully identified partners with genetically-similar HIV infections-the likely source of infection-but only 40% of index cases referred partners who presented for HIV-1 testing. Future work evaluating assisted partner notification may help reach susceptible partners or more people with untreated HIV-1 infections connected to acute transmission.

CLINICAL TRIALS REGISTRATION

NCT01450189.

Informed consent for HIV phylogenetic research: A case study of urban individuals living with HIV approached for enrollment in an HIV study

INTRODUCTION

Phylogenetic analyses can provide insights on HIV tansmission dynamics. Country and state-level differences in HIV criminalization and disclosure laws and advances in next generation sequencing could impact perceived study risks.

METHODS

We present study opt-out rates and the reasons provided during enrollment for a study conducted in Boston (6/2017-8/2018).

RESULTS

Of 90 patients approached to participate, 45 did not consent to participate. Reasons for not participating included unwillingness to discuss their HIV status, privacy and confidentiality concerns, disinterest, and lack of time.

CONLUSIONS

Given low participation rates and concerns related to HIV disclosure, privacy, and confidentiality, these questions remain (1) should informed consent be required for all phylogenetic analyses, including deidentified and surveillance data? (2) what additional steps can researchers take to protect the privacy of individuals, particularly in contexts where HIV is criminalized or there have been civil/criminal cases investigating HIV transmission? And (3) what role can community members play to minimize the potential risks, particularly for those most marginalized? These questions require input from both researchers and community members living with HIV/AIDS.

Migration interacts with the local transmission of HIV in developed trade areas: A molecular transmission network analysis in China

The HIV-1 epidemic is a remarkable public health concern in China, especially in developed trade areas. We aimed to investigate the interaction of migration with the local transmission network in a typical trade area, Yiwu City, the world's largest commodity distribution center. Based on 390 pol sequences from 413 participants diagnosed between 2014 and 2016, putative transmission clusters and the underlying demographic and behavioral characteristics were analyzed. Recent infection status was determined by HIV-1 limiting antigen avidity enzyme immunoassay to identify active clusters. Multiple subtypes were identified, with a predominance of CRF01_AE (47.4%) and CRF07_BC (40.8%), followed by 9 other subtypes and 8 URFs. Multivariable analyses revealed that individuals in clusters were more likely to be local residents, infected through heterosexual behaviors, and infected with CRF01_AE (P < .05). Of men who have sex with men (MSM), 81% were linked to other MSM, and only 3% were linked to heterosexual women. Of heterosexual women, 67% were linked to heterosexual men, and 11% to MSM. Yiwu residents were more likely to link to locals than that of migrants (43% vs 20%, P < .001). By contrast, local MSM and migrant MSM all had high percentages of linkage to migrant MSM (57% vs 69%, P = .069). Our findings reveal that migration promotes the dissemination and dynamic change of HIV, which are interwoven between locals and migrants. The results highlight the far-reaching influence of migrant MSM on the local HIV transmission network.

Characteristics Associated With Human Immunodeficiency Virus Transmission Networks Involving Adolescent Girls and Young Women in Human Immunodeficiency Virus Prevention Trials Network 068 Study

We combined behavioral survey data from the human immunodeficiency virus (HIV) Prevention Trials Network 068 study with phylogenetic information to determine if cluster membership was associated with characteristics of young women and their partners. Clusters were more likely to involve young women from specific villages and schools, indicating some localized transmission.Supplemental digital content is available in the text.

Understanding disclosed and cryptic HIV transmission risk via genetic analysis: what are we missing and when does it matter?

PURPOSE OF REVIEW

To discuss the recent HIV phylogenetic analyses examining HIV transmission patterns among and within risk groups.

RECENT FINDINGS

Phylodynamic analysis has recently been applied to multiple HIV outbreaks among people who inject drugs to determine whether HIV transmission is ongoing. Large-scale analyses of datasets of HIV sequences collected for drug-resistance testing provide population-level insights into transmission patterns. One focus across world regions has been to investigate whether age-disparity is a driver of HIV transmission. In sub-Saharan Africa, researchers have examined transmission between heterosexuals and MSM and between high prevalence fishing communities and inland communities. In the US and the UK, cryptic risk groups such as nondisclosed MSM and the partners of transgender women are increasingly being uncovered based on their position in densely sampled molecular transmission networks.

SUMMARY

Analysis of HIV genetic sequence can resolve viral transmission patterns between risk groups at unprecedented scales and levels of detail. Future research should focus on understanding the effect of missing data on inferences and the biases of different methods. Uncovering groups and patterns obscured from traditional epidemiolocal analyses is exciting but should not compromise the privacy of the groups in question.

Evaluation of A Phylogenetic Pipeline to Examine Transmission Networks in A Canadian HIV Cohort

Keywords: HIV; Canada; molecular phylogenetics; viral evolution; person-to-person transmission inference; transmission network; summary statistics.

Phylogenetic Clustering by Linear Integer Programming (PhyCLIP)

Subspecies nomenclature systems of pathogens are increasingly based on sequence data. The use of phylogenetics to identify and differentiate between clusters of genetically similar pathogens is particularly prevalent in virology from the nomenclature of human papillomaviruses to highly pathogenic avian influenza (HPAI) H5Nx viruses. These nomenclature systems rely on absolute genetic distance thresholds to define the maximum genetic divergence tolerated between viruses designated as closely related. However, the phylogenetic clustering methods used in these nomenclature systems are limited by the arbitrariness of setting intra and intercluster diversity thresholds. The lack of a consensus ground truth to define well-delineated, meaningful phylogenetic subpopulations amplifies the difficulties in identifying an informative distance threshold. Consequently, phylogenetic clustering often becomes an exploratory, ad hoc exercise. Phylogenetic Clustering by Linear Integer Programming (PhyCLIP) was developed to provide a statistically principled phylogenetic clustering framework that negates the need for an arbitrarily defined distance threshold. Using the pairwise patristic distance distributions of an input phylogeny, PhyCLIP parameterizes the intra and intercluster divergence limits as statistical bounds in an integer linear programming model which is subsequently optimized to cluster as many sequences as possible. When applied to the hemagglutinin phylogeny of HPAI H5Nx viruses, PhyCLIP was not only able to recapitulate the current WHO/OIE/FAO H5 nomenclature system but also further delineated informative higher resolution clusters that capture geographically distinct subpopulations of viruses. PhyCLIP is pathogen-agnostic and can be generalized to a wide variety of research questions concerning the identification of biologically informative clusters in pathogen phylogenies. PhyCLIP is freely available at http://github.com/alvinxhan/PhyCLIP, last accessed March 15, 2019.

Distinct rates and patterns of spread of the major HIV-1 subtypes in Central and East Africa

Since the ignition of the HIV-1 group M pandemic in the beginning of the 20th century, group M lineages have spread heterogeneously throughout the world. Subtype C spread rapidly through sub-Saharan Africa and is currently the dominant HIV lineage worldwide. Yet the epidemiological and evolutionary circumstances that contributed to its epidemiological expansion remain poorly understood. Here, we analyse 346 novel pol sequences from the DRC to compare the evolutionary dynamics of the main HIV-1 lineages, subtypes A1, C and D. Our results place the origins of subtype C in the 1950s in Mbuji-Mayi, the mining city of southern DRC, while subtypes A1 and D emerged in the capital city of Kinshasa, and subtypes H and J in the less accessible port city of Matadi. Following a 15-year period of local transmission in southern DRC, we find that subtype C spread at least three-fold faster than other subtypes circulating in Central and East Africa. In conclusion, our results shed light on the origins of HIV-1 main lineages and suggest that socio-historical rather than evolutionary factors may have determined the epidemiological fate of subtype C in sub-Saharan Africa.

Exploring HIV-1 Transmission Dynamics by Combining Phylogenetic Analysis and Infection Timing

HIV-1 pol sequences obtained through baseline drug resistance testing of patients newly diagnosed between 2013 and 2017 were analyzed for genetic similarity. For 927 patients the information on genetic similarity was combined with demographic data and with information on the recency of infection. Overall, 48.3% of the patients were genetically linked with 11.4% belonging to a pair and 36.9% involved in a cluster of ≥3 members. The percentage of early diagnosed (≤4 months after infection) was 28.6%. Patients of Belgian origin were more frequently involved in transmission clusters (49.7% compared to 15.3%) and diagnosed earlier (37.4% compared to 12.2%) than patients of Sub-Saharan African origin. Of the infections reported to be locally acquired, 69.5% were linked (14.1% paired and 55.4% in a cluster). Equal parts of early and late diagnosed individuals (59.9% and 52.4%, respectively) were involved in clusters. The identification of a genetically linked individual for the majority of locally infected patients suggests a high rate of diagnosis in this population. Diagnosis however is often delayed for >4 months after infection increasing the opportunities for onward transmission. Prevention of local infection should focus on earlier diagnosis and protection of the still uninfected members of sexual networks with human immunodeficiency virus (HIV)-infected members.

Ethical issues associated with HIV molecular epidemiology: a qualitative exploratory study using inductive analytic approaches

Background

HIV molecular epidemiology is increasingly recognized as a vital source of information for understanding HIV transmission dynamics. Despite extensive use of these data-intensive techniques in both research and public health settings, the ethical issues associated with this science have received minimal attention. As the discipline evolves, there is reasonable concern that existing ethical and legal frameworks and standards might lag behind the rapid methodological developments in this field. This is a follow-up on our earlier work that applied a predetermined analytical framework to examine the perspectives of a sample of scientists from the fields of epidemiology, public health, virology and bioethics on key ethical issues associated with HIV molecular epidemiology in HIV network research.

Methods

Fourteen in-depth interviews were conducted with scientists from the fields of molecular epidemiology, public health, virology and bioethics. Inductive analytical approaches were applied to identify key themes that emerged from the data.

Results

Our interviewees acknowledged the potential positive impact of molecular epidemiology in the fight against HIV. However, they were concerned that HIV phylogenetics research messages may be incorrectly interpreted if not presented at the appropriate level. There was consensus that HIV phylogenetics research presents a potential risk to privacy, but the probability and magnitude of this risk was less obvious. Although participants acknowledged the social value that could be realized from the analysis of HIV genetic sequences, there was a perceived fear that the boundaries for use of HIV sequence data were not clearly defined.

Conclusions

Our findings highlight distinct ethical issues arising from HIV molecular epidemiology. As the discipline evolves and HIV sequence data become increasingly available, it is critical to ensure that ethical standards keep pace with biomedical advancements. We argue that the ethical issues raised in this study, whether real or perceived, require further conceptual and empirical examination.

Leveraging Phylogenetics to Understand HIV Transmission and Partner Notification Networks

BACKGROUND

Partner notification is an important component of public health test and treat interventions. To enhance this essential function, we assessed the potential for molecular methods to supplement routine partner notification and corroborate HIV networks.

METHODS

All persons diagnosed with HIV infection in Wake County, NC, during 2012-2013 and their disclosed sexual partners were included in a sexual network. A data set containing HIV-1 pol sequences collected in NC during 1997-2014 from 15,246 persons was matched to HIV-positive persons in the network and used to identify putative transmission clusters. Both networks were compared.

RESULTS

The partner notification network comprised 280 index cases and 383 sexual partners and high-risk social contacts (n = 131 HIV-positive). Of the 411 HIV-positive persons in the partner notification network, 181 (44%) did not match to a HIV sequence, 61 (15%) had sequences but were not identified in a transmission cluster, and 169 (41%) were identified in a transmission cluster. More than half (59%) of transmission clusters bridged sexual network partnerships that were not recognized in the partner notification; most of these clusters were dominated by men who have sex with men.

CONCLUSIONS

Partner notification and HIV sequence analysis provide complementary representations of the existent partnerships underlying the HIV transmission network. The partner notification network components were bridged by transmission clusters, particularly among components dominated by men who have sex with men. Supplementing the partner notification network with phylogenetic data highlighted avenues for intervention.

A Short-Term Assessment of Nascent HIV-1 Transmission Clusters Among Newly Diagnosed Individuals Using Envelope Sequence-Based Phylogenetic Analyses

The identification of transmission clusters (TCs) of HIV-1 using phylogenetic analyses can provide insights into viral transmission network and help improve prevention strategies. We compared the use of partial HIV-1 envelope fragment of 1,070 bp with its loop 3 (108 bp) to determine its utility in inferring HIV-1 transmission clustering. Serum samples of recently (n = 106) and chronically (n = 156) HIV-1-infected patients with status confirmed were sequenced. HIV-1 envelope nucleotide-based phylogenetic analyses were used to infer HIV-1 TCs. Those were constructed using ClusterPickerGUI_1.2.3 considering a pairwise genetic distance of ≤10% threshold. Logistic regression analyses were used to examine the relationship between the demographic factors that were likely associated with HIV-1 clustering. Ninety-eight distinct consensus envelope sequences were subjected to phylogenetic analyses. Using a partial envelope fragment sequence, 42 sequences were grouped into 15 distinct small TCs while the V3 loop reproduces 10 clusters. The agreement between the partial envelope and the V3 loop fragments was significantly moderate with a Cohen's kappa (κ) coefficient of 0.59, p < .00001. The mean age (<38.8 years) and HIV-1 B subtype are two factors identified that were significantly associated with HIV-1 transmission clustering in the cohort, odds ratio (OR) = 0.25, 95% confidence interval (CI, 0.04–0.66), p = .002 and OR: 0.17, 95% CI (0.10–0.61), p = .011, respectively. The present study confirms that a partial fragment of the HIV-1 envelope sequence is a better predictor of transmission clustering. However, the loop 3 segment may be useful in screening purposes and may be more amenable to integration in surveillance programs.

Molecular Epidemiology of HIV-1 Subtype B Reveals Heterogeneous Transmission Risk: Implications for Intervention and Control

Background

The impact of HIV pre-exposure prophylaxis (PrEP) depends on infections averted by protecting vulnerable individuals as well as infections averted by preventing transmission by those who would have been infected if not receiving PrEP. Analysis of HIV phylogenies reveals risk factors for transmission, which we examine as potential criteria for allocating PrEP.

Methods

We analyzed 6912 HIV-1 partial pol sequences from men who have sex with men (MSM) in the United Kingdom combined with global reference sequences and patient-level metadata. Population genetic models were developed that adjust for stage of infection, global migration of HIV lineages, and changing incidence of infection through time. Models were extended to simulate the effects of providing susceptible MSM with PrEP.

Results

We found that young age <25 years confers higher risk of HIV transmission (relative risk = 2.52 [95% confidence interval, 2.32–2.73]) and that young MSM are more likely to transmit to one another than expected by chance. Simulated interventions indicate that 4-fold more infections can be averted over 5 years by focusing PrEP on young MSM.

Conclusions

Concentrating PrEP doses on young individuals can avert more infections than random allocation.

Enhanced use of phylogenetic data to inform public health approaches to HIV among men who have sex with men

The multidimensional nature and continued evolution of HIV epidemics among men who have sex with men (MSM) requires innovative intervention approaches. Strategies are needed that recognise the individual, social and structural factors driving HIV transmission; that can pinpoint networks with heightened transmission risk; and that can help target intervention in real time. HIV phylogenetics is a rapidly evolving field with strong promise for informing innovative responses to the HIV epidemic among MSM. Currently, HIV phylogenetic insights are providing new understandings of characteristics of HIV epidemics involving MSM, social networks influencing transmission, characteristics of HIV transmission clusters involving MSM, targets for antiretroviral and other prevention strategies and dynamics of emergent epidemics. Maximising the potential of HIV phylogenetics for HIV responses among MSM will require attention to key methodological challenges and ethical considerations, as well as resolving key implementation and scientific questions. Enhanced and integrated use of HIV surveillance, sociobehavioural and phylogenetic data resources are becoming increasingly critical for informing public health approaches to HIV among MSM.

HIV-TRACE (TRAnsmission Cluster Engine): a Tool for Large Scale Molecular Epidemiology of HIV-1 and Other Rapidly Evolving Pathogens

In modern applications of molecular epidemiology, genetic sequence data are routinely used to identify clusters of transmission in rapidly evolving pathogens, most notably HIV-1. Traditional 'shoe-leather' epidemiology infers transmission clusters by tracing chains of partners sharing epidemiological connections (e.g., sexual contact). Here, we present a computational tool for identifying a molecular transmission analog of such clusters: HIV-TRACE (TRAnsmission Cluster Engine). HIV-TRACE implements an approach inspired by traditional epidemiology, by identifying chains of partners whose viral genetic relatedness imply direct or indirect epidemiological connections. Molecular transmission clusters are constructed using codon-aware pairwise alignment to a reference sequence followed by pairwise genetic distance estimation among all sequences. This approach is computationally tractable and is capable of identifying HIV-1 transmission clusters in large surveillance databases comprising tens or hundreds of thousands of sequences in near real time, that is, on the order of minutes to hours. HIV-TRACE is available at www.hivtrace.org and from www.github.com/veg/hivtrace, along with the accompanying result visualization module from www.github.com/veg/hivtrace-viz. Importantly, the approach underlying HIV-TRACE is not limited to the study of HIV-1 and can be applied to study outbreaks and epidemics of other rapidly evolving pathogens.

Genetic Cluster Analysis for HIV Prevention

PURPOSE OF REVIEW

This review summarizes the use of genetic similarity clusters to understand HIV transmission and inform prevention efforts.

RECENT FINDINGS

Recent emphases include the development of real-time cluster identification in order to interrupt transmission chains, the use of clusters to estimate rates of transmission along the HIV care cascade, and the extension of cluster analyses to understand transmission in the generalized epidemics of sub-Saharan Africa. Importantly, this recent empirical work has been accompanied by theoretical work that elucidates the processes that underlie HIV genetic similarity clusters; multiple studies suggest that clusters are not necessarily enriched with individuals with high transmission rates, but rather can reflect variation in sampling times within a population, with individuals sampled early in infection more likely to cluster. Analyses of genetic similarity clusters have great promise to inform HIV epidemiology and prevention. Future emphases should include the collection of additional sequence data from underrepresented populations, such as those in sub-Saharan Africa, and further development and evaluation of clustering methods.

Using Contact Patterns to Inform HIV Interventions in Persons Who Inject Drugs in Northern Vietnam

OBJECTIVES

Population mixing patterns can greatly inform allocation of HIV prevention interventions such as treatment as prevention or preexposure prophylaxis. Characterizing contact patterns among subgroups can help identify the specific combinations of contact expected to result in the greatest number of new infections.

SETTING

Baseline data from an intervention to reduce HIV-related risk behaviors in male persons who inject drugs (PWID) in the Northern Vietnamese province of Thai Nguyen were used for the analysis.

METHODS

Egocentric network data were provided by PWID who reported any drug-injection equipment sharing in the previous 3 months. Age-dependent mixing was assessed to explore its epidemiological implications on risk of HIV transmission risk (among those HIV-infected) and HIV acquisition risk (among those not infected) in PWID.

RESULTS

A total of 1139 PWID collectively reported 2070 equipment-sharing partnerships in the previous 3 months. Mixing by age identified the 30-34 and 35-39 years age groups as the groups from whom the largest number of new infections was transmitted, making them primary targets for treatment as prevention. Among the uninfected, 25-29, 30-35, and 35-39 years age groups had the highest HIV acquisition rate, making them the primary targets for preexposure prophylaxis.

CONCLUSIONS

Collection and analysis of contact patterns in PWID is feasible and can greatly inform infectious disease dynamics and targeting of appropriate interventions. Results presented also provide much needed empirical data on mixing to improve mathematical models of disease transmission in this population.

Estimating Effects of HIV Sequencing Data Completeness on Transmission Network Patterns and Detection of Growing HIV Transmission Clusters

HIV nucleotide sequence data can identify clusters of persons with genetically similar strains suggesting transmission. We simulated the effect of lowered data completeness, defined by the percent of persons with diagnosed HIV with a reported sequence, on transmission patterns and detection of growing HIV transmission clusters. We analyzed HIV surveillance data for persons with HIV diagnosed during 2008-2014 who resided in Michigan or Washington. We calculated genetic distances, constructed the inferred transmission network for each jurisdiction, and compared transmission network characteristics and detection of growing transmission clusters in the full dataset with artificially reduced datasets. Simulating lower levels of completeness resulted in decreased percentages of persons linked to a cluster from high completeness (full dataset) to low completeness (5%) (Michigan: 54%-18%; Washington, 46%-16%). Patterns of transmission between certain populations remained robust as data completeness level was reduced. As data completeness was artificially decreased, sensitivity of cluster detection substantially diminished in both states. In Michigan, sensitivity decreased from 100% with the full dataset, to 62% at 50% completeness and 21% at 25% completeness. In Washington, sensitivity decreased from 100% with the full dataset, to 71% at 50% completeness and 29% at 25% completeness. Lower sequence data completeness limits the ability to detect clusters that may benefit from investigation; however, inferences can be made about transmission patterns even with low data completeness, given sufficient numbers. Data completeness should be prioritized, as lack of or delays in detection of transmission clusters could result in additional infections.

Inferring HIV-1 transmission networks and sources of epidemic spread in Africa with deep-sequence phylogenetic analysis

To prevent new infections with human immunodeficiency virus type 1 (HIV-1) in sub-Saharan Africa, UNAIDS recommends targeting interventions to populations that are at high risk of acquiring and passing on the virus. Yet it is often unclear who and where these 'source' populations are. Here we demonstrate how viral deep-sequencing can be used to reconstruct HIV-1 transmission networks and to infer the direction of transmission in these networks. We are able to deep-sequence virus from a large population-based sample of infected individuals in Rakai District, Uganda, reconstruct partial transmission networks, and infer the direction of transmission within them at an estimated error rate of 16.3% [8.8-28.3%]. With this error rate, deep-sequence phylogenetics cannot be used against individuals in legal contexts, but is sufficiently low for population-level inferences into the sources of epidemic spread. The technique presents new opportunities for characterizing source populations and for targeting of HIV-1 prevention interventions in Africa.

Increasing importance of European lineages in seeding the hepatitis C virus subtype 1a epidemic in Spain

BackgroundReducing the burden of the hepatitis C virus (HCV) requires large-scale deployment of intervention programmes, which can be informed by the dynamic pattern of HCV spread. In Spain, ongoing transmission of HCV is mostly fuelled by people who inject drugs (PWID) infected with subtype 1a (HCV1a).AimOur aim was to map how infections spread within and between populations, which could help formulate more effective intervention programmes to halt the HCV1a epidemic in Spain.MethodsEpidemiological links between HCV1a viruses from a convenience sample of 283 patients in Spain, mostly PWID, collected between 2014 and 2016, and 1,317, 1,291 and 1,009 samples collected abroad between 1989 and 2016 were reconstructed using sequences covering the NS3, NS5A and NS5B genes. To efficiently do so, fast maximum likelihood-based tree estimation was coupled to a flexible Bayesian discrete phylogeographic inference method.ResultsThe transmission network structure of the Spanish HCV1a epidemic was shaped by continuous seeding of HCV1a into Spain, almost exclusively from North America and European countries. The latter became increasingly relevant and have dominated in recent times. Export from Spain to other countries in Europe was also strongly supported, although Spain was a net sink for European HCV1a lineages. Spatial reconstructions showed that the epidemic in Spain is diffuse, without large, dominant within-country networks.ConclusionTo boost the effectiveness of local intervention efforts, concerted supra-national strategies to control HCV1a transmission are needed, with a strong focus on the most important drivers of ongoing transmission, i.e. PWID and other high-risk populations.

Insights on transmission of HIV from phylogenetic analysis to locally optimize HIV prevention strategies

PURPOSE OF REVIEW

Phylogenetic analysis can identify transmission networks by clustering genetically related HIV genotypes that are routinely collected. In this study, we will review phylogenetic insights gained on transmission of HIV and phylogenetically optimized HIV prevention strategies.

RECENT FINDINGS

Phylogenetic analysis reports that HIV transmission varies by geographical region and by route of transmission. In high-income countries, HIV is predominantly transmitted between recently infected MSM who live in the same country. In rural Uganda, transmission of HIV is frequently between different communities. Age-discrepant transmission has been reported across the world. Four studies have used phylogenetic optimization of HIV prevention. Three studies predict that immediate treatment after diagnosis would have prevented 19-42% of infections, and that preexposure prophylaxis would have prevented 66% of infections. One phylogenetic study guided a public health response to an actively ongoing HIV outbreak. Phylogenetic clustering requires a dense sample of patients and small time-gaps between infection and diagnosis.

SUMMARY

Phylogenetic analysis can be an important tool to identify a local strategy that prevents most infections. Future studies that use phylogenetic analysis for optimizing HIV prevention strategies should also include cost-effectiveness so that the most cost-effective prevention method is identified.

Experts' Perspectives on Key Ethical Issues Associated With HIV Phylogenetics as Applied in HIV Transmission Dynamics Research

The use of phylogenetics in HIV molecular epidemiology has considerably increased our ability to understand the origin, spread, and characteristics of HIV epidemics. Despite its potential to advance knowledge on HIV transmission dynamics, the ethical issues associated with HIV molecular epidemiology have received minimal attention. In-depth interviews were conducted with scientists from diverse backgrounds to explore their perspectives on ethical issues associated with phylogenetic analysis of HIV genetic data as applied to HIV transmission dynamics studies. The Emanuel framework was used as the analytical framework. Favorable risk-benefit ratio and informed consent were the most invoked ethical principles and fair participant selection the least. Fear of loss of privacy and disclosure of HIV transmission were invariably cited as key ethical concerns. As HIV sequence data become increasingly available, comprehensive guidelines should be developed to guide its access, sharing and use, cognizant of the potential harms that may result.

Growth of HIV-1 Molecular Transmission Clusters in New York City

Background

HIV-1 genetic sequences can be used to infer viral transmission history and dynamics. Throughout the United States, HIV-1 sequences from drug resistance testing are reported to local public health departments.

Methods

We investigated whether inferred HIV transmission network dynamics can identify individuals and clusters of individuals most likely to give rise to future HIV cases in a surveillance setting. We used HIV-TRACE, a genetic distance-based clustering tool, to infer molecular transmission clusters from HIV-1 pro/RT sequences from 65736 people in the New York City surveillance registry. Logistic and LASSO regression analyses were used to identify correlates of clustering and cluster growth, respectively. We performed retrospective transmission network analyses to evaluate individual- and cluster-level prioritization schemes for identifying parts of the network most likely to give rise to new cases in the subsequent year.

Results

Individual-level prioritization schemes predicted network growth better than random targeting. Across the 3600 inferred molecular transmission clusters, previous growth dynamics were superior predictors of future transmission cluster growth compared to individual-level prediction schemes. Cluster-level prioritization schemes considering previous cluster growth relative to cluster size further improved network growth predictions.

Conclusions

Prevention efforts based on HIV molecular epidemiology may improve public health outcomes in a US surveillance setting.

Ethical considerations in global HIV phylogenetic research

Phylogenetic analysis of pathogens is an increasingly powerful way to reduce the spread of epidemics, including HIV. As a result, phylogenetic approaches are becoming embedded in public health and research programmes, as well as outbreak responses, presenting unique ethical, legal, and social issues that are not adequately addressed by existing bioethics literature. We formed a multidisciplinary working group to explore the ethical issues arising from the design of, conduct in, and use of results from HIV phylogenetic studies, and to propose recommendations to minimise the associated risks to both individuals and groups. We identified eight key ethical domains, within which we highlighted factors that make HIV phylogenetic research unique. In this Review, we endeavoured to provide a framework to assist researchers, public health practitioners, and funding institutions to ensure that HIV phylogenetic studies are designed, done, and disseminated in an ethical manner. Our conclusions also have broader relevance for pathogen phylogenetics.

Transmission dynamics among participants initiating antiretroviral therapy upon diagnosis of early acute HIV-1 infection in Thailand

OBJECTIVE

To assess transmission characteristics in a predominantly MSM cohort initiating antiretroviral therapy (ART) immediately following diagnosis of acute HIV-1infection (AHI).

METHODS

A longitudinal study (2009-2017) was performed in participants with AHI (n = 439) attending a single clinic in Bangkok. Plasma samples obtained prior to ART were used to obtain HIV-1 pol sequences and combined with clinical and epidemiologic data to assess transmission dynamics (cluster formation and size) using phylogenetic analysis. Clusters were estimated using maximum likelihood, genetic distance of 1.5% and visual inspection. The potential transmitter(s) in a cluster was determined using time to viral suppression and interview data.

RESULTS

The cohort was predominantly MSM (93%) and infected with HIV-1 CRF01_AE (87%). Medians (ranges) for age and viral load prior to ART were 26 (18-70) years and 5.9 (2.5-8.2) log10 HIV-1 RNA copies/ml. Median time from history of HIV-1 exposure to diagnosis was 19 (3-61) days. Viral suppression was observed in 388 of 412 (94%) participants at a median time of 12 weeks following ART. Twenty-six clusters with median cluster size of 2 (2-5) representing 62 of 439 (14%) participants were observed. Younger age was associated with cluster formation: median 28 versus 30 years for unique infections (P = 0.01). A potential transmitter was identified in 11 of 26 (42%) clusters.

CONCLUSION

Despite high rates of viral suppression following diagnosis and treatment of AHI within a cohort of young Thai MSM, HIV-1 transmission continued, reflecting the need to expand awareness and treatment access to the entire MSM population.

Introduction of Phylodynamics for Controlling the HIV/AIDS Epidemic in Korea

As over 1000 new cases of HIV/AIDS occur in Korea annually, preventive health programs against HIV/AIDS are urgently needed. Since phylodynamic studies have been suggested as a way to understand how infectious diseases are transmitted and evolve, phylodynamic inferences can be a useful tool for HIV/AIDS research. In particular, phylodynamic models are helpful for dating the origins of an epidemic and estimating its basic reproduction number. Thus, the introduction of phylodynamics would be a highly valuable step towards controlling the HIV/AIDS epidemic in Korea.

HIV-1 diversity among young women in rural South Africa: HPTN 068

Background

South Africa has one of the highest rates of HIV-1 (HIV) infection world-wide, with the highest rates among young women. We analyzed the molecular epidemiology and evolutionary history of HIV in young women attending high school in rural South Africa.

Methods

Samples were obtained from the HPTN 068 randomized controlled trial, which evaluated the effect of cash transfers for school attendance on HIV incidence in women aged 13–20 years (Mpumalanga province, 2011–2015). Plasma samples from HIV-infected participants were analyzed using the ViroSeq HIV-1 Genotyping assay. Phylogenetic analysis was performed using 200 pol gene study sequences and 2,294 subtype C reference sequences from South Africa. Transmission clusters were identified using Cluster Picker and HIV-TRACE, and were characterized using demographic and other epidemiological data. Phylodynamic analyses were performed using the BEAST software.

Results

The study enrolled 2,533 young women who were followed through their expected high school graduation date (main study); some participants had a post-study assessment (follow-up study). Two-hundred-twelve of 2,533 enrolled young women had HIV infection. HIV pol sequences were obtained for 94% (n = 201/212) of the HIV-infected participants. All but one of the sequences were HIV-1 subtype C; the non-C subtype sequence was excluded from further analysis. Median pairwise genetic distance between the subtype C sequences was 6.4% (IQR: 5.6–7.2). Overall, 26% of study sequences fell into 21 phylogenetic clusters with 2–6 women per cluster. Thirteen (62%) clusters included women who were HIV-infected at enrollment. Clustering was not associated with study arm, demographic or other epidemiological factors. The estimated date of origin of HIV subtype C in the study population was 1958 (95% highest posterior density [HPD]: 1931–1980), and the median estimated substitution rate among study pol sequences was 1.98x10^-3 (95% HPD: 1.15x10^-3–2.81x10^-3) per site per year.

Conclusions

Phylogenetic analysis suggests that multiple HIV subtype C sublineages circulate among school age girls in South Africa. There were no substantive differences in the molecular epidemiology of HIV between control and intervention arms in the HPTN 068 trial.

Tracking external introductions of HIV using phylodynamics reveals a major source of infections in rural KwaZulu-Natal, South Africa

Despite increasing access to antiretrovirals, HIV incidence in rural KwaZulu-Natal remains among the highest ever reported in Africa. While many epidemiological factors have been invoked to explain such high incidence, widespread human mobility and viral movement suggest that transmission between communities may be a major source of new infections. High cross-community transmission rates call into question how effective increasing the coverage of antiretroviral therapy locally will be at preventing new infections, especially if many new cases arise from external introductions. To help address this question, we use a phylodynamic model to reconstruct epidemic dynamics and estimate the relative contribution of local transmission versus external introductions to overall incidence in KwaZulu-Natal from HIV-1 phylogenies. By comparing our results with population-based surveillance data, we show that we can reliably estimate incidence from viral phylogenies once viral movement in and out of the local population is accounted for. Our analysis reveals that early epidemic dynamics were largely driven by external introductions. More recently, we estimate that 35 per cent (95% confidence interval: 20-60%) of new infections arise from external introductions. These results highlight the growing need to consider larger-scale regional transmission dynamics when designing and testing prevention strategies.

Molecular evolution methods to study HIV-1 epidemics

Nucleotide sequences of HIV isolates are obtained routinely to evaluate the presence of resistance mutations to antiretroviral drugs. But, beyond their clinical use, these and other viral sequences include a wealth of information that can be used to better understand and characterize the epidemiology of HIV in relevant populations. In this review, we provide a brief overview of the main methods used to analyze HIV sequences, the data bases where reference sequences can be obtained, and some caveats about the possible applications for public health of these analyses, along with some considerations about their limitations and correct usage to derive robust and reliable conclusions.

Factors Driving the HIV Epidemic in Southern Africa

The HIV pandemic has disproportionately impacted sub-Saharan Africa and Southern Africa in particular. The concurrent presence of overlapping epidemic drivers likely underpins how and why the HIV epidemic is so explosive in this region, with implications for understanding approaches to reduce transmission. In this review, we discuss the relative contribution and interaction between epidemic drivers in the Southern African context, including factors both distally and proximally associated with the likelihood and degree of exposure to HIV and factors that increase the probability of transmission when exposure occurs. In particular, we focus on young women as a key population in need of HIV prevention and highlight factors that increase their risk on several levels.