Detailed molecular epidemiologic characterization of HIV-1 infection in Bulgaria reveals broad diversity and evolving phylodynamics

Transmitted HIV Drug Resistance in Bulgaria Occurs in Clusters of Individuals from Different Transmission Groups and Various Subtypes (2012-2020)

Transmitted HIV drug resistance in Bulgaria was first reported in 2015 using data from 1988-2011. We determined the prevalence of surveillance drug resistance mutations (SDRMs) and HIV-1 genetic diversity in Bulgaria during 2012-2020 using polymerase sequences from 1053 of 2010 (52.4%) antiretroviral therapy (ART)-naive individuals. Sequences were analyzed for DRM using the WHO HIV SDRM list implemented in the calculated population resistance tool at Stanford University. Genetic diversity was inferred using automated subtyping tools and phylogenetics. Cluster detection and characterization was performed using MicrobeTrace. The overall rate of SDRMs was 5.7% (60/1053), with 2.2% having resistance to nucleoside reverse transcriptase inhibitors (NRTIs), 1.8% to non-nucleoside reverse transcriptase inhibitors (NNRTIs), 2.1% to protease inhibitors (PIs), and 0.4% with dual-class SDRMs. We found high HIV-1 diversity, with the majority being subtype B (60.4%), followed by F1 (6.9%), CRF02_AG (5.2%), A1 (3.7%), CRF12_BF (0.8%), and other subtypes and recombinant forms (23%). Most (34/60, 56.7%) of the SDRMs were present in transmission clusters of different subtypes composed mostly of male-to-male sexual contact (MMSC), including a 14-member cluster of subtype B sequences from 12 MMSC and two males reporting heterosexual contact; 13 had the L90M PI mutation and one had the T215S NRTI SDRM. We found a low SDRM prevalence amid high HIV-1 diversity among ART-naive patients in Bulgaria during 2012-2020. The majority of SDRMs were found in transmission clusters containing MMSC, indicative of onward spread of SDRM in drug-naive individuals. Our study provides valuable information on the transmission dynamics of HIV drug resistance in the context of high genetic diversity in Bulgaria, for the development of enhanced prevention strategies to end the epidemic.

Molecular Epidemiological Analysis of the Origin and Transmission Dynamics of the HIV-1 CRF01_AE Sub-Epidemic in Bulgaria

HIV-1 subtype CRF01_AE is the second most predominant strain in Bulgaria, yet little is known about the molecular epidemiology of its origin and transmissibility. We used a phylodynamics approach to better understand this sub-epidemic by analyzing 270 HIV-1 polymerase (pol) sequences collected from persons diagnosed with HIV/AIDS between 1995 and 2019. Using network analyses at a 1.5% genetic distance threshold (d), we found a large 154-member outbreak cluster composed mostly of persons who inject drugs (PWID) that were predominantly men. At d = 0.5%, which was used to identify more recent transmission, the large cluster dissociated into three clusters of 18, 12, and 7 members, respectively, five dyads, and 107 singletons. Phylogenetic analysis of the Bulgarian sequences with publicly available global sequences showed that CRF01_AE likely originated from multiple Asian countries, with Vietnam as the likely source of the outbreak cluster between 1988 and 1990. Our findings indicate that CRF01_AE was introduced into Bulgaria multiple times since 1988, and infections then rapidly spread among PWID locally with bridging to other risk groups and countries. CRF01_AE continues to spread in Bulgaria as evidenced by the more recent large clusters identified at d = 0.5%, highlighting the importance of public health prevention efforts in the PWID communities.

Molecular Epidemiology of the HIV-1 Subtype B Sub-Epidemic in Bulgaria

HIV-1 subtype B is the predominant strain in Bulgaria, yet little is known about the molecular epidemiology of these infections, including its origin and transmissibility. We used a phylodynamics approach by combining and analyzing 663 HIV-1 polymerase (pol) sequences collected from persons diagnosed with HIV/AIDS between 1988-2018 and associated epidemiologic data to better understand this sub-epidemic in Bulgaria. Using network analyses at a 1.5% genetic distance threshold (d) we found several large phylogenetic clusters composed mostly of men who have sex with men (MSM) and male heterosexuals (HET). However, at d = 0.5%, used to identify more recent transmission, the largest clusters dissociated to become smaller in size. The majority of female HET and persons with other transmission risks were singletons or pairs in the network. Phylogenetic analysis of the Bulgarian pol sequences with publicly available global sequences showed that subtype B was likely introduced into Bulgaria from multiple countries, including Israel and several European countries. Our findings indicate that subtype B was introduced into Bulgaria multiple times since 1988 and then infections rapidly spread among MSM and non-disclosed MSM. These high-risk behaviors continue to spread subtype B infection in Bulgaria as evidenced by the large clusters at d = 0.5%. Relatively low levels of antiretroviral drug resistance were observed in our study. Prevention strategies should continue to include increased testing and linkage to care and treatment, as well as expanded outreach to the MSM communities.

The prevalence of hepatitis B and C co-infections among people with HIV-1 in Bulgaria: 2010–2015

Aim: To evaluate hepatitis B virus (HBV) and hepatitis C virus (HCV) among individuals with HIV/AIDS in Bulgaria diagnosed between 2010 and 2015. Materials & methods: A total of 1158 individuals were diagnosed with HIV/AIDS during the study period. Different transmission groups were tested with ELISA and real-time PCR for HBV and HCV markers. Results: Hepatitis B surface antigen and hepatitis C virus antiboby were found in 9.3 and 23.2% of the tested. HBV DNA and HCV RNA has been found in 47.4 and 69.6%. Hepatitis B and C co-infections were predominant in multiple risk behavior groups, including people who inject drugs, men who have sex with men, prisoners and Roma individuals. Conclusion: HIV prevalence in Bulgaria is low but the rates of hepatitis B and C co-infections among these patients fall within the upper range reported in Europe.

HIV-1 genetic diversity and demographic characteristics in Bulgaria

HIV-1 strain diversity in Bulgaria is extensive and includes contributions from nearly all major subtypes and the Circulating Recombinant Forms (CRF): 01_AE, 02_AG, and 05_DF. Prior to this study, HIV-1 sequence information from Bulgaria has been based solely on the pro-RT gene, which represent less than 15% of the viral genome. To further characterize HIV-1 in Bulgaria, assess participant risk behaviors, and strengthen knowledge of circulating strains in the region, the study “Genetic Subtypes of HIV-1 in Bulgaria (RV240)” was conducted. This study employed the real time-PCR based Multi-region Hybridization Assay (MHA) B/non-B and HIV-1 sequencing to survey 215 of the approximately 1100 known HIV-1 infected Bulgarian adults (2008–2009) and determine if they were infected with subtype B HIV-1. The results indicated a subtype B prevalence of 40% and demonstrate the application of the MHA B/non-B in an area containing broad HIV-1 strain diversity. Within the assessed risk behaviors, the proportion of subtype B infection was greatest in men who have sex with men and lowest among those with drug use risk factors. During this study, 15 near full-length genomes and 22 envelope sequences were isolated from study participants. Phylogenetic analysis shows the presence of subtypes A1, B, C, F1, and G, CRF01_AE, CRF02_AG, CRF05_DF, and one unique recombinant form (URF). These sequences also show the presence of two strain groups containing participants with similar risk factors. Previous studies in African and Asian cohorts have shown that co-circulation of multiple subtypes can lead to viral recombination within super-infected individuals and the emergence of new URFs. The low prevalence of URFs in the presence of high subtype diversity in this study, may be the result of successful infection prevention and control programs. Continued epidemiological monitoring and support of infection prevention programs will help maintain control of the HIV-1 epidemic in Bulgaria.

Earlier Initiation of Antiretroviral Treatment Coincides With an Initial Control of the HIV-1 Sub-Subtype F1 Outbreak Among Men-Having-Sex-With-Men in Flanders, Belgium

Human immunodeficiency virus type 1 (HIV-1) non-B subtype infections occurred in Belgium since the 1980s, mainly amongst migrants and heterosexuals, whereas subtype B predominated in men-having-sex-with-men (MSM). In the last decade, the diagnosis of F1 sub-subtype in particular has increased substantially, which prompted us to perform a detailed reconstruction of its epidemiological history. To this purpose, the Belgian AIDS Reference Laboratories collected HIV-1 pol sequences from all sub-subtype F1-infected patients for whom genotypic drug resistance testing was requested as part of routine clinical follow-up. This data was complemented with HIV-1 pol sequences from countries with a high burden of F1 infections or a potential role in the global origin of sub-subtype F1. The molecular epidemiology of the Belgian subtype F1 epidemic was investigated using Bayesian phylogenetic inference and transmission dynamics were characterized based on birth-death models. F1 sequences were retained from 297 patients diagnosed and linked to care in Belgium between 1988 and 2015. Phylogenetic inference indicated that among the 297 Belgian F1 sequences, 191 belonged to a monophyletic group that mainly contained sequences from people likely infected in Belgium (OR 26.67, 95% CI 9.59–74.15), diagnosed in Flanders (OR 7.28, 95% CI 4.23–12.53), diagnosed at a recent stage of infection (OR 7.19, 95% CI 2.88-17.95) or declared to be MSM (OR 34.8, 95% CI 16.0–75.6). Together with a Spanish clade, this Belgian clade was embedded in the genetic diversity of Brazilian subtype F1 strains and most probably emerged after one or only a few migration events from Brazil to the European continent before 2002. The origin of the Belgian outbreak was dated back to 2002 (95% higher posterior density 2000–2004) and birth-death models suggested that its extensive growth had been controlled (R_e < 1) by 2012, coinciding with a time period where delay in antiretroviral treatment initiation substantially declined. In conclusion, phylogenetic reconstruction of the Belgian HIV-1 sub-subtype F1 epidemic illustrates the introduction and substantial dissemination of viral strains in a geographically restricted risk group that was most likely controlled by effective treatment as prevention.

A 28-Year History of HIV-1 Drug Resistance and Transmission in Washington, DC

Washington, DC consistently has one of the highest annual rates of new HIV-1 diagnoses in the United States over the last 10 years. To guide intervention and prevention strategies to combat DC HIV infection, it is helpful to understand HIV transmission dynamics in a historical context. Toward this aim, we conducted a retrospective study (years 1987–2015) of 3,349 HIV pol sequences (1,026 bp) from 1,995 individuals living in the DC area belonging to three different cohorts. We coupled HIV sequence data with clinical information (sex, risk factor, race/ethnicity, viral load, subtype, anti-retroviral regimen) to identify circulating drug resistant mutations (DRM) and transmission clusters and assess their persistence over time. Of the transmission clusters identified in the DC area, 78.0 and 31.7% involved MSM and heterosexuals, respectively. The longest spread of time for a single cluster was 5 years (2007–2012) using a distance-based network inference approach and 27 years (1987–2014) using a maximum likelihood phylogenetic approach. We found eight subtypes and nine recombinants. Genetic diversity increased steadily over time with a slight peak in 2009 and remained constant thereafter until 2015. Nucleotide diversity also increased over time while relative genetic diversity (BEAST) remained relatively steady over the last 28 years with slight increases since 2000 in subtypes B and C. Sequences from individuals on drug therapy contained the highest total number of DRMs (1,104–1,600) and unique DRMs (63–97) and the highest proportion (>20%) of resistant individuals. Heterosexuals (43.94%), MSM (40.13%), and unknown (44.26%) risk factors showed similar prevalence of DRMs, while injection drug users had a lower prevalence (33.33%). Finally, there was a 60% spike in the number of codons with DRMs between 2007 and 2010. Past patterns of HIV transmission and DRM accumulation over time described here will help to predict future efficacy of ART drugs based on DRMs persisting over time and identify risk groups of interest for prevention and intervention efforts within the DC population. Our results show how longitudinal data can help to understand the temporal dynamics of HIV-1 at the local level.

Origin and Spread of HIV-1 Subtype B Among Heterosexual Individuals in Bulgaria

Human immunodeficiency virus (HIV) was originally introduced in Bulgaria through heterosexual transmission (HET) and later transferred to other vulnerable groups along with numerous more recent introductions from outside Bulgaria. To define the diversity, origins, and dynamics of the HIV-1 subtypes prevalent in HET population in Bulgaria, we applied phylogenetic and phylodynamic analyses using polymerase (pol) sequences from HET individuals to infer the spatiotemporal evolutionary history of the HIV-1 epidemic in this population in Bulgaria. High genetic diversity was found, including 13 different HIV-1 subtypes: 45.7% subtype B, 19.9% CRF01_AE, 7.5% CRF02_AG, 7.5% sub-subtypes A1 and A6, 7.1% subtype C, 5.3% subtype F1, 4.0% URFs, 1.2% CRF05_DF, 0.6% subtype G, 0.3% CRF04_cpx, 0.3% CRF29_BF, 0.3% CRF14_BG, and 0.3% subtype H. The estimated root of the subtype B in the phylogenetic tree dated back to the year 1980 largely due to multiple introductions of subtype B from outside the country. Several significant clades have been identified highlighting six different main epidemic entrances of subtype B dating from 1989 to 2007. The Bayesian skyline plot showed two different exponential growth periods starting in the 1980s to 1990 followed by a constant phase up to about 2008, with another exponential growth period from 2008 to the year 2012. The migration analysis identified dynamic pattern of gene flow and demonstrated that many HET probably acquired the infection abroad (14.6%), while only (6.6%) of non-HET were infected outside country. The phylogenetic analysis showed an intermixing between sequences from Bulgarians with sequences from other countries, suggesting different HIV introduction in this country followed by the internal spread through local transmission networks.

A national study of the molecular epidemiology of HIV-1 in Australia 2005-2012

Introduction

Rates of new HIV-1 diagnoses are increasing in Australia, with evidence of an increasing proportion of non-B HIV-1 subtypes reflecting a growing impact of migration and travel. The present study aims to define HIV-1 subtype diversity patterns and investigate possible HIV-1 transmission networks within Australia.

Methods

The Australian Molecular Epidemiology Network (AMEN) HIV collaborating sites in Western Australia, South Australia, Victoria, Queensland and western Sydney (New South Wales), provided baseline HIV-1 partial pol sequence, age and gender information for 4,873 patients who had genotypes performed during 2005–2012. HIV-1 phylogenetic analyses utilised MEGA V6, with a stringent classification of transmission pairs or clusters (bootstrap ≥98%, genetic distance ≤1.5% from at least one other sequence in the cluster).

Results

HIV-1 subtype B represented 74.5% of the 4,873 sequences (WA 59%, SA 68.4%, w-Syd 73.8%, Vic 75.6%, Qld 82.1%), with similar proportion of transmission pairs and clusters found in the B and non-B cohorts (23% vs 24.5% of sequences, p = 0.3). Significantly more subtype B clusters were comprised of ≥3 sequences compared with non-B clusters (45.0% vs 24.0%, p = 0.021) and significantly more subtype B pairs and clusters were male-only (88% compared to 53% CRF01_AE and 17% subtype C clusters). Factors associated with being in a cluster of any size included; being sequenced in a more recent time period (p<0.001), being younger (p<0.001), being male (p = 0.023) and having a B subtype (p = 0.02). Being in a larger cluster (>3) was associated with being sequenced in a more recent time period (p = 0.05) and being male (p = 0.008).

Conclusion

This nationwide HIV-1 study of 4,873 patient sequences highlights the increased diversity of HIV-1 subtypes within the Australian epidemic, as well as differences in transmission networks associated with these HIV-1 subtypes. These findings provide epidemiological insights not readily available using standard surveillance methods and can inform the development of effective public health strategies in the current paradigm of HIV prevention in Australia.

High Rate of Hepatitis B and C Coinfections Among People Living with HIV-1 in Bulgaria: 2010-2014

In a representative nationwide study, we have determined the prevalence of hepatitis B virus (HBV) and hepatitis C virus (HCV) coinfections among HIV-positive patients diagnosed during the period 2010-2014 in Bulgaria. Despite a relatively low rate of new HIV diagnoses, the rates of hepatitis B and C coinfections among these patients fell within the upper range reported in Europe. HBsAg and HCV antibodies (Ab) were found in 10.4% and 25.6% of the tested individuals, respectively. Importantly, high rates of active hepatitis infections were confirmed by detection of HBV DNA in 51.1% and HCV RNA in 78.1% of the tested individuals. Hepatitis coinfections affected mostly high risk groups and persons with multiple risk behavior, including people who inject drugs, men who have sex with men, prisoners, and Roma people.

The global spread of HIV-1 subtype B epidemic

Human immunodeficiency virus type 1 (HIV-1) was discovered in the early 1980s when the virus had already established a pandemic. For at least three decades the epidemic in the Western World has been dominated by subtype B infections, as part of a sub-epidemic that traveled from Africa through Haiti to United States. However, the pattern of the subsequent spread still remains poorly understood. Here we analyze a large dataset of globally representative HIV-1 subtype B strains to map their spread around the world over the last 50years and describe significant spread patterns. We show that subtype B travelled from North America to Western Europe in different occasions, while Central/Eastern Europe remained isolated for the most part of the early epidemic. Looking with more detail in European countries we see that the United Kingdom, France and Switzerland exchanged viral isolates with non-European countries than with European ones. The observed pattern is likely to mirror geopolitical landmarks in the post-World War II era, namely the rise and the fall of the Iron Curtain and the European colonialism. In conclusion, HIV-1 spread through specific migration routes which are consistent with geopolitical factors that affected human activities during the last 50years, such as migration, tourism and trade. Our findings support the argument that epidemic control policies should be global and incorporate political and socioeconomic factors.

Molecular epidemiology of HIV-1 infection in Europe: An overview

Human Immunodeficiency Virus type 1 (HIV-1) is characterised by vast genetic diversity. Globally circulating HIV-1 viruses are classified into distinct phylogenetic strains (subtypes, sub-subtypes) and several recombinant forms. Here we describe the characteristics and evolution of European HIV-1 epidemic over time through a review of published literature and updated queries of existing HIV-1 sequence databases. HIV-1 in Western and Central Europe was introduced in the early-1980s in the form of subtype B, which is still the predominant clade. However, in Eastern Europe (Former Soviet Union (FSU) countries and Russia) the predominant strain, introduced into Ukraine in the mid-1990s, is subtype A (A_FSU) with transmission mostly occurring in People Who Inject Drugs (PWID). In recent years, the epidemic is evolving towards a complex tapestry with an increase in the prevalence of non-B subtypes and recombinants in Western and Central Europe. Non-B epidemics are mainly associated with immigrants, heterosexuals and females but more recently, non-B clades have also spread amongst groups where non-B strains were previously absent - non-immigrant European populations and amongst men having sex with men (MSM). In some countries, non-B clades have spread amongst the native population, for example subtype G in Portugal and subtype A in Greece, Albania and Cyprus. Romania provides a unique case where sub-subtype F1 has predominated throughout the epidemic. In contrast, HIV-1 epidemic in FSU countries remains more homogeneous with A_FSU clade predominating in all countries. The differences between the evolution of the Western epidemic and the Eastern epidemic may be attributable to differences in transmission risk behaviours, lifestyle and the patterns of human mobility. The study of HIV-1 epidemic diversity provides a useful tool by which we can understand the history of the pandemic in addition to allowing us to monitor the spread and growth of the epidemic over time.

Phylodynamics of the major HIV-1 CRF02_AG African lineages and its global dissemination

The HIV-1 CRF02_AG clade is the most prevalent HIV variant in West and West-Central Africa and its detection outside Africa is increasingly common. Little is known, however, about the number and phylodynamics of major CRF02_AG lineages circulating worldwide. To this end, a total of 3170 HIV-1 CRF02_AG-like pol sequences isolated around the world, over a period of 25years (1989 to 2013), were analyzed using Maximum Likelihood and Bayesian coalescent-based methods. Our results suggest that most of the current CRF02_AG diversity comes from the dissemination of a few founder strains out of Central Africa into West Africa and Cameroon between the late 1960s and the middle 1980s. The CRF02_AG strain introduced into West Africa established a large regional epidemic with low phylogeographic structure. This strain was also successfully disseminated out of the West African region and originated at least three large secondary outbreaks in Cameroon at around the late 1970s, in the former Soviet Union (FSU) countries at around the late 1990s, and in Bulgaria/Germany at around the early 2000s. The CRF02_AG African lineages introduced into Cameroon remained mostly restricted to this country and its neighbors. Demographic reconstructions indicate that major CRF02_AG clades circulating in Africa exhibited a decline in growth rate since the middle 1980s/1990s, whereas CRF02_AG clades in Europe and the FSU countries continue to grow exponentially until the middle to late 2000s. Substantial differences in the median estimated growth rate of the same CRF02_AG clade circulating in different regions (0.63-2.00year^-1), and of different CRF02_AG clades circulating in the same country (0.41-0.75year^-1) were observed. Thus, the cause of the epidemic outcome of the different HIV-1 CRF02_AG lineages is probably multifactorial.

Occupational HIV Transmission Among Male Adult Film Performers - Multiple States, 2014

In 2014, the California Department of Public Health was notified by a local health department of a diagnosis of acute human immunodeficiency virus (HIV) infection* and rectal gonorrhea in a male adult film industry performer, aged 25 years (patient A). Patient A had a 6-day history of rash, fever, and sore throat suggestive of acute retroviral syndrome at the time of examination. He was informed of his positive HIV and gonorrhea test results 6 days after his examination. Patient A had a negative HIV-1 RNA qualitative nucleic acid amplification test (NAAT)(†) 10 days before symptom onset. This investigation found that during the 22 days between the negative NAAT and being informed of his positive HIV test results, two different production companies directed patient A to have condomless sex with a total of 12 male performers. Patient A also provided contact information for five male non-work-related sexual partners during the month before and after his symptom onset. Patient A had additional partners during this time period for which no locating information was provided. Neither patient A nor any of his interviewed sexual partners reported taking HIV preexposure prophylaxis (PrEP). Contact tracing and phylogenetic analysis of HIV sequences amplified from pretreatment plasma revealed that a non-work-related partner likely infected patient A, and that patient A likely subsequently infected both a coworker during the second film production and a non-work-related partner during the interval between his negative test and receipt of his positive HIV results. Adult film performers and production companies, medical providers, and all persons at risk for HIV should be aware that testing alone is not sufficient to prevent HIV transmission. Condom use provides additional protection from HIV and sexually transmitted infections (STIs). Performers and all persons at risk for HIV infection in their professional and personal lives should discuss the use of PrEP with their medical providers.

Low HIV-1 transmitted drug resistance in Bulgaria against a background of high clade diversity

OBJECTIVES

To determine transmitted drug resistance (TDR) and HIV-1 genetic diversity in Bulgaria.

METHODS

The prevalence of TDR and HIV-1 subtypes was determined in 305/1446 (21.1%) persons newly diagnosed with HIV/AIDS from 1988 to 2011. TDR mutations (TDRMs) in protease and reverse transcriptase were defined using the WHO HIV drug mutation list. Phylogenetic analysis was used to infer polymerase (pol) genotype.

RESULTS

TDRMs were found in 16/305 (5.2%) persons, 11 (3.6%) with resistance to NRTIs, 5 (1.6%) with resistance to NNRTIs and 3 (0.9%) with resistance to PIs. Dual-class TDRMs were found in three (1.0%) patients and one statistically supported cluster of TDRMs comprising two individuals with subtype B infection. TDRMs were found in 10 heterosexuals, 4 MSM and two intravenous drug users. Phylogenetic analyses identified high HIV-1 diversity consisting of mostly subtype B (44.6%), subtype C (3.3%), sub-subtype A1 (2.6%), sub-subtype F1 (2.3%), sub-subtype A-like (3.6%), subtype G (0.3%), CRF14_BG (1.6%), CRF05_DF (1.3%), CRF03_AB (0.3%) and unique recombinant forms (1.3%).

CONCLUSIONS

We found a low prevalence of TDR against a background of high HIV-1 genetic diversity among antiretroviral-naive patients in Bulgaria. Our results provide baseline data on TDR and support continued surveillance of high-risk populations in Bulgaria to better target treatment and prevention efforts.

Molecular tools for studying HIV transmission in sexual networks

PURPOSE OF REVIEW

Phylogenetics is frequently used for studies of population-based HIV transmission. The purpose of this review is to highlight the current utilities and limitations of phylogenetics in HIV epidemiological research from sample collection through to data analysis.

RECENT FINDINGS

Studies of HIV phylogenies can provide critical information about HIV epidemics that are otherwise difficult to obtain through traditional study design such as transmission of drug-resistant virus, mixing between demographic groups, and rapidity of viral spread within populations. However, recent results from empirical and theoretical studies of HIV phylogenies challenge some of the key assumptions and interpretations from phylogenetic studies. Recent findings include lack of transmission bottlenecks in MSM and injection drug use epidemics, evidence for preferential transmission of HIV virus in heterosexual epidemics, and limited evidence that tree topologies correlate directly with underlying network structures. Other challenges include a lack of a standardized definition for a phylogenetic transmission cluster and biased or sparse sampling of HIV transmission networks.

SUMMARY

Phylogenetics is an important tool for HIV research, and offers opportunities to understand critical aspects of the HIV epidemic. Like all epidemiological research, the methods used and interpretation of results from phylogenetic studies should be made cautiously with careful consideration.