Using phylogenetic analysis to trace HIV-1 migration among western European injecting drug users seroconverting from 1984 to 1997

Recent and Rapid Transmission of HIV Among People Who Inject Drugs in Scotland Revealed Through Phylogenetic Analysis

Background

Harm reduction has dramatically reduced HIV incidence among people who inject drugs (PWID). In Glasgow, Scotland, <10 infections/year have been diagnosed among PWID since the mid-1990s. However, in 2015 a sharp rise in diagnoses was noted among PWID; many were subtype C with 2 identical drug-resistant mutations and some displayed low avidity, suggesting the infections were linked and recent.

Methods

We collected Scottish pol sequences and identified closely related sequences from public databases. Genetic linkage was ascertained among 228 Scottish, 1820 UK, and 524 global sequences. The outbreak cluster was extracted to estimate epidemic parameters.

Results

All 104 outbreak sequences originated from Scotland and contained E138A and V179E. Mean genetic distance was <1% and mean time between transmissions was 6.7 months. The average number of onward transmissions consistently exceeded 1, indicating that spread was ongoing.

Conclusions

In contrast to other recent HIV outbreaks among PWID, harm reduction services were not clearly reduced in Scotland. Nonetheless, the high proportion of individuals with a history of homelessness (45%) suggests that services were inadequate for those in precarious living situations. The high prevalence of hepatitis C (>90%) is indicative of sharing of injecting equipment. Monitoring the epidemic phylogenetically in real time may accelerate public health action.

The application of HIV molecular epidemiology to public health

HIV is responsible for one of the largest viral pandemics in human history. Despite a concerted global response for prevention and treatment, the virus persists. Thus, urgent public health action, utilizing novel interventions, is needed to prevent future transmission events, critical to eliminating HIV. For public health planning to prove effective and successful, we need to understand the dynamics of regional epidemics and to intervene appropriately. HIV molecular epidemiology tools as implemented in phylogenetic, phylodynamic and phylogeographic analyses have proven to be powerful tools in public health planning across many studies. Numerous applications with HIV suggest that molecular methods alone or in combination with mathematical modelling can provide inferences about the transmission dynamics, critical epidemiological parameters (prevalence, incidence, effective number of infections, Re, generation times, time between infection and diagnosis), or the spatiotemporal characteristics of epidemics. Molecular tools have been used to assess the impact of an intervention and outbreak investigation which are of great public health relevance. In some settings, molecular sequence data may be more readily available than HIV surveillance data, and can therefore allow for molecular analyses to be conducted more easily. Nonetheless, classic methods have an integral role in monitoring and evaluation of public health programmes, and should supplement emerging techniques from the field of molecular epidemiology. Importantly, molecular epidemiology remains a promising approach in responding to viral diseases.

An epidemiological investigation to reconstruct a probable human immunodeficiency virus-1 transmission network: a case report

BACKGROUND

Recently published studies have highlighted the importance of phylogenetic and phylodynamic analyses in supporting epidemiological investigations to reconstruct the transmission network of human immunodeficiency virus. Here, we report a case of sexual transmission of human immunodeficiency virus type 1 between a man and a woman that marks once more the importance of a tightened collaboration between phylogeny and epidemiology.

CASE PRESENTATION

We describe a case of human immunodeficiency virus type 1 subtype B transmission in a stable Caucasian heterosexual couple. The man was 30 years old and the woman was 21 years old at the time of their presentation to the Department of Public Health and Infectious Diseases of the University of Rome "Sapienza". The couple reported a history of drug abuse.

CONCLUSION

Phylogenetic analysis is a powerful technique that if properly used can prove valuable in research investigations. In the case presented here, a phylogenetic analysis alongside epidemiological evidence allowed us to determine the most probable source of the human immunodeficiency virus infection. The dated tree allowed us to date the transmission event, the time point, and the direction of transmission based on the phylogeny, which agreed with the presumptive time of infection determined from clinical history-taking.

Distinct HIV type 1 strains in different risk groups and the absence of new infections by drug-resistant strains in Lithuania

To analyze HIV-1 genotypes in Lithuania and the transmission of drug-resistant viruses, HIV-1 sequences were obtained from 138 individuals, who were diagnosed as HIV-1 infected in 1990-2008 and represented all major risk groups. Subtype A strains, dominating in the former Soviet Union (90% of cases), were found in 60% of individuals, followed by subtype B (22%) and CRF03_AB (12%) strains. The remaining 7% of the strains included variants belonging to subtype C, CRF01_AE, CRF02_AG, more complex recombinant forms, and strains that could not be reliably genotyped. Analysis of virus genotypes per risk group revealed the circulation of distinct HIV-1 strains in different risk groups: subtype A viruses were present in 82% of injecting drug users (IDUs), but less than a half of heterosexually infected individuals and cases with unknown transmission route, and none of men having sex with men (MSM). We observed no mutations causing drug resistance among 27 newly diagnosed HIV-1 cases.

Transmission of risk-group specific HIV-1 strains among Dutch drug users for more than 20 years and their replacement by nonspecific strains after switching to low-harm drug practices

OBJECTIVES

To characterize HIV-1 epidemiological networks of men having sex with men (MSM) and drug users (DUs) in the Netherlands for >30 years.

DESIGN AND METHODS

Previously, we demonstrated different origin of the HIV-1 epidemics in Dutch MSM and DUs. To achieve the study objectives, risk group-specific genetic markers in the pol gene were examined in 315 participants of the Amsterdam Cohort Studies on HIV/AIDS who were registered as HIV-1 infected in 1981-2011.

RESULTS

Phylogenetic analysis demonstrated circulation of distinct virus strains in the 2 networks, with 98% of viruses of MSM clustering together and apart from strains of 73% DUs. Nine genetic markers that significantly distinguished virus strains specific for DUs were identified, of which 3 were ≥90% conserved. Over the total observation period, only 6% of viruses (4 of MSM and 14 of DUs) clustered with those of the other risk group. Among these sequences, the 3 most conserved genetic markers of that other risk group were 87% conserved.All 4 cases of DU-specific viruses among MSM occurred in 1980s-early 1990s. Viruses nonspecific for DUs were causing new infections among DUs at the rate of 20% till 2002 and replaced DU-specific strains among new infections thereafter, coinciding with switching of DUs to low-harm drug practices.

CONCLUSIONS

Dutch MSM and DUs have remained separate epidemiological networks for decades, despite their geographical and behavioral overlap. Switching to low-harm drug practices among DUs resulted in new infections caused by HIV-1 strains originating from other risk groups.

The HIV Epidemic: High-Income Countries

The HIV epidemic in higher-income nations is driven by receptive anal intercourse, injection drug use through needle/syringe sharing, and, less efficiently, vaginal intercourse. Alcohol and noninjecting drug use increase sexual HIV vulnerability. Appropriate diagnostic screening has nearly eliminated blood/blood product-related transmissions and, with antiretroviral therapy, has reduced mother-to-child transmission radically. Affected subgroups have changed over time (e.g., increasing numbers of Black and minority ethnic men who have sex with men). Molecular phylogenetic approaches have established historical links between HIV strains from central Africa to those in the United States and thence to Europe. However, Europe did not just receive virus from the United States, as it was also imported from Africa directly. Initial introductions led to epidemics in different risk groups in Western Europe distinguished by viral clades/sequences, and likewise, more recent explosive epidemics linked to injection drug use in Eastern Europe are associated with specific strains. Recent developments in phylodynamic approaches have made it possible to obtain estimates of sequence evolution rates and network parameters for epidemics.

Dynamics of two separate but linked HIV-1 CRF01_AE outbreaks among injection drug users in Stockholm, Sweden, and Helsinki, Finland

Detailed phylogenetic analyses were performed to characterize an HIV-1 outbreak among injection drug users (IDUs) in Stockholm, Sweden, in 2006. This study investigated the source and dynamics of HIV-1 spread during the outbreak as well as associated demographic and clinical factors. Seventy Swedish IDUs diagnosed during 2004 to 2007 were studied. Demographic, clinical, and laboratory data were collected, and the V3 region of the HIV-1 envelope gene was sequenced to allow detailed phylogenetic analyses. The results showed that the Stockholm outbreak was caused by a CRF01_AE variant imported from Helsinki, Finland, around 2003, which was quiescent until the outbreak started in 2006. Local Swedish subtype B variants continued to spread at a lower rate. The number of new CRF01_AE cases over a rooted phylogenetic tree accurately reflected the transmission dynamics and showed a temporary increase, by a factor of 12, in HIV incidence during the outbreak. Virus levels were similar in CRF01_AE and subtype B infections, arguing against differences in contagiousness. Similarly, there were no major differences in other baseline characteristics. Instead, the outbreak in Stockholm (and Helsinki) was best explained by an introduction of HIV into a standing network of previously uninfected IDUs. The combination of phylogenetics and epidemiological data creates a powerful tool for investigating outbreaks of HIV and other infectious diseases that could improve surveillance and prevention.

Population dynamics of HIV-1 subtype B in a cohort of men-having-sex-with-men in Rome, Italy

A recent increase in HIV diagnoses among men-having-sex-with-men (MSM) has been shown by surveillance data from Europe and Italy, and new approaches to inferring viral population dynamics from heterochronously sampled gene sequences have been developed. The aim of this study was to reconstruct the epidemiological history of HIV-1 subtype B in a homogeneous group of Italian MSM using a coalescent-based Bayesian framework. A total of 125 HIV-1 subtype B pol sequences were analyzed using Bayesian methods and a relaxed molecular clock to reconstruct their dated phylogeny and estimate population dynamics. At least 10 epidemiological clusters of 3-9 isolates were identified: half including the largest clades originated in the early 1990s and the other half radiated from 1999. Demographic analysis showed that the HIV epidemic grew in accordance with a logistic model characterized by a rapid exponential increase in the effective number of infections (r = 1.54 year) starting from the early 1980s and reaching a plateau 10 years later. Our data suggest that the HIV B epidemic entered our MSM population through multiple transmission chains about 20 years later than in other Western European country. Epidemiological clusters originating in the early 2000s suggest a recent re-emergence of HIV in Italian MSM.

Tracing the HIV-1 subtype B mobility in Europe: a phylogeographic approach

BACKGROUND

The prevalence and the origin of HIV-1 subtype B, the most prevalent circulating clade among the long-term residents in Europe, have been studied extensively. However the spatial diffusion of the epidemic from the perspective of the virus has not previously been traced.

RESULTS

In the current study we inferred the migration history of HIV-1 subtype B by way of a phylogeography of viral sequences sampled from 16 European countries and Israel. Migration events were inferred from viral phylogenies by character reconstruction using parsimony. With regard to the spatial dispersal of the HIV subtype B sequences across viral phylogenies, in most of the countries in Europe the epidemic was introduced by multiple sources and subsequently spread within local networks. Poland provides an exception where most of the infections were the result of a single point introduction. According to the significant migratory pathways, we show that there are considerable differences across Europe. Specifically, Greece, Portugal, Serbia and Spain, provide sources shedding HIV-1; Austria, Belgium and Luxembourg, on the other hand, are migratory targets, while for Denmark, Germany, Italy, Israel, Norway, the Netherlands, Sweden, Switzerland and the UK we inferred significant bidirectional migration. For Poland no significant migratory pathways were inferred.

CONCLUSION

Subtype B phylogeographies provide a new insight about the geographical distribution of viral lineages, as well as the significant pathways of virus dispersal across Europe, suggesting that intervention strategies should also address tourists, travellers and migrants.

Multiple HIV-1 introductions into the Swedish intravenous drug user population

In 2001, an increase of HIV-1 diagnoses among intravenous drug users (IVDU) was reported in Sweden. In nearby countries, Finland, Russia and the Baltic states, recent outbreaks had been described. Since there was a concern that these outbreaks would carry over to Sweden a study was initiated to determine the factors leading to the Swedish increase of HIV-1 diagnosed IVDUs. HIV-1 env V3 sequences were obtained from 97 patients enrolled in ongoing epidemiological studies encompassing the years 1987--2004 with a focus on 2001--2002. The sequences were used for maximum likelihood and Bayesian inference of the molecular epidemiology. Among the virus spreading in 2001--2002, we found that four different subtypes/CRFs were present in the Swedish IVDU population (A, B, CRF01_AE and CRF06_cpx). Subtype B constituted 85% of the infections, established by 12 independent introductions into the IVDU population. The worrisome increase in 2001 was mainly not a result of import of the outbreaks in nearby countries, but rather a higher detection rate of secondary cases due to efficient epidemiological tracing of the generally slow spread of established forms of subtype B in the IVDU community. However, a few of the non-subtype B cases were linked to the outbreaks in Finland, Estonia and Latvia. Because HIV-1 outbreaks can easily be exported from one country to another amongst IVDUs, this prompts continued surveillance in the Baltic Sea Region.

An outbreak of HIV-1 subtype G among Italian injecting drug users

We describe an outbreak of subtype G among injecting drug users (IDU) in northern Italy newly infected with HIV. We analysed pol gene sequences from samples of 139 individuals from different risk groups. Non-B subtypes were more frequently detected among IDU than in homosexual or heterosexual contacts. All G subtypes but one were found among IDU. The phylogenetic analysis indicated that the outbreak was of monophyletic origin and was caused by HIV-1 strains similar to those from western Africa.

Genetic and phylogenetic evolution of HIV-1 in a low subtype heterogeneity epidemic: the Italian example

The Human Immunodeficiency Virus type 1 (HIV-1) is classified into genetic groups, subtypes and sub-subtypes which show a specific geographic distribution pattern. The HIV-1 epidemic in Italy, as in most of the Western Countries, has traditionally affected the Intra-venous drug user (IDU) and Homosexual (Homo) risk groups and has been sustained by the genetic B subtype. In the last years, however, the HIV-1 transmission rate among heterosexuals has dramatically increased, becoming the prevalent transmission route. In fact, while the traditional risk groups have high levels of knowledge and avoid high-risk practices, the heterosexuals do not sufficiently perceive the risk of HIV-1 infection. This misperception, linked to the growing number of immigrants from non-Western Countries, where non-B clades and circulating recombinant forms (CRFs) are prevalent, is progressively introducing HIV-1 variants of non-B subtype in the Italian epidemic. This is in agreement with reports from other Western European Countries.

In this context, the Italian HIV-1 epidemic is still characterized by low subtype heterogeneity and represents a paradigmatic example of the European situation. The continuous molecular evolution of the B subtype HIV-1 isolates, characteristic of a long-lasting epidemic, together with the introduction of new subtypes as well as recombinant forms may have significant implications for diagnostic, treatment, and vaccine development. The study and monitoring of the genetic evolution of the HIV-1 represent, therefore, an essential strategy for controlling the local as well as global HIV-1 epidemic and for developing efficient preventive and therapeutic strategies.

Coinfection with HIV and hepatitis C virus among injection drug users in southern China

BACKGROUND

We sought to examine coinfection with human immunodeficiency virus (HIV) and hepatitis C virus (HCV) among injection drug users (IDUs) in Guangxi, China.

METHODS

A longitudinal cohort of IDUs (547 subjects) was established to study risk factors for bloodborne infections. At each visit, participants completed questionnaires defining demographic characteristics, patterns of drug use, and sexual behaviors. Blood samples were collected and analyzed for the presence and genotype of HIV and HCV.

RESULTS

Coinfection with HIV and HCV was found in 17.6% of the IDUs. HCV was present in 95.1% of HIV-positive and 70.4% of HIV-negative heroin users. The prevalence of HIV in HCV-positive and HCV-negative heroin users was 23.4% and 3.6%, respectively. Multivariate logistic regression analysis revealed that sexual activity during the past 6 months and duration of injection drug use were significantly associated with coinfection with HIV and HCV. The main circulating HCV genotypes included 6a (38%), 3b (37%), and 1a (19%), whereas genotypes 6e (4%), 3a (2%), and 1b (1%) were present in only a few IDUs. Multiple HCV genotypes were present at each study site and did not segregate by HIV status or subtype.

CONCLUSIONS

HCV is highly prevalent in IDUs throughout southern China. In Guangxi, HIV infections are the result of parenteral and sexual transmission, and, therefore, all IDUs are at high risk of coinfection with HIV and HCV. Molecular tracking of HCV may be a more sensitive predictor of the future spread of the HIV-1 epidemic than is HIV subtyping. This study emphasizes that, without implementation of injection prevention and primary substance abuse programs in China, the extent and effect of coinfection with HIV and HCV will only increase.

High prevalence of the CCR5Delta32 HIV-resistance mutation among Estonian HIV type 1-infected individuals

The aim of this survey was to investigate human immunodeficiency virus type 1 (HIV-1) coreceptor, chemokine receptor 5 (CCR5), polymorphism among Estonian HIV-1-infected individuals. Homozygous CCR5Delta32 genotypes have been associated with resistance to HIV-1 infection; however, inconsistent evidence exists as to whether a single copy of a mutant allele among heterozygotes confers protection from HIV-1 infection. In an Estonian population the frequency of the CCR5Delta32 allele has been found to be among the greatest observed to date. Ironically, Estonia is concomitantly characterized by a very high HIV-1 prevalence. We compared the allele frequencies in a healthy control population to the HIV-positive group. The frequency of heterozygous individuals did not differ significantly between the HIV-positive group and the control population. Allele frequencies were analyzed among different risk groups as well as groups with different HIV genetic backgrounds. We did not find a difference between CCR5Delta32 allele frequencies among intravenous drug users (IDUs) and sexually infected persons. Likewise, the distribution of CCR5Delta32 allele frequencies among patients infected with different subtypes did not differ while data from "pure" subtypes A, B, and CRF06_cpx were pooled and evaluated against unique recombinant forms.

HIV epidemiology in the Northwestern Federal District of Russia: dominance of HIV type 1 subtype A

A rapidly advancing epidemic of HIV-1 infection has been documented among injecting drug users (IDUs) in Russia. The Northwestern Federal District was the first of the seven Russian Federal Districts involved in a drug-related HIV epidemic through an outbreak in Kaliningrad in 1996. The Northwestern Federal District has a high HIV prevalence rate having reached 252 per 100,000 by the end of 2003, exceeding the Russian average (180) by 1.4 times. The epidemic peaked in 2001. Since then the annual number of new cases has decreased, probably reflecting saturation among at least some IDU populations. However, at the same time, the heterosexual spread of HIV has become more prominent. To study the genetic epidemiology of HIV-1, samples were collected from 150 individuals covering a wide geographical area and different transmission groups in the Northwestern Federal District. Phylogenetic analysis revealed that an Eastern European subtype A HIV-1 strain similar to those reported earlier among IDUs in other regions of Russia accounted for 80% of HIV-1 infections and was the predominant subtype in six out of seven administrative territories studied both among IDUs and heterosexually infected persons. As an exception to the dominant role of the Eastern European subtype A strain, the CRF03-AB strain was found to be dominant in the city of Cherepovets located in the north central European Russian territory of Vologda Oblast. This is the first report of the CRF03-AB strain causing an outbreak outside the Kaliningrad region.

Molecular testing of multiple HIV-1 transmissions in a criminal case

OBJECTIVE

To test the a priori hypothesis of HIV-1 transmission from one suspect to six recipients in a criminal case.

METHODS

Partial pol and/or env sequences were obtained for at least two samples of the suspect and the victims. Appropriate local controls were sampled based on epidemiological and subtype criteria. Phylogenetic testing was performed using different reconstruction methods.

RESULTS

Phylogenetic analyses consistently inferred a monophyletic cluster for the suspect and victim samples in both genome regions. This was highly supported by parametric and non-parametric bootstrapping techniques. Moreover, the controls most closely related to the suspect-victim cluster had a similar geographical origin to the suspect.

CONCLUSIONS

Taking into account the limitations on the conclusions that can be drawn from molecular investigations we could infer that our molecular data is consistent with a scenario of multiple HIV transmission between suspect and victims.

Molecular epidemiology of HTLV-2 infection among intravenous drug users in Spain

BACKGROUND

HTLV-2 infection is prevalent among intravenous drug users (IDUs), which often are coinfected with HIV-1. Two main subtypes, HTLV-2a and -2b, have been described among European IDUs, with a characteristically geographic distribution: IIa in northern Europe and IIb in the south. Although frequent traveling within Europe might have favoured a wide dissemination of diverse HTLV-2 variants, HTLV-2b was so far the only subtype reported in Spain.

OBJECTIVE

To investigate the current molecular epidemiology of HTLV-2 in Spain.

STUDY DESIGN

Twenty-six new HTLV-2 strains isolated from IDUs over the last 11 years in different geographic regions were examined. The HTLV LTR region (620bp) was sequenced and phylogenetic analyses were performed. Endonuclease restriction sites were examined to further characterize the HTLV-2 subgroup.

RESULTS

All 26 individuals infected with HTLV-2 showed viruses belonging to the HTLV-2 b4 clade and were closely related to the previously reported HTLV-2 Spanish and Italian IDU isolates.

CONCLUSIONS

HTLV-2 subtype b4 continues to be the only HTLV-2 subgroup recognized so far in Spain and no introduction of other HTLV-2 variants has occurred over the last 11 years.

HIV-1 co/super-infection in intravenous drug users

BACKGROUND

The frequency of HIV-1 co/super-infection is unknown despite their implications for public health and vaccine development. This issue was addressed during an epidemic of both CRF11 and B subtype among intravenous drug users (IVDUs).

METHODS

Bulk sequencing of reverse transcriptase, protease and C2V3 regions and subtype-specific nested polymerase chain reaction (PCR) in plasma and proviral DNA were performed using baseline and follow-up samples collected in recently infected IVDUs between 1998-2002 and in IVDUs with chronic infection living in the same area and presenting an unexpected rise of viremia (> 1 log10).

RESULTS

In 58 recently infected patients, three B/CRF-11 co-infections, 25 B, 28 CRF-11 and two other subtypes were detected at baseline. In the three co-infected patients, both CRF-11 and B were detected in plasma and proviral DNA and persisted during follow-up. B- and CFR-11-specific PCR performed on follow-up samples of 40 of 58 recently infected patients (median follow-up, 14.5 months) revealed a transient B super-infection in a patient initially infected by CRF-11. Five of 156 chronic IVDUs (total follow-up: 346 years) had an unexpected rise of viremia. In two of them, aviremic without treatment for years after an initial B infection, a symptomatic CRF-11 super-infection occurred and was associated with high viral load and a fall of CD4 cell count.

CONCLUSIONS

In recently infected IVDUs, co-infection B/CRF-11 is relatively frequent (5%). In chronically infected IVDUs super-infection may be transient and may occur in patients controlling efficiently HIV infection by the initial strain.

Spreading of HIV-1 subtype G and envB/gagG recombinant strains among injecting drug users in Lisbon, Portugal

We have evaluated the genetic diversity of HIV-1 strains infecting injecting drug users (IDUs) in Lisbon, Portugal. Heteroduplex mobility assay and/or phylogenetic analysis revealed that env (C2V3C3 or gp41) subtype B is present in 63.7% of the 135 viral samples studied, followed by subtypes G (23.7%), A (6.7%), F (5.2%), and D (0.7%). Similar analysis of gag (p24/p7) performed on 91 of the specimens demonstrated that 49.5% of the infections were caused by subtype G viruses; other gag subtypes identified were B (39.5%), F (3.3%), A and D (1.1.% each), and the recombinant circulating form CRF02_AG (5.5%). Discordant env/gag sub-types were detected in 34.1% of the strains and may reflect the presence of dual infections and/or recombinant viruses. The presumptive B/G recombinant form was highly predominant (21 of 31). The genetic pattern of HIV-1 subtype B and G strains is suggestive of multiple introductions and recombination episodes and of a longstanding presence of both subtypes in the country. C2V3C3 amino acid sequences from IDU-derived subtype G viruses presented highly significant signatures, which distinguish the variants from this transmission group. The unusually high prevalence of subtype G sequences (34.1%), independent of the geographic origin of the infected individuals, makes this IDU HIV-1 epidemic unique.

Travel and the spread of HIV-1 genetic variants

HIV-1 comprises three groups, the main (M group), O (outlier) and N (non-M, non-O). The M group, divided into 11 subtypes, is responsible for the global HIV-1 pandemic. Recombination between M subtypes has resulted in the generation of multiple circulating recombinant forms (CRFs) consisting of mosaic lineages. Most subtypes and CRFs are represented in Africa, whereas predominance of one or a few subtypes was reported initially elsewhere. This finding reflects the African origin of the epidemic. In western countries, where the B subtype is predominant, there is a steep increase in non B-subtypes and CRFs, while new recombinants emerge worldwide. Travellers contribute to the spread of HIV-1 genetic diversity worldwide, and in the developing world migration of rural populations and civil war are additional contributing factors. The spreading of HIV-1 variants has implications for diagnostic, treatment, and vaccine development.