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

Inferring transmission heterogeneity using virus genealogies: Estimation and targeted prevention

Spread of HIV typically involves uneven transmission patterns where some individuals spread to a large number of individuals while others to only a few or none. Such transmission heterogeneity can impact how fast and how much an epidemic spreads. Further, more efficient interventions may be achieved by taking such transmission heterogeneity into account. To address these issues, we developed two phylogenetic methods based on virus sequence data: 1) to generally detect if significant transmission heterogeneity is present, and 2) to pinpoint where in a phylogeny high-level spread is occurring. We derive inference procedures to estimate model parameters, including the amount of transmission heterogeneity, in a sampled epidemic. We show that it is possible to detect transmission heterogeneity under a wide range of simulated situations, including incomplete sampling, varying levels of heterogeneity, and including within-host genetic diversity. When evaluating real HIV-1 data from different epidemic scenarios, we found a lower level of transmission heterogeneity in slowly spreading situations and a higher level of heterogeneity in data that included a rapid outbreak, while R0 and Sackin's index (overall tree shape statistic) were similar in the two scenarios, suggesting that our new method is able to detect transmission heterogeneity in real data. We then show by simulations that targeted prevention, where we pinpoint high-level spread using a coalescence measurement, is efficient when sequence data are collected in an ongoing surveillance system. Such phylogeny-guided prevention is efficient under both single-step contact tracing as well as iterative contact tracing as compared to random intervention.

Recent increased identification and transmission of HIV-1 unique recombinant forms in Sweden

A temporal increase in non-B subtypes has earlier been described in Sweden by us and we hypothesized that this increased viral heterogeneity may become a hotspot for the development of more complex and unique recombinant forms (URFs) if the epidemics converge. In the present study, we performed subtyping using four automated tools and phylogenetic analysis by RAxML of pol gene sequences (n = 5246) and HIV-1 near full-length genome (HIV-NFLG) sequences (n = 104). A CD4⁺ T-cell decline trajectory algorithm was used to estimate time of HIV infection. Transmission clusters were identified using the family-joining method. The analysis of HIV-NFLG and pol gene described 10.6% (11/104) and 2.6% (137/5246) of the strains as URFs, respectively. An increasing trend of URFs was observed in recent years by both approaches (p = 0·0082; p < 0·0001). Transmission cluster analysis using the pol gene of all URFs identified 14 clusters with two to eight sequences. Larger transmission clusters of URFs (BF1 and 01B) were observed among MSM who mostly were sero-diagnosed in recent time. Understanding the increased appearance and transmission of URFs in recent years could have importance for public health interventions and the use of HIV-NFLG would provide better statistical support for such assessments.

Inference of Transmission Network Structure from HIV Phylogenetic Trees

Phylogenetic inference is an attractive means to reconstruct transmission histories and epidemics. However, there is not a perfect correspondence between transmission history and virus phylogeny. Both node height and topological differences may occur, depending on the interaction between within-host evolutionary dynamics and between-host transmission patterns. To investigate these interactions, we added a within-host evolutionary model in epidemiological simulations and examined if the resulting phylogeny could recover different types of contact networks. To further improve realism, we also introduced patient-specific differences in infectivity across disease stages, and on the epidemic level we considered incomplete sampling and the age of the epidemic. Second, we implemented an inference method based on approximate Bayesian computation (ABC) to discriminate among three well-studied network models and jointly estimate both network parameters and key epidemiological quantities such as the infection rate. Our ABC framework used both topological and distance-based tree statistics for comparison between simulated and observed trees. Overall, our simulations showed that a virus time-scaled phylogeny (genealogy) may be substantially different from the between-host transmission tree. This has important implications for the interpretation of what a phylogeny reveals about the underlying epidemic contact network. In particular, we found that while the within-host evolutionary process obscures the transmission tree, the diversification process and infectivity dynamics also add discriminatory power to differentiate between different types of contact networks. We also found that the possibility to differentiate contact networks depends on how far an epidemic has progressed, where distance-based tree statistics have more power early in an epidemic. Finally, we applied our ABC inference on two different outbreaks from the Swedish HIV-1 epidemic.

Dynamic Variation in Sexual Contact Rates in a Cohort of HIV-Negative Gay Men

Human immunodeficiency virus (HIV) transmission models that include variability in sexual behavior over time have shown increased incidence, prevalence, and acute-state transmission rates for a given population risk profile. This raises the question of whether dynamic variation in individual sexual behavior is a real phenomenon that can be observed and measured. To study this dynamic variation, we developed a model incorporating heterogeneity in both between-person and within-person sexual contact patterns. Using novel methodology that we call iterated filtering for longitudinal data, we fitted this model by maximum likelihood to longitudinal survey data from the Centers for Disease Control and Prevention's Collaborative HIV Seroincidence Study (1992-1995). We found evidence for individual heterogeneity in sexual behavior over time. We simulated an epidemic process and found that inclusion of empirically measured levels of dynamic variation in individual-level sexual behavior brought the theoretical predictions of HIV incidence into closer alignment with reality given the measured per-act probabilities of transmission. The methods developed here provide a framework for quantifying variation in sexual behaviors that helps in understanding the HIV epidemic among gay men.

Global Dispersal Pattern of HIV Type 1 Subtype CRF01_AE: A Genetic Trace of Human Mobility Related to Heterosexual Sexual Activities Centralized in Southeast Asia

BACKGROUND

Human immunodeficiency virus type 1 (HIV-1) subtype CRF01_AE originated in Africa and then passed to Thailand, where it established a major epidemic. Despite the global presence of CRF01_AE, little is known about its subsequent dispersal pattern.

METHODS

We assembled a global data set of 2736 CRF01_AE sequences by pooling sequences from public databases and patient-cohort studies. We estimated viral dispersal patterns, using statistical phylogeographic analysis run over bootstrap trees estimated by the maximum likelihood method.

RESULTS

We show that Thailand has been the source of viral dispersal to most areas worldwide, including 17 of 20 sampled countries in Europe. Japan, Singapore, Vietnam, and other Asian countries have played a secondary role in the viral dissemination. In contrast, China and Taiwan have mainly imported strains from neighboring Asian countries, North America, and Africa without any significant viral exportation.

DISCUSSION

The central role of Thailand in the global spread of CRF01_AE can be probably explained by the popularity of Thailand as a vacation destination characterized by sex tourism and by Thai emigration to the Western world. Our study highlights the unique case of CRF01_AE, the only globally distributed non-B clade whose global dispersal did not originate in Africa.

Transmission of hepatitis C virus among intravenous drug users in the Uppsala region of Sweden

Background

Epidemiology and transmission patterns of hepatitis C virus (HCV) are important subjects as we enter a new era of treatment with directly acting antivirals (DAAs). The highest prevalence of HCV in developed countries is found among intravenous drug users (IDUs), where unsafe needle sharing practices provide the main route of infection. Efforts to prohibit the continuous spread of HCV among these groups have been initiated by the community services and health care providers. Our goal was to understand how HCV was transmitted among IDUs within a limited population group. We provide a retrospective study (2005–2007) of the HCV transmission patterns in a population of IDUs in the Uppsala region of Sweden.

Method

Eighty-two serum samples were collected from IDUs in Uppsala County. Our reverse transcription nested polymerase chain reaction (RT-nested PCR) and sequencing method enabled a comprehensive genetic analysis for a broad spectrum of genotypes of two relatively conserved regions, NS5B and NS3, that encodes for the viral polymerase and protease, respectively. HCV RNA in serum samples was amplified and sequenced with in-house primers. Sequence similarities between individuals and subgroups were analyzed with maximum likelihood (ML) phylogenetic trees. Published HCV reference sequences from other geographic regions and countries were also included for clarity.

Results

Phylogenetic analysis was possible for 59 NS5B (72%) and 29 NS3 (35%) sequences from Uppsala patients. Additionally, we also included 15 NS3 sequences from Örebro patients, making a total of 44 NS3 sequences for the analysis. By analyzing the NS3 sequences, two transmission sets were found between the IDUs (>98% sequence identity), with one set consisting of two individuals and another set consisting of three individuals. In addition, the phylogenetic analysis done with our serum samples displayed clusters that distinguished them from the reference sequences.

Conclusion

Our method seems to enable us to trace the HCV transmission between IDUs. Furthermore, the method is fairly independent of the time of infection because the method uses relatively conserved HCV sequence regions (i.e. NS5B and NS3).

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.

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.